Next Article in Journal
Omicron BA.2 Lineage, the “Stealth” Variant: Is It Truly a Silent Epidemic? A Literature Review
Previous Article in Journal
Partial Endothelial Nitric Oxide Synthase Deficiency Exacerbates Cognitive Deficit and Amyloid Pathology in the APPswe/PS1ΔE9 Mouse Model of Alzheimer’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols

1
Department of Cardiology, Angiology and Pneumology, Heidelberg University Hospital, 69120 Heidelberg, Germany
2
German Center for Cardiovascular Research (DZHK), Partner Site Heidelberg/Mannheim, 69120 Heidelberg, Germany
3
Institute of Physiology and Pathophysiology, Heidelberg University, 69120 Heidelberg, Germany
4
EMBL Imaging Centre, European Molecular Biology Laboratory (EMBL), Meyerhofstraße 1, 69117 Heidelberg, Germany
5
Department of Molecular Neurobiology, Max-Planck-Institute for Medical Research, Jahnstraße 29, 69120 Heidelberg, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(13), 7318; https://doi.org/10.3390/ijms23137318
Submission received: 8 June 2022 / Revised: 27 June 2022 / Accepted: 28 June 2022 / Published: 30 June 2022
(This article belongs to the Topic Stem Cell Differentiation and Applications)

Abstract

:
Current protocols for the differentiation of human-induced pluripotent stem cells (hiPSC) into cardiomyocytes only generate a small amount of cardiac pacemaker cells. In previous work, we reported the generation of high amounts of cardiac pacemaker cells by co-culturing hiPSC with mouse visceral endoderm-like (END2) cells. However, potential medical applications of cardiac pacemaker cells generated according to this protocol, comprise an incalculable xenogeneic risk. We thus aimed to establish novel protocols maintaining the differentiation efficiency of the END2 cell-based protocol, yet eliminating the use of END2 cells. Three protocols were based on the activation and inhibition of the Wingless/Integrated (Wnt) signaling pathway, supplemented either with retinoic acid and the Wnt activator CHIR99021 (protocol B) or with the NODAL inhibitor SB431542 (protocol C) or with a combination of all three components (protocol D). An additional fourth protocol (protocol E) was used, which was originally developed by the manufacturer STEMCELL Technologies for the differentiation of hiPSC or hESC into atrial cardiomyocytes. All protocols (B, C, D, E) were compared to the END2 cell-based protocol A, serving as reference, in terms of their ability to differentiate hiPSC into cardiac pacemaker cells. Our analysis revealed that protocol E induced upregulation of 12 out of 15 cardiac pacemaker-specific genes. For comparison, reference protocol A upregulated 11, while protocols B, C and D upregulated 9, 10 and 8 cardiac pacemaker-specific genes, respectively. Cells differentiated according to protocol E displayed intense fluorescence signals of cardiac pacemaker-specific markers and showed excellent rate responsiveness to adrenergic and cholinergic stimulation. In conclusion, we characterized four novel and END2 cell-independent protocols for the differentiation of hiPSC into cardiac pacemaker cells, of which protocol E was the most efficient.

1. Introduction

Diseases that impair the proper function of the human sinus node (hSAN), are currently treated by cardiac pacemaker devices. However, they bear a lot of disadvantages like lead dysfunction, the need for frequent battery change, perioperative infections requiring total system explantation and the lack of response to the autonomic nervous system [1,2]. Therefore, the implementation of new therapeutic strategies, which bypass the side effects of cardiac pacemaker devices, is paramount. In this context, the development of a biological pacemaker is very promising [3] and can be achieved by approaches consisting of cell [4,5] and genetic therapeutic strategies [5]. Due to the fact that genetic strategies bear the risk of oncogenicity [6,7], cell strategies may be the most auspicious approach. The discovery of human-induced pluripotent stem cells (hiPSC) [8] provided access to cell therapeutic strategies, especially in the field of regenerative medicine [9]. Based on their pluripotency, hiPSC can be differentiated into all types of somatic cells [8]. Protocols for the differentiation of hiPSC or human embryonic stem cells (hESC) into cardiomyocytes have already been established [10,11,12,13,14,15,16].
However, although particularly efficient regarding the differentiation of hiPSC or hESC into cardiomyocytes of ventricular phenotype, the amount of cardiac pacemaker cells received remains, at approximately 10%, very low [17]. In this regard, we established a pacemaker-specific protocol that co-cultures hiPSC with mouse visceral endoderm-like cells (END2 cells) followed by growth in fetal bovine serum (FBS)-enriched culture medium to promote differentiation toward cardiac pacemaker cells [4]. Importantly, we increased the actual yield of cardiac pacemaker cells to 63.4% without the need for recurring genetic manipulations [4]. However, the use of mouse END2 cells bears a xenogeneic risk, which might raise safety issues regarding future applications in humans [18] as well as the restricted scalability of an END2 cell-based co-culture system [19]. Recently, it was shown that inhibition of NODAL by the small molecule SB431542 during deactivation of the Wingless/Integrated (Wnt) signaling pathway by IWR-1 increases the fraction of hiPSC-derived cardiac pacemaker-like cells [20]. Similar results were obtained by simultaneous application of SB431542 and the Wnt inhibitor IWP-2 [21]. Moreover, a second application of the Wnt activator CHIR99021 as well as administration of retinoic acid (RA) after Wnt inhibition by IWR-1 promoted differentiation of hiPSC toward a more atrial- and epicardial-like fate with a reduced proportion of ventricular-like cardiomyocytes and correlated with high expression levels of T-box transcription factor 18 (Tbx18) [22], known to induce the development of cardiac pacemaker tissue [23,24]. Furthermore, the manufacturer STEMCELL Technologies recently established a protocol for the differentiation of hiPSC into cardiomyocytes with predominantly atrial-like characteristics.
Due to the limitations in future medical applications associated with the use of END2 cells [18,19], we aimed to develop new protocols for the differentiation of hiPSC into cardiac pacemaker cells that substitute END2 cell-based co-culture. In these new differentiation protocols, we used culture media and components that induce activation or inhibition of Wnt signaling [15,16,17], supplemented either with the NODAL inhibitor SB431542 or with RA and a second application of the Wnt activator CHIR99021 or combining all three components to enhance hiPSC differentiation toward cardiac pacemaker cells. The following differentiation step in FBS-enriched medium, however, remained unchanged [4]. Moreover, the protocol from the manufacturer STEMCELL Technologies was validated for its ability to differentiate hiPSC into cardiac pacemaker cells. Finally, we compared these novel protocols in terms of differentiation efficiency with the END2 cell-based pacemaker-specific protocol [4].

2. Results

2.1. Novel Cell Culture Protocols for the Differentiation of hiPSC into Cardiac Pacemaker Cells

Figure 1 gives a summary of the various protocols used for the differentiation of hiPSC into cardiac pacemaker cells. These novel protocols B–E were compared to the already published reference protocol A [4] regarding differentiation efficiency of hiPSC toward the cardiac pacemaker cell lineage. Cell differentiation by reference protocol A was initiated by a co-culture of hiPSC and END2 cells [4], whereas the novel protocols B–D were based on activation and inhibition of Wnt signaling pathways [15,16,17] supplemented either with RA and a second application of the Wnt activator CHIR99021 [22] (protocol B) or with the NODAL inhibitor SB431542 [20] (protocol C) or with all three substances (protocol D). Importantly, after initiation of spontaneous beating, cells were transferred to an FBS-enriched culture medium, identical for the protocols A–D. Differentiation protocol E consisted of the newly developed STEMdiff™ Atrial Cardiomyocyte Differentiation Kit from STEMCELL Technologies. Approximately 50% of initially seeded hiPSC were spontaneously beating if differentiated by protocols B–D. These results were similar to those achieved by reference protocol A [4]. In contrast, protocol E yielded at least 80% spontaneously beating areas.

2.2. Proliferation and Survival Analysis of Cardiac Pacemaker Cells

We used ATP-based proliferation assays to evaluate whether hiPSC-derived cardiac pacemaker cells differentiated by protocols B–E continue to differentiate. Relative light units (RLU) of cardiac pacemaker cells differentiated by protocols B-E were therefore compared to RLU of hiPSC and cardiac pacemaker cells differentiated by protocol A, which were already determined in our previous work [4]. Results of the ATP-based assays are depicted in Figure 2. Compared to hiPSC, which showed an exponential proliferation (Figure 2B) [4], hiPSC-derived cardiac pacemaker cells showed constant RLU, hence demonstrating that they did not proliferate any more during and at the end of differentiation (Figure 2A). Of note, results were similar for all novel differentiation protocols (B–E) and comparable to those of reference protocol A [4]. Moreover, as the RLU remained constant during differentiation, we can also conclude that there was no relevant cell death, thus indicating that survival of hiPSC-derived cardiac pacemaker cells was not significantly impaired.

2.3. Valuation of the New Protocols’ Differentiation Efficiency toward Cardiac Pacemaker Cells by qRT-PCR

Our aim was to estimate the efficiency of the novel protocols to differentiate hiPSC into cardiac pacemaker cells. For that purpose, two distinct references were used to evaluate the gene expression profile of cardiac pacemaker cells generated by the novel protocols B-E. As the first reference, we used the transcription profile obtained from cardiac pacemaker cells generated by protocol A, based on co-culture with END2 cells [4]. Statistically significant differences regarding transcription levels between each of the novel protocols and reference protocol A are marked by a hashtag (#) (Figure 3 and Figure 4). As a second reference, we chose the transcription profile of the human right atrium (hRA) and human sinoatrial node (hSAN) [4]. Statistically significant differences between each of the differentiation protocols and hRA are indicated by an asterisk (*) (Figure 3 and Figure 4). Of note, hSAN transcript levels were not determined, but calculated by multiplying the hRA transcript levels with a hSAN/hRA ratio [4], on the basis of data published by Chandler et al. [25].

2.3.1. Pacemaker-Specific Transcription Factors

The pacemaker-specific transcription factors T-box transcription factors 3 (Tbx3) and 18 (Tbx18) play an important role in the development of cardiac pacemaker tissue through activation of nodal-specific genetic pathways leading to SAN formation and suppression of pathways involved in the generation of working myocardium [23,26]. Cardiac pacemaker cells generated from hiPSC by protocol E showed TBX18 transcript levels similar to those in hRA and cells differentiated by reference protocol A (Figure 3B). In contrast, TBX3 transcript levels were significantly higher in cardiac pacemaker cells generated by protocol E than in hRA but lower than in cells differentiated according to reference protocol A (Figure 3A). Importantly, cardiac pacemaker cells differentiated by protocols B, C and D showed significantly lower TBX3 and TBX18 gene expression levels than hRA and cells generated by reference protocol A (Figure 3A,B). We also determined the transcript levels of the transcription factor SHOX2 and its target Bmp4, both involved in the development of cardiac pacemaker cells [27]. SHOX2 was highly transcribed in cardiac pacemaker cells generated by protocol E with levels similar to those in hRA and significantly higher than those in cells differentiated by reference protocol A (Figure 3C). In contrast, cells generated by protocols B, C and D only displayed very low gene expression levels of SHOX2 (Figure 3C). BMP4, however, was markedly transcribed in cardiac pacemaker cells generated by protocols C and D reaching levels similar to those in cells differentiated by reference protocol A and exceeding those in hRA (Figure 3D). In cells generated by protocols B and E, BMP4 levels remained relatively low though (Figure 3D).

2.3.2. Transcription Factors and Markers of the Working Myocardium

T-box transcription factor 5 (Tbx5), homeobox protein Nkx2.5 and myocyte enhancer factor 2C (Mef2c) promote differentiation into working myocardium [28]. While Tbx5 and Nkx2.5 are both involved in the development of ventricular cardiomyocytes [29], Nkx2.5 was shown to reduce SAN formation in embryonic mice [30], indicating a more suppressive role for Nkx2.5 on cardiac pacemaker-specific cell differentiation.
In this regard, we observed significantly lower expression levels of TBX5 and MEF2C in cardiac pacemaker cells generated by the new protocols than in hRA (Figure 3E,F). NKX2.5 was scarcely transcribed in cardiac pacemaker cells generated by protocols B and E, with levels similar to those observed in hRA and in cells differentiated by protocol A (Figure 3G). In contrast, cardiac pacemaker cells generated by protocol C and D revealed high NKX2.5 transcript levels (Figure 3G). The fact that the myocardial marker cardiac Troponin I (cTNI) was barely detectable in cardiac pacemaker cells generated by differentiation protocols B, C, D and E, impressively underlines the high selectivity of these protocols toward the generation of cardiac pacemaker cells (Figure 3H).

2.3.3. Ion Channels and Transporters

The funny current (If) is one of the major cardiac pacemaker currents [31] generated by hyperpolarization-activated cyclic nucleotide channels (HCN). It is importantly involved in SAN pacemaker activity [32,33] and early embryonic heart development [34,35,36]. SAN-specific HCN-isoforms are first of all HCN4 and to a lower extent HCN1 [25,36]. HCN2, by contrast, is predominantly expressed in atrial cells [25]. There is, however, some evidence that it also contributes to SAN function [37,38]. hiPSC-differentiated cardiac pacemaker cells by protocol E showed very high transcription levels of HCN1 (Figure 4A) and HCN4 (Figure 4C), clearly exceeding those of hRA and hSAN. HCN4 was even higher transcribed in cardiac pacemaker cells generated by protocol E than in cells differentiated by reference protocol A (Figure 4C). Interestingly, cells generated by protocols B, C and D showed HCN4 transcript levels exceeding those of hRA or hSAN, although they still ranged below those of cardiac pacemaker cells generated by protocols A and E (Figure 4C). This was in fact different from HCN1, where cells generated by protocols B, C and D showed lower transcript levels than hRA or hSAN (Figure 4A). HCN2, by contrast, was abundantly transcribed in cardiac pacemaker cells generated by protocols D and E, displaying significantly higher gene expression levels than hRA and cells derived from protocol A (Figure 4B).
The sodium-calcium-exchanger 1 (NCX1), also crucial for cardiac pacemaking [39], displayed high transcription levels in all cardiac pacemaker cells, independent of the applied protocols and surpassed transcript levels assessed in hRA or hSAN (Figure 4D). However, the highest NCX1 transcript levels were observed in cells generated by protocol E (Figure 4D). L-type (Cav1.2 and Cav1.3) and t-type (Cav3.1) calcium channels as part of the membrane clock system contribute importantly to the pacemaker current [40,41]. Accordingly, cardiac pacemaker cells generated by the differentiation protocols A, B, C, D and E showed transcript levels of Cav1.2, Cav1.3 and Cav3.1 similar or even higher than those observed in hRA or hSAN (Figure 4E–G). Interestingly, Cav1.2 showed the highest transcript levels in cells generated by protocol C (Figure 4E), while cells generated by protocol E led to the highest levels of Cav1.3 and Cav3.1 transcripts (Figure 4F,G).
Transcription of the major depolarizing cardiac sodium current Nav1.5, responsible for the action potential upstroke especially in ventricular-like cardiomyocytes [42], remained low, with the exception of cells generated by protocols C and D with levels comparable to those in hRA or hSAN (Figure 4H). Transcription of the inward rectifier Kir2.1 and potassium channel Kv4.3 was low in all cardiac pacemaker cells regardless of the used differentiation protocol, showing transcript levels below those in hRA and comparable to those in the hSAN (Figure 4I,K). As distinct from cardiac pacemaker cells generated by protocols C and D, we observed only minor hERG transcript levels in cells differentiated by protocols A, B and E (Figure 4J).

2.3.4. Connexins

Propagation of electrical signals within cardiac tissue is ensured by connexins (Cx) [43]. Cx45 is associated with the SAN and the conduction system [43], while Cx40 and Cx43 belong to the working myocardium [44,45]. Cx45 transcript levels were highest in cardiac pacemaker cells generated by protocol E, clearly exceeding levels in cells differentiated by reference protocol A (Figure 4L). In contrast, cells generated by protocols B, C and D showed significantly lower levels (Figure 4L). Independently of the used differentiation protocol, Cx40 was transcribed on negligible levels in all cells far below hRA or hSAN levels (Figure 4M). Similarly, ventricle-specific Cx43 was scarcely transcribed in cells generated by protocols A and E with levels similar to those in the hSAN (Figure 4N), while significant higher levels emerged in cells generated by protocols B, C and D (Figure 4N).

2.3.5. Adrenergic and Cholinergic Receptors

Adrenergic and cholinergic receptors are required for beating rate modulation of cardiac pacemaker cells [46,47]. Except for cells generated by protocols A and B, cardiac pacemaker cells revealed transcription levels of the α1a-adrenoreceptor (AR) higher than in hRA or hSAN (Figure 4O). Similar results were obtained for the β2-AR (Figure 4Q). This was different from β1-AR, where the highest transcript levels assessed in cardiac pacemaker cells generated by protocols A, B and E, still ranged below hRA or hSAN levels (Figure 4P). Cholinergic receptor CHMR2 transcription was highest in cells generated by protocols B and E and clearly surmounted hRA levels (Figure 4R).

2.4. Comparison of the Differentiation Protocols Regarding the Transcription of Cardiac Pacemaker-Specific Genes

Table 1 represents an overview of the number of cardiac pacemaker-specific genes transcribed in cells generated by the differentiation protocols A-E. Genes were marked by a “+” sign if transcript levels were similar or higher than in hRA. The highest number of cardiac pacemaker-specific genes transcribed was observed in cells generated by protocol E (12/15 genes), even surpassing the number in cells differentiated by reference protocol A (11/15 genes). Cells generated by protocol D revealed the lowest number of active cardiac pacemaker-specific genes (8/15 genes).
In conclusion, cardiac pacemaker cells generated by protocol E displayed the best profile in terms of pacemaker specificity and hence were used exclusively for further characterization including immunocytochemistry and pharmacological testing.

2.5. Immunocytochemical Analysis of Protocol E Derived Cardiac Pacemaker Cells

Confocal imaging showed pronounced fluorescence signals of the cardiac pacemaker-specific markers HCN1, HCN4, NCX1, Cav1.2, Cx45, Tbx3, Tbx18 and Shox2 (Figure 5A). By contrast, we observed only weak signals of the ventricular-like markers Nav1.5, Cx43 and Nkx2.5 (Figure 5B) and moderate signals of the atrial-like markers HCN2 and Cx40 (Figure 5C).

2.6. Spontaneous Beating Rate and Pharmacological Testing

The spontaneous beating rate of protocol E-derived cardiac pacemaker cells was assessed by video recordings and analyzed using a structural similarity approach (SSIM). Cells exhibited a spontaneous beating rate of 70.3 ± 4.8 beats per minute (bpm) (n = 30) (Figure 6A). Upon stimulation with 1 µM of the adrenergic substance isoproterenol, spontaneous beating rate significantly increased to 136.0 ± 7.3 bpm (p < 0.001) (n = 30) (Figure 6B). By contrast, treatment of cells with 1 µM of the cholinergic agent carbachol led to a significant decrease in spontaneous beating rate to 36.4 ± 2.8 bpm (p < 0.001) (n = 30) (Figure 6C). A similar effect was observed when we applied the HCN-channel inhibitor ivabradine in a concentration of 3 µM [48] that significantly decreased spontaneous beating rate to 40.4 ± 2.4 bpm (p < 0.001) (n = 30) (Figure 6D), pointing to the fact that there is no significant difference between carbachol and ivabradine regarding the ability to decrease beating rate (p = 0.357).

3. Discussion

Cardiac pacemaker devices are actually the only possibility to treat diseases affecting the hSAN. Due to their disadvantages including side effects and the missing neurohumoral regulation of rate response, the development of a biological pacemaker is of utmost importance [1,2,3]. Considering the increased oncogenic risk of a genetic therapeutic strategy [6,7], a cell-based approach appears to be very promising [4,5]. Owing to their pluripotency, hiPSC [8] represent a beneficial cell source for cell-based therapy [9]. Differentiation of hiPSC toward cardiomyocytes has already been successfully realized [10,11,12,13,14,15,16]. They consisted of the induction of mesoderm formation using either Wnt activation by CHIR99021, a small molecule inhibiting glycogen synthase kinase (GSK3) [15,16] or activation of the bone morphogenetic protein (BMP) and ACTIVIN A pathways [11,12,21]. After mesoderm formation, further differentiation into cardiomyocytes was achieved by inhibition of the Wnt pathways using the Wnt inhibitors IWP-2 or IWR-1 [15,16,21] or by culturing mesodermal cells after embryoid body formation in monocultures [11,12,14]. However, following these protocols, only 10% of the cells developed a pacemaker-like phenotype [17]. Overexpression of pacemaker-specific genes like TBX3 [49] or SHOX2 [50] after transfer into mouse pluripotent stem cells or of TBX18 in cardiomyocytes of the working myocardium [6,7] resulted in a higher number of pacemaker-like cells, but, raised safety issues due to the increased risk of oncogenicity [21]. By contrast, in previous work, we were able to establish a protocol yielding a higher amount of cardiac pacemaker cells without the need for recurring genetic modifications [4].
Our protocol mainly consisted of two steps. In a first step, co-culture of hiPSC with END2 cells was performed to induce spontaneously beating cell clusters [10]. In a second step, beating cell clusters were transferred to an FBS-enriched cell culture medium to promote further differentiation toward a cardiac pacemaker-like phenotype [51]. Following our protocol, 63.4% of generated cells had a cardiac pacemaker-like phenotype [4].
However, the xenogeneic origin of END2 cells raises serious safety issues regarding future applications of cardiac pacemaker cells in a human clinical setting [18]. In addition, an END2 cell-mediated differentiation is less scalable than a differentiation by defined factors [19]. Therefore, we aimed to establish novel cell culture conditions that yield an increased differentiation efficiency comparable to our previous protocol [4] and substitute the use of END2 cells.
We wanted to replace the first step of our previous protocol, namely the co-culture of hiPSC and END2 cells by xenogeneic-free small molecules involved in the activation and inhibition of the Wnt signaling pathway [15,16,17]. The second step, however, remained unchanged and consisted of the maturation of spontaneous beating cell clusters in FBS-enriched cell culture medium as previously described [4]. It has been recently shown that inhibition of NODAL by the small molecule SB431542 promotes the differentiation of hiPSC into cardiac pacemaker cells [20,21]. The addition of the NODAL inhibitor SB431542 to the Wnt inhibitor IWR-1 led to a 1.4–2.4-fold higher amount of hiPSC-derived cardiac pacemaker-specific cells compared to solely treatment with IWR-1 [20]. Furthermore, it was shown that administration of retinoic acid (RA), as well as a second application of the Wnt activator CHIR99021 after cell treatment by the Wnt inhibitor IWR-1, induced an atrial- and epicardial-like differentiation of hiPSC [21].
In consequence, our novel protocols B-D were based on activation and inhibition of Wnt signaling, supplemented either with RA and a second application of the Wnt activator CHIR99021 (protocol B) or with the NODAL inhibitor SB431542 (protocol C) or with all three substances (protocol D). Moreover, we wanted to examine the capability of the STEMdiff™ Atrial Cardiomyocyte Differentiation Kit from STEMCELL Technologies to differentiate hiPSC into cardiac pacemaker cells (protocol E), while our previously published END2 cell-based protocol for cardiac pacemaker-specific differentiation served as a reference (protocol A) [4].
To evaluate the differentiation efficiency of the protocols B-E, we compared the transcription profile of cardiac pacemaker cells generated by protocols B-E to that of cardiac pacemaker cells generated by reference protocol A [4]. Importantly, only spontaneously beating areas of 56-day-old cells were used for analysis.
Although all protocols induced the upregulation of cardiac pacemaker-specific genes upon differentiation (Table 1), we observed significant differences regarding their differentiation efficiency. The administration of the NODAL inhibitor SB431542 (protocol C) was more efficient than a second application of the Wnt inhibitor CHIR99021 in combination with RA (protocol B) and led to an upregulation of 10 cardiac pacemaker-specific genes, whereas protocol B only upregulated 9 out of the 15 genes analyzed. Of note, protocol D displayed the worst performance with only eight genes upregulated, thus indicating that treatment with a combination of the Wnt inhibitor IWR-1 and NODAL inhibitor SB431542 followed by a combination of Wnt activator CHIR99021 and RA caused no enhanced differentiation efficiency toward cardiac pacemaker cells.
Despite being commercialized as a kit for atrial-like differentiation of hiPSC, protocol E led to significant overexpression of most of the cardiac pacemaker-specific genes analyzed, namely 12 out of 15. For comparison, reference protocol A [4] induced upregulation of 11 cardiac pacemaker-specific genes. Importantly, the pacemaker-specific transcription factor Shox2, which contributes to the development of cardiac pacemakers by upregulation of the pro-pacemaking transcription factors Tbx3 and Tbx18 [50,52], showed very high transcription levels in cardiac pacemaker cells generated by protocol E, which were clearly different from cells generated by the other protocols A, B, C and D. Similarly, transcription of Tbx18, which is critically involved in the formation of the sinus node head [53] showed highest levels in cells differentiated by protocols A and E, whereas Tbx3, important for the pacemaker-specific differentiation within cardiac tissue [53,54,55], showed highest levels in cells generated by protocol A, followed by cells differentiated by protocol E, thus clearly surpassing hRA levels. In contrast, much lower levels were observed in cells generated by protocols B-D.
Consistent with the transcription profiles, immunocytochemical analysis of cardiac pacemaker cells generated by protocol E displayed intense fluorescence signals of Tbx3 and Tbx18, whereas expression of Nkx2.5, known to suppress hSAN formation [30] by downregulating Shox2 activity [56], was nearly undetectable. In conclusion, high SHOX2 transcript levels, as well as intense Shox2 protein signals detected in cardiac pacemaker cells generated by protocol E, strongly indicate that this culture condition favors the differentiation of hiPSC toward a cardiac pacemaker cell lineage.
It could be shown that Shox2 leads to an upregulation of the cardiac pacemaker-specific ion channel HCN4 [57], which next to HCN1, is responsible for the generation of the spontaneous depolarizing pacemaker current If in sinoatrial cells [25,32,33,36]. Similarly, Tbx3 is known to activate the gene expression of HCN4 and of the cardiac pacemaker-specific calcium channel Cav3.1 [55]. Conversely, Tbx3 suppresses the expression of myocardial genes like Cx40, Cx43 and SCN5A [55]. Moreover, it could be shown that overexpression of Tbx18 in hiPSC induced upregulation of HCN4 and downregulation of Kir2.1 [24]. The non-pacemaker-specific, ventricular-like transcription factor Nkx2.5, by contrast, suppresses HCN4 as well as TBX3 [58] and upregulates the atrial marker Cx40 [59]. In line with this, protocol E led to pronounced transcription of the pacemaker-specific ion channels HCN1 and HCN4. These findings were confirmed by confocal analysis showing high protein levels of HCN1 and HCN4 in cells generated by protocol E. Similarly, we observed a strongly increased transcription for cardiac pacemaker-specific ion channels like NCX1 [39], Cav1.2, Cav1.3 and Cav3.1 [40,41] in these cells, substantiated by intense immunocytochemical signals. In contrast to cells produced by reference protocol A and nearly lacking HCN2 transcripts, protocol-E-induced cells showed much higher levels of the atrial marker HCN2 [25]. Although Cx40 [44,45] was higher transcribed in cells generated by protocol E than in cells differentiated by protocol A, its gene expression levels ranged, however, far below those of hRA or hSAN. Transcription of ventricular and non-pacemaker markers such as cTnI, SCN5A [42], Cx43 [44,45] and Kir2.1 [60] remained relatively low too. Particularly, the very low transcription as well as the weak immunocytochemical signals of the ventricular-specific markers Nav1.5 [42] and Cx43 [44,45] underlined that culturing hiPSC following protocol E excludes a ventricular-like differentiation of hiPSC. Moreover, the appearance of only moderate fluorescence signals of the atrial markers HCN2 and Cx40 confirmed these observations. In contrast, the cardiac pacemaker-specific marker Cx45 was highly detected on mRNA as well as on protein levels. These findings clearly indicate that protocol E preferentially upregulates the expression of cardiac pacemaker-specific markers, even more efficiently than reference protocol A [4]. Moreover, the fact that the cardiac pacemaker-specific relationship between transcription factors and their target genes is maintained in cardiac pacemaker cells differentiated by protocol E, emphasizes once more the high potential of protocol E to differentiate hiPSC into cardiac pacemaker cells.
To assess the maturity of hiPSC-derived cardiac pacemaker cells and to exclude any proliferative activity during differentiation, we performed ATP-based cell proliferation assays. In contrast to hiPSC, which showed exponential growth, no significant cell proliferation or cell death could be recorded for cardiac pacemaker cells, independently of the used differentiation protocol. These findings are insofar important because they underline the stability of hiPSC-derived cardiac pacemaker cells, which is a prerequisite for a functioning biological pacemaker [61].
The rate responsiveness of sinoatrial cells to adrenergic and cholinergic substances is paramount for a proper cardiac pacemaker function [46,47]. Cardiac pacemaker cells generated by protocol E not only displayed high gene expression levels of adrenergic and cholinergic receptors but also revealed appropriate modulation by the adrenergic substance isoproterenol and the cholinergic substance carbachol. Stimulation with 1 µM isoproterenol led to a significant increase in spontaneous beating rate from 70.3 ± 4.8 bpm to 136.0 ± 7.3 bpm (p < 0.001) indicating an excellent adrenergic response. Conversely, treatment with 1 µM carbachol resulted in a significant decrease in spontaneous beating rate from 70.3 ± 4.8 bpm to 36.4 ± 2.8 bpm (p < 0.001), hence demonstrating an adequate response of cardiac pacemaker cells generated by protocol E to cholinergic agents. In addition, ivabradine has been shown to specifically inhibit If in cardiac pacemaker cells at a concentration of 3 µM [48]. Equally, in cardiac pacemaker cells generated by protocol E, 3 µM ivabradine lowered the spontaneous beating rate from 65.7 ± 1.3 bpm to 40.4 ± 2.4 bpm (p < 0.001). Of note, no significant changes in beating rate could be observed in hiPSC-derived cardiomyocytes, even if high concentrations of ivabradine up to 9 µM [62] were applied.

4. Materials and Methods

4.1. Differentiation of hiPSC by Novel Cell Culture Protocols

hiPSC generated from human dermal fibroblasts (Catalog number: DF-F-ZB, BioCat GmbH, Heidelberg, Germany) [4] were used for the various differentiation assays. They were cultured in mTeSR™1 medium (STEMCELL Technologies Inc., Vancouver, BC, Canada) and grown on six-well plates (Greiner Bio-One GmbH, Frickenhausen, Germany), coated with Corning® Matrigel® hESC-Qualified Matrix (Corning, Inc., Corning, NY, USA) according to the manufacturer’s instructions. We developed four novel cell culture protocols (protocols B-E) and compared them to the END2 cell-based differentiation protocol (protocol A) [4] in terms of their ability to differentiate hiPSC into cardiac pacemaker cells. A schematic overview of the differentiation protocols A-E is depicted in Figure 1. To induce cardiac pacemaker-specific differentiation of hiPSC according to cell culture protocols B-D, hiPSC colonies of 90% confluency were mechanically dissociated and cultivated in RPMI 1640 medium (Thermo Fisher Scientific Inc., Waltham, MA, USA) containing B-27 supplement without insulin (Thermo Fisher Scientific Inc., Waltham, MA, USA) from days 0–6. On days 0 + 1, cell culture medium was supplemented with 6 µM of the Wnt activator CHIR 99021 (Sigma-Aldrich, St. Louis, MO, USA), and on days 3 + 4, the Wnt inhibitor IWR-1 was added at a concentration of 5 µM (Sigma-Aldrich, St. Louis, MO, USA). From day 7 to the start of spontaneous cell beating (days 10–12) cells were cultured in RPMI 1640 medium containing B-27 supplement (including insulin) (Thermo Fisher Scientific Inc., Waltham, MA, USA). From day 13 to day 56, spontaneously beating cells were cultured in FBS-enriched medium as previously described [4]. Cell culture protocols B-D differed from each other by following modifications: protocol B, 1 µM of retinoic acid (RA, Sigma-Aldrich) and 5 µM of CHIR 99021 were added on days 5 + 6; protocol C, 5 µM of the NODAL inhibitor SB431542 (Tocris Bioscience, Bristol, UK) were added on days 3 + 4; protocol D, 1 µM of RA and 5 µM of CHIR 99021 were added on days 5 + 6 in addition to 5 µM of SB431542 on days 3 + 4; protocol E used the STEMdiff™ Atrial Cardiomyocyte Differentiation Kit from STEMCELL Technologies (STEMCELL Technologies Inc., Vancouver, BC, Canada) to induce differentiation of hiPSC. For each novel protocol (B–E) differentiation assays were performed 16 times (n = 16 six-well plates).

4.2. Cell Proliferation Assays

To assess whether hiPSC continues to proliferate during differentiation, we performed adenosine triphosphate (ATP)-based cell proliferation assays using the BioFix® Lumi ATP assay (Macherey-Nagel, Düren, Germany) according to the manufacturer’s instructions. Therefore, we analyzed cell proliferation of differentiated hiPSC (dhiPSC) cultured according to protocols B-E (n = 4 six-well plates for each protocol) on days 14, 28, 42 and 56 after the start of differentiation. For each ATP-based assay, we used dhiPSC from one well of a six-well plate, and each experiment was repeated six times.

4.3. RNA Isolation and cDNA Synthesis

On day 56, total RNA was isolated from dhiPSC cultured according to protocols B–E (n = 4 six-well plates for each protocol) by Trizol reagent (Thermo Fisher Scientific Inc.,) following the manufacturer’s instructions. Importantly, only spontaneously beating areas of dhiPSC, subsequently termed cardiac pacemaker cells, were used for RNA analysis. cDNA synthesis was realized by reverse transcription according to the manufacturer’s protocol (Superscript II; Thermo Fisher Scientific Inc.).

4.4. Quantitative Real-Time Polymerase Chain Reaction

For quantitative real-time polymerase chain reaction (qRT-PCR), we used an ABS 7500 Realtime PCR System (Thermo Fisher Scientific Inc.) according to the manufacturer’s instructions. In detail, each well of a 96-well optical detection plate was loaded with 0.5 µL cDNA, 5 µL TaqMan Fast Universal Master Mix (Thermo Fisher Scientific Inc.), 0.55 µL 6-carboxyfluorescein (FAM)-labeled pre-designed TaqMan probes and primers (TaqMan Gene Expression Assays, Thermo Fisher Scientific Inc.) and 3.95 µL distilled RNAse-free water (Thermo Fisher Scientific Inc.). The used TaqMan probes and primers (Thermo Fisher Scientific Inc.) for the characterization of cardiac pacemaker cells are listed in Table 2. Transcript levels of analyzed genes were compared to a reference consisting of predesigned primers and probes for the housekeeping genes glyceraldehyde 3-phosphate dehydrogenase (GAPDH), hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1) and beta-actin (ACTB) (all from Thermo Fisher Scientific Inc.). Assessment of gene expression relative to the combined reference of GAPDH, HPRT1 and ACTB was calculated using a modified threshold cycle (CT) quantification method as described elsewhere [63]. All PCR reactions were carried out in triplicates, and the results were calculated as the average of the triplicates. The relative abundance of gene transcripts in cardiac pacemaker cells generated by protocols B–E was compared to the already published relative abundance of gene expression in cardiac pacemaker cells differentiated by protocol A [4]. Finally, the relative abundance of gene expressions in cardiac pacemaker cells generated by protocols A–E was compared to the relative abundance of gene transcripts in the human right atrium (hRA) or human sinoatrial node (hSAN), which has already been published [4]. hSAN transcript levels were determined [4] by multiplying hRA transcript levels with a previously published hRA and hSAN gene expression ratio [25].

4.5. Immunocytochemistry of hiPSC-Derived Cardiac Pacemaker Cells

dhiPSC differentiated by protocol E were used on day 56 of cell culture for immunocytochemical analysis. For dissociation into single cells, 56-day-old spontaneous beating areas of dhiPSC cultured on six-well plates (n = 4) in Maintenance medium (STEMCELL Technologies) were treated with STEMdiff™ Cardiomyocyte Dissociation Kit (STEMCELL Technologies Inc.) according to the manufacturer’s instructions. Dissociated cardiac pacemaker cells were seeded on Matrigel-coated (Corning® Matrigel® hESC-Qualified Matrix) 0.5 cm diameter coverslips (Menzel-Gläser, Gerhard Menzel GmbH, Braunschweig, Germany) in 24-well plates (Greiner Bio-One) (n = 4) and cultured for one day in STEMdiff™ Cardiomyocyte Support Medium (component of STEMdiff™ Cardiomyocyte Dissociation Kit). After 24 h, cell culture medium was switched to Maintenance medium. For immunocytochemistry, dissociated cardiac pacemaker cells were treated with 0.5% TritonTM X-100 (Sigma-Aldrich)/PBS solution for 15 min, followed by 0.1 M glycine (Merck)/PBS solution for 60 min before being incubated in blocking solution (PBS with 2% bovine serum albumin; Sigma-Aldrich) for 2–3 h. Afterward, cells were incubated with primary antibodies, diluted 1:200 in blocking solution, overnight at 4 °C. Subsequent primary antibodies were used: anti-HCN1 (IgG, rabbit, Alomone Labs, Ltd., Jerusalem, Israel), anti-HCN2 (IgG, rabbit, Alomone Labs, Ltd.), anti-HCN4 (IgG, rabbit, Alomone Labs, Ltd.), anti-NCX1 (IgG, rabbit, Alomone Labs, Ltd.), anti-Cav1.2 (IgG, mouse, Abcam, Cambridge, UK), anti-Nav1.5 (IgG, rabbit, Alomone Labs, Ltd.), anti-Cx40 (IgG, rabbit, Alomone Labs, Ltd.), anti-Cx43 (IgG, rabbit, Alomone Labs, Ltd.), anti-Cx45 (IgG, mouse, Abcam) and anti-Tbx3 (IgG, mouse, Abcam), anti-Tbx18 (IgG, mouse, Santa Cruz Biotechnology Inc., Dallas, TX, USA), anti-Shox2 (IgG, mouse, Abcam) and anti-Nkx2.5 (IgG, rabbit, Abcam). The following day, cardiac pacemaker cells were washed three times with PBS and incubated with secondary antibodies, diluted 1:200 in blocking solution, at 4 °C for 4 h. Secondary antibodies included MFP555 goat anti-rabbit IgG (H + L) (MoBiTec GmbH, Göttingen, Germany) and MFP488 goat anti-mouse IgG (H + L) (MoBiTec GmbH). After three washes with PBS, nuclei were counterstained with 1 μg/mL 4′,6-diamidino-2-phenylindole (DAPI, Sigma-Aldrich, St. Louis, MO, USA) for 30 min. Afterward, cells were washed three times with PBS and mounted with Citifluor glycerol/PBS-solution AF1 (Agar Scientific Ltd., Stansted, UK). They were visualized by confocal microscopy using the Olympus confocal microscope FV3000 (Olympus, Tokyo, Japan). Editing and analysis of confocal images were performed with ImageJ software (Version 1.53) (National Institutes of Health, Bethesda, MD, USA).

4.6. Assessment of Spontaneous Beating Rate of Cardiac Pacemaker Cells by Movement Frequency Analysis with Structural Similarity Approach

To assess the spontaneous beating rate of cardiac pacemaker cells differentiated by protocol E (n = 24 wells of six-well plates), cell video recording was performed using the Leica DMi1 inverted microscope equipped with the camera Leica MC120 HD (Leica Camera AG, Wetzlar, Germany). Video analysis was realized by applying ImageJ software (Version 1.53). To measure and analyze the movement (beating) frequency of spontaneous beating cardiac pacemaker cells from the movie, we used a structural similarity approach implemented as a custom Fiji (Laboratory for Optical and Computational Instrumentation, University of Wisconsin-Madison, Madison, WI, USA) Groovy script. The structural similarity measurement was inspired by Wang et al. [64] and slightly modified to better fit the brightness and contrast and signal-to-noise ratio of our data. More specifically, we took a given time frame, or the time average of the whole movie as our reference frame. Then we computed the pair-wise similarity (SSIM) value, between the reference frame and each time frame of the whole movie. The pair-wise similarity score can be computed either using the whole field of view of the movie or a selected area as region-of-interest (ROI).

4.7. Pharmacological Testing

To analyze the influence of pharmacological substances on the spontaneous beating rate, we treated cardiac pacemaker cells differentiated by protocol E with 1 µM isoproterenol hydrochloride (Sigma-Aldrich) (n = 6 wells of a six-well plate), 1 µM carbachol (Sigma-Aldrich) (n = 6 wells of a six-well plate) and 3 µM ivabradine hydrochloride (Sigma-Aldrich) (n = 6 wells of a six-well plate).

4.8. Statistical Analysis

For each protocol (B-E), hiPSC differentiation was performed on 16 six-well plates. Four six-well plates were used for RNA extraction and subsequent qRT-PCR analysis and further four six-well plates for the cell proliferation assays. For protocol E, four six-well plates of dhiPSC were used for immunocytochemistry and further four six-well plates for pharmacological testing and video analysis. For each well, five video recordings were performed.
GraphPad Prism (Version 9.2.0, GraphPad Prism Software, Inc., San Diego, CA, USA) was used for statistical analysis. Data of multiple groups were compared by one-way ANOVA, followed by a Tukey post hoc test. Moreover, a two-tailed unpaired Student’s t-test was performed for pairwise comparison. Differences were judged significant at a level p < 0.05. Arithmetic mean as well as the standard error of the mean (SEM) were used for data presentation.

5. Conclusions

The best cardiac pacemaker-specific transcription profile was attained in cells generated by protocol E providing high levels of gene expression in 12 of 15 cardiac pacemaker-specific genes, while reference protocol A only induced upregulation of 11 cardiac pacemaker-specific genes [4]. Moreover, except for HCN2 and Cx40, we observed no relevant transcriptional induction of genes representative of atrial- or ventricular-like markers in cells generated by protocol E.
Although developed for atrial-like differentiation of hiPSC, protocol E led to a significant induction of a high number of genes essentially required for the differentiation of hiPSCs to cardiac pacemaker cells (HCN4, HCN1, Cav1.2, Cav1.3, Cav3.1, TBX3, TBX18 and SHOX2) with the absolute advantage to obviate the need for END2 cells.

Author Contributions

Conceptualization, F.F.D.; methodology, F.F.D., N.D.U. and Z.H.; software, F.F.D. and Z.H.; validation, F.F.D.; formal analysis, F.F.D.; investigation, F.F.D.; resources, F.F.D. and P.A.S.; data curation, F.F.D.; writing—original draft preparation, F.F.D.; writing—review and editing, F.F.D., N.D.U., Z.H., M.K., R.R., N.F. and P.A.S.; visualization, F.F.D.; supervision, P.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fondation Coeur-Daniel Wagner, Fondation de Luxembourg (F.F.D.) and by the Deutsche Forschungsgemeinschaft (project number: 398004441) (P.A.S.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank Simone Bauer for her excellent technical work as well as the advanced light microscope facility of the EMBL Heidelberg, Germany for supplying the confocal microscope. For the publication fee we acknowledge financial support by “Deutsche Forschungsgemeinschaft” within the funding programme “Open Access Publikationskosten” as well as by Heidelberg University.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. Rosen, M.R.; Robinson, R.B.; Brink, P.R.; Cohen, I.S. The road to biological pacing. Nat. Rev. Cardiol. 2011, 8, 656–666. [Google Scholar] [CrossRef] [PubMed]
  2. Rosen, M.R. Gene therapy and biological pacing. N. Engl. J. Med. 2014, 371, 1158–1159. [Google Scholar] [CrossRef] [PubMed]
  3. Cingolani, E.; Goldhaber, J.I.; Marbán, E. Next-generation pacemakers: From small devices to biological pacemakers. Nat. Rev. Cardiol. 2018, 15, 139–150. [Google Scholar] [CrossRef] [PubMed]
  4. Schweizer, P.A.; Darche, F.F.; Ullrich, N.D.; Geschwill, P.; Greber, B.; Rivinius, R.; Seyler, C.; Müller-Decker, K.; Draguhn, A.; Utikal, J.; et al. Subtype-specific differentiation of cardiac pacemaker cell clusters from human induced pluripotent stem cells. Stem Cell Res. Ther. 2017, 8, 229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cho, H.C.; Marbán, E. Biological therapies for cardiac arrhythmias: Can genes and cells replace drugs and devices? Circ. Res. 2010, 106, 674–685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kapoor, N.; Liang, W.; Marbán, E.; Cho, H.C. Direct conversion of quiescent cardiomyocytes to pacemaker cells by expression of Tbx18. Nat. Biotechnol. 2013, 31, 54–62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Hu, Y.F.; Dawkins, J.F.; Cho, H.C.; Marbán, E.; Cingolani, E. Biological pacemaker created by minimally invasive somatic reprogramming in pigs with complete heart block. Sci. Transl. Med. 2014, 6, 245ra94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Takahashi, K.; Yamanaka, S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell 2006, 126, 663–676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Matsa, E.; Burridge, P.W.; Wu, J.C. Human stem cells for modeling heart disease and for drug discovery. Sci. Transl. Med. 2014, 6, 239ps6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Mummery, C.; Ward-van Oostwaard, D.; Doevendans, P.; Spijker, R.; van den Brink, S.; Hassink, R.; van der Heyden, M.; Opthof, T.; Pera, M.; de la Riviere, A.B.; et al. Differentiation of human embryonic stem cells to cardiomyocytes: Role of coculture with visceral endoderm-like cells. Circulation 2003, 107, 2733–2740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Laflamme, M.A.; Chen, K.Y.; Naumova, A.V.; Muskheli, V.; Fugate, J.A.; Dupras, S.K.; Reinecke, H.; Xu, C.; Hassanipour, M.; Police, S.; et al. Cardiomyocytes derived from human embryonic stem cells in pro-survival factors enhance function of infarcted rat hearts. Nat. Biotechnol. 2007, 25, 1015–1024. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, L.; Soonpaa, M.H.; Adler, E.D.; Roepke, T.K.; Kattman, S.J.; Kennedy, M.; Henckaerts, E.; Bonham, K.; Abbott, G.W.; Linden, R.M.; et al. Human cardiovascular progenitor cells develop from a KDR+ embryonic-stemcell-derived population. Nature 2008, 453, 524–528. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, J.; Wilson, G.F.; Soerens, A.G.; Koonce, C.H.; Yu, J.; Palecek, S.P.; Thomson, J.A.; Kamp, T.J. Functional cardiomyocytes derived from human induced pluripotent stem cells. Circ. Res. 2009, 104, e30–e41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zwi, L.; Caspi, O.; Arbel, G.; Huber, I.; Gepstein, A.; Park, I.H.; Gepstein, L. Cardiomyocyte differentiation of human induced pluripotent stem cells. Circulation 2009, 120, 1513–1523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lian, X.; Hsiao, C.; Wilson, G.; Zhu, K.; Hazeltine, L.B.; Azarin, S.M.; Raval, K.K.; Zhang, J.; Kamp, T.J.; Palecek, S.P. Robust cardiomyocyte differentiation from human pluripotent stem cells via temporal modulation of canonical Wnt signaling. Proc. Natl. Acad. Sci. USA 2012, 109, E1848–E1857. [Google Scholar] [CrossRef] [Green Version]
  16. Lian, X.; Zhang, J.; Azarin, S.M.; Zhu, K.; Hazeltine, L.B.; Bao, X.; Hsiao, C.; Kamp, T.J.; Palecek, S.P. Directed cardiomyocyte differentiation from human pluripotent stem cells by modulating Wnt/beta-catenin signaling under fully defined conditions. Nat. Protoc. 2013, 8, 162–175. [Google Scholar] [CrossRef] [Green Version]
  17. Yechikov, S.; Copaciu, R.; Gluck, J.M.; Deng, W.; Chiamvimonvat, N.; Chan, J.W.; Lieu, D.K. Same-single-cell analysis of pacemaker-specific markers in human induced pluripotent stem cell-derived cardiomyocyte subtypes classified by electrophysiology. Stem Cells 2016, 34, 2670–2680. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, X.Q.; Graichen, R.; Soo, S.Y.; Balakrishnan, T.; Rahmat, S.N.; Sieh, S.; Tham, S.C.; Freund, C.; Moore, J.; Mummery, C.; et al. Chemically defined medium supporting cardiomyocyte differentiation of human embryonic stem cells. Differentiation 2008, 76, 958–970. [Google Scholar] [CrossRef]
  19. Graichen, R.; Xu, X.; Braam, S.R.; Balakrishnan, T.; Norfiza, S.; Sieh, S.; Soo, S.Y.; Tham, S.C.; Mummery, C.; Colman, A.; et al. Enhanced cardiomyogenesis of human embryonic stem cells by a small molecular inhibitor of p38 MAPK. Differentiation 2008, 76, 357–370. [Google Scholar] [CrossRef]
  20. Yechikov, S.; Kao, H.K.J.; Chang, C.W.; Pretto, D.; Zhang, X.D.; Sun, Y.H.; Smithers, R.; Sirish, P.; Nolta, J.A.; Chan, J.W.; et al. NODAL inhibition promotes differentiation of pacemaker-like cardiomyocytes from human induced pluripotent stem cells. Stem Cell Res. 2020, 49, 102043. [Google Scholar] [CrossRef]
  21. Protze, S.I.; Liu, J.; Nussinovitch, U.; Ohana, L.; Backx, P.H.; Gepstein, L.; Keller, G.M. Sinoatrial node cardiomyocytes derived from human pluripotent cells function as a biological pacemaker. Nat. Biotechnol. 2017, 35, 56–68. [Google Scholar] [CrossRef] [PubMed]
  22. Zhao, J.; Cao, H.; Tian, L.; Huo, W.; Zhai, K.; Wang, P.; Ji, G.; Ma, Y. Efficient Differentiation of TBX18+/WT1+ Epicardial-Like Cells from Human Pluripotent Stem Cells Using Small Molecular Compounds. Stem Cells Dev. 2017, 26, 528–540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Christoffels, V.M.; Mommersteeg, M.T.; Trowe, M.O.; Prall, O.W.; de Gier-de Vries, C.; Soufan, A.T.; Bussen, M.; Schuster-Gossler, K.; Harvey, R.P.; Moorman, A.F.; et al. Formation of the venous pole of the heart from an Nkx2-5-negative precursor population requires Tbx18. Circ. Res. 2006, 98, 1555–1563. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, W.; Zhao, H.; Quan, D.; Tang, Y.; Wang, X.; Huang, C. Tbx18 promoted the conversion of human-induced pluripotent stem cell-derived cardiomyocytes into sinoatrial node-like pacemaker cells. Cell Biol. Int. 2022, 46, 403–414. [Google Scholar] [CrossRef]
  25. Chandler, N.J.; Greener, I.D.; Tellez, J.O.; Inada, S.; Musa, H.; Molenaar, P.; Difrancesco, D.; Baruscotti, M.; Longhi, R.; Anderson, R.H.; et al. Molecular architecture of the human sinus node: Insights into the function of the cardiac pacemaker. Circulation 2009, 119, 1562–1575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sizarov, A.; Devalla, H.D.; Anderson, R.H.; Passier, R.; Christoffels, V.M.; Moorman, A.F. Molecular analysis of patterning of conduction tissues in the developing human heart. Circ. Arrhythm. Electrophysiol. 2011, 4, 532–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Puskaric, S.; Schmitteckert, S.; Mori, A.D.; Glaser, A.; Schneider, K.U.; Bruneau, B.G.; Blaschke, R.J.; Steinbeisser, H.; Rappold, G. Shox2 mediates Tbx5 activity by regulating Bmp4 in the pacemaker region of the developing heart. Hum. Mol. Genet. 2010, 19, 4625–4633. [Google Scholar] [CrossRef] [Green Version]
  28. Olson, E.N. Gene regulatory networks in the evolution and development of the heart. Science 2006, 313, 1922–1927. [Google Scholar] [CrossRef] [Green Version]
  29. Moskowitz, I.P.; Kim, J.B.; Moore, M.L.; Wolf, C.M.; Peterson, M.A.; Shendure, J.; Nobrega, M.A.; Yokota, Y.; Berul, C.; Izumo, S.; et al. A molecular pathway including Id2, Tbx5, and Nkx2-5 required for cardiac conduction system development. Cell 2007, 129, 1365–1376. [Google Scholar] [CrossRef] [Green Version]
  30. Espinoza-Lewis, R.A.; Liu, H.; Sun, C.; Chen, C.; Jiao, K.; Chen, Y. Ectopic expression of Nkx2.5 suppresses the formation of the sinoatrial node in mice. Dev. Biol. 2011, 356, 359–369. [Google Scholar] [CrossRef] [Green Version]
  31. DiFrancesco, D.; Tortora, P. Direct activation of cardiac pacemaker channels by intracellular cyclic AMP. Nature 1991, 351, 145–147. [Google Scholar] [CrossRef]
  32. DiFrancesco, D. The role of the funny current in pacemaker activity. Circ. Res. 2010, 106, 434–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Baruscotti, M.; Bucchi, A.; Viscomi, C.; Mandelli, G.; Consalez, G.; Gnecchi-Rusconi, T.; Montano, N.; Casali, K.R.; Micheloni, S.; Barbuti, A.; et al. Deep bradycardia and heart block caused by inducible cardiac-specific knockout of the pacemaker channel gene Hcn4. Proc. Natl. Acad. Sci. USA 2011, 108, 1705–1710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Später, D.; Abramczuk, M.K.; Buac, K.; Zangi, L.; Stachel, M.W.; Clarke, J.; Sahara, M.; Ludwig, A.; Chien, K.R. A HCN4+ cardiomyogenic progenitor derived from the first heart field and human pluripotent stem cells. Nat. Cell. Biol. 2013, 15, 1098–1106. [Google Scholar] [CrossRef] [PubMed]
  35. Schweizer, P.A.; Schröter, J.; Greiner, S.; Haas, J.; Yampolsky, P.; Mereles, D.; Buss, S.J.; Seyler, C.; Bruehl, C.; Draguhn, A.; et al. The symptom complex of familial sinus node dysfunction and myocardial noncompaction is associated with mutations in the HCN4 channel. J. Am. Coll. Cardiol. 2014, 64, 757–767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Schweizer, P.A.; Yampolsky, P.; Malik, R.; Thomas, D.; Zehelein, J.; Katus, H.A.; Koenen, M. Transcription profiling of HCN-channel isotypes throughout mouse cardiac development. Basic Res. Cardiol. 2009, 104, 621–629. [Google Scholar] [CrossRef] [Green Version]
  37. Ludwig, A.; Budde, T.; Stieber, J.; Moosmang, S.; Wahl, C.; Holthoff, K.; Langebartels, A.; Wotjak, C.; Munsch, T.; Zong, X.; et al. Absence epilepsy and sinus dysrhythmia in mice lacking the pacemaker channel HCN2. EMBO J. 2003, 22, 216–224. [Google Scholar] [CrossRef] [Green Version]
  38. Li, N.; Csepe, T.A.; Hansen, B.J.; Dobrzynski, H.; Higgins, R.S.; Kilic, A.; Mohler, P.J.; Janssen, P.M.; Rosen, M.R.; Biesiadecki, B.J.; et al. Molecular Mapping of Sinoatrial Node HCN Channel Expression in the Human Heart. Circ. Arrhythm. Electrophysiol. 2015, 8, 1219–1227. [Google Scholar] [CrossRef] [Green Version]
  39. Herrmann, S.; Lipp, P.; Wiesen, K.; Stieber, J.; Nguyen, H.; Kaiser, E.; Ludwig, A. The cardiac sodium-calcium exchanger NCX1 is a key player in the initiation and maintenance of a stable heart rhythm. Cardiovasc. Res. 2013, 99, 780–788. [Google Scholar] [CrossRef]
  40. Maltsev, V.A.; Lakatta, E.G. Synergism of coupled subsarcolemmal Ca2+ clocks and sarcolemmal voltage clocks confers robust and flexible pacemaker function in a novel pacemaker cell model. Am. J. Physiol. Heart Circ. Physiol. 2009, 296, H594–H615. [Google Scholar] [CrossRef] [Green Version]
  41. Lakatta, E.G.; Maltsev, V.A.; Vinogradova, T.M. A coupled SYSTEM of intracellular Ca2+ clocks and surface membrane voltage clocks controls the timekeeping mechanism of the heart’s pacemaker. Circ. Res. 2010, 106, 659–673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Rook, M.B.; Evers, M.M.; Vos, M.A.; Bierhuizen, M.F. Biology of cardiac sodium channel Nav1.5 expression. Cardiovasc. Res. 2012, 93, 12–23. [Google Scholar] [CrossRef] [Green Version]
  43. Jansen, J.A.; van Veen, T.A.; de Bakker, J.M.; van Rijen, H.V. Cardiac connexins and impulse propagation. J. Mol. Cell. Cardiol. 2010, 48, 76–82. [Google Scholar] [CrossRef] [PubMed]
  44. Gaborit, N.; Le Bouter, S.; Szuts, V.; Varro, A.; Escande, D.; Nattel, S.; Demolombe, S. Regional and tissue specific transcript signatures of ion channel genes in the non-diseased human heart. J. Physiol. 2007, 582, 675–693. [Google Scholar] [CrossRef] [PubMed]
  45. Fedorov, V.V.; Schuessler, R.B.; Hemphill, M.L.; Ambrosi, C.M.; Chang, R.; Voloshina, A.S.; Brown, K.; Hucker, W.J.; Efimov, I.R. Structural and functional evidence for discrete exit pathways that connect the canine sinoatrial node and atria. Circ. Res. 2009, 104, 915–923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhang, H.; Butters, T.; Adeniran, I.; Higham, J.; Holden, A.V.; Boyett, M.R.; Hancox, J.C. Modeling the chronotropic effect of isoprenaline on rabbit sinoatrial node. Front. Physiol. 2012, 3, 241. [Google Scholar] [CrossRef] [Green Version]
  47. Behar, J.; Ganesan, A.; Zhang, J.; Yaniv, Y. The Autonomic Nervous System Regulates the Heart Rate through cAMP-PKA Dependent and Independent Coupled-Clock Pacemaker Cell Mechanisms. Front. Physiol. 2016, 7, 419. [Google Scholar] [CrossRef] [Green Version]
  48. Yaniv, Y.; Maltsev, V.A.; Ziman, B.D.; Spurgeon, H.A.; Lakatta, E.G. The “funny” current (I(f)) inhibition by ivabradine at membrane potentials encompassing spontaneous depolarization in pacemaker cells. Molecules 2012, 17, 8241–8254. [Google Scholar] [CrossRef]
  49. Jung, J.J.; Husse, B.; Rimmbach, C.; Krebs, S.; Stieber, J.; Steinhoff, G.; Dendorfer, A.; Franz, W.M.; David, R. Programming and isolation of highly pure physiologically and pharmacologically functional sinus-nodal bodies from pluripotent stem cells. Stem Cell Rep. 2014, 2, 592–605. [Google Scholar] [CrossRef] [Green Version]
  50. Ionta, V.; Liang, W.; Kim, E.H.; Rafie, R.; Giacomello, A.; Marbán, E.; Cho, H.C. SHOX2 overexpression favors differentiation of embryonic stem cells into cardiac pacemaker cells, improving biological pacing ability. Stem Cell Rep. 2015, 4, 129–142. [Google Scholar] [CrossRef] [Green Version]
  51. Scavone, A.; Capilupo, D.; Mazzocchi, N.; Crespi, A.; Zoia, S.; Campostrini, G.; Bucchi, A.; Milanesi, R.; Baruscotti, M.; Benedetti, S.; et al. Embryonic stem cell-derived CD166+ precursors develop into fully functional sinoatrial-like cells. Circ. Res. 2013, 113, 389–398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Espinoza-Lewis, R.A.; Yu, L.; He, F.; Liu, H.; Tang, R.; Shi, J.; Sun, X.; Martin, J.F.; Wang, D.; Yang, J.; et al. Shox2 is essential for the differentiation of cardiac pacemaker cells by repressing Nkx2-5. Dev. Biol. 2009, 327, 376–385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wiese, C.; Grieskamp, T.; Airik, R.; Mommersteeg, M.T.; Gardiwal, A.; de Gier-de Vries, C.; Schuster-Gossler, K.; Moorman, A.F.; Kispert, A.; Christoffels, V.M. Formation of the sinus node head and differentiation of sinus node myocardium are independently regulated by Tbx18 and Tbx3. Circ. Res. 2009, 104, 388–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hoogaars, W.M.; Tessari, A.; Moorman, A.F.; de Boer, P.A.; Hagoort, J.; Soufan, A.T.; Campione, M.; Christoffels, V.M. The transcriptional repressor Tbx3 delineates the developing central conduction system of the heart. Cardiovasc. Res. 2004, 62, 489–499. [Google Scholar] [CrossRef] [PubMed]
  55. Hoogaars, W.M.; Engel, A.; Brons, J.F.; Verkerk, A.O.; de Lange, F.J.; Wong, L.Y.; Bakker, M.L.; Clout, D.E.; Wakker, V.; Barnett, P.; et al. Tbx3 controls the sinoatrial node gene program and imposes pacemaker function on the atria. Genes Dev. 2007, 21, 1098–1112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Ye, W.; Wang, J.; Song, Y.; Yu, D.; Sun, C.; Liu, C.; Chen, F.; Zhang, Y.; Wang, F.; Harvey, R.P.; et al. A common Shox2-Nkx2-5 antagonistic mechanism primes the pacemaker cell fate in the pulmonary vein myocardium and sinoatrial node. Development 2015, 142, 2521–2532. [Google Scholar]
  57. Blaschke, R.J.; Hahurij, N.D.; Kuijper, S.; Just, S.; Wisse, L.J.; Deissler, K.; Maxelon, T.; Anastassiadis, K.; Spitzer, J.; Hardt, S.E.; et al. Targeted mutation reveals essential functions of the homeodomain transcription factor Shox2 in sinoatrial and pacemaking development. Circulation 2007, 115, 1830–1838. [Google Scholar] [CrossRef] [Green Version]
  58. Mommersteeg, M.T.; Hoogaars, W.M.; Prall, O.W.; de Gier-de Vries, C.; Wiese, C.; Clout, D.E.; Papaioannou, V.E.; Brown, N.A.; Harvey, R.P.; Moorman, A.F.; et al. Molecular pathway for the localized formation of the sinoatrial node. Circ. Res. 2007, 100, 3541–3562. [Google Scholar] [CrossRef] [Green Version]
  59. Stennard, F.A.; Harvey, R.P. T-box transcription factors and their roles in regulatory hierarchies in the developing heart. Development 2005, 132, 4897–4910. [Google Scholar] [CrossRef] [Green Version]
  60. Miake, J.; Marbán, E.; Nuss, H.B. Biological pacemaker created by gene transfer. Nature 2002, 419, 132–133. [Google Scholar] [CrossRef]
  61. Rosen, M.R.; Brink, P.R.; Cohen, I.S.; Robinson, R.B. Cardiac pacing: From biological to electronic … to biological? Circ. Arrhythm. Electrophysiol. 2008, 1, 54–61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Kim, J.J.; Yang, L.; Lin, B.; Zhu, X.; Sun, B.; Kaplan, A.D.; Bett, G.C.; Rasmusson, R.L.; London, B.; Salama, G. Mechanism of automaticity in cardiomyocytes derived from human induced pluripotent stem cells. J. Mol. Cell. Cardiol. 2015, 81, 81–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Harrell, M.D.; Harbi, S.; Hoffman, J.F.; Zavadil, J.; Coetzee, W.A. Large-scale analysis of ion channel gene expression in the mouse heart during perinatal development. Physiol. Genom. 2007, 28, 273–283. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE Trans Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic overview of the novel differentiation protocols. (A): Reference protocol A based on a co-culture of hiPSC and END2 cells. (BD): Novel protocols (BD) based on activation (6 µM CHIR 99021 on days 0 + 1) and inhibition (5 µM IWR-1 on days 3 + 4) of Wingless/Integrated (Wnt) signaling with following modifications for each protocol: (B): Additional administration of 1 µM RA and 5 µM CHIR 99021 on days 5 + 6. (C): Additional administration of 5 µM of the NODAL inhibitor SB431542 on days 3 + 4. (D): Additional administration of 5 µM of the NODAL inhibitor SB431542 on days 3 + 4 as well as of 1 µM RA and 5 µM CHIR 99021 on days 5 + 6. (E): Protocol of the manufacturer STEMCELL Technologies consisting of three differentiation compounds (AC), followed by culture in Maintenance medium. Abbreviations: hiPSC = human induced pluripotent stem cells, END2 cells = mouse visceral endoderm-like cells, FBS = fetal bovine serum, d = day, RPMI = RPMI-1640, B-27 = B-27 supplement, w/o = without, CHIR = CHIR 99021, IWR = IWR-1, RA = retinoic acid, SB = SB431542, Diff. A, B and C = differentiation compounds A, B and C.
Figure 1. Schematic overview of the novel differentiation protocols. (A): Reference protocol A based on a co-culture of hiPSC and END2 cells. (BD): Novel protocols (BD) based on activation (6 µM CHIR 99021 on days 0 + 1) and inhibition (5 µM IWR-1 on days 3 + 4) of Wingless/Integrated (Wnt) signaling with following modifications for each protocol: (B): Additional administration of 1 µM RA and 5 µM CHIR 99021 on days 5 + 6. (C): Additional administration of 5 µM of the NODAL inhibitor SB431542 on days 3 + 4. (D): Additional administration of 5 µM of the NODAL inhibitor SB431542 on days 3 + 4 as well as of 1 µM RA and 5 µM CHIR 99021 on days 5 + 6. (E): Protocol of the manufacturer STEMCELL Technologies consisting of three differentiation compounds (AC), followed by culture in Maintenance medium. Abbreviations: hiPSC = human induced pluripotent stem cells, END2 cells = mouse visceral endoderm-like cells, FBS = fetal bovine serum, d = day, RPMI = RPMI-1640, B-27 = B-27 supplement, w/o = without, CHIR = CHIR 99021, IWR = IWR-1, RA = retinoic acid, SB = SB431542, Diff. A, B and C = differentiation compounds A, B and C.
Ijms 23 07318 g001
Figure 2. ATP-based cell proliferation assays. (A): hiPSC-derived cardiac pacemaker cells differentiated by protocols A (violet), B (red), C (dark blue), D (yellow) and E (blue). (B): hiPSC. Data are provided as means ± SEM. * p < 0.05, ** p < 0.01. Abbreviations: hiPSC = human induced pluripotent stem cells, RLU = relative light units, ATP = adenosine triphosphate.
Figure 2. ATP-based cell proliferation assays. (A): hiPSC-derived cardiac pacemaker cells differentiated by protocols A (violet), B (red), C (dark blue), D (yellow) and E (blue). (B): hiPSC. Data are provided as means ± SEM. * p < 0.05, ** p < 0.01. Abbreviations: hiPSC = human induced pluripotent stem cells, RLU = relative light units, ATP = adenosine triphosphate.
Ijms 23 07318 g002
Figure 3. Gene expression of pacemaker-specific and myocardium-related transcription factors in cardiac pacemaker cells. Abundance of mRNA transcripts (in%) relative to gene expression of the housekeeping genes GAPDH/HPRT1/ACTB. (AD): Pacemaker-specific transcription factors. (EG): Myocardium-related transcription factors. (H): Myocardial marker troponin I (cTnI). Data are provided as means ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, Diff. prot. A–E versus hRA. # p < 0.05, ## p < 0.01, ### p < 0.001, Diff. prot. B–E versus Diff. prot. A. Abbreviations: hRA: human right atrium, hSAN: human sinoatrial node, Diff. prot. = differentiation protocol.
Figure 3. Gene expression of pacemaker-specific and myocardium-related transcription factors in cardiac pacemaker cells. Abundance of mRNA transcripts (in%) relative to gene expression of the housekeeping genes GAPDH/HPRT1/ACTB. (AD): Pacemaker-specific transcription factors. (EG): Myocardium-related transcription factors. (H): Myocardial marker troponin I (cTnI). Data are provided as means ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, Diff. prot. A–E versus hRA. # p < 0.05, ## p < 0.01, ### p < 0.001, Diff. prot. B–E versus Diff. prot. A. Abbreviations: hRA: human right atrium, hSAN: human sinoatrial node, Diff. prot. = differentiation protocol.
Ijms 23 07318 g003
Figure 4. Gene expression of ion channels, transporters, connexins, adrenergic and cholinergic receptors in cardiac pacemaker cells. Abundance of mRNA transcripts (in%) relative to gene expression of the housekeeping genes GAPDH/HPRT1/ACTB. (AC): HCN channels. (D): Na-Ca exchanger NCX1. (EG): calcium channels. (H): sodium channel SCN5A. (IK): potassium channels. (L–N): connexins. (OR): adrenergic and cholinergic receptors. Data are provided as means ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, Diff. prot. A–E versus hRA. # p < 0.05, ## p < 0.01, ### p < 0.001, Diff. prot. B–E versus Diff. prot. A. Abbreviations: hRA: human right atrium, hSAN: human sinoatrial node, Diff. prot. = differentiation protocol.
Figure 4. Gene expression of ion channels, transporters, connexins, adrenergic and cholinergic receptors in cardiac pacemaker cells. Abundance of mRNA transcripts (in%) relative to gene expression of the housekeeping genes GAPDH/HPRT1/ACTB. (AC): HCN channels. (D): Na-Ca exchanger NCX1. (EG): calcium channels. (H): sodium channel SCN5A. (IK): potassium channels. (L–N): connexins. (OR): adrenergic and cholinergic receptors. Data are provided as means ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, Diff. prot. A–E versus hRA. # p < 0.05, ## p < 0.01, ### p < 0.001, Diff. prot. B–E versus Diff. prot. A. Abbreviations: hRA: human right atrium, hSAN: human sinoatrial node, Diff. prot. = differentiation protocol.
Ijms 23 07318 g004
Figure 5. Immunocytochemical analysis of cardiac pacemaker cells differentiated by protocol E. Signals were visualized by confocal microscopy: (A): Pacemaker-specific markers. (B): Ventricular- specific markers. (C): Atrial-specific markers. Nuclei were counterstained with 4′,6-diamidino-2-phenylindole (DAPI). Scale bars = 10 μm.
Figure 5. Immunocytochemical analysis of cardiac pacemaker cells differentiated by protocol E. Signals were visualized by confocal microscopy: (A): Pacemaker-specific markers. (B): Ventricular- specific markers. (C): Atrial-specific markers. Nuclei were counterstained with 4′,6-diamidino-2-phenylindole (DAPI). Scale bars = 10 μm.
Ijms 23 07318 g005
Figure 6. Spontaneous beating rate of cardiac pacemaker cells differentiated by protocol E. Cell video recording was performed using an inverted microscope equipped with a camera. Video analysis was realized by applying ImageJ software (Version 1.53) and structural similarity approach. (A): Basal spontaneous beating rate. (B): Spontaneous beating rate upon stimulation with 1 µM isoproterenol. (C): Spontaneous beating rate upon treatment with 1 µM carbachol. (D): Spontaneous beating rate upon treatment with 3 µM ivabradine. Abbreviation: SSIM: pair-wise similarity value.
Figure 6. Spontaneous beating rate of cardiac pacemaker cells differentiated by protocol E. Cell video recording was performed using an inverted microscope equipped with a camera. Video analysis was realized by applying ImageJ software (Version 1.53) and structural similarity approach. (A): Basal spontaneous beating rate. (B): Spontaneous beating rate upon stimulation with 1 µM isoproterenol. (C): Spontaneous beating rate upon treatment with 1 µM carbachol. (D): Spontaneous beating rate upon treatment with 3 µM ivabradine. Abbreviation: SSIM: pair-wise similarity value.
Ijms 23 07318 g006
Table 1. Assessment of pacemaker-specific gene expression in cardiac pacemaker cells.
Table 1. Assessment of pacemaker-specific gene expression in cardiac pacemaker cells.
Diff. Prot. ADiff. Prot. BDiff. Prot. CDiff. Prot. DDiff. Prot. E
TBX3++++
TBX18++
SHOX2+
BMP4+++
HCN1++
HCN4+++++
NCX1+++++
Cav1.2++++
Cav1.3+++++
Cav3.1+++++
Cx45
α1a-AR++++
β1-AR
β2-AR+++++
CHRM2+++++
Total11/159/1510/158/1512/15
Table 2. Overview of the used TaqMan probes and primers.
Table 2. Overview of the used TaqMan probes and primers.
Analyzed GenePrimer Accession Number (TaqMan Probe and Primer)
AbbreviationWritten-Out Form
GAPDHglyceraldehyde-3-phosphate dehydrogenaseHs99999905_m1
HPRT1hypoxanthine phosphoribosyltransferase 1Hs01003270_m1
ACTBactin betaHs010606665_m1
TBX3T-box transcription factor 3Hs00195612_m1
TBX18T-box transcription factor 18Hs01385458_m1
SHOX2short stature homeobox 2Hs00243203_m1
BMP4bone morphogenetic protein 4Hs00370078_m1
TBX5T-box transcription factor 5Hs01052563_m1
MEF2Cmyocyte enhancer factor 2CHs01554602_g1
NKX2.5NK2 homeobox 5Hs00231763_m1
cTnItroponin I3, cardiac typeHs01036382_g1
HCN1hyperpolarization activated cyclic nucleotide gated potassium channel 1Hs00395037_m1
HCN2hyperpolarization activated cyclic nucleotide gated potassium channel 2Hs00606903_m1
HCN4hyperpolarization activated cyclic nucleotide gated potassium channel 4Hs00175760_m1
NCX1sodium/calcium exchanger proteinHs01062258_m1
Cav1.2calcium voltage-gated channel subunit alpha1 CHs00167681_m1
Cav1.3calcium voltage-gated channel subunit alpha1 DHs00167753_m1
Cav3.1calcium voltage-gated channel subunit alpha1 GHs00367969_m1
SCN5Asodium voltage-gated channel alpha subunit 5Hs00165693_m1
Kir2.1potassium inwardly rectifying channel subfamily J member 2Hs00542478_m1
hERGpotassium voltage-gated channel subfamily H member 2Hs00265315_m1
Kv4.3potassium voltage-gated channel subfamily D member 3Hs00542597_m1
Cx45gap junction protein gamma 1Hs01087407_s1
Cx40gap junction protein alpha 5Hs99999170_s1
Cx43gap junction protein alpha 1Hs00748445_s1
α1a-ARadrenoceptor alpha 1AHs00169124_m1
β1-ARadrenoceptor beta 1Hs02330048_s1
β2-ARadrenoceptor beta 2Hs00240532_s1
CHRM2cholinergic receptor muscarinic 2Hs00265208_s1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Darche, F.F.; Ullrich, N.D.; Huang, Z.; Koenen, M.; Rivinius, R.; Frey, N.; Schweizer, P.A. Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols. Int. J. Mol. Sci. 2022, 23, 7318. https://doi.org/10.3390/ijms23137318

AMA Style

Darche FF, Ullrich ND, Huang Z, Koenen M, Rivinius R, Frey N, Schweizer PA. Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols. International Journal of Molecular Sciences. 2022; 23(13):7318. https://doi.org/10.3390/ijms23137318

Chicago/Turabian Style

Darche, Fabrice F., Nina D. Ullrich, Ziqiang Huang, Michael Koenen, Rasmus Rivinius, Norbert Frey, and Patrick A. Schweizer. 2022. "Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols" International Journal of Molecular Sciences 23, no. 13: 7318. https://doi.org/10.3390/ijms23137318

APA Style

Darche, F. F., Ullrich, N. D., Huang, Z., Koenen, M., Rivinius, R., Frey, N., & Schweizer, P. A. (2022). Improved Generation of Human Induced Pluripotent Stem Cell-Derived Cardiac Pacemaker Cells Using Novel Differentiation Protocols. International Journal of Molecular Sciences, 23(13), 7318. https://doi.org/10.3390/ijms23137318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop