- Original Paper
- Open Access
Mitogen-activated protein kinase p38 modulates pacemaker ion channels differentiation in P19-derived pluripotent cells
The Journal of Physiological Sciences volume 70, Article number: 39 (2020)
Signal regulators during early cardiogenetic differentiation for the cellular automaticity are largely unknown. Our investigations were designed to clarify the role of transcription factors and their modulators in P19-derived cardiomyocytes to the expression of cardiac pacemaker ion channels. Transcription factors Csx/Nkx2.5 and GATA4 but not MEF2C were markedly inhibited by p38 MAP kinase inhibition in a distinct manner; expression but not phosphorylation of GATA4 was reduced by inhibition of p38 MAP kinase actions. In the presence of an ERK1/2,5 inhibitor PD98059 or a JNK MAP kinase inhibitor SP600125, P19 cells successfully differentiated into cardiomyocytes displaying spontaneous beatings with expression of three types of pacemaker ion channels. We demonstrate that acquisition of cellular automaticity and the expression of pacemaker ion channels are regulated by the transcription factors, Csx/Nkx2.5 and GATA4, through intracellular signals including p38 MAP kinase in the process of P19-derived pluripotent cells differentiation into cardiomyocytes.
The generation of pacemaker ionic channels is regarded as one of the most important factors responsible for the spontaneous beating of myocytes. To date, many studies have demonstrated several ionic mechanisms for the generation of slow diastolic depolarization in action potentials at the sino-atrial node in the heart [1, 2], including the L-type Ca2+ channel current (ICa.L), the T-type Ca2+ channel current (ICa.T) , the hyperpolarization-activated inward current (If) , the rapidly activating delayed rectifier K+ current (IKr) , and the sustained inward current (Ist) .
In mammalian embryonal differentiation, morphogenesis of the cardiovascular system is initiated in 20 gestational days in humans and 7.5 gestational days in mice, followed by the formation of a simple tubular heart and generation of spontaneous beating . The generation of specific contractile proteins and ion channel proteins in earlier cardiac differentiation is triggered by cardiac-specific gene expression. It has recently been shown that the expression of cardiac genes is regulated by several specific transcription factors during embryonal heart differentiation [8, 9]. Csx/Nkx2.5, GATA4, and MEF2C are considered to be cardiac-specific transcription factors that play the critical roles in the early development of the heart, and serve as useful molecular markers to identify cardiac inductive signals from other tissues or germ layers .
The homeobox gene Csx/Nkx2.5, a vertebrate homologue of Drosophila tinman, is one of the earliest known markers of mesoderm that give rise to myocardium, and its expression persists throughout the myocardium in the fully formed heart [7, 8]. GATA4 serves as a cardiac-specific member of the GATA family of zinc finger transcription factors which is detected very early in the cardiogenetic area and persists later in the developing heart . MEF2C belongs to the MEF2 family of MADS-box, and is expressed throughout the heart during mouse embryogenesis which controls cardiac morphogenesis and myogenesis [8, 12]. In the stimulation of cellular proliferation factors, including BMP, ET-1, and CT-1, the differentiation of cardiomyocytes is generated via the activation of intranuclear transcription factors such as Csx/Nkx2.5, GATA4, and MEF2C [7, 10]. However, in the early stage of embryonal development, signal transduction systems from the cellular membrane receptor signals into intranuclear transcription factors, subsequently activate the cardiac-specific proteins including ion channel expression, are poorly understood. On the other hand, the MAP kinase pathways, ERK1/2, ERK5, JNK, and p38 superfamilies, are major signaling systems by which the cells transfer extracellular signals into intracellular responses . It is known that MAP kinase performs the key roles in the differentiation of cardiomyocytes and in morphogenesis .
P19CL6 cells, clonal derivatives of mouse pluripotent P19 embryonal carcinoma cells, efficiently differentiate into spontaneously beating cardiomyocytes under adherent conditions when treated with 1% dimethyl sulfoxide (Me2SO) [10, 15]. It is widely accepted that P19 cells are excellent models for studying the development of cardiogenesis and cardiac gene expression in more than 100 studies [10, 15,16,17,18]. In this context, we focused on the molecular mechanisms of the expression of cardiac ionic channels, in order to assess the possible interaction of MAP kinase and specific cardiac transcription factors in P19-derived cardiomyocytes. We show that the expression of pacemaker ion channels are regulated by the transcription factors, Csx/Nkx2.5 and GATA4, through intracellular signals including p38 MAP kinase in the process of P19 cells differentiation into cardiomyocytes with automaticity.
Materials and methods
Cell culture and differentiation
P19CL6 cells were cultured with α-MEM containing 1% dimethyl sulfoxide (Me2SO) to induce differentiation into spontaneously beating cardiomyocytes as described previously [10, 15], with or without one of the following MAPK inhibitors: PD98059 (PD), an inhibitor of ERK1/2,5; SP600125 (SP), an inhibitor of JNK; SB203580 (SB), a specific inhibitor of the p38-pathway. For long-term inhibition of MAP kinase signals, 10 μM concentration of each inhibitor was applied [17, 18]. Days of differentiation were numbered consecutively after the first day of 1% Me2SO treatment, considered as day 0.
Single cardiomyocyte preparation
Spontaneously beating single cardiomyocytes were chosen for electrophysiological recordings by using of whole-cell current clamp and voltage clamp experiments [4, 19, 20]. All electrophysiological measurements were performed at 36 °C. Data were acquired at a sampling rate of 10 kHz, filtered at 2 kHz, and analyzed off-line by using the Pulsefit (HEKA) and Sigma9.0 software (SPSS Inc. Chicago, IL, USA).
Solutions and chemicals
For measuring membrane potentials, the bath solution (normal Tyrode solution) contained (in mM) NaCl 140, KCl 5.4, CaCl2 1.8, MgCl2 1, HEPES 10, and glucose 10, with the pH adjusted to 7.4 with NaOH. The pipette solution contained (in mM) KCl 140, MgCl2 2, Mg-ATP 5, creatine phosphate (disodium salt) 5, EGTA 0.05, and HEPES 10, with the pH adjusted to 7.2 with KOH. For measuring ICa, the bath solution contained the normal Tyrode solution supplemented with 30 μM tetrodotoxin to eliminate Na+ current. The pipette solution contained (in mM) CsCl 130, MgCl2 2, Mg-ATP 2, Na2-GTP 2, EGTA 10, and HEPES 5, with the pH adjusted to 7.2 with CsOH. For measuring If, the bath solution contained the normal Tyrode solution supplemented with (in mM) BaCl2 1, CdCl2 0.4, tetrodotoxin 0.03, 4,4′-diisothiocyanostilbene-2,2′-disulfonic acid (DIDS) 0.1, and 4-aminopyridine (4-AP) 1, to eliminate inward rectifier K+ current, Ca2+ current (T- and L-type), Na+ current, and transient outward K+ current (Ito), respectively. The pipette solution contained (in mM) KCl 140, MgCl2 1, EGTA 10, HEPES 5, and Mg-ATP 5, with the pH adjusted to 7.2 with KOH. PD98059, SP600125, and SB203580 were purchased from Calbiochem (La Jolla, CA, USA). All other chemicals were purchased from Wako Pure Chemical Industries (Osaka, Japan).
Quantitative real-time RT-PCR
Total RNA was extracted from P19-derived cardiomyocytes using Isogen (Nippongene, Tokyo, Japan). The cDNA was synthesized from 1 μg of total RNA using Transcriptor First Strand cDNA Synthesis Kit (Roche Molecular System Inc., Alameda, CA, USA). The real-time PCR was performed on Light Cycler (Roche) using the FastStart DNA Master SYBR Green I (Roche) as a detection reagent. For Csx/Nkx2.5, primer F sequenced 5′-CAG TGG AGC TGG ACA AAG CC-3′, and primer R sequenced 5′-TAG CGA CGG TTC TGG AAC CA-3′ were used. For GATA4, primer F sequenced 5′-GCC TGT ATG TAA TGC CTG CG-3′, and primer R sequenced 5′-CCG AGC AGG AAT TTG AAG AGG -3′ were used. For MEF2C, primer F sequenced 5′-GTA TGT CTC CTG GTG TAA CA-3′, and primer R sequenced 5′-GGA TAT CCT CCC ATT CCT TG-3′ were used. As an internal control in each reaction, mouse GAPDH primer F of 5′-CCA AGG TCA TCC ATG ACA AC-3′, and primer R of 5′-TTA CTC CTT GGA GGC CAT GT-3′ were used. The expression of each target mRNA relative to GAPDH under experimental and control conditions was calculated based on the threshold cycle (CT) as r = 2−Δ(ΔCT), where ΔCT = CT target − CT GAPDH and Δ(ΔCT) = ΔCT experimental − ΔCT control. Gel electrophoresis was performed to confirm the correct size of amplification and the absence of unspecific bands.
Western blot analysis
P19-derived cardiomyocytes were lysed in cold cell lysis buffer. The lysates were clarified by centrifugation at 12 000g for 15 min at 4 °C. The supernatants were used for phospho-p38 and total p38 Western blot analysis.
Analysis of data
All data are presented as mean ± S.E. One-way ANOVA followed by a Bonferroni post hoc test was used for multiple comparisons. Differences were considered significant when p values were less than 0.05.
Action potentials in differentiating P19-derived cells
Spontaneous beating activity is thought of as the unique character of the cardiomyocyte, distinct from other cells. We investigated mechanisms for the spontaneous beating activity of P19-derived cardiomyocytes by means of electrical excitability or action potentials. P19-derived cardiomyocytes displayed spontaneous beating, which was initialized on differentiated day 10 with very slow frequency. The beating rates reached 89.2 ± 17.4 bpm (n = 15) on days 14 to 17 (Table 1). Action potentials (APs) measured by current-clamp mode indicated a maximum diastolic potential (MDP) of − 58.6 ± 2.5 mV (n = 15) on days 14 to 17 (Fig. 1a). On the other hand, in the absence of Me2SO, P19 cells did not differentiate into cardiomyocytes, thus demonstrating no spontaneous beating activity (Table 1). The membrane potentials (MP) were − 52.8 ± 8.5 mV (n = 11) on differentiated day 0 (data not shown) and − 51.9 ± 6.7 mV (n = 14) on days 14 to 17.
To clarify the intracellular signals in the cardiogenetic process, MAP kinase-dependent signal transduction systems were explored in accordance with the acquisition of cellular excitability. We first examined the possible contribution of ERK1/2- and ERK5-dependent cascades to the development of action potential configurations. PD98059 (PD) has been widely used as a selective inhibitor of ERK1/2,5, which belong to MAP kinase kinase superfamilies. We recorded the action potentials of P19-derived cells on differentiation days 14 to 17, the cells of which were treated with PD throughout the cardiogenetic process from day 0. Similar to the control condition (1% Me2SO without MAPK kinase inhibitor), P19-derived cardiomyocytes displayed spontaneous beating of 83 ± 12 bpm, MDP of − 57.9 ± 1.7 mV (Fig. 1b), APD50 of 105 ± 9 ms, and APD75 of 299 ± 19 ms in the presence of PD (Table 1). To study the possible interaction of the JNK pathway, another branch of MAP kinase cascade, we used SP600125 (SP) in order to selectively abolish the JNK MAP kinase-dependent signal cascade. We found that action potentials were distinctly formed (Fig. 1c), and their parameters of APs were similar to those observed under control conditions (Table 1). In the presence of SB203580 (SB), an inhibitor of p38, P19 cells were found to be quiescent without spontaneous beating during the observation period of 30 days. The measured MDP (− 49.0 ± 3.3 mV, n = 11) was statistically identical to that of undifferentiated P19 cells (− 51.9 ± 6.7 mV, n = 14) on days 14 to 17 (Fig. 1d and Table 1). To identify the detailed electrophysiological development of P19-derived cells in the presence of PD, SP, or SB, we initially focused on parameters of membrane potentials: APD50, APD75, and MDP (Table 1). The APD50 in Me2SO + PD-treated cells were significantly shorter, and maximum diastolic potentials in Me2SO + SB treated cells were significantly depolarized in comparison with that in Me2SO treated cells. Other parameters in PD and SP groups were all statistically identical to those in Me2SO-treated cells.
Electrophysiological characters of pacemaker channel currents in differentiated P19 cells
We investigated the development of membrane currents focusing on pacemaker channel currents such as L-type Ca2+ current (ICa.L), T-type Ca2+ current (ICa.T), and hyperpolarization-activated inward current (If). Expression of three types of pacemaker currents was identified in the P19-derived cardiomyocytes on differentiation days 14 to 17 (Fig. 2a). However, under the differentiating condition from day 0 to day 7, these three types of ion channels were not detected (data not shown), which was consistent with quiescence of the cells on those days. In the presence of PD and SP, ICa.L, ICa.T and If were all distinctly formed in function (Fig. 2b, c), which were nearly identical to those under control conditions (Fig. 2a). In the presence of SB, these ionic currents were not expressed in function (Fig. 2d). Subsequently, we investigated ICa.L in terms of current density and voltage dependency. All of the ICa.L values recorded, except those in the presence of SB, demonstrated a similar current (I)–voltage (V) relationship (Fig. 2e). The maximum peak currents were obtained at the test potential of 0 mV. The maximum conductance for the L-type Ca2+ channel (Gmax, Ca.L) were 0.96 ± 0.05 nS pF−1 (n = 19) under control conditions, 1.52 ± 0.06 nS pF−1 (n = 17) in the presence of PD, and 0.95 ± 0.06 nS pF−1 (n = 11) in the presence of SP (Fig. 2h), whereas ICa.L was completely abolished by SB (Fig. 2e, h, Table 2). In the presence of PD, ICa.L and Gmax, Ca.L were significantly larger than those in control by unknown reasons.
The I–V relationship of ICa.T indicated that the maximum conductance for the T-type Ca2+ channel (Gmax, Ca.T) obtained at − 30 mV was 0.64 ± 0.04 nS pF−1 (n = 17) under control conditions, and was not affected by the presence of PD or SP, that is to say, Gmax, Ca.T and the maximum peak current of ICa.T under the effect of PD or SP were highly similar to those under control conditions, whereas ICa.T was completely abolished by SB (Fig. 2f, i, Table 2).
If densities at − 120 mV were − 26.9 ± 4.7 pA pF−1 (n = 37) under control conditions, − 28.5 ± 3.9 pA pF−1 (n = 16) in the presence of PD, and − 37.4 ± 6.3 pA pF−1 (n = 10) in the presence of SP, indicating no statistically different significance between control group and PD group. Conductance curves for If indicate that half-activated potentials (V1/2) were − 94.5 ± 3.4 mV under control conditions, − 90.1 ± 5.6 mV in the presence of PD, and − 89.4 ± 5.3 mV in the presence of SP (Fig. 2g). If and Gmax,f were significantly larger in SP group (Fig. 2g, j). However, in the presence of SB, If displayed substantially no inward currents like to those of ICa.L and ICa.T. Interestingly, under the effect of SB, P19 cells were retarded in development in terms of membrane capacitance (Table 2), indicating a significant role of p38 in cell differentiation and growth.
Transcription factor changes in differentiated P19 cells
In order to identify the signals to initiate ion channel expression in terms of transcription factor activation, we performed RT-PCR analyses for the detection of the mRNA encoding cardiac-specific transcription factors Csx/Nkx2.5, GATA4, and MEF2C which may be responsible for the expression of ion channels. Three cardiac-specific transcription factors appeared as early as day 1 in the differentiation stage, and their intensive expression were successively detected in the presence of Me2SO (Fig. 3a). Importantly, expression of Csx/Nkx2.5, GATA4 and MEF2C were quantitatively identical in amount and nearly stable from day 4 to day 14 (Fig. 3b, c).
We also found that under the effect of an ERK1/2,5-inhibitor PD, the expression of the three specific transcription factors was appreciably increased, with the same time course (Fig. 3d–f) as that under control conditions (Fig. 3a–c). Simultaneously, expression of cardiac-specific transcription factors under the effect of a JNK MAP kinase inhibitor SP was almost identical to those under control conditions in terms of time course and total RNA molecules (Fig. 3g–i). Therefore, JNK MAPK kinase was identified as not modifying the cardiogenetic process in terms of cardiac pacemaker ion channels expression. Csx/Nkx2.5 mRNA was slowly and suppressively expressed by a p38 inhibitor SB. Importantly, GATA4 mRNA was declined on differentiation day 4 or later by SB (Fig. 3l). On the other hand, MEF2C mRNA expression was unaffected by a p38 inhibitor SB. These results indicate that the expression of cardiac pacemaker ion channels and the beating property of myocytes are distinctly controlled by specific transcription factors, Csx/Nkx2.5 and GATA4, activated by the p38 MAP kinase pathway.
Expression of transcription factors were quantitatively compared in the presence of MAP kinase inhibitors on the differentiation day 14 (Fig. 4). Expression of Csx/Nkx2.5 was halved and expression of GATA4 was highly inhibited in the presence of SB (Fig. 4a, b). On the other hand, the expression levels of Csx/Nkx2.5, GATA4 and MEF2C were all significantly increased by PD (Fig. 4a–c). Taken together, p38 MAP kinase-dependent signal transduction and classical MAP kinase 1/2,5 signal transduction contrarily modulate transcription factors Csx/Nkx2.5 and GATA4 in the process of cardiogenetic differentiation.
P38 protein changes in differentiated P19 cells
In order to identify roles of p38 and phospho-p38 proteins in cell differentiation, we performed Western blot analysis. Although total p38 protein was unchanged, phospho-p38 protein levels, or activated form of p38, were increased 1.8-fold on the differentiation day 4 in the presence of Me2SO (Fig. 5a). Importantly, phospho-p38 expression was quantitatively unchanged, and nearly stable from day 4 to day 14, which are comparable to changes of the transcription factors (Fig. 3c). On the differentiation day 14, phospho-p38 expression was highly inhibited in the presence of an inhibitor of p38 phosphorylation SB, although it was unaffected by PD and SP (Fig. 5b).
p38-dependent expression of Csx/Nkx2.5 and GATA4
Unlike Csx/Nkx2.5, GATA4 is regulated to enhance its transcriptional potency by the MAP kinase signal cascade through direct phosphorylation. To further confirm and specify the role of p38 MAP kinase, protein expression of Csx/Nkx2.5 and GATA4 was examined with or without the action of p38 MAP kinase. Expression of Csx/Nkx2.5 and GATA4 proteins were significantly reduced by a p38 MAP kinase inhibitor SB in a concentration-dependent manner (Fig. 6a, b). Although expression of phospho-GATA4 was similarly inhibited by SB, decreases of phospho-GATA4 by SB were identical to those in total GATA4 as indicated by phospho-GATA4/GATA4 ratio in Fig. 6d, which suggests that p38 MAP kinase-dependent phosphorylation acts not to activate GATA4 transcriptional potency but to increase GATA4 transcription.
The electrical and mechanical mechanisms governing the precise and highly organized actions in the heart are extremely complex, requiring coordinated neural and humoral factors in the healthy and pathological conditions [21,22,23,24,25,26,27,28,29,30,31]. Such sophisticated coordination also depends on the developmental changes in the steps of cardiogenesis including the generation of pacemaker potentials [1,2,3,4,5,6, 32,33,34]. Our present study suggests the possible mechanism of the intracellular signal pathways that regulates expression of cardiac-specific transcription factors responsible for the generation of pacemaker potentials/ion channels in P19-derived cardiomyocytes.
Actions of the MAP kinase pathway on differentiation of cardiomyocytes
During the cardiogenetic differentiation of P19, three distinct MAP kinases, namely, JNK, p38, and ERK1/2, are coordinately activated, and contribute to both proliferation and differentiation . Our results demonstrate that P19-derived cells exert cardiomyocyte-like electrophysiological features of spontaneous beating activity with distinct expression of pacemaker ion channels under the inhibition of either ERK1/2,5 or JNK MAP kinase. These observations are consistent with a previous study demonstrating a dispensable role of “classical” MAP kinase pathway of ERK1/2 on cardiomyocyte differentiation in mouse ES cells . Based on the current study, inhibition of p38 MAP kinase absolutely disrupted the cardiogenetic process of P19 cells in the following ways: (1) spontaneous beating activity was absent during the observation period of 30 days; (2) cellular membrane potentials were identical to those of undifferentiated P19 cells; and (3) current density of ICa.L, ICa.T and If was substantially nil. Therefore, our results imply that p38 MAP kinase plays a pivotal role as the intracellular signal mediator generating pacemaker ionic channels during the process of P19-derived cardiomyocyte differentiation.
The possible relationships between Csx/Nkx2.5 and GATA4
Several critical transcription factors, including Csx/Nkx2.5, GATA4 and MEF2C, have been shown to regulate cardiac-specific genes, with no single transcription factor responsible for the differentiation of lateral plate mesoderm into cardiac cells . Recent studies have indicated the functional cooperation of Csx/Nkx2.5 and GATA4 for cardiogenetic processes [10, 14]. Based on our experiments, the interaction of Csx/Nkx2.5 with GATA4 in cardiogenesis is likely due to four possible mechanisms:
Csx/Nkx2.5-dependent activation of GATA4. Several studies have elucidated the critical role of Csx/Nkx2.5 during specific cardiac differentiation and have identified direct downstream targets for Csx/Nkx2.5, namely, ANP, cardiac α-actin, MHC, and others [9, 10]. The expression of Csx/Nkx2.5 was detected in cardiac crescent formation at embryonic day 7.5 (E7.5), prior to the expression of GATA4 at E8 in the formation of the linear heart tubes . Under the inhibition of p38 MAP kinase, expression of Csx/Nkx2.5 was slightly maintained on day 4 or later, while expression of GATA4 was highly eliminated (Fig. 3). Importantly, Csx/Nkx2.5, as a cofactor of GATA4, has been shown to be able to recruit GATA4 to cardiac gene promoters .
GATA4-dependent activation of Csx/Nkx2.5. It is notable that the extracellular signal of BMP is transformed via the interaction of Smad1/4 and GATA4 into Csx/Nkx2.5 activation . Under the inhibition of p38 MAP kinase, the expression of GATA4 preceded Csx/Nkx2.5 on day 2 (Fig. 3k), indicating that GATA4 may activate the expression of Csx/Nkx2.5. Our results disclosed that on differentiation day 14, GATA4 was nearly completely inhibited while expression of Csx/Nkx2.5 persisted. It is therefore suggested that Csx/Nkx2.5 activation by GATA4 could be sustained for a short period of time.
Synergistic interaction of GATA4 and Csx/Nkx2.5. Several recent studies have demonstrated that both homeodomain of Csx/Nkx2.5 and DNA-binding domain of MEF2C served as cofactors that interact with the C-terminal domain of GATA4, generating a complex for synergistic transactivation of genes . The most likely correct speculation is that expression of GATA4 plays a crucial role during differentiation, and that p38 MAP kinase predominantly modulates expression of GATA4, while the expression of Csx/Nkx2.5 is stimulated by dual pathways of GATA4 and other intracellular factors.
Independent interaction of GATA4 with Csx/Nkx2.5. The simplest explanation of our results is that GATA4 and Csx/Nkx2.5 independently regulate downstream targets and that both factors are required to activate all related cardiac ion channel genes expression needed for differentiation and cell beating. The expression of GATA4 exhibited high sensitivity to a p38 MAP kinase inhibitor, while the expression of Csx/Nkx2.5 was mediated, in a part, by the p38 MAP kinase pathway, possibly regulated by other intracellular factors such as SRF and others .
MAPK unaffected MEF2C activation
The expression of MEF2C is regulated by many intracellular and intranuclear transcription factors. GATA4 may physically interact with MEF2C, thus regulating cardiac gene expression . In P19CL6noggin cells, the overexpression of Csx/Nkx2.5 induced the expression of MEF2C, which suggests that expression of MEF2C was positively regulated by Csx/Nkx2.5 . However, our current study could be considered partly inconsistent with these previous findings in that the expression of MEF2C was almost identical to those under control conditions utilizing the p38 MAP kinase inhibitor, whereas Csx/Nkx2.5 and GATA4 were notably abolished. During cardiogenesis, other intracellular pathways such as PI3-kinase pathway  could be a candidate for the MEF2C activator.
p38 MAP kinase modulates the expression of GATA4
Several studies have disclosed that protein kinases, including ERK1/2  and p38 MAP kinase , catalyze Ser-105-specific phosphorylation of GATA4 in cardiomyocytes. Based on our findings in Fig. 3e, ERK1/2,5 and p38 MAP kinase may be distinct in function in terms of cardiogenetic differentiation in P19 cells. Consistent with our study, in the context of rat BNP promoter, p38, but not other ERK-dependent transcriptions, is responsible for GATA4 binding site activation in the promoter region, suggesting that the activation of p38 MAP kinase is necessary for ET1-induced GATA4 activation that leads to BNP gene expression in neonatal rat cardiomyocytes . Interestingly, a recent study demonstrated that p38 kinase was moderately activated by the ERK1/2,5 inhibitor PD in neonatal rat cardiomyocytes . Consistent with these studies, our work demonstrates that the expression of GATA4 and Csx/Nkx2.5 was augmented under the influence of ERK1/2,5 inhibition by PD (Fig. 4), supporting our hypothesis that p38 MAP kinase-dependent activation of GATA4 and Csx/Nkx2.5 is crucial for cardiogenesis in P19 cells. This is supportive of our electrophysiological determination that ICa.L and Gmax, Ca.L with PD were larger than those under control conditions (Fig. 2e, h). On the other hand, inhibition of the JNK pathway had no effect on ET1-induced GATA4 activation . Moreover, overexpression of MEKK1, an upstream kinase of JNK, independently induced BNP promoter activation on GATA4 actions . Consistently, treatment of P19 with a JNK inhibitor had no influence on the differentiation or expression of pacemaker ion channels in this study.
Although we have successfully identified that activation of transcription factors Csx/Nkx2.5 and GATA4 is a trigger for pacemaker ion channels expression including Cav3.2-T-type Ca2+ channel regulation , we need further investigation to clarify these possible association with other transcription factors such as Tbx3 and Tbx18, which reportedly form sino-atrial node in the early cardiogenesis . Although the generation of cardiac pacemaker potentials relies on a complex interplay between different types of currents carried by ion channels/transporters, we only focused on the three major ion channels on this study. Moreover, we were unable to identify molecules responsible for the automaticity in P19-derived cardiomyocytes. In cardiac pacemaker cells, the spontaneous action potential and resulting contraction is brought about through time- and voltage-dependent changes in the gating of ionic channels. The alternative activation of inward and outward current systems produces a feedback loop which generates the oscillation of the membrane potential. A large number of ionic currents in single sinus node cells have been identified, and, using the conductance and gating parameters of the individual current components, it is possible to reproduce the pacemaker activity by computer simulation. For instance, an example of the pacemaker model simulation, in which the action potential and ionic currents in sinus node cells were reproduced using the Kyoto Model developed by Sarai et al. . According to their study, more than 13 ionic currents have been incorporated. The currents include L-(ICaL) and T-type Ca2+ currents (ICaT), the delayed rectifier K+ current (IK), Na+-K+ pump current (INaK), muscarinic K+ current (IK,ACh), hyperpolarization-activated cation current (If), background cation current and a dihydropyridine-sensitive sustained inward current (Ist), and others. It is generally accepted that the rising phase of the action potential is due to ICaL and that the repolarization is caused by the inactivation of ICaL and the activation of IK. As the membrane hyperpolarizes, IK begins to deactivate, which in turn triggers a gradual depolarization after the maximum diastolic potential. During the later part of the pacemaker potential, several inward currents such as If, ICaT and Ist participate in driving the membrane towards the threshold of the following action potential. Among them, Cav1.3-L-type Ca2+ channel and HCN2/4-If channel could be major contributors [46,47,48], although we were unable to confirm their molecules in this P19-derived cardiomyocyte. These are the description on the classical “membrane clock”. Recent experimental evidence based on confocal cell imaging and supported by numerical modeling, however, shows that initiation of cardiac impulse is a more complex phenomenon and involves yet another clock that resides under the sarcolemma . This clock is the sarcoplasmic reticulum (SR): it generates spontaneous, but precisely timed, rhythmic, submembrane, local Ca2+ releases that appear not at the maximum diastolic potentials but during the late diastolic depolarization. We could not confirm molecules for the Ca2+ clock either in this investigation. Nevertheless, transcriptional regulation of other ion channels/transporters such as K+ channels, Na+ channels, Cl− channels and Na+–Ca2+ exchanger by Csx/Nkx2.5 and GATA4 needs to be elucidated. In addition, pacemaker rhythm in this cardiogenetic model is much slower than the beating ratio of mouse sino-atrial node, which requires further investigations to identify molecules needed to tune and mature pacemaking cardiomyocytes.
Our results demonstrate that the transcription factors Csx/Nkx2.5 and GATA4 association with p38 MAP kinase but no other MAP kinases regulates transcription of pacemaker ion channels contributing to spontaneous beating activity during the process of cardiac myocyte differentiation.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available in the Oita University School of Medicine repository. They are also available from the corresponding author on reasonable request.
- If :
Hyperpolarization-activated inward current
- ICa.L :
L-type Ca2+ channel current
- ICa.T :
T-type Ca2+ current
Extracellular signal-regulated kinase
- IKr :
Rapidly activating delayed rectifier K+ current
- Ist :
Sustained inward current
- Ito :
Transient outward K+ current
Beat per minute
Maximum diastolic potentials
- APD50 :
The duration of 50% repolarization of action potentials
- APD75 :
The duration of 75% repolarization of action potentials
- V1/2 :
- Cm :
Cellular membrane capacitance
- Gmax :
The maximum chord conductance
Irisawa H, Brown HF, Giles W (1993) Cardiac pacemaking in the sinoatrial node. Physiol Rev 73:197–227
Cho HS, Takano M, Noma A (2003) The electrophysiological properties of spontaneously beating pacemaker cells isolated from mouse sinoatrial node. J Physiol 550:169–180
Hagiwara N, Irisawa H, Kameyama M (1988) Contribution of two types of calcium currents to the pacemaker potentials of rabbit sino-atrial node cells. J Physiol 395:233–253
DiFrancesco D (1991) The contribution of the “pacemaker” current (If) to generation of spontaneous activity in rabbit sino-atrial node myocytes. J Physiol 434:23–40
Boyett MR, Honjo H, Kodama I (2000) The sinoatrial node, a heterogeneous pacemaker structure. Cardiovasc Res 47:658–687
Mitsuiye T, Shinagawa Y, Noma A (2000) Sustained inward current during pacemaker depolarization in mammalian sinoatrial node cells. Circ Res 87:88–91
Kirby ML (2002) Molecular embryogenesis of the heart. Pediatr Dev Pathol 5:516–543
Bruneau BG (2002) Transcriptional regulation of vertebrate cardiac morphogenesis. Circ Res 90:509–519
Akazawa H, Komuro I (2003) Roles of cardiac transcription factors in cardiac hypertrophy. Circ Res 92:1079–1088
Monzen K, Shiojima I, Hiroi Y, Kudoh S, Oka T, Takimoto E, Hayashi D, Hosoda T, Habara-Ohkubo A, Nakaoka T, Fujita T, Yazaki Y, Komuro I (1999) Bone morphogenetic proteins induce cardiomyocyte differentiation through the mitogen-activated protein kinase kinase kinase TAK1 and cardiac transcription factors Csx/Nkx2.5 and GATA-4. Mol Cell Biol 19:7096–7105
Pikkarainen S, Tokola H, Kerkela R, Ruskoaho H (2004) GATA transcription factors in the developing and adult heart. Cardiovasc Res 63:196–207
Lin Q, Schwarz J, Bucana C, Olson EN (1997) Control of mouse cardiac morphogenesis and myogenesis by transcription factor MEF2C. Science 276:1404–1407
Chang L, Karin M (2001) Mammalian MAP kinase signalling cascades. Nature 410:37–40
Eriksson M, Leppa S (2002) Mitogen-activated protein kinases and activator protein 1 are required for proliferation and cardiomyocyte differentiation of P19 embryonal carcinoma cells. J Biol Chem 277:15992–16001
Habara-Ohkubo A (1999) Differentiation of beating cardiac muscle cells from a derivative of P19 embryonal carcinoma cells. Cell Struct Funct 21:101–110
van der Heyden MA, Defize LH (2003) Twenty one years of P19 cells: what an embryonal carcinoma cell line taught us about cardiomyocyte differentiation. Cardiovasc Res 58:292–302
Yan L, Tang Q, Shen D, Peng S, Zheng Q, Guo H, Jiang M, Deng W (2008) SOCS-1 inhibits TNF-α-induced cardiomyocyte apoptosis via ERK1/2 pathway activation. Inflammation 31:180–188
Yin F, Wang YY, Du JH, Li C, Lu ZZ, Han C, Zhang YY (2006) Noncanonical cAMP pathway and p38 MAPK mediate β2-adrenergic receptor-induced IL-6 production in neonatal mouse cardiac fibroblasts. J Mol Cell Cardiol 40:384–393
Masuda K, Takanari H, Morishima M, Ma F, Wang Y, Takahashi N, Ono K (2018) Testosterone-mediated upregulation of delayed rectifier potassium channel in cardiomyocytes causes abbreviation of QT intervals in rats. J Physiol Sci 68:759–767
Liu Y, Wang P, Ma F, Zheng M, Liu G, Kume S, Kurokawa T, Ono K (2019) Asparagine-linked glycosylation modifies voltage-dependent gating properties of CaV3.1-T-type Ca2 + channel. J Physiol Sci 69:335–343
Garza MA, Wason EA, Cruger JR, Chung E, Zhang JQ (2019) Strength training attenuates post-infarct cardiac dysfunction and remodeling. J Physiol Sci 69:523–530
Lookin O, Protsenko Y (2019) The lack of slow force response in failing rat myocardium: role of stretch-induced modulation of Ca-TnC kinetics. J Physiol Sci 69:345–357
Li M, Zheng C, Kawada T, Inagaki M, Uemura K, Sugimachi M (2019) Chronic vagal nerve stimulation exerts additional beneficial effects on the beta-blocker-treated failing heart. J Physiol Sci 69:295–303
Tashiro M, Inoue H, Konishi M (2019) Modulation of Mg2 + influx and cytoplasmic free Mg2 + concentration in rat ventricular myocytes. J Physiol Sci 69:97–102
Balakina-Vikulova NA, Panfilov A, Solovyova O, Katsnelson LB (2020) Mechano-calcium and mechano-electric feedbacks in the human cardiomyocyte analyzed in a mathematical model. J Physiol Sci 70(1):12. https://doi.org/10.1186/s12576-020-00741-6
Yuan JQ, Yuan Y, Pan SS, Cai K (2020) Altered expression levels of autophagy-associated proteins during exercise preconditioning indicate the involvement of autophagy in cardioprotection against exercise-induced myocardial injury. J Physiol Sci 70(1):10. https://doi.org/10.1186/s12576-020-00738-1
Obata K, Morita H, Takaki M (2020) Mechanism underlying the negative inotropic effect in rat left ventricle in hyperthermia: the role of TRPV1. J Physiol Sci 70(1):4. https://doi.org/10.1186/s12576-020-00734-5
Kawai E, Takeda R, Ota A, Morita E, Imai D, Suzuki Y, Yokoyama H, Ueda SY, Nakahara H, Miyamoto T, Okazaki K (2020) Increase in diastolic blood pressure induced by fragrance inhalation of grapefruit essential oil is positively correlated with muscle sympathetic nerve activity. J Physiol Sci 70(1):2. https://doi.org/10.1186/s12576-020-00733-6
Fu Q, Ogoh S (2019) Sex differences in baroreflex function in health and disease. J Physiol Sci 69(6):851–859. https://doi.org/10.1007/s12576-019-00727-z
Nakahara H, Kawada T, Ueda SY, Kawai E, Yamamoto H, Sugimachi M, Miyamoto T (2019) Acupoint dependence of depressor and bradycardic responses elicited by manual acupuncture stimulation in humans. J Physiol Sci 69(6):1077–1084. https://doi.org/10.1007/s12576-019-00728-y
Koba S, Hisatome I, Watanabe T (2019) Augmented fear bradycardia in rats with heart failure. J Physiol Sci 69(6):875–883. https://doi.org/10.1007/s12576-019-00698-1
Oshita K, Kozasa Y, Nakagawa Y, Kuwabara Y, Kuwahara K, Nakagawa T, Nakashima N, Hiraki T, Takano M (2019) Overexpression of the HCN2 channel increases the arrhythmogenicity induced by hypokalemia. J Physiol Sci 69:653–660
Kawada T, Shimizu S, Uemura K, Hayama Y, Yamamoto H, Shishido T, Nishikawa T, Sugimachi M (2019) Ivabradine preserves dynamic sympathetic control of heart rate despite inducing significant bradycardia in rats. J Physiol Sci 69:211–222
Takagi D, Okamoto Y, Ohba T, Yamamoto H, Ono K (2020) Comparative study of hyperpolarization-activated currents in pulmonary vein cardiomyocytes isolated from rat, guinea pig, and rabbit. J Physiol Sci 70(1):6. https://doi.org/10.1186/s12576-020-00736-3
Davidson SM, Morange M (2000) Hsp25 and the p38 MAPK pathway are involved in differentiation of cardiomyocytes. Dev Biol 218:146–160
Sepulveda JL, Vlahopoulos S, Iyer D, Belaguli N, Schwartz RJ (2002) Combinatorial expression of GATA4, Nkx2-5, and serum response factor directs early cardiac gene activity. J Biol Chem 277:25775–25782
Brown CO 3rd, Chi X, Garcia-Gras E, Shirai M, Feng XH, Schwartz RJ (2004) The cardiac determination factor, Nkx2-5, is activated by mutual cofactors GATA-4 and Smad1/4 via a novel upstream enhancer. J Biol Chem 279:10659–10669
Dodou E, Verzi MP, Anderson JP, Xu SM, Black BL (2004) Mef2c is a direct transcriptional target of ISL1 and GATA factors in the anterior heart field during mouse embryonic development. Developmen 131:3931–3942
Liang Q, Wiese RJ, Bueno OF, Dai YS, Markham BE, Molkentin JD (2001) The transcription factor GATA4 is activated by extracellular signal-regulated kinase 1- and 2-mediated phosphorylation of serine 105 in cardiomyocytes. Mol Cell Biol 21:7460–7469
Charron F, Tsimiklis G, Arcand M, Robitaille L, Liang Q, Molkentin JD, Meloche S, Nemer M (2001) Tissue-specific GATA factors are transcriptional effectors of the small GTPase RhoA. Genes Dev 15:2702–2719
Kerkela R, Pikkarainen S, Majalahti-Palviainen T, Tokola H, Ruskoaho H (2002) Distinct roles of mitogen-activated protein kinase pathways in GATA-4 transcription factor-mediated regulation of B-type natriuretic peptide gene. J Biol Chem 277:13752–13760
Clerk A, Michael A, Sugden PH (1998) Stimulation of the p38 mitogen-activated protein kinase pathway in neonatal rat ventricular myocytes by the G protein-coupled receptor agonists, endothelin-1 and phenylephrine: a role in cardiac myocyte hypertrophy? J Cell Biol 142:523–535
Wang Y, Morishima M, Zheng M, Uchino T, Mannen K, Takahashi A, Nakaya Y, Komuro I, Ono K (2007) Transcription factors Csx/Nkx2.5 and GATA4 distinctly regulate expression of Ca2 + channels in neonatal rat heart. J Mol Cell Cardiol 42:1045–1053
Hoogaars WM, Engel A, Brons JF, Verkerk AO, de Lange FJ, Wong LY, Bakker ML, Clout DE, Wakker V, Barnett P, Ravesloot JH, Moorman AF, Verheijck EE, Christoffels VM (2007) Tbx3 controls the sinoatrial node gene program and imposes pacemaker function on the atria. Genes Dev 21:1098–1112
Sarai N, Matsuoka S, Kuratomi S, Ono K, Noma A (2003) Role of individual ionic current systems in the SA node hypothesized by a model study. Jpn J Physiol 53:125–134
Mangoni ME, Traboulsie A, Leoni AL, Couette B, Marger L, Le Quang K, Kupfer E, Cohen-Solal A, Vilar J, Shin HS, Escande D, Charpentier F, Nargeot J, Lory P (2006) Bradycardia and slowing of the atrioventricular conduction in mice lacking CaV3.1/1G T-type calcium channels. Circ Res. 98:1422–1430
Mangoni ME, Couette B, Bourinet E, Platzer J, Reimer D, Striessnig J, Nargeot J (2003) Functional role of L-type Cav1.3 Ca2 + channels in cardiac pacemaker activity. Proc Natl Acad Sci USA 100:5543–5548
Baig SM, Koschak A, Lieb A, Gebhart M, Dafinger C, Nürnberg G, Ali A, Ahmad I, Sinnegger-Brauns MJ, Brandt N, Engel J, Mangoni ME, Farooq M, Khan HU, Nürnberg P, Striessnig J, Bolz HJ (2011) Loss of Ca(v)1.3 (CACNA1D) function in a human channelopathy with bradycardia and congenital deafness. Nat Neurosci 14:77–84
Adachi T, Shibata S, Okamoto Y, Sato S, Fujisawa S, Ohba T, Ono K (2013) The mechanism of increased postnatal heart rate and sinoatrial node pacemaker activity in mice. J Physiol Sci 63:133–146
We thank Dr. Jianliang Wu and Mrs. Moriyama for their technical assistance. This study was supported by the KAKEN grants #21590934 (KO). P19CL6 cells were generous gift from Dr. Komuro (University of Tokyo).
This study was supported in part by Grants-in-aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (KAKEN, Nos. 15590759, 17590775) to K.O.
Ethics approval and consent to participate
All procedures performed in this study were in accordance with the guidelines on the use and care of laboratory animals as put forward by the Physiological Society of Japan and under the control of the Ethics Committee of Animal Care and Experimentation, Oita University School of Medicine, Japan.
Consent for publication
This study does not involve human subject. Consent for publication is not required for this article.
All authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zheng, M., Kang, L., Uchino, T. et al. Mitogen-activated protein kinase p38 modulates pacemaker ion channels differentiation in P19-derived pluripotent cells. J Physiol Sci 70, 39 (2020). https://doi.org/10.1186/s12576-020-00766-x
- p38MAP kinase