During embryogenesis, paracrine signaling between tissues in close proximity contributes to the determination of their respective cell fate(s) and development into functional organs. Organoids are in vitro models that mimic organ formation and cellular heterogeneity, but lack the paracrine input of surrounding tissues. Here, we describe a human multilineage iPSC-derived organoid that recapitulates cooperative cardiac and gut development and maturation, with extensive cellular and structural complexity in both tissues. We demonstrate that the presence of endoderm tissue (gut/intestine) in the organoids contributed to the development of cardiac tissue features characteristic of stages after heart tube formation, including cardiomyocyte expansion, compartmentalization, enrichment of atrial/nodal cells, myocardial compaction, and fetal-like functional maturation. Overall, this study demonstrates the ability to generate and mature cooperative tissues originating from different germ lineages within a single organoid model, an advance that will further the examination of multi-tissue interactions during development, physiological maturation, and disease.

Failure of septation of the interventricular septum (IVS) is the most common congenital heart defect (CHD), but mechanisms for patterning the IVS are largely unknown. We show that a Tbx5+/Mef2cAHF+ progenitor lineage forms a compartment boundary bisecting the IVS. This coordinated population originates at a first- and second heart field interface, subsequently forming a morphogenetic nexus. Ablation of Tbx5+/Mef2cAHF+ progenitors cause IVS disorganization, right ventricular hypoplasia and mixing of IVS lineages. Reduced dosage of the CHD transcription factor TBX5 disrupts boundary position and integrity, resulting in ventricular septation defects (VSDs) and patterning defects, including Slit2 and Ntn1 misexpression. Reducing NTN1 dosage partly rescues cardiac defects in Tbx5 mutant embryos. Loss of Slit2 or Ntn1 causes VSDs and perturbed septal lineage distributions. Thus, we identify essential cues that direct progenitors to pattern a compartment boundary for proper cardiac septation, revealing new mechanisms for cardiac birth defects.

During the phylotypic period, embryos from different genera show similar gene expression patterns, implying common regulatory mechanisms. Here we set out to identify enhancers involved in the initial events of cardiogenesis, which occurs during the phylotypic period. We isolate early cardiac progenitor cells from zebrafish embryos and characterize 3838 open chromatin regions specific to this cell population. Of these regions, 162 overlap with conserved non-coding elements (CNEs) that also map to open chromatin regions in human. Most of the zebrafish conserved open chromatin elements tested drive gene expression in the developing heart. Despite modest sequence identity, human orthologous open chromatin regions recapitulate the spatial temporal expression patterns of the zebrafish sequence, potentially providing a basis for phylotypic gene expression patterns. Genome-wide, we discover 5598 zebrafish-human conserved open chromatin regions, suggesting that a diverse repertoire of ancient enhancers is established prior to organogenesis and the phylotypic period.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

The mechanisms by which transcription factor haploinsufficiency alters the epigenetic and transcriptional landscape in human cells to cause disease are unknown. Here, we utilized human induced pluripotent stem cell (iPSC)-derived endothelial cells (ECs) to show that heterozygous nonsense mutations in NOTCH1 that cause aortic valve calcification disrupt the epigenetic architecture, resulting in derepression of latent pro-osteogenic and -inflammatory gene networks. Hemodynamic shear stress, which protects valves from calcification in vivo, activated anti-osteogenic and anti-inflammatory networks in NOTCH1(+/+), but not NOTCH1(+/-), iPSC-derived ECs. NOTCH1 haploinsufficiency altered H3K27ac at NOTCH1-bound enhancers, dysregulating downstream transcription of more than 1,000 genes involved in osteogenesis, inflammation, and oxidative stress. Computational predictions of the disrupted NOTCH1-dependent gene network revealed regulatory nodes that, when modulated, restored the network toward the NOTCH1(+/+) state. Our results highlight how alterations in transcription factor dosage affect gene networks leading to human disease and reveal nodes for potential therapeutic intervention.

Cranial malformations are a significant cause of perinatal morbidity and mortality. Iroquois homeobox transcription factors (IRX) are expressed early in bone tissue formation and facilitate patterning and mineralization of the skeleton. Mice lacking Irx5 appear grossly normal, suggesting that redundancy within the Iroquois family. However, global loss of both Irx3 and Irx5 in mice leads to significant skeletal malformations and embryonic lethality from cardiac defects. Here, we study the bone-specific functions of Irx3 and Irx5 using Osx-Cre to drive osteoblast lineage-specific deletion of Irx3 in Irx5(-/-) mice. Although we found that the Osx-Cre transgene alone could also affect craniofacial mineralization, newborn Irx3 (flox/flox) /Irx5(-/-)/Osx-Cre (+) mice displayed additional mineralization defects in parietal, interparietal, and frontal bones with enlarged sutures and reduced calvarial expression of osteogenic genes. Newborn endochondral long bones were largely unaffected, but we observed marked reductions in 3-4-week old bone mineral content of Irx3 (flox/flox) /Irx5(-/-)/Osx-Cre (+) mice. Our findings indicate that IRX3 and IRX5 can work together to regulate mineralization of specific cranial bones. Our results also provide insight into the causes of the skeletal changes and mineralization defects seen in Hamamy syndrome patients carrying mutations in IRX5.

Large-scale quantitative analysis of transcriptional co-expression has been used to dissect regulatory networks and to predict the functions of new genes discovered by genome sequencing in model organisms such as yeast. Although the idea that tissue-specific expression is indicative of gene function in mammals is widely accepted, it has not been objectively tested nor compared with the related but distinct strategy of correlating gene co-expression as a means to predict gene function.

Ubiquitously expressed CTCF is involved in numerous cellular functions, such as organizing chromatin into TAD structures. In contrast, its paralog, CTCFL, is normally only present in the testis. However, it is also aberrantly expressed in many cancers. While it is known that shared and unique zinc finger sequences in CTCF and CTCFL enable CTCFL to bind competitively to a subset of CTCF binding sites as well as its own unique locations, the impact of CTCFL on chromosome organization and gene expression has not been comprehensively analyzed in the context of CTCF function. Using an inducible complementation system, we analyze the impact of expressing CTCFL and CTCF-CTCFL chimeric proteins in the presence or absence of endogenous CTCF to clarify the relative and combined contribution of CTCF and CTCFL to chromosome organization and transcription.

Heart failure with reduced ejection fraction (HFrEF) is associated with high mortality, highlighting an urgent need for new therapeutic strategies. As stress-activated cardiac signaling cascades converge on the nucleus to drive maladaptive gene programs, interdicting pathological transcription is a conceptually attractive approach for HFrEF therapy. Here, we demonstrate that CDK7/12/13 are critical regulators of transcription activation in the heart that can be pharmacologically inhibited to improve HFrEF. CDK7/12/13 inhibition using the first-in-class inhibitor THZ1 or RNAi blocks stress-induced transcription and pathologic hypertrophy in cultured rodent cardiomyocytes. THZ1 potently attenuates adverse cardiac remodeling and HFrEF pathogenesis in mice and blocks cardinal features of disease in human iPSC-derived cardiomyocytes. THZ1 suppresses Pol II enrichment at stress-transactivated cardiac genes and inhibits a specific pathologic gene program in the failing mouse heart. These data identify CDK7/12/13 as druggable regulators of cardiac gene transactivation during disease-related stress, suggesting that HFrEF features a critical dependency on transcription that can be therapeutically exploited.

Using four-dimensional whole-embryo light sheet imaging with improved and accessible computational tools, we longitudinally reconstruct early murine cardiac development at single-cell resolution. Nascent mesoderm progenitors form opposing density and motility gradients, converting the temporal birth sequence of gastrulation into a spatial anterolateral-to-posteromedial arrangement. Migrating precardiac mesoderm does not strictly preserve cellular neighbor relationships, and spatial patterns only become solidified as the cardiac crescent emerges. Progenitors undergo a mesenchymal-to-epithelial transition, with a first heart field (FHF) ridge apposing a motile juxta-cardiac field (JCF). Anchored along the ridge, the FHF epithelium rotates the JCF forward to form the initial heart tube, along with push-pull morphodynamics of the second heart field. In Mesp1 mutants that fail to make a cardiac crescent, mesoderm remains highly motile but directionally incoherent, resulting in density gradient inversion. Our practicable live embryo imaging approach defines spatial origins and behaviors of cardiac progenitors and identifies their unanticipated morphological transitions.

Pancreatic cancer is characterized by extensive resistance to conventional therapies, making clinical management a challenge. Here we map the epigenetic dependencies of cancer stem cells, cells that preferentially evade therapy and drive progression, and identify SWI/SNF complex member SMARCD3 as a regulator of pancreatic cancer cells. Although SWI/SNF subunits often act as tumor suppressors, we show that SMARCD3 is amplified in cancer, enriched in pancreatic cancer stem cells and upregulated in the human disease. Diverse genetic mouse models of pancreatic cancer and stage-specific Smarcd3 deletion reveal that Smarcd3 loss preferentially impacts established tumors, improving survival especially in context of chemotherapy. Mechanistically, SMARCD3 acts with FOXA1 to control lipid and fatty acid metabolism, programs associated with therapy resistance and poor prognosis in cancer. These data identify SMARCD3 as an epigenetic modulator responsible for establishing the metabolic landscape in aggressive pancreatic cancer cells and a potential target for new therapies.

Heart development is exquisitely sensitive to the precise temporal regulation of thousands of genes that govern developmental decisions during differentiation. However, we currently lack a detailed understanding of how chromatin and gene expression patterns are coordinated during developmental transitions in the cardiac lineage. Here, we interrogated the transcriptome and several histone modifications across the genome during defined stages of cardiac differentiation. We find distinct chromatin patterns that are coordinated with stage-specific expression of functionally related genes, including many human disease-associated genes. Moreover, we discover a novel preactivation chromatin pattern at the promoters of genes associated with heart development and cardiac function. We further identify stage-specific distal enhancer elements and find enriched DNA binding motifs within these regions that predict sets of transcription factors that orchestrate cardiac differentiation. Together, these findings form a basis for understanding developmentally regulated chromatin transitions during lineage commitment and the molecular etiology of congenital heart disease.

Lysine acetylation regulates transcription by targeting histones and nonhistone proteins. Here we report that the central regulator of transcription, RNA polymerase II, is subject to acetylation in mammalian cells. Acetylation occurs at eight lysines within the C-terminal domain (CTD) of the largest polymerase subunit and is mediated by p300/KAT3B. CTD acetylation is specifically enriched downstream of the transcription start sites of polymerase-occupied genes genome-wide, indicating a role in early stages of transcription initiation or elongation. Mutation of lysines or p300 inhibitor treatment causes the loss of epidermal growth-factor-induced expression of c-Fos and Egr2, immediate-early genes with promoter-proximally paused polymerases, but does not affect expression or polymerase occupancy at housekeeping genes. Our studies identify acetylation as a new modification of the mammalian RNA polymerase II required for the induction of growth factor response genes.

Direct reprogramming of adult somatic cells into alternative cell types has been shown for several lineages. We previously showed that GATA4, MEF2C, and TBX5 (GMT) directly reprogrammed nonmyocyte mouse heart cells into induced cardiomyocyte-like cells (iCMs) in vitro and in vivo. However, GMT alone appears insufficient in human fibroblasts, at least in vitro. Here, we show that GMT plus ESRRG and MESP1 induced global cardiac gene-expression and phenotypic shifts in human fibroblasts derived from embryonic stem cells, fetal heart, and neonatal skin. Adding Myocardin and ZFPM2 enhanced reprogramming, including sarcomere formation, calcium transients, and action potentials, although the efficiency remained low. Human iCM reprogramming was epigenetically stable. Furthermore, we found that transforming growth factor β signaling was important for, and improved the efficiency of, human iCM reprogramming. These findings demonstrate that human fibroblasts can be directly reprogrammed toward the cardiac lineage, and lay the foundation for future refinements in vitro and in vivo.

KMT2D, which encodes a histone H3K4 methyltransferase, has been implicated in human congenital heart disease in the context of Kabuki syndrome. However, its role in heart development is not understood. Here, we demonstrate a requirement for KMT2D in cardiac precursors and cardiomyocytes during cardiogenesis in mice. Gene expression analysis revealed downregulation of ion transport and cell cycle genes, leading to altered calcium handling and cell cycle defects. We further determined that myocardial Kmt2d deletion led to decreased H3K4me1 and H3K4me2 at enhancers and promoters. Finally, we identified KMT2D-bound regions in cardiomyocytes, of which a subset was associated with decreased gene expression and decreased H3K4me2 in mutant hearts. This subset included genes related to ion transport, hypoxia-reoxygenation and cell cycle regulation, suggesting that KMT2D is important for these processes. Our findings indicate that KMT2D is essential for regulating cardiac gene expression during heart development primarily via H3K4 di-methylation.

The interplay between different levels of gene regulation in modulating developmental transcriptional programs, such as histone modifications and chromatin remodeling, is not well understood. Here, we show that the chromatin remodeling factor Brg1 is required for enhancer activation in mesoderm induction. In an embryonic stem cell-based directed differentiation assay, the absence of Brg1 results in a failure of cardiomyocyte differentiation and broad deregulation of lineage-specific gene expression during mesoderm induction. We find that Brg1 co-localizes with H3K27ac at distal enhancers and is required for robust H3K27 acetylation at distal enhancers that are activated during mesoderm induction. Brg1 is also required to maintain Polycomb-mediated repression of non-mesodermal developmental regulators, suggesting cooperativity between Brg1 and Polycomb complexes. Thus, Brg1 is essential for modulating active and repressive chromatin states during mesoderm lineage commitment, in particular the activation of developmentally important enhancers. These findings demonstrate interplay between chromatin remodeling complexes and histone modifications that, together, ensure robust and broad gene regulation during crucial lineage commitment decisions.

Maintenance of vascular integrity is required for embryogenesis and organ homeostasis. However, the gene expression programs that stabilize blood vessels are poorly understood. Here, we show that the histone methyltransferase Ezh2 maintains integrity of the developing vasculature by repressing a transcriptional program that activates expression of Mmp9. Inactivation of Ezh2 in developing mouse endothelium caused embryonic lethality with compromised vascular integrity and increased extracellular matrix degradation. Genome-wide approaches showed that Ezh2 targets Mmp9 and its activators Fosl1 and Klf5. In addition, we uncovered Creb3l1 as an Ezh2 target that directly activates Mmp9 gene expression in the endothelium. Furthermore, genetic inactivation of Mmp9 rescued vascular integrity defects in Ezh2-deficient embryos. Thus, epigenetic repression of Creb3l1, Fosl1, Klf5 and Mmp9 by Ezh2 in endothelial cells maintains the integrity of the developing vasculature, potentially linking this transcriptional network to diseases with compromised vascular integrity.

Tbx20 is a member of the T-box transcription factor family expressed in the forming hearts of vertebrate and invertebrate embryos. We report here analysis of Tbx20 expression during murine cardiac development and assessment of DNA-binding and transcriptional properties of Tbx20 isoforms. Tbx20 was expressed in myocardium and endocardium, including high levels in endocardial cushions. cDNAs generated by alternative splicing encode at least four Tbx20 isoforms, and Tbx20a uniquely carried strong transactivation and transrepression domains in its C terminus. Isoforms with an intact T-box bound specifically to DNA sites resembling the consensus brachyury half site, although with less avidity compared with the related factor, Tbx5. Tbx20 physically interacted with cardiac transcription factors Nkx2-5, GATA4, and GATA5, collaborating to synergistically activate cardiac gene expression. Among cardiac GATA factors, there was preferential synergy with GATA5, implicated in endocardial differentiation. In Xenopus embryos, enforced expression of Tbx20a, but not Tbx20b, led to induction of mesodermal and endodermal lineage markers as well as cell migration, indicating that the long Tbx20a isoform uniquely bears functional domains that can alter gene expression and developmental behaviour in an in vivo context. We propose that Tbx20 plays an integrated role in the ancient myogenic program of the heart, and has been additionally coopted during evolution of vertebrates for endocardial cushion development.

In the nascent mesoderm, levels of Brachyury (TBXT) expression must be precisely regulated to ensure cells exit the primitive streak and pattern the anterior-posterior axis, but how this varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dose reduction during early human gastrulation using human induced pluripotent stem cell (hiPSC)-based models of gastrulation and mesoderm differentiation. Multiomic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprised of WT, TBXT heterozygous (TBXT-Het), or TBXT null (TBXT-KO) hiPSCs reveal that varying TBXT dosage does not compromise a cell's ability to differentiate into nascent mesoderm, but that the loss of TBXT significantly delays the temporal progression of the epithelial to mesenchymal transition (EMT). This delay is dependent on TBXT dose, as cells heterozygous for TBXT proceed with EMT at an intermediate pace relative to WT or TBXT-KO. By differentiating iPSCs of the allelic series into nascent mesoderm in a monolayer format, we further illustrate that TBXT dose directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that EMT progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

Transcriptional regulation of thousands of genes instructs complex morphogenetic and molecular events for heart development. Cardiac transcription factors choreograph gene expression at each stage of differentiation by interacting with cofactors, including chromatin-modifying enzymes, and by binding to a constellation of regulatory DNA elements. Here, we present salient examples relevant to cardiovascular development and heart disease, and review techniques that can sharpen our understanding of cardiovascular biology. We discuss the interplay between cardiac transcription factors, cis-regulatory elements, and chromatin as dynamic regulatory networks, to orchestrate sequential deployment of the cardiac gene expression program.