Wnt dose escalation during the exit from pluripotency identifies tranilast as a regulator of cardiac mesoderm

Wnt signaling is a critical determinant of cell lineage development. This study used Wnt dose-dependent induction programs to gain insights into molecular regulation of stem cell differentiation. We performed single-cell RNA sequencing of hiPSCs responding to a dose escalation protocol with Wnt agonist CHIR-99021 during the exit from pluripotency to identify cell types and genetic activity driven by Wnt stimulation. Results of activated gene sets and cell types were used to build a multiple regression model that predicts the efficiency of cardiomyocyte differentiation. Cross-referencing Wnt-associated gene expression profiles to the Connectivity Map database, we identified the small-molecule drug, tranilast. We found that tranilast synergistically activates Wnt signaling to promote cardiac lineage differentiation, which we validate by in vitro analysis of hiPSC differentiation and in vivo analysis of developing quail embryos. Our study provides an integrated workflow that links experimental datasets, prediction models, and small-molecule databases to identify drug-like compounds that control cell differentiation.

High-content fluorescence bioassay investigates pore formation, ion channel modulation and cell membrane lysis induced by venoms

Venoms comprise highly sophisticated bioactive molecules modulating ion channels, receptors, coagulation factors, and the cellular membranes. This array of targets and bioactivities requires advanced high-content bioassays to facilitate the development of novel envenomation treatments and biotechnological and pharmacological agents. In response to the existing gap in venom research, we developed a cutting-edge fluorescence-based high-throughput and high-content cellular assay. This assay enables the simultaneous identification of prevalent cellular activities induced by venoms such as membrane lysis, pore formation, and ion channel modulation. By integrating intracellular calcium with extracellular nucleic acid measurements, we have successfully distinguished these venom mechanisms within a single cellular assay. Our high-content bioassay was applied across three cell types exposed to venom components representing lytic, ion pore-forming or ion channel modulator toxins. Beyond unveiling distinct profiles for these action mechanisms, we found that the pore-forming latrotoxin α-Lt1a prefers human neuroblastoma to kidney cells and cardiomyocytes, while the lytic bee peptide melittin is not selective. Furthermore, evaluation of snake venoms showed that Elapid species induced rapid membrane lysis, while Viper species showed variable to no activity on neuroblastoma cells. These findings underscore the ability of our high-content bioassay to discriminate between clades and interspecific traits, aligning with clinical observations at venom level, beyond discriminating among ion pore-forming, membrane lysis and ion channel modulation. We hope our research will expedite the comprehension of venom biology and the diversity of toxins that elicit cytotoxic, cardiotoxic and neurotoxic effects, and assist in identifying venom components that hold the potential to benefit humankind.

HOPX-associated molecular programs control cardiomyocyte cell states underpinning cardiac structure and function

Genomic regulation of cardiomyocyte differentiation is central to heart development and function. This study uses genetic loss-of-function human-induced pluripotent stem cell-derived cardiomyocytes to evaluate the genomic regulatory basis of the non-DNA-binding homeodomain protein HOPX. We show that HOPX interacts with and controls cardiac genes and enhancer networks associated with diverse aspects of heart development. Using perturbation studies in vitro, we define how upstream cell growth and proliferation control HOPX transcription to regulate cardiac gene programs. We then use cell, organoid, and zebrafish regeneration models to demonstrate that HOPX-regulated gene programs control cardiomyocyte function in development and disease. Collectively, this study mechanistically links cell signaling pathways as upstream regulators of HOPX transcription to control gene programs underpinning cardiomyocyte identity and function.

A village in a dish model system for population-scale hiPSC studies

The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and different conditions. Human induced pluripotent stem cells are uniquely suited to study these context-dependent effects but cell lines from hundreds or thousands of individuals are required. Village cultures, where multiple induced pluripotent stem lines are cultured and differentiated in a single dish, provide an elegant solution for scaling induced pluripotent stem experiments to the necessary sample sizes required for population-scale studies. Here, we show the utility of village models, demonstrating how cells can be assigned to an induced pluripotent stem line using single-cell sequencing and illustrating that the genetic, epigenetic or induced pluripotent stem line-specific effects explain a large percentage of gene expression variation for many genes. We demonstrate that village methods can effectively detect induced pluripotent stem line-specific effects, including sensitive dynamics of cell states.

New Drug Targets and Preclinical Modelling Recommendations for Treating Acute Myocardial Infarction

Acute myocardial infarction (AMI) is the leading cause of morbidity and mortality worldwide and the primary underlying risk factor for heart failure. Despite decades of research and clinical trials, there are no drugs currently available to prevent organ damage from acute ischaemic injuries of the heart. In order to address the increasing global burden of heart failure, drug, gene, and cell-based regeneration technologies are advancing into clinical testing. In this review we highlight the burden of disease associated with AMI and the therapeutic landscape based on market analyses. New studies revealing the role of acid-sensitive cardiac ion channels and other proton-gated ion channels in cardiac ischaemia are providing renewed interest in pre- and post-conditioning agents with novel mechanisms of action that may also have implications for gene- and cell-based therapeutics. Furthermore, we present guidelines that couple new cell technologies and data resources with traditional animal modelling pipelines to help de-risk drug candidates aimed at treating AMI. We propose that improved preclinical pipelines and increased investment in drug target identification for AMI is critical to stem the increasing global health burden of heart failure.

Vascular cells improve functionality of human cardiac organoids

Crosstalk between cardiac cells is critical for heart performance. Here we show that vascular cells within human cardiac organoids (hCOs) enhance their maturation, force of contraction, and utility in disease modeling. Herein we optimize our protocol to generate vascular populations in addition to epicardial, fibroblast, and cardiomyocyte cells that self-organize into in-vivo-like structures in hCOs. We identify mechanisms of communication between endothelial cells, pericytes, fibroblasts, and cardiomyocytes that ultimately contribute to cardiac organoid maturation. In particular, (1) endothelial-derived LAMA5 regulates expression of mature sarcomeric proteins and contractility, and (2) paracrine platelet-derived growth factor receptor β (PDGFRβ) signaling from vascular cells upregulates matrix deposition to augment hCO contractile force. Finally, we demonstrate that vascular cells determine the magnitude of diastolic dysfunction caused by inflammatory factors and identify a paracrine role of endothelin driving dysfunction. Together this study highlights the importance and role of vascular cells in organoid models.

Dynamic chromatin organization and regulatory interactions in human endothelial cell differentiation

Vascular endothelial cells are a mesoderm-derived lineage with many essential functions, including angiogenesis and coagulation. The gene-regulatory mechanisms underpinning endothelial specialization are largely unknown, as are the roles of chromatin organization in regulating endothelial cell transcription. To investigate the relationships between chromatin organization and gene expression, we induced endothelial cell differentiation from human pluripotent stem cells and performed Hi-C and RNA-sequencing assays at specific time points. Long-range intrachromosomal contacts increase over the course of differentiation, accompanied by widespread heteroeuchromatic compartment transitions that are tightly associated with transcription. Dynamic topologically associating domain boundaries strengthen and converge on an endothelial cell state, and function to regulate gene expression. Chromatin pairwise point interactions (DNA loops) increase in frequency during differentiation and are linked to the expression of genes essential to vascular biology. Chromatin dynamics guide transcription in endothelial cell development and promote the divergence of endothelial cells from cardiomyocytes.

Organization of gene programs revealed by unsupervised analysis of diverse gene-trait associations

Genome wide association studies provide statistical measures of gene–trait associations that reveal how genetic variation influences phenotypes. This study develops an unsupervised dimensionality reduction method called UnTANGLeD (Unsupervised Trait Analysis of Networks from Gene Level Data) which organizes 16,849 genes into discrete gene programs by measuring the statistical association between genetic variants and 1,393 diverse complex traits. UnTANGLeD reveals 173 gene clusters enriched for protein–protein interactions and highly distinct biological processes governing development, signalling, disease, and homeostasis. We identify diverse gene networks with robust interactions but not associated with known biological processes. Analysis of independent disease traits shows that UnTANGLeD gene clusters are conserved across all complex traits, providing a simple and powerful framework to predict novel gene candidates and programs influencing orthogonal disease phenotypes. Collectively, this study demonstrates that gene programs co-ordinately orchestrating cell functions can be identified without reliance on prior knowledge, providing a method for use in functional annotation, hypothesis generation, machine learning and prediction algorithms, and the interpretation of diverse genomic data.

Temporal perturbation of histone deacetylase activity reveals a requirement for HDAC1-3 in mesendoderm cell differentiation

Histone deacetylases (HDACs) are a class of enzymes that control chromatin state and influence cell fate. We evaluated the chromatin accessibility and transcriptome dynamics of zinc-containing HDACs during cell differentiation in vitro coupled with chemical perturbation to identify the role of HDACs in mesendoderm cell fate specification. Single-cell RNA sequencing analyses of HDAC expression during human pluripotent stem cell (hPSC) differentiation in vitro and mouse gastrulation in vivo reveal a unique association of HDAC1 and -3 with mesendoderm gene programs during exit from pluripotency. Functional perturbation with small molecules reveals that inhibition of HDAC1 and -3, but not HDAC2, induces mesoderm while impeding endoderm and early cardiac progenitor specification. These data identify unique biological functions of the structurally homologous enzymes HDAC1-3 in influencing hPSC differentiation from pluripotency toward mesendodermal and cardiac progenitor populations.

Cavin4 interacts with Bin1 to promote T-tubule formation and stability in developing skeletal muscle

The cavin proteins are essential for caveola biogenesis and function. Here, we identify a role for the muscle-specific component, Cavin4, in skeletal muscle T-tubule development by analyzing two vertebrate systems, mouse and zebrafish. In both models, Cavin4 localized to T-tubules, and loss of Cavin4 resulted in aberrant T-tubule maturation. In zebrafish, which possess duplicated cavin4 paralogs, Cavin4b was shown to directly interact with the T-tubule-associated BAR domain protein Bin1. Loss of both Cavin4a and Cavin4b caused aberrant accumulation of interconnected caveolae within the T-tubules, a fragmented T-tubule network enriched in Caveolin-3, and an impaired Ca2+ response upon mechanical stimulation. We propose a role for Cavin4 in remodeling the T-tubule membrane early in development by recycling caveolar components from the T-tubule to the sarcolemma. This generates a stable T-tubule domain lacking caveolae that is essential for T-tubule function.

Integrating single-cell genomics pipelines to discover mechanisms of stem cell differentiation

Pluripotent stem cells underpin a growing sector that leverages their differentiation potential for research, industry, and clinical applications. This review evaluates the landscape of methods in single-cell transcriptomics that are enabling accelerated discovery in stem cell science. We focus on strategies for scaling stem cell differentiation through multiplexed single-cell analyses, for evaluating molecular regulation of cell differentiation using new analysis algorithms, and methods for integration and projection analysis to classify and benchmark stem cell derivatives against in vivo cell types. By discussing the available methods, comparing their strengths, and illustrating strategies for developing integrated analysis pipelines, we provide user considerations to inform their implementation and interpretation.

scGPS: Determining Cell States and Global Fate Potential of Subpopulations

Finding cell states and their transcriptional relatedness is a main outcome from analysing single-cell data. In developmental biology, determining whether cells are related in a differentiation lineage remains a major challenge. A seamless analysis pipeline from cell clustering to estimating the probability of transitions between cell clusters is lacking. Here, we present Single Cell Global fate Potential of Subpopulations (scGPS) to characterise transcriptional relationship between cell states. scGPS decomposes mixed cell populations in one or more samples into clusters (SCORE algorithm) and estimates pairwise transitioning potential (scGPS algorithm) of any pair of clusters. SCORE allows for the assessment and selection of stable clustering results, a major challenge in clustering analysis. scGPS implements a novel approach, with machine learning classification, to flexibly construct trajectory connections between clusters. scGPS also has a feature selection functionality by network and modelling approaches to find biological processes and driver genes that connect cell populations. We applied scGPS in diverse developmental contexts and show superior results compared to a range of clustering and trajectory analysis methods. scGPS is able to identify the dynamics of cellular plasticity in a user-friendly workflow, that is fast and memory efficient. scGPS is implemented in R with optimised functions using C++ and is publicly available in Bioconductor.

Therapeutic Inhibition of Acid Sensing Ion Channel 1a Recovers Heart Function After Ischemia-Reperfusion Injury

Background: Ischemia-reperfusion injury (IRI) is one of the major risk factors implicated in morbidity and mortality associated with cardiovascular disease. During cardiac ischemia, the build-up of acidic metabolites results in decreased intracellular and extracellular pH that can reach as low as 6.0-6.5. The resulting tissue acidosis exacerbates ischemic injury and significantly impacts cardiac function. Methods: We used genetic and pharmacological methods to investigate the role of acid sensing ion channel 1a (ASIC1a) in cardiac IRI at the cellular and whole organ level. Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) as well as ex vivo and in vivo models of IRI were used to test the efficacy of ASIC1a inhibitors as pre- and post-conditioning therapeutic agents. Results: Analysis of human complex trait genetics indicate that variants in the ASIC1 genetic locus are significantly associated with cardiac and cerebrovascular ischemic injuries. Using hiPSC-CMs in vitro and murine ex vivo heart models, we demonstrate that genetic ablation of ASIC1a improves cardiomyocyte viability after acute IRI. Therapeutic blockade of ASIC1a using specific and potent pharmacological inhibitors recapitulates this cardioprotective effect. We used an in vivo model of myocardial infarction (MI) and two models of ex vivo donor heart procurement and storage as clinical models to show that ASIC1a inhibition improves post-IRI cardiac viability. Use of ASIC1a inhibitors as pre- or post-conditioning agents provided equivalent cardioprotection to benchmark drugs, including the sodium-hydrogen exchange inhibitor zoniporide. At the cellular and whole organ level, we show that acute exposure to ASIC1a inhibitors has no impact on cardiac ion channels regulating baseline electromechanical coupling and physiological performance. Conclusions: Collectively, our data provide compelling evidence for a novel pharmacological strategy involving ASIC1a blockade as a cardioprotective therapy to improve the viability of hearts subjected to IRI.

Elp2 mutations perturb the epitranscriptome and lead to a complex neurodevelopmental phenotype

Intellectual disability (ID) and autism spectrum disorder (ASD) are the most common neurodevelopmental disorders and are characterized by substantial impairment in intellectual and adaptive functioning, with their genetic and molecular basis remaining largely unknown. Here, we identify biallelic variants in the gene encoding one of the Elongator complex subunits, ELP2, in patients with ID and ASD. Modelling the variants in mice recapitulates the patient features, with brain imaging and tractography analysis revealing microcephaly, loss of white matter tract integrity and an aberrant functional connectome. We show that the Elp2 mutations negatively impact the activity of the complex and its function in translation via tRNA modification. Further, we elucidate that the mutations perturb protein homeostasis leading to impaired neurogenesis, myelin loss and neurodegeneration. Collectively, our data demonstrate an unexpected role for tRNA modification in the pathogenesis of monogenic ID and ASD and define Elp2 as a key regulator of brain development.

Single cell eQTL analysis identifies cell type-specific genetic control of gene expression in fibroblasts and reprogrammed induced pluripotent stem cells

BACKGROUND

The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution. By integrating single cell RNA-sequencing (scRNA-seq) and population genetics, we apply a framework in which to evaluate cell type-specific effects of genetic variation on gene expression.

RESULTS

Here, we perform scRNA-seq on 64,018 fibroblasts from 79 donors and map expression quantitative trait loci (eQTLs) at the level of individual cell types. We demonstrate that the majority of eQTLs detected in fibroblasts are specific to an individual cell subtype. To address if the allelic effects on gene expression are maintained following cell reprogramming, we generate scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type-specific eQTLs in iPSCs and show that the eQTLs in fibroblasts almost entirely disappear during reprogramming.

CONCLUSIONS

This work provides an atlas of how genetic variation influences gene expression across cell subtypes and provides evidence for patterns of genetic architecture that lead to cell type-specific eQTL effects.

Conserved Epigenetic Regulatory Logic Infers Genes Governing Cell Identity

Determining genes that orchestrate cell differentiation in development and disease remains a fundamental goal of cell biology. This study establishes a genome-wide metric based on the gene-repressive trimethylation of histone H3 at lysine 27 (H3K27me3) across hundreds of diverse cell types to identify genetic regulators of cell differentiation. We introduce a computational method, TRIAGE, which uses discordance between gene-repressive tendency and expression to identify genetic drivers of cell identity. We apply TRIAGE to millions of genome-wide single-cell transcriptomes, diverse omics platforms, and eukaryotic cells and tissue types. Using a wide range of data, we validate the performance of TRIAGE in identifying cell-type-specific regulatory factors across diverse species including human, mouse, boar, bird, fish, and tunicate. Using CRISPR gene editing, we use TRIAGE to experimentally validate RNF220 as a regulator of Ciona cardiopharyngeal development and SIX3 as required for differentiation of endoderm in human pluripotent stem cells. A record of this paper’s transparent peer review process is included in the Supplemental Information.

Perturbing genes controlling cell decisions have major implications in development or disease. However, identifying key regulatory genes from the thousands expressed in a cell is challenging. TRIAGE is a computational method that distills patterns of epigenetic repression across diverse cell types to infer regulatory genes using input gene expression data from any cell type. Demonstrating its utility, we combine single-cell RNA-seq and TRIAGE to identify and experimentally confirm novel regulators of heart development in evolutionarily distant species.

Genotype-free demultiplexing of pooled single-cell RNA-seq

A variety of methods have been developed to demultiplex pooled samples in a single cell RNA sequencing (scRNA-seq) experiment which either require hashtag barcodes or sample genotypes prior to pooling. We introduce scSplit which utilizes genetic differences inferred from scRNA-seq data alone to demultiplex pooled samples. scSplit also enables mapping clusters to original samples. Using simulated, merged, and pooled multi-individual datasets, we show that scSplit prediction is highly concordant with demuxlet predictions and is highly consistent with the known truth in cell-hashing dataset. scSplit is ideally suited to samples without external genotype information and is available at: https://github.com/jon-xu/scSplit.

Genetic Lineage Tracing of Sca-1+ Cells Reveals Endothelial but Not Myogenic Contribution to the Murine Heart

BACKGROUND

The adult mammalian heart displays a cardiomyocyte turnover rate of ≈1%/y throughout postnatal life and after injuries such as myocardial infarction (MI), but the question of which cell types drive this low level of new cardiomyocyte formation remains contentious. Cardiac-resident stem cells marked by stem cell antigen-1 (Sca-1, gene name Ly6a) have been proposed as an important source of cardiomyocyte renewal. However, the in vivo contribution of endogenous Sca-1⁺ cells to the heart at baseline or after MI has not been investigated.

METHODS

Here we generated Ly6a gene-targeted mice containing either a constitutive or an inducible Cre recombinase to perform genetic lineage tracing of Sca-1⁺ cells in vivo.

RESULTS

We observed that the contribution of endogenous Sca-1⁺ cells to the cardiomyocyte population in the heart was <0.005% throughout all of cardiac development, with aging, or after MI. In contrast, Sca-1⁺ cells abundantly contributed to the cardiac vasculature in mice during physiological growth and in the post-MI heart during cardiac remodeling. Specifically, Sca-1 lineage-traced endothelial cells expanded postnatally in the mouse heart after birth and into adulthood. Moreover, pulse labeling of Sca-1⁺ cells with an inducible Ly6a-MerCreMer allele also revealed a preferential expansion of Sca-1 lineage-traced endothelial cells after MI injury in the mouse.

CONCLUSIONS

Cardiac-resident Sca-1⁺ cells are not significant contributors to cardiomyocyte renewal in vivo. However, cardiac Sca-1⁺ cells represent a subset of vascular endothelial cells that expand postnatally with enhanced responsiveness to pathological stress in vivo.

Single-Cell Transcriptomic Analysis of Cardiac Differentiation from Human PSCs Reveals HOPX-Dependent Cardiomyocyte Maturation

Cardiac differentiation of human pluripotent stem cells (hPSCs) requires orchestration of dynamic gene regulatory networks during stepwise fate transitions but often generates immature cell types that do not fully recapitulate properties of their adult counterparts, suggesting incomplete activation of key transcriptional networks. We performed extensive single-cell transcriptomic analyses to map fate choices and gene expression programs during cardiac differentiation of hPSCs and identified strategies to improve in vitro cardiomyocyte differentiation. Utilizing genetic gain- and loss-of-function approaches, we found that hypertrophic signaling is not effectively activated during monolayer-based cardiac differentiation, thereby preventing expression of HOPX and its activation of downstream genes that govern late stages of cardiomyocyte maturation. This study therefore provides a key transcriptional roadmap of in vitro cardiac differentiation at single-cell resolution, revealing fundamental mechanisms underlying heart development and differentiation of hPSC-derived cardiomyocytes.

Single-cell RNA-seq of human induced pluripotent stem cells reveals cellular heterogeneity and cell state transitions between subpopulations

Heterogeneity of cell states represented in pluripotent cultures has not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC-CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method and, through this, identified four subpopulations distinguishable on the basis of their pluripotent state, including a core pluripotent population (48.3%), proliferative (47.8%), early primed for differentiation (2.8%), and late primed for differentiation (1.1%). For each subpopulation, we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets composed of 165 unique genes that denote the specific pluripotency states; using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to threefold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations and support our conclusions with results from two orthogonal pseudotime trajectory methods.

Chromatin and Transcriptional Analysis of Mesoderm Progenitor Cells Identifies HOPX as a Regulator of Primitive Hematopoiesis

We analyzed chromatin dynamics and transcriptional activity of human embryonic stem cell (hESC)-derived cardiac progenitor cells (CPCs) and KDR⁺/CD34⁺ endothelial cells generated from different mesodermal origins. Using an unbiased algorithm to hierarchically rank genes modulated at the level of chromatin and transcription, we identified candidate regulators of mesodermal lineage determination. HOPX, a non-DNA-binding homeodomain protein, was identified as a candidate regulator of blood-forming endothelial cells. Using HOPX reporter and knockout hESCs, we show that HOPX regulates blood formation. Loss of HOPX does not impact endothelial fate specification but markedly reduces primitive hematopoiesis, acting at least in part through failure to suppress Wnt/β-catenin signaling. Thus, chromatin state analysis permits identification of regulators of mesodermal specification, including a conserved role for HOPX in governing primitive hematopoiesis.

Inducible CRISPR genome editing platform in naive human embryonic stem cells reveals JARID2 function in self-renewal

To easily edit the genome of naïve human embryonic stem cells (hESC), we introduced a dual cassette encoding an inducible Cas9 into the AAVS1 site of naïve hESC (iCas9). The iCas9 line retained karyotypic stability, expression of pluripotency markers, differentiation potential, and stability in 5iLA and EPS pluripotency conditions. The iCas9 line induced efficient homology-directed repair (HDR) and non-homologous end joining (NHEJ) based mutations through CRISPR-Cas9 system. We utilized the iCas9 line to study the epigenetic regulator, PRC2 in early human pluripotency. The PRC2 requirement distinguishes between early pluripotency stages, however, what regulates PRC2 activity in these stages is not understood. We show reduced H3K27me3 and pluripotency markers in JARID2 2iL-I-F hESC mutants, indicating JARID2 requirement in maintenance of hESC 2iL-I-F state. These data suggest that JARID2 regulates PRC2 in 2iL-I-F state and the lack of PRC2 function in 5iLA state may be due to lack of sufficient JARID2 protein.

Subtle modifications to oxytocin produce ligands that retain potency and improved selectivity across species

Oxytocin and vasopressin mediate various physiological functions that are important for osmoregulation, reproduction, cardiovascular function, social behavior, memory, and learning through four G protein-coupled receptors that are also implicated in high-profile disorders. Targeting these receptors is challenging because of the difficulty in obtaining ligands that retain selectivity across rodents and humans for translational studies. We identified a selective and more stable oxytocin receptor (OTR) agonist by subtly modifying the pharmacophore framework of human oxytocin and vasopressin. [Se-Se]-oxytocin-OH displayed similar potency to oxytocin but improved selectivity for OTR, an effect that was retained in mice. Centrally infused [Se-Se]-oxytocin-OH potently reversed social fear in mice, confirming that this action was mediated by OTR and not by V1a or V1b vasopressin receptors. In addition, [Se-Se]-oxytocin-OH produced a more regular contraction pattern than did oxytocin in a preclinical labor induction and augmentation model using myometrial strips from cesarean sections. [Se-Se]-oxytocin-OH had no activity in human cardiomyocytes, indicating a potentially improved safety profile and therapeutic window compared to those of clinically used oxytocin. In conclusion, [Se-Se]-oxytocin-OH is a novel probe for validating OTR as a therapeutic target in various biological systems and is a promising new lead for therapeutic development. Our medicinal chemistry approach may also be applicable to other peptidergic signaling systems with similar selectivity issues.

Generating high-purity cardiac and endothelial derivatives from patterned mesoderm using human pluripotent stem cells

Human pluripotent stem cells (hPSCs) provide a valuable model for the study of human development and a means to generate a scalable source of cells for therapeutic applications. This protocol specifies cell fate efficiently into cardiac and endothelial lineages from hPSCs. The protocol takes 2 weeks to complete and requires experience in hPSC culture and differentiation techniques. Building on lessons taken from early development, this monolayer-directed differentiation protocol uses different concentrations of activin A and bone morphogenetic protein 4 (BMP4) to polarize cells into mesodermal subtypes that reflect mid-primitive-streak cardiogenic mesoderm and posterior-primitive-streak hemogenic mesoderm. This differentiation platform provides a basis for generating distinct cardiovascular progenitor populations that enable the derivation of cardiomyocytes and functionally distinct endothelial cell (EC) subtypes from cardiogenic versus hemogenic mesoderm with high efficiency without cell sorting. ECs derived from cardiogenic and hemogenic mesoderm can be matured into >90% CD31⁺/VE-cadherin⁺ definitive ECs. To test the functionality of ECs at different stages of differentiation, we provide methods for assaying the blood-forming potential and de novo lumen-forming activity of ECs. To our knowledge, this is the first protocol that provides a common platform for directed differentiation of cardiomyocytes and endothelial subtypes from hPSCs. This protocol yields endothelial differentiation efficiencies exceeding those of previously published protocols. Derivation of these cell types is a critical step toward understanding the basis of disease and generating cells with therapeutic potential.

The developmental origins and lineage contributions of endocardial endothelium

Endocardial development involves a complex orchestration of cell fate decisions that coordinate with endoderm formation and other mesodermal cell lineages. Historically, investigations into the contribution of endocardium in the developing embryo was constrained to the heart where these cells give rise to the inner lining of the myocardium and are a major contributor to valve formation. In recent years, studies have continued to elucidate the complexities of endocardial fate commitment revealing a much broader scope of lineage potential from developing endocardium. These studies cover a wide range of species and model systems and show direct contribution or fate potential of endocardium giving rise to cardiac vasculature, blood, fibroblast, and cardiomyocyte lineages. This review focuses on the marked expansion of knowledge in the area of endocardial fate potential. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Inhibition of β-catenin signaling respecifies anterior-like endothelium into beating human cardiomyocytes

During vertebrate development, mesodermal fate choices are regulated by interactions between morphogens such as activin/nodal, BMPs and Wnt/β-catenin that define anterior-posterior patterning and specify downstream derivatives including cardiomyocyte, endothelial and hematopoietic cells. We used human embryonic stem cells to explore how these pathways control mesodermal fate choices in vitro. Varying doses of activin A and BMP4 to mimic cytokine gradient polarization in the anterior-posterior axis of the embryo led to differential activity of Wnt/β-catenin signaling and specified distinct anterior-like (high activin/low BMP) and posterior-like (low activin/high BMP) mesodermal populations. Cardiogenic mesoderm was generated under conditions specifying anterior-like mesoderm, whereas blood-forming endothelium was generated from posterior-like mesoderm, and vessel-forming CD31⁺ endothelial cells were generated from all mesoderm origins. Surprisingly, inhibition of β-catenin signaling led to the highly efficient respecification of anterior-like endothelium into beating cardiomyocytes. Cardiac respecification was not observed in posterior-derived endothelial cells. Thus, activin/BMP gradients specify distinct mesodermal subpopulations that generate cell derivatives with unique angiogenic, hemogenic and cardiogenic properties that should be useful for understanding embryogenesis and developing therapeutics.

Summary: The manipulation of signals that control embryonic patterning allows human pluripotent stem cells to be differentiated into endothelial subpopulations with distinct haematopoietic, angiogenic and cardiogenic potential.

Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development

Background

The pervasive expression of circular RNA is a recently discovered feature of gene expression in highly diverged eukaryotes, but the functions of most circular RNAs are still unknown. Computational methods to discover and quantify circular RNA are essential. Moreover, discovering biological contexts where circular RNAs are regulated will shed light on potential functional roles they may play.

Results

We present a new algorithm that increases the sensitivity and specificity of circular RNA detection by discovering and quantifying circular and linear RNA splicing events at both annotated and un-annotated exon boundaries, including intergenic regions of the genome, with high statistical confidence. Unlike approaches that rely on read count and exon homology to determine confidence in prediction of circular RNA expression, our algorithm uses a statistical approach. Using our algorithm, we unveiled striking induction of general and tissue-specific circular RNAs, including in the heart and lung, during human fetal development. We discover regions of the human fetal brain, such as the frontal cortex, with marked enrichment for genes where circular RNA isoforms are dominant.

Conclusions

The vast majority of circular RNA production occurs at major spliceosome splice sites; however, we find the first examples of developmentally induced circular RNAs processed by the minor spliceosome, and an enriched propensity of minor spliceosome donors to splice into circular RNA at un-annotated, rather than annotated, exons. Together, these results suggest a potentially significant role for circular RNA in human development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0690-5) contains supplementary material, which is available to authorized users.

Cardiac development in zebrafish and human embryonic stem cells is inhibited by exposure to tobacco cigarettes and e-cigarettes

Background

Maternal smoking is a risk factor for low birth weight and other adverse developmental outcomes.

Objective

We sought to determine the impact of standard tobacco cigarettes and e-cigarettes on heart development in vitro and in vivo.

Methods

Zebrafish (Danio rerio) were used to assess developmental effects in vivo and cardiac differentiation of human embryonic stem cells (hESCs) was used as a model for in vitro cardiac development.

Results

In zebrafish, exposure to both types of cigarettes results in broad, dose-dependent developmental defects coupled with severe heart malformation, pericardial edema and reduced heart function. Tobacco cigarettes are more toxic than e-cigarettes at comparable nicotine concentrations. During cardiac differentiation of hESCs, tobacco smoke exposure results in a delayed transition through mesoderm. Both types of cigarettes decrease expression of cardiac transcription factors in cardiac progenitor cells, suggesting a persistent delay in differentiation. In definitive human cardiomyocytes, both e-cigarette- and tobacco cigarette-treated samples showed reduced expression of sarcomeric genes such as MLC2v and MYL6. Furthermore, tobacco cigarette-treated samples had delayed onset of beating and showed low levels and aberrant localization of N-cadherin, reduced myofilament content with significantly reduced sarcomere length, and increased expression of the immature cardiac marker smooth muscle alpha-actin.

Conclusion

These data indicate a negative effect of both tobacco cigarettes and e-cigarettes on heart development in vitro and in vivo. Tobacco cigarettes are more toxic than E-cigarettes and exhibit a broader spectrum of cardiac developmental defects.

Methods for assessing the electromechanical integration of human pluripotent stem cell-derived cardiomyocyte grafts

Cardiomyocytes derived from human pluripotent stem cells show tremendous promise for the replacement of myocardium and contractile function lost to infarction. However, until recently, no methods were available to directly determine whether these stem cell-derived grafts actually couple with host myocardium and fire synchronously following transplantation in either intact or injured hearts. To resolve this uncertainty, our group has developed techniques for the intravital imaging of hearts engrafted with stem cell-derived cardiomyocytes that have been modified to express the genetically encoded protein calcium sensor, GCaMP. When combined with the simultaneously recorded electrocardiogram, this protocol allows one to make quantitative assessments as to the presence and extent of host-graft electrical coupling as well as the timing and pattern of graft activation. As described here, this system has been employed to investigate the electromechanical integration of human embryonic stem cell-derived cardiomyocytes in a guinea pig model of cardiac injury, but analogous approaches should be applicable to other human graft cell types and animal models.

Transmembrane protein 88: a Wnt regulatory protein that specifies cardiomyocyte development

Genetic regulation of the cell fate transition from lateral plate mesoderm to the specification of cardiomyocytes requires suppression of Wnt/β-catenin signaling, but the mechanism for this is not well understood. By analyzing gene expression and chromatin dynamics during directed differentiation of human embryonic stem cells (hESCs), we identified a suppressor of Wnt/β-catenin signaling, transmembrane protein 88 (TMEM88), as a potential regulator of cardiovascular progenitor cell (CVP) specification. During the transition from mesoderm to the CVP, TMEM88 has a chromatin signature of genes that mediate cell fate decisions, and its expression is highly upregulated in advance of key cardiac transcription factors in vitro and in vivo. In early zebrafish embryos, tmem88a is expressed broadly in the lateral plate mesoderm, including the bilateral heart fields. Short hairpin RNA targeting of TMEM88 during hESC cardiac differentiation increases Wnt/β-catenin signaling, confirming its role as a suppressor of this pathway. TMEM88 knockdown has no effect on NKX2.5 or GATA4 expression, but 80% of genes most highly induced during CVP development have reduced expression, suggesting adoption of a new cell fate. In support of this, analysis of later stage cell differentiation showed that TMEM88 knockdown inhibits cardiomyocyte differentiation and promotes endothelial differentiation. Taken together, TMEM88 is crucial for heart development and acts downstream of GATA factors in the pre-cardiac mesoderm to specify lineage commitment of cardiomyocyte development through inhibition of Wnt/β-catenin signaling.

Diastolic dysfunction and thin filament dysregulation resulting from excitation-contraction uncoupling in a mouse model of restrictive cardiomyopathy

Restrictive cardiomyopathy (RCM) has been linked to mutations in the thin filament regulatory protein cardiac troponin I (cTnI). As the pathogenesis of RCM from genotype to clinical phenotype is not fully understood, transgenic (Tg) mice were generated with cardiac specific expression of an RCM-linked missense mutation (R193H) in cTnI. R193H Tg mouse hearts with 15% stoichiometric replacement had smaller hearts and significantly elevated end diastolic pressures (EDP) in vivo. Using a unique carbon microfiber-based assay, membrane intact R193H adult cardiac myocytes generated higher passive tensions across a range of physiologic sarcomere lengths resulting in significant Ca(2+) independent cellular diastolic tone that was manifest in vivo as elevated organ-level EDP. Sarcomere relaxation and Ca(2+) decay was uncoupled in isolated R193H Tg adult myocytes due to the increase in myofilament Ca(2+) sensitivity of tension, decreased passive compliance of the sarcomere, and adaptive in vivo changes in which phospholamban (PLN) content was decreased. Further evidence of Ca(2+) and mechanical uncoupling in R193H Tg myocytes was demonstrated by the biphasic response of relaxation to increased pacing frequency versus the negative staircase seen with Ca(2+) decay. In comparison, non-transgenic myocyte relaxation closely paralleled the accelerated Ca(2+) decay. Ca(2+) transient amplitude was also significantly blunted in R193H Tg myocytes despite normal mechanical shortening resulting in myocyte hypercontractility when compared to non-transgenics. These results identify for the first time that a single point mutation in cTnI, R193H, directly causes elevated EDP due to a myocyte intrinsic loss of compliance independent of Ca(2+) cycling or altered cardiac morphology. The compound influence of impaired relaxation and elevated EDP represents a clinically severe form of diastolic dysfunction similar to the hemodynamic state documented in RCM patients.

Cardiac disease in mucopolysaccharidosis type I attributed to catecholaminergic and hemodynamic deficiencies

Cardiac dysfunction is a common cause of death among pediatric patients with mutations in the lysosomal hydrolase α-l-iduronidase (IDUA) gene, which causes mucopolysaccharidosis type I (MPS-I). The purpose of this study was to analyze adrenergic regulation of cardiac hemodynamic function in MPS-I. An analysis of murine heart function was performed using conductance micromanometry to assess in vivo cardiac hemodynamics. Although MPS-I (IDUA(-/-)) mice were able to maintain normal cardiac output and ejection fraction at baseline, this cohort had significantly compromised systolic and diastolic function compared with IDUA(+/-) control mice. During dobutamine infusion MPS-I mice did not significantly increase cardiac output from baseline, indicative of blunted cardiac reserve. Autonomic tone, measured functionally by β-blockade, indicated that MPS-I mice required catecholaminergic stimulation to maintain baseline hemodynamics. Survival analysis showed mortality only among MPS-I mice. Linear regression analysis revealed that heightened end-systolic volume in the resting heart is significantly correlated with susceptibility to mortality in MPS-I hearts. This study reveals that cardiac remodeling in the pathology of MPS-I involves heightened adrenergic tone at the expense of cardiac reserve with cardiac decompensation predicted on the basis of increased baseline systolic volumes.

Influence of genetic background on ex vivo and in vivo cardiac function in several commonly used inbred mouse strains

Inbred mouse strains play a critical role in biomedical research. Genetic homogeneity within inbred strains and their general amenability to genetic manipulation have made them an ideal resource for dissecting the physiological function(s) of individual genes. However, the inbreeding that makes inbred mice so useful also results in genetic divergence between them. This genetic divergence is often unaccounted for but may be a confounding factor when comparing studies that have utilized distinct inbred strains. Here, we compared the cardiac function of C57BL/6J mice to seven other commonly used inbred mouse strains: FVB/NJ, DBA/2J, C3H/HeJ, BALB/cJ, 129X1/SvJ, C57BL/10SnJ, and 129S1/SvImJ. The assays used to compare cardiac function were the ex vivo isolated Langendorff heart preparation and in vivo real-time hemodynamic analysis using conductance micromanometry. We report significant strain-dependent differences in cardiac function between C57BL/6J and other commonly used inbred strains. C57BL/6J maintained better cardiac function than most inbred strains after ex vivo ischemia, particularly compared with 129S1/SvImJ, 129X1/SvJ, and C57BL/10SnJ strains. However, during in vivo acute hypoxia 129X1/SvJ and 129S1/SvImJ maintained relatively normal cardiac function, whereas C57BL/6J animals showed dramatic cardiac decompensation. Additionally, C3H/HeJ showed rapid and marked cardiac decompensation in response to esmolol infusion compared with effects of other strains. These findings demonstrate the complex effects of genetic divergence between inbred strains on cardiac function. These results may help inform analysis of gene ablation or transgenic studies and further demonstrate specific quantitative traits that could be useful in discovery of genetic modifiers relevant to cardiac health and disease.

Aesthetic cardiology: adipose-derived stem cells for myocardial repair

Stem cell biology has increasingly gained scientific and public interest in recent years. In particular, the use of stem cells for treatment of heart disease has been strongly pursued within the scientific and medical communities. Significant effort has gone into the use of adult tissue-derived stem cells for cardiac repair including bone marrow, blood, and cardiac-derived cell populations. Significant interest in this area has been balanced by the difficulties of understanding stem cells, cardiac injury, and the amalgamation of these areas of investigation in translational medicine. Recent studies have emerged on adipose-derived stem cells which show the potential for cardiac lineage development in vitro and may have application in cell-mediated in vivo therapy for the diseased heart. This review provides a summary of current findings within the field of adipose-derived stem cell biology regarding their cardiac differentiation potential.

Pathogenic peptide deviations support a model of adaptive evolution of chordate cardiac performance by troponin mutations

In cardiac muscle, the troponin (cTn) complex is a key regulator of myofilament calcium sensitivity because it serves as a molecular switch required for translating myocyte calcium fluxes into sarcomeric contraction and relaxation. Studies of several species suggest that ectotherm chordates have myofilaments with heightened calcium responsiveness. However, genetic polymorphisms in cTn that cause increased myofilament sensitivity to activating calcium in mammals result in cardiac disease including arrhythmias, diastolic dysfunction, and increased susceptibility to sudden cardiac death. We hypothesized that specific residue modifications in the regulatory arm of troponin I (TnI) were critical in mediating the observed decrease in myofilament calcium sensitivity within the mammalian taxa. We performed large-scale phylogenetic analysis, atomic resolution molecular dynamics simulations and modeling, and computational alanine scanning. This study provides evidence that a His to Ala substitution within mammalian cardiac TnI (cTnI) reduced the thermodynamic potential at the interface between cTnI and cardiac TnC (cTnC) in the calcium-saturated state by disrupting a strong intermolecular electrostatic interaction. This key residue modification reduced myofilament calcium sensitivity by making cTnI molecularly untethered from cTnC. To meet the requirements for refined mammalian adult cardiac performance, we propose that compensatory evolutionary pressures favored mutations that enhanced the relaxation properties of cTn by decreasing its sensitivity to activating calcium.

Artificial selection for whole animal low intrinsic aerobic capacity co-segregates with hypoxia-induced cardiac pump failure

Oxygen metabolism is a strong predictor of the general health and fitness of an organism. In this study, we hypothesized that a divergence in intrinsic aerobic fitness would co-segregate with susceptibility for cardiovascular dysfunction. To test this hypothesis, cardiac function was assessed in rats specifically selected over nineteen generations for their low (LCR) and high (HCR) intrinsic aerobic running capacity. As an integrative marker of native aerobic capacity, run time to exhaustion between LCR and HCR rats had markedly diverged by 436% at generation nineteen of artificial selection. In vivo assessment of baseline cardiac function by echocardiography and catheter-based conductance micromanometry showed no marked difference in cardiac performance. However, when challenged by exposure to acute hypoxia, cardiac pump failure occurred significantly earlier in LCR rats compared to HCR animals. Acute cardiac decompensation in LCR rats was exclusively due to the development of intractable irregular ventricular contractions. Analysis of isolated cardiac myocytes showed significantly slower sarcomeric relaxation and delayed kinetics of calcium cycling in LCR myocytes compared to HCR myocytes. This study also revealed that artificial selection for low native aerobic capacity is a novel pathologic stimulus that results in myosin heavy chain isoform switching in the heart as shown by increased levels of beta-MHC in LCR rats. Together, these results provide evidence that alterations in sub-cellular calcium handling and MHC isoform composition are associated with susceptibility to compensatory cardiac remodeling and hypoxia induced pump failure in animals with low intrinsic aerobic capacity.

Single histidine button in cardiac troponin I sustains heart performance in response to severe hypercapnic respiratory acidosis in vivo

Intracellular acidosis is a profound negative regulator of myocardial performance. We hypothesized that titrating myofilament calcium sensitivity by a single histidine substituted cardiac troponin I (A164H) would protect the whole animal physiological response to acidosis in vivo. To experimentally induce severe hypercapnic acidosis, mice were exposed to a 40% CO(2) challenge. By echocardiography, it was found that systolic function and ventricular geometry were maintained in cTnI A164H transgenic (Tg) mice. By contrast, non-Tg (Ntg) littermates experienced rapid and marked cardiac decompensation during this same challenge. For detailed hemodymanic assessment, Millar pressure-conductance catheterization was performed while animals were treated with a beta-blocker, esmolol, during a severe hypercapnic acidosis challenge. Survival and load-independent measures of contractility were significantly greater in Tg vs. Ntg mice. This assay showed that Ntg mice had 100% mortality within 5 min of acidosis. By contrast, systolic and diastolic function were protected in Tg mice during acidosis, and they had 100% survival. This study shows that, independent of any beta-adrenergic compensation, myofilament-based molecular manipulation of inotropy by histidine-modified troponin I maintains cardiac inotropic and lusitropic performance and markedly improves survival during severe acidosis in vivo.

Single histidine-substituted cardiac troponin I confers protection from age-related systolic and diastolic dysfunction

AIMS

Contractile dysfunction associated with myocardial ischaemia is a significant cause of morbidity and mortality in the elderly. Strategies to protect the aged heart from ischaemia-mediated pump failure are needed. We hypothesized that troponin I-mediated augmentation of myofilament calcium sensitivity would protect cardiac function in aged mice.

METHODS AND RESULTS

To address this, we investigated transgenic (Tg) mice expressing a histidine-substituted form of adult cardiac troponin I (cTnI A164H), which increases myofilament calcium sensitivity in a pH-dependent manner. Serial echocardiography revealed that Tg hearts showed significantly improved systolic function at 4 months, which was sustained for 2 years based on ejection fraction and velocity of circumferential fibre shortening. Age-related diastolic dysfunction was also attenuated in Tg mice as assessed by Doppler measurements of the mitral valve inflow and lateral annulus Doppler tissue imaging. During acute hypoxia, cardiac contractility significantly improved in aged Tg mice made evident by increased stroke volume, end systolic pressure, and +dP/dt compared with non-transgenic mice.

CONCLUSION

This study shows that increasing myofilament function by means of a pH-responsive histidine button engineered into cTnI results in enhanced baseline heart function in Tg mice over their lifetime, and during acute hypoxia improves survival in aged mice by maintaining cardiac contractility.

Non-canonical Wnt signaling enhances differentiation of Sca1+/c-kit+ adipose-derived murine stromal vascular cells into spontaneously beating cardiac myocytes

Recent reports have described a stem cell population termed stromal vascular cells (SVCs) derived from the stromal vascular fraction of adipose tissue, which are capable of intrinsic differentiation into spontaneously beating cardiomyocytes in vitro. The objective of this study was to further define the cardiac lineage differentiation potential of SVCs in vitro and to establish methods for enriching SVC-derived beating cardiac myocytes. SVCs were isolated from the stromal vascular fraction of murine adipose tissue. Cells were cultured in methylcellulose-based murine stem cell media. Analysis of SVC-derived beating myocytes included Western blot and calcium imaging. Enrichment of acutely isolated SVCs was carried out using antibody-tagged magnetic nanoparticles, and pharmacologic manipulation of Wnt and cytokine signaling. Under initial media conditions, spontaneously beating SVCs expressed both cardiac developmental and adult protein isoforms. Functionally, this specialized population can spontaneously contract and pace under field stimulation and shows the presence of coordinated calcium transients. Importantly, this study provides evidence for two independent mechanisms of enriching the cardiac differentiation of SVCs. First, this study shows that differentiation of SVCs into cardiac myocytes is augmented by non-canonical Wnt agonists, canonical Wnt antagonists, and cytokines. Second, SVCs capable of cardiac lineage differentiation can be enriched by selection for stem cell-specific membrane markers Sca1 and c-kit. Adipose-derived SVCs are a unique population of stem cells that show evidence of cardiac lineage development making them a potential source for stem cell-based cardiac regeneration studies.