Towards personalized clinical in-silico modeling of atrial anatomy and electrophysiology

Impact of effective refractory period personalization on arrhythmia vulnerability in patient-specific atrial computer models

AIMS

The effective refractory period (ERP) is one of the main electrophysiological properties governing arrhythmia, yet ERP personalization is rarely performed when creating patient-specific computer models of the atria to inform clinical decision-making. This study evaluates the impact of integrating clinical ERP measurements into personalized in silico models on arrhythmia vulnerability.

METHODS AND RESULTS

Clinical ERP measurements were obtained in seven patients from multiple locations in the atria. Atrial geometries from the electroanatomical mapping system were used to generate personalized anatomical atrial models. The Courtemanche M. et al. cellular model was adjusted to reproduce patient-specific ERP. Four modeling approaches were compared: homogeneous (A), heterogeneous (B), regional (C), and continuous (D) ERP distributions. Non-personalized approaches (A and B) were based on literature data, while personalized approaches (C and D) were based on patient measurements. Modeling effects were assessed on arrhythmia vulnerability and tachycardia cycle length, with sensitivity analysis on ERP measurement uncertainty. Mean vulnerability was 3.4 ± 4.0%, 7.7 ± 3.4%, 9.0 ± 5.1%, and 7.0 ± 3.6% for scenarios A-D, respectively. Mean tachycardia cycle length was 167.1 ± 12.6 ms, 158.4 ± 27.5 ms, 265.2 ± 39.9 ms, and 285.9 ± 77.3 ms for scenarios A-D, respectively. Incorporating perturbations to the measured ERP in the range of 2, 5, 10, 20, and 50 ms changed the vulnerability of the model to 5.8 ± 2.7%, 6.1 ± 3.5%, 6.9 ± 3.7%, 5.2 ± 3.5%, and 9.7 ± 10.0%, respectively.

CONCLUSION

Increased ERP dispersion had a greater effect on re-entry dynamics than on vulnerability. Inducibility was higher in personalized scenarios compared with scenarios with uniformly reduced ERP; however, this effect was reversed when incorporating fibrosis informed by low-voltage areas. Effective refractory period measurement uncertainty up to 20 ms slightly influenced vulnerability. Electrophysiological personalization of atrial in silico models appears essential and requires confirmation in larger cohorts.

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.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

Unsupervised Classification of Atrial Electrograms for Electroanatomic Mapping of Human Persistent Atrial Fibrillation

OBJECTIVE

Ablation treatment for persistent atrial fibrillation (persAF) remains challenging due to the absence of a 'ground truth' for atrial substrate characterization and the presence of multiple mechanisms driving the arrhythmia. We implemented an unsupervised classification to identify clusters of atrial electrograms (AEGs) with similar patterns, which were then validated by AEG-derived markers.

METHODS

956 bipolar AEGs were collected from 11 persAF patients. CARTO variables (Biosense Webster; ICL, ACI and SCI) were used to create a 3D space, and subsequently used to perform an unsupervised classification with k-means. The characteristics of the identified groups were investigated using nine AEG-derived markers: sample entropy (SampEn), dominant frequency, organization index (OI), determinism, laminarity, recurrence rate (RR), peak-to-peak (PP) amplitude, cycle length (CL), and wave similarity (WS).

RESULTS

Five AEG classes with distinct characteristics were identified (F = 582, P<0.0001). The presence of fractionation increased from class 1 to 5, as reflected by the nine markers. Class 1 (25%) included organized AEGs with high WS, determinism, laminarity, and RR, and low SampEn. Class 5 (20%) comprised fractionated AEGs with in low WS, OI, determinism, laminarity, and RR, and in high SampEn. Classes 2 (12%), 3 (13%) and 4 (30%) suggested different degrees of AEG organization.

CONCLUSIONS

Our results expand and reinterpret the criteria used for automated AEG classification. The nine markers highlighted electrophysiological differences among the five classes found by the k-means, which could provide a more complete characterization of persAF substrate during ablation target identification in future clinical studies.

Epicardial Fibrosis Explains Increased Endo-Epicardial Dissociation and Epicardial Breakthroughs in Human Atrial Fibrillation

Background

Atrial fibrillation (AF) is accompanied by progressive epicardial fibrosis, dissociation of electrical activity between the epicardial layer and the endocardial bundle network, and transmural conduction (breakthroughs). However, causal relationships between these phenomena have not been demonstrated yet. Our goal was to test the hypothesis that epicardial fibrosis suffices to increase endo–epicardial dissociation (EED) and breakthroughs (BT) during AF.

Methods

We simulated the effect of fibrosis in the epicardial layer on EED and BT in a detailed, high-resolution, three-dimensional model of the human atria with realistic electrophysiology. The model results were compared with simultaneous endo–epicardial mapping in human atria. The model geometry, specifically built for this study, was based on MR images and histo-anatomical studies. Clinical data were obtained in four patients with longstanding persistent AF (persAF) and three patients without a history of AF.

Results

The AF cycle length (AFCL), conduction velocity (CV), and EED were comparable in the mapping studies and the simulations. EED increased from 24.1 ± 3.4 to 56.58 ± 6.2% (p < 0.05), and number of BTs per cycle from 0.89 ± 0.55 to 6.74 ± 2.11% (p < 0.05), in different degrees of fibrosis in the epicardial layer. In both mapping data and simulations, EED correlated with prevalence of BTs. Fibrosis also increased the number of fibrillation waves per cycle in the model.

Conclusion

A realistic 3D computer model of AF in which epicardial fibrosis was increased, in the absence of other pathological changes, showed increases in EED and epicardial BT comparable to those in longstanding persAF. Thus, epicardial fibrosis can explain both phenomena.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia

In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.

Universal atrial coordinates applied to visualisation, registration and construction of patient specific meshes

•

We introduce a coordinate system for the atria based on anatomical landmarks.

•

We construct the coordinates from solutions to Laplace’s equation.

•

We demonstrate the mapping of both scalar and vector data between different atria.

•

The coordinate system was used for registration and 2D visualisation of multimodal data.

•

Patient specific meshes with atrial structures and fibre direction were constructed using just five landmark points.

Integrating spatial information about atrial physiology and anatomy in a single patient from multimodal datasets, as well as generalizing these data across patients, requires a common coordinate system. In the atria, this is challenging due to the complexity and variability of the anatomy. We aimed to develop and validate a Universal Atrial Coordinate (UAC) system for the following applications: combination and assessment of multimodal data; comparison of spatial data across patients; 2D visualization; and construction of patient specific geometries to test mechanistic hypotheses. Left and right atrial LGE-MRI data were segmented and meshed. Two coordinates were calculated for each atrium by solving Laplace’s equation, with boundary conditions assigned using five landmark points. The coordinate system was used to map spatial information between atrial meshes, including scalar fields measured using different mapping modalities, and atrial anatomic structures and fibre directions from a reference geometry. Average error in point transfer from a source mesh to a destination mesh and back again was less than 0.1 mm for the left atrium and 0.02 mm for the right atrium. Patient specific meshes were constructed using the coordinate system and phase singularity density maps from arrhythmia simulations were visualised in 2D. In conclusion, we have developed a universal atrial coordinate system allowing automatic registration of imaging and electroanatomic mapping data, 2D visualisation, and patient specific model creation.

Simulated P wave morphology in the presence of endo-epicardial activation delay

Aims

Evidences of asynchrony between epicardial and endocardial activation in the atrial wall have been reported. We used a computer model of the atria and torso to investigate the consequences of such activation delay on P wave morphology, while controlling for P wave duration.

Methods and results

We created 390 models of the atria based on the same geometry. These models differed by atrial wall thickness (from 2 to 3 mm), transmural coupling, and tissue conductivity in the endocardial and epicardial layers. Among them, 18 were in baseline, 186 had slower conduction in the epicardium layer and 186 in the endocardial layer. Conduction properties were adjusted in such a way that total activation time was the same in all models. P waves on a 16-lead system were simulated during sinus rhythm. Activation maps were similar in all cases. Endo-epicardial delay varied between -5.5 and 5.5 ms vs. 0 ± 0.5 ms in baseline. All P waves had the same duration but variability in their morphology was observed. With slower epicardial conduction, P wave amplitude was reduced by an average of 20% on leads V3-V5 and P wave area decreased by 50% on leads V1-V2 and by 40% on lead V3. Reversed, lower magnitude effects were observed with slower endocardial conduction.

Conclusion

An endo-epicardial delay of a few milliseconds is sufficient to significantly alter P wave morphology, even if the activation map remains the same.

Characterization of human persistent atrial fibrillation electrograms using recurrence quantification analysis

Atrial fibrillation (AF) is regarded as a complex arrhythmia, with one or more co-existing mechanisms, resulting in an intricate structure of atrial activations. Fractionated atrial electrograms (AEGs) were thought to represent arrhythmogenic tissue and hence have been suggested as targets for radiofrequency ablation. However, current methods for ablation target identification have resulted in suboptimal outcomes for persistent AF (persAF) treatment, possibly due to the complex spatiotemporal dynamics of these mechanisms. In the present work, we sought to characterize the dynamics of atrial tissue activations from AEGs collected during persAF using recurrence plots (RPs) and recurrence quantification analysis (RQA). 797 bipolar AEGs were collected from 18 persAF patients undergoing pulmonary vein isolation (PVI). Automated AEG classification (normal vs. fractionated) was performed using the CARTO criteria (Biosense Webster). For each AEG, RPs were evaluated in a phase space estimated following Takens' theorem. Seven RQA variables were obtained from the RPs: recurrence rate; determinism; average diagonal line length; Shannon entropy of diagonal length distribution; laminarity; trapping time; and Shannon entropy of vertical length distribution. The results show that the RQA variables were significantly affected by PVI, and that the variables were effective in discriminating normal vs. fractionated AEGs. Additionally, diagonal structures associated with deterministic behavior were still present in the RPs from fractionated AEGs, leading to a high residual determinism, which could be related to unstable periodic orbits and suggesting a possible chaotic behavior. Therefore, these results contribute to a nonlinear perspective of the spatiotemporal dynamics of persAF.

A work flow to build and validate patient specific left atrium electrophysiology models from catheter measurements

Biophysical models of the atrium provide a physically constrained framework for describing the current state of an atrium and allow predictions of how that atrium will respond to therapy. We propose a work flow to simulate patient specific electrophysiological heterogeneity from clinical data and validate the resulting biophysical models. In 7 patients, we recorded the atrial anatomy with an electroanatomical mapping system (St Jude Velocity); we then applied an S1-S2 electrical stimulation protocol from the coronary sinus (CS) and the high right atrium (HRA) whilst recording the activation patterns using a PentaRay catheter with 10 bipolar electrodes at 12 ± 2 sites across the atrium. Using only the activation times measured with a PentaRay catheter and caused by a stimulus applied in the CS with a remote catheter we fitted the four parameters for a modified Mitchell-Schaeffer model and the tissue conductivity to the recorded local conduction velocity restitution curve and estimated local effective refractory period. Model parameters were then interpolated across each atrium. The fitted model recapitulated the S1-S2 activation times for CS pacing giving a correlation ranging between 0.81 and 0.98. The model was validated by comparing simulated activations times with the independently recorded HRA pacing S1-S2 activation times, giving a correlation ranging between 0.65 and 0.96. The resulting work flow provides the first validated cohort of models that capture clinically measured patient specific electrophysiological heterogeneity.

Efficient parameterization of cardiac action potential models using a genetic algorithm

Finding appropriate values for parameters in mathematical models of cardiac cells is a challenging task. Here, we show that it is possible to obtain good parameterizations in as little as 30-40 s when as many as 27 parameters are fit simultaneously using a genetic algorithm and two flexible phenomenological models of cardiac action potentials. We demonstrate how our implementation works by considering cases of "model recovery" in which we attempt to find parameter values that match model-derived action potential data from several cycle lengths. We assess performance by evaluating the parameter values obtained, action potentials at fit and non-fit cycle lengths, and bifurcation plots for fidelity to the truth as well as consistency across different runs of the algorithm. We also fit the models to action potentials recorded experimentally using microelectrodes and analyze performance. We find that our implementation can efficiently obtain model parameterizations that are in good agreement with the dynamics exhibited by the underlying systems that are included in the fitting process. However, the parameter values obtained in good parameterizations can exhibit a significant amount of variability, raising issues of parameter identifiability and sensitivity. Along similar lines, we also find that the two models differ in terms of the ease of obtaining parameterizations that reproduce model dynamics accurately, most likely reflecting different levels of parameter identifiability for the two models.

Atrial Electrogram Fractionation Distribution before and after Pulmonary Vein Isolation in Human Persistent Atrial Fibrillation-A Retrospective Multivariate Statistical Analysis

Purpose: Complex fractionated atrial electrograms (CFAE)-guided ablation after pulmonary vein isolation (PVI) has been used for persistent atrial fibrillation (persAF) therapy. This strategy has shown suboptimal outcomes due to, among other factors, undetected changes in the atrial tissue following PVI. In the present work, we investigate CFAE distribution before and after PVI in patients with persAF using a multivariate statistical model.

Methods: 207 pairs of atrial electrograms (AEGs) were collected before and after PVI respectively, from corresponding LA regions in 18 persAF patients. Twelve attributes were measured from the AEGs, before and after PVI. Statistical models based on multivariate analysis of variance (MANOVA) and linear discriminant analysis (LDA) have been used to characterize the atrial regions and AEGs.

Results: PVI significantly reduced CFAEs in the LA (70 vs. 40%; P < 0.0001). Four types of LA regions were identified, based on the AEGs characteristics: (i) fractionated before PVI that remained fractionated after PVI (31% of the collected points); (ii) fractionated that converted to normal (39%); (iii) normal prior to PVI that became fractionated (9%) and; (iv) normal that remained normal (21%). Individually, the attributes failed to distinguish these LA regions, but multivariate statistical models were effective in their discrimination (P < 0.0001).

Conclusion: Our results have unveiled that there are LA regions resistant to PVI, while others are affected by it. Although, traditional methods were unable to identify these different regions, the proposed multivariate statistical model discriminated LA regions resistant to PVI from those affected by it without prior ablation information.

Three-dimensional Integrated Functional, Structural, and Computational Mapping to Define the Structural "Fingerprints" of Heart-Specific Atrial Fibrillation Drivers in Human Heart Ex Vivo

BACKGROUND

Structural remodeling of human atria plays a key role in sustaining atrial fibrillation (AF), but insufficient quantitative analysis of human atrial structure impedes the treatment of AF. We aimed to develop a novel 3-dimensional (3D) structural and computational simulation analysis tool that could reveal the structural contributors to human reentrant AF drivers.

METHODS AND RESULTS

High-resolution panoramic epicardial optical mapping of the coronary-perfused explanted intact human atria (63-year-old woman, chronic hypertension, heart weight 608 g) was conducted during sinus rhythm and sustained AF maintained by spatially stable reentrant AF drivers in the left and right atrium. The whole atria (107×61×85 mm³) were then imaged with contrast-enhancement MRI (9.4 T, 180×180×360-μm³ resolution). The entire 3D human atria were analyzed for wall thickness (0.4-11.7 mm), myofiber orientations, and transmural fibrosis (36.9% subendocardium; 14.2% midwall; 3.4% subepicardium). The 3D computational analysis revealed that a specific combination of wall thickness and fibrosis ranges were primarily present in the optically defined AF driver regions versus nondriver tissue. Finally, a 3D human heart-specific atrial computer model was developed by integrating 3D structural and functional mapping data to test AF induction, maintenance, and ablation strategies. This 3D model reproduced the optically defined reentrant AF drivers, which were uninducible when fibrosis and myofiber anisotropy were removed from the model.

CONCLUSIONS

Our novel 3D computational high-resolution framework may be used to quantitatively analyze structural substrates, such as wall thickness, myofiber orientation, and fibrosis, underlying localized AF drivers, and aid the development of new patient-specific treatments.

Minimizing discordances in automated classification of fractionated electrograms in human persistent atrial fibrillation

Ablation of persistent atrial fibrillation (persAF) targeting complex fractionated atrial electrograms (CFAEs) detected by automated algorithms has produced conflicting outcomes in previous electrophysiological studies. We hypothesize that the differences in these algorithms could lead to discordant CFAE classifications by the available mapping systems, giving rise to potential disparities in CFAE-guided ablation. This study reports the results of a head-to-head comparison of CFAE detection performed by NavX (St. Jude Medical) versus CARTO (Biosense Webster) on the same bipolar electrogram data (797 electrograms) from 18 persAF patients. We propose revised thresholds for both primary and complementary indices to minimize the differences in CFAE classification performed by either system. Using the default thresholds [NavX: CFE-Mean ≤ 120 ms; CARTO: ICL ≥ 7], NavX classified 70 % of the electrograms as CFAEs, while CARTO detected 36 % (Cohen’s kappa κ ≈ 0.3, P < 0.0001). Using revised thresholds found using receiver operating characteristic curves [NavX: CFE-Mean ≤ 84 ms, CFE-SD ≤ 47 ms; CARTO: ICL ≥ 4, ACI ≤ 82 ms, SCI ≤ 58 ms], NavX classified 45 %, while CARTO detected 42 % (κ ≈ 0.5, P < 0.0001). Our results show that CFAE target identification is dependent on the system and thresholds used by the electrophysiological study. The thresholds found in this work counterbalance the differences in automated CFAE classification performed by each system. This could facilitate comparisons of CFAE ablation outcomes guided by either NavX or CARTO in future works.

Optimization of catheter ablation of atrial fibrillation: insights gained from clinically-derived computer models

Atrial fibrillation (AF) is the most common heart rhythm disturbance, and its treatment is an increasing economic burden on the health care system. Despite recent intense clinical, experimental and basic research activity, the treatment of AF with current antiarrhythmic drugs and catheter/surgical therapies remains limited. Radiofrequency catheter ablation (RFCA) is widely used to treat patients with AF. Current clinical ablation strategies are largely based on atrial anatomy and/or substrate detected using different approaches, and they vary from one clinical center to another. The nature of clinical ablation leads to ambiguity regarding the optimal patient personalization of the therapy partly due to the fact that each empirical configuration of ablation lines made in a patient is irreversible during one ablation procedure. To investigate optimized ablation lesion line sets, in silico experimentation is an ideal solution. 3D computer models give us a unique advantage to plan and assess the effectiveness of different ablation strategies before and during RFCA. Reliability of in silico assessment is ensured by inclusion of accurate 3D atrial geometry, realistic fiber orientation, accurate fibrosis distribution and cellular kinetics; however, most of this detailed information in the current computer models is extrapolated from animal models and not from the human heart. The predictive power of computer models will increase as they are validated with human experimental and clinical data. To make the most from a computer model, one needs to develop 3D computer models based on the same functionally and structurally mapped intact human atria with high spatial resolution. The purpose of this review paper is to summarize recent developments in clinically-derived computer models and the clinical insights they provide for catheter ablation.

Three-dimensional cardiac computational modelling: methods, features and applications

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

Big heart data: advancing health informatics through data sharing in cardiovascular imaging

The burden of heart disease is rapidly worsening due to the increasing prevalence of obesity and diabetes. Data sharing and open database resources for heart health informatics are important for advancing our understanding of cardiovascular function, disease progression and therapeutics. Data sharing enables valuable information, often obtained at considerable expense and effort, to be reused beyond the specific objectives of the original study. Many government funding agencies and journal publishers are requiring data reuse, and are providing mechanisms for data curation and archival. Tools and infrastructure are available to archive anonymous data from a wide range of studies, from descriptive epidemiological data to gigabytes of imaging data. Meta-analyses can be performed to combine raw data from disparate studies to obtain unique comparisons or to enhance statistical power. Open benchmark datasets are invaluable for validating data analysis algorithms and objectively comparing results. This review provides a rationale for increased data sharing and surveys recent progress in the cardiovascular domain. We also highlight the potential of recent large cardiovascular epidemiological studies enabling collaborative efforts to facilitate data sharing, algorithms benchmarking, disease modeling and statistical atlases.

Heterogeneous and anisotropic integrative model of pulmonary veins: computational study of arrhythmogenic substrate for atrial fibrillation

Mechanisms underlying the genesis of re-entrant substrate for the most common cardiac arrhythmia, atrial fibrillation (AF), are not well understood. In this study, we develop a multi-scale three-dimensional computational model that integrates cellular electrophysiology of the left atrium (LA) and pulmonary veins (PVs) with the respective tissue geometry and fibre orientation. The latter is reconstructed in unique detail from high-resolution (approx. 70 μm) contrast micro-computed tomography data. The model is used to explore the mechanisms of re-entry initiation and sustenance in the PV region, regarded as the primary source of high-frequency electrical activity in AF. Simulations of the three-dimensional model demonstrate that an initial break-down of normal electrical excitation wave-fronts can be caused by the electrical heterogeneity between the PVs and LA. High tissue anisotropy is then responsible for the slow conduction and generation of a re-entrant circuit near the PVs. Evidence of such circuits has been seen clinically in AF patients. Our computational study suggests that primarily the combination of electrical heterogeneity and conduction anisotropy between the PVs and LA tissues leads to the generation of a high-frequency (approx. 10 Hz) re-entrant source near the PV sleeves, thus providing new insights into the arrhythmogenic mechanisms of excitation waves underlying AF.