Modeling ventricular repolarization: effects of transmural and apex-to-base heterogeneities in action potential durations

Heterogeneities in the densities of membrane ionic currents of myocytes cause regional variations in action potential duration (APD) at various intramural depths and along the apico-basal and circumferential directions in the left ventricle. This work extends our previous study of cartesian slabs to ventricular walls shaped as an ellipsoidal volume and including both transmural and apex-to-base APD heterogeneities. Our 3D simulation study investigates the combined effect on repolarization sequences and APD distributions of: (a) the intrinsic APD heterogeneity across the wall and along the apex-to-base direction, and (b) the electrotonic currents that modulate the APDs when myocytes are embedded in a ventricular wall with fiber rotation and orthotropic anisotropy. Our findings show that: (i) the transmural and apex-to-base heterogeneities have only a weak influence on the repolarization patterns on myocardial layers parallel to the epicardium; (ii) the patterns of APD distribution on the epicardial surface are mostly affected by the apex-to-base heterogeneities and do not reveal the APD transmural heterogeneity; (iii) the transmural heterogeneity is clearly discernible in both repolarization and APD patterns only on transmural sections; (iv) the apex-to-base heterogeneity is clearly discernible only in APD patterns on layers parallel to the epicardium. Thus, in our orthotropic ellipsoidal wall, the complex 3D electrotonic modulation of APDs does not fully mix the effects of the transmural and apex-to-base heterogeneity. The intrinsic spatial heterogeneity of the APDs is unmasked in the modulated APD patterns only in the appropriate transmural or intramural sections. These findings are independent of the stimulus location (epicardial, endocardial) and of Purkinje involvement.     

Repolarization of canine ventricular myocardium under the supraventricular rhythm

The ventricular myocardium is characterized by heterogeneity of activation-recovery interval durations. The transmural ARI gradients are present in the right ventricular apex (ARIs monotonically decreased as one moved from the endocardium to the epicardium), and in the left ventricular base (repolarization in the subepicardial layers was significantly shorter than that in the midmyo cardial layers whereas subendocardial ARIs did not differ from the others). The repolarization pattern of these myocardial regions is governed by the distribution of ARIs. In the apical left ventricular and basal right ventricular areas, no significant transmural differences in the repolarization durations were found. The repolarization pattern of these myocardial regions is governed by the activation sequence. In the right ventricle, ARIs were significantly longer at the base and shorter at the apex. In contrast, in the left ventricle, the apical ARIs were prolonged whereas the basal ARIs were abbreviated. The apex-to-base sequence of myocardial repolarization seems to depend on apex-to-base gradient of activation-recovery intervals durations.     

Effect of activation sequence on transmural patterns of repolarization and action potential duration in rabbit ventricular myocardium

Although transmural heterogeneity of action potential duration (APD) is established in single cells isolated from different tissue layers, the extent to which it produces transmural gradients of repolarization in electrotonically coupled ventricular myocardium remains controversial. The purpose of this study was to examine the relative contribution of intrinsic cellular gradients of APD and electrotonic influences to transmural repolarization in rabbit ventricular myocardium. Transmural optical mapping was performed in left ventricular wedge preparations from eight rabbits. Transmural patterns of activation, repolarization, and APD were recorded during endocardial and epicardial stimulation. Experimental results were compared with modeled data during variations in electrotonic coupling. A transmural gradient of APD was evident during endocardial stimulation, which reflected differences previously seen in isolated cells, with the longest APD at the endocardium and the shortest at the epicardium (endo: 165 ± 5 vs. epi: 147 ± 4 ms; P < 0.05). During epicardial stimulation, this gradient reversed (epi: 162 ± 4 vs. endo: 148 ± 6 ms; P < 0.05). In both activation sequences, transmural repolarization followed activation and APD shortened along the activation path such that significant transmural gradients of repolarization did not occur. This correlation between transmural activation time and APD was recapitulated in simulations and varied with changes in intercellular coupling, confirming that it is mediated by electrotonic current flow between cells. These data suggest that electrotonic influences are important in determining the transmural repolarization sequence in rabbit ventricular myocardium and that they are sufficient to overcome intrinsic differences in the electrophysiological properties of the cells across the ventricular wall.     

Transmural and apicobasal gradients in repolarization contribute to T-wave genesis in human surface ECG

The cellular basis of the T-wave morphology of surface ECG remains controversial in clinical cardiology. We examined the effect of action potential duration (APD) distribution on T-wave morphology using a realistic model of the human ventricle and torso. We developed a finite-element model of the ventricle consisting of ～26 million elements, including the conduction system, each implemented with the ion current model of cardiomyocytes. This model was embedded in a torso model with distinct organ structures to obtain the standard ECG leads. The APD distribution was changed in the transmural direction by locating the M cells in either the endocardial or epicardial region. We also introduced apicobasal gradients by modifying the ion channel parameters. Both the transmural gradient (with M cells on the endocardial side) and the apicobasal gradient produced positive T waves, although a very large gradient was required for the apicobasal gradient. By contrast, T waves obtained with the transmural gradient were highly symmetric and, therefore, did not represent the true physiological state. Only combination of the transmural and the moderate apicobasal gradients produced physiological T waves in surface ECG. Positive T waves in surface ECG mainly originated from the transmural distribution of APD with M cells on the endocardial side, although the apicobasal gradient was also required to attain the physiological waveform.     

Electrophysiology of Heart Failure Using a Rabbit Model: From the Failing Myocyte to Ventricular Fibrillation

Heart failure is a leading cause of death, yet its underlying electrophysiological (EP) mechanisms are not well understood. In this study, we use a multiscale approach to analyze a model of heart failure and connect its results to features of the electrocardiogram (ECG). The heart failure model is derived by modifying a previously validated electrophysiology model for a healthy rabbit heart. Specifically, in accordance with the heart failure literature, we modified the cell EP by changing both membrane currents and calcium handling. At the tissue level, we modeled the increased gap junction lateralization and lower conduction velocity due to downregulation of Connexin 43. At the biventricular level, we reduced the apex-to-base and transmural gradients of action potential duration (APD). The failing cell model was first validated by reproducing the longer action potential, slower and lower calcium transient, and earlier alternans characteristic of heart failure EP. Subsequently, we compared the electrical wave propagation in one dimensional cables of healthy and failing cells. The validated cell model was then used to simulate the EP of heart failure in an anatomically accurate biventricular rabbit model. As pacing cycle length decreases, both the normal and failing heart develop T-wave alternans, but only the failing heart shows QRS alternans (although moderate) at rapid pacing. Moreover, T-wave alternans is significantly more pronounced in the failing heart. At rapid pacing, APD maps show areas of conduction block in the failing heart. Finally, accelerated pacing initiated wave reentry and breakup in the failing heart. Further, the onset of VF was not observed with an upregulation of SERCA, a potential drug therapy, using the same protocol. The changes introduced at the cell and tissue level have increased the failing heart’s susceptibility to dynamic instabilities and arrhythmias under rapid pacing. However, the observed increase in arrhythmogenic potential is not due to a steepening of the restitution curve (not present in our model), but rather to a novel blocking mechanism.     

Heterogeneous connexin43 expression produces electrophysiological heterogeneities across ventricular wall

Recently we found that electrophysiological (EP) heterogeneities between subepicardial and midmyocardial cells can form a substrate for reentrant ventricular arrhythmias. However, cell-to-cell coupling through gap junctions is expected to attenuate transmural heterogeneities between cell types spanning the ventricular wall. Because connexin43 (Cx43) is the principal ventricular gap junction protein, we hypothesized that transmural EP heterogeneities are in part produced by heterogeneous Cx43 expression across the ventricular wall. The left ventricles of eight dogs were sectioned to expose the transmural surface. To determine whether heterogeneous Cx43 expression influenced EP function, high-resolution transmural optical mapping of the arterially perfused canine wedge preparation was used to measure transmural conduction velocity (thetaTM), dV/dt(max), transmural space constant (lambdaTM), and transmural gradients of action potential duration (APD). Relative Cx43 expression, quantified by confocal immunofluorescence, was significantly lower (by 24 +/- 17%; P < 0.05) in subepicardial compared with deeper layers. Importantly, reduced subepicardial Cx43 was associated with transmural heterogeneities of EP function evidenced by selectively reduced subepicardial thetaTM (by 18 +/- 9%; P < 0.05) compared with deeper layers. In subepicardial regions, dV/dt(max) was fastest (by 19 +/- 15%) and lambdaTM was smallest (by 18.1 +/- 2%), which suggests that conduction slowing was attributable to localized uncoupling rather than reduced excitability. The maximum transmural APD gradients occurred in the same regions where Cx43 expression was lowest; this suggests that Cx43 expression patterns served to maintain APD gradients across the transmural wall. These data demonstrate that heterogeneous Cx43 expression is closely associated with functionally significant EP heterogeneities across the transmural wall. Therefore, Cx43 expression patterns can potentially contribute to arrhythmic substrates that are dependent on transmural electrophysiological heterogeneities.     

Modeling of IK1 mutations in human left ventricular myocytes and tissue

Elucidation of the cellular basis of arrhythmias in ion channelopathy disorders is complicated by the inherent difficulties in studying human cardiac tissue. Thus we used a computer modeling approach to study the mechanisms of cellular dysfunction induced by mutations in inward rectifier potassium channel (K(ir))2.1 that cause Andersen-Tawil syndrome (ATS). ATS is an autosomal dominant disorder associated with ventricular arrhythmias that uncommonly degenerate into the lethal arrhythmia torsade de pointes. We simulated the cellular and tissue effects of a potent disease-causing mutation D71V K(ir)2.1 with mathematical models of human ventricular myocytes and a bidomain model of transmural conduction. The D71V K(ir)2.1 mutation caused significant action potential duration prolongation in subendocardial, midmyocardial, and subepicardial myocytes but did not significantly increase transmural dispersion of repolarization. Simulations of the D71V mutation at shorter cycle lengths induced stable action potential alternans in midmyocardial, but not subendocardial or subepicardial cells. The action potential alternans was manifested as an abbreviated QRS complex in the transmural ECG, the result of action potential propagation failure in the midmyocardial tissue. In addition, our simulations of D71V mutation recapitulate several key ECG features of ATS, including QT prolongation, T-wave flattening, and QRS widening. Thus our modeling approach faithfully recapitulates several features of ATS and provides a mechanistic explanation for the low frequency of torsade de pointes arrhythmia in ATS.     

Effects of ectopic pacing on repolarization of the chicken left ventricle

Effects of ectopic pacing on left ventricular repolarization were studied in six anesthetized open-chest chickens. In each animal, unipolar electrograms were acquired from as many as 98 sites with 14 plunge needles (seven transmural locations between epicardium and endocardium in each needle). Activation-recovery intervals (ARIs), corrected to the cycle length, were used for estimating repolarization. At baseline, the nonuniform ARI distribution in the left ventricle resulted in the apicobasal differences being greater than the transmural gradient. Nonuniform ARI prolongation caused by ectopic pacing resulted in decreasing the transmural repolarization gradient and increasing the differences in the apex-to-base direction. The basal, but not apical transmural differences contributed to the total left ventricular transmural gradient. The total left ventricular apicobasal gradient was contributed by the apicobasal differences in mid-myocardial and subendocardial layers more than in subepicardial ones. Thus, in in situ chicken hearts, the transmural and apicobasal ARI gradients exist within the left ventricle with the shortest ARIs in the basal subepicardium and the longest ARIs in the subendocardium of the apical and middle parts of the left ventricle. Apicobasal compared to transmural heterogeneity of local repolarization properties contributes more to the total left ventricular repolarization gradient.     

Torsion of the human left ventricle: experimental analysis and computational modeling

We set a twofold investigation: we assess left ventricular (LV) rotation and twist in the human heart through 3D-echocardiographic speckle tracking, and use representative experimental data as benchmark with respect to numerical results obtained by solving our mechanical model of the LV. We aim at new insight into the relationships between myocardial contraction patterns and the overall behavior at the scale of the whole organ. It is concluded that torsional rotation is sensitive to transmural gradients of contractility which is assumed linearly related to action potential duration (APD). Pressure-volume loops and other basic strain measures are not affected by these gradients. Therefore, realistic torsional behavior of human LV may indeed correspond to the electrophysiological and functional differences between endocardial and epicardial cells recently observed in non-failing hearts. Future investigations need now to integrate the mechanical model proposed here with minimal models of human ventricular APD to drive excitation-contraction coupling transmurally.     

Subepicardial phase 0 block and discontinuous transmural conduction underlie right precordial ST-segment elevation by a SCN5A loss-of-function mutation

Two mechanisms are generally proposed to explain right precordial ST-segment elevation in Brugada syndrome: 1) right ventricular (RV) subepicardial action potential shortening and/or loss of dome causing transmural dispersion of repolarization; and 2) RV conduction delay. Here we report novel mechanistic insights into ST-segment elevation associated with a Na(+) current (I(Na)) loss-of-function mutation from studies in a Dutch kindred with the COOH-terminal SCN5A variant p.Phe2004Leu. The proband, a man, experienced syncope at age 22 yr and had coved-type ST-segment elevations in ECG leads V1 and V2 and negative T waves in V2. Peak and persistent mutant I(Na) were significantly decreased. I(Na) closed-state inactivation was increased, slow inactivation accelerated, and recovery from inactivation delayed. Computer-simulated I(Na)-dependent excitation was decremental from endo- to epicardium at cycle length 1,000 ms, not at cycle length 300 ms. Propagation was discontinuous across the midmyocardial to epicardial transition region, exhibiting a long local delay due to phase 0 block. Beyond this region, axial excitatory current was provided by phase 2 (dome) of the M-cell action potentials and depended on L-type Ca(2+) current ("phase 2 conduction"). These results explain right precordial ST-segment elevation on the basis of RV transmural gradients of membrane potentials during early repolarization caused by discontinuous conduction. The late slow-upstroke action potentials at the subepicardium produce T-wave inversion in the computed ECG waveform, in line with the clinical ECG.     

Dynamical and cellular electrophysiological mechanisms of ECG changes during ischaemia

The interpretation of normal and pathological electrocardiographic (ECG) patterns in terms of the underlying cellular and tissue electrophysiology is rudimentary, as the existing theories rely on geometrical aspects. We relate effects of sub-endocardial ischaemia on the ST-segment depression in ECG to patterns of transmural action potential propagation in a one-dimensional virtual ventricular wall. Our computational study exposes two electrophysiological mechanisms of ST depression: dynamic-predominantly positive spatial gradients in the membrane potential during abnormal repolarization of the wall, produced by action potential duration changes in the ischaemic region; and static-a negative spatial gradient of the resting membrane potential between the normal and ischaemic regions. Hyperkalaemia is the major contributor to both these mechanisms at the cellular level. These results complement simulations of the effects of cardiac geometry on the ECG, and dissect spatio-temporal and cellular electrophysiological mechanisms of ST depression seen in sub-endocardial ischaemia.     

Effects of transmural electrical heterogeneities and electrotonic interactions on the dispersion of cardiac repolarization and action potential duration: A simulation study

It has been shown in the literature that myocytes isolated from the ventricular walls at various intramural depths have different action potential durations (APDs). When these myocytes are embedded in the ventricular wall, their inhomogeneous properties affect the sequence of repolarization and the actual distribution of the APDs in the entire wall. In this article, we implement a mathematical model to simulate the combined effect of (a) the non-homogeneous intrinsic membrane properties (in particular the non-homogeneous APDs) and (b) the electrotonic currents that modulate the APDs when the myocytes are embedded in the ventricular myocardium. In particular, we study the effect of (a) and (b) on the excitation and repolarization sequences and on the distribution of APDs in the ventricles. We implement a Monodomain tissue representation that includes orthotropic anisotropy, transmural fiber rotation and homogeneous or heterogeneous transmural intrinsic membrane properties, modeled according to the phase I Luo-Rudy membrane ionic model. Three-dimensional simulations are performed in a cartesian slab with a parallel finite element solver employing structured isoparametric trilinear finite elements in space and a semi-implicit adaptive method in time. Simulations of excitation and repolarization sequences elicited by epicardial or endocardial pacing show that in a homogeneous slab the repolarization pathways approximately follow the activation sequence. Conversely, in the heterogeneous cases considered in this study, we observed two repolarization wavefronts that started from the epi and the endocardial faces respectively and collided in the thickness of the wall and in one case an additional repolarization wave starting from an intramural site. Introducing the heterogeneities along the transmural epi-endocardial direction affected both the repolarization sequence and the APD dispersion, but these effects were clearly discernible only in transmural planes. By contrast, in planes parallel to epi- and endocardium the APD distribution remained remarkably similar to that observed in the homogeneous model. Therefore, the patterns of the repolarization sequence and APD dispersion on the epicardial surface (or any other intramural surface parallel to it) do not reveal the uniform transmural heterogeneity.     

面向短QT综合征的药物作用建模与仿真研究

短QT综合征(short QT syndrome,SQTS)是以心电图QT间期、心室和心房不应期明显缩短为主要显性特征,并伴有晕厥、高发心源性猝死(sudden cardiac death,SCD)和恶性心律失常风险的一类遗传性心肌离子通道病.据目前资料信息,关于SQTS致病机理的报道比较多,而对SQTS药物治疗的报道罕见.为了揭示在SQTS下的药物作用,本文通过计算机仿真构建人体心室细胞和组织的药物作用模型,利用该模型,从亚细胞、细胞、组织三个尺度,模拟SQT1、SQT2和SQT3下的普罗帕酮药物作用过程,并仿真心电图的变化情况.仿真结果表明:在SQT1下普罗帕酮延长了动作电位时程(action potential duration,APD)和心电图QT间期,并降低T波幅值;相反,在SQT2和SQT3下普罗帕酮缩短了APD和QT间期.计算使用药物前后细胞间膜电压和APD空间离散度的变化,定量分析了普罗帕酮降低T波振幅的原因.总之,对SQT1,普罗帕酮有效;对SQT2和SQT3,普罗帕酮没有改变其致心律失常的危险.仿真结果为普罗帕酮用于临床治疗SQTS提供理论参考.     

Nonequilibrium Arrhythmic States and Transitions in a Mathematical Model for Diffuse Fibrosis in Human Cardiac Tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to fibroblasts via a gap-junctional conductance , reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.     

Ventricular hypertrophy amplifies transmural repolarization dispersion and induces early afterdepolarization

The effects of left ventricular hypertrophy (LVH) on the generation of phase 2 early afterdepolarization (EAD) and transmural dispersion of repolarization (TDR) were assessed using arterially perfused rabbit ventricular wedge preparations. Transmembrane action potentials from epicardium, subendocardium, and endocardium were simultaneously recorded together with a transmural ECG. Transmural action potential duration (APD) was also mapped. LVH (renovascular hypertension model) produced significant prolongation in ventricular APD and QT interval. Preferential APD prolongation in subendocardium and endocardium was associated with a marked increase in TDR. Phase 2 EADs were generated from subendocardium or endocardium in all LVH rabbits (15 of 15) in the absence of APD prolonging agents at basic cycle lengths of 2,000-4,000 ms. Phase 2 EAD could produce "R on T" extrasystoles, initiating polymorphic ventricular tachycardia (VT). This study provides the first direct evidence from intracellular recordings that phase 2 EAD could be generated from rabbit intact hypertrophied LV wall in the absence of APD prolonging agents, resulting in R on T extrasystoles capable of initiating polymorphic VT under enhanced TDR.     

Propagation of normal beats and re-entry in a computational model of ventricular cardiac tissue with regional differences in action potential shape and duration

There is substantial experimental evidence from studies using both intact tissue and isolated single cells to support the existence of different cell types within the ventricular wall of the heart, each possessing different electrical properties. However other studies have failed to find these differences, and instead support the idea that electrical coupling in vivo between regions with different cell types smoothes out differences in action potential shape and duration. In this study we have used a computational model of electrical activation in heterogenous 2D and 3D cardiac tissue to investigate the propagation of both normal beats and arrhythmias. We used the Luo–Rudy dynamic model for guinea pig ventricular cells, with simplified Ca²⁺ handling and transmural heterogeneity in I_Ks and I_to. With normal cell-to-cell coupling, a layer of M cells was not necessary for the formation of an upright T wave in the simulated electrocardiogram, and the amplitude and configuration of the T wave was not greatly affected by the thickness and configuration of the M cell layer. Transmural gradients in repolarisation pushed re-entrant waves with an intramural filament towards either the base or the apex of the ventricles, and caused transient break up of re-entry with a transmural filament.     

Chemical ablation of the Purkinje system causes early termination and activation rate slowing of long-duration ventricular fibrillation in dogs

Endocardial mapping has suggested that Purkinje fibers may play a role in the maintenance of long-duration ventricular fibrillation (LDVF). To determine the influence of Purkinje fibers on LDVF, we chemically ablated the Purkinje system with Lugol solution and recorded endocardial and transmural activation during LDVF. Dog hearts were isolated and perfused, and the ventricular endocardium was exposed and treated with Lugol solution (n = 6) or normal Tyrode solution as a control (n = 6). The left anterior papillary muscle endocardium was mapped with a 504-electrode (21 x 24) plaque with electrodes spaced 1 mm apart. Transmural activation was recorded with a six-electrode plunge needle on each side of the plaque. Ventricular fibrillation (VF) was induced, and perfusion was halted. LDVF spontaneously terminated sooner in Lugol-ablated hearts than in control hearts (4.9 +/- 1.5 vs. 9.2 +/- 3.2 min, P = 0.01). After termination of VF, both the control and Lugol hearts were typically excitable, but only short episodes of VF could be reinduced. Endocardial activation rates were similar during the first 2 min of LDVF for Lugol-ablated and control hearts but were significantly slower in Lugol hearts by 3 min. In control hearts, the endocardium activated more rapidly than the epicardium after 4 min of LDVF with wave fronts propagating most often from the endocardium to epicardium. No difference in transmural activation rate or wave front direction was observed in Lugol hearts. Ablation of the subendocardium hastens VF spontaneous termination and alters VF activation sequences, suggesting that Purkinje fibers are important in the maintenance of LDVF.     

Sequence of ventricular repolarization under body cooling in the frog

Lowering the temperature is known to prolong the repolarization of cardiomyocytes. However, whether the prolongation of action potentials is uniform throughout the myocardium, and whether this prolongation is obvious in ECG, remains unclear. Ventricular repolarization sequences and body surface potential distributions were studied in 20 frogs Rana temporaria using epicardial and body surface potential mapping. An apex-to-base ventricular repolarization sequence corresponded to the distribution of local repolarization durations was demonstrated at the temperature of 18 degrees C. The body surface potential distribution during the ST-T complex was characterized by the cranial negative and caudal positive potential areas. Under the body cooling to 10 degrees C, repolarization prolonged to a greater extent at the apex that resulted in the base-to-apex repolarization sequence, which, in turn, caused an inversion in the body surface potential distribution with cranial portion of the body being positive and caudal portion being negative.     

Electrophysiology and anatomy of embryonic rabbit hearts before and after septation

Mechanisms of cardiac pacemaking and conduction system (CPCS) development are difficult to study, in part because of the absence of models that are physiologically similar to humans in which we can label the entire CPCS. Investigations of the adult rabbit heart have provided insight into normal and abnormal cardiac conduction. The adult and the embryonic rabbit have an endogenous marker of the entire cardiac conduction system, neurofilament 160 (NF-160). Previous work suggested that ventricular septation correlates with critical phases in avian CPCS development, in contrast to the mouse CPCS. Combining high-resolution optical mapping with immunohistochemical analysis of the embryonic rabbit heart, we investigated the significance of ventricular septation in patterning the rabbit embryonic conduction system. We hypothesized that 1) completion of ventricular septation does not correlate with changes in the ventricular activation sequence in rabbit embryos and 2) CPCS anatomy determines the activation sequence of the embryonic heart. We found that preseptated (days 11-13, n = 13) and postseptated (day 15, n = 5) hearts had similar "apex-to-base" ventricular excitation. PR intervals were not significantly different in either group. CPCS anatomy revealed continuity of the NF-160-positive tract connecting the presumptive sinoatrial node, atrioventricular (AV) junction, and ventricular conduction system. The presence of collagen in the AV junction coincided with the appearance of an AV interval. We conclude that the apex-to-base ventricular activation sequence in the rabbit embryo is present before completion of ventricular septation. CPCS anatomy reflects global cardiac activation as demonstrated by high-resolution optical mapping.