Identifying the onset of the QRS complex is an important step for localizing the site of origin (SOO) of premature ventricular complexes (PVCs) and the exit site of Ventricular Tachycardia (VT). However, identifying the QRS onset is challenging due to signal noise, baseline wander, motion artifact, and muscle artifact. Furthermore, in VT, QRS onset detection is especially difficult due to the overlap with repolarization from the prior beat. In this study, 7706 captured bipolar pacing beats (Stim-QRS < 40 ms) pooled from 384 anatomically widely dispersed pacing sites of 15 patients were used for an attention-based Swin-Unet neural network. We also utilized a self-supervised pretraining technique using 88253 unannotated ECG records. The algorithm correctly identified most of the onsets for ECG signals with bipolar pacing-site ECG dataset, achieving a sensitivity of 0.958 and a 1.924 ± 4.275 milliseconds prediction error. Our algorithm also achieved a prediction error of 1.518 ± 8.702 milliseconds for the QT Database (QTDB), and a prediction error of 1.333 ± 7.575 milliseconds for the Lobachevsky University Electrocardiography Database (LUDB) public datasets. We also achieved high inter-dataset performance, which supports the practical performance of the method, with a sensitivity of 0.927 for QTDB and a sensitivity of 0.981 for LUDB. The AI model achieves accurate onset detection in paced ECGs with spike-removed inputs, providing a controlled, high-fidelity training setting for future efforts in generalizing to VT ECGs. The use of self-supervised pretraining further improves the detector's accuracy, showcasing the applicability of the approach and using unannotated ECG signals for downstream tasks.

Electrocardiogram (ECG) recordings are affected by the heart's three-dimensional orientation within the thorax, i.e., the anatomical axis. Various cardiac conditions can cause the anatomical axis to shift and/or alter the pattern of electrical activation, leading to changes in the electrical axis. Nevertheless, there remains a lack of a formal, population-level study of the interplay between the cardiac anatomical and electrical axes and the factors that affect them. In this context, this study aimed to: (1) propose standardised definitions for the cardiac anatomical and electrical axes, (2) characterise their population-wide interplay in healthy conditions, (3) evaluate the impact of hypertension on their distribution and (4) identify associations with phenotypical and disease characteristics. Using cardiac magnetic resonance images and 12-lead ECGs from ~39,000 UK Biobank participants, patient-specific, paired biventricular geometries and vectorcardiograms were constructed. Five anatomical and four electrical axis definitions were computed, with the optimal pair of definitions selected based on their mutual alignment in 3D space within 28,000 healthy subjects. Accordingly, the anatomical axis was defined as the vector from the apex to the spatial centre of the four valves, and the electrical axis as the direction of the maximum QRS dipole. Mean angular separation in 3D, [Formula: see text], was 145.0° ± 16.8° in the healthy cohort. The electrical axes exhibited a much larger variability, and strong evidence of anatomical-electrical coupling was identified. Increasing BMI notably affected the anatomical axis, rotating the heart more horizontally-a pattern mirrored by the electrical axis. Both axes were also significantly influenced by sex and, to a lesser extent, age. The axes were then studied in the sub-cohort of ~3,500 UK BioBank participants with primary hypertension, where a similar rotational pattern as that with increasing BMI was revealed. Finally, phenome-wide association studies in the 39,000 participants reveal associations between the axes angular metrics and phenotypes signalling an increased afterload, and an association to hypertension among other clinical conditions. These findings underscore the complex anatomical-electrical interplay and highlight the potential of cardiac axes biomarkers for an improved clinical ECG interpretation and disease characterisation.

Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.

[This corrects the article DOI: 10.1371/journal.pcbi.1008851.].

Cardiac anatomy plays a crucial role in determining cardiac function. However, there is a poor understanding of how specific and localised anatomical changes affect different cardiac functional outputs. In this work, we test the hypothesis that in a statistical shape model (SSM), the modes that are most relevant for describing anatomy are also most important for determining the output of cardiac electromechanics simulations. We made patient-specific four-chamber heart meshes (n = 20) from cardiac CT images in asymptomatic subjects and created a SSM from 19 cases. Nine modes captured 90% of the anatomical variation in the SSM. Functional simulation outputs correlated best with modes 2, 3 and 9 on average (R = 0.49 ± 0.17, 0.37 ± 0.23 and 0.34 ± 0.17 respectively). We performed a global sensitivity analysis to identify the different modes responsible for different simulated electrical and mechanical measures of cardiac function. Modes 2 and 9 were the most important for determining simulated left ventricular mechanics and pressure-derived phenotypes. Mode 2 explained 28.56 ± 16.48% and 25.5 ± 20.85, and mode 9 explained 12.1 ± 8.74% and 13.54 ± 16.91% of the variances of mechanics and pressure-derived phenotypes, respectively. Electrophysiological biomarkers were explained by the interaction of 3 ± 1 modes. In the healthy adult human heart, shape modes that explain large portions of anatomical variance do not explain equivalent levels of electromechanical functional variation. As a result, in cardiac models, representing patient anatomy using a limited number of modes of anatomical variation can cause a loss in accuracy of simulated electromechanical function.

Atrial fibrillation (AF) is a supraventricular tachyarrhythmia characterized by complete absence of coordinated atrial contraction and is associated with an increased morbidity and mortality. Personalized computational modeling provides a novel framework for integrating and interpreting the role of atrial electrophysiology (EP) including the underlying anatomy and microstructure in the development and sustenance of AF. Coronary computed tomography angiography data were segmented using a statistics-based approach and the smoothed voxel representations were discretized into high-resolution tetrahedral finite element (FE) meshes. To estimate the complex left atrial myofiber architecture, individual fiber fields were generated according to morphological data on the endo- and epicardial surfaces based on local solutions of Laplace's equation and transmurally interpolated to tetrahedral elements. The influence of variable transmural microstructures was quantified through EP simulations on 3 patients using 5 different fiber interpolation functions. Personalized geometrical models included the heterogeneous thickness distribution of the left atrial myocardium and subsequent discretization led to high-fidelity tetrahedral FE meshes. The novel algorithm for automated incorporation of the left atrial fiber architecture provided a realistic estimate of the atrial microstructure and was able to qualitatively capture all important fiber bundles. Consistent maximum local activation times were predicted in EP simulations using individual transmural fiber interpolation functions for each patient suggesting a negligible effect of the transmural myofiber architecture on EP. The established modeling pipeline provides a robust framework for the rapid development of personalized model cohorts accounting for detailed anatomy and microstructure and facilitates simulations of atrial EP.

Knowledge of atrial wall thickness (AWT) has the potential to provide important information for patient stratification and the planning of interventions in atrial arrhythmias. To date, information about AWT has only been acquired in post-mortem or poor-contrast computed tomography (CT) studies, providing limited coverage and highly variable estimates of AWT. We present a novel contrast agent-free MRI sequence for imaging AWT and use it to create personalized AWT maps and a biatrial atlas. A novel black-blood phase-sensitive inversion recovery protocol was used to image ten volunteers and, as proof of concept, two atrial fibrillation patients. Both atria were manually segmented to create subject-specific AWT maps using an average of nearest neighbors approach. These were then registered non-linearly to generate an AWT atlas. AWT was 2.4 ± 0.7 and 2.7 ± 0.7 mm in the left and right atria, respectively, in good agreement with post-mortem and CT data, where available. AWT was 2.6 ± 0.7 mm in the left atrium of a patient without structural heart disease, similar to that of volunteers. In a patient with structural heart disease, the AWT was increased to 3.1 ± 1.3 mm. We successfully designed an MRI protocol to non-invasively measure AWT and create the first whole-atria AWT atlas. The atlas can be used as a reference to study alterations in thickness caused by atrial pathology. The protocol can be used to acquire personalized AWT maps in a clinical setting and assist in the treatment of atrial arrhythmias.