Shox2 is essential for the differentiation of cardiac pacemaker cells by repressing Nkx2-5

Odz4 upregulates SAN-specific genes to promote differentiation into cardiac pacemaker-like cells

Cardiac arrhythmias stemming from abnormal sinoatrial node (SAN) function can lead to sudden death. Developing a biological pacemaker device for treating sick sinus syndrome (SSS) could offer a potential cure. Understanding SAN differentiation is crucial, yet its regulatory mechanism remains unclear. We reanalyzed published RNA-seq data and identified Odz4 as a SAN-specific candidate. In situ hybridization revealed Odz4 expression in the cardiac crescent and throughout the cardiac conduction system (CCS). To assess the role of Odz4 in CCS differentiation, we utilized a Tet-Off inducible system for its intracellular domain (ICD). Embryonic bodies (EBs) exogenously expressing Odz4-ICD exhibited an increased propensity to develop into pacemaker-like cells with enhanced automaticity and upregulated expression of SAN-specific genes. CellChat and GO analyses unveiled SAN-specific enrichment of ligand-receptor sets, especially Ptn-Ncl, and extracellular matrix components in the group exogenously expressing Odz4-ICD. Our findings underscore the significance of Odz4 in SAN development and offer fresh insights into biological pacemaker establishment.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

A new paradigm for generating high-quality cardiac pacemaker cells from mouse pluripotent stem cells

Cardiac biological pacing (BP) is one of the future directions for bradyarrhythmias intervention. Currently, cardiac pacemaker cells (PCs) used for cardiac BP are mainly derived from pluripotent stem cells (PSCs). However, the production of high-quality cardiac PCs from PSCs remains a challenge. Here, we developed a cardiac PC differentiation strategy by adopting dual PC markers and simulating the developmental route of PCs. First, two PC markers, Shox2 and Hcn4, were selected to establish Shox2:EGFP; Hcn4:mCherry mouse PSC reporter line. Then, by stepwise guiding naïve PSCs to cardiac PCs following naïve to formative pluripotency transition and manipulating signaling pathways during cardiac PCs differentiation, we designed the FSK method that increased the yield of SHOX2+; HCN4+ cells with typical PC characteristics, which was 12 and 42 folds higher than that of the embryoid body (EB) and the monolayer M10 methods respectively. In addition, the in vitro cardiac PCs differentiation trajectory was mapped by single-cell RNA sequencing (scRNA-seq), which resembled in vivo PCs development, and ZFP503 was verified as a key regulator of cardiac PCs differentiation. These PSC-derived cardiac PCs have the potential to drive advances in cardiac BP technology, help with the understanding of PCs (patho)physiology, and benefit drug discovery for PC-related diseases as well.

Transcriptional regulation of the postnatal cardiac conduction system heterogeneity

The cardiac conduction system (CCS) is a network of specialized cardiomyocytes that coordinates electrical impulse generation and propagation for synchronized heart contractions. Although the components of the CCS, including the sinoatrial node, atrioventricular node, His bundle, bundle branches, and Purkinje fibers, were anatomically discovered more than 100 years ago, their molecular constituents and regulatory mechanisms remain incompletely understood. Here, we demonstrate the transcriptomic landscape of the postnatal mouse CCS at a single-cell resolution with spatial information. Integration of single-cell and spatial transcriptomics uncover region-specific markers and zonation patterns of expression. Network inference shows heterogeneous gene regulatory networks across the CCS. Notably, region-specific gene regulation is recapitulated in vitro using neonatal mouse atrial and ventricular myocytes overexpressing CCS-specific transcription factors, Tbx3 and/or Irx3. This finding is supported by ATAC-seq of different CCS regions, Tbx3 ChIP-seq, and Irx motifs. Overall, this study provides comprehensive molecular profiles of the postnatal CCS and elucidates gene regulatory mechanisms contributing to its heterogeneity.

Segregation of morphogenetic regulatory function of Shox2 from its cell fate guardian role in sinoatrial node development

Shox2 plays a vital role in the morphogenesis and physiological function of the sinoatrial node (SAN), the primary cardiac pacemaker, manifested by the formation of a hypoplastic SAN and failed differentiation of pacemaker cells in Shox2 mutants. Shox2 and Nkx2-5 are co-expressed in the developing SAN and regulate the fate of the pacemaker cells through a Shox2-Nkx2-5 antagonistic mechanism. Here we show that simultaneous inactivation of Nkx2-5 in the SAN of Shox2 mutants (dKO) rescued the pacemaking cell fate but not the hypoplastic defects, indicating uncoupling of SAN cell fate determination and morphogenesis. Single-cell RNA-seq revealed that the presumptive SAN cells of Shox2-/- mutants failed to activate pacemaking program but remained in a progenitor state preceding working myocardium, while both wildtype and dKO SAN cells displayed normal pacemaking cell fate with similar cellular state. Shox2 thus acts as a safeguard but not a determinant to ensure the pacemaking cell fate through the Shox2-Nkx2-5 antagonistic mechanism, which is segregated from its morphogenetic regulatory function in SAN development.

Unsupervised ensemble-based phenotyping enhances discoverability of genes related to left-ventricular morphology

Recent genome-wide association studies have successfully identified associations between genetic variants and simple cardiac morphological parameters derived from cardiac magnetic resonance images. However, the emergence of large databases, including genetic data linked to cardiac magnetic resonance facilitates the investigation of more nuanced patterns of cardiac shape variability than those studied so far. Here we propose a framework for gene discovery coined unsupervised phenotype ensembles. The unsupervised phenotype ensemble builds a redundant yet highly expressive representation by pooling a set of phenotypes learnt in an unsupervised manner, using deep learning models trained with different hyperparameters. These phenotypes are then analysed via genome-wide association studies, retaining only highly confident and stable associations across the ensemble. We applied our approach to the UK Biobank database to extract geometric features of the left ventricle from image-derived three-dimensional meshes. We demonstrate that our approach greatly improves the discoverability of genes that influence left ventricle shape, identifying 49 loci with study-wide significance and 25 with suggestive significance. We argue that our approach would enable more extensive discovery of gene associations with image-derived phenotypes for other organs or image modalities.

Canonical Wnt signaling directs the generation of functional human PSC-derived atrioventricular canal cardiomyocytes in bioprinted cardiac tissues

The creation of a functional 3D bioprinted human heart remains challenging, largely due to the lack of some crucial cardiac cell types, including the atrioventricular canal (AVC) cardiomyocytes, which are essential to slow down the electrical impulse between the atrium and ventricle. By utilizing single-cell RNA sequencing analysis and a 3D bioprinting technology, we discover that stage-specific activation of canonical Wnt signaling creates functional AVC cardiomyocytes derived from human pluripotent stem cells. These cardiomyocytes display morphological characteristics and express molecular markers of AVC cardiomyocytes, including transcription factors TBX2 and MSX2. When bioprinted in prefabricated cardiac tissues, these cardiomyocytes successfully delay the electrical impulse, demonstrating their capability of functioning as the AVC cardiomyocytes in vitro. Thus, these findings not only identify canonical Wnt signaling as a key regulator of the AVC cardiomyocyte differentiation in vitro, but, more importantly, provide a critical cellular source for the biofabrication of a functional human heart.

The cardiac conduction system: History, development, and disease

The heart is the first organ to form during embryonic development, establishing the circulatory infrastructure necessary to sustain life and enable downstream organogenesis. Critical to the heart's function is its ability to initiate and propagate electrical impulses that allow for the coordinated contraction and relaxation of its chambers, and thus, the movement of blood and nutrients. Several specialized structures within the heart, collectively known as the cardiac conduction system (CCS), are responsible for this phenomenon. In this review, we discuss the discovery and scientific history of the mammalian cardiac conduction system as well as the key genes and transcription factors implicated in the formation of its major structures. We also describe known human diseases related to CCS development and explore existing challenges in the clinical context.

Cellular heterogeneity of pluripotent stem cell-derived cardiomyocyte grafts is mechanistically linked to treatable arrhythmias

Preclinical data have confirmed that human pluripotent stem cell-derived cardiomyocytes (PSC-CMs) can remuscularize the injured or diseased heart, with several clinical trials now in planning or recruitment stages. However, because ventricular arrhythmias represent a complication following engraftment of intramyocardially injected PSC-CMs, it is necessary to provide treatment strategies to control or prevent engraftment arrhythmias (EAs). Here, we show in a porcine model of myocardial infarction and PSC-CM transplantation that EAs are mechanistically linked to cellular heterogeneity in the input PSC-CM and resultant graft. Specifically, we identify atrial and pacemaker-like cardiomyocytes as culprit arrhythmogenic subpopulations. Two unique surface marker signatures, signal regulatory protein α (SIRPA)+CD90-CD200+ and SIRPA+CD90-CD200-, identify arrhythmogenic and non-arrhythmogenic cardiomyocytes, respectively. Our data suggest that modifications to current PSC-CM-production and/or PSC-CM-selection protocols could potentially prevent EAs. We further show that pharmacologic and interventional anti-arrhythmic strategies can control and potentially abolish these arrhythmias.

Development of the Cardiac Conduction System

The electrical impulses that coordinate the sequential, rhythmic contractions of the atria and ventricles are initiated and tightly regulated by the specialized tissues of the cardiac conduction system. In the mature heart, these impulses are generated by the pacemaker cardiomyocytes of the sinoatrial node, propagated through the atria to the atrioventricular node where they are delayed and then rapidly propagated to the atrioventricular bundle, right and left bundle branches, and finally, the peripheral ventricular conduction system. Each of these specialized components arise by complex patterning events during embryonic development. This chapter addresses the origins and transcriptional networks and signaling pathways that drive the development and maintain the function of the cardiac conduction system.

The multi-lineage transcription factor ISL1 controls cardiomyocyte cell fate through interaction with NKX2.5

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription.

D-galactose causes sinoatrial node dysfunction: from phenotype to mechanism

With the population aging, age-related sinoatrial node dysfunction (SND) has been on the rise. Sinoatrial node (SAN) degeneration is an important factor for the age-related SND development. However, there is no suitable animal modeling method in this field. Here, we investigated whether D-galactose could induce SAN degeneration and explored the associated mechanism. In vivo, twelve C57BL/6 mice were divided into Control and D-galactose group to receive corresponding treatments. Senescence was confirmed by analyzing the hair and weight; cardiac function was evaluated through echocardiography, cerebral blood flux and serum-BNP; the SAN function was evaluated by electrocardiogram; fibrotic change was evaluated by Masson's trichrome staining and oxidative stress was assessed through DHE staining and serum indicators. Mechanism was verified through immunofluorescence-staining and Western blotting. In vitro, mouse-atrial-myocytes were treated with D-galactose, and edaravone was utilized as the ROS scavenger. Senescence, oxidative stress, proliferation ability and mechanism were verified through various methods, and intuitive evidence was obtained through electrophysiological assay. Finally, we concluded that D-galactose can be used to induce age-related SND, in which oxidative stress plays a key role, causing PITX2 ectopic expression and downregulates SHOX2 expression, then through the downstream GATA4/NKX2-5 axis, results in pacing-related ion channels dysfunction, and hence SND development.

Risk stratification of atrial fibrillation and stroke using single nucleotide polymorphism and circulating biomarkers

BACKGROUND

Atrial fibrillation (AF) is the most common sustained arrhythmia, and it causes a high rate of complications such as stroke. It is known that AF begins as paroxysmal form and gradually progresses to persistent form, and sometimes it is difficult to identify paroxysmal AF (PAF) before having stroke. The aim of this study is to evaluate the risk of PAF and stroke using genetic analysis and circulating biomarkers.

MATERIALS AND METHODS

A total of 600 adult subjects were enrolled (300 from PAF and control groups). Peripheral blood was drawn to identify the genetic variation and biomarkers. Ten single nucleotide polymorphisms (SNPs) were analyzed, and circulating cell-free DNA (cfDNA) was measured from plasma. Four microRNAs (miR-99a-5p, miR-192-5p, miR-214-3p, and miR-342-5p) were quantified in serum using quantitative RT-PCR.

RESULTS

Genotyping identified 4 single nucleotide polymorphisms (SNPs) that were significantly associated with AF (rs6817105, rs3807989, rs10824026, and rs2106261), and the genetic risk score using 4 SNPs showed the area under the curve (AUC) of 0.631. Circulating miRNAs and cfDNA did not show significant differences between PAF and control groups. The concentration of cfDNA was significantly higher in patients with a history of stroke, and the AUC was 0.950 to estimate the association with stroke.

CONCLUSION

The risk of AF could be assessed by genetic risk score. Furthermore, the risk of stroke might be evaluated by plasma cfDNA level.

Local tissue mechanics control cardiac pacemaker cell embryonic patterning

Cardiac pacemaker cells (CPCs) initiate the electric impulses that drive the rhythmic beating of the heart. CPCs reside in a heterogeneous, ECM-rich microenvironment termed the sinoatrial node (SAN). Surprisingly, little is known regarding the biochemical composition or mechanical properties of the SAN, and how the unique structural characteristics present in this region of the heart influence CPC function remains poorly understood. Here, we have identified that SAN development involves the construction of a "soft" macromolecular ECM that specifically encapsulates CPCs. In addition, we demonstrate that subjecting embryonic CPCs to substrate stiffnesses higher than those measured in vivo results in loss of coherent electrical oscillation and dysregulation of the HCN4 and NCX1 ion channels required for CPC automaticity. Collectively, these data indicate that local mechanics play a critical role in maintaining the embryonic CPC function while also quantitatively defining the range of material properties that are optimal for embryonic CPC maturation.

Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish

Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.

Conversion of Unmodified Stem Cells to Pacemaker Cells by Overexpression of Key Developmental Genes

Arrhythmias of the heart are currently treated by implanting electronic pacemakers and defibrillators. Unmodified adipose tissue-derived stem cells (ASCs) have the potential to differentiate into all three germ layers but have not yet been tested for the generation of pacemaker and Purkinje cells. We investigated if-based on overexpression of dominant conduction cell-specific genes in ASCs-biological pacemaker cells could be induced. Here we show that by overexpression of certain genes that are active during the natural development of the conduction system, the differentiation of ASCs to pacemaker and Purkinje-like cells is feasible. Our study revealed that the most effective procedure consisted of short-term upregulation of gene combinations SHOX2-TBX5-HCN2, and to a lesser extent SHOX2-TBX3-HCN2. Single-gene expression protocols were ineffective. Future clinical implantation of such pacemaker and Purkinje cells, derived from unmodified ASCs of the same patient, could open up new horizons for the treatment of arrythmias.

A potential immunotherapeutic and prognostic biomarker for multiple tumors including glioma: SHOX2

BACKGROUND

Short stature homeobox 2 (SHOX2) is significant gene in the development and progression of multiple types of tumors. Nonetheless, the biological role of SHOX2 within pan-cancer datasets has not been investigated. Thus, comprehensive bioinformatics analyses of pan-cancer datasets were conducted to explore how SHOX2 regulates tumorigenesis.

METHODS

A variety of tumor datasets and online analytical tools, including SangerBox, TIMER2, LinkedOmic, GEPIA2 and cBioPortal, were applied to explore SHOX2 expression in various tumors. To ascertain the connections between SHOX2 expression and genetic alterations, SHOX2-related genes and tumor immunity, the pan-cancer datasets were examined. In vitro assays were applied to verify the biological functions of SHOX2 in glioma cells via CCK-8, wound healing, Transwell and colony formation assays.

RESULTS

Analyses found that SHOX2 was overexpressed in multiple cancer types. SHOX2 expression level was significantly correlated with isocitrate dehydrogenase (IDH), 1p/19q, O⁶-methylguanine DNA methyltransferase (MGMT) status and new types of glioma patients. High mRNA expression levels of SHOX2 were associated with a poor prognosis in multiple tumor patients. KEGG enrichment analysis showed that SHOX2-related genes were associated with cell cycle and DNA damage repair. Genetic alterations of SHOX2 were identified in multiple types of cancers, including duplications and deep mutations. Immune analysis showed that SHOX2 was closely correlated with the tumor mutation burden (TMB), microsatellite instability (MSI), neoantigen and neoantigens and immune checkpoint (ICP) in a variety of tumors and could influence the immunotherapy sensitivity of cancers. CCK-8, wound healing, Transwell and colony formation experiments showed that SHOX2 knockdown inhibited glioma cell proliferation, migration, invasion and colony formation abilities.

CONCLUSION

SHOX2 was overexpressed in multiple cancer types in TCGA cohort. SHOX2 knockdown inhibited glioma cell proliferation, migration and colony formation ability. Our study showed that SHOX2 may be an immunotherapeutic and promising prognostic biomarker in certain types of tumors.

Sinus node dysfunction: current understanding and future directions

The sinoatrial node (SAN) is the primary pacemaker of the heart. Normal SAN function is crucial in maintaining proper cardiac rhythm and contraction. Sinus node dysfunction (SND) is due to abnormalities within the SAN, which can affect the heartbeat frequency, regularity, and the propagation of electrical pulses through the cardiac conduction system. As a result, SND often increases the risk of cardiac arrhythmias. SND is most commonly seen as a disease of the elderly given the role of degenerative fibrosis as well as other age-dependent changes in its pathogenesis. Despite the prevalence of SND, current treatment is limited to pacemaker implantation, which is associated with substantial medical costs and complications. Emerging evidence has identified various genetic abnormalities that can cause SND, shedding light on the molecular underpinnings of SND. Identification of these molecular mechanisms and pathways implicated in the pathogenesis of SND is hoped to identify novel therapeutic targets for the development of more effective therapies for this disease. In this review article, we examine the anatomy of the SAN and the pathophysiology and epidemiology of SND. We then discuss in detail the most common genetic mutations correlated with SND and provide our perspectives on future research and therapeutic opportunities in this field.

A transcription factor atlas of directed differentiation

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

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.

High cardiomyocyte diversity in human early prenatal heart development

Cardiomyocytes play key roles during cardiogenesis, but have poorly understood features, especially in prenatal stages. Here, we characterized human prenatal cardiomyocytes, 6.5-7 weeks post-conception, by integrating single-cell RNA sequencing, spatial transcriptomics, and ligand-receptor interaction information. Using a computational workflow developed to dissect cell type heterogeneity, localize cell types, and explore their molecular interactions, we identified eight types of developing cardiomyocyte, more than double compared to the ones identified in the Human Developmental Cell Atlas. These have high variability in cell cycle activity, mitochondrial content, and connexin gene expression, and are differentially distributed in the ventricles, including outflow tract, and atria, including sinoatrial node. Moreover, cardiomyocyte ligand-receptor crosstalk is mainly with non-cardiomyocyte cell types, encompassing cardiogenesis-related pathways. Thus, early prenatal human cardiomyocytes are highly heterogeneous and develop unique location-dependent properties, with complex ligand-receptor crosstalk. Further elucidation of their developmental dynamics may give rise to new therapies.

Molecular and Functional Remodeling of Superior and Inferior SAN in a Rat Model of HCM

BACKGROUND

Recently, our laboratory presented functional and molecular evidence for the presence of 2 competing sinoatrial node (SAN) pacemakers in healthy human and rat hearts. Anatomically localized near the superior vena cava and inferior vena cava, the superior and inferior SANs (sSAN and iSAN, respectively) preferentially control fast and slow normal heart rates. However, only 1 dominant pacemaker, primarily the sSAN, was functional in the failing rat heart with hypertrophic cardiomyopathy.

OBJECTIVES

This study aimed to determine the transcriptional basis of functional silencing of 1 of 2 dominant pacemakers in failing rat hearts.

METHODS

Ascending aortic constriction was performed on 1-week-old male Sprague-Dawley rat pups to induce left ventricular hypertrophy and heart failure. The dominant pacemaker was anatomically mapped in adult (10-12 weeks old) healthy and failing rat hearts using optical mapping in isolated right atrial tissue preparations. RNA sequencing was used to identify regional sSAN/iSAN gene expression differences between healthy and failing rat hearts.

RESULTS

In all failing rat hearts optically mapped in this study (n = 4), only the sSAN pacemaker was functional, while the iSAN was silent. Compared to healthy rat hearts, a total of 3,640 genes were downregulated, and 4,518 genes were upregulated in failing rat hearts. The functional quiescence of the iSAN in these failing rat hearts may be explained by their downregulation of sodium, potassium, and calcium ion channels as well as their downregulation of specific structural genes, including ankyrin, titin, and myosin heavy chain. Moreover, the iSAN showed predominant downregulation of several key transcription factors such as Tbx5, Tbx3, Shox2, and Smad9.

CONCLUSIONS

Pressure-overload-induced heart failure resulted in significant downregulation of critical transcription factors, ion channels, and structural transcripts of the iSAN, which could explain the functional silencing of the iSAN in failing rat hearts.

Molecular and electrophysiological evaluation of human cardiomyocyte subtypes to facilitate generation of composite cardiac models

Paucity of physiologically relevant cardiac models has limited the widespread application of human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes in drug development. Here, we performed comprehensive characterization of hiPSC-derived cardiomyocyte subtypes from 2D and 3D cultures and established a novel 3D model to study impulse initiation and propagation. Directed differentiation approaches were used to generate sinoatrial nodal (SANCM), atrial (ACM) and ventricular cardiomyocytes (VCM). Single cell RNA sequencing established that the protocols yield distinct cell populations in line with expected identities, which was also confirmed by electrophysiological characterization. In 3D EHT cultures of all subtypes, we observed prominent expression of stretch-responsive genes such as NPPA. Response to rate modulating drugs noradrenaline, carbachol and ivabradine were comparable in single cells and EHTs. Differences in the speed of impulse propagation between the subtypes were more pronounced in EHTs compared with 2D monolayers owing to a progressive increase in conduction velocities in atrial and ventricular cardiomyocytes, in line with a more mature phenotype. In a novel binary EHT model of pacemaker-atrial interface, the SANCM end of the tissue consistently paced the EHTs under baseline conditions, which was inhibited by ivabradine. Taken together, our data provide comprehensive insights into molecular and electrophysiological properties of hiPSC-derived cardiomyocyte subtypes, facilitating the creation of next generation composite cardiac models for drug discovery, disease modeling and cell-based regenerative therapies.

A single cell transcriptional roadmap of human pacemaker cell differentiation

Each heartbeat is triggered by the sinoatrial node (SAN), the primary pacemaker of the heart. Studies in animal models have revealed that pacemaker cells share a common progenitor with the (pro)epicardium, and that the pacemaker cardiomyocytes further diversify into ‘transitional’, ‘tail’, and ‘head’ subtypes. However, the underlying molecular mechanisms, especially of human pacemaker cell development, are poorly understood. Here, we performed single cell RNA sequencing (scRNA-seq) and trajectory inference on human induced pluripotent stem cells (hiPSCs) differentiating to SAN-like cardiomyocytes (SANCMs) to construct a roadmap of transcriptional changes and lineage decisions. In differentiated SANCM, we identified distinct clusters that closely resemble different subpopulations of the in vivo SAN. Moreover, the presence of a side population of proepicardial cells suggested their shared ontogeny with SANCM, as also reported in vivo. Our results demonstrate that the divergence of SANCM and proepicardial lineages is determined by WNT signaling. Furthermore, we uncovered roles for TGFβ and WNT signaling in the branching of transitional and head SANCM subtypes, respectively. These findings provide new insights into the molecular processes involved in human pacemaker cell differentiation, opening new avenues for complex disease modeling in vitro and inform approaches for cell therapy-based regeneration of the SAN.

NRF-2/HO-1 Pathway-Mediated SHOX2 Activation Is a Key Switch for Heart Rate Acceleration by Yixin-Fumai Granules

Population aging has led to increased sick sinus syndrome (SSS) incidence; however, no effective and safe medical therapy has been reported thus far. Yixin-Fumai granules (YXFMs), a Chinese medicine granule designed for bradyarrhythmia treatment, can effectively increase SSS patients' heart rate. Senescence-induced sinoatrial node (SAN) degeneration is an important part of SSS pathogenesis, and older people often show high levels of oxidative stress; reactive oxygen species (ROS) accumulation in the SAN causes abnormal SAN pacing or conduction functions. The current study observed the protective effects of YXFMs on senescent SAN and explored the relationship between the NRF-2/HO-1 pathway, SHOX2, and T-type calcium channels. We selected naturally senescent C57BL/6 mice with bradycardia to simulate SSS; electrocardiography, Masson's trichrome staining, and DHE staining were used to assess SAN function and tissue damage. Immunofluorescence staining and Western blotting were used to assay related proteins. In vitro, we treated human-induced pluripotent stem cell-derived atrial myocytes (hiPSC-AMs) and mouse atrial myocyte-derived cell line HL-1 with D-galactose to simulate senescent SAN-pacemaker cells. CardioExcyte96 was used to evaluate the pulsatile function of the hiPSC-AMs, and the mechanism was verified by DCFH-DA, immunofluorescence staining, RT-qPCR, and Western blotting. The results demonstrated that YXFMs effectively inhibited senescence-induced SAN hypofunction, and this effect possibly originated from scavenging of ROS and promotion of NRF-2, SHOX2, and T-type calcium channel expression. In vitro experiment results indicated that ML385, si-SHOX2, LDN193189, and Mibefradil reversed YXFMs' effects. Moreover, we, for the first time, found that ROS accumulation may hinder SHOX2 expression; YXFMs can activate SHOX2 through the NRF-2/HO-1 pathway-mediated ROS scavenging and then regulate CACNA1G through the SHOX2/BMP4/GATA4/NKX2-5 axis, improve T-type calcium channel function, and ameliorate the SAN dysfunction. Finally, through network pharmacology and molecular docking, we screened for the most stable YXFMs compound that docks to NRF-2, laying the foundation for future studies.

devCellPy is a machine learning-enabled pipeline for automated annotation of complex multilayered single-cell transcriptomic data

A major informatic challenge in single cell RNA-sequencing analysis is the precise annotation of datasets where cells exhibit complex multilayered identities or transitory states. Here, we present devCellPy a highly accurate and precise machine learning-enabled tool that enables automated prediction of cell types across complex annotation hierarchies. To demonstrate the power of devCellPy, we construct a murine cardiac developmental atlas from published datasets encompassing 104,199 cells from E6.5-E16.5 and train devCellPy to generate a cardiac prediction algorithm. Using this algorithm, we observe a high prediction accuracy (>90%) across multiple layers of annotation and across de novo murine developmental data. Furthermore, we conduct a cross-species prediction of cardiomyocyte subtypes from in vitro-derived human induced pluripotent stem cells and unexpectedly uncover a predominance of left ventricular (LV) identity that we confirmed by an LV-specific TBX5 lineage tracing system. Together, our results show devCellPy to be a useful tool for automated cell prediction across complex cellular hierarchies, species, and experimental systems.

Recombinant Adeno-Associated Viral Vector-Mediated Gene Transfer of hTBX18 Generates Pacemaker Cells from Ventricular Cardiomyocytes

Sinoatrial node dysfunction can manifest as bradycardia, leading to symptoms of syncope and sudden cardiac death. Electronic pacemakers are the current standard of care but are limited due to a lack of biological chronotropic control, cost of revision surgeries, and risk of lead- and device-related complications. We therefore aimed to develop a biological alternative to electronic devices by using a clinically relevant gene therapy vector to demonstrate conversion of cardiomyocytes into sinoatrial node-like cells in an in vitro context. Neonatal rat ventricular myocytes were transduced with recombinant adeno-associated virus vector 6 encoding either hTBX18 or green fluorescent protein and maintained for 3 weeks. At the endpoint, qPCR, Western blot analysis and immunocytochemistry were used to assess for reprogramming into pacemaker cells. Cell morphology and Arclight action potentials were imaged via confocal microscopy. Compared to GFP, hTBX18-transduced cells showed that hTBX18, HCN4 and Cx45 were upregulated. Cx43 was significantly downregulated, while sarcomeric α-actinin remained unchanged. Cardiomyocytes transduced with hTBX18 acquired the tapering morphology of native pacemaker cells, as compared to the block-like, striated appearance of ventricular cardiomyocytes. Analysis of the action potentials showed phase 4 depolarization and a significant decrease in the APD50 of the hTBX18-transduced cells. We have demonstrated that rAAV-hTBX18 gene transfer to ventricular myocytes results in morphological, molecular, physiological, and functional changes, recapitulating the pacemaker phenotype in an in vitro setting. The generation of these induced pacemaker-like cells using a clinically relevant vector opens new prospects for biological pacemaker development.

RA signaling pathway combined with Wnt signaling pathway regulates human-induced pluripotent stem cells (hiPSCs) differentiation to sinus node-like cells

BACKGROUND

The source of SAN is debated among researchers. Many studies have shown that RA and Wnt signaling are involved in heart development. In this study, we investigated the role of retinoic acid (RA) and Wnt signaling in the induction of sinus node-like cells.

METHODS

The experimental samples were divided into four groups: control group (CHIR = 0), CHIR = 3, RA + CHIR = 0 andRA + CHIR = 3. After 20 days of differentiation, Western blot, RT-qPCR, immunofluorescence and flow cytometry were performed to identify sinus node-like cells. Finally, whole-cell patch clamp technique was used to record pacing funny current and action potential (AP) in four groups.

RESULTS

The best intervention method used in our experiment was RA = 0.25 µmol/L D5-D9 + CHIR = 3 µmol/L D5-D7. Results showed that CHIR can increase the expression of ISL-1 and TBX3, while RA mainly elevated Shox2. Immunofluorescence assay and flow cytometry further illustrated that combining RA with CHIR can induce sinus node-like cells (CTNT+Shox2+Nkx2.5-). Moreover, CHIR might reduce the frequency of cell beats, but in conjunction with RA could partly compensate for this side effect. Whole cell patch clamps were able to record funny current and the typical sinus node AP in the experimental group, which did not appear in the control group.

CONCLUSIONS

Combining RA with Wnt signaling within a specific period can induce sinus node-like cells.

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

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.

Comparative Study of Transcriptome in the Hearts Isolated from Mice, Rats, and Humans

The heart is a significant organ in mammalian life, and the heartbeat mechanism has been an essential focus of science. However, few studies have focused on species differences. Accordingly, challenges remain in studying genes that have universal functions across species and genes that determine species differences. Here, we analyzed transcriptome data in mouse, rat, and human atria, ventricles, and sinoatrial nodes (SA) obtained from different platforms and compared them by calculating specificity measure (SPM) values in consideration of species differences. Among the three heart regions, the species differences in SA were the greatest, and we searched for genes that determined the essential characteristics of SA, which was SHOX2 in our criteria. The SPM value of SHOX2 was prominently high across species. Similarly, by calculating SPM values, we identified 3 atrial-specific, 11 ventricular-specific, and 17 SA-specific markers. Ontology analysis identified 70 cardiac region- and species-specific ontologies. These results suggest that reanalyzing existing data by calculating SPM values may identify novel tissue-specific genes and species-dependent gene expression. This study identified the importance of SHOX2 as an SA-specific transcription factor, a novel cardiac regional marker, and species-dependent ontologies.

Preferential Expression of Ca2+-Stimulable Adenylyl Cyclase III in the Supraventricular Area, Including Arrhythmogenic Pulmonary Vein of the Rat Heart

Ectopic excitability in pulmonary veins (PVs) is the major cause of atrial fibrillation. We previously reported that the inositol trisphosphate receptor in rat PV cardiomyocytes cooperates with the Na+-Ca2+ exchanger to provoke ectopic automaticity in response to norepinephrine. Here, we focused on adenylyl cyclase (AC) as another effector of norepinephrine stimulation. RT-PCR, immunohistochemistry, and Western blotting revealed that the abundant expression of Ca2+-stimulable AC3 was restricted to the supraventricular area, including the PVs. All the other AC isotypes hardly displayed any region-specific expressions. Immunostaining of isolated cardiomyocytes showed an enriched expression of AC3 along the t-tubules in PV myocytes. The cAMP-dependent response of L-type Ca2+ currents in the PV and LA cells is strengthened by the 0.1 mM intracellular Ca2+ condition, unlike in the ventricular cells. The norepinephrine-induced automaticity of PV cardiomyocytes was reversibly suppressed by 100 µM SQ22536, an adenine-like AC inhibitor. These findings suggest that the specific expression of AC3 along t-tubules may contribute to arrhythmogenic automaticity in rat PV cardiomyocytes.

A dual SHOX2:GFP; MYH6:mCherry knockin hESC reporter line for derivation of human SAN-like cells

The sinoatrial node (SAN) is the primary pacemaker of the heart. The human SAN is poorly understood due to limited primary tissue access and limitations in robust in vitro derivation methods. We developed a dual SHOX2:GFP; MYH6:mCherry knockin human embryonic stem cell (hESC) reporter line, which allows the identification and purification of SAN-like cells. Using this line, we performed several rounds of chemical screens and developed an efficient strategy to generate and purify hESC-derived SAN-like cells (hESC-SAN). The derived hESC-SAN cells display molecular and electrophysiological characteristics of bona fide nodal cells, which allowed exploration of their transcriptional profile at single-cell level. In sum, our dual reporter system facilitated an effective strategy for deriving human SAN-like cells, which can potentially be used for future disease modeling and drug discovery.

GOing Forward With the Cardiac Conduction System Using Gene Ontology

The cardiac conduction system (CCS) comprises critical components responsible for the initiation, propagation, and coordination of the action potential. Aberrant CCS development can cause conduction abnormalities, including sick sinus syndrome, accessory pathways, and atrioventricular and bundle branch blocks. Gene Ontology (GO; http://geneontology.org/) is an invaluable global bioinformatics resource which provides structured, computable knowledge describing the functions of gene products. Many gene products are known to be involved in CCS development; however, this information is not comprehensively captured by GO. To address the needs of the heart development research community, this study aimed to describe the specific roles of proteins reported in the literature to be involved with CCS development and/or function. 14 proteins were prioritized for GO annotation which led to the curation of 15 peer-reviewed primary experimental articles using carefully selected GO terms. 152 descriptive GO annotations, including those describing sinoatrial node and atrioventricular node development were created and submitted to the GO Consortium database. A functional enrichment analysis of 35 key CCS development proteins confirmed that this work has improved the in-silico interpretation of this CCS dataset. This work may improve future investigations of the CCS with application of high-throughput methods such as genome-wide association studies analysis, proteomics, and transcriptomics.

Human-induced pluripotent stem cells in cardiovascular research: current approaches in cardiac differentiation, maturation strategies, and scalable production

Manifestations of cardiovascular diseases (CVDs) in a patient or a population differ based on inherent biological makeup, lifestyle, and exposure to environmental risk factors. These variables mean that therapeutic interventions may not provide the same benefit to every patient. In the context of CVDs, human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer an opportunity to model CVDs in a patient-specific manner. From a pharmacological perspective, iPSC-CM models can serve as go/no-go tests to evaluate drug safety. To develop personalized therapies for early diagnosis and treatment, human-relevant disease models are essential. Hence, to implement and leverage the utility of iPSC-CMs for large-scale treatment or drug discovery, it is critical to (i) carefully evaluate the relevant limitations of iPSC-CM differentiations, (ii) establish quality standards for defining the state of cell maturity, and (iii) employ techniques that allow scalability and throughput with minimal batch-to-batch variability. In this review, we briefly describe progress made with iPSC-CMs in disease modelling and pharmacological testing, as well as current iPSC-CM maturation techniques. Finally, we discuss current platforms for large-scale manufacturing of iPSC-CMs that will enable high-throughput drug screening applications.

Modeling Human Cardiac Arrhythmias: Insights from Zebrafish

Cardiac arrhythmia, or irregular heart rhythm, is associated with morbidity and mortality and is described as one of the most important future public health challenges. Therefore, developing new models of cardiac arrhythmia is critical for understanding disease mechanisms, determining genetic underpinnings, and developing new therapeutic strategies. In the last few decades, the zebrafish has emerged as an attractive model to reproduce in vivo human cardiac pathologies, including arrhythmias. Here, we highlight the contribution of zebrafish to the field and discuss the available cardiac arrhythmia models. Further, we outline techniques to assess potential heart rhythm defects in larval and adult zebrafish. As genetic tools in zebrafish continue to bloom, this model will be crucial for functional genomics studies and to develop personalized anti-arrhythmic therapies.

PRRX1 Loss-of-Function Mutations Underlying Familial Atrial Fibrillation

Background

Atrial fibrillation (AF) is the most common form of clinical cardiac dysrhythmia responsible for thromboembolic cerebral stroke, congestive heart failure, and death. Aggregating evidence highlights the strong genetic basis of AF. Nevertheless, AF is of pronounced genetic heterogeneity, and in an overwhelming majority of patients, the genetic determinants underpinning AF remain elusive.

Methods and Results

By genome‐wide screening with polymorphic microsatellite markers and linkage analysis in a 4‐generation Chinese family affected with autosomal‐dominant AF, a novel locus for AF was mapped to chromosome 1q24.2–q25.1, a 3.20‐cM (≈4.19 Mbp) interval between markers D1S2851 and D1S218, with the greatest 2‐point logarithm of odds score of 4.8165 for the marker D1S452 at recombination fraction=0.00. Whole‐exome sequencing and bioinformatics analyses showed that within the mapping region, only the mutation in the paired related homeobox 1 (PRRX1) gene, NM_022716.4:c.319C>T;(p.Gln107*), cosegregated with AF in the family. In addition, sequencing analyses of PRRX1 in another cohort of 225 unrelated patients with AF revealed a new mutation, NM_022716.4:c.437G>T; (p.Arg146Ile), in a patient. The 2 mutations were absent in 908 control subjects. Biological analyses in HeLa cells demonstrated that the 2 mutants had significantly diminished transactivation on the target genes ISL1 and SHOX2 and markedly decreased ability to bind the promoters of ISL1 and SHOX2 (2 genes causally linked to AF), although with normal intracellular distribution.

Conclusions

This study first indicates that PRRX1 loss‐of‐function mutations predispose to AF, which provides novel insight into the molecular pathogenesis underpinning AF, implying potential implications for precisive prophylaxis and management of AF.

Inherited and Acquired Rhythm Disturbances in Sick Sinus Syndrome, Brugada Syndrome, and Atrial Fibrillation: Lessons from Preclinical Modeling

Rhythm disturbances are life-threatening cardiovascular diseases, accounting for many deaths annually worldwide. Abnormal electrical activity might arise in a structurally normal heart in response to specific triggers or as a consequence of cardiac tissue alterations, in both cases with catastrophic consequences on heart global functioning. Preclinical modeling by recapitulating human pathophysiology of rhythm disturbances is fundamental to increase the comprehension of these diseases and propose effective strategies for their prevention, diagnosis, and clinical management. In silico, in vivo, and in vitro models found variable application to dissect many congenital and acquired rhythm disturbances. In the copious list of rhythm disturbances, diseases of the conduction system, as sick sinus syndrome, Brugada syndrome, and atrial fibrillation, have found extensive preclinical modeling. In addition, the electrical remodeling as a result of other cardiovascular diseases has also been investigated in models of hypertrophic cardiomyopathy, cardiac fibrosis, as well as arrhythmias induced by other non-cardiac pathologies, stress, and drug cardiotoxicity. This review aims to offer a critical overview on the effective ability of in silico bioinformatic tools, in vivo animal studies, in vitro models to provide insights on human heart rhythm pathophysiology in case of sick sinus syndrome, Brugada syndrome, and atrial fibrillation and advance their safe and successful translation into the cardiology arena.

Spatiotemporal expression of long noncoding RNA Moshe modulates heart cell lineage commitment

The roles of long non-coding RNA (LncRNA) have been highlighted in various development processes including congenital heart defects (CHD). Here, we characterized the molecular function of LncRNA, Moshe (1010001N08ik-203), one of the Gata6 antisense transcripts located upstream of Gata6, which is involved in both heart development and the most common type of congenital heart defect, atrial septal defect (ASD). During mouse embryonic development, Moshe was first detected during the cardiac mesoderm stage (E8.5 to E9.5) where Gata6 is expressed and continues to increase at the atrioventricular septum (E12.5), which is involved in ASD. Functionally, the knock-down of Moshe during cardiogenesis caused significant repression of Nkx2.5 in cardiac progenitor stages and resulted in the increase in major SHF lineage genes, such as cardiac transcriptional factors (Isl1, Hand2, Tbx2), endothelial-specific genes (Cd31, Flk1, Tie1, vWF), a smooth muscle actin (a-Sma) and sinoatrial node-specific genes (Shox2, Tbx18). Chromatin Isolation by RNA Purification showed Moshe activates Nkx2.5 gene expression via direct binding to its promoter region. Of note, Moshe was conserved across species, including human, pig and mouse. Altogether, this study suggests that Moshe is a heart-enriched lncRNA that controls a sophisticated network of cardiogenesis by repressing genes in SHF via Nkx2.5 during cardiac development and may play an important role in ASD.

Implementing Biological Pacemakers: Design Criteria for Successful

Each heartbeat that pumps blood throughout the body is initiated by an electrical impulse generated in the sinoatrial node (SAN). However, a number of disease conditions can hamper the ability of the SAN's pacemaker cells to generate consistent action potentials and maintain an orderly conduction path, leading to arrhythmias. For symptomatic patients, current treatments rely on implantation of an electronic pacing device. However, complications inherent to the indwelling hardware give pause to categorical use of device therapy for a subset of populations, including pediatric patients or those with temporary pacing needs. Cellular-based biological pacemakers, derived in vitro or in situ, could function as a therapeutic alternative to current electronic pacemakers. Understanding how biological pacemakers measure up to the SAN would facilitate defining and demonstrating its advantages over current treatments. In this review, we discuss recent approaches to creating biological pacemakers and delineate design criteria to guide future progress based on insights from basic biology of the SAN. We emphasize the need for long-term efficacy in vivo via maintenance of relevant proteins, source-sink balance, a niche reflective of the native SAN microenvironment, and chronotropic competence. With a focus on such criteria, combined with delivery methods tailored for disease indications, clinical implementation will be attainable.

Inducible cardiomyocyte injury within the atrioventricular conduction system uncovers latent regenerative capacity in mice

The cardiac conduction system (CCS) ensures regular contractile function, and injury to any of its components can cause cardiac dysrhythmia. Although all cardiomyocytes (CMs) originate from common progenitors, the CCS is composed of biologically distinct cell types with unique functional and developmental characteristics. In contrast to ventricular cardiomyocytes, which continue to proliferate after birth, most CCS cells terminally exit the cell cycle during fetal development. Although the CCS should thus provide a poor substrate for postnatal injury repair, its regenerative capacity remains untested. Here, we describe a genetic system for ablating CMs that reside within the atrioventricular conduction system (AVCS). Adult mouse AVCS ablation resulted in regenerative failure characterized by persistent atrioventricular conduction defects and contractile dysfunction. In contrast, AVCS injury in neonatal mice led to recovery in a subset of these mice, thus providing evidence for CCS plasticity. Furthermore, CM proliferation did not appear to completely account for the observed functional recovery, suggesting that mechanisms regulating recovery from dysrhythmia are likely to be distinct from cardiac regeneration associated with ventricular injury. Taken together, we anticipate that our results will motivate further mechanistic studies of CCS plasticity and enable the exploration of rhythm restoration as an alternative therapeutic strategy.

Genomic and physiological analyses of the zebrafish atrioventricular canal reveal molecular building blocks of the secondary pacemaker region

The atrioventricular canal (AVC) is the site where key structures responsible for functional division between heart regions are established, most importantly, the atrioventricular (AV) conduction system and cardiac valves. To elucidate the mechanism underlying AVC development and function, we utilized transgenic zebrafish line sqet31Et expressing EGFP in the AVC to isolate this cell population and profile its transcriptome at 48 and 72 hpf. The zebrafish AVC transcriptome exhibits hallmarks of mammalian AV node, including the expression of genes implicated in its development and those encoding connexins forming low conductance gap junctions. Transcriptome analysis uncovered protein-coding and noncoding transcripts enriched in AVC, which have not been previously associated with this structure, as well as dynamic expression of epithelial-to-mesenchymal transition markers and components of TGF-β, Notch, and Wnt signaling pathways likely reflecting ongoing AVC and valve development. Using transgenic line Tg(myl7:mermaid) encoding voltage-sensitive fluorescent protein, we show that abolishing the pacemaker-containing sinoatrial ring (SAR) through Isl1 loss of function resulted in spontaneous activation in the AVC region, suggesting that it possesses inherent automaticity although insufficient to replace the SAR. The SAR and AVC transcriptomes express partially overlapping species of ion channels and gap junction proteins, reflecting their distinct roles. Besides identifying conserved aspects between zebrafish and mammalian conduction systems, our results established molecular hallmarks of the developing AVC which underlies its role in structural and electrophysiological separation between heart chambers. This data constitutes a valuable resource for studying AVC development and function, and identification of novel candidate genes implicated in these processes.

E2A ablation enhances proportion of nodal-like cardiomyocytes in cardiac-specific differentiation of human embryonic stem cells

BACKGROUND

Human sinoatrial cardiomyocytes are essential building blocks for cell therapies of conduction system disorders. However, current differentiation protocols for deriving nodal cardiomyocytes from human pluripotent stem cells (hPSCs) are very inefficient.

METHODS

By employing the hPSCs to cardiomyocyte (CM) in vitro differentiation system and generating E2A-knockout hESCs using CRISPR/Cas9 gene editing technology, we analyze the functions of E2A in CM differentiation.

FINDINGS

We found that knockout of the transcription factor E2A substantially increased the proportion of nodal-like cells in hESC-derived CMs. The E2A ablated CMs displayed smaller cell size, increased beating rates, weaker contractile force, and other functional characteristics similar to sinoatrial node (SAN) cells. Transcriptomic analyses indicated that ion channel-encoding genes were up-regulated in E2A ablated CMs. E2A directly bounded to the promoters of genes key to SAN development via conserved E-box motif, and promoted their expression. Unexpect enhanced activity of NOTCH pathway after E2A ablation could also facilate to induct ventricle workingtype CMs reprogramming into SAN-like cells.

INTERPRETATION

Our study revealed a new role for E2A during directed cardiac differentiation of hESCs and may provide new clues for enhancing induction efficiency of SAN-like cardiomyocytes from hPSCs in the future.

FUNDING

This work was supported by the NSFC (No.82070391, N.S.; No.81870175 and 81922006, P.L.), the National Key R&D Program of China (2018YFC2000202, N.S.; 2017YFA0103700, P.L.), the Haiju program of National Children's Medical Center EK1125180102, and Innovative research team of high-level local universities in Shanghai and a key laboratory program of the Education Commission of Shanghai Municipality (ZDSYS14005).

SHOX2 cooperates with STAT3 to promote breast cancer metastasis through the transcriptional activation of WASF3

BACKGROUND

Metastasis is most often the root cause of cancer-related death. Human short stature homeobox 2 (SHOX2), a homeodomain transcription factor, is a novel inducer of epithelial-to-mesenchymal transition in breast cancer cells, though its exact role and underlying mechanisms in metastasis are not well understood.

METHODS

TCGA analysis was performed to identify the clinical relevance of SHOX2 in breast cancer. Gene depletion was achieved by short hairpin RNA and small interfering RNA. Molecular regulations and alterations were assessed by Western blotting, immunoprecipitation, immunohistochemistry, qRT-PCR, chromatin immunoprecipitation coupled with qPCR (ChIP-qPCR), and ChIP/re-ChIP. The impact of SHOX2 signaling on tumor growth and metastasis was evaluated in orthotopic breast tumor mice.

RESULTS

The expression level of SHOX2 is strongly associated with poor distant metastasis-free survival in breast cancer patients and inactivation of SHOX2 suppresses breast tumor growth and metastasis in mice. In breast cancer cells, SHOX2 directly activates Wiskott-Aldridge syndrome protein family member 3 (WASF3), a metastasis-promoting gene, at the transcriptional level, leading to a significant increase in metastatic potential. Mechanistically, SHOX2 activates signal transducer and activator of transcription 3 (STAT3) and recruits it to the WASF3 promoter, where STAT3 cooperates with SHOX2 to form a functional immunocomplex to promote WASF3 transcriptional activity in breast cancer cells. WASF3 knockdown abrogates SHOX2-induced metastasis, but not SHOX2-dependent tumorigenesis.

CONCLUSIONS

These findings provide a critical link between the SHOX2-STAT3-WASF3 signaling axis and metastasis and suggest that the targeting of this signaling node may represent a valuable alternative strategy for combating breast cancer metastasis.

Zebrafish Models in Therapeutic Research of Cardiac Conduction Disease

Malfunction in the cardiac conduction system (CCS) due to congenital anomalies or diseases can cause cardiac conduction disease (CCD), which results in disturbances in cardiac rhythm, leading to syncope and even sudden cardiac death. Insights into development of the CCS components, including pacemaker cardiomyocytes (CMs), atrioventricular node (AVN) and the ventricular conduction system (VCS), can shed light on the pathological and molecular mechanisms underlying CCD, provide approaches for generating human pluripotent stem cell (hPSC)-derived CCS cells, and thus improve therapeutic treatment for such a potentially life-threatening disorder of the heart. However, the cellular and molecular mechanisms controlling CCS development remain elusive. The zebrafish has become a valuable vertebrate model to investigate early development of CCS components because of its unique features such as external fertilization, embryonic optical transparency and the ability to survive even with severe cardiovascular defects during development. In this review, we highlight how the zebrafish has been utilized to dissect the cellular and molecular mechanisms of CCS development, and how the evolutionarily conserved developmental mechanisms discovered in zebrafish could be applied to directing the creation of hPSC-derived CCS cells, therefore providing potential therapeutic strategies that may contribute to better treatment for CCD patients.

Transcriptional and Epigenetic Landscape of Cardiac Pacemaker Cells: Insights Into Cellular Specialization in the Sinoatrial Node

Cardiac pacemaker cells differentiate and functionally specialize early in embryonic development through activation of critical gene regulatory networks. In general, cellular specification and differentiation require that combinations of cell type-specific transcriptional regulators activate expression of key effector genes by binding to DNA regulatory elements including enhancers and promoters. However, because genomic DNA is tightly packaged by histones that must be covalently modified in order to render DNA regulatory elements and promoters accessible for transcription, the process of development and differentiation is intimately connected to the epigenetic regulation of chromatin accessibility. Although the difficulty of obtaining sufficient quantities of pure populations of pacemaker cells has limited progress in this field, the advent of low-input genomic technologies has the potential to catalyze a rapid growth of knowledge in this important area. The goal of this review is to outline the key transcriptional networks that control pacemaker cell development, with particular attention to our emerging understanding of how chromatin accessibility is modified and regulated during pacemaker cell differentiation. In addition, we will discuss the relevance of these findings to adult sinus node function, sinus node diseases, and origins of genetic variation in heart rhythm. Lastly, we will outline the current challenges facing this field and promising directions for future investigation.

PDGFRα-Signaling Is Dispensable for the Development of the Sinoatrial Node After Its Fate Commitment

Palate-derived growth factor receptor α (Pdgfrα) signaling has been reported to play important roles in the cardiac development. A previous study utilizing Pdgfrα conventional knockout mice reported hypoplasia of the sinus venous myocardium including the sinoatrial node (SAN) accompanied by increased expression of Nkx2.5. This mouse line embryos die by E11.5 due to embryonic lethality, rendering them difficult to investigate the details. To elucidate the underlying mechanism, in this study, we revisited this observation by generation of specific ablation of Pdgfrα in the SAN by Shox2-Cre at E9.5, using a Shox2-Cre;Pdgfrα ^flox/flox conditional mouse line. Surprisingly, we found that resultant homozygous mutant mice did not exhibit any malformation in SAN morphology as compared to their wild-type littermates. Further analysis revealed the normal cardiac function in adult mutant mice assessed by the record of heart rate and electrocardiogram and unaltered expression of Nkx2.5 in the E13.5 SAN of Pdgfrα conditional knockout mice. Our results unambiguously demonstrate that Pdgfrα is dispensable for SAN development after its fate commitment in mice.

Hedgehog signalling controls sinoatrial node development and atrioventricular cushion formation

Smoothened is a key receptor of the hedgehog pathway, but the roles of Smoothened in cardiac development remain incompletely understood. In this study, we found that the conditional knockout of Smoothened from the mesoderm impaired the development of the venous pole of the heart and resulted in hypoplasia of the atrium/inflow tract (IFT) and a low heart rate. The blockage of Smoothened led to reduced expression of genes critical for sinoatrial node (SAN) development in the IFT. In a cardiac cell culture model, we identified a Gli2–Tbx5–Hcn4 pathway that controls SAN development. In the mutant embryos, the endocardial-to-mesenchymal transition (EndMT) in the atrioventricular cushion failed, and Bmp signalling was downregulated. The addition of Bmp2 rescued the EndMT in mutant explant cultures. Furthermore, we analysed Gli2⁺ scRNAseq and Tbx5^−/− RNAseq data and explored the potential genes downstream of hedgehog signalling in posterior second heart field derivatives. In conclusion, our study reveals that Smoothened-mediated hedgehog signalling controls posterior cardiac progenitor commitment, which suggests that the mutation of Smoothened might be involved in the aetiology of congenital heart diseases related to the cardiac conduction system and heart valves.

Generation of cardiomyocytes by atrioventricular node cells in long-term cultures

Turnover of cardiac pacemaker cells may occur during the lifetime of the body, and we recently raised the hypothesis that specialized cardiac cells have in common the potential to generate cardiomyocytes from fibroblasts. To examine this hypothesis, we analyzed the ability of atrioventricular node cells (AVNCs) to generate functional cardiomyocytes in long-term culture. AVNCs were isolated from adult guinea pig hearts and cultured for up to three weeks. Under phase-contrast microscopic observation over time, it was found that within a week, a number of fibroblasts gathered around the AVNCs and formed cell clusters, and thereafter the cell clusters started to beat spontaneously. The nascent cell clusters expanded their area gradually by three weeks in culture and expressed specific cardiac genes and proteins. Maturation of newly formed cardiomyocytes seems to be slow in cultures of AVNCs compared with those of sinoatrial node cells. Stimulation of muscarinic receptors with acetylcholine induced a beating rate decrease which was blocked by atropine, and activation of adenylate cyclase activity with forskolin increased the beat rate, while stimulation of beta adrenoceptors by isoproterenol had no effect. These results indicate that AVNCs form a cluster of cells with properties of functional cardiomyocytes and provide evidence to support the hypothesis.

•

Atrioventricular node cells formed spontaneously beating cell clusters in culture.

•

The cell clusters expressed specific cardiac genes and proteins.

•

The properties of the cell clusters depended on their spontaneous beating rate.

•

Specialized cardiac cells may in common have the ability to generate cardiomyocytes.

Altered Expression of Candidate Genes in Mayer-Rokitansky-Küster-Hauser Syndrome May Influence Vaginal Keratinocytes Biology: A Focus on Protein Kinase X

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is a rare and complex disease defined by congenital aplasia of the vagina and uterus in 46,XX women, often associated with kidney and urinary tract anomalies. The aetiopathogenesis of MRKH syndrome is still largely unknown. Herein, we investigated the role of selected candidate genes in the aetiopathogenesis of MRKH syndrome, with a focus on PRKX, which encodes for protein kinase X. Through RT-qPCR analyses performed on vaginal dimple samples from patients, and principal component analysis (PCA), we highlighted a phenotype-related expression pattern of PRKX, MUC1, HOXC8 and GREB1L in MRKH patients. By using an in vitro approach, we proved that PRKX ectopic overexpression in a cell model of vaginal keratinocytes promotes cell motility through epithelial-to-mesenchymal transition (EMT) activation, a fundamental process in urogenital tract morphogenesis. Moreover, our findings showed that PRKX upregulation in vaginal keratinocytes is able to affect transcriptional levels of HOX genes, implicated in urinary and genital tract development. Our study identified the dysregulation of PRKX expression as a possible molecular cause for MRKH syndrome. Moreover, we propose the specific role of PRKX in vaginal keratinocyte biology as one of the possible mechanisms underlying this complex disease.