Single-cell multi-modal integrative analyses highlight functional dynamic gene regulatory networks directing human cardiac development.

Illuminating the precise stepwise genetic programs directing cardiac development provides insights into the mechanisms of congenital heart disease and strategies for cardiac regenerative therapies. Here, we integrate in vitro and in vivo human single-cell multi-omic studies with high-throughput functional genomic screening to reveal dynamic, cardiac-specific gene regulatory networks (GRNs) and transcriptional regulators during human cardiomyocyte development. Interrogating developmental trajectories reconstructed from single-cell data unexpectedly reveal divergent cardiomyocyte lineages with distinct gene programs based on developmental signaling pathways. High-throughput functional genomic screens identify key transcription factors from inferred GRNs that are functionally relevant for cardiomyocyte lineages derived from each pathway. Notably, we discover a critical heat shock transcription factor 1 (HSF1)-mediated cardiometabolic GRN controlling cardiac mitochondrial/metabolic function and cell survival, also observed in fetal human cardiomyocytes. Overall, these multi-modal genomic studies enable the systematic discovery and validation of coordinated GRNs and transcriptional regulators controlling the development of distinct human cardiomyocyte populations.

Spatially organized cellular communities form the developing human heart.

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

Cell Layers: uncovering clustering structure in unsupervised single-cell transcriptomic analysis.

Motivation: Unsupervised clustering of single-cell transcriptomics is a powerful method for identifying cell populations. Static visualization techniques for single-cell clustering only display results for a single resolution parameter. Analysts will often evaluate more than one resolution parameter but then only report one. Results: We developed Cell Layers, an interactive Sankey tool for the quantitative investigation of gene expression, co-expression, biological processes and cluster integrity across clustering resolutions. Cell Layers enhances the interpretability of single-cell clustering by linking molecular data and cluster evaluation metrics, providing novel insight into cell populations. Availability and implementation: https://github.com/apblair/CellLayers.

Human-gained heart enhancers are associated with species-specific cardiac attributes.

The heart, a vital organ which is first to develop, has adapted its size, structure and function in order to accommodate the circulatory demands for a broad range of animals. Although heart development is controlled by a relatively conserved network of transcriptional/chromatin regulators, how the human heart has evolved species-specific features to maintain adequate cardiac output and function remains to be defined. Here, we show through comparative epigenomic analysis the identification of enhancers and promoters that have gained activity in humans during cardiogenesis. These cis-regulatory elements (CREs) are associated with genes involved in heart development and function, and may account for species-specific differences between human and mouse hearts. Supporting these findings, genetic variants that are associated with human cardiac phenotypic/disease traits, particularly those differing between human and mouse, are enriched in human-gained CREs. During early stages of human cardiogenesis, these CREs are also gained within genomic loci of transcriptional regulators, potentially expanding their role in human heart development. In particular, we discovered that gained enhancers in the locus of the early human developmental regulator ZIC3 are selectively accessible within a subpopulation of mesoderm cells which exhibits cardiogenic potential, thus possibly extending the function of ZIC3 beyond its conserved left-right asymmetry role. Genetic deletion of these enhancers identified a human gained enhancer that was required for not only ZIC3 and early cardiac gene expression at the mesoderm stage but also cardiomyocyte differentiation. Overall, our results illuminate how human gained CREs may contribute to human-specific cardiac attributes, and provide insight into how transcriptional regulators may gain cardiac developmental roles through the evolutionary acquisition of enhancers.

Cardiac cell type-specific gene regulatory programs and disease risk association.

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.

Canonical Wnt5b Signaling Directs Outlying Nkx2.5+ Mesoderm into Pacemaker Cardiomyocytes.

Pacemaker cardiomyocytes that create the sinoatrial node are essential for the initiation and maintenance of proper heart rhythm. However, illuminating developmental cues that direct their differentiation has remained particularly challenging due to the unclear cellular origins of these specialized cardiomyocytes. By discovering the origins of pacemaker cardiomyocytes, we reveal an evolutionarily conserved Wnt signaling mechanism that coordinates gene regulatory changes directing mesoderm cell fate decisions, which lead to the differentiation of pacemaker cardiomyocytes. We show that in zebrafish, pacemaker cardiomyocytes derive from a subset of Nkx2.5+ mesoderm that responds to canonical Wnt5b signaling to initiate the cardiac pacemaker program, including activation of pacemaker cell differentiation transcription factors Isl1 and Tbx18 and silencing of Nkx2.5. Moreover, applying these developmental findings to human pluripotent stem cells (hPSCs) notably results in the creation of hPSC-pacemaker cardiomyocytes, which successfully pace three-dimensional bioprinted hPSC-cardiomyocytes, thus providing potential strategies for biological cardiac pacemaker therapy.

Transcriptionally active HERV-H retrotransposons demarcate topologically associating domains in human pluripotent stem cells.

Chromatin architecture has been implicated in cell type-specific gene regulatory programs, yet how chromatin remodels during development remains to be fully elucidated. Here, by interrogating chromatin reorganization during human pluripotent stem cell (hPSC) differentiation, we discover a role for the primate-specific endogenous retrotransposon human endogenous retrovirus subfamily H (HERV-H) in creating topologically associating domains (TADs) in hPSCs. Deleting these HERV-H elements eliminates their corresponding TAD boundaries and reduces the transcription of upstream genes, while de novo insertion of HERV-H elements can introduce new TAD boundaries. The ability of HERV-H to create TAD boundaries depends on high transcription, as transcriptional repression of HERV-H elements prevents the formation of boundaries. This ability is not limited to hPSCs, as these actively transcribed HERV-H elements and their corresponding TAD boundaries also appear in pluripotent stem cells from other hominids but not in more distantly related species lacking HERV-H elements. Overall, our results provide direct evidence for retrotransposons in actively shaping cell type- and species-specific chromatin architecture.

Combinatorial interactions of genetic variants in human cardiomyopathy.

Dilated cardiomyopathy (DCM) is a leading cause of morbidity and mortality worldwide; yet how genetic variation and environmental factors impact DCM heritability remains unclear. Here, we report that compound genetic interactions between DNA sequence variants contribute to the complex heritability of DCM. By using genetic data from a large family with a history of DCM, we discovered that heterozygous sequence variants in the TROPOMYOSIN 1 (TPM1) and VINCULIN (VCL) genes cose-gregate in individuals affected by DCM. In vitro studies of patient-derived and isogenic human-pluripotent-stem-cell-derived cardio-myocytes that were genome-edited via CRISPR to create an allelic series of TPM1 and VCL variants revealed that cardiomyocytes with both TPM1 and VCL variants display reduced contractility and sarcomeres that are less organized. Analyses of mice genetically engineered to harbour these human TPM1 and VCL variants show that stress on the heart may also influence the variable penetrance and expressivity of DCM-associated genetic variants in vivo. We conclude that compound genetic variants can interact combinatorially to induce DCM, particularly when influenced by other disease-provoking stressors.

mRNAs and lncRNAs intrinsically form secondary structures with short end-to-end distances.

The 5' and 3' termini of RNA play important roles in many cellular processes. Using Förster resonance energy transfer (FRET), we show that mRNAs and lncRNAs have an intrinsic propensity to fold in the absence of proteins into structures in which the 5' end and 3' end are ≤7 nm apart irrespective of mRNA length. Computational estimates suggest that the inherent proximity of the ends is a universal property of most mRNA and lncRNA sequences. Only guanosine-depleted RNA sequences with low sequence complexity are unstructured and exhibit end-to-end distances expected for the random coil conformation of RNA. While the biological implications remain to be explored, short end-to-end distances could facilitate the binding of protein factors that regulate translation initiation by bridging mRNA 5' and 3' ends. Furthermore, our studies provide the basis for measuring, computing and manipulating end-to-end distances and secondary structure in RNA in research and biotechnology.

Cell-Surface Marker Signature for Enrichment of Ventricular Cardiomyocytes Derived from Human Embryonic Stem Cells.

To facilitate understanding of human cardiomyocyte (CM) subtype specification, and the study of ventricular CM biology in particular, we developed a broadly applicable strategy for enrichment of ventricular cardiomyocytes (VCMs) derived from human embryonic stem cells (hESCs). A bacterial artificial chromosome transgenic H9 hESC line in which GFP expression was driven by the human ventricular-specific myosin light chain 2 (MYL2) promoter was generated, and screened to identify cell-surface markers specific for MYL2-GFP-expressing VCMs. A CD77+/CD200- cell-surface signature facilitated isolation of >97% cardiac troponin I-positive cells from H9 hESC differentiation cultures, with 65% expressing MYL2-GFP. This study provides a tool for VCM enrichment when using some, but not all, human pluripotent stem cell lines. Tools generated in this study can be utilized toward understanding CM subtype specification, and enriching for VCMs for therapeutic applications.

EF-G Activation by Phosphate Analogs.

Elongation factor G (EF-G) is a universally conserved translational GTPase that promotes the translocation of tRNA and mRNA through the ribosome. EF-G binds to the ribosome in a GTP-bound form and subsequently catalyzes GTP hydrolysis. The contribution of the ribosome-stimulated GTP hydrolysis by EF-G to tRNA/mRNA translocation remains debated. Here, we show that while EF-G•GDP does not stably bind to the ribosome and induce translocation, EF-G•GDP in complex with phosphate group analogs BeF3(-) and AlF4(-) promotes the translocation of tRNA and mRNA. Furthermore, the rates of mRNA translocation induced by EF-G in the presence of GTP and a non-hydrolyzable analog of GTP, GDP•BeF3(-) are similar. Our results are consistent with the model suggesting that GTP hydrolysis is not directly coupled to mRNA/tRNA translocation. Hence, GTP binding is required to induce the activated, translocation-competent conformation of EF-G while GTP hydrolysis triggers EF-G release from the ribosome.

Movement of elongation factor G between compact and extended conformations.

Previous structural studies suggested that ribosomal translocation is accompanied by large interdomain rearrangements of elongation factor G (EF-G). Here, we follow the movement of domain IV of EF-G relative to domain II of EF-G using ensemble and single-molecule Förster resonance energy transfer. Our results indicate that ribosome-free EF-G predominantly adopts a compact conformation that can also, albeit infrequently, transition into a more extended conformation in which domain IV moves away from domain II. By contrast, ribosome-bound EF-G predominantly adopts an extended conformation regardless of whether it is interacting with pretranslocation ribosomes or with posttranslocation ribosomes. Our data suggest that ribosome-bound EF-G may also occasionally sample at least one more compact conformation. GTP hydrolysis catalyzed by EF-G does not affect the relative stability of the observed conformations in ribosome-free and ribosome-bound EF-G. Our data support a model suggesting that, upon binding to a pretranslocation ribosome, EF-G moves from a compact to a more extended conformation. This transition is not coupled to but likely precedes both GTP hydrolysis and mRNA/tRNA translocation.

Following movement of domain IV of elongation factor G during ribosomal translocation.

Translocation of mRNA and tRNAs through the ribosome is catalyzed by a universally conserved elongation factor (EF-G in prokaryotes and EF-2 in eukaryotes). Previous studies have suggested that ribosome-bound EF-G undergoes significant structural rearrangements. Here, we follow the movement of domain IV of EF-G, which is critical for the catalysis of translocation, relative to protein S12 of the small ribosomal subunit using single-molecule FRET. We show that ribosome-bound EF-G adopts distinct conformations corresponding to the pre- and posttranslocation states of the ribosome. Our results suggest that, upon ribosomal translocation, domain IV of EF-G moves toward the A site of the small ribosomal subunit and facilitates the movement of peptidyl-tRNA from the A to the P site. We found no evidence of direct coupling between the observed movement of domain IV of EF-G and GTP hydrolysis. In addition, our results suggest that the pretranslocation conformation of the EF-G-ribosome complex is significantly less stable than the posttranslocation conformation. Hence, the structural rearrangement of EF-G makes a considerable energetic contribution to promoting tRNA translocation.

Human biliverdin reductase-based peptides activate and inhibit glucose uptake through direct interaction with the kinase domain of insulin receptor.

Insulin binding changes conformation of the insulin receptor kinase (IRK) domain and initiates glucose uptake through the insulin, IGF-1, phosphatidyl inositol 3-kinase (PI3K), and MAPK pathways; human biliverdin reductase (hBVR) is an IRK substrate and pathway effector. This is the first report on hBVR peptide-mediated IRK activation and conformational change. (290)KYCCSRK, which increased IRK V(max) without changing K(m), stimulated glucose uptake and potentiated insulin and IGF-1 stimulation in 4 cell lines. KYCCSRK in native hBVR was necessary for the hBVR and IRK cross-activation. Peptide treatment also activated PI3K downstream effectors, Akt and ERK, phosphorylation, and Elk transcriptional activity. In cells transfected with CMV-regulated EGFP-VP-peptide plasmid, C(292)→A mutant did not stimulate glucose uptake; K(296)→A decreased uptake and kinase activity. KEDQYMKMTV, corresponding to hBVR's SH2-binding domain, was a potent inhibitor of glucose uptake and IRK. The mechanism of action of peptides was examined using cells expressing IRK (aa 988-1263) activated by coexpressed KYCCSRK. Three active cys-mutants of IRK, with fluorophore coupled to cysteines, C(1056), C(1138), or C(1234), were examined for changes in fluorescence emission spectra in the presence of peptides. KYCCSRK and KEDQYMKMTV bound to different sites in IRK. The findings identify novel agents for activating or inhibiting insulin signaling and offer a new approach for treatment of type 2 diabetes and hypoglycemia.