A computer model of normal conduction in the human atria

Three-dimensional impulse propagation in myocardium: arrhythmogenic mechanisms at the tissue level

Impulse propagation in the heart depends on the excitability of individual cardiomyocytes, impulse transmission between adjacent myocytes, and the 3-dimensional arrangement of those cells. Here, we review the role of each of these factors in normal and aberrant cardiac electric activation, with particular emphasis on the effects of 3-dimensional myocyte architecture at the tissue scale. The analysis draws on findings from in vivo and in vitro experiments, as well as biophysically based computer models that have been used to integrate and interpret these experimental data. It indicates that discontinuous arrangement of myocytes and extracellular connective tissue at the tissue scale can give rise to current source-to-sink mismatch, spatiotemporal distribution of refractoriness, and rate-sensitive electric instability, which contribute to the initiation and maintenance of reentrant cardiac arrhythmia. This exacerbates the risk of rhythm disturbance associated with heart disease. We conclude that structure-based, multiscale computer models that incorporate accurate information about local cellular electric activity provide a powerful platform for investigating the basis of reentrant cardiac arrhythmia. However, it is important that these models capture key features of structure and related electric function at the tissue scale.

Theoretical considerations for mapping activation in human cardiac fibrillation

Defining mechanisms for cardiac fibrillation is challenging because, in contrast to other arrhythmias, fibrillation exhibits complex non-repeatability in spatiotemporal activation but paradoxically exhibits conserved spatial gradients in rate, dominant frequency, and electrical propagation. Unlike animal models, in which fibrillation can be mapped at high spatial and temporal resolution using optical dyes or arrays of contact electrodes, mapping of cardiac fibrillation in patients is constrained practically to lower resolutions or smaller fields-of-view. In many animal models, atrial fibrillation is maintained by localized electrical rotors and focal sources. However, until recently, few studies had revealed localized sources in human fibrillation, so that the impact of mapping constraints on the ability to identify rotors or focal sources in humans was not described. Here, we determine the minimum spatial and temporal resolutions theoretically required to detect rigidly rotating spiral waves and focal sources, then extend these requirements for spiral waves in computer simulations. Finally, we apply our results to clinical data acquired during human atrial fibrillation using a novel technique termed focal impulse and rotor mapping (FIRM). Our results provide theoretical justification and clinical demonstration that FIRM meets the spatio-temporal resolution requirements to reliably identify rotors and focal sources for human atrial fibrillation.

Dominant frequency and organization index maps in a realistic three-dimensional computational model of atrial fibrillation

AIMS

To study, using simulation, the spectral characteristics of different patterns of atrial fibrillation (AF) at high spatial resolution. Dominant frequency (DF) and organization index (OI) maps have been used to approximate the location of the focal source of high frequency during AF events.

METHODS AND RESULTS

A realistic three-dimensional model of the human atria that includes fibre orientation, electrophysiological heterogeneity, and anisotropy was implemented. The cellular model was modified to simulate electrical remodelling. More than 43 000 electrograms were calculated on the surface, and were processed to reconstitute the DF and OI maps. Atrial fibrillation episodes were triggered by a source of transitory and of continuous activity (both with a cycle length of 130 ms) in five different locations. The maps obtained during the AF events triggered by transitory foci did not show areas with high DF or OI values. When continuous foci were applied, the DF maps show ample zones with high values in the atrium where the focus was applied; while OI maps display smaller areas with high values, always within the areas of high DF and, in three of five locations, this high-value area was located at the site of focus application and at the nearby area. In the other two locations, the area presenting the highest OI values is small and located at the site of focus application, which allowed its precise localization.

CONCLUSION

Organization index maps provide a better approximation than DF maps for the localization of ectopic sources of high frequency and continuous activity during episodes of simulated AF in remodelled tissue.

Fully automated initiation of simulated episodes of atrial arrhythmias

AIMS

To develop computational tools for automatically initiating a large number of independent episodes of atrial arrhythmias in electro-anatomical computer models of the atria and therefore facilitating the design of in silico experiments.

METHODS AND RESULTS

A biophysical model of the atria was constructed from segmented medical images of the human atria of a patient with atrial fibrillation (AF). A set of 40 initial conditions were generated based on a priori knowledge about wavefront propagation and the number and location of reentries (1-6 randomly distributed over the atrial epicardium). Simulations were run from each of these initial conditions in three substrates representing different forms of AF dynamics (stable rotors; multiple unstable meandering wavelets; and wavelets broken by repolarization heterogeneities). To demonstrate the applicability of the initiation method for testing clinical of therapeutic interventions, the channel I(Kr) was blocked after 2 s of simulation and its effect on the number of functional reentries was documented. The use of pre-computed initial conditions enabled to successfully generate episodes of simulated AF in each substrate. Blockade of I(Kr) channel prolonged action potential duration, resulting in a reduction of the number of functional reentries. In the substrate with unstable spiral waves, the effect was sufficiently large to terminate AF in about two-thirds of the cases. In the two other substrates, the effect was minor.

CONCLUSION

These new simulation tools may help investigate in computer models therapeutic interventions in different substrates in order to identify substrate-specific optimal therapy.

A simplified 3D model of whole heart electrical activity and 12-lead ECG generation

We present a computationally efficient three-dimensional bidomain model of torso-embedded whole heart electrical activity, with spontaneous initiation of activation in the sinoatrial node, incorporating a specialized conduction system with heterogeneous action potential morphologies throughout the heart. The simplified geometry incorporates the whole heart as a volume source, with heart cavities, lungs, and torso as passive volume conductors. We placed four surface electrodes at the limbs of the torso: V_R, V_L, V_F and V_GND and six electrodes on the chest to simulate the Einthoven, Goldberger-augmented and precordial leads of a standard 12-lead system. By placing additional seven electrodes at the appropriate torso positions, we were also able to calculate the vectorcardiogram of the Frank lead system. Themodel was able to simulate realistic electrocardiogram (ECG) morphologies for the 12 standard leads, orthogonal X, Y, and Z leads, as well as the vectorcardiogram under normal and pathological heart states. Thus, simplified and easy replicable 3D cardiac bidomain model offers a compromise between computational load and model complexity and can be used as an investigative tool to adjust cell, tissue, and whole heart properties, such as setting ischemic lesions or regions of myocardial infarction, to readily investigate their effects on whole ECG morphology.

A three-dimensional finite element model of human atrial anatomy: new methods for cubic Hermite meshes with extraordinary vertices

High-order cubic Hermite finite elements have been valuable in modeling cardiac geometry, fiber orientations, biomechanics, and electrophysiology, but their use in solving three-dimensional problems has been limited to ventricular models with simple topologies. Here, we utilized a subdivision surface scheme and derived a generalization of the "local-to-global" derivative mapping scheme of cubic Hermite finite elements to construct bicubic and tricubic Hermite models of the human atria with extraordinary vertices from computed tomography images of a patient with atrial fibrillation. To an accuracy of 0.6 mm, we were able to capture the left atrial geometry with only 142 bicubic Hermite finite elements, and the right atrial geometry with only 90. The left and right atrial bicubic Hermite meshes were G1 continuous everywhere except in the one-neighborhood of extraordinary vertices, where the mean dot products of normals at adjacent elements were 0.928 and 0.925. We also constructed two biatrial tricubic Hermite models and defined fiber orientation fields in agreement with diagrammatic data from the literature using only 42 angle parameters. The meshes all have good quality metrics, uniform element sizes, and elements with aspect ratios near unity, and are shared with the public. These new methods will allow for more compact and efficient patient-specific models of human atrial and whole heart physiology.

Simulation of biatrial conduction via different pathways during sinus rhythm with a detailed human atrial model

In order to better understand biatrial conduction, investigate various conduction pathways, and compare the differences between isotropic and anisotropic conductions in human atria, we present a simulation study of biatrial conduction with known/assumed conduction pathways using a recently developed human atrial model. In addition to known pathways: (1) Bachmann's bundle (BB), (2) limbus of fossa ovalis (LFO), and (3) coronary sinus (CS), we also hypothesize that there exist two fast conduction bundles that connect the crista terminalis (CT), LFO, and CS. Our simulation demonstrates that use of these fast conduction bundles results in a conduction pattern consistent with experimental data. The comparison of isotropic and anisotropoic conductions in the BB case showed that the atrial working muscles had small effect on conduction time and conduction speed, although the conductivities assigned in anisotropic conduction were two to four times higher than the isotropic conduction. In conclusion, we suggest that the hypothesized intercaval bundles play a significant role in the biatrial conduction and that myofiber orientation has larger effects on the conduction system than the atrial working muscles. This study presents readers with new insights into human atrial conduction.

A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship

The most common sustained cardiac arrhythmias in humans are atrial tachyarrhythmias, mainly atrial fibrillation. Areas of complex fractionated atrial electrograms and high dominant frequency have been proposed as critical regions for maintaining atrial fibrillation; however, there is a paucity of data on the relationship between the characteristics of electrograms and the propagation pattern underlying them. In this study, a realistic 3D computer model of the human atria has been developed to investigate this relationship. The model includes a realistic geometry with fiber orientation, anisotropic conductivity and electrophysiological heterogeneity. We simulated different tachyarrhythmic episodes applying both transient and continuous ectopic activity. Electrograms and their dominant frequency and organization index values were calculated over the entire atrial surface. Our simulations show electrograms with simple potentials, with little or no cycle length variations, narrow frequency peaks and high organization index values during stable and regular activity as the observed in atrial flutter, atrial tachycardia (except in areas of conduction block) and in areas closer to ectopic activity during focal atrial fibrillation. By contrast, cycle length variations and polymorphic electrograms with single, double and fragmented potentials were observed in areas of irregular and unstable activity during atrial fibrillation episodes. Our results also show: (1) electrograms with potentials without negative deflection related to spiral or curved wavefronts that pass over the recording point and move away, (2) potentials with a much greater proportion of positive deflection than negative in areas of wave collisions, (3) double potentials related with wave fragmentations or blocking lines and (4) fragmented electrograms associated with pivot points. Our model is the first human atrial model with realistic fiber orientation used to investigate the relationship between different atrial arrhythmic propagation patterns and the electrograms observed at more than 43000 points on the atrial surface.

Personalization of atrial anatomy and electrophysiology as a basis for clinical modeling of radio-frequency ablation of atrial fibrillation

Multiscale cardiac modeling has made great advances over the last decade. Highly detailed atrial models were created and used for the investigation of initiation and perpetuation of atrial fibrillation. The next challenge is the use of personalized atrial models in clinical practice. In this study, a framework of simple and robust tools is presented, which enables the generation and validation of patient-specific anatomical and electrophysiological atrial models. Introduction of rule-based atrial fiber orientation produced a realistic excitation sequence and a better correlation to the measured electrocardiograms. Personalization of the global conduction velocity lead to a precise match of the measured P-wave duration. The use of a virtual cohort of nine patient and volunteer models averaged out possible model-specific errors. Intra-atrial excitation conduction was personalized manually from left atrial local activation time maps. Inclusion of LE-MRI data into the simulations revealed possible gaps in ablation lesions. A fast marching level set approach to compute atrial depolarization was extended to incorporate anisotropy and conduction velocity heterogeneities and reproduced the monodomain solution. The presented chain of tools is an important step towards the use of atrial models for the patient-specific AF diagnosis and ablation therapy planing.

Application of micro-computed tomography with iodine staining to cardiac imaging, segmentation, and computational model development

Micro-computed tomography (micro-CT) has been widely used to generate high-resolution 3-D tissue images from small animals nondestructively, especially for mineralized skeletal tissues. However, its application to the analysis of soft cardiovascular tissues has been limited by poor inter-tissue contrast. Recent ex vivo studies have shown that contrast between muscular and connective tissue in micro-CT images can be enhanced by staining with iodine. In the present study, we apply this novel technique for imaging of cardiovascular structures in canine hearts. We optimize the method to obtain high-resolution X-ray micro-CT images of the canine atria and its distinctive regions-including the Bachmann's bundle, atrioventricular node, pulmonary arteries and veins-with clear inter-tissue contrast. The imaging results are used to reconstruct and segment the detailed 3-D geometry of the atria. Structure tensor analysis shows that the arrangement of atrial fibers can also be characterized using the enhanced micro-CT images, as iodine preferentially accumulates within the muscular fibers rather than in connective tissues. This novel technique can be particularly useful in nondestructive imaging of 3-D cardiac architectures from large animals and humans, due to the combination of relatively high speed ( ~ 1 h/per scan of the large canine heart) and high voxel resolution (36 μm) provided. In summary, contrast micro-CT facilitates fast and nondestructive imaging and segmenting of detailed 3-D cardiovascular geometries, as well as measuring fiber orientation, which are crucial in constructing biophysically detailed computational cardiac models.

An eikonal-diffusion solver and its application to the interpolation and the simulation of reentrant cardiac activations

Electrical propagation of the cardiac impulse in the myocardium can be described by the eikonal-diffusion equation. This equation governs the field of activation times in a domain where conduction properties are specified. This approach has been applied to knowledge-based interpolation of sparse measurements of activation times and to the creation of initial conditions for detailed ionic models of cardiac propagation. This paper presents the mathematical basis, matrix formulation, and compact Matlab implementation of an iterative finite-element solver (triangular meshes) for the eikonal-diffusion equation extended to reentrant activations, which automatically identifies the period of reentry and computes the resulting isochrones. An iterative algorithm is designed to perform Laplacian interpolation of reentrant activation maps to be used as initial estimate for the eikonal-diffusion solver. The performance of the algorithm is analyzed in test-case geometries (ventricular slice and simplified atrial surface model).

Human atrial fibrillation: insights from computational electrophysiological models

Computational electrophysiology has proven useful to investigate the mechanisms of cardiac arrhythmias at various spatial scales, from isolated myocytes to the whole heart. This article reviews how mathematical modeling has aided our understanding of human atrial myocyte electrophysiology to study the contribution of structural and electrical remodeling to human atrial fibrillation. Potential new avenues of investigation and model development are suggested.

An image-based model of the whole human heart with detailed anatomical structure and fiber orientation

Many heart anatomy models have been developed to study the electrophysiological properties of the human heart. However, none of them includes the geometry of the whole human heart. In this study, an anatomically detailed mathematical model of the human heart was firstly reconstructed from the computed tomography images. In the reconstructed model, the atria consisted of atrial muscles, sinoatrial node, crista terminalis, pectinate muscles, Bachmann's bundle, intercaval bundles, and limbus of the fossa ovalis. The atrioventricular junction included the atrioventricular node and atrioventricular ring, and the ventricles had ventricular muscles, His bundle, bundle branches, and Purkinje network. The epicardial and endocardial myofiber orientations of the ventricles and one layer of atrial myofiber orientation were then measured. They were calculated using linear interpolation technique and minimum distance algorithm, respectively. To the best of our knowledge, this is the first anatomically-detailed human heart model with corresponding experimentally measured fibers orientation. In addition, the whole heart excitation propagation was simulated using a monodomain model. The simulated normal activation sequence agreed well with the published experimental findings.

Modeling defibrillation of the heart: approaches and insights

Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices (ICDs), constitutes the most important means of combating sudden cardiac death. While ICD therapy has proved to be efficient and reliable, defibrillation is a traumatic experience. Thus, research on defibrillation mechanisms, particularly aimed at lowering defibrillation voltage, remains an important topic. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart is the most promising approach to achieve this goal. The aim of this paper is to assess the current state-of-the-art in ventricular defibrillation modeling, focusing on both numerical modeling approaches and major insights that have been obtained using defibrillation models, primarily those of realistic ventricular geometry. The paper showcases the contributions that modeling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modeling of the heart (i.e., cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level and has the potential to contribute to the betterment of the clinical practice of defibrillation.

Patient-specific identification of optimal ubiquitous electrocardiogram (U-ECG) placement using a three-dimensional model of cardiac electrophysiology

A bipolar mini-ECG for ubiquitous healthcare (U-ECG) has been introduced, and various studies using the U-ECG device are in progress. Because it uses two electrodes within a small torso surface area, the design of the U-ECG must be suitable for detecting ECG signals. Using a 3-D model of cardiac electrophysiology, we have developed a simulation method for identifying the optimal placement of U-ECG electrodes on the torso surface. We simulated the heart-torso model to obtain a body surface potential map and ECG waveforms, which were compared with the empirical data. Using this model, we determined the optimal placement of the two U-ECG electrodes, spaced 5 cm apart, for detecting the P, R, and T waves. The ECG data, obtained using the optimal U-ECG placement for a specific wave, showed a clear shape for the target wave, but equivocal shapes for the other waves. The present study provides an efficient simulation method to identify the optimal attachment position and direction of the U-ECG electrodes on the surface of the torso.

Modeling atrial arrhythmias: impact on clinical diagnosis and therapies

Atrial arrhythmias are the most frequent sustained rhythm disorders in humans and often lead to severe complications such as heart failure and stroke. Despite the important insights provided by animal models into the mechanisms of atrial arrhythmias, direct translation of experimental findings to new therapies in patients has not been straightforward. With the advances in computer technology, large-scale electroanatomical computer models of the atria that integrate information from the molecular to organ scale have reached a level of sophistication that they can be used to interpret the outcome of experimental and clinical studies and aid in the rational design of therapies. This paper reviews the state-of-the-art of computer models of the electrical dynamics of the atria and discusses the evolving role of simulation in assisting the clinical diagnosis and treatment of atrial arrhythmias.

Computational modeling of the human atrial anatomy and electrophysiology

This review article gives a comprehensive survey of the progress made in computational modeling of the human atria during the last 10 years. Modeling the anatomy has emerged from simple "peanut"-like structures to very detailed models including atrial wall and fiber direction. Electrophysiological models started with just two cellular models in 1998. Today, five models exist considering e.g. details of intracellular compartments and atrial heterogeneity. On the pathological side, modeling atrial remodeling and fibrotic tissue are the other important aspects. The bridge to data that are measured in the catheter laboratory and on the body surface (ECG) is under construction. Every measurement can be used either for model personalization or for validation. Potential clinical applications are briefly outlined and future research perspectives are suggested.

Clinical mapping approach to diagnose electrical rotors and focal impulse sources for human atrial fibrillation

INTRODUCTION

The perpetuating mechanisms for human atrial fibrillation (AF) remain undefined. Localized rotors and focal beat sources may sustain AF in elegant animal models, but there has been no direct evidence for localized sources in human AF using traditional methods. We developed a clinical computational mapping approach, guided by human atrial tissue physiology, to reveal sources of human AF.

METHODS AND RESULTS

In 49 AF patients referred for ablation (62 ± 9 years; 30 persistent), we defined repolarization dynamics using monophasic action potentials (MAPs) and recorded AF activation from 64-pole basket catheters in left atrium and, in n = 20 patients, in both atria. Careful positioning of basket catheters was required for optimal mapping. AF electrograms at 64-128 electrodes were combined with repolarization and conduction dynamics to construct spatiotemporal AF maps. We observed sustained sources in 47/49 patients, in the form of electrical rotors (n = 57) and focal beats (n = 11) that controlled local atrial activation with peripheral wavebreak (fibrillatory conduction). Patients with persistent AF had more sources than those with paroxysmal AF (2.1 ± 1.0 vs 1.5 ± 0.8, P = 0.02), related to shorter cycle length (163 ± 19 milliseconds vs 187 ± 25 milliseconds, P < 0.001). Approximately one-quarter of sources lay in the right atrium.

CONCLUSIONS

Physiologically guided computational mapping revealed sustained electrical rotors and repetitive focal beats during human AF for the first time. These localized sources were present in 96% of AF patients, and controlled AF activity. These results provide novel mechanistic insights into human AF and lay the foundation for mechanistically tailored approaches to AF ablation.

Effects of electrical and structural remodeling on atrial fibrillation maintenance: a simulation study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength--electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely.

Estimating the time scale and anatomical location of atrial fibrillation spontaneous termination in a biophysical model

Due to their transient nature, spontaneous terminations of atrial fibrillation (AF) are difficult to investigate. Apparently, confounding experimental findings about the time scale of this phenomenon have been reported, with values ranging from 1 s to 1 min. We propose a biophysical modeling approach to study the mechanisms of spontaneous termination in two models of AF with different levels of dynamical complexity. 8 s preceding spontaneous terminations were studied and the evolution of cycle length and wavefront propagation were documented to assess the time scale and anatomical location of the phenomenon. Results suggest that termination mechanisms are dependent on the underlying complexity of AF. During simulated AF of low complexity, the total process of spontaneous termination lasted 3,200 ms and was triggered in the left atrium 800 ms earlier than in the right atrium. The last fibrillatory activity was observed more often in the right atrium. These asymmetric termination mechanisms in both time and space were not observed during spontaneous terminations of complex AF simulations, which showed less predictable termination patterns lasting only 1,600 ms. This study contributes to the interpretation of previous clinical observations, and illustrates how computer modeling provides a complementary approach to study the mechanisms of cardiac arrhythmias.

3D virtual human atria: A computational platform for studying clinical atrial fibrillation

Despite a vast amount of experimental and clinical data on the underlying ionic, cellular and tissue substrates, the mechanisms of common atrial arrhythmias (such as atrial fibrillation, AF) arising from the functional interactions at the whole atria level remain unclear. Computational modelling provides a quantitative framework for integrating such multi-scale data and understanding the arrhythmogenic behaviour that emerges from the collective spatio-temporal dynamics in all parts of the heart. In this study, we have developed a multi-scale hierarchy of biophysically detailed computational models for the human atria--the 3D virtual human atria. Primarily, diffusion tensor MRI reconstruction of the tissue geometry and fibre orientation in the human sinoatrial node (SAN) and surrounding atrial muscle was integrated into the 3D model of the whole atria dissected from the Visible Human dataset. The anatomical models were combined with the heterogeneous atrial action potential (AP) models, and used to simulate the AP conduction in the human atria under various conditions: SAN pacemaking and atrial activation in the normal rhythm, break-down of regular AP wave-fronts during rapid atrial pacing, and the genesis of multiple re-entrant wavelets characteristic of AF. Contributions of different properties of the tissue to mechanisms of the normal rhythm and arrhythmogenesis were investigated. Primarily, the simulations showed that tissue heterogeneity caused the break-down of the normal AP wave-fronts at rapid pacing rates, which initiated a pair of re-entrant spiral waves; and tissue anisotropy resulted in a further break-down of the spiral waves into multiple meandering wavelets characteristic of AF. The 3D virtual atria model itself was incorporated into the torso model to simulate the body surface ECG patterns in the normal and arrhythmic conditions. Therefore, a state-of-the-art computational platform has been developed, which can be used for studying multi-scale electrical phenomena during atrial conduction and AF arrhythmogenesis. Results of such simulations can be directly compared with electrophysiological and endocardial mapping data, as well as clinical ECG recordings. The virtual human atria can provide in-depth insights into 3D excitation propagation processes within atrial walls of a whole heart in vivo, which is beyond the current technical capabilities of experimental or clinical set-ups.

Propagation of the Sinus Impulse Into the Koch Triangle and Localization, Timing, and Origin of the Multicomponent Potentials Recorded in This Area

Background—

The presence of a conduction block at the level of the Koch triangle (KT) and the origin of the multicomponent potentials inside this area are controversial issues. We investigated the propagation of the sinus impulse into the KT and the characteristics of multicomponent potentials recorded in that area in patients with and without atrioventricular nodal reentrant tachycardia (AVNRT).

Methods and Results—

Thirty-two patients (16 with AVNRT, 16 without AVNRT) underwent a sinus rhythm electroanatomic mapping of the right atrium (RA). Conduction velocities in the RA and in the KT were evaluated quantitatively on activation maps and qualitatively on isochronal and propagation maps. The presence, location, and timing of different types of multicomponent potentials were evaluated. A mean of 149±44 points were sampled in the RA, whereas a mean of 79±21 points were collected inside the KT. Propagation block at the level of crista terminalis was not found in any patient, whereas slow conduction inside the KT was found in all (median conduction velocity, 122 cm/s [110 to 135 cm/s] outside KT versus 60 cm/s [48 to 75 cm/s] inside KT; P <0.0001). Jackman potentials were identified inside KT in almost all the patients and were invariably found on the line of collision between the wavefronts activating the KT in opposite directions.

Conclusions—

No conduction block was detected inside the KT in patients with and without AVNRT. Conduction slowing was demonstrated during propagation of the sinus impulse inside the KT. The genesis of the Jackman potential may be related to the collision of the wavefronts activating KT in opposite directions.

Image-based models of cardiac structure in health and disease

Computational approaches to investigating the electromechanics of healthy and diseased hearts are becoming essential for the comprehensive understanding of cardiac function. In this article, we first present a brief review of existing image-based computational models of cardiac structure. We then provide a detailed explanation of a processing pipeline which we have recently developed for constructing realistic computational models of the heart from high resolution structural and diffusion tensor (DT) magnetic resonance (MR) images acquired ex vivo. The presentation of the pipeline incorporates a review of the methodologies that can be used to reconstruct models of cardiac structure. In this pipeline, the structural image is segmented to reconstruct the ventricles, normal myocardium, and infarct. A finite element mesh is generated from the segmented structural image, and fiber orientations are assigned to the elements based on DTMR data. The methods were applied to construct seven different models of healthy and diseased hearts. These models contain millions of elements, with spatial resolutions in the order of hundreds of microns, providing unprecedented detail in the representation of cardiac structure for simulation studies.

Optimizing local capture of atrial fibrillation by rapid pacing: study of the influence of tissue dynamics

While successful termination by pacing of organized atrial tachycardias has been observed in patients, rapid pacing of AF can induce a local capture of the atrial tissue but in general no termination. The purpose of this study was to perform a systematic evaluation of the ability to capture AF by rapid pacing in a biophysical model of the atria with different dynamics in terms of conduction velocity (CV) and action potential duration (APD). Rapid pacing was applied during 30 s at five locations on the atria, for pacing cycle lengths in the range 60-110% of the mean AF cycle length (AFCL(mean)). Local AF capture could be achieved using rapid pacing at pacing sites located distal to major anatomical obstacles. Optimal pacing cycle lengths were found in the range 74-80% AFCL(mean) (capture window width: 14.6 ± 3% AFCL(mean)). An increase/decrease in CV or APD led to a significant shrinking/stretching of the capture window. Capture did not depend on AFCL, but did depend on the atrial substrate as characterized by an estimate of its wavelength, a better capture being achieved at shorter wavelengths. This model-based study suggests that a proper selection of the pacing site and cycle length can influence local capture results and that atrial tissue properties (CV and APD) are determinants of the response to rapid pacing.

Incorporating histology into a 3D microscopic computer model of myocardium to study propagation at a cellular level

We introduce a 3D model of cardiac tissue to study at a microscopic level the relationship between tissue morphology and propagation of depolarization. Unlike the classical bidomain approach, in which tissue properties are described by the apparent conductivity of the tissue, in this "microdomain" approach, we included histology by modeling the actual shape of the intracellular and extracellular spaces that contain spatially distributed gap-junctions and membranes. The histological model of the tissue was generated by a computer algorithm that can be tuned to model different histological changes. For healthy tissue, the model predicted a realistic conduction velocity of 0.42 m/s based solely on the parameters derived from histology. A comparison with a brick-shaped, simplified model showed that conduction depended to a moderate extent on the shape of myocytes; a comparison with a one-dimensional bidomain model with the same overall shape and structure showed that the apparent conductivity of the tissue can be used to create an equivalent bidomain model. In summary, the microdomain approach offers a means of directly incorporating structural and functional parameters into models of cardiac activation and propagation and thus provides a valuable bridge between the cellular and tissue domains in the myocardium.

Structure specific models of electrical function in the right atrial appendage

Atrial fibrillation is the most common cardiac arrhythmia. It is prevalent in the elderly and contributes to mortality in congestive heart failure. Computer models of atrial electrical activation that incorporate realistic structure provide a means of investigating the mechanisms that initiate and maintain reentrant atrial arrhythmia. We are extending computational and experimental techniques that are well established in our laboratory to develop a detailed structure-based model of atrial electrical function. The 3D geometry of the atria and veno-atrial junctions is reconstructed from magnetic resonance images and detailed structure from specific regions of the atria is acquired using a semi-automated extended-volume imaging system. As an example of our approach, we present a reconstruction of the pig right atrial appendage (RAA), including the pectinate muscles (PM) and crista terminalis (CT). The RAA was embedded in wax and the block surface was serially etched, stained and imaged, then removed using an ultramiller to produce a uniformly-spaced image stack. Tissue was segmented and connected voxels were selected using a 3D region-growing algorithm to construct RAA geometry. Electrical activity has been modeled on this structure using the Courtemanche atrial cell activation model. A bidomain formulation was used employing a grid-based finite element solver. The RAA was activated by applying a stimulus (150 microA/mm3, 5 ms) to the 27 grid points at the top of the CT. Despite the complex structure of the PM, RAA activation was relatively uniform.

Mechanisms of transition from normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of tissue heterogeneity and anisotropy

Experimental evidence suggests that regional differences in action potential (AP) morphology can provide a substrate for initiation and maintenance of reentrant arrhythmias in the right atrium (RA), but the relationships between the complex electrophysiological and anatomical organization of the RA and the genesis of reentry are unclear. In this study, a biophysically detailed three-dimensional computer model of the right atrial tissue was constructed to study the role of tissue heterogeneity and anisotropy in arrhythmogenesis. The model of Lindblad et al. for a rabbit atrial cell was modified to incorporate experimental data on regional differences in several ionic currents (primarily, I(Na), I(CaL), I(K1), I(to), and I(sus)) between the crista terminalis and pectinate muscle cells. The modified model was validated by its ability to reproduce the AP properties measured experimentally. The anatomical model of the rabbit RA (including tissue geometry and fiber orientation) was based on a recent histological reconstruction. Simulations with the resultant electrophysiologically and anatomically detailed three-dimensional model show that complex organization of the RA tissue causes breakdown of regular AP conduction patterns at high pacing rates (>11.75 Hz): as the AP in the crista terminalis cells is longer, and electrotonic coupling transverse to fibers of the crista terminalis is weak, high-frequency pacing at the border between the crista terminalis and pectinate muscles results in a unidirectional conduction block toward the crista terminalis and generation of reentry. Contributions of the tissue heterogeneity and anisotropy to reentry initiation mechanisms are quantified by measuring action potential duration (APD) gradients at the border between the crista terminalis and pectinate muscles: the APD gradients are high in areas where both heterogeneity and anisotropy are high, such that intrinsic APD differences are not diminished by electrotonic interactions. Thus, our detailed computer model reconstructs complex electrical activity in the RA, and provides new insights into the mechanisms of transition from focal atrial tachycardia into reentry.

Automatically generated, anatomically accurate meshes for cardiac electrophysiology problems

Significant advancements in imaging technology and the dramatic increase in computer power over the last few years broke the ground for the construction of anatomically realistic models of the heart at an unprecedented level of detail. To effectively make use of high-resolution imaging datasets for modeling purposes, the imaged objects have to be discretized. This procedure is trivial for structured grids. However, to develop generally applicable heart models, unstructured grids are much preferable. In this study, a novel image-based unstructured mesh generation technique is proposed. It uses the dual mesh of an octree applied directly to segmented 3-D image stacks. The method produces conformal, boundary-fitted, and hexahedra-dominant meshes. The algorithm operates fully automatically with no requirements for interactivity and generates accurate volume-preserving representations of arbitrarily complex geometries with smooth surfaces. The method is very well suited for cardiac electrophysiological simulations. In the myocardium, the algorithm minimizes variations in element size, whereas in the surrounding medium, the element size is grown larger with the distance to the myocardial surfaces to reduce the computational burden. The numerical feasibility of the approach is demonstrated by discretizing and solving the monodomain and bidomain equations on the generated grids for two preparations of high experimental relevance, a left ventricular wedge preparation, and a papillary muscle.

A tissue-specific model of reentry in the right atrial appendage

INTRODUCTION

Atrial fibrillation is prevalent in the elderly and contributes to mortality in congestive heart failure. Development of computer models of atrial electrical activation that incorporate realistic structures provides a means of investigating the mechanisms that initiate and maintain reentrant atrial arrhythmia. As a step toward this, we have developed a model of the right atrial appendage (RAA) including detailed geometry of the pectinate muscles (PM) and crista terminalis (CT) with high spatial resolution, as well as complete fiber architecture.

METHODS AND RESULTS

Detailed structural images of a pig RAA were acquired using a semiautomated extended-volume imaging system. The generally accepted anisotropic ratio of 10:1 was adopted in the computer model. To deal with the regional action potential duration heterogeneity in the RAA, a Courtemanche cell model and a Luo-Rudy cell model were used for the CT and PM, respectively. Activation through the CT and PM network was adequately reproduced with acceptable accuracy using reduced-order computer models. Using a train of reducing cycle length stimuli applied to a CT/PM junction, we observed functional block both parallel with and perpendicular to the axis of the CT.

CONCLUSION

With stimulation from the CT at the junction of a PM, we conclude: (a) that conduction block within the CT is due to a reduced safety factor; and (b) that unidirectional block and reentry within the CT is due to its high anisotropy. Regional differences in effective refractive period do not explain the observed conduction block.

A morphologically realistic shell model of atrial propagation and ablation

A three dimensional morphologically realistic model of atrial propagation is developed, based on the male Visible Human dataset and the Fitzhugh-Nagumo equations for cardiac excitation. The atrial shell geometry incorporates eleven different anatomical structures, including the pulmonary veins, and the septum, Bachmann's bundle and coronary sinus as interatrial conduction pathways. Although the model utilizes a simplified cellular model of cardiac excitation it is able to reproduce a variety of electrophysiological phenomena including: autorhythmicity of the sinoatrial node and its ability to excite surrounding atrial tissue, spiral re-entrant wavefronts, ectopic beats originating in the PV and their termination by circumferential ablation of the PV. The model is an important tool to quantitatively study atrial excitation under normal conditions and during atrial fibrillation.

Selective site pacing: rationale and practical application

Although it has become traditional to place permanent pacemaker leads at the right ventricular apex and right atrial appendage, pacing from these locations poorly mimics normal physiology. A growing evidence base shows that right ventricular apical pacing results in ventricular dyssynchrony and various adverse effects. Provocative data from early trials suggest that pacing from alternate sites in the right ventricle--His bundle pacing, para-Hisian pacing, septal right ventricular outflow tract pacing, and right ventricular midseptal pacing--may lead to improved results. Similarly, early data suggest that right atrial pacing near Bachmann's bundle may lead to superior outcomes when compared with pacing from the right atrial appendage. Several large-scale, randomized clinical trials are now under way to establish the future role of selective site pacing.

Properties of two human atrial cell models in tissue: restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.

Cardiac electrical dynamics: maximizing dynamical heterogeneity

The relationships between key features of the cardiac electrical activity, such as electrical restitution, discordant alternans, wavebreak, and reentry, and the onset of ventricular tachyarrhythmias have been characterized extensively under the condition of constant rapid pacing. However, it is unlikely that this scenario applies directly to the clinical situation, where the induction of ventricular tachycardia (VT) typically is associated with the interruption of normal cardiac rhythm by several premature beats. To address this issue, we have developed a general theory to explain why specific patterns of premature stimuli increase dynamic heterogeneity of repolarization and precipitate conduction block. The theory predicts that conduction block is caused by (1) creation of a spatial gradient in diastolic interval (DI) by waves traveling at slightly different velocities (ie, conduction velocity dispersion) and (2) amplification of the spatial gradient in DI over subsequent action potentials, secondary to a strong dependence of action potential duration on the preceding DI (ie, a steep action potential duration restitution function). Tests of this theory have been conducted in computer models of homogeneous tissue, where increased spatial dispersion of repolarization during premature stimulation can be attributed solely to the development of dynamical heterogeneity, and in a canine model exhibiting spontaneously occurring VT and sudden death. Our results thus far indicate that the probability of inducing ventricular fibrillation (VF) in the animal model is highest for those sequences predicted to cause conduction block in the computer model. An understanding of the mechanisms underlying these observations will help to identify key electrical phenomena in the onset of VT and fibrillation. Drug and electrical therapies can then be improved by targeting these specific phenomena.

Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence

Ventricular fibrillation is a lethal arrhythmia characterized by multiple wavelets usually starting from a single or figure-of-eight re-entrant circuit. Understanding the factors regulating vulnerability to the re-entry is essential for developing effective therapeutic strategies to prevent ventricular fibrillation. In this study, we investigated how pre-existing tissue heterogeneities and electrical restitution properties affect the initiation of re-entry by premature extrastimuli in two-dimensional cardiac tissue models. We studied two pacing protocols for inducing re-entry following the "sinus" rhythm (S1) beat: (1) a single premature (S2) extrastimulus in heterogeneous tissue; (2) two premature extrastimuli (S2 and S3) in homogeneous tissue. In the first case, the vulnerable window of re-entry is determined by the spatial dimension and extent of the heterogeneity, and is also affected by electrical restitution properties and the location of the premature stimulus. The vulnerable window first increases as the action potential duration (APD) difference between the inside and outside of the heterogeneous region increases, but then decreases as this difference increases further. Steeper APD restitution reduces the vulnerable window of re-entry. In the second case, electrical restitution plays an essential role. When APD restitution is flat, no re-entry can be induced. When APD restitution is steep, re-entry can be induced by an S3 over a range of S1S2 intervals, which is also affected by conduction velocity restitution. When APD restitution is even steeper, the vulnerable window is reduced due to collision of the spiral tips.

Reduction of P-Wave Duration and Successful Pulmonary Vein Isolation in Patients with Atrial Fibrillation

INTRODUCTION

We hypothesize that successful pulmonary vein (PV) isolation can shorten the P-wave duration in patients with atrial fibrillation (AF).

METHODS AND RESULTS

We recorded magnified surface electrocardiogram (ECG) and P-wave signal-averaged ECG using 12 electrode leads before and after 31 PV isolation procedures in 27 patients with AF. The patients were followed for 16 +/- 4 months. Repeat ablation studies documented failed PV isolation in seven patients with AF recurrences. At baseline, the maximal P-wave duration in patients without AF recurrence (161 +/- 7 msec) was slightly shorter than that in patients with AF recurrence (168 +/- 10 msec, P < 0.05). After ablation, patients without recurrence showed a significant reduction of P-wave duration from 161 +/- 7 msec to 151 +/- 8 msec (P < 0.0001). In contrast, no change of P-wave duration was noted in patients with recurrences. These findings were confirmed with signal averaged ECG of the P-waves. Three-dimensional (3-D) computer simulation using an atrial cell model showed that elimination of the muscle sleeves inside the PV resulted in a shortening of the P-wave duration and change of the terminal portion of the P-wave morphology.

CONCLUSIONS

A significant shortening of P-wave duration by P-wave signal-averaged ECG can be used as an indicator for successful PV isolation. These findings suggest that activation of the PV muscle sleeves may be an important component of the terminal portion of the P-wave on surface ECG.

Myocardial segment-specific model generation for simulating the electrical action of the heart

BACKGROUND

Computer models of the electrical and mechanical actions of the heart, solved on geometrically realistic domains, are becoming an increasingly useful scientific tool. Construction of these models requires detailed measurement of the microstructural features which impact on the function of the heart. Currently a few generic cardiac models are in use for a wide range of simulation problems, and contributions to publicly accessible databases of cardiac structures, on which models can be solved, remain rare. This paper presents to-date the largest database of porcine left ventricular segment microstructural architecture, for use in both electrical and mechanical simulation.

METHODS

Cryosectioning techniques were used to reconstruct the myofibre and myosheet orientations in tissue blocks of size ~15 x 15 x 15 mm, taken from the mid-anterior left ventricular freewall, of seven hearts. Tissue sections were gathered on orthogonal planes, and the angles of intersection of myofibres and myosheets with these planes determined automatically with a gradient intensity based algorithm. These angles were then combined to provide a description of myofibre and myosheet variation throughout the tissue, in a form able to be input to biophysically based computational models of the heart.

RESULTS

Several microstructural features were common across all hearts. Myofibres rotated through 141 +/- 18 degrees (mean +/- SD) from epicardium to endocardium, in near linear fashion. In the outer two-thirds of the wall sheet angles were predominantly negative, however, in the inner one-third an abrupt change in sheet angle, with reversal in sign, was seen in six of the seven hearts. Two distinct populations of sheets with orthogonal orientations often co-existed, usually with one population dominating. The utility of the tissue structures was demonstrated by simulating the passive and active electrical responses of two of the tissue blocks to current injection. Distinct patterns of electrical response were obtained in the two tissue blocks, illustrating the importance of testing model based predictions on a variety of tissue architectures.

CONCLUSION

This study significantly expands the set of geometries on which models of cardiac function can be solved.

Simulation of atrial electrophysiology and body surface potentials for normal and abnormal rhythm

The effect of different atrial electrical activation sequences (sinus rhythm and atrial flutter circling in the right atrium) on the body surface potentials is investigated in this study. A realistic volume conductor model consisting of atria, lungs, chest and blood masses is generated from image stacks recorded by magnetic resonance imaging. The electrical sources-the transmembrane potentials-within the atrial volumetric model are simulated for different atrial rhythms employing a cellular automaton capable of considering different parameters depending on the specific properties of the tissues. The potentials on the torso surface are computed applying the finite element method for solving the differential equations derived from the bidomain theory. Both the simulated atrial activation patterns and the computed torso potentials for atrial sinus rhythm and atrial flutter are in qualitatively and quantitatively good agreement with data measured in humans. The simulation of body surface potentials generated by different electrical activation sequences in the atria or ventricles allows testing and assessing noninvasive imaging of cardiac electrophysiology, as both the potentials on the body surface and the reference activation in the heart are available.

Fibrillatory conduction in branching atrial tissue--Insight from volumetric and monolayer computer models

Increased local load in branching atrial tissue (muscle fibers and bundle insertions) influences wave propagation during atrial fibrillation (AF). This computer model study reveals two principal phenomena: if the branching is distant from the driving rotor (>19 mm), the load causes local slowing of conduction or wavebreaks. If the driving rotor is close to the branching, the increased load causes first a slow drift of the rotor towards the branching. Finally, the rotor anchors, and a stable, repeatable pattern of activation can be observed. Variation of the bundle geometry from a cylindrical, volumetric structure to a flat strip of a comparable load in a monolayer model changed the local activation sequence in the proximity of the bundle. However, the global behavior and the basic effects are similar in all models. Wavebreaks in branching tissue contribute to the chaotic nature of AF (fibrillatory conduction). The stabilization (anchoring) of driving rotors by branching tissue might contribute to maintain sustained AF.

A biophysical model of atrial fibrillation to define the appropriate ablation pattern in modified maze

OBJECTIVE

The surgical Maze III procedure remains the gold standard in treating atrial fibrillation (AF); however due to clinical difficulties and higher risks, less invasive ablation alternatives are clinically investigated. The present study aims to define more efficient ablation patterns of the modified maze procedure using a biophysical model of human atria with chronic AF.

METHODS

A three-dimensional model of human atria was developed using both MRI-imaging and a one-layer cellular model reproducing experimentally observed atrial cellular properties. Sustained AF could be induced by a burst-pacing protocol. Ablation lines were implemented in rendering the cardiac cells non-conductive, mimicking transmural lines. Lines were progressively implemented respectively around pulmonary veins (PV), left atrial appendage (LAA), left atrial isthmus (LAI), cavo-tricuspid isthmus (CTI), and intercaval lines (SIVC) in the computer model, defining the following patterns: P1=PV, P2=P1+LAA, P3=P2+LAI, P4=P3+CTI, P5=P3+SIVC, P6=P5+CTI. Forty simulations were done for each pattern and proportion of sinus rhythm (SR) conversion and time-to-AF termination (TAFT) were assessed.

RESULTS

The most efficient patterns are P5, P6, and Maze III with 100% success. The main difference is expressed in decreasing mean TAFT with a correlation coefficient R=-0.8. There is an inflexion point for 100% success rate at a 7.5s TAFT, meaning that no additional line is mandatory beyond pattern P5.

CONCLUSIONS

Our biophysical model suggests that Maze III could be simplified in his right atrial pattern to a single line joining both vena cavae. This has to be confirmed in clinical settings.

Specific electrocardiographic markers of P-wave morphology in interatrial block

BACKGROUND

Interatrial block (IAB; P waves >/= 110 milliseconds), the conduction delay between the right (RA) and left atrium (LA), is depicted on the electrocardiogram (ECG) as prolonged, often bifid ("notched"), P waves with distinguishable RA and LA components. Although electrophysiologic (EP) studies give some insight on how RA and LA components are depicted on the surface ECG in normal conduction, few if at all any, have conclusively demonstrated this correlation with IAB. Using existing EP knowledge, we investigated if such P-wave markers on bedside ECGs exist in IAB and appraised their utility in IAB recognition.

METHODS

We reviewed the medical records of patients admitted to a general hospital from December 1, 2004, to December 15, 2004. Of those, 151 patients had been admitted for nonacute presentations and were screened with 12-lead ECGs. Thirty-eight ECGs were excluded for nonsinus and paced rhythms, severe motion artifact, errors in lead placement, absence of adequate patient identification, and duplicate patient admissions after discharge. The remaining 113 ECGs were then evaluated for IAB. Sixty-three patients who did not have IAB formed the control (group A), whereas of the remaining 50 patients with IAB, 24 who had past ECGs for comparison formed group B1 and 26 without past ECGs formed group B2. Groups were compared for common clinical comorbidities, whereas sensitivity and specificity were calculated for significant P-wave markers. P values were also calculated, with a value of <.05 considered significant.

RESULTS

Clinical characteristics of patients in all groups were statistically comparable. Overall, almost all P waves in patients with IAB (groups 1 and 2) appeared "notched" (94%, P < .0001; sensitivity, 75%; specificity, 94% for IAB recognition; positive predictive value, 94%). P-wave RA components were commonly depicted as "domes," whereas their LA counterparts formed "spikes" (48%, P < .0001; sensitivity 96%; specificity, 70% for IAB recognition). When groups B1 and B2 were compared with increased accuracy, more P waves in group B1 were noted to have notches and had easily discernible RA and LA components; often, the RA duration is longer than the LA duration. In addition, more "dome-and-spike" complexes could be determined when past ECGs were present for comparison. These markers could be found on any bedside ECG lead in IAB but were predominant on leads II and V3 to V6.

CONCLUSIONS

Specific noninvasive surface markers such as P-wave "dome-and-spike" complexes and "notches" in any lead (predominantly leads II and V3-V6) on the bedside ECG could alert clinicians to measure P waves and so identify IAB.