Complex Interdependence Regulates Heterotypic Transcription Factor Distribution and Coordinates Cardiogenesis

Context-dependent T-BOX transcription factor family: from biology to targeted therapy

T-BOX factors belong to an evolutionarily conserved family of transcription factors. T-BOX factors not only play key roles in growth and development but are also involved in immunity, cancer initiation, and progression. Moreover, the same T-BOX molecule exhibits different or even opposite effects in various developmental processes and tumor microenvironments. Understanding the multiple roles of context-dependent T-BOX factors in malignancies is vital for uncovering the potential of T-BOX-targeted cancer therapy. We summarize the physiological roles of T-BOX factors in different developmental processes and their pathological roles observed when their expression is dysregulated. We also discuss their regulatory roles in tumor immune microenvironment (TIME) and the newly arising questions that remain unresolved. This review will help in systematically and comprehensively understanding the vital role of the T-BOX transcription factor family in tumor physiology, pathology, and immunity. The intention is to provide valuable information to support the development of T-BOX-targeted therapy.

The Impact of Natriuretic Peptides on Heart Development, Homeostasis, and Disease

During mammalian heart development, the clustered genes encoding peptide hormones, Natriuretic Peptide A (NPPA; ANP) and B (NPPB; BNP), are transcriptionally co-regulated and co-expressed predominately in the atrial and ventricular trabecular cardiomyocytes. After birth, expression of NPPA and a natural antisense transcript NPPA-AS1 becomes restricted to the atrial cardiomyocytes. Both NPPA and NPPB are induced by cardiac stress and serve as markers for cardiovascular dysfunction or injury. NPPB gene products are extensively used as diagnostic and prognostic biomarkers for various cardiovascular disorders. Membrane-localized guanylyl cyclase receptors on many cell types throughout the body mediate the signaling of the natriuretic peptide ligands through the generation of intracellular cGMP, which interacts with and modulates the activity of cGMP-activated kinase and other enzymes and ion channels. The natriuretic peptide system plays a fundamental role in cardio-renal homeostasis, and its potent diuretic and vasodilatory effects provide compensatory mechanisms in cardiac pathophysiological conditions and heart failure. In addition, both peptides, but also CNP, have important intracardiac actions during heart development and homeostasis independent of the systemic functions. Exploration of the intracardiac functions may provide new leads for the therapeutic utility of natriuretic peptide-mediated signaling in heart diseases and rhythm disorders. Here, we review recent insights into the regulation of expression and intracardiac functions of NPPA and NPPB during heart development, homeostasis, and disease.

Molecular regulation of PPARγ/RXRα signaling by the novel cofactor ZFP407

Cofactors interacting with PPARγ can regulate adipogenesis and adipocyte metabolism by modulating the transcriptional activity and selectivity of PPARγ signaling. ZFP407 was previously demonstrated to regulate PPARγ target genes such as GLUT4, and its overexpression improved glucose homeostasis in mice. Here, using a series of molecular assays, including protein-interaction studies, mutagenesis, and ChIP-seq, ZFP407 was found to interact with the PPARγ/RXRα protein complex in the nucleus of adipocytes. Consistent with this observation, ZFP407 ChIP-seq peaks significantly overlapped with PPARγ ChIP-seq peaks, with more than half of ZFP407 peaks overlapping with PPARγ peaks. Transcription factor binding motifs enriched in these overlapping sites included CTCF, RARα/RXRγ, TP73, and ELK1, which regulate cellular development and function within adipocytes. Site-directed mutagenesis of frequent PPARγ phosphorylation or SUMOylation sites did not prevent its regulation by ZFP407, while mutagenesis of ZFP407 domains potentially necessary for RXR and PPARγ binding abrogated any impact of ZFP407 on PPARγ activity. These data suggest that ZFP407 controls the activity of PPARγ, but does so independently of post-translational modifications, likely by direct binding, establishing ZFP407 as a newly identified PPARγ cofactor. In addition, ZFP407 ChIP-seq analyses identified regions that did not overlap with PPARγ peaks. These non-overlapping peaks were significantly enriched for the transcription factor binding motifs of TBX19, PAX8, HSF4, and ZKSCAN3, which may contribute to the PPARγ-independent functions of ZFP407 in adipocytes and other cell types.

Early life stage exposure to fenbuconazole causes multigenerational cardiac developmental defects in zebrafish and potential reasons

With the increasing use of triazole fungicides in agriculture, triazole pesticides have aroused great concern about their toxicity and ecological risk. The current study investigated the impairments of embryonic exposure to fenbuconazole (FBZ) on cardiac transgenerational toxicity and related mechanisms. The fertilized eggs were exposed to 5, 50 and 500 ng/L FBZ for 72 h, and the larvae were then raised to adulthood in clean water. The adult fish were mated with unexposed fish to produce maternal and paternal F1 and F2 embryos, respectively. The results showed that increased arrhythmia were observed in F0, F1 and F2 larvae. Transcriptome sequencing indicated that the pathway of adrenergic signaling in cardiomyocytes was enriched in F0 and F2 larvae. In both F0 and F1 adult zebrafish hearts, ADRB2 protein expression decreased, and transcription of genes related to cardiac development and Ca2+ homeostasis was downregulated. These alterations might cause cardiac developmental defects. Significantly decreased protein levels of H3K9Ac and H3K14Ac might be linked with the downregulation in transcription of cardiac development genes. Protein‒protein interaction analysis exhibited that the pathway affecting the heart was well inherited in the paternal line. These results provide new ideas for the analysis and prevention of congenital heart disease.

Single-nucleotide variant calling in single-cell sequencing data with Monopogen

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Generation of human excitatory forebrain neurons by cooperative binding of proneural NGN2 and homeobox factor EMX1

Generation of defined neuronal subtypes from human pluripotent stem cells remains a challenge. The proneural factor NGN2 has been shown to overcome experimental variability observed by morphogen-guided differentiation and directly converts pluripotent stem cells into neurons, but their cellular heterogeneity has not been investigated yet. Here, we found that NGN2 reproducibly produces three different kinds of excitatory neurons characterized by partial coactivation of other neurotransmitter programs. We explored two principle approaches to achieve more precise specification: prepatterning the chromatin landscape that NGN2 is exposed to and combining NGN2 with region-specific transcription factors. Unexpectedly, the chromatin context of regionalized neural progenitors only mildly altered genomic NGN2 binding and its transcriptional response and did not affect neurotransmitter specification. In contrast, coexpression of region-specific homeobox factors such as EMX1 resulted in drastic redistribution of NGN2 including recruitment to homeobox targets and resulted in glutamatergic neurons with silenced nonglutamatergic programs. These results provide the molecular basis for a blueprint for improved strategies for generating a plethora of defined neuronal subpopulations from pluripotent stem cells for therapeutic or disease-modeling purposes.

TEA domain transcription factor 1(TEAD1) induces cardiac fibroblasts cells remodeling through BRD4/Wnt4 pathway

Cardiac fibroblasts (CFs) are the primary cells tasked with depositing and remodeling collagen and significantly associated with heart failure (HF). TEAD1 has been shown to be essential for heart development and homeostasis. However, fibroblast endogenous TEAD1 in cardiac remodeling remains incompletely understood. Transcriptomic analyses revealed consistently upregulated cardiac TEAD1 expression in mice 4 weeks after transverse aortic constriction (TAC) and Ang-II infusion. Further investigation revealed that CFs were the primary cell type expressing elevated TEAD1 levels in response to pressure overload. Conditional TEAD1 knockout was achieved by crossing TEAD1-floxed mice with CFs- and myofibroblasts-specific Cre mice. Echocardiographic and histological analyses demonstrated that CFs- and myofibroblasts-specific TEAD1 deficiency and treatment with TEAD1 inhibitor, VT103, ameliorated TAC-induced cardiac remodeling. Mechanistically, RNA-seq and ChIP-seq analysis identified Wnt4 as a novel TEAD1 target. TEAD1 has been shown to promote the fibroblast-to-myofibroblast transition through the Wnt signalling pathway, and genetic Wnt4 knockdown inhibited the pro-transformation phenotype in CFs with TEAD1 overexpression. Furthermore, co-immunoprecipitation combined with mass spectrometry, chromatin immunoprecipitation, and luciferase assays demonstrated interaction between TEAD1 and BET protein BRD4, leading to the binding and activation of the Wnt4 promoter. In conclusion, TEAD1 is an essential regulator of the pro-fibrotic CFs phenotype associated with pathological cardiac remodeling via the BRD4/Wnt4 signalling pathway.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

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.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

Antigen perception in T cells by long-term Erk and NFAT signaling dynamics

Immune system threat detection hinges on T cells' ability to perceive varying peptide-major histocompatibility complex (pMHC) antigens. As the Erk and NFAT pathways link T cell receptor engagement to gene regulation, their signaling dynamics may convey information about pMHC inputs. To test this idea, we developed a dual reporter mouse strain and a quantitative imaging assay that, together, enable simultaneous monitoring of Erk and NFAT dynamics in live T cells over day-long timescales as they respond to varying pMHC inputs. Both pathways initially activate uniformly across various pMHC inputs but diverge only over longer (9+ h) timescales, enabling independent encoding of pMHC affinity and dose. These late signaling dynamics are decoded via multiple temporal and combinatorial mechanisms to generate pMHC-specific transcriptional responses. Our findings underscore the importance of long timescale signaling dynamics in antigen perception and establish a framework for understanding T cell responses under diverse contexts.

Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies have mapped >90% of CVD-associated variants within the noncoding genome, which can alter the function of regulatory proteins, such as transcription factors (TFs). However, due to the overwhelming number of single-nucleotide polymorphisms (SNPs) (>500,000) in genome-wide association studies, prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1535 CVD-associated SNPs within TF footprints and putative cardiac enhancers plus 14,218 variants in linkage disequilibrium with genotype-dependent gene expression in cardiac tissues. Using ChIP-seq data from two cardiac TFs (NKX2-5 and TBX5) in human-induced pluripotent stem cell-derived cardiomyocytes, we trained a large-scale gapped k-mer SVM model to identify CVD-associated SNPs that altered NKX2-5 and TBX5 binding. The model was tested by scoring human heart TF genomic footprints within putative enhancers and measuring in vitro binding through electrophoretic mobility shift assay. Five variants predicted to alter NKX2-5 (rs59310144, rs6715570, and rs61872084) and TBX5 (rs7612445 and rs7790964) binding were prioritized for in vitro validation based on the magnitude of the predicted change in binding and are in cardiac tissue eQTLs. All five variants altered NKX2-5 and TBX5 DNA binding. We present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro analysis.

Stomatal cell fate commitment via transcriptional and epigenetic control: Timing is crucial

The formation of stomata presents a compelling model system for comprehending the initiation, proliferation, commitment and differentiation of de novo lineage-specific stem cells. Precise, timely and robust cell fate and identity decisions are crucial for the proper progression and differentiation of functional stomata. Deviations from this precise specification result in developmental abnormalities and nonfunctional stomata. However, the molecular underpinnings of timely cell fate commitment have just begun to be unravelled. In this review, we explore the key regulatory strategies governing cell fate commitment, emphasizing the distinctions between embryonic and postembryonic stomatal development. Furthermore, the interplay of transcription factors and cell cycle machineries is pivotal in specifying the transition into differentiation. We aim to synthesize recent studies utilizing single-cell as well as cell-type-specific transcriptomics, epigenomics and chromatin accessibility profiling to shed light on how master-regulatory transcription factors and epigenetic machineries mutually influence each other to drive fate commitment and maintenance.

c-JUN is a barrier in hESC to cardiomyocyte transition

Loss of c-JUN leads to early mouse embryonic death, possibly because of a failure to develop a normal cardiac system. How c-JUN regulates human cardiomyocyte cell fate remains unknown. Here, we used the in vitro differentiation of human pluripotent stem cells into cardiomyocytes to study the role of c-JUN. Surprisingly, the knockout of c-JUN improved cardiomyocyte generation, as determined by the number of TNNT2+ cells. ATAC-seq data showed that the c-JUN defect led to increased chromatin accessibility on critical regulatory elements related to cardiomyocyte development. ChIP-seq data showed that the knockout c-JUN increased RBBP5 and SETD1B expression, leading to improved H3K4me3 deposition on key genes that regulate cardiogenesis. The c-JUN KO phenotype could be copied using the histone demethylase inhibitor CPI-455, which also up-regulated H3K4me3 levels and increased cardiomyocyte generation. Single-cell RNA-seq data defined three cell branches, and knockout c-JUN activated more regulons that are related to cardiogenesis. In summary, our data demonstrated that c-JUN could regulate cardiomyocyte cell fate by modulating H3K4me3 modification and chromatin accessibility and shed light on how c-JUN regulates heart development in humans.

Prioritizing Cardiovascular Disease-Associated Variants Altering NKX2-5 Binding through an Integrative Computational Approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies (GWAS) have mapped > 90% of CVD-associated variants within the non-coding genome, which can alter the function of regulatory proteins, like transcription factors (TFs). However, due to the overwhelming number of GWAS single nucleotide polymorphisms (SNPs) (>500,000), prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1,535 CVD-associated SNPs that occur within human heart footprints/enhancers and 9,309 variants in linkage disequilibrium (LD) with differential gene expression profiles in cardiac tissue. Using hiPSC-CM ChIP-seq data from NKX2-5 and TBX5, two cardiac TFs essential for proper heart development, we trained a large-scale gapped k-mer SVM (LS-GKM-SVM) predictive model that can identify binding sites altered by CVD-associated SNPs. The computational predictive model was tested by scoring human heart footprints and enhancers in vitro through electrophoretic mobility shift assay (EMSA). Three variants (rs59310144, rs6715570, and rs61872084) were prioritized for in vitro validation based on their eQTL in cardiac tissue and LS-GKM-SVM prediction to alter NKX2-5 DNA binding. All three variants altered NKX2-5 DNA binding. In summary, we present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro experimental analysis.

Modeling a variant of unknown significance in the Drosophila ortholog of the human cardiogenic gene NKX2.5

Sequencing of human genome samples has unearthed genetic variants for which functional testing is necessary to validate their clinical significance. We used the Drosophila system to analyze a variant of unknown significance in the human congenital heart disease gene NKX2.5 (also known as NKX2-5). We generated an R321N allele of the NKX2.5 ortholog tinman (tin) to model a human K158N variant and tested its function in vitro and in vivo. The R321N Tin isoform bound poorly to DNA in vitro and was deficient in activating a Tin-dependent enhancer in tissue culture. Mutant Tin also showed a significantly reduced interaction with a Drosophila T-box cardiac factor named Dorsocross1. We generated a tinR321N allele using CRISPR/Cas9, for which homozygotes were viable and had normal heart specification, but showed defects in the differentiation of the adult heart that were exacerbated by further loss of tin function. We propose that the human K158N variant is pathogenic through causing a deficiency in DNA binding and a reduced ability to interact with a cardiac co-factor, and that cardiac defects might arise later in development or adult life.

Cardiac competence of the paraxial head mesoderm fades concomitant with a shift towards the head skeletal muscle programme

The vertebrate head mesoderm provides the heart, the great vessels, some smooth and most head skeletal muscle, in addition to parts of the skull. It has been speculated that the ability to generate cardiac and smooth muscle is the evolutionary ground-state of the tissue. However, whether indeed the entire head mesoderm has generic cardiac competence, how long this may last, and what happens as cardiac competence fades, is not clear. Bone morphogenetic proteins (Bmps) are known to promote cardiogenesis. Using 41 different marker genes in the chicken embryo, we show that the paraxial head mesoderm that normally does not engage in cardiogenesis has the ability to respond to Bmp for a long time. However, Bmp signals are interpreted differently at different time points. Up to early head fold stages, the paraxial head mesoderm is able to read Bmps as signal to engage in the cardiac programme; the ability to upregulate smooth muscle markers is retained slightly longer. Notably, as cardiac competence fades, Bmp promotes the head skeletal muscle programme instead. The switch from cardiac to skeletal muscle competence is Wnt-independent as Wnt caudalises the head mesoderm and also suppresses Msc-inducing Bmp provided by the prechordal plate, thus suppressing both the cardiac and the head skeletal muscle programmes. Our study for the first time suggests a specific transition state in the embryo when cardiac competence is replaced by skeletal muscle competence. It sets the stage to unravel the cardiac-skeletal muscle antagonism that is known to partially collapse in heart failure.

Antigen perception in T cells by long-term Erk and NFAT signaling dynamics

Immune system threat detection hinges on T cells' ability to perceive varying peptide major-histocompatibility complex (pMHC) antigens. As the Erk and NFAT pathways link T cell receptor engagement to gene regulation, their signaling dynamics may convey information about pMHC inputs. To test this idea, we developed a dual reporter mouse strain and a quantitative imaging assay that, together, enable simultaneous monitoring of Erk and NFAT dynamics in live T cells over day-long timescales as they respond to varying pMHC inputs. Both pathways initially activate uniformly across various pMHC inputs, but diverge only over longer (9+ hrs) timescales, enabling independent encoding of pMHC affinity and dose. These late signaling dynamics are decoded via multiple temporal and combinatorial mechanisms to generate pMHC-specific transcriptional responses. Our findings underscore the importance of long timescale signaling dynamics in antigen perception, and establish a framework for understanding T cell responses under diverse contexts.

SIGNIFICANCE STATEMENT

To counter diverse pathogens, T cells mount distinct responses to varying peptide-major histocompatibility complex ligands (pMHCs). They perceive the affinity of pMHCs for the T cell receptor (TCR), which reflects its foreignness, as well as pMHC abundance. By tracking signaling responses in single living cells to different pMHCs, we find that T cells can independently perceive pMHC affinity vs dose, and encode this information through the dynamics of Erk and NFAT signaling pathways downstream of the TCR. These dynamics are jointly decoded by gene regulatory mechanisms to produce pMHC-specific activation responses. Our work reveals how T cells can elicit tailored functional responses to diverse threats and how dysregulation of these responses may lead to immune pathologies.

CHD-associated enhancers shape human cardiomyocyte lineage commitment

Enhancers orchestrate gene expression programs that drive multicellular development and lineage commitment. Thus, genetic variants at enhancers are thought to contribute to developmental diseases by altering cell fate commitment. However, while many variant-containing enhancers have been identified, studies to endogenously test the impact of these enhancers on lineage commitment have been lacking. We perform a single-cell CRISPRi screen to assess the endogenous roles of 25 enhancers and putative cardiac target genes implicated in genetic studies of congenital heart defects (CHDs). We identify 16 enhancers whose repression leads to deficient differentiation of human cardiomyocytes (CMs). A focused CRISPRi validation screen shows that repression of TBX5 enhancers delays the transcriptional switch from mid- to late-stage CM states. Endogenous genetic deletions of two TBX5 enhancers phenocopy epigenetic perturbations. Together, these results identify critical enhancers of cardiac development and suggest that misregulation of these enhancers could contribute to cardiac defects in human patients.

Tbx5 maintains atrial identity by regulating an atrial enhancer network

Understanding how the atrial and ventricular chambers of the heart maintain their distinct identity is a prerequisite for treating chamber-specific diseases. Here, we selectively inactivated the transcription factor Tbx5 in the atrial working myocardium of the neonatal mouse heart to show that it is required to maintain atrial identity. Atrial Tbx5 inactivation downregulated highly chamber specific genes such as Myl7 and Nppa , and conversely, increased the expression of ventricular identity genes including Myl2 . Using combined single nucleus transcriptome and open chromatin profiling, we assessed genomic accessibility changes underlying the altered atrial identity expression program, identifying 1846 genomic loci with greater accessibility in control atrial cardiomyocytes compared to KO aCMs. 69% of the control-enriched ATAC regions were bound by TBX5, demonstrating a role for TBX5 in maintaining atrial genomic accessibility. These regions were associated with genes that had higher expression in control aCMs compared to KO aCMs, suggesting they act as TBX5-dependent enhancers. We tested this hypothesis by analyzing enhancer chromatin looping using HiChIP and found 510 chromatin loops that were sensitive to TBX5 dosage. Of the loops enriched in control aCMs, 73.7% contained anchors in control-enriched ATAC regions. Together, these data demonstrate a genomic role for TBX5 in maintaining the atrial gene expression program by binding to atrial enhancers and preserving tissue-specific chromatin architecture of atrial enhancers.

The long noncoding RNA Charme supervises cardiomyocyte maturation by controlling cell differentiation programs in the developing heart

Long noncoding RNAs (lncRNAs) are emerging as critical regulators of heart physiology and disease, although the studies unveiling their modes of action are still limited to few examples. We recently identified pCharme, a chromatin-associated lncRNA whose functional knockout in mice results in defective myogenesis and morphological remodeling of the cardiac muscle. Here, we combined Cap-Analysis of Gene Expression (CAGE), single-cell (sc)RNA sequencing, and whole-mount in situ hybridization analyses to study pCharme cardiac expression. Since the early steps of cardiomyogenesis, we found the lncRNA being specifically restricted to cardiomyocytes, where it assists the formation of specific nuclear condensates containing MATR3, as well as important RNAs for cardiac development. In line with the functional significance of these activities, pCharme ablation in mice results in a delayed maturation of cardiomyocytes, which ultimately leads to morphological alterations of the ventricular myocardium. Since congenital anomalies in myocardium are clinically relevant in humans and predispose patients to major complications, the identification of novel genes controlling cardiac morphology becomes crucial. Our study offers unique insights into a novel lncRNA-mediated regulatory mechanism promoting cardiomyocyte maturation and bears relevance to Charme locus for future theranostic applications.

An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5^RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5^RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

Deciphering Transcriptional Networks during Human Cardiac Development

Human heart development is governed by transcription factor (TF) networks controlling dynamic and temporal gene expression alterations. Therefore, to comprehensively characterize these transcriptional regulations, day-to-day transcriptomic profiles were generated throughout the directed cardiac differentiation, starting from three distinct human- induced pluripotent stem cell lines from healthy donors (32 days). We applied an expression-based correlation score to the chronological expression profiles of the TF genes, and clustered them into 12 sequential gene expression waves. We then identified a regulatory network of more than 23,000 activation and inhibition links between 216 TFs. Within this network, we observed previously unknown inferred transcriptional activations linking IRX3 and IRX5 TFs to three master cardiac TFs: GATA4, NKX2-5 and TBX5. Luciferase and co-immunoprecipitation assays demonstrated that these five TFs could (1) activate each other's expression; (2) interact physically as multiprotein complexes; and (3) together, finely regulate the expression of SCN5A, encoding the major cardiac sodium channel. Altogether, these results unveiled thousands of interactions between TFs, generating multiple robust hypotheses governing human cardiac development.

Toward improved understanding of cardiac development and congenital heart disease: The advent of cardiac organoids

Human cardiac organoid systems hold significant promise for mechanistic studies of early heart morphogenesis and an improved understanding of congenital cardiac disease.

During the past decade, we have witnessed remarkable progress in genome editing technology, stem cell research, and bioengineering. The fundamental basic research discoveries accelerate rapidly into clinical translation, paving the way for myocardial regeneration, better understanding of the structural heart disease, and bioengineering of heart structures and even entire hearts. The new horizon is vast and diverse, ranging from creating universal stem cell biobanking to genome edited heart xenotransplantation. Herein, a group of experts from Duke University discuss the state of the art and the possible influence of cardiac organoids on our understanding of structural heart disease. It may not be immediately clear now in what practical ways this technology will be translated into our daily work, yet the current progress in bioengineering will likely have a very significant influence on our surgical practice.

Igor E. Konstantinov, MD, PhD, FRACS

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells–that also express NANOG–to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

The authors show that the transcription factors NANOG and GATA6 co-bind the same enhancers in common progenitors before divergent epiblast and primitive endoderm lineages emerge. This may help maintain plasticity at early stages and facilitate bifurcation into distinct lineages

Identification and functional analysis of genetic variants of ISL1 gene promoter in human atrial septal defects

BACKGROUND

Atrial septal defect (ASD) is a common type of congenital heart disease. A gene promoter plays pivotal role in the disease development. This study was designed to investigate the pathological role of variants of the ISL1 gene promoter region in ASD patients.

METHODS

Total DNA extracted from 625 subjects, including 332 ASD patients and 293 healthy controls, was sequenced to identify variants in the promoter region of ISL1 gene. Further functional analyses of the variants were performed with dual luciferase reporter assay and electrophoretic mobility shift assay (EMSA). All possible binding sites of transcription factor affected by the identified variants were predicted using the JASPAR database.

RESULTS

Four variants in the ISL1 gene promoter were found only in patients with ASD by sequencing. Three of the four variants [g.4923 G > C (rs541081886), g.5079 A > G (rs1371835943) and g.5309 G > A (rs116222082)] significantly decreased the transcriptional activities compared with the wild-type ISL1 gene promoter (p < 0.05). The EMSA revealed that these variants [g.4923 G > C (rs541081886), g.5079 A > G (rs1371835943) and g.5309 G > A (rs116222082)] in the ISL1 gene promoter affected the number and affinity of binding sites of transcription factors. Further analysis with the online JASPAR database demonstrated that a cluster of putative binding sites for transcription factors may be altered by these variants.

CONCLUSIONS

These sequence variants identified from the promoter region of ISL1 gene in ASD patients are probably involved in the development of ASD by affecting the transcriptional activity and altering ISL1 levels. Therefore, these findings may provide new insights into the molecular etiology and potential therapeutic strategy of ASD.

Dynamic chromatin accessibility deploys heterotypic cis/trans-acting factors driving stomatal cell-fate commitment

Chromatin architecture and transcription factor (TF) binding underpin cell-fate specification during development, but their mutual regulatory relationships remain unclear. Here we report an atlas of dynamic chromatin landscapes during stomatal cell-lineage progression, in which sequential cell-state transitions are governed by lineage-specific bHLH TFs. Major reprogramming of chromatin accessibility occurs at the proliferation-to-differentiation transition. We discover novel co-cis regulatory elements (CREs) signifying the early precursor stage, BBR/BPC (GAGA) and bHLH (E-box) motifs, where master-regulatory bHLH TFs, SPEECHLESS and MUTE, consecutively bind to initiate and terminate the proliferative state, respectively. BPC TFs complex with MUTE to repress SPEECHLESS expression through a local deposition of repressive histone marks. We elucidate the mechanism by which cell-state-specific heterotypic TF complexes facilitate cell-fate commitment by recruiting chromatin modifiers via key co-CREs.

Cardiac Development in the Presence of Cadmium: An in Vitro Study Using Human Embryonic Stem Cells and Cardiac Organoids

BACKGROUND

Exposure to cadmium (Cd) is associated with cardiovascular diseases. Maternal Cd exposure is a significant risk factor for congenital heart disease. However, mechanisms of Cd on developmental cardiotoxicity are not well defined.

OBJECTIVES

We evaluated the effects of Cd on the different stages (mesoderm, cardiac induction, cardiac function) of cardiac development using an early embryo development in vitro model and two- or three-dimensional (2- or 3D) cardiomyocyte and cardiac organoid formation models mimicking early cardiac development.

METHODS

Embryonic stem cells (ESCs) form 3D aggregates, called embryoid bodies, that recapitulate events involved with early embryogenesis (e.g., germ layer formation). This model was used for early germ layer formation and signaling pathway identification. The 2D cardiomyocyte differentiation from the NKX2-5eGFP/w human ESCs model was used to explore the effects of Cd exposure on cardiomyocyte formation and to model mesoderm differentiation and cardiac induction, allowing us to explore different developmental windows of Cd toxicity. The 3D cardiac organoid model was used in evaluating the effects of Cd exposure on contractility and cardiac development.

RESULTS

Cd (0.6μM; 110 ppb) lowered the differentiation of embryoid bodies to mesoderm via suppression of Wnt/β-catenin-signaling pathways. During early mesoderm induction, the mesoderm-associated transcription factors MESP1 and EOMES showed a transient up-regulation, which decreased later in the cardiac induction stage. Cd (0.15μM) lowered mesoderm formation and cardiac induction through suppression of the transcription factors and mesoderm marker genes HAND1, SNAI2, HOPX, and the cardiac-specific genes NKX2-5, GATA4, troponin T, and alpha-actinin. In addition, Cd-induced histone modifications for both gene activation (H3K4me3) and repression (H3K27me3), which play vital roles in regulating mesoderm commitment markers. The effects of Cd inhibition on cardiomyocyte differentiation were confirmed in 3D cardiac organoids.

DISCUSSION

In conclusion, using a human ESC-derived 2D/3D in vitro differentiation model system and cardiac organoids, we demonstrated that low-dose Cd suppressed mesoderm formation through mesoderm gene histone modification, thus inhibiting cardiomyocyte differentiation and cardiac induction. The studies provide valuable insights into cellular events and molecular mechanisms associated with Cd-induced congenital heart disease. https://doi.org/10.1289/EHP11208.

Regulation of cardiac ion channels by transcription factors: Looking for new opportunities of druggable targets for the treatment of arrhythmias

Cardiac electrical activity is governed by different ion channels that generate action potentials. Acquired or inherited abnormalities in the expression and/or function of ion channels usually result in electrophysiological changes that can cause cardiac arrhythmias. Transcription factors (TFs) control gene transcription by binding to specific DNA sequences adjacent to target genes. Linkage analysis, candidate-gene screening within families, and genome-wide association studies have linked rare and common genetic variants in the genes encoding TFs with genetically-determined cardiac arrhythmias. Besides its critical role in cardiac development, recent data demonstrated that they control cardiac electrical activity through the direct regulation of the expression and function of cardiac ion channels in adult hearts. This narrative review summarizes some studies showing functional data on regulation of the main human atrial and ventricular Na⁺, Ca²⁺, and K⁺ channels by cardiac TFs such as Pitx2c, Tbx20, Tbx5, Zfhx3, among others. The results have improved our understanding of the mechanisms regulating cardiac electrical activity and may open new avenues for therapeutic interventions in cardiac acquired or inherited arrhythmias through the identification of TFs as potential drug targets. Even though TFs have for a long time been considered as 'undruggable' targets, advances in structural biology have led to the identification of unique pockets in TFs amenable to be targeted with small-molecule drugs or peptides that are emerging as novel therapeutic drugs.

The TH1 cell lineage-determining transcription factor T-bet suppresses TH2 gene expression by redistributing GATA3 away from TH2 genes

Lineage-determining transcription factors (LD-TFs) drive the differentiation of progenitor cells into a specific lineage. In CD4+ T cells, T-bet dictates differentiation of the TH1 lineage, whereas GATA3 drives differentiation of the alternative TH2 lineage. However, LD-TFs, including T-bet and GATA3, are frequently co-expressed but how this affects LD-TF function is not known. By expressing T-bet and GATA3 separately or together in mouse T cells, we show that T-bet sequesters GATA3 at its target sites, thereby removing GATA3 from TH2 genes. This redistribution of GATA3 is independent of GATA3 DNA binding activity and is instead mediated by the T-bet DNA binding domain, which interacts with the GATA3 DNA binding domain and changes GATA3's sequence binding preference. This mechanism allows T-bet to drive the TH1 gene expression program in the presence of GATA3. We propose that redistribution of one LD-TF by another may be a common mechanism that could explain how specific cell fate choices can be made even in the presence of other transcription factors driving alternative differentiation pathways.

CHD4 is recruited by GATA4 and NKX2-5 to repress noncardiac gene programs in the developing heart

The nucleosome remodeling and deacetylase (NuRD) complex is one of the central chromatin remodeling complexes that mediates gene repression. NuRD is essential for numerous developmental events, including heart development. Clinical and genetic studies have provided direct evidence for the role of chromodomain helicase DNA-binding protein 4 (CHD4), the catalytic component of NuRD, in congenital heart disease (CHD), including atrial and ventricular septal defects. Furthermore, it has been demonstrated that CHD4 is essential for mammalian cardiomyocyte formation and function. A key unresolved question is how CHD4/NuRD is localized to specific cardiac target genes, as neither CHD4 nor NuRD can directly bind DNA. Here, we coupled a bioinformatics-based approach with mass spectrometry analyses to demonstrate that CHD4 interacts with the core cardiac transcription factors GATA4, NKX2-5, and TBX5 during embryonic heart development. Using transcriptomics and genome-wide occupancy data, we characterized the genomic landscape of GATA4, NKX2-5, and TBX5 repression and defined the direct cardiac gene targets of the GATA4-CHD4, NKX2-5-CHD4, and TBX5-CHD4 complexes. These data were used to identify putative cis-regulatory elements controlled by these complexes. We genetically interrogated two of these silencers in vivo: Acta1 and Myh11 We show that deletion of these silencers leads to inappropriate skeletal and smooth muscle gene misexpression, respectively, in the embryonic heart. These results delineate how CHD4/NuRD is localized to specific cardiac loci and explicates how mutations in the broadly expressed CHD4 protein lead to cardiac-specific disease states.

Pioneer factors as master regulators of the epigenome and cell fate

Pioneer factors are transcription factors with the unique ability to initiate opening of closed chromatin. The stability of cell identity relies on robust mechanisms that maintain the epigenome and chromatin accessibility to transcription factors. Pioneer factors counter these mechanisms to implement new cell fates through binding of DNA target sites in closed chromatin and introduction of active-chromatin histone modifications, primarily at enhancers. As master regulators of enhancer activation, pioneers are thus crucial for the implementation of correct cell fate decisions in development, and as such, they hold tremendous potential for therapy through cellular reprogramming. The power of pioneer factors to reshape the epigenome also presents an Achilles heel, as their misexpression has major pathological consequences, such as in cancer. In this Review, we discuss the emerging mechanisms of pioneer factor functions and their roles in cell fate specification, cellular reprogramming and cancer.

Transcription factor protein interactomes reveal genetic determinants in heart disease

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

PRDM16 Is a Compact Myocardium-Enriched Transcription Factor Required to Maintain Compact Myocardial Cardiomyocyte Identity in Left Ventricle

BACKGROUND

Left ventricular noncompaction cardiomyopathy (LVNC) was discovered half a century ago as a cardiomyopathy with excessive trabeculation and a thin ventricular wall. In the decades since, numerous studies have demonstrated that LVNC primarily has an effect on left ventricles (LVs) and is often associated with LV dilation and dysfunction. However, in part because of the lack of suitable mouse models that faithfully mirror the selective LV vulnerability in patients, mechanisms underlying the susceptibility of LVs to dilation and dysfunction in LVNC remain unknown. Genetic studies have revealed that deletions and mutations in PRDM16 (PR domain-containing 16) cause LVNC, but previous conditional Prdm16 knockout mouse models do not mirror the LVNC phenotype in patients, and the underlying molecular mechanisms by which PRDM16 deficiency causes LVNC are still unclear.

METHODS

Prdm16 cardiomyocyte-specific knockout (Prdm16cKO) mice were generated and analyzed for cardiac phenotypes. RNA sequencing and chromatin immunoprecipitation deep sequencing were performed to identify direct transcriptional targets of PRDM16 in cardiomyocytes. Single-cell RNA sequencing in combination with spatial transcriptomics was used to determine cardiomyocyte identity at the single-cell level.

RESULTS

Cardiomyocyte-specific ablation of Prdm16 in mice caused LV-specific dilation and dysfunction, as well as biventricular noncompaction, which fully recapitulated LVNC in patients. PRDM16 functioned mechanistically as a compact myocardium-enriched transcription factor that activated compact myocardial genes while repressing trabecular myocardial genes in LV compact myocardium. Consequently, Prdm16cKO LV compact myocardial cardiomyocytes shifted from their normal transcriptomic identity to a transcriptional signature resembling trabecular myocardial cardiomyocytes or neurons. Chamber-specific transcriptional regulation by PRDM16 was attributable in part to its cooperation with LV-enriched transcription factors Tbx5 and Hand1.

CONCLUSIONS

These results demonstrate that disruption of proper specification of compact cardiomyocytes may play a key role in the pathogenesis of LVNC. They also shed light on underlying mechanisms of the LV-restricted transcriptional program governing LV chamber growth and maturation, providing a tangible explanation for the susceptibility of LV in a subset of LVNC cardiomyopathies.

Genomic enhancers in cardiac development and disease

The Human Genome Project marked a major milestone in the scientific community as it unravelled the ~3 billion bases that are central to crucial aspects of human life. Despite this achievement, it only scratched the surface of understanding how each nucleotide matters, both individually and as part of a larger unit. Beyond the coding genome, which comprises only ~2% of the whole genome, scientists have realized that large portions of the genome, not known to code for any protein, were crucial for regulating the coding genes. These large portions of the genome comprise the 'non-coding genome'. The history of gene regulation mediated by proteins that bind to the regulatory non-coding genome dates back many decades to the 1960s. However, the original definition of 'enhancers' was first used in the early 1980s. In this Review, we summarize benchmark studies that have mapped the role of cardiac enhancers in disease and development. We highlight instances in which enhancer-localized genetic variants explain the missing link to cardiac pathogenesis. Finally, we inspire readers to consider the next phase of exploring enhancer-based gene therapy for cardiovascular disease.

Nkx2-3 induces autophagy inhibiting proliferation and migration of vascular smooth muscle cells via AMPK/mTOR signaling pathway

Vascular remodeling and restenosis are common complications after percutaneous coronary intervention. Excessive proliferation and migration of vascular smooth muscle cells (VSMCs) play important roles in intimal hyperplasia-induced vascular restenosis. NK2 Homeobox 3 (Nkx2-3), a critical member of Nkx family, is involved in tissue differentiation and organ development. However, the role of Nkx2-3 in VSMCs proliferation and migration remains unknown. In this study, we used carotid balloon injury model and platelet-derived growth factor-BB (PDGF)-treated VSMCs as in vivo and in vitro experimental models. EdU assay and CCK-8 assay were used to detect cell proliferation. Migration was measured by scratch test. Hematoxylin and eosin staining and immunohistochemistry staining were used to evaluate the intimal hyperplasia. The autophagy level was detected by fluorescent mRFP-GFP-LC3 in vitro and by transmission electron microscopy in vivo. It was shown that Nkx2-3 was upregulated both in balloon injured carotid arteries and PDGF-stimulated VSMCs. Adenovirus-mediated Nkx2-3 overexpression inhibited intimal hyperplasia after balloon injury, and suppressed VSMCs proliferation and migration induced by PDGF. Conversely, silencing of Nkx2-3 by small interfering RNA exaggerated proliferation and migration of VSMCs. Furthermore, we found that Nkx2-3 enhanced autophagy level, while the autophagy inhibitor 3-MA eliminated the inhibitory effect of Nkx2-3 on VSMCs proliferation and migration both in vivo and in vitro. Moreover, Nkx2-3 promoted autophagy in VSMCs by activating the AMP-activated protein kinase/mammalian target of rapamycin (AMPK/mTOR) signaling pathway. These results demonstrated for the first time that Nkx2-3 inhibited VSMCs proliferation and migration through AMPK/mTOR-mediated autophagy.

Control of cardiomyocyte differentiation timing by intercellular signaling pathways

Congenital Heart Disease (CHD), malformations of the heart present at birth, is the most common class of life-threatening birth defect (Hoffman (1995) [1], Gelb (2004) [2], Gelb (2014) [3]). A major research challenge is to elucidate the genetic determinants of CHD and mechanistically link CHD ontogeny to a molecular understanding of heart development. Although the embryonic origins of CHD are unclear in most cases, dysregulation of cardiovascular lineage specification, patterning, proliferation, migration or differentiation have been described (Olson (2004) [4], Olson (2006) [5], Srivastava (2006) [6], Dunwoodie (2007) [7], Bruneau (2008) [8]). Cardiac differentiation is the process whereby cells become progressively more dedicated in a trajectory through the cardiac lineage towards mature cardiomyocytes. Defects in cardiac differentiation have been linked to CHD, although how the complex control of cardiac differentiation prevents CHD is just beginning to be understood. The stages of cardiac differentiation are highly stereotyped and have been well-characterized (Kattman et al. (2011) [9], Wamstad et al. (2012) [10], Luna-Zurita et al. (2016) [11], Loh et al. (2016) [12], DeLaughter et al. (2016) [13]); however, the developmental and molecular mechanisms that promote or delay the transition of a cell through these stages have not been as deeply investigated. Tight temporal control of progenitor differentiation is critically important for normal organ size, spatial organization, and cellular physiology and homeostasis of all organ systems (Raff et al. (1985) [14], Amthor et al. (1998) [15], Kopan et al. (2014) [16]). This review will focus on the action of signaling pathways in the control of cardiomyocyte differentiation timing. Numerous signaling pathways, including the Wnt, Fibroblast Growth Factor, Hedgehog, Bone Morphogenetic Protein, Insulin-like Growth Factor, Thyroid Hormone and Hippo pathways, have all been implicated in promoting or inhibiting transitions along the cardiac differentiation trajectory. Gaining a deeper understanding of the mechanisms controlling cardiac differentiation timing promises to yield insights into the etiology of CHD and to inform approaches to restore function to damaged hearts.

Functional analysis of novel genetic variants of NKX2-5 associated with nonsyndromic congenital heart disease

NKX2-5, a master cardiac regulatory transcription factor was the first known genetic cause of congenital heart diseases (CHDs). To further investigate its role in CHD pathogenesis, we performed mutational screening of 285 CHD probands and 200 healthy controls. Five coding sequence variants were identified in six CHD cases (2.1%), including three in the N-terminal region (p.A61G, p.R95L, and p.E131K) and one each in homeodomain (HD) (p.A148E) and tyrosine-rich domain (p.P247A). Variant-p.A148E showed tertiary structure changes and differential DNA binding affinity of mutant compared to wild type. Two N-terminal variants-p.A61G and p.E131K along with HD variant p.A148E demonstrated significantly reduced transcriptional activity of Nppa and Actc1 promoters in dual luciferase promoter assay supported by their reduced expression in qRT-PCR. Nonetheless, variant p.R95L affected the synergy of NKX2-5 with serum response factor and TBX5 leading to significantly decreased Actc1 promoter activity depicting a distinctive role of this region. The aberrant expression of other target genes-Irx4, Mef2c, Bmp10, Myh6, Myh7, and Myocd is also observed in response to NKX2-5 variants, possibly due to the defective gene regulatory network. Severely impaired downstream promoter activities and abnormal expression of target genes due to N-terminal variants supports the emerging role of this region during cardiac-developmental pathways.

Intra- and intercellular signaling pathways associated with drug-induced cardiac pathophysiology

Cardiac physiology and homeostasis are maintained by the interaction of multiple cell types, via both intra- and intercellular signaling pathways. Perturbations in these signaling pathways induced by oncology therapies can reduce cardiac function, ultimately leading to heart failure. As cancer survival increases, related cardiovascular complications are becoming increasingly prevalent, thus identifying the perturbations and cell signaling drivers of cardiotoxicity is increasingly important. Here, we discuss the homotypic and heterotypic cellular interactions that form the basis of intra- and intercellular cardiac signaling pathways, and how oncological agents disrupt these pathways, leading to heart failure. We also highlight the emerging systems biology techniques that can be applied, enabling a deeper understanding of the intra- and intercellular signaling pathways across multiple cell types associated with cardiovascular toxicity.

iPSC Modeling of RBM20-Deficient DCM Identifies Upregulation of RBM20 as a Therapeutic Strategy

Recent advances in induced pluripotent stem cell (iPSC) technology and directed differentiation of iPSCs into cardiomyocytes (iPSC-CMs) make it possible to model genetic heart disease in vitro. We apply CRISPR/Cas9 genome editing technology to introduce three RBM20 mutations in iPSCs and differentiate them into iPSC-CMs to establish an in vitro model of RBM20 mutant dilated cardiomyopathy (DCM). In iPSC-CMs harboring a known causal RBM20 variant, the splicing of RBM20 target genes, calcium handling, and contractility are impaired consistent with the disease manifestation in patients. A variant (Pro633Leu) identified by exome sequencing of patient genomes displays the same disease phenotypes, thus establishing this variant as disease causing. We find that all-trans retinoic acid upregulates RBM20 expression and reverts the splicing, calcium handling, and contractility defects in iPSC-CMs with different causal RBM20 mutations. These results suggest that pharmacological upregulation of RBM20 expression is a promising therapeutic strategy for DCM patients with a heterozygous mutation in RBM20.

Briganti et al. use iPSC and CRISPR/Cas9 to create a model of RBM20-deficient dilated cardiomyopathy (DCM) that recapitulates mRNA splicing and contractile defects of the disease. They evaluate pharmacological upregulation of RBM20 as a therapeutic strategy. All-trans retinoic acid upregulates RBM20 expression and ameliorates the in vitro hallmarks of disease.

Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors

Meis1 and Meis2 are homeodomain transcription factors that regulate organogenesis through cooperation with Hox proteins. Elimination of Meis genes after limb induction has shown their role in limb proximo-distal patterning; however, limb development in the complete absence of Meis function has not been studied. Here, we report that Meis1/2 inactivation in the lateral plate mesoderm of mouse embryos leads to limb agenesis. Meis and Tbx factors converge in this function, extensively co-binding with Tbx to genomic sites and co-regulating enhancers of Fgf10, a critical factor in limb initiation. Limbs with three deleted Meis alleles show proximal-specific skeletal hypoplasia and agenesis of posterior skeletal elements. This failure in posterior specification results from an early role of Meis factors in establishing the limb antero-posterior prepattern required for Shh activation. Our results demonstrate roles for Meis transcription factors in early limb development and identify their involvement in previously undescribed interaction networks that regulate organogenesis.

Variant Intronic Enhancer Controls SCN10A-short Expression and Heart Conduction

BACKGROUND

Genetic variants in SCN10A, encoding the neuronal voltage-gated sodium channel Na_V1.8, are strongly associated with atrial fibrillation, Brugada syndrome, cardiac conduction velocities, and heart rate. The cardiac function of SCN10A has not been resolved, however, and diverging mechanisms have been proposed. Here, we investigated the cardiac expression of SCN10A and the function of a variant-sensitive intronic enhancer previously linked to the regulation of SCN5A, encoding the major essential cardiac sodium channel Na_V1.5.

METHODS

The expression of SCN10A was investigated in mouse and human hearts. With the use of CRISPR/Cas9 genome editing, the mouse intronic enhancer was disrupted, and mutant mice were characterized by transcriptomic and electrophysiological analyses. The association of genetic variants at SCN5A-SCN10A enhancer regions and gene expression were evaluated by genome-wide association studies single-nucleotide polymorphism mapping and expression quantitative trait loci analysis.

RESULTS

We found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene (Scn10a-short). Transcription occurs from an intronic enhancer-promoter complex, whereas full-length Scn10a transcript was undetectable in the human and mouse heart. Expression quantitative trait loci analysis revealed that the genetic variants in linkage disequilibrium with genetic variant rs6801957 in the intronic enhancer associate with SCN10A transcript levels in the heart. Genetic modification of the enhancer in the mouse genome led to reduced cardiac Scn10a-short expression in atria and ventricles, reduced cardiac sodium current in atrial cardiomyocytes, atrial conduction slowing and arrhythmia, whereas the expression of Scn5a, the presumed enhancer target gene, remained unaffected. In patch-clamp transfection experiments, expression of Scn10a-short-encoded Na_V1.8-short increased Na_V1.5-mediated sodium current. We propose that noncoding genetic variation modulates transcriptional regulation of Scn10a-short in cardiomyocytes that impacts Na_V1.5-mediated sodium current and heart rhythm.

CONCLUSIONS

Genetic variants in and around SCN10A modulate enhancer function and expression of a cardiac-specific SCN10A-short transcript. We propose that noncoding genetic variation modulates transcriptional regulation of a functional C-terminal portion of Na_V1.8 in cardiomyocytes that impacts on Na_V1.5 function, cardiac conduction velocities, and arrhythmia susceptibility.

Logic and lineage impacts on functional transcription factor deployment for T-cell fate commitment

Transcription factors are the major agents that read the regulatory sequence information in the genome to initiate changes in expression of specific genes, both in development and in physiological activation responses. Their actions depend on site-specific DNA binding and are largely guided by their individual DNA target sequence specificities. However, their action is far more conditional in a real developmental context than would be expected for simple reading of local genomic DNA sequence, which is common to all cells in the organism. They are constrained by slow-changing chromatin states and by interactions with other transcription factors, which affect their occupancy patterns of potential sites across the genome. These mechanisms lead to emergent discontinuities in function even for transcription factors with minimally changing expression. This is well revealed by diverse lineages of blood cells developing throughout life from hematopoietic stem cells, which use overlapping combinations of transcription factors to drive strongly divergent gene regulation programs. Here, using development of T lymphocytes from hematopoietic multipotent progenitor cells as a focus, recent evidence is reviewed on how binding specificity and dynamics, transcription factor cooperativity, and chromatin state changes impact the effective regulatory functions of key transcription factors including PU.1, Runx1, Notch-RBPJ, and Bcl11b.

GATA-targeted compounds modulate cardiac subtype cell differentiation in dual reporter stem cell line

Background

Pharmacological modulation of cell fate decisions and developmental gene regulatory networks holds promise for the treatment of heart failure. Compounds that target tissue-specific transcription factors could overcome non-specific effects of small molecules and lead to the regeneration of heart muscle following myocardial infarction. Due to cellular heterogeneity in the heart, the activation of gene programs representing specific atrial and ventricular cardiomyocyte subtypes would be highly desirable. Chemical compounds that modulate atrial and ventricular cell fate could be used to improve subtype-specific differentiation of endogenous or exogenously delivered progenitor cells in order to promote cardiac regeneration.

Methods

Transcription factor GATA4-targeted compounds that have previously shown in vivo efficacy in cardiac injury models were tested for stage-specific activation of atrial and ventricular reporter genes in differentiating pluripotent stem cells using a dual reporter assay. Chemically induced gene expression changes were characterized by qRT-PCR, global run-on sequencing (GRO-seq) and immunoblotting, and the network of cooperative proteins of GATA4 and NKX2-5 were further explored by the examination of the GATA4 and NKX2-5 interactome by BioID. Reporter gene assays were conducted to examine combinatorial effects of GATA-targeted compounds and bromodomain and extraterminal domain (BET) inhibition on chamber-specific gene expression.

Results

GATA4-targeted compounds 3i-1000 and 3i-1103 were identified as differential modulators of atrial and ventricular gene expression. More detailed structure-function analysis revealed a distinct subclass of GATA4/NKX2-5 inhibitory compounds with an acetyl lysine-like domain that contributed to ventricular cells (%Myl2-eGFP+). Additionally, BioID analysis indicated broad interaction between GATA4 and BET family of proteins, such as BRD4. This indicated the involvement of epigenetic modulators in the regulation of GATA-dependent transcription. In this line, reporter gene assays with combinatorial treatment of 3i-1000 and the BET bromodomain inhibitor (+)-JQ1 demonstrated the cooperative role of GATA4 and BRD4 in the modulation of chamber-specific cardiac gene expression.

Conclusions

Collectively, these results indicate the potential for therapeutic alteration of cell fate decisions and pathological gene regulatory networks by GATA4-targeted compounds modulating chamber-specific transcriptional programs in multipotent cardiac progenitor cells and cardiomyocytes. The compound scaffolds described within this study could be used to develop regenerative strategies for myocardial regeneration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-021-02259-z.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

Enhancer grammar in development, evolution, and disease: dependencies and interplay

Each language has standard books describing that language's grammatical rules. Biologists have searched for similar, albeit more complex, principles relating enhancer sequence to gene expression. Here, we review the literature on enhancer grammar. We introduce dependency grammar, a model where enhancers encode information based on dependencies between enhancer features shaped by mechanistic, evolutionary, and biological constraints. Classifying enhancers based on the types of dependencies may identify unifying principles relating enhancer sequence to gene expression. Such rules would allow us to read the instructions for development within genomes and pinpoint causal enhancer variants underlying disease and evolutionary changes.

How transcription factors drive choice of the T cell fate

Recent evidence has elucidated how multipotent blood progenitors transform their identities in the thymus and undergo commitment to become T cells. Together with environmental signals, a core group of transcription factors have essential roles in this process by directly activating and repressing specific genes. Many of these transcription factors also function in later T cell development, but control different genes. Here, we review how these transcription factors work to change the activities of specific genomic loci during early intrathymic development to establish T cell lineage identity. We introduce the key regulators and highlight newly emergent insights into the rules that govern their actions. Whole-genome deep sequencing-based analysis has revealed unexpectedly rich relationships between inherited epigenetic states, transcription factor-DNA binding affinity thresholds and influences of given transcription factors on the activities of other factors in the same cells. Together, these mechanisms determine T cell identity and make the lineage choice irreversible.

Transcription factors regulate early T cell development via redeployment of other factors: Functional dynamics of constitutively required factors in cell fate decisions

Establishment of cell lineage identity from multipotent progenitors is controlled by cooperative actions of lineage-specific and stably expressed transcription factors, combined with input from environmental signals. Lineage-specific master transcription factors activate and repress gene expression by recruiting consistently expressed transcription factors and chromatin modifiers to their target loci. Recent technical advances in genome-wide and multi-omics analysis have shed light on unexpected mechanisms that underlie more complicated actions of transcription factors in cell fate decisions. In this review, we discuss functional dynamics of stably expressed and continuously required factors, Notch and Runx family members, throughout developmental stages of early T cell development in the thymus. Pre- and post-commitment stage-specific transcription factors induce dynamic redeployment of Notch and Runx binding genomic regions. Thus, together with stage-specific transcription factors, shared transcription factors across distinct developmental stages regulate acquisition of T lineage identity.

Modeling Human TBX5 Haploinsufficiency Predicts Regulatory Networks for Congenital Heart Disease

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.