Defining human cardiac transcription factor hierarchies using integrated single-cell heterogeneity analysis

The HH-GLI2-CKS1B network regulates the proliferation-to-maturation transition of cardiomyocytes

Cardiomyocyte (CM) proliferation and maturation are highly linked processes, however, the extent to which these processes are controlled by a single signaling axis is unclear. Here, we show the previously undescribed role of Hedgehog (HH)-GLI2-CKS1B cascade in regulation of the toggle between CM proliferation and maturation. Here we show downregulation of GLI-signaling in adult human CM, adult murine CM, and in late-stage hiPSC-CM leading to their maturation. In early-stage hiPSC-CM, inhibition of HH- or GLI-proteins enhanced CM maturation with increased maturation indices, increased calcium handling, and transcriptome. Mechanistically, we identified CKS1B, as a new effector of GLI2 in CMs. GLI2 binds the CKS1B promoter to regulate its expression. CKS1B overexpression in late-stage hiPSC-CMs led to increased proliferation with loss of maturation in CMs. Next, analysis of datasets of patients with heart disease showed a significant enrichment of GLI2-signaling in patients with ischemic heart failure (HF) or dilated-cardiomyopathy (DCM) disease, indicating operational GLI2-signaling in the stressed heart. Thus, the Hh-GLI2-CKS1B axis regulates the proliferation-maturation transition and provides targets to enhance cardiac tissue engineering and regenerative therapies.

LMNA -Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-derived iPSC Differentiation Support Cell type and Lineage-specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA -Related Dilated Cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C ( LMNA ) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. Molecular mechanisms of disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA -Related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (4 Patient and 8 Control) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for Cardiac Progenitors to Cardiomyocytes (CM) and Epicardium-Derived Cells (EPDC). Data integration and comparative analyses of Patient and Control cells found cell type and lineage differentially expressed genes (DEG) with enrichment to support pathway dysregulation. Top DEG and enriched pathways included: 10 ZNF genes and RNA polymerase II transcription in Pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CM; LMNA and epigenetic regulation and DDIT4 and mTORC1 signaling in EPDC. Top DEG also included: XIST and other X-linked genes, six imprinted genes: SNRPN , PWAR6 , NDN , PEG10 , MEG3 , MEG8 , and enriched gene sets in metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

OneSC: A computational platform for recapitulating cell state transitions

Computational modelling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a wet lab. Recent advancements in single-cell RNA sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico 'synthetic' cells that faithfully mimic the temporal trajectories. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

Engineered Gold and Silica Nanoparticle-Incorporated Hydrogel Scaffolds for Human Stem Cell-Derived Cardiac Tissue Engineering

Electrically conductive biomaterials and nanomaterials have demonstrated great potential in the development of functional and mature cardiac tissues. In particular, gold nanomaterials have emerged as promising candidates due to their biocompatibility and ease of fabrication for cardiac tissue engineering utilizing rat- or stem cell-derived cardiomyocytes (CMs). However, despite significant advancements, it is still not clear whether the enhancement in cardiac tissue function is primarily due to the electroconductivity features of gold nanoparticles or the structural changes of the scaffold resulting from the addition of these nanoparticles. To address this question, we developed nanoengineered hydrogel scaffolds comprising gelatin methacrylate (GelMA) embedded with either electrically conductive gold nanorods (GNRs) or nonconductive silica nanoparticles (SNPs). This enabled us to simultaneously assess the roles of electrically conductive and nonconductive nanomaterials in the functionality and fate of human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). Our studies revealed that both GNR- and SNP-incorporated hydrogel scaffolds exhibited excellent biocompatibility and similar cardiac cell attachment. Although the expression of sarcomere alpha-actinin did not significantly differ among the conditions, a more organized sarcomere structure was observed within the GNR-embedded hydrogels compared to the nonconductive nanoengineered scaffolds. Furthermore, electrical coupling was notably improved in GNR-embedded scaffolds, as evidenced by the synchronous calcium flux and enhanced calcium transient intensity. While we did not observe a significant difference in the gene expression profile of human cardiac tissues formed on the conductive GNR- and nonconductive SNP-incorporated hydrogels, we noticed marginal improvements in the expression of some calcium and structural genes in the nanomaterial-embedded hydrogel groups as compared to the control condition. Given that the cardiac tissues formed atop the nonconductive SNP-based scaffolds (used as the control for conductivity) also displayed similar levels of gene expression as compared to the conductive hydrogels, it suggests that the electrical conductivity of nanomaterials (i.e., GNRs) may not be the sole factor influencing the function and fate of hiPSC-derived cardiac tissues when cells are cultured atop the scaffolds. Overall, our findings provide additional insights into the role of electrically conductive gold nanoparticles in regulating the functionalities of hiPSC-CMs.

Tachycardia-induced metabolic rewiring as a driver of contractile dysfunction

Prolonged tachycardia-a risk factor for cardiovascular morbidity and mortality-can induce cardiomyopathy in the absence of structural disease in the heart. Here, by leveraging human patient data, a canine model of tachycardia and engineered heart tissue generated from human induced pluripotent stem cells, we show that metabolic rewiring during tachycardia drives contractile dysfunction by promoting tissue hypoxia, elevated glucose utilization and the suppression of oxidative phosphorylation. Mechanistically, a metabolic shift towards anaerobic glycolysis disrupts the redox balance of nicotinamide adenine dinucleotide (NAD), resulting in increased global protein acetylation (and in particular the acetylation of sarcoplasmic/endoplasmic reticulum Ca2+-ATPase), a molecular signature of heart failure. Restoration of NAD redox by NAD+ supplementation reduced sarcoplasmic/endoplasmic reticulum Ca2+-ATPase acetylation and accelerated the functional recovery of the engineered heart tissue after tachycardia. Understanding how metabolic rewiring drives tachycardia-induced cardiomyopathy opens up opportunities for therapeutic intervention.

CD151 expression marks atrial- and ventricular- differentiation from human induced pluripotent stem cells

Current differentiation protocols for human induced pluripotent stem cells (hiPSCs) produce heterogeneous cardiomyocytes (CMs). Although chamber-specific CM selection using cell surface antigens enhances biomedical applications, a cell surface marker that accurately distinguishes between hiPSC-derived atrial CMs (ACMs) and ventricular CMs (VCMs) has not yet been identified. We have developed an approach for obtaining functional hiPSC-ACMs and -VCMs based on CD151 expression. For ACM differentiation, we found that ACMs are enriched in the CD151low population and that CD151 expression is correlated with the expression of Notch4 and its ligands. Furthermore, Notch signaling inhibition followed by selecting the CD151low population during atrial differentiation leads to the highly efficient generation of ACMs as evidenced by gene expression and electrophysiology. In contrast, for VCM differentiation, VCMs exhibiting a ventricular-related gene signature and uniform action potentials are enriched in the CD151high population. Our findings enable the production of high-quality ACMs and VCMs appropriate for hiPSC-derived chamber-specific disease models and other applications.

Unsupervised analysis of whole transcriptome data from human pluripotent stem cells cardiac differentiation

The main objective of the present work was to highlight differences and similarities in gene expression patterns between different pluripotent stem cell cardiac differentiation protocols, using a workflow based on unsupervised machine learning algorithms to analyse the transcriptome of cells cultured as a 2D monolayer or as 3D aggregates. This unsupervised approach effectively allowed to portray the transcriptomic changes that occurred throughout the differentiation processes, with a visual representation of the entire transcriptome. The results allowed to corroborate previously reported data and also to unveil new gene expression patterns. In particular, it was possible to identify a correlation between low cardiomyocyte differentiation efficiencies and the early expression of a set of non-mesodermal genes, which can be further explored as predictive markers of differentiation efficiency. The workflow here developed can also be applied to analyse other stem cell differentiation transcriptomic datasets, envisaging future clinical implementation of cellular therapies.

Modeling drug-induced mitochondrial toxicity with human primary cardiomyocytes

Mitochondrial toxicity induced by therapeutic drugs is a major contributor for cardiotoxicity, posing a serious threat to pharmaceutical industries and patients' lives. However, mitochondrial toxicity testing is not incorporated into routine cardiac safety screening procedures. To accurately model native human cardiomyocytes, we comprehensively evaluated mitochondrial responses of adult human primary cardiomyocytes (hPCMs) to a nucleoside analog, remdesivir (RDV). Comparison of their response to human pluripotent stem cell-derived cardiomyocytes revealed that the latter utilized a mitophagy-based mitochondrial recovery response that was absent in hPCMs. Accordingly, action potential duration was elongated in hPCMs, reflecting clinical incidences of RDV-induced QT prolongation. In a screen for mitochondrial protectants, we identified mitochondrial ROS as a primary mediator of RDV-induced cardiotoxicity. Our study demonstrates the utility of hPCMs in the detection of clinically relevant cardiac toxicities, and offers a framework for hPCM-based high-throughput screening of cardioprotective agents.

Single-cell RNA sequencing reveals the transcriptional heterogeneity of Tbx18-positive cardiac cells during heart development

The T-box family transcription factor 18 (Tbx18) has been found to play a critical role in regulating the development of the mammalian heart during the primary stages of embryonic development while the cellular heterogeneity and landscape of Tbx18-positive (Tbx18+) cardiac cells remain incompletely characterized. Here, we analyzed prior published single-cell RNA sequencing (scRNA-seq) mouse heart data to explore the heterogeneity of Tbx18+ cardiac cell subpopulations and provide a comprehensive transcriptional landscape of Tbx18+ cardiac cells during their development. Bioinformatic analysis methods were utilized to identify the heterogeneity between cell groups. Based on the gene expression characteristics, Tbx18+ cardiac cells can be classified into a minimum of two distinct cell populations, namely fibroblast-like cells and cardiomyocytes. In terms of temporal heterogeneity, these cells exhibit three developmental stages, namely the MEM stage, ML_P0 stage, and P stage Tbx18+ cardiac cells. Furthermore, Tbx18+ cardiac cells encompass several cell types, including cardiac progenitor-like cells, cardiomyocytes, and epicardial/stromal cells, as determined by specific transcriptional regulatory networks. The scRNA-seq results revealed the involvement of extracellular matrix (ECM) signals and epicardial epithelial-to-mesenchymal transition (EMT) in the development of Tbx18+ cardiac cells. The utilization of a lineage-tracing model served to validate the crucial function of Tbx18 in the differentiation of cardiac cells. Consequently, these findings offer a comprehensive depiction of the cellular heterogeneity within Tbx18+ cardiac cells.

Cardiac Transcription Factors and Regulatory Networks

Cardiac development is a fine-tuned process governed by complex transcriptional networks, in which transcription factors (TFs) interact with other regulatory layers. In this chapter, we introduce the core cardiac TFs including Gata, Hand, Nkx2, Mef2, Srf, and Tbx. These factors regulate each other's expression and can also act in a combinatorial manner on their downstream targets. Their disruption leads to various cardiac phenotypes in mice, and mutations in humans have been associated with congenital heart defects. In the second part of the chapter, we discuss different levels of regulation including cis-regulatory elements, chromatin structure, and microRNAs, which can interact with transcription factors, modulate their function, or are downstream targets. Finally, examples of disturbances of the cardiac regulatory network leading to congenital heart diseases in human are provided.

JAK2 as a surface marker for enrichment of human pluripotent stem cells-derived ventricular cardiomyocytes

BACKGROUND

Human pluripotent stem cell (hPSC)-derived cardiomyocytes (CMs) hold great promise for cardiac disease modelling, drug discovery and regenerative medicine. Despite the advancement in various differentiation protocols, the heterogeneity of the generated population composed of diverse cardiac subtypes poses a significant challenge to their practical applications. Mixed populations of cardiac subtypes can compromise disease modelling and drug discovery, while transplanting them may lead to undesired arrhythmias as they may not integrate and synchronize with the host tissue's contractility. It is therefore crucial to identify cell surface markers that could enable high purity of ventricular CMs for subsequent applications.

METHODS

By exploiting the fact that immature CMs expressing myosin light chain 2A (MLC2A) will gradually express myosin light chain 2 V (MLC2V) protein as they mature towards ventricular fate, we isolated signal regulatory protein alpha (SIRPA)-positive CMs expressing intracellular MLC2A or MLC2V using MARIS (method for analysing RNA following intracellular sorting). Subsequently, RNA sequencing analysis was performed to examine the gene expression profile of MLC2A + and MLC2V + sorted CMs. We identified genes that were significantly up-regulated in MLC2V + samples to be potential surface marker candidates for ventricular specification. To validate these surface markers, we performed immunostaining and western blot analysis to measure MLC2A and MLC2V protein expressions in SIRPA + CMs that were either positive or negative for the putative surface markers, JAK2 (Janus kinase 2) or CD200. We then characterized the electrophysiological properties of surface marker-sorted CMs, using fluo-4 AM, a green-fluorescent calcium indicator, to measure the cellular calcium transient at the single cell level. For functional validation, we investigated the response of the surface marker-sorted CMs to vernakalant, an atrial-selective anti-arrhythmic agent.

RESULTS

In this study, while JAK2 and CD200 were identified as potential surface markers for the purification of ventricular-like CMs, the SIRPA+/JAK2+ population showed a higher percentage of MLC2V-expressing cells (~ 90%) compared to SIRPA+/CD200+ population (~ 75%). SIRPA+/JAK2+ sorted CMs exhibited ventricular-like electrophysiological properties, including slower beating rate, slower calcium depolarization and longer calcium repolarization duration. Importantly, vernakalant had limited to no significant effect on the calcium repolarization duration of SIRPA+/JAK2+ population, indicating their enrichment for ventricular-like CMs.

CONCLUSION

Our study lays the groundwork for the identification of cardiac subtype surface markers that allow purification of cardiomyocyte sub-populations. Our findings suggest that JAK2 can be employed as a cell surface marker for enrichment of hPSC-derived ventricular-like CMs.

SARS-CoV-2 spike receptor-binding domain is internalized and promotes protein ISGylation in human induced pluripotent stem cell-derived cardiomyocytes

Although an increased risk of myocarditis has been observed after vaccination with mRNA encoding severe acute respiratory syndrome coronavirus 2 spike protein, its underlying mechanism has not been elucidated. This study investigated the direct effects of spike receptor-binding domain (S-RBD) on human cardiomyocytes differentiated from induced pluripotent stem cells (iPSC-CMs). Immunostaining experiments using ACE2 wild-type (WT) and knockout (KO) iPSC-CMs treated with purified S-RBD demonstrated that S-RBD was bound to ACE2 and internalized into the subcellular space in the iPSC-CMs, depending on ACE2. Immunostaining combined with live cell imaging using a recombinant S-RBD fused to the superfolder GFP (S-RBD-sfGFP) demonstrated that S-RBD was bound to the cell membrane, co-localized with RAB5A, and then delivered from the endosomes to the lysosomes in iPSC-CMs. Quantitative PCR array analysis followed by single cell RNA sequence analysis clarified that S-RBD-sfGFP treatment significantly upregulated the NF-kβ pathway-related gene (CXCL1) in the differentiated non-cardiomyocytes, while upregulated interferon (IFN)-responsive genes (IFI6, ISG15, and IFITM3) in the matured cardiomyocytes. S-RBD-sfGFP treatment promoted protein ISGylation, an ISG15-mediated post-translational modification in ACE2-WT-iPSC-CMs, which was suppressed in ACE2-KO-iPSC-CMs. Our experimental study demonstrates that S-RBD is internalized through the endolysosomal pathway, which upregulates IFN-responsive genes and promotes ISGylation in the iPSC-CMs.

Targeted activation of human ether-à-go-go-related gene channels rescues electrical instability induced by the R56Q+/- long QT syndrome variant

AIMS

Long QT syndrome type 2 (LQTS2) is associated with inherited variants in the cardiac human ether-à-go-go-related gene (hERG) K+ channel. However, the pathogenicity of hERG channel gene variants is often uncertain. Using CRISPR-Cas9 gene-edited hiPSC-derived cardiomyocytes (hiPSC-CMs), we investigated the pathogenic mechanism underlying the LQTS-associated hERG R56Q variant and its phenotypic rescue by using the Type 1 hERG activator, RPR260243.

METHODS AND RESULTS

The above approaches enable characterization of the unclear causative mechanism of arrhythmia in the R56Q variant (an N-terminal PAS domain mutation that primarily accelerates channel deactivation) and translational investigation of the potential for targeted pharmacologic manipulation of hERG deactivation. Using perforated patch clamp electrophysiology of single hiPSC-CMs, programmed electrical stimulation showed that the hERG R56Q variant does not significantly alter the mean action potential duration (APD90). However, the R56Q variant increases the beat-to-beat variability in APD90 during pacing at constant cycle lengths, enhances the variance of APD90 during rate transitions, and increases the incidence of 2:1 block. During paired S1-S2 stimulations measuring electrical restitution properties, the R56Q variant was also found to increase the variability in rise time and duration of the response to premature stimulations. Application of the hERG channel activator, RPR260243, reduces the APD variance in hERG R56Q hiPSC-CMs, reduces the variability in responses to premature stimulations, and increases the post-repolarization refractoriness.

CONCLUSION

Based on our findings, we propose that the hERG R56Q variant leads to heterogeneous APD dynamics, which could result in spatial dispersion of repolarization and increased risk for re-entry without significantly affecting the average APD90. Furthermore, our data highlight the antiarrhythmic potential of targeted slowing of hERG deactivation gating, which we demonstrate increases protection against premature action potentials and reduces electrical heterogeneity in hiPSC-CMs.

PITX2 gain-of-function mutation associated with atrial fibrillation alters mitochondrial activity in human iPSC atrial-like cardiomyocytes

Atrial fibrillation (AF) is the most common cardiac arrhythmia worldwide; however, the underlying causes of AF initiation are still poorly understood, particularly because currently available models do not allow in distinguishing the initial causes from maladaptive remodeling that induces and perpetuates AF. Lately, the genetic background has been proven to be important in the AF onset. iPSC-derived cardiomyocytes, being patient- and mutation-specific, may help solve this diatribe by showing the initial cell-autonomous changes underlying the development of the disease. Transcription factor paired-like homeodomain 2 (PITX2) has been identified as a key regulator of atrial development/differentiation, and the PITX2 genomic locus has the highest association with paroxysmal AF. PITX2 influences mitochondrial activity, and alterations in either its expression or function have been widely associated with AF. In this work, we investigate the activity of mitochondria in iPSC-derived atrial cardiomyocytes (aCMs) obtained from a young patient (24 years old) with paroxysmal AF, carrying a gain-of-function mutation in PITX2 (rs138163892) and from its isogenic control (CTRL) in which the heterozygous point mutation has been reverted to WT. PITX2 aCMs show a higher mitochondrial content, increased mitochondrial activity, and superoxide production under basal conditions when compared to CTRL aCMs. However, increasing mitochondrial workload by FCCP or β-adrenergic stimulation allows us to unmask mitochondrial defects in PITX2 aCMs, which are incapable of responding efficiently to the higher energy demand, determining ATP deficiency.

Distinct mononuclear diploid cardiac subpopulation with minimal cell-cell communications persists in embryonic and adult mammalian heart

A small proportion of mononuclear diploid cardiomyocytes (MNDCMs), with regeneration potential, could persist in adult mammalian heart. However, the heterogeneity of MNDCMs and changes during development remains to be illuminated. To this end, 12 645 cardiac cells were generated from embryonic day 17.5 and postnatal days 2 and 8 mice by single-cell RNA sequencing. Three cardiac developmental paths were identified: two switching to cardiomyocytes (CM) maturation with close CM-fibroblast (FB) communications and one maintaining MNDCM status with least CM-FB communications. Proliferative MNDCMs having interactions with macrophages and non-proliferative MNDCMs (non-pMNDCMs) with minimal cell-cell communications were identified in the third path. The non-pMNDCMs possessed distinct properties: the lowest mitochondrial metabolisms, the highest glycolysis, and high expression of Myl4 and Tnni1. Single-nucleus RNA sequencing and immunohistochemical staining further proved that the Myl4+Tnni1+ MNDCMs persisted in embryonic and adult hearts. These MNDCMs were mapped to the heart by integrating the spatial and single-cell transcriptomic data. In conclusion, a novel non-pMNDCM subpopulation with minimal cell-cell communications was unveiled, highlighting the importance of microenvironment contribution to CM fate during maturation. These findings could improve the understanding of MNDCM heterogeneity and cardiac development, thus providing new clues for approaches to effective cardiac regeneration.

Hey2 enhancer activity defines unipotent progenitors for left ventricular cardiomyocytes in juxta-cardiac field of early mouse embryo

The cardiac crescent is the first structure of the heart and contains progenitor cells of the first heart field, which primarily differentiate into left ventricular cardiomyocytes. The interface between the forming cardiac crescent and extraembryonic tissue is known as the juxta-cardiac field (JCF), and progenitor cells in this heart field contribute to the myocardium of the left ventricle and atrioventricular canal as well as the epicardium. However, it is unclear whether there are progenitor cells that differentiate specifically into left ventricular cardiomyocytes. We have previously demonstrated that an enhancer of the gene encoding the Hey2 bHLH transcriptional repressor is activated in the ventricular myocardium during mouse embryonic development. In this study, we aimed to investigate the characteristics of cardiomyocyte progenitor cells and their cell lineages by analyzing Hey2 enhancer activity at the earliest stages of heart formation. We found that the Hey2 enhancer initiated its activity prior to cardiomyocyte differentiation within the JCF. Hey2 enhancer-active cells were present rostrally to the Tbx5-expressing region at the early phase of cardiac crescent formation and differentiated exclusively into left ventricular cardiomyocytes in a lineage distinct from the Tbx5-positive lineage. By the late phase of cardiac crescent formation, Hey2 enhancer activity became significantly overlapped with Tbx5 expression in cells that contribute to the left ventricular myocardium. Our study reveals that a population of unipotent progenitor cells for left ventricular cardiomyocytes emerge in the JCF, providing further insight into the mode of cell type diversification during early cardiac development.

Activation of AMPK promotes cardiac differentiation by stimulating the autophagy pathway

Autophagy, a critical catabolic process for cell survival against different types of stress, has a role in the differentiation of various cells, such as cardiomyocytes. Adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK) is an energy-sensing protein kinase involved in the regulation of autophagy. In addition to its direct role in regulating autophagy, AMPK can also influence other cellular processes by regulating mitochondrial function, posttranslational acetylation, cardiomyocyte metabolism, mitochondrial autophagy, endoplasmic reticulum stress, and apoptosis. As AMPK is involved in the control of various cellular processes, it can influence the health and survival of cardiomyocytes. This study investigated the effects of an AMPK inducer (Metformin) and an autophagy inhibitor (Hydroxychloroquine) on the differentiation of human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs). The results showed that autophagy was upregulated during cardiac differentiation. Furthermore, AMPK activation increased the expression of CM-specific markers in hPSC-CMs. Additionally, autophagy inhibition impaired cardiomyocyte differentiation by targeting autophagosome-lysosome fusion. These results indicate the significance of autophagy in cardiomyocyte differentiation. In conclusion, AMPK might be a promising target for the regulation of cardiomyocyte generation by in vitro differentiation of pluripotent stem cells.

Recent advances in regulating the proliferation or maturation of human-induced pluripotent stem cell-derived cardiomyocytes

In the last decade, human-induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM)-based cell therapy has drawn broad attention as a potential therapy for treating injured hearts. However, mass production of hiPSC-CMs remains challenging, limiting their translational potential in regenerative medicine. Therefore, multiple strategies including cell cycle regulators, small molecules, co-culture systems, and epigenetic modifiers have been used to improve the proliferation of hiPSC-CMs. On the other hand, the immaturity of these proliferative hiPSC-CMs could lead to lethal arrhythmias due to their limited ability to functionally couple with resident cardiomyocytes. To achieve functional maturity, numerous methods such as prolonged culture, biochemical or biophysical stimulation, in vivo transplantation, and 3D culture approaches have been employed. In this review, we summarize recent approaches used to promote hiPSC-CM proliferation, and thoroughly review recent advances in promoting hiPSC-CM maturation, which will serve as the foundation for large-scale production of mature hiPSC-CMs for future clinical applications.

APLNR marks a cardiac progenitor derived with human induced pluripotent stem cells

Cardiomyocytes can be readily derived from human induced pluripotent stem cell (hiPSC) lines, yet its efficacy varies across different batches of the same and different hiPSC lines. To unravel the inconsistencies of in vitro cardiac differentiation, we utilized single cell transcriptomics on hiPSCs undergoing cardiac differentiation and identified cardiac and extra-cardiac lineages throughout differentiation. We further identified APLNR as a surface marker for in vitro cardiac progenitors and immunomagnetically isolated them. Differentiation of isolated in vitro APLNR+ cardiac progenitors derived from multiple hiPSC lines resulted in predominantly cardiomyocytes accompanied with cardiac mesenchyme. Transcriptomic analysis of differentiating in vitro APLNR+ cardiac progenitors revealed transient expression of cardiac progenitor markers before further commitment into cardiomyocyte and cardiac mesenchyme. Analysis of in vivo human and mouse embryo single cell transcriptomic datasets have identified APLNR expression in early cardiac progenitors of multiple lineages. This platform enables generation of in vitro cardiac progenitors from multiple hiPSC lines without genetic manipulation, which has potential applications in studying cardiac development, disease modelling and cardiac regeneration.

Insights into Endothelin Receptors in Pulmonary Hypertension

Pulmonary hypertension (PH) is a disease which affects the cardiopulmonary system; it is defined as a mean pulmonary artery pressure (mPAP) > 20 mmHg as measured by right heart catheterization at rest, and is caused by complex and diverse mechanisms. In response to stimuli such as hypoxia and ischemia, the expression and synthesis of endothelin (ET) increase, leading to the activation of various signaling pathways downstream of it and producing effects such as the induction of abnormal vascular proliferation during the development of the disease. This paper reviews the regulation of endothelin receptors and their pathways in normal physiological processes and disease processes, and describes the mechanistic roles of ET receptor antagonists that are currently approved and used in clinical studies. Current clinical researches on ET are focused on the development of multi-target combinations and novel delivery methods to improve efficacy and patient compliance while reducing side effects. In this review, future research directions and trends of ET targets are described, including monotherapy and precision medicine.

Systematic discovery of transcription factors that improve hPSC-derived cardiomyocyte maturation via temporal analysis of bioengineered cardiac tissues

Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) have the potential to become powerful tools for disease modeling, drug testing, and transplantation; however, their immaturity limits their applications. Transcription factor (TF) overexpression can improve hPSC-CM maturity, but identifying these TFs has been elusive. Toward this, we establish here an experimental framework for systematic identification of maturation enhancing factors. Specifically, we performed temporal transcriptome RNAseq analyses of progressively matured hPSC-derived cardiomyocytes across 2D and 3D differentiation systems and further compared these bioengineered tissues to native fetal and adult-derived tissues. These analyses revealed 22 TFs whose expression did not increase in 2D differentiation systems but progressively increased in 3D culture systems and adult mature cell types. Individually overexpressing each of these TFs in immature hPSC-CMs identified five TFs (KLF15, ZBTB20, ESRRA, HOPX, and CAMTA2) as regulators of calcium handling, metabolic function, and hypertrophy. Notably, the combinatorial overexpression of KLF15, ESRRA, and HOPX improved all three maturation parameters simultaneously. Taken together, we introduce a new TF cocktail that can be used in solo or in conjunction with other strategies to improve hPSC-CM maturation and anticipate that our generalizable methodology can also be implemented to identify maturation-associated TFs for other stem cell progenies.

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.

Single-cell transcriptomics for the assessment of cardiac disease

Cardiovascular disease is the leading cause of death globally. An advanced understanding of cardiovascular disease mechanisms is required to improve therapeutic strategies and patient risk stratification. State-of-the-art, large-scale, single-cell and single-nucleus transcriptomics facilitate the exploration of the cardiac cellular landscape at an unprecedented level, beyond its descriptive features, and can further our understanding of the mechanisms of disease and guide functional studies. In this Review, we provide an overview of the technical challenges in the experimental design of single-cell and single-nucleus transcriptomics studies, as well as a discussion of the type of inferences that can be made from the data derived from these studies. Furthermore, we describe novel findings derived from transcriptomics studies for each major cardiac cell type in both health and disease, and from development to adulthood. This Review also provides a guide to interpreting the exhaustive list of newly identified cardiac cell types and states, and highlights the consensus and discordances in annotation, indicating an urgent need for standardization. We describe advanced applications such as integration of single-cell data with spatial transcriptomics to map genes and cells on tissue and define cellular microenvironments that regulate homeostasis and disease progression. Finally, we discuss current and future translational and clinical implications of novel transcriptomics approaches, and provide an outlook of how these technologies will change the way we diagnose and treat heart disease.

Multimodal charting of molecular and functional cell states via in situ electro-sequencing

Paired mapping of single-cell gene expression and electrophysiology is essential to understand gene-to-function relationships in electrogenic tissues. Here, we developed in situ electro-sequencing (electro-seq) that combines flexible bioelectronics with in situ RNA sequencing to stably map millisecond-timescale electrical activity and profile single-cell gene expression from the same cells across intact biological networks, including cardiac and neural patches. When applied to human-induced pluripotent stem-cell-derived cardiomyocyte patches, in situ electro-seq enabled multimodal in situ analysis of cardiomyocyte electrophysiology and gene expression at the cellular level, jointly defining cell states and developmental trajectories. Using machine-learning-based cross-modal analysis, in situ electro-seq identified gene-to-electrophysiology relationships throughout cardiomyocyte development and accurately reconstructed the evolution of gene expression profiles based on long-term stable electrical measurements. In situ electro-seq could be applicable to create spatiotemporal multimodal maps in electrogenic tissues, potentiating the discovery of cell types and gene programs responsible for electrophysiological function and dysfunction.

Trajectory reconstruction identifies dysregulation of perinatal maturation programs in pluripotent stem cell-derived cardiomyocytes

A limitation in the application of pluripotent stem cell-derived cardiomyocytes (PSC-CMs) is the failure of these cells to achieve full functional maturity. The mechanisms by which directed differentiation differs from endogenous development, leading to consequent PSC-CM maturation arrest, remain unclear. Here, we generate a single-cell RNA sequencing (scRNA-seq) reference of mouse in vivo CM maturation with extensive sampling of previously difficult-to-isolate perinatal time periods. We subsequently generate isogenic embryonic stem cells to create an in vitro scRNA-seq reference of PSC-CM-directed differentiation. Through trajectory reconstruction, we identify an endogenous perinatal maturation program that is poorly recapitulated in vitro. By comparison with published human datasets, we identify a network of nine transcription factors (TFs) whose targets are consistently dysregulated in PSC-CMs across species. Notably, these TFs are only partially activated in common ex vivo approaches to engineer PSC-CM maturation. Our study can be leveraged toward improving the clinical viability of PSC-CMs.

Single-cell transcriptomes in the heart: when every epigenome counts

The response of an organ to stimuli emerges from the actions of individual cells. Recent cardiac single-cell RNA-sequencing studies of development, injury, and reprogramming have uncovered heterogeneous populations even among previously well-defined cell types, raising questions about what level of experimental resolution corresponds to disease-relevant, tissue-level phenotypes. In this review, we explore the biological meaning behind this cellular heterogeneity by undertaking an exhaustive analysis of single-cell transcriptomics in the heart (including a comprehensive, annotated compendium of studies published to date) and evaluating new models for the cardiac function that have emerged from these studies (including discussion and schematics that depict new hypotheses in the field). We evaluate the evidence to support the biological actions of newly identified cell populations and debate questions related to the role of cell-to-cell variability in development and disease. Finally, we present emerging epigenomic approaches that, when combined with single-cell RNA-sequencing, can resolve basic mechanisms of gene regulation and variability in cell phenotype.

Tissue-embedded stretchable nanoelectronics reveal endothelial cell-mediated electrical maturation of human 3D cardiac microtissues

Clinical translation of stem cell therapies for heart disease requires electrical integration of transplanted cardiomyocytes. Generation of electrically matured human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) is critical for electrical integration. Here, we found that hiPSC-derived endothelial cells (hiPSC-ECs) promoted the expression of selected maturation markers in hiPSC-CMs. Using tissue-embedded stretchable mesh nanoelectronics, we achieved a long-term stable map of human three-dimensional (3D) cardiac microtissue electrical activity. The results revealed that hiPSC-ECs accelerated the electrical maturation of hiPSC-CMs in 3D cardiac microtissues. Machine learning-based pseudotime trajectory inference of cardiomyocyte electrical signals further revealed the electrical phenotypic transition path during development. Guided by the electrical recording data, single-cell RNA sequencing identified that hiPSC-ECs promoted cardiomyocyte subpopulations with a more mature phenotype, and multiple ligand-receptor interactions were up-regulated between hiPSC-ECs and hiPSC-CMs, revealing a coordinated multifactorial mechanism of hiPSC-CM electrical maturation. Collectively, these findings show that hiPSC-ECs drive hiPSC-CM electrical maturation via multiple intercellular pathways.

Deep learning in spatial transcriptomics: Learning from the next next-generation sequencing

Spatial transcriptomics (ST) technologies are rapidly becoming the extension of single-cell RNA sequencing (scRNAseq), holding the potential of profiling gene expression at a single-cell resolution while maintaining cellular compositions within a tissue. Having both expression profiles and tissue organization enables researchers to better understand cellular interactions and heterogeneity, providing insight into complex biological processes that would not be possible with traditional sequencing technologies. Data generated by ST technologies are inherently noisy, high-dimensional, sparse, and multi-modal (including histological images, count matrices, etc.), thus requiring specialized computational tools for accurate and robust analysis. However, many ST studies currently utilize traditional scRNAseq tools, which are inadequate for analyzing complex ST datasets. On the other hand, many of the existing ST-specific methods are built upon traditional statistical or machine learning frameworks, which have shown to be sub-optimal in many applications due to the scale, multi-modality, and limitations of spatially resolved data (such as spatial resolution, sensitivity, and gene coverage). Given these intricacies, researchers have developed deep learning (DL)-based models to alleviate ST-specific challenges. These methods include new state-of-the-art models in alignment, spatial reconstruction, and spatial clustering, among others. However, DL models for ST analysis are nascent and remain largely underexplored. In this review, we provide an overview of existing state-of-the-art tools for analyzing spatially resolved transcriptomics while delving deeper into the DL-based approaches. We discuss the new frontiers and the open questions in this field and highlight domains in which we anticipate transformational DL applications.

Impaired Human Cardiac Cell Development due to NOTCH1 Deficiency

BACKGROUND

NOTCH1 pathogenic variants are implicated in multiple types of congenital heart defects including hypoplastic left heart syndrome, where the left ventricle is underdeveloped. It is unknown how NOTCH1 regulates human cardiac cell lineage determination and cardiomyocyte proliferation. In addition, mechanisms by which NOTCH1 pathogenic variants lead to ventricular hypoplasia in hypoplastic left heart syndrome remain elusive.

METHODS

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 genome editing was utilized to delete NOTCH1 in human induced pluripotent stem cells. Cardiac differentiation was carried out by sequential modulation of WNT signaling, and NOTCH1 knockout and wild-type differentiating cells were collected at day 0, 2, 5, 10, 14, and 30 for single-cell RNA-seq.

RESULTS

Human NOTCH1 knockout induced pluripotent stem cells are able to generate functional cardiomyocytes and endothelial cells, suggesting that NOTCH1 is not required for mesoderm differentiation and cardiovascular development in vitro. However, disruption of NOTCH1 blocks human ventricular-like cardiomyocyte differentiation but promotes atrial-like cardiomyocyte generation through shortening the action potential duration. NOTCH1 deficiency leads to defective proliferation of early human cardiomyocytes, and transcriptomic analysis indicates that pathways involved in cell cycle progression and mitosis are downregulated in NOTCH1 knockout cardiomyocytes. Single-cell transcriptomic analysis reveals abnormal cell lineage determination of cardiac mesoderm, which is manifested by the biased differentiation toward epicardial and second heart field progenitors at the expense of first heart field progenitors in NOTCH1 knockout cell populations.

CONCLUSIONS

NOTCH1 is essential for human ventricular-like cardiomyocyte differentiation and proliferation through balancing cell fate determination of cardiac mesoderm and modulating cell cycle progression. Because first heart field progenitors primarily contribute to the left ventricle, we speculate that pathogenic NOTCH1 variants lead to biased differentiation of first heart field progenitors, blocked ventricular-like cardiomyocyte differentiation, and defective cardiomyocyte proliferation, which collaboratively contribute to left ventricular hypoplasia in hypoplastic left heart syndrome.

Identification of SALL4 Expressing Islet-1+ Cardiovascular Progenitor Cell Clones

The utilization of cardiac progenitor cells (CPCs) has been shown to induce favorable regenerative effects. While there are various populations of endogenous CPCs in the heart, there is no consensus regarding which population is ideal for cell-based regenerative therapy. Early-stage progenitor cells can be differentiated into all cardiovascular lineages, including cardiomyocytes and endothelial cells. Identifying an Islet-1+ (Isl-1+) early-stage progenitor population with enhanced stemness, multipotency and differentiation potential would be beneficial for the development of novel regenerative therapies. Here, we investigated the transcriptome of human neonatal Isl-1+ CPCs. Isl-1+ human neonatal CPCs exhibit enhanced stemness properties and were found to express Spalt-like transcription factor 4 (SALL4). SALL4 plays a role in embryonic development as well as proliferation and expansion of hematopoietic progenitor cells. SALL4, SOX2, EpCAM and TBX5 are co-expressed in the majority of Isl-1+ clones isolated from neonatal patients. The pre-mesendodermal transcript TFAP2C was identified in select Isl-1, SALL4, SOX2, EpCAM and TBX5 expressing clones. The ability to isolate and expand pre-mesendodermal stage cells from human patients is a novel finding that holds potential value for applications in regenerative medicine.

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.

Single-cell and spatial transcriptomics: Advances in heart development and disease applications

Current transcriptomics technologies, including bulk RNA-seq, single-cell RNA sequencing (scRNA-seq), single-nucleus RNA-sequencing (snRNA-seq), and spatial transcriptomics (ST), provide novel insights into the spatial and temporal dynamics of gene expression during cardiac development and disease processes. Cardiac development is a highly sophisticated process involving the regulation of numerous key genes and signaling pathways at specific anatomical sites and developmental stages. Exploring the cell biological mechanisms involved in cardiogenesis also contributes to congenital heart disease research. Meanwhile, the severity of distinct heart diseases, such as coronary heart disease, valvular disease, cardiomyopathy, and heart failure, is associated with cellular transcriptional heterogeneity and phenotypic alteration. Integrating transcriptomic technologies in the clinical diagnosis and treatment of heart diseases will aid in advancing precision medicine. In this review, we summarize applications of scRNA-seq and ST in the cardiac field, including organogenesis and clinical diseases, and provide insights into the promise of single-cell and spatial transcriptomics in translational research and precision medicine.

Transcription Factors - the Essence of Heart Regeneration: A Potential Novel Therapeutic Strategy

Myocardial cell injury and following sequelae are the primary reasons for death globally. Unfortunately, myocardiocytes in adults have limited regeneration capacity. Therefore, the generation of neo myocardiocytes from non-myocardial cells is a surrogate strategy. Transcription factors (TFs) can be recruited to achieve this tremendous goal. Transcriptomic analyses have suggested that GATA, Mef2c, and Tbx5 (GMT cocktail) are master TFs to transdifferentiate/reprogram cell linage of fibroblasts, somatic cells, mesodermal cells into myocardiocytes. However, adding MESP1, MYOCD, ESRRG, and ZFPM2 TFs induces the generation of more efficient and physiomorphological features for induced myocardiocytes. Moreover, the same cocktail of transcription factors can induce the proliferation and differentiation of induced/pluripotent stem cells into myocardial cells. Amelioration of impaired myocardial cells involves the activation of healing transcription factors, which are induced by inflammation mediators; IL6, tumor growth factor β, and IL22. Transcription factors regulate the cellular and subcellular physiology of myocardiocytes to include mitotic cell cycling regulation, karyokinesis and cytokinesis, hypertrophic growth, adult sarcomeric contractile protein gene expression, fatty acid metabolism, and mitochondrial biogenesis and maturation. Cell therapy by transcription factors can be applied to cardiogenesis and ameliorating impaired cardiocytes. Transcription factors are the cornerstone in cell differentiation.

Surfaceome mapping of primary human heart cells with CellSurfer uncovers cardiomyocyte surface protein LSMEM2 and proteome dynamics in failing hearts

Cardiac cell surface proteins are drug targets and useful biomarkers for discriminating among cellular phenotypes and disease states. Here we developed an analytical platform, CellSurfer, that enables quantitative cell surface proteome (surfaceome) profiling of cells present in limited quantities, and we apply it to isolated primary human heart cells. We report experimental evidence of surface localization and extracellular domains for 1,144 N-glycoproteins, including cell-type-restricted and region-restricted glycoproteins. We identified a surface protein specific for healthy cardiomyocytes, LSMEM2, and validated an anti-LSMEM2 monoclonal antibody for flow cytometry and imaging. Surfaceome comparisons among pluripotent stem cell derivatives and their primary counterparts highlighted important differences with direct implications for drug screening and disease modeling. Finally, 20% of cell surface proteins, including LSMEM2, were differentially abundant between failing and non-failing cardiomyocytes. These results represent a rich resource to advance development of cell type and organ-specific targets for drug delivery, disease modeling, immunophenotyping and in vivo imaging.

Human-Induced Pluripotent Stem Cell-Based Differentiation of Cardiomyocyte Subtypes for Drug Discovery and Cell Therapy

Drug attrition rates have increased over the past few years, accompanied with growing costs for the pharmaceutical industry and consumers. Lack of in vitro models connecting the results of toxicity screening assays with clinical outcomes accounts for this high attrition rate. The emergence of cardiomyocytes derived from human pluripotent stem cells provides an amenable source of cells for disease modeling, drug discovery, and cardiotoxicity screening. Functionally similar to to embryonic stem cells, but with fewer ethical concerns, induced pluripotent stem cells (iPSCs) can recapitulate patient-specific genetic backgrounds, which would be a huge revolution for personalized medicine. The generated iPSC-derived cardiomyocytes (iPSC-CMs) represent different subtypes including ventricular-, atrial-, and nodal-like cardiomyocytes. Purifying these subtypes for chamber-specific drug screening presents opportunities and challenges. In this chapter, we discuss the strategies for the purification of iPSC-CMs, the use of iPSC-CMs for drug discovery and cardiotoxicity test, and the current limitations of iPSC-CMs that should be overcome for wider and more precise cardiovascular applications.

Integrative single-cell analysis of cardiogenesis identifies developmental trajectories and non-coding mutations in congenital heart disease

To define the multi-cellular epigenomic and transcriptional landscape of cardiac cellular development, we generated single-cell chromatin accessibility maps of human fetal heart tissues. We identified eight major differentiation trajectories involving primary cardiac cell types, each associated with dynamic transcription factor (TF) activity signatures. We contrasted regulatory landscapes of iPSC-derived cardiac cell types and their in vivo counterparts, which enabled optimization of in vitro differentiation of epicardial cells. Further, we interpreted sequence based deep learning models of cell-type-resolved chromatin accessibility profiles to decipher underlying TF motif lexicons. De novo mutations predicted to affect chromatin accessibility in arterial endothelium were enriched in congenital heart disease (CHD) cases vs. controls. In vitro studies in iPSCs validated the functional impact of identified variation on the predicted developmental cell types. This work thus defines the cell-type-resolved cis-regulatory sequence determinants of heart development and identifies disruption of cell type-specific regulatory elements in CHD.

Strategy of Patient-Specific Therapeutics in Cardiovascular Disease Through Single-Cell RNA Sequencing

Recently, single cell RNA sequencing (scRNA-seq) technology has enabled the discovery of novel or rare subtypes of cells and their characteristics. This technique has advanced unprecedented biomedical research by enabling the profiling and analysis of the transcriptomes of single cells at high resolution and throughput. Thus, scRNA-seq has contributed to recent advances in cardiovascular research by the generation of cell atlases of heart and blood vessels and the elucidation of mechanisms involved in cardiovascular development and diseases. This review summarizes the overall workflow of the scRNA-seq technique itself and key findings in the cardiovascular development and diseases based on the previous studies. In particular, we focused on how the single-cell sequencing technology can be utilized in clinical field and precision medicine to treat specific diseases.

Single-cell transcriptomic profiling reveals specific maturation signatures in human cardiomyocytes derived from LMNB2-inactivated induced pluripotent stem cells

Mammalian cardiomyocyte maturation entails phenotypic and functional optimization during the late fetal and postnatal phases of heart development, both processes driven and coordinated by complex gene regulatory networks. Cardiomyocytes derived from human induced pluripotent stem cells (iPSCs) are heterogenous and immature, barely resembling their adult in vivo counterparts. To characterize relevant developmental programs and maturation states during human iPSC-cardiomyocyte differentiation, we performed single-cell transcriptomic sequencing, which revealed six cardiomyocyte subpopulations, whose heterogeneity was defined by cell cycle and maturation states. Two of those subpopulations were characterized by a mature, non-proliferative transcriptional profile. To further investigate the proliferation-maturation transition in cardiomyocytes, we induced loss-of-function of LMNB2, which represses cell cycle progression in primary cardiomyocytes in vivo. This resulted in increased maturation in LMNB2-inactivated cardiomyocytes, characterized by transcriptional profiles related to myofibril structure and energy metabolism. Furthermore, we identified maturation signatures and maturational trajectories unique for control and LMNB2-inactivated cardiomyocytes. By comparing these datasets with single-cell transcriptomes of human fetal hearts, we were able to define spatiotemporal maturation states in human iPSC-cardiomyocytes. Our results provide an integrated approach for comparing in vitro-differentiated cardiomyocytes with their in vivo counterparts and suggest a strategy to promote cardiomyocyte maturation.

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.

Fluorescence Spectroscopy of Low-Level Endogenous β-adrenergic Receptor Expression at the Plasma Membrane of Differentiating Human iPSC-Derived Cardiomyocytes

The potential of human-induced pluripotent stem cells (hiPSCs) to be differentiated into cardiomyocytes (CMs) mimicking adult CMs functional morphology, marker genes and signaling characteristics has been investigated since over a decade. The evolution of the membrane localization of CM-specific G protein-coupled receptors throughout differentiation has received, however, only limited attention to date. We employ here advanced fluorescent spectroscopy, namely linescan Fluorescence Correlation Spectroscopy (FCS), to observe how the plasma membrane abundance of the β1- and β2-adrenergic receptors (β1/2-ARs), labelled using a bright and photostable fluorescent antagonist, evolves during the long-term monolayer culture of hiPSC-derived CMs. We compare it to the kinetics of observed mRNA levels in wildtype (WT) hiPSCs and in two CRISPR/Cas9 knock-in clones. We conduct these observations against the backdrop of our recent report of cell-to-cell expression variability, as well as of the subcellular localization heterogeneity of β-ARs in adult CMs.

Time-regulated transcripts with the potential to modulate human pluripotent stem cell-derived cardiomyocyte differentiation

Background

Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are a promising disease model, even though hiPSC-CMs cultured for extended periods display an undifferentiated transcriptional landscape. MiRNA–target gene interactions contribute to fine-tuning the genetic program governing cardiac maturation and may uncover critical pathways to be targeted.

Methods

We analyzed a hiPSC-CM public dataset to identify time-regulated miRNA–target gene interactions based on three logical steps of filtering. We validated this process in silico using 14 human and mouse public datasets, and further confirmed the findings by sampling seven time points over a 30-day protocol with a hiPSC-CM clone developed in our laboratory. We then added miRNA mimics from the top eight miRNAs candidates in three cell clones in two different moments of cardiac specification and maturation to assess their impact on differentiation characteristics including proliferation, sarcomere structure, contractility, and calcium handling.

Results

We uncovered 324 interactions among 29 differentially expressed genes and 51 miRNAs from 20,543 transcripts through 120 days of hiPSC-CM differentiation and selected 16 genes and 25 miRNAs based on the inverse pattern of expression (Pearson R-values < − 0.5) and consistency in different datasets. We validated 16 inverse interactions among eight genes and 12 miRNAs (Person R-values < − 0.5) during hiPSC-CMs differentiation and used miRNAs mimics to verify proliferation, structural and functional features related to maturation. We also demonstrated that miR-124 affects Ca²⁺ handling altering features associated with hiPSC-CMs maturation.

Conclusion

We uncovered time-regulated transcripts influencing pathways affecting cardiac differentiation/maturation axis and showed that the top-scoring miRNAs indeed affect primarily structural features highlighting their role in the hiPSC-CM maturation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-03138-x.

Acoustic Fabrication of Living Cardiomyocyte-based Hybrid Biorobots

Organized assemblies of cells have demonstrated promise as bioinspired actuators and devices; still, the fabrication of such "biorobots" has predominantly relied on passive assembly methods that reduce design capabilities. To address this, we have developed a strategy for the rapid formation of functional biorobots composed of live cardiomyocytes. We employ tunable acoustic fields to facilitate the efficient aggregation of millions of cells into high-density macroscopic architectures with directed cell orientation and enhanced cell-cell interaction. These biorobots can perform actuation functions both through naturally occurring contraction-relaxation cycles and through external control with chemical and electrical stimuli. We demonstrate that these biorobots can be used to achieve controlled actuation of a soft skeleton and pumping of microparticles. The biocompatible acoustic assembly strategy described here should prove generally useful for cellular manipulation in the context of tissue engineering, soft robotics, and other applications.

Advances in application of single-cell RNA sequencing in cardiovascular research

Single-cell RNA sequencing (scRNA-seq) provides high-resolution information on transcriptomic changes at the single-cell level, which is of great significance for distinguishing cell subtypes, identifying stem cell differentiation processes, and identifying targets for disease treatment. In recent years, emerging single-cell RNA sequencing technologies have been used to make breakthroughs regarding decoding developmental trajectories, phenotypic transitions, and cellular interactions in the cardiovascular system, providing new insights into cardiovascular disease. This paper reviews the technical processes of single-cell RNA sequencing and the latest progress based on single-cell RNA sequencing in the field of cardiovascular system research, compares single-cell RNA sequencing with other single-cell technologies, and summarizes the extended applications and advantages and disadvantages of single-cell RNA sequencing. Finally, the prospects for applying single-cell RNA sequencing in the field of cardiovascular research are discussed.

Dissecting the Heterogeneity of Human Thoracic Aortic Aneurysms Using Single-Cell Transcriptomics

Thoracic aortic aneurysm is a life-threatening condition caused by weakening of the thoracic aorta wall, often developing silently until dissection or rupture occurs. Despite substantial efforts in the past decade, there have been no significant therapeutic advances to prevent or clinically manage diverse forms of thoracic aortic aneurysm and dissection with the only effective treatment being surgical repair. There is an urgent need to understand intra- and inter-aneurysmal heterogeneity underlying thoracic aortic aneurysm and dissection pathogenesis. The human aortic wall consists of many cell types and exhibits significant regional heterogeneity. High-throughput single-cell RNA sequencing has emerged as the principal tool to reveal the complexity in human tissues and clinical specimens. Recent single-cell RNA sequencing studies of different aortic cell populations both in vivo and in vitro began to dissect this complexity and have provided valuable information. In this review, we summarize these findings and discuss the potential applications of single-cell transcriptomics and related high-content technologies in human thoracic aortic aneurysm and dissection research, as well as the challenges associated with it.

Corin: A Key Mediator in Sodium Homeostasis, Vascular Remodeling, and Heart Failure

Simple Summary

Atrial natriuretic peptide (ANP) is an important hormone that regulates many physiological and pathological processes, including electrolyte and body fluid balance, blood volume and pressure, cardiac channel activity and function, inflammatory response, lipid metabolism, and vascular remodeling. Corin is a transmembrane serine protease that activates ANP. Variants in the CORIN gene are associated with cardiovascular disease, including hypertension, cardiac hypertrophy, atrial fibrillation, heart failure, and preeclampsia. The current data indicate a key role of corin-mediated ANP production and signaling in the maintenance of cardiovascular homeostasis. In this review, we discuss the latest findings regarding the molecular and cellular mechanisms underlying the role of corin in sodium homeostasis, uterine spiral artery remodeling, and heart failure.

Abstract

Atrial natriuretic peptide (ANP) is a crucial element of the cardiac endocrine function that promotes natriuresis, diuresis, and vasodilation, thereby protecting normal blood pressure and cardiac function. Corin is a type II transmembrane serine protease that is highly expressed in the heart, where it converts the ANP precursor to mature ANP. Corin deficiency prevents ANP activation and causes hypertension and heart disease. In addition to the heart, corin is expressed in other tissues, including those of the kidney, skin, and uterus, where corin-mediated ANP production and signaling act locally to promote sodium excretion and vascular remodeling. These results indicate that corin and ANP function in many tissues via endocrine and autocrine mechanisms. In heart failure patients, impaired natriuretic peptide processing is a common pathological mechanism that contributes to sodium and body fluid retention. In this review, we discuss most recent findings regarding the role of corin in non-cardiac tissues, including the kidney and skin, in regulating sodium homeostasis and body fluid excretion. Moreover, we describe the molecular mechanisms underlying corin and ANP function in supporting orderly cellular events in uterine spiral artery remodeling. Finally, we assess the potential of corin-based approaches to enhance natriuretic peptide production and activity as a treatment of heart failure.

A novel Cbx1, PurB, and Sp3 complex mediates long-term silencing of tissue- and lineage-specific genes

miRNA-based cellular fate reprogramming offers an opportunity to investigate the mechanisms of long-term gene silencing. To further understand how genes are silenced in a tissue-specific manner, we leveraged our miRNA-based method of reprogramming fibroblasts into cardiomyocytes. Through screening approaches, we identified three proteins that were downregulated during reprogramming of fibroblasts into cardiomyocytes: heterochromatin protein Cbx1, transcriptional activator protein PurB, and transcription factor Sp3. We show that knockdown of Cbx1, PurB, and Sp3 was sufficient to induce cardiomyocyte gene expression in fibroblasts. Similarly, gene editing to ablate Cbx1, PurB, and Sp3 expression induced fibroblasts to convert into cardiomyocytes in vivo. Furthermore, high-throughput DNA sequencing and coimmunoprecipitation experiments indicated that Cbx1, PurB, and Sp3 also bound together as a complex and were necessary to localize nucleosomes to cardiomyocyte genes on the chromosome. Finally, we found that the expression of these genes led to nucleosome modification via H3K27me3 (trimethylated histone-H3 lysine-27) deposition through an interaction with the polycomb repressive PRC2 complex. In summary, we conclude that Cbx1, PurB, and Sp3 control cell fate by actively repressing lineage-specific genes.

Computational profiling of hiPSC-derived heart organoids reveals chamber defects associated with NKX2-5 deficiency

Heart organoids have the potential to generate primary heart-like anatomical structures and hold great promise as in vitro models for cardiac disease. However, their properties have not yet been fully studied, which hinders their wide spread application. Here we report the development of differentiation systems for ventricular and atrial heart organoids, enabling the study of heart diseases with chamber defects. We show that our systems generate chamber-specific organoids comprising of the major cardiac cell types, and we use single cell RNA sequencing together with sample multiplexing to characterize the cells we generate. To that end, we developed a machine learning label transfer approach leveraging cell type, chamber, and laterality annotations available for primary human fetal heart cells. We then used this model to analyze organoid cells from an isogeneic line carrying an Ebstein’s anomaly associated genetic variant in NKX2-5, and we successfully recapitulated the disease’s atrialized ventricular defects. In summary, we have established a workflow integrating heart organoids and computational analysis to model heart development in normal and disease states.

A human cardiac organoid system, coupled with single cell RNA sequencing and machine learning for transcriptional phenotyping, was developed. This allowed investigation of a genetic variant associated with Ebstein’s Anomaly, a congenital heart disease with chamber defects.

Deciphering Cardiac Biology and Disease by Single-Cell Transcriptomic Profiling

By detecting minute molecular changes in hundreds to millions of single cells, single-cell RNA sequencing allows for the comprehensive characterization of the diversity and dynamics of cells in the heart. Our understanding of the heart has been transformed through the recognition of cellular heterogeneity, the construction of regulatory networks, the building of lineage trajectories, and the mapping of intercellular crosstalk. In this review, we introduce cardiac progenitors and their transcriptional regulation during embryonic development, highlight cellular heterogeneity and cell subtype functions in cardiac health and disease, and discuss insights gained from the study of pluripotent stem-cell-derived cardiomyocytes.

Uncompensated mitochondrial oxidative stress underlies heart failure in an iPSC-derived model of congenital heart disease

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease with 30% mortality from heart failure (HF) in the first year of life, but the cause of early HF remains unknown. Induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CM) from patients with HLHS showed that early HF is associated with increased apoptosis, mitochondrial respiration defects, and redox stress from abnormal mitochondrial permeability transition pore (mPTP) opening and failed antioxidant response. In contrast, iPSC-CM from patients without early HF showed normal respiration with elevated antioxidant response. Single-cell transcriptomics confirmed that early HF is associated with mitochondrial dysfunction accompanied with endoplasmic reticulum (ER) stress. These findings indicate that uncompensated oxidative stress underlies early HF in HLHS. Importantly, mitochondrial respiration defects, oxidative stress, and apoptosis were rescued by treatment with sildenafil to inhibit mPTP opening or TUDCA to suppress ER stress. Together these findings point to the potential use of patient iPSC-CM for modeling clinical heart failure and the development of therapeutics.