The SAN is a subtle biological system that comprises multiple subdomains and cell types with distinct functions, including pacemaker cells (PCs) and several non-PCs
[3]. Some previous studies that used genetically deficient mice (
TBX3, TBX5, TBX18, SHOX2, NKX2-5, and
PITX2) provided initial insights into SAN development and function
[4,5,6,7,8][4][5][6][7][8]. Unfortunately, these studies are often limited by the small size and heterogeneous organization of the SAN, providing insufficient knowledge on gene regulation of the SAN
[9]. Bulk RNA sequencing (RNA-seq) provides transcriptome data for the study of the SAN, further advancing SAN research. However, given the lack of specific molecular markers and the small size of the SAN, transcriptome data are often distorted and affected by other cells
[10]. Despite significant progress in SAN research,
ourthe understanding of its cellular features, differentiation trajectories, and regulatory mechanisms is incomplete. Some obstacles in investigating the molecular features of the SAN include (1) the low total number of PCs in the heart, (2) the complex three-dimensional anatomy of the SAN, (3) difficulties in separating PCs from the surrounding myocardium, (4) lack of specific molecular markers, and (5) significant cell heterogeneity in SAN tissues
[10,11][10][11].
2. SAN Structure
The human SAN is a compact, slightly elongated, banana-shaped structure located at the junction of the lateral superior vena cava (SVC) and right atrium (RA). Functional and structural mapping of the human SAN has revealed that it can be distinguished into head, center, and tail compartments. The head is usually closer to the epicardium, and the tail is typically tilted toward the endocardium
[15,16,17,18,19][13][14][15][16][17]. Different regions of the SAN pacemaker clusters exhibit diverse functional and molecular features, which are activated through distinct cellular mechanisms
[20,21,22][18][19][20] (
Figure 1A).
In mice, the SAN artery, autonomic nerve fibers, telocytes, monocytes, macrophages, adipocytes, and fibroblasts were observed via electron microscopy
[23][21]. High populations of fibroblasts in the SAN maintain the source–sink balance and ensure normal electrical conduction
[3]. The heterogeneous structure of the SAN makes studying its function challenging. Development of the scRNA-seq technology has aided in examining these heterogeneous and rare cells.
2.1. Three-Dimensional Structure of the SAN
In adult mice, the SAN head/center, which is densely packed with clusters of PCs, is the leading pacemaker region, extends superiorly, and wraps the right SVC
[24][22]. The SAN also extends inferiorly toward the inferior vena cava (IVC), forming the tail
[25,26][23][24]. In mice, the SAN can be genetically divided into a
TBX18+/NKX2-5− head domain and a
TBX18−/NKX2-5+ tail domain
[7] (
Figure 1B,C).
TBX18 deletion in the SAN head of mice has been shown to result in a normal heart rate
[7], whereas
SHOX2 ablation in the SAN tail has led to severe SSS
[27][25]. These findings highlight the different roles of different parts of the SAN in pacemaking and conduction.
Understanding the heterogeneity of the SAN complex is necessary to fully characterize its pacing and conduction functions. To further explore the complex 3D anatomy of the SAN, Goodyer et al. used scRNA-seq in E16.5 mice
[10]. The cells underwent unsupervised clustering via t-distributed stochastic neighbor embedding (t-SNE). Cluster 9 showed significant enrichment of nodal markers. Further experiments identified two distinct subclusters within cluster 9, i.e., compact SAN (cSAN) and transitional cells (Tz). cSAN showed high levels of established SAN markers with low or no known atrial myocardium (AM) marker expression
[10]. Further investigation of the cellular heterogeneity of the cSAN by Goodyer et al. revealed two distinct clusters in the cSAN, consistent with the previously recognized head and tail subdomains
[10]. Both clusters showed high expressions of the nodal markers
HCN4 and
SHOX2. The head cluster showed increased expression of
TBX18 and decreased expression of
NKX2-5, whereas the tail cluster showed decreased expression of
TBX18 and increased expression of
NKX2-5 [7], consistent with previous reports
[28][26]. A comparison of these two clusters revealed differentially expressed genes not previously reported.
IGFBP5, PDE1A, and
VSNL1 were abundant in the head, whereas
SMPX, ALDH1B1, and
SLC22A1 were abundant in the tail. These findings may enable the examination of different functional characteristics of the head and tail
[10].
To further explore the different functional and molecular bases of the SAN head and tail, Li et al. (2019) examined the action potentials (APs) of the SAN’s head and tail separately using a patch clamp and found that the SAN tail has unique APs that mediate the SAN head and A. By performing scRNA-seq assays on SAN regions highlighted by GFP isolated from E13.5 SHOX2
Cre/+ in R26RmTmG mice, they also verified that NKX2-5 plays an important role in the development and function of the SAN tail. They found four cardiomyocyte subpopulations, one of which expressed both SAN- and AM-specific genes, consistent with the properties of the SAN tail. They also revealed that the SAN tail has unique electrophysiological properties and transcriptomic features
[29][27].
According to Li et al.’s integrated NIOM-3D structural study, the human SAN is functionally insulated/discontinued from the atria by fibrotic tissue, fat, and myofibers, with the exception of 2–5 discrete sinoatrial conduction pathways (SACPs)
[19,30][17][28]. SACPs transmit electrical impulses to the RA
[31,32][29][30]. SACP structure includes tracts of myofibers with transitional cells on the SAN border that merge with atrial myobundles
[33][31].
The human fetal SAN contains predominantly PCs and limited collagen and transitional cells. From approximately 2 weeks after birth, Tz gradually increase in number and intertwine with the collagen framework
[34][32]. Importantly, Tz have been implicated in SSS
[33][31]. However, they remain poorly understood due to identification and isolation challenges. Goodyer et al. provided transcriptomic data for Tz
[10]. In their study, Tz expressed lower levels of nodal markers and were enriched for AM markers
GJA5 and
SCN5A [10]. However, no new markers were specific to Tz, and it remains uncertain as to how Tz enter the SAN architecture.
2.2. SAN Microenvironment
The unique microenvironment of the SAN maintains its functional characteristics. However, the formation and function of the SAN microenvironment remain poorly understood. To determine the cellular processes that contribute to the formation of the SAN microenvironment, Bressan et al. performed RNA-seq on chicken embryos earlier than E3.5 and found that factors associated with epithelial-to-mesenchymal transition and fibrosis were upregulated in the pacemaker region
[35][33]. Whole-mount in situ hybridization and microsurgery demonstrated that proepicardium-derived mesenchymal cells (including fibroblasts) contributed to the formation of the SAN microenvironment and that the integration of mesenchymal cells is essential to protecting PCs from the electrical “load” of the adjacent AM
[35][33]. However, how mesenchymal cells interact with PCs and regulate physiological functions in the SAN requires further investigation.
Recently, Chou et al. compared
TBX18-induced PCs and ventricular myocardium via RNA-seq and showed the predominance of glucose metabolism and glycolysis in the GO analysis
[36][34].
TBX18-induced PCs co-culturing with fibroblasts activated the integrin-dependent mitogen-activated protein kinase-E2F1 signal through cell–cell contact and induced Aldoc expression in PCs
[36][34]. Their results provided new insights into the link between fibroblasts and PCs and a new way to mediate SSS therapy by regulating Aldoc.
The human SAN comprises 35%–55% fibrotic tissue
[37][35]; cardiac diseases may further increase fibrosis within the SAN, potentially impairing automaticity and the conduction of electricity to the atria. Recently, Li et al. described the transcriptional characteristics of fibroblasts in the SANs of patients with heart failure but did not further elaborate on the connection between fibroblasts and PCs
[38][36].
In addition to fibroblasts, many other cell types exist in the SAN
[23][21]. The interactions among these cells and their roles in maintaining sinus rhythm have not been well elucidated. Single-cell transcriptomics is a promising tool for investigating cell–cell interaction networks
[39][37], which is also the direction of
ourthe research group in the future.
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
[26[24][28],
30], 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
[22][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
[40][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][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
[41][39], which underlies the funny current (If), an essential factor for the maintenance of sinus rhythm
[42][40]. In humans,
HCN4 is expressed in PCs and RA, and so cannot be used for identifying human PCs as the sole marker
[43][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
[44,45][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
[46,47][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
[40][38]. Recent transcriptome studies have shown significant changes in SAN-related genes when
ISL1 is ablated in the SAN
[46][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
[46][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
[46][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
[40][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
[48][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.
[46][44], several
SHOX2 target genes were discovered, namely,
CAV1, FKBP10, IGFBP5, MCF2l, and
NR2F2, which were validated in mouse and zebrafish models
[48][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
[49][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
[50,51][48][49]. Single-cell transcriptomics has revealed the role of
PITX2 in cardiac development and left–right cellular specification
[52][50]. The CM-RA1 cluster with a SAN transcriptional signature was more abundant in
PITX2 mutants than that in controls
[52][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
[53][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
[40][38]. Some TFs have been used to direct stem cell differentiation toward SANLPCs
[54][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[10][53],
55], neurotrimin (
NTM)
[10], copine 5 (
CPNE5)
[10], regulator of G-protein signaling type 6 (
RGS6)
[10[10][54],
56], arachidonate 8-lipoxygenase (
ALOX8)
[55][53], sodium–hydrogen exchange regulatory cofactor 2 (
SLC9A3R2)
[55][53], and fibronectin leucine-rich transmembrane protein 3 (
FLRT3)
[57][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
[55][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
[14][56]. Although its function is unknown,
VSNL1 expression has been detected in the venous pole region of the developing heart
[58][57]. Moreover,
VSNL1 deficiency reduces the heart rate in human-induced pluripotent stem cell-derived cardiomyocytes and mice
[14][56].
FLRT3, a cell-autonomous regulator of the adherens junction of PCs, can mediate the transmission of electrical activity by regulating gap junctions
[57][55].