Investigating the Influence of Parameters on Electrogram Morphology Using Gaussian Process Emulators

Manisha Sahota1, Fernando Campos2, Shuang Qian1, John Whitaker1, Steven E Williams3, Matthijs Cluitmans4, Steven Niederer5, Martin Bishop1
1King's College London, 2UNESP, 3University of Edinburgh, 4Maastricht University, 5Imperial College London


Abstract

Background: Intracardiac electrograms (EGMs) are commonly used to assess cardiac arrhythmias to gain information about structural alterations of underlying tissue. Computational models offer a systematic approach to study how parameters such as myocardial tissue thickness influence EGM morphology. However, global sensitivity analyses are challenging to perform due to the high computational cost of such multi-parametric and nonlinear models.

Objective: To evaluate the impact of myocardial tissue thickness, electrode size and conduction slowing on the amplitude and minimum temporal derivative of the unipolar EGM voltage (dV/dtmin), identifying key parameters that influence the EGM morphology.

Methods: A realistic computational model of healthy myocardial tissue and bath was generated with a highly conductive electrode on the endocardial surface. Bidomain simulations were performed with a planar wavefront perpendicular to the electrode. The impact of myocardial tissue thickness, electrode size and conduction slowing via altering fibrosis density, sodium channel conductance (gNa) and intracellular (σi) and extracellular (σe) conductivities on EGM and action potential (AP) morphology were investigated by training Gaussian process emulators (GPEs) to perform a Sobol variance-based global sensitivity analysis (GSA).

Results: The trained emulators achieved a median R2 test score of 0.91 and ISE of 97.5. GSA revealed a strong influence of fibrosis density and σi/σe on EGM amplitude, with reduced fibrosis density and increased σi/σe generating EGMs with larger amplitudes. Both σi/σe and gNa exerted similar effects on dV/dtmin, though gNa more significantly influenced AP morphology, as expected. Myocardial tissue thickness heterogeneity most greatly influenced the EGM amplitude but had minimal association with dV/dtmin. Electrode width had the greatest impact on dV/dtmin, though electrode length and height had negligible influence.

Conclusion: This study demonstrates the efficacy of GPEs as a robust approach to perform a computationally efficient GSA, gaining deeper insights into the link between cell and tissue properties and intracardiac EGMs.