Human Sinoatrial Node Research Articles

Profiling Reduced Expression of Contractile and Mitochondrial mRNAs in the Human Sinoatrial Node vs. Right Atrium and Predicting Their Suppressed Expression by Transcription Factors and/or microRNAs.

(1) Background: The sinus node (SN) is the main pacemaker of the heart. It is characterized by pacemaker cells that lack mitochondria and contractile elements. We investigated the possibility that transcription factors (TFs) and microRNAs (miRs) present in the SN can regulate gene expression that affects SN morphology and function. (2) Methods: From human next-generation sequencing data, a list of mRNAs that are expressed at lower levels in the SN compared with the right atrium (RA) was compiled. The mRNAs were then classified into contractile, mitochondrial or glycogen mRNAs using bioinformatic software, RStudio and Ingenuity Pathway Analysis. The mRNAs were combined with TFs and miRs to predict their interactions. (3) Results: From a compilation of the 1357 mRNAs, 280 contractile mRNAs and 198 mitochondrial mRNAs were identified to be expressed at lower levels in the SN compared with RA. TFs and miRs were shown to interact with contractile and mitochondrial function-related mRNAs. (4) Conclusions: In human SN, TFs (MYCN, SOX2, NUPR1 and PRDM16) mainly regulate mitochondrial mRNAs (COX5A, SLC25A11 and NDUFA8), while miRs (miR-153-3p, miR-654-5p, miR-10a-5p and miR-215-5p) mainly regulate contractile mRNAs (RYR2, CAMK2A and PRKAR1A). TF and miR-mRNA interactions provide a further understanding of the complex molecular makeup of the SN and potential therapeutic targets for cardiovascular treatments.

Caveolar Compartmentalization of Pacemaker Signaling is Required for Stable Rhythmicity of Sinus Nodal Cells and is Disrupted in Heart Failure.

Heart rhythm relies on complex interactions between electrogenic membrane proteins and intracellular Ca2+ signaling in sinoatrial node (SAN) myocytes; however, mechanisms underlying the functional organization of proteins involved in SAN pacemaking and its structural foundation remain elusive. Caveolae are nanoscale, plasma membrane pits that compartmentalize various ion channels and transporters, including those involved in SAN pacemaking, via binding with the caveolin-3 scaffolding protein, but the precise role of caveolae in cardiac pacemaker function is unknown. Our objective was to determine the role of caveolae in SAN pacemaking and dysfunction (SND). Biochemical co-purification, in vivo electrocardiogram monitoring, ex vivo optical mapping, in vitro confocal Ca2+ imaging, and immunofluorescent and electron microscopy analyses were performed in wild type, cardiac-specific caveolin-3 knockout, and 8-weeks post-myocardial infarction heart failure (HF) mice. SAN tissue samples from donor human hearts were used for biochemical studies. We utilized a novel 3-dimensional single SAN cell mathematical model to determine the functional outcomes of protein nanodomain-specific localization and redistribution in SAN pacemaking. In both mouse and human SANs, caveolae compartmentalized HCN4, Cav1.2, Cav1.3, Cav3.1 and NCX1 proteins within discrete pacemaker signalosomes via direct association with caveolin-3. This compartmentalization positioned electrogenic sarcolemmal proteins near the subsarcolemmal sarcoplasmic reticulum (SR) membrane and ensured fast and robust activation of NCX1 by subsarcolemmal local SR Ca2+ release events (LCRs), which diffuse across ~15-nm subsarcolemmal cleft. Disruption of caveolae led to the development of SND via suppression of pacemaker automaticity through a 50% decrease of the L-type Ca2+ current, a negative shift of the HCN current (I f) activation curve, and a 40% reduction of Na+/Ca2+-exchanger function, along with ~2.3-times widening of the sarcolemma-SR distance. These changes significantly decreased the SAN depolarizing force, both during diastolic depolarization and upstroke phase, leading to bradycardia, sinus pauses, recurrent development of SAN quiescence, and significant increase in heart rate lability. Computational modeling, supported by biochemical studies, identified NCX1 redistribution to extra-caveolar membrane as the primary mechanism of SAN pauses and quiescence due to the impaired ability of NCX1 to be effectively activated by LCRs and trigger action potentials. HF remodeling mirrored caveolae disruption leading to NCX1-LCR uncoupling and SND. SAN pacemaking is driven by complex protein interactions within a nanoscale caveolar pacemaker signalosome. Disruption of caveolae leads to SND, potentially demonstrating a new dimension of SAN remodeling and providing a newly recognized target for therapy.

Differentiation of Sinoatrial-like Cardiomyocytes as a Biological Pacemaker Model.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are widely used for disease modeling and pharmacological screening. However, their application has mainly focused on inherited cardiopathies affecting ventricular cardiomyocytes, leading to extensive knowledge on generating ventricular-like hiPSC-CMs. Electronic pacemakers, despite their utility, have significant disadvantages, including lack of hormonal responsiveness, infection risk, limited battery life, and inability to adapt to changes in heart size. Therefore, developing an in vitro multiscale model of the human sinoatrial node (SAN) pacemaker using hiPSC-CM and SAN-like cardiomyocyte differentiation protocols is essential. This would enhance the understanding of SAN-related pathologies and support targeted therapies. Generating SAN-like cardiomyocytes offers the potential for biological pacemakers and specialized conduction tissues, promising significant benefits for patients with conduction system defects. This review focuses on arrythmias related to pacemaker dysfunction, examining protocols' advantages and drawbacks for generating SAN-like cardiomyocytes from hESCs/hiPSCs, and discussing therapeutic approaches involving their engraftment in animal models.

Abstract We045: Development of SA Node Organoids: A Multicellular Approach to Biological Pacemaking

The rising prevalence of sinoatrial (SA) node dysfunction, attributed to the aging population, is emerging as a critical health issue. The SA node, a natural organoid within the heart, comprises the pacemaker cells embedded within a connective tissue matrix, insulated electrically from the surrounding atrial muscle cells. This isolated architecture, along with a heterogeneous cellular composition, is crucial for the SA node's robust pacemaking and electrical conduction. Developing methods to reconstruct this crucial heterogeneous structure could offer a long-term solution for individuals with SA node dysfunction. However, the use of current biological pacemakers, which typically involve a single cell type, has often led to unstable heart rhythms during one-month in vivo integration, challenging their long-term clinical application as biological pacemakers. To address existing SA node model limitations, we developed a SA node organoid that reproduces the human SA node's multicellular tissue structure. Initially, we developed a protocol for generating cardiac SA node organoids by employing a lentivirus to deliver dox-inducible genes for key pacemaker transcription factors, specifically TBX-18 or Shox-2. We mixed wild-type human-induced pluripotent stem cells (hiPSCs) with those modified to express TBX-18 or Shox-2. These cell mixtures were aggregated, embedded in a diluted Matrigel, and then exposed to a cardiac differentiation protocol, activating TBX-18 or Shox-2 at specific times (D1 or D7) via Dox administration. By D9, these organoids began to exhibit contraction, which lasted beyond 60 days, with notable upregulation of two transcription factors and pacemaker-specific genes following Dox treatment. Subsequently, we integrated the SA node organoid into 2D hiPSC-derived atrial muscle tissue constructs, maintaining geometric isolation within the cardiac tissue to enhance source-sink mismatch. This isolation allowed the SA node organoid to successfully trigger activity in adjacent dormant cardiomyocytes, preserving the rhythm for over 60 days. These findings could redefine the role of tissue architecture and cell diversity in the SA node's pacemaking function, potentially leading to a high-fidelity, tissue-level biological pacemaker as a new treatment strategy for SA node dysfunctions.

Morphology of human sinoatrial node and its surrounding right atrial muscle in the global obesity pandemic-does fat matter?

The sinus node (SN) is the main pacemaker site of the heart, located in the upper right atrium at the junction of the superior vena cava and right atrium. The precise morphology of the SN in the human heart remains relatively unclear especially the SN microscopical anatomy in the hearts of aged and obese individuals. In this study, the histology of the SN with surrounding right atrial (RA) muscle was analyzed from young non-obese, aged non-obese, aged obese and young obese individuals. The impacts of aging and obesity on fibrosis, apoptosis and cellular hypertrophy were investigated in the SN and RA. Moreover, the impact of obesity on P wave morphology in ECG was also analyzed to determine the speed and conduction of the impulse generated by the SN. Human SN/RA specimens were dissected from 23 post-mortem hearts (preserved in 4% formaldehyde solution), under Polish local ethical rules. The SN/RA tissue blocks were embedded in paraffin and histologically stained with Masson's Trichrome. High and low-magnification images were taken, and analysis was done for appropriate statistical tests on Prism (GraphPad, USA). 12-lead ECGs from 14 patients under Polish local ethical rules were obtained. The P wave morphologies from lead II, lead III and lead aVF were analyzed. Compared to the surrounding RA, the SN in all four groups has significantly more connective tissue (P ≤ 0.05) (young non-obese individuals, aged non-obese individuals, aged obese individuals and young obese individuals) and significantly smaller nodal cells (P ≤ 0.05) (young non-obese individuals, aged non-obese individuals, aged obese individuals, young obese individuals). In aging, overall, there was a significant increase in fibrosis, apoptosis, and cellular hypertrophy in the SN (P ≤ 0.05) and RA (P ≤ 0.05). Obesity did not further exacerbate fibrosis but caused a further increase in cellular hypertrophy (SN P ≤ 0.05, RA P ≤ 0.05), especially in young obese individuals. However, there was more infiltrating fat within the SN and RA bundles in obesity. Compared to the young non-obese individuals, the young obese individuals showed decreased P wave amplitude and P wave slope in aVF lead. Aging and obesity are two risk factors for extensive fibrosis and cellular hypertrophy in SN and RA. Obesity exacerbates the morphological alterations, especially hypertrophy of nodal and atrial myocytes. These morphological alterations might lead to functional alterations and eventually cause cardiovascular diseases, such as SN dysfunction, atrial fibrillation, bradycardia, and heart failure.

Mechanistic insight into the functional role of human sinoatrial node conduction pathways and pacemaker compartments heterogeneity: A computer model analysis.

The sinoatrial node (SAN), the primary pacemaker of the heart, is responsible for the initiation and robust regulation of sinus rhythm. 3D mapping studies of the ex-vivo human heart suggested that the robust regulation of sinus rhythm relies on specialized fibrotically-insulated pacemaker compartments (head, center and tail) with heterogeneous expressions of key ion channels and receptors. They also revealed up to five sinoatrial conduction pathways (SACPs), which electrically connect the SAN with neighboring right atrium (RA). To elucidate the role of these structural-molecular factors in the functional robustness of human SAN, we developed comprehensive biophysical computer models of the SAN based on 3D structural, functional and molecular mapping of ex-vivo human hearts. Our key finding is that the electrical insulation of the SAN except SACPs, the heterogeneous expression of If, INa currents and adenosine A1 receptors (A1R) across SAN pacemaker-conduction compartments are required to experimentally reproduce observed SAN activation patterns and important phenomena such as shifts of the leading pacemaker and preferential SACP. In particular, we found that the insulating border between the SAN and RA, is required for robust SAN function and protection from SAN arrest during adenosine challenge. The heterogeneity in the expression of A1R within the human SAN compartments underlies the direction of pacemaker shift and preferential SACPs in the presence of adenosine. Alterations of INa current and fibrotic remodelling in SACPs can significantly modulate SAN conduction and shift the preferential SACP/exit from SAN. Finally, we show that disease-induced fibrotic remodeling, INa suppression or increased adenosine make the human SAN vulnerable to pacing-induced exit blocks and reentrant arrhythmia. In summary, our computer model recapitulates the structural and functional features of the human SAN and can be a valuable tool for investigating mechanisms of SAN automaticity and conduction as well as SAN arrhythmia mechanisms under different pathophysiological conditions.

A gain-of-function HCN4 mutant in the HCN domain is responsible for inappropriate sinus tachycardia in a Spanish family.

In a family with inappropriate sinus tachycardia (IST), we identified a mutation (p.V240M) of the hyperpolarization-activated cyclic nucleotide-gated type 4 (HCN4) channel, which contributes to the pacemaker current (If) in human sinoatrial node cells. Here, we clinically study fifteen family members and functionally analyze the p.V240M variant. Macroscopic (IHCN4) and single-channel currents were recorded using patch-clamp in cells expressing human native (WT) and/or p.V240M HCN4 channels. All p.V240M mutation carriers exhibited IST that was accompanied by cardiomyopathy in adults. IHCN4 generated by p.V240M channels either alone or in combination with WT was significantly greater than that generated by WT channels alone. The variant, which lies in the N-terminal HCN domain, increased the single-channel conductance and opening frequency and probability of HCN4 channels. Conversely, it did not modify the channel sensitivity for cAMP and ivabradine or the level of expression at the membrane. Treatment with ivabradine based on functional data reversed the IST and the cardiomyopathy of the carriers. In computer simulations, the p.V240M gain-of-function variant increases If and beating rate and thus explains the IST of the carriers. The results demonstrate the importance of the unique HCN domain in HCN4, which stabilizes the channels in the closed state.

Abstract P1125: Development Of A Multicellular Sinoatrial Node Organoid For Modeling Robust Pacemaking And Conduction

The prevalence and healthcare burden of sinoatrial (SA) node dysfunction have been continuously increasing due to the growth in the elderly population, emphasizing the need for more detailed studies of SA node functions to enable effective therapy for treating and preventing SA node dysfunction. However, current SA node or pacemaker models are limited to isolated single cell-type cells or cell clusters, leaving a gap in the ability to model autonomous cardiac contraction. In contrast, the SA node is a natural organoid with elaborate insulated architecture and heterogeneous cellular composition. Furthermore, current single cell-type pacemakers have worsened heart rhythm stability during one-month in vivo integration, limiting their application as clinically viable biological pacemakers that can generate robust pacemaking and conduction. To address the current limitations of SA node models, we developed a three-dimensional multicellular SA node organoid by reproducing the human SA node's multicellular tissue structure and fail-safe mechanisms. First, inspired by the SA node's partial electrical insulation, we structurally isolated a group of cardiomyocytes with a single exit pathway that allowed for electrical connection to neighboring cardiomyocytes. We found that the geometrically isolated cardiomyocytes successfully activated neighboring quiescent cardiomyocytes without any pacemaker-related gene expression, and they maintained the rhythm for three months. Second, the overexpression of HCN4, a key pacemaker channel, in these isolated cardiomyocytes further improved heart rhythm variability down to the millisecond level, suggesting the interplay between pacemaker gene expression and nodal architecture for long-term, high-fidelity pacing. Finally, we successfully fabricated and integrated the SA node organoid into our 3D pacemaker organoid-myocardium platform, and the integrated organoid dominated rhythm control in 3D cardiac muscle constructs for 30 days. These studies will help define if tissue-level architecture and multicellular compositions mediate SA node's robust pacemaking and conduction and may reveal a high-fidelity tissue-level biological pacemaker as a novel therapeutic strategy for SA node dysfunctions.

BS-452759-2 HIGH-RESOLUTION 3D MULTIMODALITY IMAGING TO DISTINGUISH INTRAMURAL STRUCTURE OF HUMAN SINOATRIAL NODE VS EPI AND ENDO ATRIAL EXITS

Defining and targeting the human sinoatrial node (SAN) is challenging since clinical electrode mapping can visualize only surface epi/endocardial (Epi/Endo) exits of SAN activation through conduction pathways and not the leading pacemaker within the 3D intramural structure. To distinguish the functional-structural features of the human SAN pacemaker vs Epi-Endo exits by 3D multimodality imaging. Simultaneous EpiEndo dual-sided near infrared transmural optical mapping (NIOM, 300-500μm2 resolution) was conducted on coronary perfused explanted human atria (n=5, 45±20 y.o.). SAN pacemaker and exit sites were defined by NIOM activation maps. Histology-validated contrast-enhanced MRI (CE-MRI, ∼100μm3 resolution) was used to quantitatively analyze the 3D atrial structure of the SAN pacemaker vs exit sites, including wall thickness, transmural fiber orientation and fibrosis content in the sub-Epi, intramural, and sub-Endo layers. All ex-vivo preparations had stable sinus rhythm (SR, 77±14 bpm) at baseline with leading pacemakers in the intramural SANpacemaker compartments. In total, 24 SR recordings at baseline and autonomic stimulations) were analyzed. NIOM detected 29 exits on Epi and 32 exits on Endo. The majority of Epi or Endo earliest exits (88%) were distributed on the CT while secondary exits were located along the CT border. The most distant superior, middle, and inferior exits were 14.1±2.4mm, 8.8±6.3mm, and 12.3±3.1mm away from the SAN pacemakers respectively. 3D quantitative structural analysis revealed that transmural wall thickness is significantly less across the SAN pacemaker than atrial exit sites, but the fibrosis content of intramural and sub-Endo layers is significantly higher in the SAN than in exit sites. There are transmural myofibers that were defined primarily at the exit sites along CT but not within SAN. This study provides a quantitative atlas of the 3D architecture of the human SAN, and reveals that SAN pacemakers can be distinguished from both Epi/Endo exits by wall thickness, intramural/sub-Endo fibrosis content, and transmural fiber orientations. The new 3D transmural imaging analysis approach can be integrated with clinical electrode mapping, which may help to define SAN pacemaker and guide patient-specific targeted ablation of atrial SAN and CT arrhythmias.

Human Sinoatrial Node Pacemaker Activity: Role of the Slow Component of the Delayed Rectifier K+ Current, IKs.

The pacemaker activity of the sinoatrial node (SAN) has been studied extensively in animal species but is virtually unexplored in humans. Here we assess the role of the slowly activating component of the delayed rectifier K+ current (IKs) in human SAN pacemaker activity and its dependence on heart rate and β-adrenergic stimulation. HEK-293 cells were transiently transfected with wild-type KCNQ1 and KCNE1 cDNA, encoding the α- and β-subunits of the IKs channel, respectively. KCNQ1/KCNE1 currents were recorded both during a traditional voltage clamp and during an action potential (AP) clamp with human SAN-like APs. Forskolin (10 µmol/L) was used to increase the intracellular cAMP level, thus mimicking β-adrenergic stimulation. The experimentally observed effects were evaluated in the Fabbri-Severi computer model of an isolated human SAN cell. Transfected HEK-293 cells displayed large IKs-like outward currents in response to depolarizing voltage clamp steps. Forskolin significantly increased the current density and significantly shifted the half-maximal activation voltage towards more negative potentials. Furthermore, forskolin significantly accelerated activation without affecting the rate of deactivation. During an AP clamp, the KCNQ1/KCNE1 current was substantial during the AP phase, but relatively small during diastolic depolarization. In the presence of forskolin, the KCNQ1/KCNE1 current during both the AP phase and diastolic depolarization increased, resulting in a clearly active KCNQ1/KCNE1 current during diastolic depolarization, particularly at shorter cycle lengths. Computer simulations demonstrated that IKs reduces the intrinsic beating rate through its slowing effect on diastolic depolarization at all levels of autonomic tone and that gain-of-function mutations in KCNQ1 may exert a marked bradycardic effect during vagal tone. In conclusion, IKs is active during human SAN pacemaker activity and has a strong dependence on heart rate and cAMP level, with a prominent role at all levels of autonomic tone.

Development of a new histological identification method of human sinoatrial node suitable for immunohistochemical study.

Histological identification of the human sinoatrial node (SAN) remains a challenge. Conventional identification methods, such as Lev's method, have certain limitations. The aim of our study was to develop a new histological identification method that could properly identify the sinoatrial node, applicable to the immunohistochemical study of intra-nodal structures. Thirty-nine human autopsied hearts were included in this study. The cases included 23 men and 16 women ranging in age from 20 to 99years. The sinoatrial area from eight control samples was cut in the vertical section using the conventional Lev's method. In our new method, called the "En face one-block method," the sinoatrial node was cut in "En face" at the junction of the right border of the right appendage and superior vena cava, placed in one long cassette, and serially cut using a microtome. Immunostaining was performed using primary antibodies against CD31, podoplanin (D2-40), S-100, and other proteins. The average area of the SAN on the slide glass in our new method was 32.2 mm2, which was significantly larger than that (3.59 mm2) of the control samples by Lev's method. The SAN area was positively correlated with age (r = 0.357; p = 0.026), especially in women (r = 0.626; p = 0.0095). The SAN group had significantly lower percentage of CD31-positive blood capillaries, higher percentage of podoplanin-positive lymphatic channels, and S-100-positive peripheral nerves. We successfully developed a novel cutting method applicable to immunohistochemical studies, with which we could provide a bird's-eye view of the sinoatrial nodes.

Coupling and heterogeneity modulate pacemaking capability in healthy and diseased two-dimensional sinoatrial node tissue models.

Both experimental and modeling studies have attempted to determine mechanisms by which a small anatomical region, such as the sinoatrial node (SAN), can robustly drive electrical activity in the human heart. However, despite many advances from prior research, important questions remain unanswered. This study aimed to investigate, through mathematical modeling, the roles of intercellular coupling and cellular heterogeneity in synchronization and pacemaking within the healthy and diseased SAN. In a multicellular computational model of a monolayer of either human or rabbit SAN cells, simulations revealed that heterogenous cells synchronize their discharge frequency into a unique beating rhythm across a wide range of heterogeneity and intercellular coupling values. However, an unanticipated behavior appeared under pathological conditions where perturbation of ionic currents led to reduced excitability. Under these conditions, an intermediate range of intercellular coupling (900-4000 MΩ) was beneficial to SAN automaticity, enabling a very small portion of tissue (3.4%) to drive propagation, with propagation failure occurring at both lower and higher resistances. This protective effect of intercellular coupling and heterogeneity, seen in both human and rabbit tissues, highlights the remarkable resilience of the SAN. Overall, the model presented in this work allowed insight into how spontaneous beating of the SAN tissue may be preserved in the face of perturbations that can cause individual cells to lose automaticity. The simulations suggest that certain degrees of gap junctional coupling protect the SAN from ionic perturbations that can be caused by drugs or mutations.

Abstract 15525: Resolving Human Sinoatrial Node Leading Pacemaker vs Atrial Early Activation Sites by 3D High Resolution Imaging to Guide Ablation Treatment of Atrial Arrhythmia

Introduction: Surgical/ablation treatment for atrial arrhythmias such as anatomic sinoatrial node (SAN) tachycardia and atrial fibrillation either target or preserve the SAN. However, defining SAN is challenging because clinical electrophysiological (EP) mapping can visualize only exits of SAN activation as early activation sites (EAS) and not the leading pacemaker (LP) within the 3D intramural SAN. Methods: High resolution (300-900μm 2 ) including epi/endocardial (Epi/Endo) dual sided near infrared optical mapping (NIOM) was conducted on coronary perfused explanted human atria (n=26, 19-69 y.o.). SAN was defined by optical action potential morphologies as the region of slow diastolic depolarization (slow upstroke) preceding atrial excitation. Serial histology and contrast enhanced MRI (CE-MRI, 100μm 3 resolution) were used to define the 3D fibrotic structure of the SAN pacemaker complex. Results: During sinus rhythm (SR) (83±22 bpm), the LP was primarily located in the SAN center. Electrical impulses exit from SAN to atria through 1-2 sinoatrial conduction pathways (SACP), leading to discrete EAS along crista terminalis (CT), preferentially from lateral superior/middle SACPs. The distance between SAN LP and EAS varied from 3.5-23 mm. Heterogeneous atrial wall thickness, fibrosis content and myofiber orientation along the SAN and CT (5-15mm thick) regions led to complex intramural conduction from SAN LP to EAS and substantial Epi-Endo activation dyssynchrony. Conclusions: EAS visualized on both Epi/Endo mapping mainly distributed along thick CT, but not on the surface projection of SAN, due to intramural fiber orientation of preferential SACPs. The higher fibrotic content in the human SAN than CT detected by CE-MRI, integrated with EP mapping can be helpful to accurately define SAN structure and pathways for reentrant SAN arrhythmias ablations.

Characterization of the pace-and-drive capacity of the human sinoatrial node: A 3D in silico study

The sinoatrial node (SAN) is a complex structure that spontaneously depolarizes rhythmically (“pacing”) and excites the surrounding non-automatic cardiac cells (“drive”) to initiate each heart beat. However, the mechanisms by which the SAN cells can activate the large and hyperpolarized surrounding cardiac tissue are incompletely understood. Experimental studies demonstrated the presence of an insulating border that separates the SAN from the hyperpolarizing influence of the surrounding myocardium, except at a discrete number of sinoatrial exit pathways (SEPs). We propose a highly detailed 3D model of the human SAN, including 3D SEPs to study the requirements for successful electrical activation of the primary pacemaking structure of the human heart. A total of 788 simulations investigate the ability of the SAN to pace and drive with different heterogeneous characteristics of the nodal tissue (gradient and mosaic models) and myocyte orientation. A sigmoidal distribution of the tissue conductivity combined with a mosaic model of SAN and atrial cells in the SEP was able to drive the right atrium (RA) at varying rates induced by gradual If block. Additionally, we investigated the influence of the SEPs by varying their number, length, and width. SEPs created a transition zone of transmembrane voltage and ionic currents to enable successful pace and drive. Unsuccessful simulations showed a hyperpolarized transmembrane voltage (−66 mV), which blocked the L-type channels and attenuated the sodium-calcium exchanger. The fiber direction influenced the SEPs that preferentially activated the crista terminalis (CT). The location of the leading pacemaker site (LPS) shifted toward the SEP-free areas. LPSs were located closer to the SEP-free areas (3.46 ± 1.42 mm), where the hyperpolarizing influence of the CT was reduced, compared with a larger distance from the LPS to the areas where SEPs were located (7.17± 0.98 mm). This study identified the geometrical and electrophysiological aspects of the 3D SAN-SEP-CT structure required for successful pace and drive in silico.

SHOX2 refines the identification of human sinoatrial nodal cell population in the invitro cardiac differentiation.

IntroductionDysfunction of the sinoatrial node (SAN) cells causes arrhythmias, and many patients require artificial cardiac pacemaker implantation. However, the mechanism of impaired SAN automaticity remains unknown, and the generation of human SAN cells in vitro may provide a platform for understanding the pathogenesis of SAN dysfunction. The short stature homeobox 2 (SHOX2) and hyperpolarization-activated cyclic nucleotide-gated cation channel 4 (HCN4) genes are specifically expressed in SAN cells and are important for SAN development and automaticity. In this study, we aimed to purify and characterize human SAN-like cells in vitro, using HCN4 and SHOX2 as SAN markers.MethodsWe developed an HCN4-EGFP/SHOX2-mCherry dual reporter cell line derived from human induced pluripotent stem cells (hiPSCs), and HCN4 and SHOX2 gene expressions were visualized using the fluorescent proteins EGFP and mCherry, respectively. The dual reporter cell line was established using an HCN4-EGFP bacterial artificial chromosome-based semi-knock-in system and a CRISPR-Cas9-dependent knock-in system with a SHOX2-mCherry targeting vector. Flow cytometry, RT-PCR, and whole-cell patch-clamp analyses were performed to identify SAN-like cells.ResultsFlow cytometry analysis and cell sorting isolated HCN4-EGFP single-positive (HCN4+/SHOX2-) and HCN4-EGFP/SHOX2-mCherry double-positive (HCN4+/SHOX2+) cells. RT-PCR analyses showed that SAN-related genes were enriched within the HCN4+/SHOX2+ cells. Further, electrophysiological analyses showed that approximately 70% of the HCN4+/SHOX2+ cells exhibited SAN-like electrophysiological characteristics, as defined by the action potential parameters of the maximum upstroke velocity and action potential duration.ConclusionsThe HCN4-EGFP/SHOX2-mCherry dual reporter hiPSC system developed in this study enabled the enrichment of SAN-like cells within a mixed HCN4+/SHOX2+ population of differentiating cardiac cells. This novel cell line is useful for the further enrichment of human SAN-like cells. It may contribute to regenerative medicine, for example, biological pacemakers, as well as testing for cardiotoxic and chronotropic actions of novel drug candidates.

SARS-CoV-2 Infection Induces Ferroptosis of Sinoatrial Node Pacemaker Cells.

Background:Increasing evidence suggests that cardiac arrhythmias are frequent clinical features of coronavirus disease 2019 (COVID-19). Sinus node damage may lead to bradycardia. However, it is challenging to explore human sinoatrial node (SAN) pathophysiology due to difficulty in isolating and culturing human SAN cells. Embryonic stem cells (ESCs) can be a source to derive human SAN-like pacemaker cells for disease modeling.Methods:We used both a hamster model and human ESC (hESC)–derived SAN-like pacemaker cells to explore the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on the pacemaker cells of the heart. In the hamster model, quantitative real-time polymerase chain reaction and immunostaining were used to detect viral RNA and protein, respectively. We then created a dual knock-in SHOX2:GFP;MYH6:mCherry hESC reporter line to establish a highly efficient strategy to derive functional human SAN-like pacemaker cells, which was further characterized by single-cell RNA sequencing. Following exposure to SARS-CoV-2, quantitative real-time polymerase chain reaction, immunostaining, and RNA sequencing were used to confirm infection and determine the host response of hESC-SAN–like pacemaker cells. Finally, a high content chemical screen was performed to identify drugs that can inhibit SARS-CoV-2 infection, and block SARS-CoV-2–induced ferroptosis.Results:Viral RNA and spike protein were detected in SAN cells in the hearts of infected hamsters. We established an efficient strategy to derive from hESCs functional human SAN-like pacemaker cells, which express pacemaker markers and display SAN-like action potentials. Furthermore, SARS-CoV-2 infection causes dysfunction of human SAN-like pacemaker cells and induces ferroptosis. Two drug candidates, deferoxamine and imatinib, were identified from the high content screen, able to block SARS-CoV-2 infection and infection-associated ferroptosis.Conclusions:Using a hamster model, we showed that primary pacemaker cells in the heart can be infected by SARS-CoV-2. Infection of hESC-derived functional SAN-like pacemaker cells demonstrates ferroptosis as a potential mechanism for causing cardiac arrhythmias in patients with COVID-19. Finally, we identified candidate drugs that can protect the SAN cells from SARS-CoV-2 infection.

Mechanisms of Sinoatrial Node Dysfunction in Heart Failure With Preserved Ejection Fraction.

The ability to increase heart rate during exercise and other stressors is a key homeostatic feature of the sinoatrial node (SAN). When the physiological heart rate response is blunted, chronotropic incompetence limits exercise capacity, a common problem in patients with heart failure with preserved ejection fraction (HFpEF). Despite its clinical relevance, the mechanisms of chronotropic incompetence remain unknown. Dahl salt-sensitive rats fed a high-salt diet and C57Bl6 mice fed a high-fat diet and an inhibitor of constitutive nitric oxide synthase (Nω-nitro-L-arginine methyl ester [L-NAME]; 2-hit) were used as models of HFpEF. Myocardial infarction was created to induce HF with reduced ejection fraction. Rats and mice fed with a normal diet or those that had a sham surgery served as respective controls. A comprehensive characterization of SAN function and chronotropic response was conducted by in vivo, ex vivo, and single-cell electrophysiologic studies. RNA sequencing of SAN was performed to identify transcriptomic changes. Computational modeling of biophysically-detailed human HFpEF SAN was created. Rats with phenotypically-verified HFpEF exhibited limited chronotropic response associated with intrinsic SAN dysfunction, including impaired β-adrenergic responsiveness and an alternating leading pacemaker within the SAN. Prolonged SAN recovery time and reduced SAN sensitivity to isoproterenol were confirmed in the 2-hit mouse model. Adenosine challenge unmasked conduction blocks within the SAN, which were associated with structural remodeling. Chronotropic incompetence and SAN dysfunction were also found in rats with HF with reduced ejection fraction. Single-cell studies and transcriptomic profiling revealed HFpEF-related alterations in both the "membrane clock" (ion channels) and the "Ca2+ clock" (spontaneous Ca2+ release events). The physiologic impairments were reproduced in silico by empirically-constrained quantitative modeling of human SAN function. Chronotropic incompetence and SAN dysfunction were seen in both models of HF. We identified that intrinsic abnormalities of SAN structure and function underlie the chronotropic response in HFpEF.

Hypoglycaemia combined with mild hypokalaemia reduces the heart rate and causes abnormal pacemaker activity in a computational model of a human sinoatrial cell.

Low blood glucose, hypoglycaemia, has been implicated as a possible contributing factor to sudden cardiac death (SCD) in people with diabetes but it is challenging to investigate in clinical studies. We hypothesized the effects of hypoglycaemia on the sinoatrial node (SAN) in the heart to be a candidate mechanism and adapted a computational model of the human SAN action potential developed by Fabbri et al., to investigate the effects of hypoglycaemia on the pacemaker rate. Using Latin hypercube sampling, we combined the effects of low glucose (LG) on the human ether-a-go-go-related gene channel with reduced blood potassium, hypokalaemia, and added sympathetic and parasympathetic stimulus. We showed that hypoglycaemia on its own causes a small decrease in heart rate but there was also a marked decrease in heart rate when combined with hypokalaemia. The effect of the sympathetic stimulus was diminished, causing a smaller increase in heart rate, with LG and hypokalaemia compared to normoglycaemia. By contrast, the effect of the parasympathetic stimulus was enhanced, causing a greater decrease in heart rate. We therefore demonstrate a potential mechanistic explanation for hypoglycaemia-induced bradycardia and show that sinus arrest is a plausible mechanism for SCD in people with diabetes.

Altered microRNA and mRNA profiles during heart failure in the human sinoatrial node

Heart failure (HF) is frequently accompanied with the sinoatrial node (SAN) dysfunction, which causes tachy-brady arrhythmias and increased mortality. MicroRNA (miR) alterations are associated with HF progression. However, the transcriptome of HF human SAN, and its role in HF-associated remodeling of ion channels, transporters, and receptors responsible for SAN automaticity and conduction impairments is unknown. We conducted comprehensive high-throughput transcriptomic analysis of pure human SAN primary pacemaker tissue and neighboring right atrial tissue from human transplanted HF hearts (n = 10) and non-failing (nHF) donor hearts (n = 9), using next-generation sequencing. Overall, 47 miRs and 832 mRNAs related to multiple signaling pathways, including cardiac diseases, tachy-brady arrhythmias and fibrosis, were significantly altered in HF SAN. Of the altered miRs, 27 are predicted to regulate mRNAs of major ion channels and neurotransmitter receptors which are involved in SAN automaticity (e.g. HCN1, HCN4, SLC8A1) and intranodal conduction (e.g. SCN5A, SCN8A) or both (e.g. KCNJ3, KCNJ5). Luciferase reporter assays were used to validate interactions of miRs with predicted mRNA targets. In conclusion, our study provides a profile of altered miRs in HF human SAN, and a novel transcriptome blueprint to identify molecular targets for SAN dysfunction and arrhythmia treatments in HF.

Reproducibility study of the Fabbri et al. 2017 model of the human sinus node action potential

The sinoatrial node (SAN) is the natural pacemaker of the mammalian heart. It has been the subject of several mathematical studies, aimed at reproducing its electrical response under normal sinus rhythms, as well as under various conditions. Such studies were traditionally done using data from rabbit SAN cells. More recently, human SAN cell data have become available, resulting in the publication of a human SAN cell model (Fabbri et al., 2017), along with its CellML version. Here, we used the CellML file provided by the model authors, together with some SED-ML files and Python scripts that we created to reproduce the main results of the aforementioned modeling study. EDITOR'S NOTE (v2): this article and its OMEX archive are republished with technical changes made to the corresponding Python scripts to remove a run-time error message displayed when executing each simulation.