A new algorithm to diagnose atrial ectopic origin from multi lead ECG systems - insights from 3D virtual human atria and torso

Rapid atrial arrhythmias such as atrial fibrillation (AF) predispose to ventricular arrhythmias, sudden cardiac death and stroke. Identifying the origin of atrial ectopic activity from the electrocardiogram (ECG) can help to diagnose the early onset of AF in a cost-effective manner. The complex and rapid atrial electrical activity during AF makes it difficult to obtain detailed information on atrial activation using the standard 12-lead ECG alone. Compared to conventional 12-lead ECG, more detailed ECG lead configurations may provide further information about spatio-temporal dynamics of the body surface potential (BSP) during atrial excitation. We apply a recently developed 3D human atrial model to simulate electrical activity during normal sinus rhythm and ectopic pacing. The atrial model is placed into a newly developed torso model which considers the presence of the lungs, liver and spinal cord. A boundary element method is used to compute the BSP resulting from atrial excitation. Elements of the torso mesh corresponding to the locations of the placement of the electrodes in the standard 12-lead and a more detailed 64-lead ECG configuration were selected. The ectopic focal activity was simulated at various origins across all the different regions of the atria. Simulated BSP maps during normal atrial excitation (i.e. sinoatrial node excitation) were compared to those observed experimentally (obtained from the 64-lead ECG system), showing a strong agreement between the evolution in time of the simulated and experimental data in the P-wave morphology of the ECG and dipole evolution. An algorithm to obtain the location of the stimulus from a 64-lead ECG system was developed. The algorithm presented had a success rate of 93%, meaning that it correctly identified the origin of atrial focus in 75/80 simulations, and involved a general approach relevant to any multi-lead ECG system. This represents a significant improvement over previously developed algorithms

Trigger versus Substrate: Multi-Dimensional Modulation of QT-Prolongation Associated Arrhythmic Dynamics by a hERG Channel Activator

Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca²⁺-dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca²⁺ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on INa. The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval

Effects of Heart Rate and Ventricular Wall Thickness on Non-invasive Mapping: An in silico Study

Background: Non-invasive cardiac mapping—also known as Electrocardiographic imaging (ECGi)—is a novel, painless and relatively economic method to map the electrical activation and repolarization patterns of the heart, providing a valuable tool for early identification and diagnosis of conduction abnormalities and arrhythmias. Moreover, the ability to obtain information on cardiac electrical activity non-invasively using ECGi provides the potential for a priori information to guide invasive surgical procedures, improving success rates, and reducing procedure time. Previous studies have shown the influence of clinical variables, such as heart rate, heart size, endocardial wall, and body composition on surface electrocardiogram (ECG) measurements. The influence of clinical variables on the ECG variability has provided information on cardiovascular control and its abnormalities in various pathologies. However, the effects of such clinical variables on the Body Surface Potential (BSP) and ECGi maps have yet to be systematically investigated. Methods: In this study we investigated the effects of heart size, intracardiac thickness, and heart rate on BSP and ECGi maps using a previously-developed 3D electrophysiologically-detailed ventricles-torso model. The inverse solution was solved using the three different Tikhonov regularization methods. Results: Through comparison of multiple measures of error/accuracy on the ECGi reconstructions, our results showed that using different heart geometries to solve the forward and inverse problems produced a larger estimated focal excitation location. An increase of ~2 mm in the Euclidean distance error was observed for an increase in the heart size. However, the estimation of the location of focal activity was still able to be obtained. Similarly, a Euclidean distance increase was observed when the order of regularization was reduced. For the case of activation maps reconstructed at the same ectopic focus location but different heart rates, an increase in the errors and Euclidean distance was observed when the heart rate was increased. Conclusions: Non-invasive cardiac mapping can still provide useful information about cardiac activation patterns for the cases when a different geometry is used for the inverse problem compared to the one used for the forward solution; rapid pacing rates can induce order-dependent errors in the accuracy of reconstruction

Generalized Toda Theory from Six Dimensions and the Conifold

Recently, a physical derivation of the Alday-Gaiotto-Tachikawa correspondence has been put forward. A crucial role is played by the complex Chern-Simons theory arising in the 3d-3d correspondence, whose boundary modes lead to Toda theory on a Riemann surface. We explore several features of this derivation and subsequently argue that it can be extended to a generalization of the AGT correspondence. The latter involves codimension two defects in six dimensions that wrap the Riemann surface. We use a purely geometrical description of these defects and find that the generalized AGT setup can be modeled in a pole region using generalized conifolds. Furthermore, we argue that the ordinary conifold clarifies several features of the derivation of the original AGT correspondence.Comment: 27+2 pages, 3 figure

ECG Imaging to Detect the Site of Ventricular Ischemia Using Torso Electrodes: A Computational Study

Electrocardiography provides some information useful for ischemic diagnosis. However, more recently there has been substantial growth in the area of ECG imaging, which by solving the inverse problem of electrocardiography aims to produce high-resolution mapping of the electrical and magnetic dynamics of the heart. Most inverse studies use the full resolution of the body surface potential (BSP) to reconstruct the epicardial potentials, however using a limited number of torso electrodes to interpolate the BSP is more clinically relevant and has an important effect on the reconstruction which must be quantified. A circular ischemic lesion on the right ventricle lateral wall 27 mm in radius is reconstructed using three Tikhonov methods along with 6 different electrode configurations ranging from 32 leads to 1,024 leads. The 2nd order Tikhonov solution performed the most accurately (~80% lesion identified) followed by the 1st (~50% lesion identified) and then the 0 order Tikhonov solution performed the worst with a maximum of ~30% lesion identified regardless of how many leads were used. With an increasing number of leads the solution produces less error, and the error becomes more localised around the lesion for all three regularisation methods. In noisy conditions, the relative performance gap of the 1st and 2nd order Tikhonov solutions was reduced, and determining an accurate regularisation parameter became relatively more difficult. Lesions located on the left ventricle walls were also able to be identified but comparatively to the right ventricle lateral wall performed marginally worse with lesions located on the interventricular septum being able to be indicated by the reconstructions but not successfully identified against the error. The quality of reconstruction was found to decrease as the lesion radius decreased, with a lesion radius of <20 mm becoming difficult to correctly identify against the error even when using >512 torso electrodes

Novel non-invasive algorithm to identify the origins of re-entry and ectopic foci in the atria from 64-lead ECGs: A computational study.

Atrial tachy-arrhytmias, such as atrial fibrillation (AF), are characterised by irregular electrical activity in the atria, generally associated with erratic excitation underlain by re-entrant scroll waves, fibrillatory conduction of multiple wavelets or rapid focal activity. Epidemiological studies have shown an increase in AF prevalence in the developed world associated with an ageing society, highlighting the need for effective treatment options. Catheter ablation therapy, commonly used in the treatment of AF, requires spatial information on atrial electrical excitation. The standard 12-lead electrocardiogram (ECG) provides a method for non-invasive identification of the presence of arrhythmia, due to irregularity in the ECG signal associated with atrial activation compared to sinus rhythm, but has limitations in providing specific spatial information. There is therefore a pressing need to develop novel methods to identify and locate the origin of arrhythmic excitation. Invasive methods provide direct information on atrial activity, but may induce clinical complications. Non-invasive methods avoid such complications, but their development presents a greater challenge due to the non-direct nature of monitoring. Algorithms based on the ECG signals in multiple leads (e.g. a 64-lead vest) may provide a viable approach. In this study, we used a biophysically detailed model of the human atria and torso to investigate the correlation between the morphology of the ECG signals from a 64-lead vest and the location of the origin of rapid atrial excitation arising from rapid focal activity and/or re-entrant scroll waves. A focus-location algorithm was then constructed from this correlation. The algorithm had success rates of 93% and 76% for correctly identifying the origin of focal and re-entrant excitation with a spatial resolution of 40 mm, respectively. The general approach allows its application to any multi-lead ECG system. This represents a significant extension to our previously developed algorithms to predict the AF origins in association with focal activities

Chronology of the sedimentary processes during the postglacial sea level rise in two estuaries of the Algarve coast, Southern Portugal

Four profiles of estuarine sediments obtained from boreholes drilled in the Algarve, Southern Portugal were studied in order to reconstruct the process of sediment accumulation driven by the postglacial sea level rise. In addition to the sedimentological analysis, the Foraminifera Index of Marine Influence (FIMI) permitted assessment of the nature and organization of sedimentary facies in the BelicheeGuadiana and Gilão-Almargem estuaries. The Beliche- Guadiana CM5 and Almargem G2 profiles accumulated in a sheltered environment, with the former presenting an almost continuous record of the sea level rise since ca 13 000 cal yr BP. The G1 and G3 profiles from the Gilão-Almargem area represent a more discontinuous record of the last 8000 years, which accumulated in the more dynamic environment of an outer estuary. The integration of all radiocarbon ages of dated levels, led to an estimate of sediment accumulation rates. Assuming a constant position of the sediment surface with respect to the tidal range and a negligible compaction of sediment, the sea level rose at the rate of 7 mm yr ^-1 in the period from 13 000 to 7500 cal yr BP. This process slowed down to ca 0.9 mm yr 1 from 7500 cal yr BP until the present. The marked historical change in the rate of sediment accumulation in these estuaries also occurred with the accumulation of organic matter and is, therefore, important data for global biogeochemical models of carbon. The main obstacle to obtain higher temporal resolution of the sedimentary processes was the intense anaerobic respiration of organic matter via sulphate reduction, which did not allow any accumulation of peat and, furthermore, led to erasure of the palaeontological record by acid formed from the subsequent oxidation of sulphides.FORMOSE- Sources and Retention of Organic Matter in the Estuarine Zones, PRAXIS XXI program of Portuguese Science and Technology Foundation and project MEGASIG, INTERREG IIIa program of the European Union