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.

Valvular Interstitial Cells (VICs) are a common substrate for congenital and adult heart disease yet the signaling mechanisms governing their formation during early valvulogenesis are incompletely understood. We developed an unbiased strategy to identify genes important in endocardial epithelial-to-mesenchymal transformation (EMT) using a spatial transcriptional profile. Endocardial cells overlaying the cushions of the atrioventricular canal (AVC) and outflow tract (OFT) undergo an EMT to yield VICs. RNA sequencing (RNA-seq) analysis of gene expression between AVC, OFT, and ventricles (VEN) isolated from chick and mouse embryos at comparable stages of development (chick HH18; mouse E11.0) was performed. EMT occurs in the AVC and OFT cushions, but not VEN at this time. 198 genes in the chick (n=1) and 105 genes in the mouse (n=2) were enriched 2-fold in the cushions. Gene regulatory networks (GRN) generated from cushion-enriched gene lists confirmed TGFβ as a nodal point and identified NF-κB as a potential node. To reveal previously unrecognized regulators of EMT four candidate genes, Hapln1, Id1, Foxp2, and Meis2, and a candidate pathway, NF-κB, were selected. In vivo spatial expression of each gene was confirmed by in situ hybridization and a functional role for each in endocardial EMT was determined by siRNA knockdown in a collagen gel assay. Our spatial-transcriptional profiling strategy yielded gene lists which reflected the known biology of the system. Further analysis accurately identified and validated previously unrecognized novel candidate genes and the NF-κB pathway as regulators of endocardial cell EMT in vitro.

Mutations in the lamin A/C (LMNA) gene, which encodes nuclear membrane proteins, cause a variety of human conditions including dilated cardiomyopathy (DCM) with associated cardiac conduction system disease. To investigate mechanisms responsible for electrophysiologic and myocardial phenotypes caused by dominant human LMNA mutations, we performed longitudinal evaluations in heterozygous Lmna(+/-) mice. Despite one normal allele, Lmna(+/-) mice had 50% of normal cardiac lamin A/C levels and developed cardiac abnormalities. Conduction system function was normal in neonatal Lmna(+/-) mice but, by 4 weeks of age, atrioventricular (AV) nodal myocytes had abnormally shaped nuclei and active apoptosis. Telemetric and in vivo electrophysiologic studies in 10-week-old Lmna(+/-) mice showed AV conduction defects and both atrial and ventricular arrhythmias, analogous to those observed in humans with heterozygous LMNA mutations. Isolated myocytes from 12-month-old Lmna(+/-) mice exhibited impaired contractility. In vivo cardiac studies of aged Lmna(+/-) mice revealed DCM; in some mice this occurred without overt conduction system disease. However, neither histopathology nor serum CK levels indicated skeletal muscle pathology. These data demonstrate cardiac pathology due to heterozygous Lmna mutations reflecting a 50% reduction in lamin protein levels. Lamin haploinsufficiency caused early-onset programmed cell death of AV nodal myocytes and progressive electrophysiologic disease. While lamin haploinsufficiency was better tolerated by non-conducting myocytes, ultimately, these too succumbed to diminished lamin levels leading to dilated cardiomyopathy, which presumably arose independently from conduction system disease.

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.

Impulses generated by a multicellular, bioelectric signaling center termed the sinoatrial node (SAN) stimulate the rhythmic contraction of the heart. The SAN consists of a network of electrochemically oscillating pacemaker cells encased in a heterogeneous connective tissue microenvironment. Although the cellular composition of the SAN has been a point of interest for more than a century, the biological processes that drive the tissue-level assembly of the cells within the SAN are unknown. Here, we demonstrate that the SAN's structural features result from a developmental process during which mesenchymal cells derived from a multipotent progenitor structure, the proepicardium, integrate with and surround pacemaker myocardium. This process actively remodels the forming SAN and is necessary for sustained electrogenic signal generation and propagation. Collectively, these findings provide experimental evidence for how the microenvironmental architecture of the SAN is patterned and demonstrate that proper cellular arrangement is critical for cardiac pacemaker biorhythmicity.

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 contribution of somatic mosaicism, or genetic mutations arising after oocyte fertilization, to congenital heart disease (CHD) is not well understood. Further, the relationship between mosaicism in blood and cardiovascular tissue has not been determined.

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.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell–derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

ALPK3 encodes α-kinase 3, a muscle-specific protein of unknown function. ALPK3 loss-of-function variants cause cardiomyopathy with distinctive clinical manifestations in both children and adults, but the molecular functions of ALPK3 remain poorly understood.

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.