TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids

In the nascent mesoderm, TBXT expression must be precisely regulated to ensure that cells exit the primitive streak and pattern the anterior-posterior axis, but how varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dosage reduction during early human gastrulation using human induced pluripotent stem cell models of gastrulation and mesoderm differentiation. Multi-omic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprising wild-type, TBXT heterozygous or TBXT null human induced pluripotent stem cells reveal that varying TBXT dosage does not compromise the ability of a cell to differentiate into nascent mesoderm, but instead directly influences the temporal progression of the epithelial-to-mesenchymal transition with wild type transitioning first, followed by TBXT heterozygous and then TBXT null. By differentiating cells into nascent mesoderm in a monolayer format, we further illustrate that TBXT dosage directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that epithelial-to-mesenchymal transition progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

A disrupted compartment boundary underlies abnormal cardiac patterning and congenital heart defects

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.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development

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. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development

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.

TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids

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.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

A Mesp1-dependent developmental breakpoint in transcriptional and epigenomic specification of early cardiac precursors

Transcriptional networks governing cardiac precursor cell (CPC) specification are incompletely understood owing, in part, to limitations in distinguishing CPCs from non-cardiac mesoderm in early gastrulation. We leveraged detection of early cardiac lineage transgenes within a granular single-cell transcriptomic time course of mouse embryos to identify emerging CPCs and describe their transcriptional profiles. Mesp1, a transiently expressed mesodermal transcription factor, is canonically described as an early regulator of cardiac specification. However, we observed perdurance of CPC transgene-expressing cells in Mesp1 mutants, albeit mislocalized, prompting us to investigate the scope of the role of Mesp1 in CPC emergence and differentiation. Mesp1 mutant CPCs failed to robustly activate markers of cardiomyocyte maturity and crucial cardiac transcription factors, yet they exhibited transcriptional profiles resembling cardiac mesoderm progressing towards cardiomyocyte fates. Single-cell chromatin accessibility analysis defined a Mesp1-dependent developmental breakpoint in cardiac lineage progression at a shift from mesendoderm transcriptional networks to those necessary for cardiac patterning and morphogenesis. These results reveal Mesp1-independent aspects of early CPC specification and underscore a Mesp1-dependent regulatory landscape required for progression through cardiogenesis.

Graded mesoderm assembly governs cell fate and morphogenesis of the early mammalian heart

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.

Smarcd3 is an epigenetic modulator of the metabolic landscape in pancreatic ductal adenocarcinoma

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.

Transcription Factor GATA4 Regulates Cell Type-Specific Splicing Through Direct Interaction With RNA in Human Induced Pluripotent Stem Cell-Derived Cardiac Progenitors

BACKGROUND

GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors.

METHODS

We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4's novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors.

RESULTS

We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms.

CONCLUSIONS

This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

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

Motivation

Unsupervised clustering of single-cell transcriptomics is a powerful method for identifying cell populations. Static visualization techniques for single-cell clustering only display results for a single resolution parameter. Analysts will often evaluate more than one resolution parameter but then only report one.

Results

We developed Cell Layers, an interactive Sankey tool for the quantitative investigation of gene expression, co-expression, biological processes and cluster integrity across clustering resolutions. Cell Layers enhances the interpretability of single-cell clustering by linking molecular data and cluster evaluation metrics, providing novel insight into cell populations.

Availability and implementation

https://github.com/apblair/CellLayers.

Modeling congenital heart disease: lessons from mice, hPSC-based models, and organoids

In this review, Rao et al. contrast mouse models of CHD, which maintain the anatomical arrangement of the heart, and human cellular models of CHD, which are more likely to capture human-specific biology but lack anatomical structure. They also discuss the recent development of cardiac organoids, which are a promising step toward more anatomically informative human models of CHD.

Congenital heart defects (CHDs) are among the most common birth defects, but their etiology has long been mysterious. In recent decades, the development of a variety of experimental models has led to a greater understanding of the molecular basis of CHDs. In this review, we contrast mouse models of CHD, which maintain the anatomical arrangement of the heart, and human cellular models of CHD, which are more likely to capture human-specific biology but lack anatomical structure. We also discuss the recent development of cardiac organoids, which are a promising step toward more anatomically informative human models of CHD.

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.

Co-emergence of cardiac and gut tissues promotes cardiomyocyte maturation within human iPSC-derived organoids

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.

Dissecting CTCF site function in a tense HoxD locus

In this issue of Genes & Development, Amândio and colleagues (pp. 1490-1509) dissect the function of a cluster of several CTCF binding sites in the HoxD cluster by iterative deletions in mice. They found additive functions for some, and intriguingly found that some sites function as insulators, while others function as anchors for enhancer-promoter interactions. These functions vary depending on developmental context. The work provides new insights into the roles played by CTCF in regulating developmental patterns and 3D chromatin organization.

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency

[Figure: see text].

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.

WAPL maintains a cohesin loading cycle to preserve cell-type-specific distal gene regulation

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.

Molecular basis of CTCF binding polarity in genome folding

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.

Regulation of single-cell genome organization into TADs and chromatin nanodomains

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.

The developing heart: from The Wizard of Oz to congenital heart disease

The heart is an essential organ with a fascinating developmental biology. It is also one of the organs that is most often affected in human disease, either during development or in postnatal life. Over the last few decades, insights into the development of the heart have led to fundamental new concepts in gene regulation, but also to genetic and mechanistic insights into congenital heart defects. In more recent years, the lessons learned from studying heart development have been applied to interrogating regeneration of the diseased heart, exemplifying the importance of understanding the mechanistic underpinnings that lead to the development of an organ.

Author Correction: Transcriptional profiling and therapeutic targeting of oxidative stress in neuroinflammation

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Salt-inducible kinase 1 maintains HDAC7 stability to promote pathologic cardiac remodeling

Salt-inducible kinases (SIKs) are key regulators of cellular metabolism and growth, but their role in cardiomyocyte plasticity and heart failure pathogenesis remains unknown. Here, we showed that loss of SIK1 kinase activity protected against adverse cardiac remodeling and heart failure pathogenesis in rodent models and cardiomyocytes derived from human induced pluripotent stem cells. We found that SIK1 phosphorylated and stabilized histone deacetylase 7 (HDAC7) protein during cardiac stress, an event that is required for pathologic cardiomyocyte remodeling. Gain- and loss-of-function studies of HDAC7 in cultured cardiomyocytes implicated HDAC7 as a prohypertrophic signaling effector that can induce c-Myc expression, indicating a functional departure from the canonical MEF2 corepressor function of class IIa HDACs. Taken together, our findings reveal what we believe to be a previously unrecognized role for a SIK1/HDAC7 axis in regulating cardiac stress responses and implicate this pathway as a potential target in human heart failure.

Defining the relative and combined contribution of CTCF and CTCFL to genomic regulation

BACKGROUND

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.

RESULTS

We demonstrate that the N terminus of CTCF interacts with cohesin which explains the requirement for convergent CTCF binding sites in loop formation. By analyzing CTCF and CTCFL binding in tandem, we identify phenotypically distinct sites with respect to motifs, targeting to promoter/intronic intergenic regions and chromatin folding. Finally, we reveal that the N, C, and zinc finger terminal domains play unique roles in targeting each paralog to distinct binding sites to regulate transcription, chromatin looping, and insulation.

CONCLUSION

This study clarifies the unique and combined contribution of CTCF and CTCFL to chromosome organization and transcription, with direct implications for understanding how their co-expression deregulates transcription in cancer.

Transcriptional profiling and therapeutic targeting of oxidative stress in neuroinflammation

Oxidative stress is a central part of innate-immune induced neurodegeneration. However, the transcriptomic landscape of the central nervous system (CNS) innate immune cells contributing to oxidative stress is unknown, and therapies to target their neurotoxic functions are not widely available. Here, we provide the oxidative stress innate immune cell atlas in neuroinflammatory disease, and report the discovery of new druggable pathways. Transcriptional profiling of oxidative stress-producing CNS innate immune cells (Tox-seq) identified a core oxidative stress gene signature coupled to coagulation and glutathione pathway genes shared between a microglia cluster and infiltrating macrophages. Tox-seq followed by a microglia high-throughput screen (HTS) and oxidative stress gene network analysis, identified the glutathione regulating compound acivicin with potent therapeutic effects decreasing oxidative stress and axonal damage in chronic and relapsing multiple sclerosis (MS) models. Thus, oxidative stress transcriptomics identified neurotoxic CNS innate immune populations and may enable the discovery of selective neuroprotective strategies.

Minimal in vivo requirements for developmentally regulated cardiac long intergenic non-coding RNAs

Long intergenic non-coding RNAs (lincRNAs) have been implicated in gene regulation, but their requirement for development needs empirical interrogation. We computationally identified nine murine lincRNAs that have developmentally regulated transcriptional and epigenomic profiles specific to early heart differentiation. Six of the nine lincRNAs had in vivo expression patterns supporting a potential function in heart development, including a transcript downstream of the cardiac transcription factor Hand2, which we named Handlr (Hand2-associated lincRNA), Rubie and Atcayos We genetically ablated these six lincRNAs in mouse, which suggested genomic regulatory roles for four of the cohort. However, none of the lincRNA deletions led to severe cardiac phenotypes. Thus, we stressed the hearts of adult Handlr and Atcayos mutant mice by transverse aortic banding and found that absence of these lincRNAs did not affect cardiac hypertrophy or left ventricular function post-stress. Our results support roles for lincRNA transcripts and/or transcription in the regulation of topologically associated genes. However, the individual importance of developmentally specific lincRNAs is yet to be established. Their status as either gene-like entities or epigenetic components of the nucleus should be further considered.

CTCF confers local nucleosome resiliency after DNA replication and during mitosis

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.

RNA Interactions Are Essential for CTCF-Mediated Genome Organization

The function of the CCCTC-binding factor (CTCF) in the organization of the genome has become an important area of investigation, but the mechanisms by which CTCF dynamically contributes to genome organization are not clear. We previously discovered that CTCF binds to large numbers of endogenous RNAs, promoting its self-association. In this regard, we now report two independent features that disrupt CTCF association with chromatin: inhibition of transcription and disruption of CTCF-RNA interactions through mutations of 2 of its 11 zinc fingers that are not required for CTCF binding to its cognate DNA site: zinc finger 1 (ZF1) or zinc finger 10 (ZF10). These mutations alter gene expression profiles as CTCF mutants lose their ability to form chromatin loops and thus the ability to insulate chromatin domains and to mediate CTCF long-range genomic interactions. Our results point to the importance of CTCF-mediated RNA interactions as a structural component of genome organization.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

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.

Dynamic BAF chromatin remodeling complex subunit inclusion promotes temporally distinct gene expression programs in cardiogenesis

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.

Chromatin Domains Go on Repeat in Disease

TAD boundaries are insulators of genomic neighborhoods. In this issue, Sun et al. show that disease-associated tandem repeats are located to TAD boundaries and affect their insulation. The findings have important implications for TAD function and mechanisms underlying diseases such as fragile X syndrome and Huntington's disease.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication

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.

A De Novo Shape Motif Discovery Algorithm Reveals Preferences of Transcription Factors for DNA Shape Beyond Sequence Motifs

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.

Heart enhancers with deeply conserved regulatory activity are established early in zebrafish development

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.

During early embryogenesis, critical cardiac specification events occur. Here the authors isolate cardiac progenitor cells from early zebrafish embryos and characterize accessible chromatin regions specific to this cell population, finding that many of these regions overlap with conserved non-coding elements that are ortholgous to accessible chromatin regions in human.

A novel reporter allele for monitoring Dll4 expression within the embryonic and adult mouse

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.

Cardiac-enriched BAF chromatin-remodeling complex subunit Baf60c regulates gene expression programs essential for heart development and function

How chromatin-remodeling complexes modulate gene networks to control organ-specific properties is not well understood. For example, Baf60c (Smarcd3) encodes a cardiac-enriched subunit of the SWI/SNF-like BAF chromatin complex, but its role in heart development is not fully understood. We found that constitutive loss of Baf60c leads to embryonic cardiac hypoplasia and pronounced cardiac dysfunction. Conditional deletion of Baf60c in cardiomyocytes resulted in postnatal dilated cardiomyopathy with impaired contractile function. Baf60c regulates a gene expression program that includes genes encoding contractile proteins, modulators of sarcomere function, and cardiac metabolic genes. Many of the genes deregulated in Baf60c null embryos are targets of the MEF2/SRF co-factor Myocardin (MYOCD). In a yeast two-hybrid screen, we identified MYOCD as a BAF60c interacting factor; we showed that BAF60c and MYOCD directly and functionally interact. We conclude that Baf60c is essential for coordinating a program of gene expression that regulates the fundamental functional properties of cardiomyocytes.

Cooperative activation of cardiac transcription through myocardin bridging of paired MEF2 sites

Enhancers frequently contain multiple binding sites for the same transcription factor. These homotypic binding sites often exhibit synergy, whereby the transcriptional output from two or more binding sites is greater than the sum of the contributions of the individual binding sites alone. Although this phenomenon is frequently observed, the mechanistic basis for homotypic binding site synergy is poorly understood. Here, we identify a bona fide cardiac-specific Prkaa2 enhancer that is synergistically activated by homotypic MEF2 binding sites. We show that two MEF2 sites in the enhancer function cooperatively due to bridging of the MEF2C-bound sites by the SAP domain-containing co-activator protein myocardin, and we show that paired sites buffer the enhancer from integration site-dependent effects on transcription in vivo Paired MEF2 sites are prevalent in cardiac enhancers, suggesting that this might be a common mechanism underlying synergy in the control of cardiac gene expression in vivo.

ATP-dependent chromatin remodeling during mammalian development

Precise gene expression ensures proper stem and progenitor cell differentiation, lineage commitment and organogenesis during mammalian development. ATP-dependent chromatin-remodeling complexes utilize the energy from ATP hydrolysis to reorganize chromatin and, hence, regulate gene expression. These complexes contain diverse subunits that together provide a multitude of functions, from early embryogenesis through cell differentiation and development into various adult tissues. Here, we review the functions of chromatin remodelers and their different subunits during mammalian development. We discuss the mechanisms by which chromatin remodelers function and highlight their specificities during mammalian cell differentiation and organogenesis.

Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains from Genomic Compartmentalization

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.

Expandable Cardiovascular Progenitor Cells Reprogrammed from Fibroblasts

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.

Single-Cell Resolution of Temporal Gene Expression during Heart Development

Activation of complex molecular programs in specific cell lineages governs mammalian heart development, from a primordial linear tube to a four-chamber organ. To characterize lineage-specific, spatiotemporal developmental programs, we performed single-cell RNA sequencing of >1,200 murine cells isolated at seven time points spanning embryonic day 9.5 (primordial heart tube) to postnatal day 21 (mature heart). Using unbiased transcriptional data, we classified cardiomyocytes, endothelial cells, and fibroblast-enriched cells, thus identifying markers for temporal and chamber-specific developmental programs. By harnessing these datasets, we defined developmental ages of human and mouse pluripotent stem-cell-derived cardiomyocytes and characterized lineage-specific maturation defects in hearts of mice with heterozygous mutations in Nkx2.5 that cause human heart malformations. This spatiotemporal transcriptome analysis of heart development reveals lineage-specific gene programs underlying normal cardiac development and congenital heart disease.

KMT2D regulates specific programs in heart development via histone H3 lysine 4 di-methylation

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.

Highlighted article: Cardiac-specific depletion of the H3K4 methyltransferase KMT2D causes dysregulation of genes associated with cell cycle regulation, ion homeostasis and hypoxia signaling.

Loss of Iroquois homeobox transcription factors 3 and 5 in osteoblasts disrupts cranial mineralization

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.

Accelerated Evolution of Enhancer Hotspots in the Mammal Ancestor

Mammals have evolved remarkably different sensory, reproductive, metabolic, and skeletal systems. To explore the genetic basis for these differences, we developed a comparative genomics approach to scan whole-genome multiple sequence alignments to identify regions that evolved rapidly in an ancestral lineage but are conserved within extant species. This pattern suggests that ancestral changes in function were maintained in descendants. After applying this test to therian mammals, we identified 4,797 accelerated regions, many of which are noncoding and located near developmental transcription factors. We then used mouse transgenic reporter assays to test if noncoding accelerated regions are enhancers and to determine how therian-specific substitutions affect their activity in vivo. We discovered enhancers with expression specific to the therian version in brain regions involved in the hormonal control of milk ejection, uterine contractions, blood pressure, temperature, and visual processing. This work underscores the idea that changes in developmental gene expression are important for mammalian evolution, and it pinpoints candidate genes for unique aspects of mammalian biology.

Complex Interdependence Regulates Heterotypic Transcription Factor Distribution and Coordinates Cardiogenesis

Transcription factors (TFs) are thought to function with partners to achieve specificity and precise quantitative outputs. In the developing heart, heterotypic TF interactions, such as between the T-box TF TBX5 and the homeodomain TF NKX2-5, have been proposed as a mechanism for human congenital heart defects. We report extensive and complex interdependent genomic occupancy of TBX5, NKX2-5, and the zinc finger TF GATA4 coordinately controlling cardiac gene expression, differentiation, and morphogenesis. Interdependent binding serves not only to co-regulate gene expression but also to prevent TFs from distributing to ectopic loci and activate lineage-inappropriate genes. We define preferential motif arrangements for TBX5 and NKX2-5 cooperative binding sites, supported at the atomic level by their co-crystal structure bound to DNA, revealing a direct interaction between the two factors and induced DNA bending. Complex interdependent binding mechanisms reveal tightly regulated TF genomic distribution and define a combinatorial logic for heterotypic TF regulation of differentiation.

Human disease modeling reveals integrated transcriptional and epigenetic mechanisms of NOTCH1 haploinsufficiency

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.

Investigating the transcriptional control of cardiovascular development

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.

Brg1 modulates enhancer activation in mesoderm lineage commitment

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.

SUMMARY: The chromatin remodeling factor Brg1 is essential for mesoderm induction and, by modulating active and repressive chromatin states, is involved in promoting the activation of dynamic enhancers.