Physiology of Cardiac Development: From Genetics to Signaling to Therapeutic Strategies

Single-cell reconstruction and mutation enrichment analysis identifies dysregulated cardiomyocyte and endothelial cells in congenital heart disease

Congenital heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. To catalog the putative candidate CHD risk genes, we collected 16,349 variants [single-nucleotide variants (SNVs) and Indels] impacting 8,308 genes in 3,166 CHD cases for a comprehensive meta-analysis. Using American College of Medical Genetics (ACMG) guidelines, we excluded the 0.1% of benign/likely benign variants and the resulting dataset consisted of 83% predicted loss of function variants and 17% missense variants. Seventeen percent were de novo variants. A stepwise analysis identified 90 variant-enriched CHD genes, of which six (GPATCH1, NYNRIN, TCLD2, CEP95, MAP3K19, and TTC36) were novel candidate CHD genes. Single-cell transcriptome cluster reconstruction analysis on six CHD tissues and four controls revealed upregulation of the top 10 frequently mutated genes primarily in cardiomyocytes. NOTCH1 (highest number of variants) and MYH6 (highest number of recurrent variants) expression was elevated in endocardial cells and cardiomyocytes, respectively, and 60% of these gene variants were associated with tetralogy of Fallot and coarctation of the aorta, respectively. Pseudobulk analysis using the single-cell transcriptome revealed significant (P < 0.05) upregulation of both NOTCH1 (endocardial cells) and MYH6 (cardiomyocytes) in the control heart data. We observed nine different subpopulations of CHD heart cardiomyocytes of which only four were observed in the control heart. This is the first comprehensive meta-analysis combining genomics and CHD single-cell transcriptomics, identifying the most frequently mutated CHD genes, and demonstrating CHD gene heterogeneity, suggesting that multiple genes contribute to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.NEW & NOTEWORTHY Congential heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. We present a comprehensive analysis combining genomics and CHD single-cell transcriptome. Our study identifies 90 potential candidate CHD risk genes of which 6 are novel. The risk genes have heterogenous expression suggestive of multiple genes contributing to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.

The Activation of the LIMK/Cofilin Signaling Pathway via Extracellular Matrix-Integrin Interactions Is Critical for the Generation of Mature and Vascularized Cardiac Organoids

The generation of mature and vascularized human pluripotent stem cell-derived cardiac organoids (hPSC-COs) is necessary to ensure the validity of drug screening and disease modeling. This study investigates the effects of cellular aggregate (CA) stemness and self-organization on the generation of mature and vascularized hPSC-COs and elucidates the mechanisms underlying cardiac organoid (CO) maturation and vascularization. COs derived from 2-day-old CAs with high stemness (H-COs) and COs derived from 5-day-old CAs with low stemness (L-COs) were generated in a self-organized microenvironment via Wnt signaling induction. This study finds that H-COs exhibit ventricular, structural, metabolic, and functional cardiomyocyte maturation and vessel networks consisting of endothelial cells, smooth muscle cells, pericytes, and basement membranes compared to L-COs. Transcriptional profiling shows the upregulation of genes associated with cardiac maturation and vessel formation in H-COs compared with the genes in L-COs. Through experiments with LIMK inhibitors, the activation of ROCK-LIMK-pCofilin via ECM-integrin interactions leads to cardiomyocyte maturation and vessel formation in H-COs. Furthermore, the LIMK/Cofilin signaling pathway induces TGFβ/NODAL and PDGF pathway activation for the maturation and vascularization of H-COs. The study demonstrates for the first time that LIMK/Cofilin axis activation plays an important role in the generation of mature and vascularized COs.

Evaluation of the effect of boron derivatives on cardiac differentiation of mouse pluripotent stem cells

BACKGROUND

The heart is one of the first organs to form during embryonic development and has a very important place. So much that the formation of a functional heart is completed on the 55th day of human development and the 15th day of mouse development. Myocardial, endocardial and epicardial cells, which are derived from the mesoderm layer, are the cells that form the basis of the heart. Cardiac development, like other embryonic developments, is tightly controlled and regulated by various signaling pathways. The WNT signaling pathway is the most studied of these signaling pathways and the one with the clearest relationship with heart development. It is known that boron compounds and the Wnt/β-catenin pathway are highly correlated. Therefore, this study aimed to investigate the role of boron compounds in heart development as well as its effect on pluripotency of mouse embryonic stem cells for the first time in the literature.

METHODS

Toxicity of boron compounds was evaluated by using MTS analysis and obtained results were supported by morphological pictures, Trypan Blue staining and Annexin V staining. Additionally, the possible boron-related change in pluripotency of embryonic stem cells were analyzed with alkaline phosphatase activity and immunocytochemical staining of Oct4 protein as well as gene expression levels of pluripotency related OCT4, SOX2 and KLF4 genes. The alterations in the embryonic body formation capacity of mouse embryonic stem cells due to the application boron derivatives were also evaluated. Three linage differentiation was conducted to clarify the real impact of boron compounds on embryonic development. Lastly, cardiac differentiation of mESCs was investigated by using morphological pictures, cytosolic calcium measurement, gene expression and immunocytochemical analysis of cardiac differentiation related genes and in the presence of boron compounds.

RESULTS

Obtained results show that boron treatment maintains the pluripotency of embryonic stem cells at non-toxic concentrations. Additionally, endodermal, and mesodermal fate was found to be triggered after boron treatment. Also, initiation of cardiomyocyte differentiation by boron derivative treatments caused an increased gene expression levels of cardiac differentiation related TNNT2, Nkx2.5 and ISL-1 gene expression levels.

CONCLUSION

This study indicates that boron application, which is responsible for maintaining pluripotency of mESCs, can be used for increased cardiomyocyte differentiation of mESCs.

The MEF2A transcription factor interactome in cardiomyocytes

Transcriptional regulators encoded by the Myocyte Enhancer Factor 2 (MEF2) gene family play a fundamental role in cardiac development, homeostasis and pathology. Previous studies indicate that MEF2A protein-protein interactions serve as a network hub in several cardiomyocyte cellular processes. Based on the idea that interactions with regulatory protein partners underly the diverse roles of MEF2A in cardiomyocyte gene expression, we undertook a systematic unbiased screen of the MEF2A protein interactome in primary cardiomyocytes using an affinity purification-based quantitative mass spectrometry approach. Bioinformatic processing of the MEF2A interactome revealed protein networks involved in the regulation of programmed cell death, inflammatory responses, actin dynamics and stress signaling in primary cardiomyocytes. Further biochemical and functional confirmation of specific protein-protein interactions documented a dynamic interaction between MEF2A and STAT3 proteins. Integration of transcriptome level data from MEF2A and STAT3-depleted cardiomyocytes reveals that the balance between MEF2A and STAT3 activity exerts a level of executive control over the inflammatory response and cardiomyocyte cell survival and experimentally ameliorates Phenylephrine induced cardiomyocyte hypertrophy. Lastly, we identified several MEF2A/STAT3 co-regulated genes, including the MMP9 gene. Herein, we document the cardiomyocyte MEF2A interactome, which furthers our understanding of protein networks involved in the hierarchical control of normal and pathophysiological cardiomyocyte gene expression in the mammalian heart.

Adaptation of cardiomyogenesis to the generation and maturation of cardiomyocytes from human pluripotent stem cells

The advent of methods for efficient generation and cardiac differentiation of pluripotent stem cells opened new avenues for disease modelling, drug testing, and cell therapies of the heart. However, cardiomyocytes (CM) obtained from such cells demonstrate an immature, foetal-like phenotype that involves spontaneous contractions, irregular morphology, expression of embryonic isoforms of sarcomere components, and low level of ion channels. These and other features may affect cellular response to putative therapeutic compounds and the efficient integration into the host myocardium after in vivo delivery. Therefore, novel strategies to increase the maturity of pluripotent stem cell-derived CM are of utmost importance. Several approaches have already been developed relying on molecular changes that occur during foetal and postnatal maturation of the heart, its electromechanical activity, and the cellular composition. As a better understanding of these determinants may facilitate the generation of efficient protocols for in vitro acquisition of an adult-like phenotype by immature CM, this review summarizes the most important molecular factors that govern CM during embryonic development, postnatal changes that trigger heart maturation, as well as protocols that are currently used to generate mature pluripotent stem cell-derived cardiomyocytes.

A genome-wide association study of obstructive heart defects among participants in the National Birth Defects Prevention Study

Obstructive heart defects (OHDs) share common structural lesions in arteries and cardiac valves, accounting for ~25% of all congenital heart defects. OHDs are highly heritable, resulting from interplay among maternal exposures, genetic susceptibilities, and epigenetic phenomena. A genome-wide association study was conducted in National Birth Defects Prevention Study participants (N_discovery  = 3978; N_replication  = 2507), investigating the genetic architecture of OHDs using transmission/disequilibrium tests (TDT) in complete case-parental trios (N_{discovery_TDT}  = 440; N_{replication_TDT}  = 275) and case-control analyses separately in infants (N_{discovery_CCI}  = 1635; N_{replication_CCI}  = 990) and mothers (case status defined by infant; N_{discovery_CCM}  = 1703; N_{replication_CCM}  = 1078). In the TDT analysis, the SLC44A2 single nucleotide polymorphism (SNP) rs2360743 was significantly associated with OHD (p_discovery  = 4.08 × 10^-9 ; p_replication  = 2.44 × 10^-4 ). A CAPN11 SNP (rs55877192) was suggestively associated with OHD (p_discovery  = 1.61 × 10^-7 ; p_replication  = 0.0016). Two other SNPs were suggestively associated (p < 1 × 10^-6 ) with OHD in only the discovery sample. In the case-control analyses, no SNPs were genome-wide significant, and, even with relaxed thresholds ( × _discovery  < 1 × 10^-5 and p_replication  < 0.05), only one SNP (rs188255766) in the infant analysis was associated with OHDs (p_discovery  = 1.42 × 10^-6 ; p_replication  = 0.04). Additional SNPs with p_discovery  < 1 × 10^-5 were in loci supporting previous findings but did not replicate. Overall, there was modest evidence of an association between rs2360743 and rs55877192 and OHD and some evidence validating previously published findings.

Expression Signatures of Long Noncoding RNAs in Left Ventricular Noncompaction

Long noncoding RNAs have gained widespread attention in recent years for their crucial role in biological regulation. They have been implicated in a range of developmental processes and diseases including cancer, cardiovascular, and neuronal diseases. However, the role of long noncoding RNAs (lncRNAs) in left ventricular noncompaction (LVNC) has not been explored. In this study, we investigated the expression levels of lncRNAs in the blood of LVNC patients and healthy subjects to identify differentially expressed lncRNA that develop LVNC specific biomarkers and targets for developing therapies using biological pathways. We used Agilent Human lncRNA array that contains both updated lncRNAs and mRNAs probes. We identified 1,568 upregulated and 1,141 downregulated (log fold-change > 2.0) lncRNAs that are differentially expressed between LVNC and the control group. Among them, RP11-1100L3.7 and XLOC_002730 are the most upregulated and downregulated lncRNAs. Using quantitative real-time reverse transcription polymerase chain reaction (RT-QPCR), we confirmed the differential expression of three top upregulated and downregulated lncRNAs along with two other randomly picked lncRNAs. Gene Ontology (GO) and KEGG pathways analysis with these differentially expressed lncRNAs provide insight into the cellular pathway leading to LVNC pathogenesis. We also identified 1,066 upregulated and 1,017 downregulated mRNAs. Gene set enrichment analysis (GSEA) showed that G2M, Estrogen, and inflammatory pathways are enriched in differentially expressed genes (DEG). We also identified miRNA targets for these differentially expressed genes. In this study, we first report the use of LncRNA microarray to understand the pathogenesis of LVNC and to identify several lncRNA and genes and their targets as potential biomarkers.

Dinitramine induces cardiotoxicity and morphological alterations on zebrafish embryo development

Dinitramine (DN), an herbicide in the dinitroaniline family, is used in agricultural areas to prevent unwanted plant growth. Dinitroaniline herbicides inhibit cell division by preventing microtubulin synthesis. They are strongly absorbed by the soil and can contaminate groundwater; however, the mode of action of these herbicides in non-target organisms remains unclear. In this study, we examined the developmental toxicity of DN in zebrafish embryos exposed to 1.6, 3.2, and 6.4 mg/L DN, compared to embryos exposed to DMSO (control) for 96 h. Visual assessments using transgenic zebrafish (fli1:eGFP) indicated abnormal cardiac development with enlarged ventricles and atria, decreased heartbeats, and impaired cardiac function. Along with cardiac development, vessel formation and angiogenesis were suppressed through activation of the inflammatory response. In addition, exposure to 6.4 mg/L DN for 96 h induced cell death, with upregulation of genes related to apoptosis. Our results showed that DN induced morphological changes and triggered an inflammatory response and apoptotic cell death that can impair embryonic growth and survival, providing an important mechanism of DN in aquatic organisms.

Amisulbrom causes cardiovascular toxicity in zebrafish (Danio rerio)

Amisulbrom (AML), a sulfonamide fungicide used to control oomycete diseases, is regarded as a threat to aquatic species. The objective of this study was to evaluate the potential effects of AML on fish using a zebrafish model. Zebrafish embryos were exposed to 0.0075 μM, 0.075 μM, and 0.75 μM AML. AML-treated zebrafish embryos exhibited severe developmental defects, including pericardial edema, blood-clot clustering, increased hatching rates, decreased heart rates, and abnormal hemoglobin distributions. Compared with controls, key marker genes associated with cardiovascular development (i.e., nkx2.5, myh6, myh7, myl7, alas2, hbbe1, hbbe2, and gata1a) were abnormally expressed in response to AML treatment, suggesting that AML might specifically affect cardiovascular development. These results provide a valuable reference for the effects of AML on zebrafish embryos and may help to further clarify the potential risks posed by AML to aquatic ecosystems.

LEFTY-PITX2 signaling pathway is critical for generation of mature and ventricular cardiac organoids in human pluripotent stem cell-derived cardiac mesoderm cells

The generation of mature ventricular cardiomyocytes (CMs) resembling adult CMs from human pluripotent stem cells (hPSCs) is necessary for disease modeling and drug discovery. To investigate the effect of self-organizing capacity on the generation of mature cardiac organoids (COs), we generated cardiac mesoderm cell-derived COs (CMC-COs) and CM-derived COs (CM-COs) and evaluated COs. CMC-COs exhibited more organized sarcomere structures and mitochondria, well-arranged t-tubule structures, and evenly distributed intercalated discs. Increased expressions of ventricular CM, cardiac metabolic, t-tubule formation, K+ ion channel, and junctional markers were confirmed in CMC-COs. Mature ventricular-like function such as faster motion vector speed, decreased beats per min, increased peak-to-peak duration, and prolonged APD50 and APD90 were observed in CMC-COs. Transcriptional profiling revealed that extracellular matrix-integrin, focal adhesion, and LEFTY-PITX2 signaling pathways are upregulated in CMC-COs. LEFTY knockdown affected ECM-integrin-FA signaling pathways in CMC-COs. Here, we found that high self-organizing capacity of CMCs is critical for the generation of mature and ventricular COs. We also demonstrated that LEFTY-PITX2 signaling plays key roles for CM maturation and specification into ventricular-like CM subtype in CMC-COs. CMC-COs are an attractive resource for disease modeling and drug discovery.

Chemical Genetics Reveals a Role of Squalene Synthase in TGFβ Signaling and Cardiomyogenesis

KY02111 is a widely used small molecule that boosts cardiomyogenesis of the mesoderm cells derived from pluripotent stem cells, yet its molecular mechanism of action remains elusive. The present study resolves the initially perplexing effects of KY02111 on Wnt signaling and subsequently identifies squalene synthase (SQS) as a molecular target of KY02111 and its optimized version, KY-I. By disrupting the interaction of SQS with cardiac ER-membrane protein TMEM43, KY02111 impairs TGFβ signaling, but not Wnt signaling, and thereby recapitulates the clinical mutation of TMEM43 that causes arrhythmogenic right ventricular cardiomyopathy (ARVC), an inherited heart disease that involves a substitution of myocardium with fatty tissue. These findings reveal a heretofore undescribed role of SQS in TGFβ signaling and cardiomyogenesis. KY02111 may find its use in ARVC modeling as well as serve as a chemical tool for studying TGFβ/SMAD signaling.

Knockdown of DNA methyltransferase 1 reduces DNA methylation and alters expression patterns of cardiac genes in embryonic cardiomyocytes

We previously found that DNA methyltransferase 3a (DNMT3a) plays an important role in regulating embryonic cardiomyocyte gene expression, morphology, and function. In this study, we investigated the role of the most abundant DNMT in mammalian cells, DNMT1, in these processes. It is known that DNMT1 is essential for embryonic development, during which it is involved in regulating cardiomyocyte DNA methylation and gene expression. We used siRNA to knock down DNMT1 expression in primary cultures of mouse embryonic cardiomyocytes. Immunofluorescence staining and multielectrode array were, respectively, utilized to evaluate cardiomyocyte growth and electrophysiology. RNA sequencing (RNA‐Seq) and multiplex bisulfite sequencing were, respectively, performed to examine gene expression and promoter methylation. At 72 h post‐transfection, reduction of DNMT1 expression decreased the number and increased the size of embryonic cardiomyocytes. Beat frequency and the amplitude of field action potentials were decreased by DNMT1 siRNA. RNA‐Seq analysis identified 801 up‐regulated genes and 494 down‐regulated genes in the DNMT1 knockdown cells when compared to controls. Pathway analysis of the differentially expressed genes revealed pathways that were associated with cell death and survival, cell morphology, cardiac function, and cardiac disease. Alternative splicing analysis identified 929 differentially expressed exons, including 583 up‐regulated exons and 308 down‐regulated exons. Moreover, decreased methylation levels were found in the promoters of cardiac genes Myh6, Myh7, Myh7b, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb, mef2c, mef2d, Camta2, Cdkn1A, and Cdkn1C. Of these 13 genes, 6 (Myh6, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb) and 1 (Cdkn1C) had increased or decreased gene expression, respectively. Altogether, these data show that DNMT1 is important in embryonic cardiomyocytes by regulating DNA methylation, gene expression, gene splicing, and cell function.

Telomeres and Telomerase in Heart Ontogenesis, Aging and Regeneration

The main purpose of the review article is to assess the contributions of telomere length and telomerase activity to the cardiac function at different stages of development and clarify their role in cardiac disorders. It has been shown that the telomerase complex and telomeres are of great importance in many periods of ontogenesis due to the regulation of the proliferative capacity of heart cells. The review article also discusses the problems of heart regeneration and the identification of possible causes of dysfunction of telomeres and telomerase.

scaRNA1 Levels Alter Pseudouridylation in Spliceosomal RNA U2 Affecting Alternative mRNA Splicing and Embryonic Development

The heart is the first major organ to develop during embryogenesis and must receive proper spatiotemporal signaling for proper development. Failure of proper signaling between the first and second heart fields at twenty days gestation contributes to the generation of a congenital heart defect. The most common cyanotic congenital heart defect is tetralogy of Fallot (TOF) which requires surgical intervention in the first year of life. In right ventricular tissue of infants born with TOF, the levels of scaRNA1 are reduced and mRNA splicing is dysregulated. In this study, we investigate a method of quantifying pseudouridylation levels in relation to scaRNA1 levels in spliceosomal RNA U2 in three different groups of samples: right ventricular (RV) tissue of infants born with TOF versus RV tissue from normally developing infants, scaRNA1 knockdown in primary normal cardiomyocytes derived from normally developing infants, and scaRNA1 overexpression in primary cells derived from RV tissue from infants born with TOF. We hypothesize that the amount of pseudouridylation is dependent on scaRNA1 level, compromising spliceosomal function and therefore, contributing to the generation of a congenital heart defect. Our results revealed a statistically significant decrease of pseudouridylation levels in the right ventricular tissue of infants born with TOF compared to the controls. Knocking down the scaRNA1 levels in normal primary cardiomyocytes resulted in a statistically significant decrease of pseudouridylation. Finally, an overexpression of scaRNA1 in TOF primary cells resulted in an increase in pseudouridylation levels, but it did not achieve statistical significance. Our previous research provided an association between scaRNA levels, alternative splicing, and development. Here, we demonstrate that pseudouridylation levels in spliceosomal RNA U2 is dependent on the expression level of scaRNA1. Although further investigation is needed, we believe that scaRNA expression regulates biochemical modifications to spliceosomal RNAs, adjusting the fidelity of the spliceosome, allowing for controlled alternative splicing of mRNA that is important in embryonic development. If validated, this is an underappreciated mechanism that is critical for regulating proper embryonic development.

Generation of human induced pluripotent stem cell-derived cardiomyocytes in 2D monolayer and scalable 3D suspension bioreactor cultures with reduced batch-to-batch variations

Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) are promising candidates to treat myocardial infarction and other cardiac diseases. Such treatments require pure cardiomyocytes (CMs) in large quantities. Methods: In the present study we describe an improved protocol for production of hiPSC-CMs in which hiPSCs are first converted into mesodermal cells by stimulation of wingless (Wnt) signaling using CHIR99021, which are then further differentiated into CM progenitors by simultaneous inhibition of porcupine and tankyrase pathways using IWP2 and XAV939 under continuous supplementation of ascorbate during the entire differentiation procedure. Results: The protocol resulted in reproducible generation of >90% cardiac troponin T (TNNT2)-positive cells containing highly organized sarcomeres. In 2D monolayer cultures CM yields amounted to 0.5 million cells per cm2 growth area, and on average 72 million cells per 100 mL bioreactor suspension culture without continuous perfusion. The differentiation efficiency was hardly affected by the initial seeding density of undifferentiated hiPSCs. Furthermore, batch-to-batch variations were reduced by combinatorial use of ascorbate, IWP2, and XAV939. Conclusion: Combined inhibition of porcupine and tankyrase sub-pathways of Wnt signaling and continuous ascorbate supplementation, enable robust and efficient production of hiPSC-CMs.

Cardiomyogenesis Modeling Using Pluripotent Stem Cells: The Role of Microenvironmental Signaling

Pluripotent stem cells (PSC) can be used as a model to study cardiomyogenic differentiation. In vitro modeling can reproduce cardiac development through modulation of some key signaling pathways. Therefore, many studies make use of this strategy to better understand cardiomyogenesis complexity and to determine possible ways to modulate cell fate. However, challenges remain regarding efficiency of differentiation protocols, cardiomyocyte (CM) maturation and therapeutic applications. Considering that the extracellular milieu is crucial for cellular behavior control, cardiac niche studies, such as those identifying secreted molecules from adult or neonatal tissues, allow the identification of extracellular factors that may contribute to CM differentiation and maturation. This review will focus on cardiomyogenesis modeling using PSC and the elements involved in cardiac microenvironmental signaling (the secretome - extracellular vesicles, extracellular matrix and soluble factors) that may contribute to CM specification and maturation.