Heart formation requires a highly balanced network of transcriptional activation of genes. The homeodomain transcription factor, Shox2, is essential for the formation of the sinoatrial valves and for the development of the pacemaking system. The elucidation of molecular mechanisms underlying the development of pacemaker tissue has gained clinical interest as defects in its patterning can be related to atrial arrhythmias. We have analyzed putative targets of Shox2 and identified the Bmp4 gene as a direct target. Shox2 interacts directly with the Bmp4 promoter in chromatin immunoprecipitation assays and activates transcription in luciferase-reporter assays. In addition, ectopic expression of Shox2 in Xenopus embryos stimulates transcription of the Bmp4 gene, and silencing of Shox2 in cardiomyocytes leads to a reduction in the expression of Bmp4. In Tbx5(del/+) mice, a model for Holt-Oram syndrome, and Shox2(-/-) mice, we show that the T-box transcription factor Tbx5 is a regulator of Shox2 expression in the inflow tract and that Bmp4 is regulated by Shox2 in this compartment of the embryonic heart. In addition, we could show that Tbx5 acts cooperatively with Nkx2.5 to regulate the expression of Shox2 and Bmp4. This work establishes a link between Tbx5, Shox2 and Bmp4 in the pacemaker region of the developing heart and thus contributes to the unraveling of the intricate interplay between the heart-specific transcriptional machinery and developmental signaling pathways.

RPB1, the largest subunit of RNA polymerase II, contains a highly modifiable C-terminal domain (CTD) that consists of variations of a consensus heptad repeat sequence (Y1S2P3T4S5P6S7). The consensus CTD repeat motif and tandem organization represent the ancestral state of eukaryotic RPB1, but across eukaryotes CTDs show considerable diversity in repeat organization and sequence content. These differences may reflect lineage-specific CTD functions mediated by protein interactions. Mammalian CTDs contain eight non-consensus repeats with a lysine in the seventh position (K7). Posttranslational acetylation of these sites was recently shown to be required for proper polymerase pausing and regulation of two growth factor-regulated genes.

Fibro-adipogenic progenitors (FAPs) are important components of the skeletal muscle regenerative environment. Whether FAPs support muscle regeneration or promote fibro-adipogenic degeneration is emerging as a key determinant in the pathogenesis of muscular diseases, including Duchenne muscular dystrophy (DMD). However, the molecular mechanism that controls FAP lineage commitment and activity is currently unknown. We show here that an HDAC-myomiR-BAF60 variant network regulates the fate of FAPs in dystrophic muscles of mdx mice. Combinatorial analysis of gene expression microarray, genome-wide chromatin remodeling by nuclease accessibility (NA) combined with next-generation sequencing (NA-seq), small RNA sequencing (RNA-seq), and microRNA (miR) high-throughput screening (HTS) against SWI/SNF BAF60 variants revealed that HDAC inhibitors (HDACis) derepress a "latent" myogenic program in FAPs from dystrophic muscles at early stages of disease. Specifically, HDAC inhibition induces two core components of the myogenic transcriptional machinery, MYOD and BAF60C, and up-regulates the myogenic miRs (myomiRs) (miR-1.2, miR-133, and miR-206), which target the alternative BAF60 variants BAF60A and BAF60B, ultimately directing promyogenic differentiation while suppressing the fibro-adipogenic phenotype. In contrast, FAPs from late stage dystrophic muscles are resistant to HDACi-induced chromatin remodeling at myogenic loci and fail to activate the promyogenic phenotype. These results reveal a previously unappreciated disease stage-specific bipotency of mesenchimal cells within the regenerative environment of dystrophic muscles. Resolution of such bipotency by epigenetic intervention with HDACis provides a molecular rationale for the in situ reprogramming of target cells to promote therapeutic regeneration of dystrophic muscles.

The accurate and comprehensive identification of functional regulatory sequences in mammalian genomes remains a major challenge. Here we describe site-specific integration fluorescence-activated cell sorting followed by sequencing (SIF-seq), an unbiased, medium-throughput functional assay for the discovery of distant-acting enhancers. Targeted single-copy genomic integration into pluripotent cells, reporter assays and flow cytometry are coupled with high-throughput DNA sequencing to enable parallel screening of large numbers of DNA sequences. By functionally interrogating >500 kilobases (kb) of mouse and human sequence in mouse embryonic stem cells for enhancer activity we identified enhancers at pluripotency loci including NANOG. In in vitro-differentiated cardiomyocytes and neural progenitor cells, we identified cardiac enhancers and neuronal enhancers, respectively. SIF-seq is a powerful and flexible method for de novo functional identification of mammalian enhancers in a potentially wide variety of cell types.

DNA shape adds specificity to sequence motifs but has not been explored systematically outside this context. We hypothesized that DNA-binding proteins (DBPs) preferentially occupy DNA with specific structures ("shape motifs") regardless of whether or not these correspond to high information content sequence motifs. We present ShapeMF, a Gibbs sampling algorithm that identifies de novo shape motifs. Using binding data from hundreds of in vivo and in vitro experiments, we show that most DBPs have shape motifs and can occupy these in the absence of sequence motifs. This "shape-only binding" is common for many DBPs and in regions co-bound by multiple DBPs. When shape and sequence motifs co-occur, they can be overlapping, flanking, or separated by consistent spacing. Finally, DBPs within the same protein family have different shape motifs, explaining their distinct genome-wide occupancy despite having similar sequence motifs. These results suggest that shape motifs not only complement sequence motifs but also facilitate recognition of DNA beyond conventionally defined sequence motifs.

Chromatin remodeling complexes instruct cellular differentiation and lineage specific transcription. The BRG1/BRM-associated factor (BAF) complexes are important for several aspects of differentiation. We show that the catalytic subunit gene Brg1 has a specific role in cardiac precursors (CPs) to initiate cardiac gene expression programs and repress non-cardiac expression. Using immunopurification with mass spectrometry, we have determined the dynamic composition of BAF complexes during mammalian cardiac differentiation, identifying several cell-type specific subunits. We focused on the CP- and cardiomyocyte (CM)-enriched subunits BAF60c (SMARCD3) and BAF170 (SMARCC2). Baf60c and Baf170 co-regulate gene expression with Brg1 in CPs, and in CMs their loss results in broadly deregulated cardiac gene expression. BRG1, BAF60c and BAF170 modulate chromatin accessibility, to promote accessibility at activated genes while closing chromatin at repressed genes. BAF60c and BAF170 are required for proper BAF complex composition, and BAF170 loss leads to retention of BRG1 at CP-specific sites. Thus, dynamic interdependent BAF complex subunit assembly modulates chromatin states and thereby participates in directing temporal gene expression programs in cardiogenesis.

Canonical Notch signaling requires the presence of a membrane bound ligand and a corresponding transmembrane Notch receptor. Receptor engagement induces multiple proteolytic cleavage events culminating in the nuclear accumulation of the Notch intracellular domain and its binding to a transcriptional co-factor to mediate gene expression. Notch signaling networks are essential regulators of vascular patterning and angiogenesis, as well as myriad other biological processes. Delta-like 4 (Dll4) encodes the earliest Notch ligand detected in arterial cells, and is enriched in sprouting endothelial tip cells. Dll4 expression has often been inferred by proxy using a lacZ knockin reporter allele. This is problematic, as a single copy of Dll4 is haploinsufficient. Additionally, Notch activity regulates Dll4 transcription, making it unclear whether these reporter lines accurately reflect Dll4 expression. Accordingly, precisely defining Dll4 expression is essential for determining its role in development and disease. To address these limitations, we generated a novel BAC transgenic allele with a nuclear-localized β-galactosidase reporter (Dll4-BAC-nlacZ). Through a comparative analysis, we show the BAC line overcomes previous issues of haploinsufficiency, it recapitulates Dll4 expression in vivo, and allows superior visualization and imaging. As such, this novel Dll4 reporter is an important addition to the growing Notch toolkit.

The molecular mechanisms underlying folding of mammalian chromosomes remain poorly understood. The transcription factor CTCF is a candidate regulator of chromosomal structure. Using the auxin-inducible degron system in mouse embryonic stem cells, we show that CTCF is absolutely and dose-dependently required for looping between CTCF target sites and insulation of topologically associating domains (TADs). Restoring CTCF reinstates proper architecture on altered chromosomes, indicating a powerful instructive function for CTCF in chromatin folding. CTCF remains essential for TAD organization in non-dividing cells. Surprisingly, active and inactive genome compartments remain properly segregated upon CTCF depletion, revealing that compartmentalization of mammalian chromosomes emerges independently of proper insulation of TADs. Furthermore, our data support that CTCF mediates transcriptional insulator function through enhancer blocking but not as a direct barrier to heterochromatin spreading. Beyond defining the functions of CTCF in chromosome folding, these results provide new fundamental insights into the rules governing mammalian genome organization.

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

The cohesin complex has an essential role in maintaining genome organization. However, its role in gene regulation remains largely unresolved. Here we report that the cohesin release factor WAPL creates a pool of free cohesin, in a process known as cohesin turnover, which reloads it to cell-type-specific binding sites. Paradoxically, stabilization of cohesin binding, following WAPL ablation, results in depletion of cohesin from these cell-type-specific regions, loss of gene expression and differentiation. Chromosome conformation capture experiments show that cohesin turnover is important for maintaining promoter-enhancer loops. Binding of cohesin to cell-type-specific sites is dependent on the pioneer transcription factors OCT4 (POU5F1) and SOX2, but not NANOG. We show the importance of cohesin turnover in controlling transcription and propose that a cycle of cohesin loading and off-loading, instead of static cohesin binding, mediates promoter and enhancer interactions critical for gene regulation.

The genome folds into a hierarchy of three-dimensional structures within the nucleus. At the sub-megabase scale, chromosomes form topologically associating domains (TADs)1-4. However, how TADs fold in single cells is elusive. Here, we reveal TAD features inaccessible to cell population analysis by using super-resolution microscopy. TAD structures and physical insulation associated with their borders are variable between individual cells, yet chromatin intermingling is enriched within TADs compared to adjacent TADs in most cells. The spatial segregation of TADs is further exacerbated during cell differentiation. Favored interactions within TADs are regulated by cohesin and CTCF through distinct mechanisms: cohesin generates chromatin contacts and intermingling while CTCF prevents inter-TAD contacts. Furthermore, TADs are subdivided into discrete nanodomains, which persist in cells depleted of CTCF or cohesin, whereas disruption of nucleosome contacts alters their structural organization. Altogether, these results provide a physical basis for the folding of individual chromosomes at the nanoscale.

Proper regulation of gene dosage is critical for the development of the early embryo and the extraembryonic tissues that support it. Specifically, loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. In this study, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. We generate an allelic series for CDX2 in human induced pluripotent stem cells (hiPSCs) consisting of WT, heterozygous, and homozygous null CDX2 genotypes, differentiate these cells in a 2D gastruloid model, and subject these cells to multiomic single nucleus RNA and ATAC sequencing. We identify several genes that CDX2 dose-dependently regulate cytoskeletal integrity and adhesiveness in the extraembryonic mesoderm population, including regulators of the VEGF, canonical WNT, and non-canonical WNT signaling pathways. Despite these dose-dependent gene expression patterns, snATAC-seq reveals that heterozygous CDX2 expression is capable of inducing a WT-like chromatin accessibility profile, suggesting accessibility is not sufficient to drive gene expression when the CDX2 dosage is reduced. Finally, because the loss of CDX2 or TBXT phenocopy one another in vivo, we compare differentially expressed genes in our CDX2 knock-out model to those from TBXT knock-out hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability, including ANK3 and ANGPT1. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and suggest these genes may underlie the defects in vascular development and allantoic elongation seen in the absence or reduction of CDX2 in vivo.

Stem cell-based approaches to cardiac regeneration are increasingly viable strategies for treating heart failure. Generating abundant and functional autologous cells for transplantation in such a setting, however, remains a significant challenge. Here, we isolated a cell population with extensive proliferation capacity and restricted cardiovascular differentiation potentials during cardiac transdifferentiation of mouse fibroblasts. These induced expandable cardiovascular progenitor cells (ieCPCs) proliferated extensively for more than 18 passages in chemically defined conditions, with 10(5) starting fibroblasts robustly producing 10(16) ieCPCs. ieCPCs expressed cardiac signature genes and readily differentiated into functional cardiomyocytes (CMs), endothelial cells (ECs), and smooth muscle cells (SMCs) in vitro, even after long-term expansion. When transplanted into mouse hearts following myocardial infarction, ieCPCs spontaneously differentiated into CMs, ECs, and SMCs and improved cardiac function for up to 12 weeks after transplantation. Thus, ieCPCs are a powerful system to study cardiovascular specification and provide strategies for regenerative medicine in the heart.

The molecular pathway that controls cardiogenesis is temporally and spatially regulated by master transcriptional regulators such as NKX2-5, Isl1, MEF2C, GATA4, and beta-catenin. The interplay between these factors and their downstream targets are not completely understood. Here, we studied regulation of beta-catenin and GATA4 by NKX2-5 in human fetal cardiac myocytes.

We describe the genome of the western painted turtle, Chrysemys picta bellii, one of the most widespread, abundant, and well-studied turtles. We place the genome into a comparative evolutionary context, and focus on genomic features associated with tooth loss, immune function, longevity, sex differentiation and determination, and the species' physiological capacities to withstand extreme anoxia and tissue freezing.

Unraveling the gene regulatory networks that govern development and function of the mammalian heart is critical for the rational design of therapeutic interventions in human heart disease. Using the Drosophila heart as a platform for identifying novel gene interactions leading to heart disease, we found that the Rho-GTPase Cdc42 cooperates with the cardiac transcription factor Tinman/Nkx2-5. Compound Cdc42, tinman heterozygous mutant flies exhibited impaired cardiac output and altered myofibrillar architecture, and adult heart-specific interference with Cdc42 function is sufficient to cause these same defects. We also identified K(+) channels, encoded by dSUR and slowpoke, as potential effectors of the Cdc42-Tinman interaction. To determine whether a Cdc42-Nkx2-5 interaction is conserved in the mammalian heart, we examined compound heterozygous mutant mice and found conduction system and cardiac output defects. In exploring the mechanism of Nkx2-5 interaction with Cdc42, we demonstrated that mouse Cdc42 was a target of, and negatively regulated by miR-1, which itself was negatively regulated by Nkx2-5 in the mouse heart and by Tinman in the fly heart. We conclude that Cdc42 plays a conserved role in regulating heart function and is an indirect target of Tinman/Nkx2-5 via miR-1.

Adult-onset diseases can be associated with in utero events, but mechanisms for this remain unknown(1,2). The Polycomb histone methyltransferase Ezh2 stabilizes transcription by depositing repressive marks during development that persist into adulthood(3-9), but its function in postnatal organ homeostasis is unknown. We show that Ezh2 stabilizes cardiac gene expression and prevents cardiac pathology by repressing the homeodomain transcription factor gene Six1, which functions in cardiac progenitor cells but is stably silenced upon cardiac differentiation. Deletion of Ezh2 in cardiac progenitors caused postnatal myocardial pathology and destabilized cardiac gene expression with activation of Six1-dependent skeletal muscle genes. Six1 induced cardiomyocyte hypertrophy and skeletal muscle gene expression. Furthermore, genetically reducing Six1 levels rescued the pathology of Ezh2-deficient hearts. Thus, Ezh2-mediated repression of Six1 in differentiating cardiac progenitors is essential for stable gene expression and homeostasis in the postnatal heart. Our results suggest that epigenetic dysregulation in embryonic progenitor cells is a predisposing factor for adult disease and dysregulated stress responses.

Dominant mutations in the T-box transcription factor gene TBX5 cause Holt-Oram syndrome (HOS), an inherited human disease characterized by upper limb malformations and congenital heart defects (CHDs) of variable severity. We hypothesize that minor alterations in the dosage of Tbx5 directly influences severity of CHDs. Using a mouse allelic series, we show a sensitive inverse correlation between Tbx5 dosage and abnormal cardiac morphogenesis and gene expression. The CHDs found in mice harbouring a hypomorphic allele of Tbx5 (Tbx5(lox/+) mice) are less pronounced than those found in Tbx5 haploinsufficient mice (Tbx5(del/+)), and homozygous hypomorphic (Tbx5(lox/lox)) embryos have noticeably more advanced cardiac development than Tbx5 null (Tbx5(del/del)) embryos. Examination of target gene expression across the allelic series uncovers very fine sensitivity across the range of Tbx5 dosages, in which some genes respond dramatically differently to only 15% differences in Tbx5 mRNA levels. This analysis was expanded to a genome-wide level, which uncovered a Tbx5 dosage-sensitive genetic program involving a network of cardiac transcription factors, developmentally important cell-cell signaling molecules, and ion channel proteins. These results indicate an exquisite sensitivity of the developing heart to Tbx5 dosage and provide significant insight into the transcriptional and cellular mechanisms that are disrupted in CHDs.

The access of Transcription Factors (TFs) to their cognate DNA binding motifs requires a precise control over nucleosome positioning. This is especially important following DNA replication and during mitosis, both resulting in profound changes in nucleosome organization over TF binding regions. Using mouse Embryonic Stem (ES) cells, we show that the TF CTCF displaces nucleosomes from its binding site and locally organizes large and phased nucleosomal arrays, not only in interphase steady-state but also immediately after replication and during mitosis. Correlative analyses suggest this is associated with fast gene reactivation following replication and mitosis. While regions bound by other TFs (Oct4/Sox2), display major rearrangement, the post-replication and mitotic nucleosome positioning activity of CTCF is not unique: Esrrb binding regions are also characterized by persistent nucleosome positioning. Therefore, selected TFs such as CTCF and Esrrb act as resilient TFs governing the inheritance of nucleosome positioning at regulatory regions throughout the cell-cycle.

Current models propose that boundaries of mammalian topologically associating domains (TADs) arise from the ability of the CTCF protein to stop extrusion of chromatin loops by cohesin. While the orientation of CTCF motifs determines which pairs of CTCF sites preferentially stabilize loops, the molecular basis of this polarity remains unclear. By combining ChIP-seq and single molecule live imaging we report that CTCF positions cohesin, but does not control its overall binding dynamics on chromatin. Using an inducible complementation system, we find that CTCF mutants lacking the N-terminus cannot insulate TADs properly. Cohesin remains at CTCF sites in this mutant, albeit with reduced enrichment. Given the orientation of CTCF motifs presents the N-terminus towards cohesin as it translocates from the interior of TADs, these observations explain how the orientation of CTCF binding sites translates into genome folding patterns.