Figure 1. Sinoatrial node subdomains are divided based on anatomical and genetic criteria. (
A) The sinoatrial node is displayed on a schematic diagram of the heart (red, central SAN region; blue, transitional cells; yellow, sinoatrial conduction pathways; white dashed line, SAN head; black dashed line, SAN tail). Adapted with permission from
[24][28], Copyright (2005, 2017). (
B) Gene expression profile in the SAN region (expressed genes are indicated in red). Adapted with permission from
[7], Copyright (2009). (
C) Samples of action potentials (APs) are recorded from the central sinoatrial node to the periphery of the node. Adapted with permission from
[20], Copyright (2016).
3. Transcriptional Regulation and Specific Molecular Markers of the SAN
SAN formation and function are tightly regulated by a TF network that displays a dynamic and unique expression pattern in the SAN and surrounding AM
[38]. These TFs are key players in the development and differentiation of PCs and help maintain pacemaker identity and function. A previous study on SAN formation and development used a gene-deficient mouse model; their findings revealed important aspects of SAN development and differentiation
[4][5][6][7][8] (
Figure 2). However, how key developmental regulators are regulated in individual cells at specific locations during the SAN development remains unclear.
Figure 2. A gene regulatory network that controls SAN development and function. The black arrows indicate gene activation. Black lines with blunt ends represent the inhibition of gene expression by their corresponding transcription factors. The dotted lines indicate indirect or speculation effects. The right box represents the SAN program, and the blue box indicates the SAN-related gene discovered by single-cell sequencing. The left upper box represents the atrial program; the left lower box indicates the genes regulated by SHOX2. The dotted box indicates the potential transcription factors associated with the SAN.
One of the characteristics of the SAN is the high expression of the hyperpolarization-activated and cyclic-nucleotide-gated ion channel HCN4
[39], which underlies the funny current (If), an essential factor for the maintenance of sinus rhythm
[40]. In humans,
HCN4 is expressed in PCs and RA, and so cannot be used for identifying human PCs as the sole marker
[41]. CD166/Alcam, a cell surface molecule, can help identify CD166
+ pacemaker precursors from differentiating mouse embryonic stem cells. However, CD166/Alcam is not expressed in PCs; thus, its application is limited
[42][43]. Therefore, specific cell surface molecules of PCs are still being sought.
3.1. Transcription Factors
Recently, the critical regulatory role of
ISL1 in SAN development has received increasing attention
[44][45]. The specific role of
ISL1 in the SAN is unknown because of the early death of
ISL1-deficient mice and the loss of cardiomyocytes derived from ISL1 progenitor cells
[38]. Recent transcriptome studies have shown significant changes in SAN-related genes when
ISL1 is ablated in the SAN
[44].
ISL1 deficiency in mice leads to the downregulation of
TBX3,
SHOX2, and
BMP4, which are key regulators of SAN development, as well as
HCN4, HCN1, and
CACNA1G, which are responsible for the functioning of ion channels of the SAN
[44]. However, expressions of
NPPA, PITX2, NKX2.5, GJA1, GJA5, and
SCN5A were upregulated. Their study highlighted the central role of
ISL1 in establishing the PC gene program
[44].
SHOX2 plays a key role in PC fate and can activate the pacemaker gene program (
TBX3, ISL1, and
HCN4) by repressing
NKX2-5 and the working myocardial gene program
[38]. Hoffmann et al. performed transcriptome profiling of SHOX2
+/+ vs. SHOX2
−/− ESC-derived SANLPCs via RNA-seq to explore
SHOX2 pathways involved in pacemaker differentiation
[46]. Interestingly, some SAN-related genes (
TBX3, ISL1, HCN4, and
NKX2-5) were not affected by
SHOX2 deletion. On combining this finding with data obtained by Vedantham et al.
[44], several
SHOX2 target genes were discovered, namely,
CAV1, FKBP10, IGFBP5, MCF2l, and
NR2F2, which were validated in mouse and zebrafish models
[46]. Their study provided new insights into the transcriptional regulation of
SHOX2 during pacemaker development and function (
Figure 2).
The homeobox transcription factor
PITX2 is a laterality gene responsible for establishing the right- and left-body axes, asymmetric gene expression, and organ morphogenesis
[47]. Bilateral or ectopic SAN can be found in
PITX2-deficient embryos and may contribute to atrial fibrillation in adult animals with reduced
PITX2 expression
[48][49]. Single-cell transcriptomics has revealed the role of
PITX2 in cardiac development and left–right cellular specification
[50]. The CM-RA1 cluster with a SAN transcriptional signature was more abundant in
PITX2 mutants than that in controls
[50], highlighting the inhibitory effect of
PITX2 on SAN development.
In addition to the above TFs that affect SAN development and function, next-generation sequencing was recently used to identify other unexplored novel TFs in human SAN and RA and to predict interactions between key TFs and genes involved in pacemaker mechanisms
[51]. In the adult SAN, many new TFs were highly expressed (e.g.,
FOXD3, DLX2, PHOX2B, VENTX, and
SOX2). However, the role of these TFs in the adult mammalian SAN is unknown and should be explored in the future (
Figure 2).
A unique set of TFs, which are enriched in PCs, act as activators and repressors and interact with each other to determine the fate of the SAN
[38]. Some TFs have been used to direct stem cell differentiation toward SANLPCs
[52]; however, the core TF set requires further exploration to obtain highly pure SANLPCs.
3.2. Specific Molecular Markers of the SAN
In addition to the enrichment of established nodal genes, scRNA-seq revealed a host of significant novel genes not previously reported to be involved in SAN development or function, including insulin growth factor binding protein 5 (
IGFBP5)
[10], SPARC-related modular calcium-binding protein 2 (
SMOC2)
[10][53], neurotrimin (
NTM)
[10], copine 5 (
CPNE5)
[10], regulator of G-protein signaling type 6 (
RGS6)
[10][54], arachidonate 8-lipoxygenase (
ALOX8)
[53], sodium–hydrogen exchange regulatory cofactor 2 (
SLC9A3R2)
[53], and fibronectin leucine-rich transmembrane protein 3 (
FLRT3)
[55] (
Figure 2).
Van Eif et al. used CRISPR/Cas9 to produce SMOC2 frameshift mutation mice and verify the effect of
SMOC2 on SAN function. In vivo and in vitro experiments revealed that
SMOC2 inactivation had little effect on cardiac electrophysiology
[53]. Interestingly,
SMOC2 has recently been reported as a new SAN marker. The gene was subsequently validated via immunostaining or fluorescence in situ hybridization, showing
SMOC2 enrichment within the SAN compared with the surrounding AM
[10]. In a recent study, scRNA-seq of single cells from the SAN of different mammals identified a species-conserved potential SAN marker,
VSNL1 (a member of the visinin/recoverin subfamily of neuronal calcium sensor proteins), which is abundantly expressed in the SAN but barely expressed in the AM or ventricles
[56]. Although its function is unknown,
VSNL1 expression has been detected in the venous pole region of the developing heart
[57]. Moreover,
VSNL1 deficiency reduces the heart rate in human-induced pluripotent stem cell-derived cardiomyocytes and mice
[56].
FLRT3, a cell-autonomous regulator of the adherens junction of PCs, can mediate the transmission of electrical activity by regulating gap junctions
[55].