|Home | About | Journals | Submit | Contact Us | Français|
Understanding sinoatrial node (SAN) development could help in developing therapies for SAN dysfunction. However, electrophysiological investigation of SAN development remains difficult because mutant mice with SAN dysfunctions are frequently embryonically lethal. Most research on SAN development is therefore limited to immunocytochemical observations without comparable functional studies.
We applied a multi-electrode array (MEA) recording system to study SAN development in mouse hearts acutely isolated at embryonic ages (E) 8.5 to 12.5 days. Physiological heart rates were routinely restored, enabling accurate functional assessment of SAN development. We found that dominant pacemaking activity originated from the left inflow tract (LIFT) region at E8.5, but switched to the right SAN by E12.5. Combining MEA recordings and pharmacological agents, we show that intracellular calcium (Ca2+)-mediated automaticity develops early and is the major mechanism of pulse generation in the LIFT of E8.5 hearts. Later in development at E12.5, sarcolemmal ion channels develop in the SAN at a time when pacemaker channels are down regulated in the LIFT, leading to a switch in the dominant pacemaker location. Additionally, low micromolar concentrations of tetrodotoxin (TTX), a sodium channel blocker, minimally affect pacemaker rhythm at E8.5–12.5; but suppress atrial activation and reveal a tetrodotoxin-resistant SAN-atrioventricular node (internodal) pathway that mediates internodal conduction in E12.5 hearts.
Using a physiological mapping method, we demonstrate that differential mechanistic development of automaticity between the left and right inflow tract regions confers the pacemaker location switch. Moreover, a tetrodotoxin-resistant pathway mediates preferential internodal conduction in E12.5 mouse hearts.
Cardiac myocytes in the specialized cardiac conduction system display different electrical and immunocytochemical properties from surrounding working myocytes.1, 2 Mouse models with genetic manipulations are frequently used to study embryonic development of cardiac conduction system, yet functional correlates are usually lacking. Also, physiological mechanisms underlying automaticity in mature sinoatrial node (SAN) cells remain controversial and most likely involve various sarcolemmal ion channels, ion exchangers, and intracellular Ca2+ cycling machinery.2–5 Currently, the underlying mechanisms for functional development of automaticity in early pacemaking or SAN cells remains unclear. Recent research has just begun to unravel molecular pathways governing structural development of SAN.2, 6–8 Elucidation of molecular pathways directing electrophysiological development of the SAN will help in the development of genetic and cell-based therapies for treating sinus node dysfunction.9
Physiological studies of whole embryonic mouse hearts are usually challenging because of their small size and fragility. Lack of specific SAN tracking markers,10 especially at early embryonic stages, limits the interpretation of single cardiomyocyte-based studies. Moreover, mutant mice with defective SAN function are frequently embryonic-lethal and cannot be studied further.11 Therefore, most electrophysiological studies of functional development of pacemaker or conducting cells in early embryonic hearts have used voltage- or Ca2+-sensitive dyes and derived their results by optical mapping methods. These optical mapping studies reported that the earliest pacemaking activity in embryonic chick or rat hearts was usually found in the left sinoatrial primordium (termed LIFT in this manuscript); yet sometimes it originated in the right sinoatrial region.12, 13 While optical mappings have provided valuable insights toward functional development of cardiac conduction system,12–14 the dye-loading process required for optical mapping may affect electrophysiological properties of hearts15 or slow down heart rates,12,13 which could result in wandering pacemakers as observed in early studies.12,13 It is also unlikely that embryonic hearts from mutant mice with defective SAN development could survive the harsh dye-loading process. Genetically targetable or encoded Ca2+ indicators have recently been developed to avoid harsh dye-loading processes,16 yet their widespread applications remain to be demonstrated. Consequently, the mechanism underlying the switch of pacemaker location during early cardiogenesis remains to be elucidated. An electrophysiological recording method without any influence of foreign dyes may be the most physiological means to study SAN development during early cardiogenesis and should be complementary to various existing recording methods. Additionally, because most embryonic cardiomyocytes can generate automaticity,17 a system that constantly monitors dominant pacemaker location is needed so as to prevent misleading conclusions from wandering pacemakers originating from working cardiomyocytes.
Here, we applied a multi-electrode array (MEA) recording system18 with motion detection algorithms to determine the mechanism of switch in pacemaker location (referred to here as ‘pacemaker location switch’) and signal conduction between the SAN and atrioventricular node (AVN) during early cardiogenesis at the organ level.
Details of all protocols and antibodies used are provided in Online Supplements.
To establish an embryonic mouse heart preparation with physiological heart rhythm, pregnant mice (C57/B6) were sedated by repetitive i.p. injections of Avertin and ventilated with 95%O2/5% CO2 throughout the experimental period. Embryonic age of 0.5 day was designated at noon on the day a vaginal plug was detected after overnight mating. Individual embryos were dissected consecutively and embryonic hearts were then isolated with a dissecting microscope at room temperature. Isolated hearts were immediately recovered with 37°C pre-oxygenated Dulbecco’s modified eagle medium (DMEM) to obtain optimal physiological functions.
Each acutely isolated embryonic heart was placed directly onto a MEA chamber containing 60 gold electrodes (10–30 µm electrode with 20–200 µm inter-electrode spacing, Supplemental Fig. 1). All experiments were conducted at 37°C, and MEA chamber was continuously perfused with pre-oxygenated (95% O2/5% CO2) DMEM containing 1.8 mmol/L Ca2+ with or without drugs (1–2 ml/min). Spontaneous field potentials were recorded at 10kHz by the MEA data acquisition system. The sources of pharmacological agents and their applications are detailed in Online Supplements (Supplemental Figs. 2–4). The dominant pacemaker regions during experiments were constantly monitored and a switch of pacemaker location was indicated when it occurred.
For capturing the contraction sequence, each heart was imaged with a digital camera during electrophysiological recordings. The imaging data were analyzed with motion detection software provided by Dr. Martin Fink. Video imaging procedure and analysis are detailed in Online Supplement. Video imaging was also used to confirm that each electrode recorded from the same location of individual heart throughout each respective experiment (Supplemental Movie 1).
Standard techniques were used for immunostaining, cryosections, and in situ hybridization of embryonic mouse hearts.
Data are presented as the mean ± SD (standard deviation). Using the StatView program (SAS institute, NC), statistical differences were analyzed by an ANOVA procedure with post hoc Tukey/Kramer test for multiple comparisons. A statistical difference with a p-value < 0.05 was considered significant and was labeled with an asterisk (*) in all figures.
We applied the MEA recording system and video†imaging with a motion detection algorithm to analyze cardiac electrical activities and contractions of acutely isolated embryonic mouse hearts. MEA recordings of E8.5 (Fig. 1A) and E12.5 hearts (with dorsal SAN side contacting electrodes, Fig. 1B) are shown. Figures 1C–D depicts comparable field potential electrograms from various recording electrodes of respective hearts. MEA recordings illustrated a single field potential propagating across various regions of E8.5 hearts yet revealed distinct SAN, atrial and ventricular signals in E12.5 hearts. Based on the activation sequence of field potentials from various regions of embryonic hearts, the earliest location of electrical activities (termed the dominant pacemaker site) could be determined. The local activation time (LAT) of various cardiac regions is determined by the midpoint of maximal negative slope of field potentials at each recording electrode (Fig. 1E).19 At E12.5, when distinct signals with large amplitudes could be recorded, we can also transform respective local activation time of field potential signals from all electrodes to activation maps (similar to the activation map generated by the CARTO system in clinical cardiac electrophysiology)20 for a graphic demonstration of myocardial activation sequence using MATLAB software (Fig. 1F). In Figure 1F (and Supplemental Movie 2), the earliest activation (red) at E12.5 was located at the right-sided SAN (SAN map) with bilateral atrial signal propagation (A map) and double apical ventricular breakthroughs (V map).14 Video imaging with a motion detection algorithm is used in parallel to match the sidedness of earliest contraction to earliest electrical activation for confirming proper excitation-contraction coupling (Supplemental Fig. 5). This video imaging is especially useful at E8.5 when excessive manipulation of small hearts must be avoided due to their fragility. With our recording system, measured heart rates (HRs) from these isolated hearts are relatively stable for >30 minutes and are similar to measured in vivo HRs at comparable embryonic ages.21 The initially restored HRs and stable HRs at 10 minutes (min) are: 121.6±29.2 and 108.9±28.0 beats per minute (bpm) at E8.5 (n=67); 160.0±25.5 and 142.4±26.4 bpm at E10.5 (n=29); as well as 180.6±23.5 and 171.3±21.6 bpm at E12.5 (n=81), respectively (Supplemental Fig. 2). Importantly, we found that earliest pacemaking activity always started from the LIFT region at E8.5 and generally switched to the right SAN area by E12.5 (pacemaker location switch) (Fig. 1G). The respective percentages of dominant right-sided pacemaker at E8.5, 10.5, and 12.5 are 0% (n=80), 64.1% (n=39), and 88.4% (n=95). In rare cases (~4%) when the signal quality of video imaging or MEA recording was poor, we cited the origin of the dominant pacemaker as “undetermined.” Additionally, this MEA-based recording system allows efficient pharmacological investigation of pacemaking mechanisms during cardiogenesis at the organ level (see below).
Since many mechanisms contribute to automaticity of pacemaker cells, we used the MEA system to simultaneously monitor locations and rate changes of dominant pacemaker sites by pharmacological blockers of ion channels or intracellular Ca2+ ([Ca2+]i) handling proteins to elucidate developmental changes in mechanisms controlling dominant automaticity during early cardiogenesis. Dense mapping of E8.5 hearts (Fig. 1A) was used to construct Figures 2–3.
The [Ca2+]i handling machinery, including inositol 1,4,5-trisphosphate (IP3) receptor (IP3R),22 type 2 Ryanodine receptor (RyR2),23 and type 1 sodium (Na+)/ Ca2+ exchanger (NCX1),3, 11 had been implicated in the generation of pacemaker rhythms in adult mouse SAN cells. We applied various antagonists to determine which of these [Ca2+]i handling proteins contributed to automaticity at dominant pacemaker sites during early cardiogenesis. KB-R7943 (KB) at 50 µM is reported to significantly block outward current via NCX1 (the reverse mode) and modestly inhibit inward current through NCX1 (the forward mode).24 In Figures 2A–C, application of 50 µM KB minimally affected HRs of E8.5 hearts (4.7±4.9% blockade, n=12) and decreased HRs from E12.5 SAN by only 14.2±4.3% (n=14). 10 µM Ryanodine, a RyR2 blocker,25 decreased HRs of E8.5 (n=9) and E12.5 (n=16) hearts by 18.4±6.9% and 27.9±12.1% respectively. However, after Ryanodine or KB blockade, co-application of 50 µM KB and 10 µM Ryanodine for 10 minutes completely abolished HRs originating from LIFT in all E8.5 hearts tested (Zero HR at LIFT, n=8). In 3 of 8 E8.5 hearts, residual pacemaking activity could be detected only at the right inflow tract (RIFT) region (mean residual HR at RIFT is 15.9±25.0%) with arrested LIFT automaticity after exposure to these two drugs (Fig. 2B bottom left). The remaining cardiac rhythm from RIFT in these 3 hearts was slow, irregular and eventually ceased to function if an additional 10-minute recording was performed, suggesting immature RIFT pacemakers. These results indicate that [Ca2+]i signaling is the major mechanism controlling earliest left-sided automaticity in E8.5 hearts. Results in Figure 2C also showed that KB alone minimally affected HRs of E8.5 LIFT at baseline (<5% blockade) when Ryanodine receptors are fully functional, yet significantly blocked HRs (100% blockade) when Ryanodine receptors are deficient. This non-additive effect suggests that NCX1 acts as an automaticity reserve to compensate the loss of Ryanodine receptors (or vice versa) for maintaining cardiac rhythm at E8.5 LIFT. Also, this significant compensatory mechanism by various [Ca2+]i handling apparatus at E8.5 may explain the small effect of single drug, either KB (4.7±4.9% blockade) or Ryanodine alone (18.4±6.9% blockade), on HRs originating from LIFT. In contrast, at E12.5, co-application of KB and Ryanodine decreased HRs of right-sided SAN by 45.0±9.9% (n=9, close to additive blockade from both drugs: 14.2%+27.9%=42.1%) and did not cause any pacemaker switch. Thus, although both Ryanodine receptor and NCX1 contributed to automaticity of SAN at E12.5, they did not compensate each other when one is deficient. Additionally, substantial residual heart rhythm (55±9.9%, n=9) remains after the combined blockade by KB and Ryanodine at E12.5. These results indicate additional mechanisms, such as ion channels, may play increasing contributions to the automaticity at E12.5.
Furthermore, 10 µM 2-APB (2-aminoethoxydiphenyl borate), a membrane-permeable IP3R blocker,26 minimally decreased HRs of E8.5 (n=7) and E12.5 (n=7) hearts (8% and 5% respectively, Fig. 2A, C). However, after RyR blockade, 10 µM 2-APB decreased HRs of E8.5 (n=8) and E12.5 (n=7) hearts by an additional 15% and 11% respectively (Fig. 2C), indicating that IP3Rs also play a small part in automaticity reserve when Ryanodine receptors are deficient at both ages. Ratios of residual HRs relative to baseline HRs (mean ± SD in %) after drug blockade by these pharmacological agents during development are summarized in Figure 2C and its legend. The degree of HR blockade by two drugs in combination is always more than one drug alone (ANOVA, p<0.05) at either E8.5 or E12.5 (Fig. 2C). Each pairwise comparison is indicated by connecting lines with an asterisk denoting a p-value <0.05 (Fig. 2C). Of note, each drug used alone in this study did not produce similar effects on automaticity or field potential properties as other ion channel blockers or with respect to one another; rendering non-specific action from each drug unlikely.
Although [Ca2+]i signaling plays a major role in generating automaticity at the LIFT of E8.5 hearts, the switch of dominant pacemaker site to the right-sided SAN (Fig. 1G) with a concomitant increase in pacemaker rates at E12.5 (Supplemental Fig. 2) suggests that additional mechanisms of automaticity may evolve with embryonic development. For instance, hyperpolarization-activated, cyclic nucleotide-gated (HCN) cation channels, especially the type 4 isoform (HCN4), are thought to mediate pacemaker activity of mature SAN cells.4–5 ZD7288 at 1–10 µM is a specific HCN channel blocker and is used to inhibit HCN-mediated current (If) and automaticity.27 We found that 10 µM ZD7288 blocked HRs of E8.5 hearts by 27.4±10.8% (n=20) and significantly reduced HRs of E12.5 hearts by 66.0±11.0% (n=15, Fig. 3A–B). Furthermore, T-type Ca2+ channels have also been implicated in physiological and neurohormonal modulation of heart rates.28 Mibefradil at 1 µM29 and nickel at 40 µM30 are relatively specific blockers for T-type Ca2+ channels. 1 µM Mibefradil slightly decreased HRs of E8.5 hearts by 9.4±8.0% (n=12) and significantly decreased HRs of E12.5 SANs by 28.3±7.5% (n=15, Fig. 3A–B). Nickel (40 µM) produced similar patterns of blockade on HRs as 0.5 µM Mibefradil (Supplemental Fig. 6A). Results of residual HRs (relative to baseline HRs) after exposure to HCN and T-type Ca2+ channel blockers are summarized in Figure 3B. Pairwise comparisons (indicated by connecting lines) between E8.5 and E12.5 groups after each drug treatment (ANOVA, with all p-values <0.01) in Fig. 3B demonstrate that sarcolemmal ion channels play increasing contributions to automaticity at E12.5 than at E8.5.
Immunostaining with HCN4 antibodies revealed more positive staining for HCN4 channels in LIFT than RIFT regions of E8.5 hearts (Supplemental Fig. 6B–C). HCN4 immunolabeling increased at the right SAN region and was down-regulated significantly in left atria of E12.5 hearts (right panel in Supplemental Fig. 6B, and Supplemental Fig.6C). Immunostaining with CaV3.1 and CaV3.2 antibodies, two main subunits for T-type Ca2+ channels, showed minimal presence of membrane T-type Ca2+ channel subunits in either LIFT or RIFT regions of E8.5 hearts (Supplemental Fig. 6Da–c). The CaV3.1 and CaV3.2 staining was up-regulated in both atria of E12.5 hearts and some CaV3.1/CaV3.2 staining could also be observed at the right SAN region (Supplemental Fig. 6Dd–f). These results suggest that sarcolemmal ion channels play an increasing role in generating “faster” pacemaker rhythm from E8.5 to E12.5, reflecting ongoing differentiation and/or maturation. This mechanistic maturation of automaticity appears to accompany the location switch of dominant pacemaker sites (Fig. 1G) during early cardiogenesis from E8.5 to E12.5.
The relatively TTX-resistant cardiac isoform of sodium (Na+) channels, NaV1.5, and the TTX-sensitive neuronal isoform, NaV1.1, have been suggested to underlie pacemaking and conduction of adult mouse SAN.31 We used 2 µM TTX, which should block 100% of NaV1.1 and at least >50% of NaV1.5 channels, to determine the role of Na+ channels in mediating dominant automaticity and cardiac conduction during early cardiogenesis.31 2 µM TTX minimally decreased HRs of E8.5 hearts (4.4±6.1%, n=7) and slightly prolonged field potential durations (FPDs) (9%, p=NS; Fig. 3C–F). In E12.5 hearts, MEA recordings clearly demonstrated SAN signals (Figs. 1 and and3),3), and the effect of TTX on SAN automaticity could be accurately evaluated. 2 µM TTX slightly reduced HRs of sinus rhythm (6.5±5.9%) at E12.5 (n=14) but significantly prolonged sinus node field potential duration (SNd, p<0.01 with ANOVA), sinoatrial conduction interval (SAI, p<0.01), and atrial field potential duration (AFPD, p<0.01) (Fig. 3F). 2 µM TTX also led to various degrees of sinoatrial (SA) exit blocks (Supplemental Fig. 7A) and RA-LA conduction block (Supplemental Fig. 7B). Prolongation of field potential duration and conduction interval without affecting HR suggests that Na+ channels primarily contribute to electrical propagation and have minimal contribution to the dominant automaticity of E8.5 and E12.5 hearts. Values of residual HRs and various electrophysiological parameters (mean ± SD in %) after 2 µM TTX blockade are summarized in Figures 3D and F, respectively.
Most importantly, 2 µM TTX revealed a TTX-resistant sino-ventricular conducting pathway in the RA without any atrial activation in 7 of 14 E12.5 hearts (Fig. 4A). Previous reports by Rentschler and colleagues demonstrated the existence of functional AVN-bundle branch-Purkinje axis for atrioventricular conduction and ventricular activation with earliest double breakthroughs at ventricular apex of E12.5 mouse hearts.14 Color-coded activation maps (Fig. 4B) demonstrated that the sequence of ventricular activation (double breakthroughs with apical-to-base propagation pattern) are the same between normally conducted beats with atrial activation (the first beat in Fig. 4A) and sino-ventricular (Sino-V) conduction without atrial excitation (the third beat in Fig. 4A), suggesting a preferential pathway mediating Sino-V conduction in the RA of E12.5 hearts from the SAN to ventricles via the AVN-Purkinje axis (Fig. 4B and see Supplemental Movies 3–5). The unchanged morphology of ventricular FPs in all ventricular leads (two representative leads are shown in Fig. 4C) and lack of any ventricular pre-excitation during Sino-V conduction also support that this preferential pathway connects the SAN to AVN conduction axis. Furthermore, during a stable 2:1 Sino-V conduction rhythm (Fig. 5A–B), changes of ventricular intervals (V1V1) followed the changes of preceding SAN signal intervals (S1S1) in a 1:1 relationship during both normal (S-A-V) and Sino-V conducted beats (Fig. 5B–C); suggesting ventricular activations were guided by the preceding SAN rhythm (SAN-V linking). Also, SAN signals detected from electrodes surrounding the SAN area demonstrated a high-SAN-to-low-SAN (high-to-low) conducting sequence during normal and Sino-V conducted beats (Fig. 5D–E). Thus, high-to-low SAN conducting sequences with persistent SAN-V linking ruled out the possibility of coincidental junctional escape beats during the Sino-V conduction. Moreover, during TTX blockade, two of the E12.5 hearts displayed short periods of sinus arrhythmia with wide ranges of sinus cycle-length (S1S1) variations, which enabled the assessment of responses of Sino-V intervals (SVI) during Sino-V conduction to preceding sinus cycle length (S1S1). The Sino-V intervals demonstrated decremental AVN conducting properties (Fig. 5F),32 strongly supporting that a TTX-resistant internodal pathway mediates these Sino-V conducted beats.
The most plausible explanation for this TTX-resistant pathway is that the internodal pathway does not express nor relies on cardiac Na+ channel subunits (mainly NaV1.5) for signal propagation. Using NaV1.5 antibodies, we found that an internodal pathway with very low NaV1.5 expression indeed existed at the junction of RA and right sinus horn (rsh) at E12.5. Furthermore, with both HCN4 antibody labeling and in situ hybridization,33 we found that this internodal pathway completely coincided with cells stained positive for HCN4 channels in the RA of E12.5 hearts (Fig. 6A–B). Using different markers for the cardiac conduction system (such as Tbx3 and Podoplanin), other groups also demonstrated a similar immunocytochemical distinct pathway in E12.5 mouse hearts connecting SAN to the AVN.2,7 To the best of our knowledge, this study is the first electrophysiological demonstration of a functional and TTX-resistant pathway mediating internodal conduction in E12.5 mouse hearts.
We applied a MEA-based electrophysiological mapping system that enables rapid restoration of physiological HRs, thereby allowing accurate functional assessment of early pacemaker development and pharmacological investigations. No foreign dyes were used so as to avoid any non-physiological influence on the properties of dominant pacemaker cells. Importantly, a distinct SAN signal could be consistently recorded in E12.5 mouse hearts, which enables accurate determination of automaticity mechanisms in the SAN at E12.5. We find that pacemaker dominance in mice switches from the LIFT at E8.5 to the right-sided SAN region by E12.5. Intracellular Ca2+ handling proteins develop early in the LIFT at E8.5, which dominates the control of automaticity. At later developmental stages when heart rates increase, sarcolemmal ion channels evolve in the SAN to play increasing roles in faster impulse generation. Furthermore, this mechanistic maturation of SAN automaticity, combined with concomitant ion channel down-regulation at the LIFT, confers the pacemaker location switch and permits faster responses to rising electrophysiological and neurohormonal demands (summarized in Fig. 6C). With pharmacological manipulations, a clear pacemaker location switch could be demonstrated (Fig. 2B), suggesting the existence of two pacemaking regions as early as E8.5 with different degrees of maturation and mechanisms of impulse generation. Previously, Kamino’s group used a “rhythmicity gradient” to explain the switch in dominant pacemaker locations in early cardiogenesis, yet provided no mechanistic insight.12 In contrast, we show here that the pacemaker location switch during early cardiogenesis is the result of differential maturation of pacemaking mechanisms between the left and right inflow tract regions. Moreover, loss of left-sided dominance of automaticity coincides with down-regulation of HCN4 channels, suggesting inhibitory signals exist at LIFT area during cardiac development.8
Importantly, we find that compensatory mechanisms become operational as early as E8.5 to maintain heart rates and sustain lives when the dominant pacemaking mechanism is impaired. We show that a Ryanodine receptor or NCX blocker alone only slightly affects baseline heart rates at E8.5. However, small heart rate changes by a single blocker is likely due to compensatory mechanisms from additional [Ca2+]i handling proteins because co-application of two blockers abolishes left-sided pacemaker and causes the pacemaker location switch (Fig. 2B–C).
Furthermore, mechanisms underlying automaticity become increasingly complex with further cardiac development. Ion channels, such as HCN and T-type Ca2+ channels, play increasing contributions to dominant automaticity of the SAN when heart rate becomes faster at E12.5. It is not possible to report the contribution of every ion channels to automaticity development in a single study. Therefore, future studies are needed to see if other ion channels also display increasing contributions to automaticity with cardiac development. Voltage-dependent ion channels could alter membrane potential faster than [Ca2+]i handling machineries, and their increasing roles in mediating automaticity at E12.5 might be essential for developing faster heart rates to match rising metabolic demands. Nonetheless, we find that intracellular Ca2+-mediated automaticity remains functional at E12.5, suggesting many mechanisms work together to play the fail-safe task in maintaining automaticity (automaticity reserve) at later stages of cardiac development. Also, because most embryonic cardiomyocytes exhibit automaticity,17 our study indicates that simply evaluating heart rate changes by genetic manipulations or drug applications,3–5,22,23,28 without monitoring locations of dominant pacemakers, might lead to erroneous interpretations about mechanisms of dominant automaticity during early embryonic development.
Furthermore, using TTX, we find that Na+ channels play a minor role in automaticity and contribute mainly to electrical propagation of E8.5-E12.5 mouse hearts. Most importantly, TTX treatment at low micromolar concentrations blocked atrial activation and revealed a TTX-resistant internodal conducting pathway in E12.5 hearts. This preferential pathway is consistent with immunocytochemical evidence shown by others using podoplanin and Tbx3 as cardiac conduction system markers.2,7 Here, we provided the first functional evidence that such an internodal conduction pathway exists at the junction of RA/right sinus horn myocardium of E12.5 mouse hearts. However, this internodal pathway may not exist after birth since HCN4 immunostaining disappears at this internodal location in postnatal mouse hearts.33
Finally, this MEA-based recording technique is easy to perform, consistently enables the restoration of physiological heart rates, and facilitates pharmacological investigation of embryonic mouse hearts. This technique also reveals clear SAN signals for accurate electrophysiological studies and therefore could complement various optical imaging methods. In the future, this method may potentially be applied to study SAN development of mutant mice that survive to E12.5 but could not tolerate the optical dye-loading process.
Financial Support: Dr. Chen is supported by an Early Career Development Award from American College of Cardiology Foundation, and research grants from Bechtel Trusts & Foundation, California Institute of Regenerative Medicine (CIRM) (RS1-00171-1, RT1-01143, and RB2-01512), and NIH (HL105194)
Author contributionsT. Yi and H.-S.V. Chen designed research; T. Yi, J. Wen, C. Kim and H.-S.V. Chen performed research; J. Wong and W. Soussou designed the software for activation map; S. Sink, O. Taghli-Lamallem and E. Feller contributed significantly to data analysis/interpretation; M. Fink provided motion analysis software; W. Giles evaluated and commented on the manuscript; T. Yi, and H.-S.V.C. wrote the paper.