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

Defining myocardial fiber bundle architecture in atrial digital twins

A key component in developing atrial digital twins (ADT) - virtual representations of patients' atria - is the accurate prescription of myocardial fibers which are essential for the tissue characterization. Due to the difficulty of reconstructing atrial fibers from medical imaging, a widely used strategy for fiber generation in ADT relies on mathematical models. Existing methodologies utilize semi-automatic approaches, are tailored to specific morphologies, and lack rigorous validation against imaging fiber data. In this study, we introduce a novel atrial Laplace-Dirichlet-Rule-Based Method (LDRBM) for prescribing highly detailed myofiber orientations and providing robust regional annotation in bi-atrial morphologies of any complexity. The robustness of our approach is verified in eight extremely detailed bi-atrial geometries, derived from a sub-millimeter Diffusion-Tensor-Magnetic-Resonance Imaging (DTMRI) human atrial fiber dataset. We validate the LDRBM by quantitatively recreating each of the DTMRI fiber architectures: a comprehensive comparison with DTMRI ground truth data is conducted, investigating differences between electrophysiology (EP) simulations provided by either LDRBM and DTMRI fibers. Finally, we demonstrate that the novel LDRBM outperforms current state-of-the-art (LDRBMs and Universal Atrial Coordinates) fiber models, confirming the exceptional accuracy of our methodology and the critical importance of incorporating detailed fiber orientations in EP simulations. Ultimately, this work represents a fundamental step towards the development of physics-based digital twins of the human atria, establishing a new standard for prescribing fibers in ADT.

A computationally efficient anisotropic electrophysiological multiscale uterus model: From cell to organ and myometrium to abdominal surface

BACKGROUND AND OBJECTIVE

Preterm labor is a global problem affecting the health of newborns. Despite numerous studies reporting electrophysiological changes throughout pregnancy, the underlying mechanism that triggers labor remains unclear. Electrophysiological modeling can provide additional information to better understand the physiological transition from pregnancy to labor. Previous uterine electrophysiological models do not consider either the tissue thickness or fiber structure, which have both been shown to significantly impact propagation patterns.

METHODS

This paper presents a parallel computational model of the uterus using the bioengineering modeling environment OpenCMISS. This model is a multiscale anisotropic model that spans different levels from cell to organ. At the cellular level, the model utilizes a mathematical representation of uterine myocytes based on multiple ion channels. In the 3D uterine model, fiber structures are added, ranging from horizontal rings in the inner layer to vertically downward fibers in the outer layer, to more accurately depict the electrophysiological activities of the uterus. Additionally, we have developed a multilayer volume conduction model based on the boundary element method to describe the propagation of electrical signals from the myometrium to the abdominal surface.

RESULTS

Our model can not only reproduce faithfully both local non-propagated and global propagated electrical activity, but also simulate the fast wave low and fast wave high components of the electrohysterogram (EHG) on the abdominal surface. The model results support the hypothesis that the fast wave high of the EHG signal is related to uterine excitability and fast wave low is related to signal propagation. The amplitude of the simulated signal on the abdominal surface falls in the ranges of real EHG data, which is inversely proportional to the abdominal subcutaneous fat thickness, and the signal waveform highly depends on electrode position and the relative distance to the pacemaker. In addition, the propagation velocity is highly dependent on the uterus geometry and falls in the real-world data range CONCLUSIONS: Our models facilitate a better understanding of the electrophysiological changes of the uterus during pregnancy and labor, and allow for an investigation of drug effects and/or structural or anatomical abnormalities.

Effects of spatially dense adrenergic stimulation to rotor behaviour in simulated atrial sheets

Sympathetic hyperactivity via spatially dense adrenergic stimulation may create pro-arrhythmic substrates even without structural remodelling. However, the effect of sympathetic hyperactivity on arrhythmic activity, such as rotors, is unknown. Using simulations, we examined the effects of gradually increasing the spatial density of adrenergic stimulation (AS) in atrial sheets on rotors. We compared their characteristics against rotors hosted in atrial sheets with increasing spatial density of minimally conductive (MC) elements to simulate structural remodelling due to injury or disease. We generated rotors using an S1-S2 stimulation protocol. Then, we created phase maps to identify phase singularities and map their trajectory over time. We measured each rotor's duration (s), angular speed (rad/s), and spatiotemporal organization. We demonstrated that atrial sheets with increased AS spatial densities could maintain rotors longer than with MC elements (2.6 ± 0.1 s vs. 1.5 ± 0.2 s, p<0.001). Moreover, rotors have higher angular speed (70 ± 7 rads/s vs. 60 ± 15 rads/s, p<0.05) and better spatiotemporal organization (0.56 ± 0.05 vs. 0.58 ± 0.18, p<0.05) in atrial sheets with less than 25% AS elements compared to MC elements. Our findings may help elucidate electrophysiological potential alterations in atrial substrates due to sympathetic hyperactivity, particularly among individuals with autonomic derangements caused by chronic distress.

A 3D Statistical Shape Model of the Right Ventricular Outflow Tract in Pulmonary Valve Replacement Patients Post-Surgical Repair

Assessment of the right ventricular outflow tract and pulmonary arteries (RVOT) for percutaneous pulmonary valve implantation (PPVI) uses discrete measurements (diameters and lengths) from medical images. This multi-centre study identified the 3D RVOT shape features prevalent in patients late after surgical repair of congenital heart disease (CHD). A 3D RVOT statistical shape model (SSM) was computed from 81 retrospectively selected CHD patients (14.7 ± 6.8 years) who required pulmonary valve replacement late after surgical repair. A principal component analysis identified prevalent shape features (modes) within the population which were compared with standard geometric measurements (diameter, length and surface area) and between sub-groups of diagnosis, RVOT type and dysfunction. Shape mode 1 and 2 represented RVOT size and curvature and tapering and length, respectively. Shape modes 3-5 related to branch pulmonary artery calibre, conical vs. bulbous RVOTs and RVOT curvature, respectively. Tetralogy of Fallot, transannular patch type and regurgitant RVOTs were larger and straighter while conduit and stenotic types were longer and more cylindrical than other subgroups. This SSM analysed the main 3D shape features present in a population of RVOTs, exploiting the wide 3D anatomical information provided by routine imaging. This morphological information may have implications for PPVI patient selection and device design.

Leveraging 3D Atrial Geometry for the Evaluation of Atrial Fibrillation: A Comprehensive Review

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia associated with significant morbidity and mortality. Managing risk of stroke and AF burden are pillars of AF management. Atrial geometry has long been recognized as a useful measure in achieving these goals. However, traditional diagnostic approaches often overlook the complex spatial dynamics of the atria. This review explores the emerging role of three-dimensional (3D) atrial geometry in the evaluation and management of AF. Advancements in imaging technologies and computational modeling have enabled detailed reconstructions of atrial anatomy, providing insights into the pathophysiology of AF that were previously unattainable. We examine current methodologies for interpreting 3D atrial data, including qualitative, basic quantitative, global quantitative, and statistical shape modeling approaches. We discuss their integration into clinical practice, highlighting potential benefits such as personalized treatment strategies, improved outcome prediction, and informed treatment approaches. Additionally, we discuss the challenges and limitations associated with current approaches, including technical constraints and variable interpretations, and propose future directions for research and clinical applications. This comprehensive review underscores the transformative potential of leveraging 3D atrial geometry in the evaluation and management of AF, advocating for its broader adoption in clinical practice.

Novel Domain Knowledge-Encoding Algorithm Enables Label-Efficient Deep Learning for Cardiac CT Segmentation to Guide Atrial Fibrillation Treatment in a Pilot Dataset

Background: Segmenting computed tomography (CT) is crucial in various clinical applications, such as tailoring personalized cardiac ablation for managing cardiac arrhythmias. Automating segmentation through machine learning (ML) is hindered by the necessity for large, labeled training data, which can be challenging to obtain. This article proposes a novel approach for automated, robust labeling using domain knowledge to achieve high-performance segmentation by ML from a small training set. The approach, the domain knowledge-encoding (DOKEN) algorithm, reduces the reliance on large training datasets by encoding cardiac geometry while automatically labeling the training set. The method was validated in a hold-out dataset of CT results from an atrial fibrillation (AF) ablation study. Methods: The DOKEN algorithm parses left atrial (LA) structures, extracts "anatomical knowledge" by leveraging digital LA models (available publicly), and then applies this knowledge to achieve high ML segmentation performance with a small number of training samples. The DOKEN-labeled training set was used to train a nnU-Net deep neural network (DNN) model for segmenting cardiac CT in N = 20 patients. Subsequently, the method was tested in a hold-out set with N = 100 patients (five times larger than training set) who underwent AF ablation. Results: The DOKEN algorithm integrated with the nn-Unet model achieved high segmentation performance with few training samples, with a training to test ratio of 1:5. The Dice score of the DOKEN-enhanced model was 96.7% (IQR: 95.3% to 97.7%), with a median error in surface distance of boundaries of 1.51 mm (IQR: 0.72 to 3.12) and a mean centroid-boundary distance of 1.16 mm (95% CI: -4.57 to 6.89), similar to expert results (r = 0.99; p < 0.001). In digital hearts, the novel DOKEN approach segmented the LA structures with a mean difference for the centroid-boundary distances of -0.27 mm (95% CI: -3.87 to 3.33; r = 0.99; p < 0.0001). Conclusions: The proposed novel domain knowledge-encoding algorithm was able to perform the segmentation of six substructures of the LA, reducing the need for large training data sets. The combination of domain knowledge encoding and a machine learning approach could reduce the dependence of ML on large training datasets and could potentially be applied to AF ablation procedures and extended in the future to other imaging, 3D printing, and data science applications.

Statistical shape modeling of the geometric morphology of the canine femur, tibia, and patella

Bone morphometry varies among dogs of different sizes and breeds. Studying these differences may help understand the predisposition of certain breeds for specific orthopedic pathologies. This study aimed to develop a statistical shape model (SSM) of the femur, patella, and tibia of dogs without any clinical orthopeadic abnormalities to analyze and compare morphological variations based on body weight and breed. A total of 97 CT scans were collected from different facilities and divided based on breed and body weight. The 3D models of the bones were obtained and aligned to a coordinate system. The SSM was created using principal component analysis (PCA) to analyze shape variations. The study found that the first few modes of variation accounted for a significant percentage of the total variation, with size/scale being the most prominent factor. The results provide valuable insights into normal anatomical variations and can be used for future research in understanding pathological bone morphologies and developing 3D imaging algorithms in veterinary medicine.

In Silico TRials guide optimal stratification of ATrIal FIbrillation patients to Catheter Ablation and pharmacological medicaTION: the i-STRATIFICATION study

AIMS

Patients with persistent atrial fibrillation (AF) experience 50% recurrence despite pulmonary vein isolation (PVI), and no consensus is established for secondary treatments. The aim of our i-STRATIFICATION study is to provide evidence for stratifying patients with AF recurrence after PVI to optimal pharmacological and ablation therapies, through in silico trials.

METHODS AND RESULTS

A cohort of 800 virtual patients, with variability in atrial anatomy, electrophysiology, and tissue structure (low-voltage areas, LVAs), was developed and validated against clinical data from ionic currents to electrocardiogram. Virtual patients presenting AF post-PVI underwent 12 secondary treatments. Sustained AF developed in 522 virtual patients after PVI. Second ablation procedures involving left atrial ablation alone showed 55% efficacy, only succeeding in the small right atria (<60 mL). When additional cavo-tricuspid isthmus ablation was considered, Marshall-PLAN sufficed (66% efficacy) for the small left atria (<90 mL). For the bigger left atria, a more aggressive ablation approach was required, such as anterior mitral line (75% efficacy) or posterior wall isolation plus mitral isthmus ablation (77% efficacy). Virtual patients with LVAs greatly benefited from LVA ablation in the left and right atria (100% efficacy). Conversely, in the absence of LVAs, synergistic ablation and pharmacotherapy could terminate AF. In the absence of ablation, the patient's ionic current substrate modulated the response to antiarrhythmic drugs, being the inward currents critical for optimal stratification to amiodarone or vernakalant.

CONCLUSION

In silico trials identify optimal strategies for AF treatment based on virtual patient characteristics, evidencing the power of human modelling and simulation as a clinical assisting tool.

Developing Cardiac Digital Twins at Scale: Insights from Personalised Myocardial Conduction Velocity

Large-cohort studies using cardiovascular imaging and diagnostic datasets have assessed cardiac anatomy, function, and outcomes, but typically do not reveal underlying biological mechanisms. Cardiac digital twins (CDTs) provide personalized physics- and physiology-constrained in-silico representations, enabling inference of multi-scale properties tied to these mechanisms. We constructed 3464 anatomically-accurate CDTs using cardiac magnetic resonance images from UK biobank and personalised their myocardial conduction velocities (CVs) from electrocardiograms (ECG), through an automated framework. We found well-known sex-specific differences in QRS duration were fully explained by myocardial anatomy, as CV remained consistent across sexes. Conversely, significant associations of CV with ageing and increased BMI suggest myocardial tissue remodelling. Novel associations were observed with left ventricular ejection fraction and mental-health phenotypes, through a phenome-wide association study, and CV was also linked with adverse clinical outcomes. Our study highlights the utility of population-based CDTs in assessing intersubject variability and uncovering strong links with mental health.

Non-Invasive Electroanatomical Mapping: A State-Space Approach for Myocardial Current Density Estimation

Electroanatomical mapping is a method for creating a model of the electrophysiology of the human heart. Medical professionals routinely locate and ablate the site of origin of cardiac arrhythmias with invasive catheterization. Non-invasive localization takes the form of electrocardiographic (ECG) or magnetocardiographic (MCG) imaging, where the goal is to reconstruct the electrical activity of the human heart. Non-invasive alternatives to catheter electroanatomical mapping would reduce patients' risks and open new venues for treatment planning and prevention. This work introduces a new system state-based method for estimating the electrical activity of the human heart from MCG measurements. Our model enables arbitrary propagation paths and velocities. A Kalman filter optimally estimates the current densities under the given measurements and model parameters. In an outer optimization loop, these model parameters are then optimized via gradient descent. This paper aims to establish the foundation for future research by providing a detailed mathematical explanation of the algorithm. We demonstrate the feasibility of our method through a simplified one-layer simulation. Our results show that the algorithm can learn the propagation paths from the magnetic measurements. A threshold-based segmentation into healthy and pathological tissue yields a DICE score of 0.84, a recall of 0.77, and a precision of 0.93.

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.

AugmentA: Patient-specific augmented atrial model generation tool

Digital twins of patients' hearts are a promising tool to assess arrhythmia vulnerability and to personalize therapy. However, the process of building personalized computational models can be challenging and requires a high level of human interaction. We propose a patient-specific Augmented Atria generation pipeline (AugmentA) as a highly automated framework which, starting from clinical geometrical data, provides ready-to-use atrial personalized computational models. AugmentA identifies and labels atrial orifices using only one reference point per atrium. If the user chooses to fit a statistical shape model to the input geometry, it is first rigidly aligned with the given mean shape before a non-rigid fitting procedure is applied. AugmentA automatically generates the fiber orientation and finds local conduction velocities by minimizing the error between the simulated and clinical local activation time (LAT) map. The pipeline was tested on a cohort of 29 patients on both segmented magnetic resonance images (MRI) and electroanatomical maps of the left atrium. Moreover, the pipeline was applied to a bi-atrial volumetric mesh derived from MRI. The pipeline robustly integrated fiber orientation and anatomical region annotations in 38.4 ± 5.7 s. In conclusion, AugmentA offers an automated and comprehensive pipeline delivering atrial digital twins from clinical data in procedural time.

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.

Differences in atrial substrate localization using late gadolinium enhancement-magnetic resonance imaging, electrogram voltage, and conduction velocity: a cohort study using a consistent anatomical reference frame in patients with persistent atrial fibrillation

AIMS

Electro-anatomical voltage, conduction velocity (CV) mapping, and late gadolinium enhancement (LGE) magnetic resonance imaging (MRI) have been correlated with atrial cardiomyopathy (ACM). However, the comparability between these modalities remains unclear. This study aims to (i) compare pathological substrate extent and location between current modalities, (ii) establish spatial histograms in a cohort, (iii) develop a new estimated optimized image intensity threshold (EOIIT) for LGE-MRI identifying patients with ACM, (iv) predict rhythm outcome after pulmonary vein isolation (PVI) for persistent atrial fibrillation (AF).

METHODS AND RESULTS

Thirty-six ablation-naive persistent AF patients underwent LGE-MRI and high-definition electro-anatomical mapping in sinus rhythm. Late gadolinium enhancement areas were classified using the UTAH, image intensity ratio (IIR >1.20), and new EOIIT method for comparison to low-voltage substrate (LVS) and slow conduction areas <0.2 m/s. Receiver operating characteristic analysis was used to determine LGE thresholds optimally matching LVS. Atrial cardiomyopathy was defined as LVS extent ≥5% of the left atrium (LA) surface at <0.5 mV. The degree and distribution of detected pathological substrate (percentage of individual LA surface are) varied significantly (P < 0.001) across the mapping modalities: 10% (interquartile range 0-14%) of the LA displayed LVS <0.5 mV vs. 7% (0-12%) slow conduction areas <0.2 m/s vs. 15% (8-23%) LGE with the UTAH method vs. 13% (2-23%) using IIR >1.20, with most discrepancies on the posterior LA. Optimized image intensity thresholds and each patient's mean blood pool intensity correlated linearly (R2 = 0.89, P < 0.001). Concordance between LGE-MRI-based and LVS-based ACM diagnosis improved with the novel EOIIT applied at the anterior LA [83% sensitivity, 79% specificity, area under the curve (AUC): 0.89] in comparison to the UTAH method (67% sensitivity, 75% specificity, AUC: 0.81) and IIR >1.20 (75% sensitivity, 62% specificity, AUC: 0.67).

CONCLUSION

Discordances in detected pathological substrate exist between LVS, CV, and LGE-MRI in the LA, irrespective of the LGE detection method. The new EOIIT method improves concordance of LGE-MRI-based ACM diagnosis with LVS in ablation-naive AF patients but discrepancy remains particularly on the posterior wall. All methods may enable the prediction of rhythm outcomes after PVI in patients with persistent AF.

A systematic review of cardiac in-silico clinical trials

Computational models of the heart are now being used to assess the effectiveness and feasibility of interventions through in-silico clinical trials (ISCTs). As the adoption and acceptance of ISCTs increases, best practices for reporting the methodology and analysing the results will emerge. Focusing in the area of cardiology, we aim to evaluate the types of ISCTs, their analysis methods and their reporting standards. To this end, we conducted a systematic review of cardiac ISCTs over the period of 1 January 2012-1 January 2022, following the preferred reporting items for systematic reviews and meta-analysis (PRISMA). We considered cardiac ISCTs of human patient cohorts, and excluded studies of single individuals and those in which models were used to guide a procedure without comparing against a control group. We identified 36 publications that described cardiac ISCTs, with most of the studies coming from the US and the UK. In 75% of the studies, a validation step was performed, although the specific type of validation varied between the studies. ANSYS FLUENT was the most commonly used software in 19% of ISCTs. The specific software used was not reported in 14% of the studies. Unlike clinical trials, we found a lack of consistent reporting of patient demographics, with 28% of the studies not reporting them. Uncertainty quantification was limited, with sensitivity analysis performed in only 19% of the studies. In 97% of the ISCTs, no link was provided to provide easy access to the data or models used in the study. There was no consistent naming of study types with a wide range of studies that could potentially be considered ISCTs. There is a clear need for community agreement on minimal reporting standards on patient demographics, accepted standards for ISCT cohort quality control, uncertainty quantification, and increased model and data sharing.

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.

A personalized real-time virtual model of whole heart electrophysiology

Computer models capable of representing the intrinsic personal electrophysiology (EP) of the heart in silico are termed virtual heart technologies. When anatomy and EP are tailored to individual patients within the model, such technologies are promising clinical and industrial tools. Regardless of their vast potential, few virtual technologies simulating the entire organ-scale EP of all four-chambers of the heart have been reported and widespread clinical use is limited due to high computational costs and difficulty in validation. We thus report on the development of a novel virtual technology representing the electrophysiology of all four-chambers of the heart aiming to overcome these limitations. In our previous work, a model of ventricular EP embedded in a torso was constructed from clinical magnetic resonance image (MRI) data and personalized according to the measured 12 lead electrocardiogram (ECG) of a single subject under normal sinus rhythm. This model is then expanded upon to include whole heart EP and a detailed representation of the His-Purkinje system (HPS). To test the capacities of the personalized virtual heart technology to replicate standard clinical morphological ECG features under such conditions, bundle branch blocks within both the right and the left ventricles under two different conduction velocity settings are modeled alongside sinus rhythm. To ensure clinical viability, model generation was completely automated and simulations were performed using an efficient real-time cardiac EP simulator. Close correspondence between the measured and simulated 12 lead ECG was observed under normal sinus conditions and all simulated bundle branch blocks manifested relevant clinical morphological features.

Multi-Domain Variational Autoencoders for Combined Modeling of MRI-Based Biventricular Anatomy and ECG-Based Cardiac Electrophysiology

Human cardiac function is characterized by a complex interplay of mechanical deformation and electrophysiological conduction. Similar to the underlying cardiac anatomy, these interconnected physiological patterns vary considerably across the human population with important implications for the effectiveness of clinical decision-making and the accuracy of computerized heart models. While many previous works have investigated this variability separately for either cardiac anatomy or physiology, this work aims to combine both aspects in a single data-driven approach and capture their intricate interdependencies in a multi-domain setting. To this end, we propose a novel multi-domain Variational Autoencoder (VAE) network to capture combined Electrocardiogram (ECG) and Magnetic Resonance Imaging (MRI)-based 3D anatomy information in a single model. Each VAE branch is specifically designed to address the particular challenges of the respective input domain, enabling efficient encoding, reconstruction, and synthesis of multi-domain cardiac signals. Our method achieves high reconstruction accuracy on a United Kingdom Biobank dataset, with Chamfer Distances between reconstructed and input anatomies below the underlying image resolution and ECG reconstructions outperforming multiple single-domain benchmarks by a considerable margin. The proposed VAE is capable of generating realistic virtual populations of arbitrary size with good alignment in clinical metrics between the synthesized and gold standard anatomies and Maximum Mean Discrepancy (MMD) scores of generated ECGs below those of comparable single-domain approaches. Furthermore, we observe the latent space of our VAE to be highly interpretable with separate components encoding different aspects of anatomical and ECG variability. Finally, we demonstrate that the combined anatomy and ECG representation improves the performance in a cardiac disease classification task by 3.9% in terms of Area Under the Receiver Operating Characteristic (AUROC) curve over the best corresponding single-domain modeling approach.

A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms

Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research.

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.