Characterization of diverse populations of sinoatrial node cells and their proliferation potential at single nucleus resolution

BACKGROUND: Each heartbeat is initiated in the sinoatrial node (SAN), and although a recent study (GSE130710) using single nucleus RNA-seq had discovered different populations of cell types within SAN tissue, the distinct potential functions of these cell types have not been delineated. METHODS: To infer some special potential functions of different SAN cell clusters, we applied principal component analysis (PCA), t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) to the GSE130710 dataset to reduce dimensions, followed by Pseudotime trajectory and AUCell analyses, ANOVA and Hurdle statistical models, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichments to determine functional potential of cell types. Nuclear EdU immuno-labeling of SAN tissue confirmed cell type proliferation. FINDINGS: We identified elements of a coupled clock system known to drive SAN cell pacemaking within the GSE130710 sinus node myocyte cluster, which, surprisingly, manifested signals of suppressed fatty acid and nitrogen metabolism and reduced immune gene expression. Proliferation signaling was enriched in endocardial, epicardial, epithelial cells, and macrophages, in which, fatty acid and nitrogen metabolic signals were also suppressed, but immune signaling was enhanced. EdU labeling was rare in pacemaker cells but was robust in interstitial cells. INTERPRETATION: Pacemaker cells that initiate each heartbeat manifest suppressed fatty acid and nitrogen metabolism and limited immune signaling and proliferation potential. In contrast, other populations of SAN cells not directly involved in the initiation of heartbeats, manifest robust proliferation and immune potential, likely to ensure an environment required to sustain healthy SAN tissue pacemaker function.