Cell type-specific chromatin immunoprecipitation from multicellular complex samples using BiTS-ChIP

Enhancer-promoter interactions become more instructive in the transition from cell-fate specification to tissue differentiation

To regulate expression, enhancers must come in proximity to their target gene. However, the relationship between the timing of enhancer-promoter (E-P) proximity and activity remains unclear, with examples of uncoupled, anticorrelated and correlated interactions. To assess this, we selected 600 characterized enhancers or promoters with tissue-specific activity in Drosophila embryos and performed Capture-C in FACS-purified myogenic or neurogenic cells during specification and tissue differentiation. This enabled direct comparison between E-P proximity and activity transitioning from OFF-to-ON and ON-to-OFF states across developmental conditions. This showed remarkably similar E-P topologies between specified muscle and neuronal cells, which are uncoupled from activity. During tissue differentiation, many new distal interactions emerge where changes in E-P proximity reflect changes in activity. The mode of E-P regulation therefore appears to change as embryogenesis proceeds, from largely permissive topologies during cell-fate specification to more instructive regulation during terminal tissue differentiation, when E-P proximity is coupled to activation.

Chromatin gene-gene loops support the cross-regulation of genes with related function

Chromatin loops between gene pairs have been observed in diverse contexts in both flies and vertebrates. Combining high-resolution Capture-C, DNA fluorescence in situ hybridization, and genetic perturbations, we dissect the functional role of three loops between genes with related function during Drosophila embryogenesis. By mutating the loop anchor (but not the gene) or the gene (but not loop anchor), we disentangle loop formation and gene expression and show that the 3D proximity of paralogous gene loci supports their co-regulation. Breaking the loop leads to either an attenuation or enhancement of expression and perturbs their relative levels of expression and cross-regulation. Although many loops appear constitutive across embryogenesis, their function can change in different developmental contexts. Taken together, our results indicate that chromatin gene-gene loops act as architectural scaffolds that can be used in different ways in different contexts to fine-tune the coordinated expression of genes with related functions and sustain their cross-regulation.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

CTCF, BEAF-32, and CP190 are not required for the establishment of TADs in early Drosophila embryos but have locus-specific roles

The boundaries of topologically associating domains (TADs) are delimited by insulators and/or active promoters; however, how they are initially established during embryogenesis remains unclear. Here, we examined this during the first hours of Drosophila embryogenesis. DNA-FISH confirms that intra-TAD pairwise proximity is established during zygotic genome activation (ZGA) but with extensive cell-to-cell heterogeneity. Most newly formed boundaries are occupied by combinations of CTCF, BEAF-32, and/or CP190. Depleting each insulator individually from chromatin revealed that TADs can still establish, although with lower insulation, with a subset of boundaries (~10%) being more dependent on specific insulators. Some weakened boundaries have aberrant gene expression due to unconstrained enhancer activity. However, the majority of misexpressed genes have no obvious direct relationship to changes in domain-boundary insulation. Deletion of an active promoter (thereby blocking transcription) at one boundary had a greater impact than deleting the insulator-bound region itself. This suggests that cross-talk between insulators and active promoters and/or transcription might reinforce domain boundary insulation during embryogenesis.

Chromatin Immunoprecipitation Experiments from Drosophila Ovaries

Chromatin is composed of DNA and its associated proteins, and has an essential role in all cellular processes, including those taking place during Drosophila oogenesis. In order to understand the molecular basis of chromatin-based processes, such as transcription, it is essential to be able to study how and when different proteins, such as transcription factors, histones and RNA polymerases, interact with chromatin. One of the most popular methods to study this is chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Here, we describe a ChIP-seq protocol that has been optimized for Drosophila ovaries, focusing on sample preparation through preliminary data processing.

The continuum of Drosophila embryonic development at single-cell resolution

Drosophila melanogaster is a powerful, long-standing model for metazoan development and gene regulation. We profiled chromatin accessibility in almost 1 million and gene expression in half a million nuclei from overlapping windows spanning the entirety of embryogenesis. Leveraging developmental asynchronicity within embryo collections, we applied deep neural networks to infer the age of each nucleus, resulting in continuous, multimodal views of molecular and cellular transitions in absolute time. We identify cell lineages; infer their developmental relationships; and link dynamic changes in enhancer usage, transcription factor (TF) expression, and the accessibility of TFs' cognate motifs. With these data, the dynamics of enhancer usage and gene expression can be explored within and across lineages at the scale of minutes, including for precise transitions like zygotic genome activation.

Single-nucleus RNA-sequencing in pre-cellularization Drosophila melanogaster embryos

Our current understanding of the regulation of gene expression in the early Drosophila melanogaster embryo comes from observations of a few genes at a time, as with in situ hybridizations, or observation of gene expression levels without regards to patterning, as with RNA-sequencing. Single-nucleus RNA-sequencing however, has the potential to provide new insights into the regulation of gene expression for many genes at once while simultaneously retaining information regarding the position of each nucleus prior to dissociation based on patterned gene expression. In order to establish the use of single-nucleus RNA sequencing in Drosophila embryos prior to cellularization, here we look at gene expression in control and insulator protein, dCTCF, maternal null embryos during zygotic genome activation at nuclear cycle 14. We find that early embryonic nuclei can be grouped into distinct clusters according to gene expression. From both virtual and published in situ hybridizations, we also find that these clusters correspond to spatial regions of the embryo. Lastly, we provide a resource of candidate differentially expressed genes that might show local changes in gene expression between control and maternal dCTCF null nuclei with no detectable differential expression in bulk. These results highlight the potential for single-nucleus RNA-sequencing to reveal new insights into the regulation of gene expression in the early Drosophila melanogaster embryo.

Simultaneous cellular and molecular phenotyping of embryonic mutants using single-cell regulatory trajectories

Developmental progression and cellular diversity are largely driven by transcription factors (TFs); yet, characterizing their loss-of-function phenotypes remains challenging and often disconnected from their underlying molecular mechanisms. Here, we combine single-cell regulatory genomics with loss-of-function mutants to jointly assess both cellular and molecular phenotypes. Performing sci-ATAC-seq at eight overlapping time points during Drosophila mesoderm development could reconstruct the developmental trajectories of all major muscle types and reveal the TFs and enhancers involved. To systematically assess mutant phenotypes, we developed a single-nucleus genotyping strategy to process embryo pools of mixed genotypes. Applying this to four TF mutants could identify and quantify their characterized phenotypes de novo and discover new ones, while simultaneously revealing their regulatory input and mode of action. Our approach is a general framework to dissect the functional input of TFs in a systematic, unbiased manner, identifying both cellular and molecular phenotypes at a scale and resolution that has not been feasible before.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

In vivo tissue-specific chromatin profiling in Drosophila melanogaster using GFP-tagged nuclei

The chromatin landscape defines cellular identity in multicellular organisms with unique patterns of DNA accessibility and histone marks decorating the genome of each cell type. Thus, profiling the chromatin state of different cell types in an intact organism under disease or physiological conditions can provide insight into how chromatin regulates cell homeostasis in vivo. To overcome the many challenges associated with characterizing chromatin state in specific cell types, we developed an improved approach to isolate Drosophila melanogaster nuclei tagged with a GFPKASH protein. The perinuclear space-localized KASH domain anchors GFP to the outer nuclear membrane, and expression of UAS-GFPKASH can be controlled by tissue-specific Gal4 drivers. Using this protocol, we profiled chromatin accessibility using an improved version of Assay for Transposable Accessible Chromatin followed by sequencing (ATAC-seq), called Omni-ATAC. In addition, we examined the distribution of histone marks using Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and Cleavage Under Targets and Tagmentation (CUT&Tag) in adult photoreceptor neurons. We show that the chromatin landscape of photoreceptors reflects the transcriptional state of these cells, demonstrating the quality and reproducibility of our approach for profiling the transcriptome and epigenome of specific cell types in Drosophila.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

Lineage-Resolved Enhancer and Promoter Usage during a Time Course of Embryogenesis

Enhancers are essential drivers of cell states, yet the relationship between accessibility, regulatory activity, and in vivo lineage commitment during embryogenesis remains poorly understood. Here, we measure chromatin accessibility in isolated neural and mesodermal lineages across a time course of Drosophila embryogenesis. Promoters, including tissue-specific genes, are often constitutively open, even in contexts where the gene is not expressed. In contrast, the majority of distal elements have dynamic, tissue-specific accessibility. Enhancer priming appears rarely within a lineage, perhaps reflecting the speed of Drosophila embryogenesis. However, many tissue-specific enhancers are accessible in other lineages early on and become progressively closed as embryogenesis proceeds. We demonstrate the usefulness of this tissue- and time-resolved resource to definitively identify single-cell clusters, to uncover predictive motifs, and to identify many regulators of tissue development. For one such predicted neural regulator, l(3)neo38, we generate a loss-of-function mutant and uncover an essential role for neuromuscular junction and brain development.

Epigenetic Assays in Purified Cardiomyocyte Nuclei

The adult mammalian heart's potential for regeneration is very inefficient. Importantly, adult mammalian cardiomyocytes (CMs) are characterized as a cell population with very limited mitotic potential. Conversely, the neonatal mouse heart possesses a brief, yet robust, regenerative capacity within the first week of life. Cell type-specific enrichment procedures are essential for characterizing the full spectrum of epigenomic landscapes and gene regulatory networks deployed by mammalian CMs. In this chapter, we describe a protocol useful for purifying CM nuclei from mammalian cardiac tissue. Furthermore, we detail a low-input procedure suitable for the parallel genome-wide profiling of chromatin accessibility, histone modifications, and transcription factor-binding sites. The CM nuclei purified using this process are suitable for multi-omic profiling approaches.

Chromatin Immunoprecipitation

Chromatin immunoprecipitation, commonly referred to as ChIP, is a powerful technique for the evaluation of in vivo interactions of proteins with specific regions of genomic DNA. Formaldehyde is used in this technique to cross-link proteins to DNA in vivo, followed by the extraction of chromatin from cross-linked cells and tissues. Harvested chromatin is sheared and subsequently used in an immunoprecipitation incorporating antibodies specific to protein(s) of interest and thus coprecipitating and enriching the cross-linked, protein-associated DNA. The cross-linking process can be reversed, and protein-bound DNA fragments of optimal length ranging from 200 to 1000 base pairs (bp) can subsequently be purified and measured or sequenced by numerous analytical methods. In this protocol, two different fixation methods are described in detail. The first involves the standard fixation of cells and tissue by formaldehyde if the target antigen is highly abundant. The dual cross-linking procedure presented at the end includes an additional preformaldehyde cross-linking step and can be especially useful when the target protein is in low abundance or if it is indirectly associated with chromatin DNA through another protein.

A viral toolkit for recording transcription factor-DNA interactions in live mouse tissues

Transcription factors (TFs) enact precise regulation of gene expression through site-specific, genome-wide binding. Common methods for TF-occupancy profiling, such as chromatin immunoprecipitation, are limited by requirement of TF-specific antibodies and provide only end-point snapshots of TF binding. Alternatively, TF-tagging techniques, in which a TF is fused to a DNA-modifying enzyme that marks TF-binding events across the genome as they occur, do not require TF-specific antibodies and offer the potential for unique applications, such as recording of TF occupancy over time and cell type specificity through conditional expression of the TF-enzyme fusion. Here, we create a viral toolkit for one such method, calling cards, and demonstrate that these reagents can be delivered to the live mouse brain and used to report TF occupancy. Further, we establish a Cre-dependent calling cards system and, in proof-of-principle experiments, show utility in defining cell type-specific TF profiles and recording and integrating TF-binding events across time. This versatile approach will enable unique studies of TF-mediated gene regulation in live animal models.

Cell-type-specific genomics reveals histone modification dynamics in mammalian meiosis

Meiosis is the specialized cell division during which parental genomes recombine to create genotypically unique gametes. Despite its importance, mammalian meiosis cannot be studied in vitro, greatly limiting mechanistic studies. In vivo, meiocytes progress asynchronously through meiosis and therefore the study of specific stages of meiosis is a challenge. Here, we describe a method for isolating pure sub-populations of nuclei that allows for detailed study of meiotic substages. Interrogating the H3K4me3 landscape revealed dynamic chromatin transitions between substages of meiotic prophase I, both at sites of genetic recombination and at gene promoters. We also leveraged this method to perform the first comprehensive, genome-wide survey of histone marks in meiotic prophase, revealing a heretofore unappreciated complexity of the epigenetic landscape at meiotic recombination hotspots. Ultimately, this study presents a straightforward, scalable framework for interrogating the complexities of mammalian meiosis.

Meiotic DSB formation, repair and recombination occur in a continuum of substages termed leptonema, zygonema, pachynema, and diplonema. Here, authors develop a method for isolating pure sub-populations of nuclei that allows for detailed study of meiotic substages.

Tagger-A Swiss army knife for multiomics to dissect cell type-specific mechanisms of gene expression in mice

A deep understanding of how regulation of the multiple levels of gene expression in mammalian tissues give rise to complex phenotypes has been impeded by cellular diversity. A handful of techniques were developed to tag-select nucleic acids of interest in specific cell types, thereby enabling their capture. We expanded this strategy by developing the Tagger knock-in mouse line bearing a quad-cistronic transgene combining enrichment tools for nuclei, nascent RNA, translating mRNA, and mature microRNA (miRNA). We demonstrate that Tagger can capture the desired nucleic acids, enabling multiple omics approaches to be applied to specific cell types in vivo using a single transgenic mouse line.

This Methods and Resources paper describes Tagger, a knock-in mouse line bearing a quad-cistronic transgene that enables the capture of translating mRNAs, mature miRNAs, pulse-labeled total RNA, and the nucleus, all from specific cells of complex tissues.

Isolation and Quantification Brain Region-Specific and Cell Subtype-Specific Histone (De)Acetylation in Cognitive Neuroepigenetics

The acetylation of histone tails, which relaxes compact chromatin structure and enhances the accessibility of DNA to regulatory proteins, has emerged as a key mechanism for regulating gene expression. These modifications in turn play critical roles in forming long-term memories. Chromatin immunoprecipitation (ChIP) experiments have enabled the identification of specific histone modifications and the genes most closely associated with active memory formation. Problematically, however, the majority of these studies analyze diverse populations of cell homogenates obtained from the gross dissection of large brain regions. The protocol outlined here uses methods to ascribe gene-specific histone modifications (via specific antibodies and RT-qPCR) to specific cell subtypes (via specific antibodies and cell sorting) in discrete memory-related brain regions (via microdissection) to more precisely identify the role of histone acetylation and deacetylation in cognitive neuroepigenetics.

The Hox transcription factor Ubx stabilizes lineage commitment by suppressing cellular plasticity in Drosophila

During development cells become restricted in their differentiation potential by repressing alternative cell fates, and the Polycomb complex plays a crucial role in this process. However, how alternative fate genes are lineage-specifically silenced is unclear. We studied Ultrabithorax (Ubx), a multi-lineage transcription factor of the Hox class, in two tissue lineages using sorted nuclei and interfered with Ubx in mesodermal cells. We find that depletion of Ubx leads to the de-repression of genes normally expressed in other lineages. Ubx silences expression of alternative fate genes by retaining the Polycomb Group protein Pleiohomeotic at Ubx targeted genomic regions, thereby stabilizing repressive chromatin marks in a lineage-dependent manner. Our study demonstrates that Ubx stabilizes lineage choice by suppressing the multipotency encoded in the genome via its interaction with Pho. This mechanism may explain why the Hox code is maintained throughout the lifecycle, since it could set a block to transdifferentiation in adult cells.

Islands of retroelements are major components of Drosophila centromeres

Centromeres are essential chromosomal regions that mediate kinetochore assembly and spindle attachments during cell division. Despite their functional conservation, centromeres are among the most rapidly evolving genomic regions and can shape karyotype evolution and speciation across taxa. Although significant progress has been made in identifying centromere-associated proteins, the highly repetitive centromeres of metazoans have been refractory to DNA sequencing and assembly, leaving large gaps in our understanding of their functional organization and evolution. Here, we identify the sequence composition and organization of the centromeres of Drosophila melanogaster by combining long-read sequencing, chromatin immunoprecipitation for the centromeric histone CENP-A, and high-resolution chromatin fiber imaging. Contrary to previous models that heralded satellite repeats as the major functional components, we demonstrate that functional centromeres form on islands of complex DNA sequences enriched in retroelements that are flanked by large arrays of satellite repeats. Each centromere displays distinct size and arrangement of its DNA elements but is similar in composition overall. We discover that a specific retroelement, G2/Jockey-3, is the most highly enriched sequence in CENP-A chromatin and is the only element shared among all centromeres. G2/Jockey-3 is also associated with CENP-A in the sister species D. simulans, revealing an unexpected conservation despite the reported turnover of centromeric satellite DNA. Our work reveals the DNA sequence identity of the active centromeres of a premier model organism and implicates retroelements as conserved features of centromeric DNA.

De novo recruitment of Polycomb-group proteins in Drosophila embryos

Polycomb-group (PcG)-mediated transcriptional repression of target genes can be delineated into two phases. First, following initial repression of target genes by gene-specific transcription factors, PcG proteins recognize the repressed state and assume control of the genes' repression. Second, once the silenced state is established, PcG proteins may maintain repression through an indefinite number of cell cycles. Little is understood about how PcG proteins initially recognize the repressed state of target genes and the steps leading to de novo establishment of PcG-mediated repression. We describe a genetic system in which a Drosophila PcG target gene, giant (gt), is ubiquitously repressed during early embryogenesis by a maternally expressed transcription factor, and show the temporal recruitment of components of three PcG protein complexes: PhoRC, PRC1 and PRC2. We show that de novo PcG recruitment follows a temporal hierarchy in which PhoRC stably localizes at the target gene at least 1 h before stable recruitment of PRC2 and concurrent trimethylation of histone H3 at lysine 27 (H3K27me3). The presence of PRC2 and increased levels of H3K27me3 are found to precede stable binding by PRC1.

Spatially varying cis-regulatory divergence in Drosophila embryos elucidates cis-regulatory logic

Spatial patterning of gene expression is a key process in development, yet how it evolves is still poorly understood. Both cis- and trans-acting changes could participate in complex interactions, so to isolate the cis-regulatory component of patterning evolution, we measured allele-specific spatial gene expression patterns in D. melanogaster × simulans hybrid embryos. RNA-seq of cryo-sectioned slices revealed 66 genes with strong spatially varying allele-specific expression. We found that hunchback, a major regulator of developmental patterning, had reduced expression of the D. simulans allele specifically in the anterior tip of hybrid embryos. Mathematical modeling of hunchback cis-regulation suggested a candidate transcription factor binding site variant, which we verified as causal using CRISPR-Cas9 genome editing. In sum, even comparing morphologically near-identical species we identified surprisingly extensive spatial variation in gene expression, suggesting not only that development is robust to many such changes, but also that natural selection may have ample raw material for evolving new body plans via changes in spatial patterning.

Eukaryotic core promoters and the functional basis of transcription initiation

RNA polymerase II (Pol II) core promoters are specialized DNA sequences at transcription start sites of protein-coding and non-coding genes that support the assembly of the transcription machinery and transcription initiation. They enable the highly regulated transcription of genes by selectively integrating regulatory cues from distal enhancers and their associated regulatory proteins. In this Review, we discuss the defining properties of gene core promoters, including their sequence features, chromatin architecture and transcription initiation patterns. We provide an overview of molecular mechanisms underlying the function and regulation of core promoters and their emerging functional diversity, which defines distinct transcription programmes. On the basis of the established properties of gene core promoters, we discuss transcription start sites within enhancers and integrate recent results obtained from dedicated functional assays to propose a functional model of transcription initiation. This model can explain the nature and function of transcription initiation at gene starts and at enhancers and can explain the different roles of core promoters, of Pol II and its associated factors and of the activating cues provided by enhancers and the transcription factors and cofactors they recruit.

DamID profiling of dynamic Polycomb-binding sites in Drosophila imaginal disc development and tumorigenesis

Background

Tracking dynamic protein–chromatin interactions in vivo is key to unravel transcriptional and epigenetic transitions in development and disease. However, limited availability and heterogeneous tissue composition of in vivo source material impose challenges on many experimental approaches.

Results

Here we adapt cell-type-specific DamID-seq profiling for use in Drosophila imaginal discs and make FLP/FRT-based induction accessible to GAL driver-mediated targeting of specific cell lineages. In a proof-of-principle approach, we utilize ubiquitous DamID expression to describe dynamic transitions of Polycomb-binding sites during wing imaginal disc development and in a scrib tumorigenesis model. We identify Atf3 and Ets21C as novel Polycomb target genes involved in scrib tumorigenesis and suggest that target gene regulation by Atf3 and AP-1 transcription factors, as well as modulation of insulator function, plays crucial roles in dynamic Polycomb-binding at target sites. We establish these findings by DamID-seq analysis of wing imaginal disc samples derived from 10 larvae.

Conclusions

Our study opens avenues for robust profiling of small cell population in imaginal discs in vivo and provides insights into epigenetic changes underlying transcriptional responses to tumorigenic transformation.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0196-y) contains supplementary material, which is available to authorized users.

The cis-regulatory dynamics of embryonic development at single-cell resolution

Understanding how gene regulatory networks control the progressive restriction of cell fates is a long-standing challenge. Recent advances in measuring single cell gene expression are providing new insights into lineage commitment. However, the regulatory events underlying these changes remain elusive. Here we investigate the dynamics of chromatin regulatory landscapes during embryogenesis at single cell resolution. Using single cell combinatorial indexing assay for transposase accessible chromatin (sci-ATAC-seq)1, we profiled chromatin accessibility in over 20,000 single nuclei from fixed Drosophila embryos spanning three landmark embryonic stages: 2-4 hours (hrs) after egg laying (predominantly stage 5 blastoderm nuclei), when each embryo comprises ~6,000 multipotent cells; 6-8hrs (predominantly stage 10-11), to capture a midpoint in embryonic development when major lineages in the mesoderm and ectoderm are specified; and 10-12hrs (predominantly stage 13), when each of the embryo’s >20,000 cells are undergoing terminal differentiation. Our results reveal spatial heterogeneity in the usage of the regulatory genome prior to gastrulation, a feature that aligns with future cell fate, and nuclei can be temporally ordered along developmental trajectories. During mid-embryogenesis, tissue granularity emerges such that individual cell types can be inferred by their chromatin accessibility, while maintaining a signature of their germ layer of origin. The data reveal overlapping usage of regulatory elements between cells of the endoderm and non-myogenic mesoderm, suggesting a common developmental program reminiscent of the mesendoderm lineage in other species2–4. Altogether, we identify over 30,000 distal regulatory elements exhibiting tissue-specific accessibility. We validated the germ layer specificity of a subset of these predicted enhancers in transgenic embryos, achieving 90% accuracy. Overall, our results demonstrate the power of shotgun single cell profiling of embryos to resolve dynamic changes in the chromatin landscape during development, and to uncover the cis-regulatory programs of metazoan germ layers and cell types.

The degree of enhancer or promoter activity is reflected by the levels and directionality of eRNA transcription

Here, Mikhaylichenko et al. investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. The authors demonstrate that while the timing of enhancer transcription is correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers and conclude that this is likely an inherent sequence property of the elements themselves.

Gene expression is regulated by promoters, which initiate transcription, and enhancers, which control their temporal and spatial activity. However, the discovery that mammalian enhancers also initiate transcription questions the inherent differences between enhancers and promoters. Here, we investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. We show that while the timing of enhancer transcription is generally correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers. To assess how this impacts function, we developed a dual transgenic assay to simultaneously measure enhancer and promoter activities from a single element in the same embryo. Extensive transgenic analysis revealed a relationship between the direction of endogenous transcription and the ability to function as an enhancer or promoter in vivo, although enhancer RNA (eRNA) production and activity are not always strictly coupled. Some enhancers (mainly bidirectional) can act as weak promoters, producing overlapping spatio–temporal expression. Conversely, bidirectional promoters often act as strong enhancers, while unidirectional promoters generally cannot. The balance between enhancer and promoter activity is generally reflected in the levels and directionality of eRNA transcription and is likely an inherent sequence property of the elements themselves.

MEF2C transcription factor is associated with the genetic and epigenetic risk architecture of schizophrenia and improves cognition in mice

Large-scale consortia mapping the genomic risk architectures of schizophrenia provide vast amounts of molecular information, with largely unexplored therapeutic potential. We harnessed publically available information from the Psychiatric Genomics Consortium, and report myocyte enhancer factor 2C (MEF2C) motif enrichment in sequences surrounding the top scoring single-nucleotide polymorphisms within risk loci contributing by individual small effect to disease heritability. Chromatin profiling at base-pair resolution in neuronal nucleosomes extracted from prefrontal cortex of 34 subjects, including 17 cases diagnosed with schizophrenia, revealed MEF2C motif enrichment within cis-regulatory sequences, including neuron-specific promoters and superenhancers, affected by histone H3K4 hypermethylation in disease cases. Vector-induced short- and long-term Mef2c upregulation in mouse prefrontal projection neurons consistently resulted in enhanced cognitive performance in working memory and object recognition paradigms at baseline and after psychotogenic drug challenge, in conjunction with remodeling of local connectivity. Neuronal genome tagging in vivo by Mef2c-Dam adenine methyltransferase fusion protein confirmed the link between cognitive enhancement and MEF2C occupancy at promoters harboring canonical and variant MEF2C motifs. The multilayered integrative approaches presented here provide a roadmap to uncover the therapeutic potential of transcriptional regulators for schizophrenia and related disorders.

Simultaneous Transcriptional and Epigenomic Profiling from Specific Cell Types within Heterogeneous Tissues In Vivo

Epigenomic mechanisms direct distinct gene expression programs for different cell types. Various in vivo tissues have been subjected to epigenomic analysis; however, these studies have been limited by cellular heterogeneity, resulting in composite gene expression and epigenomic profiles. Here, we introduce "NuTRAP," a transgenic mouse that allows simultaneous isolation of cell-type-specific translating mRNA and chromatin from complex tissues. Using NuTRAP, we successfully characterize gene expression and epigenomic states of various adipocyte populations in vivo, revealing significant differences compared to either whole adipose tissue or in vitro adipocyte cell lines. We find that chromatin immunoprecipitation sequencing (ChIP-seq) using NuTRAP is highly efficient, scalable, and robust with even limited cell input. We further demonstrate the general utility of NuTRAP by analyzing hepatocyte-specific epigenomic states. The NuTRAP mouse is a resource that provides a powerful system for cell-type-specific gene expression and epigenomic profiling.

McEnhancer: predicting gene expression via semi-supervised assignment of enhancers to target genes

Transcriptional enhancers regulate spatio-temporal gene expression. While genomic assays can identify putative enhancers en masse, assigning target genes is a complex challenge. We devised a machine learning approach, McEnhancer, which links target genes to putative enhancers via a semi-supervised learning algorithm that predicts gene expression patterns based on enriched sequence features. Predicted expression patterns were 73–98% accurate, predicted assignments showed strong Hi-C interaction enrichment, enhancer-associated histone modifications were evident, and known functional motifs were recovered. Our model provides a general framework to link globally identified enhancers to targets and contributes to deciphering the regulatory genome.

The Drosophila embryo at single-cell transcriptome resolution

By the onset of morphogenesis, Drosophila embryos consist of about 6000 cells that express distinct gene combinations. Here, we used single-cell sequencing of precisely staged embryos and devised DistMap, a computational mapping strategy to reconstruct the embryo and to predict spatial gene expression approaching single-cell resolution. We produced a virtual embryo with about 8000 expressed genes per cell. Our interactive Drosophila Virtual Expression eXplorer (DVEX) database generates three-dimensional virtual in situ hybridizations and computes gene expression gradients. We used DVEX to uncover patterned expression of transcription factors and long noncoding RNAs, as well as signaling pathway components. Spatial regulation of Hippo signaling during early embryogenesis suggests a mechanism for establishing asynchronous cell proliferation. Our approach is suitable to generate transcriptomic blueprints for other complex tissues.

Dual functionality of cis-regulatory elements as developmental enhancers and Polycomb response elements

Here, Erceg et al. studied the occupancy of the Drosophila PhoRC during embryogenesis and revealed extensive co-occupancy at developmental enhancers. By using an established in vivo assay for Polycomb response element (PRE) activity, they show that a subset of characterized developmental enhancers can function as PREs and silence transcription in a Polycomb-dependent manner, thereby suggesting that reuse of the same elements by the PcG system may help fine-tune gene expression and ensure the timely maintenance of cell identities.

Developmental gene expression is tightly regulated through enhancer elements, which initiate dynamic spatio–temporal expression, and Polycomb response elements (PREs), which maintain stable gene silencing. These two cis-regulatory functions are thought to operate through distinct dedicated elements. By examining the occupancy of the Drosophila pleiohomeotic repressive complex (PhoRC) during embryogenesis, we revealed extensive co-occupancy at developmental enhancers. Using an established in vivo assay for PRE activity, we demonstrated that a subset of characterized developmental enhancers can function as PREs, silencing transcription in a Polycomb-dependent manner. Conversely, some classic Drosophila PREs can function as developmental enhancers in vivo, activating spatio–temporal expression. This study therefore uncovers elements with dual function: activating transcription in some cells (enhancers) while stably maintaining transcriptional silencing in others (PREs). Given that enhancers initiate spatio–temporal gene expression, reuse of the same elements by the Polycomb group (PcG) system may help fine-tune gene expression and ensure the timely maintenance of cell identities.

Mechanisms of Groucho-mediated repression revealed by genome-wide analysis of Groucho binding and activity

Background

The transcriptional corepressor Groucho (Gro) is required for the function of many developmentally regulated DNA binding repressors, thus helping to define the gene expression profile of each cell during development. The ability of Gro to repress transcription at a distance together with its ability to oligomerize and bind to histones has led to the suggestion that Gro may spread along chromatin. However, much is unknown about the mechanism of Gro-mediated repression and about the dynamics of Gro targeting.

Results

Our chromatin immunoprecipitation sequencing analysis of temporally staged Drosophila embryos shows that Gro binds in a highly dynamic manner primarily to clusters of discrete (<1 kb) segments. Consistent with the idea that Gro may facilitate communication between silencers and promoters, Gro binding is enriched at both cis-regulatory modules, as well as within the promotors of potential target genes. While this Gro-recruitment is required for repression, our data show that it is not sufficient for repression. Integration of Gro binding data with transcriptomic analysis suggests that, contrary to what has been observed for another Gro family member, Drosophila Gro is probably a dedicated repressor. This analysis also allows us to define a set of high confidence Gro repression targets. Using publically available data regarding the physical and genetic interactions between these targets, we are able to place them in the regulatory network controlling development. Through analysis of chromatin associated pre-mRNA levels at these targets, we find that genes regulated by Gro in the embryo are enriched for characteristics of promoter proximal paused RNA polymerase II.

Conclusions

Our findings are inconsistent with a one-dimensional spreading model for long-range repression and suggest that Gro-mediated repression must be regulated at a post-recruitment step. They also show that Gro is likely a dedicated repressor that sits at a prominent highly interconnected regulatory hub in the developmental network. Furthermore, our findings suggest a role for RNA polymerase II pausing in Gro-mediated repression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3589-6) contains supplementary material, which is available to authorized users.

Identification and in silico modeling of enhancers reveals new features of the cardiac differentiation network

Developmental patterning and tissue formation are regulated through complex gene regulatory networks (GRNs) driven through the action of transcription factors (TFs) converging on enhancer elements. Here, as a point of entry to dissect the poorly defined GRN underlying cardiomyocyte differentiation, we apply an integrated approach to identify active enhancers and TFs involved in Drosophila heart development. The Drosophila heart consists of 104 cardiomyocytes, representing less than 0.5% of all cells in the embryo. By modifying BiTS-ChIP for rare cells, we examined H3K4me3 and H3K27ac chromatin landscapes to identify active promoters and enhancers specifically in cardiomyocytes. These in vivo data were complemented by a machine learning approach and extensive in vivo validation in transgenic embryos, which identified many new heart enhancers and their associated TF motifs. Our results implicate many new TFs in late stages of heart development, including Bagpipe, an Nkx3.2 ortholog, which we show is essential for differentiated heart function.

Genome-wide chromatin and gene expression profiling during memory formation and maintenance in adult mice

Recent evidence suggests that the formation and maintenance of memory requires epigenetic changes. In an effort to understand the spatio-temporal extent of learning and memory-related epigenetic changes we have charted genome-wide histone and DNA methylation profiles, in two different brain regions, two cell types, and three time-points, before and after learning. In this data descriptor we provide detailed information on data generation, give insights into the rationale of experiments, highlight necessary steps to assess data quality, offer guidelines for future use of the data and supply ready-to-use code to replicate the analysis results. The data provides a blueprint of the gene regulatory network underlying short- and long-term memory formation and maintenance. This 'healthy' gene regulatory network of learning can now be compared to changes in neurological or psychiatric diseases, providing mechanistic insights into brain disorders and highlighting potential therapeutic avenues.

Genome-Wide Ultrabithorax Binding Analysis Reveals Highly Targeted Genomic Loci at Developmental Regulators and a Potential Connection to Polycomb-Mediated Regulation

Hox homeodomain transcription factors are key regulators of animal development. They specify the identity of segments along the anterior-posterior body axis in metazoans by controlling the expression of diverse downstream targets, including transcription factors and signaling pathway components. The Drosophila melanogaster Hox factor Ultrabithorax (Ubx) directs the development of thoracic and abdominal segments and appendages, and loss of Ubx function can lead for example to the transformation of third thoracic segment appendages (e.g. halters) into second thoracic segment appendages (e.g. wings), resulting in a characteristic four-wing phenotype. Here we present a Drosophila melanogaster strain with a V5-epitope tagged Ubx allele, which we employed to obtain a high quality genome-wide map of Ubx binding sites using ChIP-seq. We confirm the sensitivity of the V5 ChIP-seq by recovering 7/8 of well-studied Ubx-dependent cis-regulatory regions. Moreover, we show that Ubx binding is predictive of enhancer activity as suggested by comparison with a genome-scale resource of in vivo tested enhancer candidates. We observed densely clustered Ubx binding sites at 12 extended genomic loci that included ANTP-C, BX-C, Polycomb complex genes, and other regulators and the clustered binding sites were frequently active enhancers. Furthermore, Ubx binding was detected at known Polycomb response elements (PREs) and was associated with significant enrichments of Pc and Pho ChIP signals in contrast to binding sites of other developmental TFs. Together, our results show that Ubx targets developmental regulators via strongly clustered binding sites and allow us to hypothesize that regulation by Ubx might involve Polycomb group proteins to maintain specific regulatory states in cooperative or mutually exclusive fashion, an attractive model that combines two groups of proteins with prominent gene regulatory roles during animal development.

Cell-type-specific profiling of protein-DNA interactions without cell isolation using targeted DamID with next-generation sequencing

This protocol is an extension to: Nat. Protoc. 2, 1467-1478 (2007); doi:10.1038/nprot.2007.148; published online 7 June 2007The ability to profile transcription and chromatin binding in a cell-type-specific manner is a powerful aid to understanding cell-fate specification and cellular function in multicellular organisms. We recently developed targeted DamID (TaDa) to enable genome-wide, cell-type-specific profiling of DNA- and chromatin-binding proteins in vivo without cell isolation. As a protocol extension, this article describes substantial modifications to an existing protocol, and it offers additional applications. TaDa builds upon DamID, a technique for detecting genome-wide DNA-binding profiles of proteins, by coupling it with the GAL4 system in Drosophila to enable both temporal and spatial resolution. TaDa ensures that Dam-fusion proteins are expressed at very low levels, thus avoiding toxicity and potential artifacts from overexpression. The modifications to the core DamID technique presented here also increase the speed of sample processing and throughput, and adapt the method to next-generation sequencing technology. TaDa is robust, reproducible and highly sensitive. Compared with other methods for cell-type-specific profiling, the technique requires no cell-sorting, cross-linking or antisera, and binding profiles can be generated from as few as 10,000 total induced cells. By profiling the genome-wide binding of RNA polymerase II (Pol II), TaDa can also identify transcribed genes in a cell-type-specific manner. Here we describe a detailed protocol for carrying out TaDa experiments and preparing the material for next-generation sequencing. Although we developed TaDa in Drosophila, it should be easily adapted to other organisms with an inducible expression system. Once transgenic animals are obtained, the entire experimental procedure-from collecting tissue samples to generating sequencing libraries-can be accomplished within 5 d.

Cell-type-specific profiling of protein-DNA interactions without cell isolation using targeted DamID with next-generation sequencing

This protocol is an extension to: Nat. Protoc. 2, 1467-1478 (2007); doi:10.1038/nprot.2007.148; published online 7 June 2007The ability to profile transcription and chromatin binding in a cell-type-specific manner is a powerful aid to understanding cell-fate specification and cellular function in multicellular organisms. We recently developed targeted DamID (TaDa) to enable genome-wide, cell-type-specific profiling of DNA- and chromatin-binding proteins in vivo without cell isolation. As a protocol extension, this article describes substantial modifications to an existing protocol, and it offers additional applications. TaDa builds upon DamID, a technique for detecting genome-wide DNA-binding profiles of proteins, by coupling it with the GAL4 system in Drosophila to enable both temporal and spatial resolution. TaDa ensures that Dam-fusion proteins are expressed at very low levels, thus avoiding toxicity and potential artifacts from overexpression. The modifications to the core DamID technique presented here also increase the speed of sample processing and throughput, and adapt the method to next-generation sequencing technology. TaDa is robust, reproducible and highly sensitive. Compared with other methods for cell-type-specific profiling, the technique requires no cell-sorting, cross-linking or antisera, and binding profiles can be generated from as few as 10,000 total induced cells. By profiling the genome-wide binding of RNA polymerase II (Pol II), TaDa can also identify transcribed genes in a cell-type-specific manner. Here we describe a detailed protocol for carrying out TaDa experiments and preparing the material for next-generation sequencing. Although we developed TaDa in Drosophila, it should be easily adapted to other organisms with an inducible expression system. Once transgenic animals are obtained, the entire experimental procedure-from collecting tissue samples to generating sequencing libraries-can be accomplished within 5 d.

Strategies for targeting primate neural circuits with viral vectors

Understanding how the brain works requires understanding how different types of neurons contribute to circuit function and organism behavior. Progress on this front has been accelerated by optogenetics and chemogenetics, which provide an unprecedented level of control over distinct neuronal types in small animals. In primates, however, targeting specific types of neurons with these tools remains challenging. In this review, we discuss existing and emerging strategies for directing genetic manipulations to targeted neurons in the adult primate central nervous system. We review the literature on viral vectors for gene delivery to neurons, focusing on adeno-associated viral vectors and lentiviral vectors, their tropism for different cell types, and prospects for new variants with improved efficacy and selectivity. We discuss two projection targeting approaches for probing neural circuits: anterograde projection targeting and retrograde transport of viral vectors. We conclude with an analysis of cell type-specific promoters and other nucleotide sequences that can be used in viral vectors to target neuronal types at the transcriptional level.

Chromatin Immunoprecipitation for Analyzing Transcription Factor Binding and Histone Modifications in Drosophila

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is an invaluable technique to assess transcription factor binding and histone modifications in a genome-wide manner, an essential step towards understanding the mechanisms that govern embryonic development. Here, we provide a detailed protocol for all steps involved in generating a ChIP-seq library, starting from embryo collection, fixation, chromatin preparation, immunoprecipitation, and finally library preparation. The protocol is optimized for Drosophila embryos, but can be easily adapted for any model organism. The resulting library is suitable for sequencing on an Illumina HiSeq or MiSeq platform.

Determination of High-Resolution 3D Chromatin Organization Using Circular Chromosome Conformation Capture (4C-seq)

3D chromatin organization is essential for many aspects of transcriptional regulation. Circular Chromosome Conformation Capture followed by Illumina sequencing (4C-seq) is among the most powerful techniques to determine 3D chromatin organization. 4C-seq, like other modifications of the original 3C technique, uses the principle of "proximity ligation" to identify and quantify ten thousands of genomic interactions at a kilobase scale in a single experiment for predefined loci in the genome.In this chapter we focus on the experimental steps in the 4C-seq protocol, providing detailed descriptions on the preparation of cells, the construction of the circularized 3C library and the generation of the Illumina high throughput sequencing library. This protocol is particularly suited for the use of mammalian tissue samples, but can be used with minimal changes on circulating cells and cell lines from other sources as well. In the final section of this chapter, we provide a brief overview of data analysis approaches, accompanied by links to publicly available analysis tools.

A Genome-Wide Perspective on Metabolism

Mammals have at least 210 histologically diverse cell types (Alberts, Molecular biology of the cell. Garland Science, New York, 2008) and the number would be even higher if functional differences are taken into account. The genome in each of these cell types is differentially programmed to express the specific set of genes needed to fulfill the phenotypical requirements of the cell. Furthermore, in each of these cell types, the gene program can be differentially modulated by exposure to external signals such as hormones or nutrients. The basis for the distinct gene programs relies on cell type-selective activation of transcriptional enhancers, which in turn are particularly sensitive to modulation. Until recently we had only fragmented insight into the regulation of a few of these enhancers; however, the recent advances in high-throughput sequencing technologies have enabled the development of a large number of technologies that can be used to obtain genome-wide insight into how genomes are reprogrammed during development and in response to specific external signals. By applying such technologies, we have begun to reveal the cross-talk between metabolism and the genome, i.e., how genomes are reprogrammed in response to metabolites, and how the regulation of metabolic networks is coordinated at the genomic level.

Methylation of RNA polymerase II non-consensus Lysine residues marks early transcription in mammalian cells

Dynamic post-translational modification of RNA polymerase II (RNAPII) coordinates the co-transcriptional recruitment of enzymatic complexes that regulate chromatin states and processing of nascent RNA. Extensive phosphorylation of serine residues at the largest RNAPII subunit occurs at its structurally-disordered C-terminal domain (CTD), which is composed of multiple heptapeptide repeats with consensus sequence Y1-S2-P3-T4-S5-P6-S7. Serine-5 and Serine-7 phosphorylation mark transcription initiation, whereas Serine-2 phosphorylation coincides with productive elongation. In vertebrates, the CTD has eight non-canonical substitutions of Serine-7 into Lysine-7, which can be acetylated (K7ac). Here, we describe mono- and di-methylation of CTD Lysine-7 residues (K7me1 and K7me2). K7me1 and K7me2 are observed during the earliest transcription stages and precede or accompany Serine-5 and Serine-7 phosphorylation. In contrast, K7ac is associated with RNAPII elongation, Serine-2 phosphorylation and mRNA expression. We identify an unexpected balance between RNAPII K7 methylation and acetylation at gene promoters, which fine-tunes gene expression levels.

DNA methylation changes in plasticity genes accompany the formation and maintenance of memory

The ability to form memories is a prerequisite for an organism's behavioral adaptation to environmental changes. At the molecular level, the acquisition and maintenance of memory requires changes in chromatin modifications. In an effort to unravel the epigenetic network underlying both short- and long-term memory, we examined chromatin modification changes in two distinct mouse brain regions, two cell types and three time points before and after contextual learning. We found that histone modifications predominantly changed during memory acquisition and correlated surprisingly little with changes in gene expression. Although long-lasting changes were almost exclusive to neurons, learning-related histone modification and DNA methylation changes also occurred in non-neuronal cell types, suggesting a functional role for non-neuronal cells in epigenetic learning. Finally, our data provide evidence for a molecular framework of memory acquisition and maintenance, wherein DNA methylation could alter the expression and splicing of genes involved in functional plasticity and synaptic wiring.

MOBE-ChIP: a large-scale chromatin immunoprecipitation assay for cell type-specific studies

Cell type-specific transcriptional regulators play critical roles in the generation and maintenance of multicellularity. As they are often expressed at low levels, in vivo DNA-binding studies of these regulators by standard chromatin immunoprecipitation (ChIP) assays are technically challenging. We describe here an optimized ChIP protocol termed Maximized Objects for Better Enrichment (MOBE)-ChIP, which enhances the sensitivity of ChIP assays for detecting cell type-specific signals. The protocol, which is based on the disproportional increase of target signals over background at higher scales, uses substantially greater volume of starting materials than conventional ChIPs to achieve high signal enrichment. This technique can capture weak binding events that are ambiguous in standard ChIP assays, and is useful both in gene-specific and whole-genome analysis. This protocol has been optimized for Arabidopsis, but should be applicable to other model systems with minor modifications. The full procedure can be completed within 3 days.

Getting Down to Specifics: Profiling Gene Expression and Protein-DNA Interactions in a Cell Type-Specific Manner

The majority of multicellular organisms are comprised of an extraordinary range of cell types, with different properties and gene expression profiles. Understanding what makes each cell type unique and how their individual characteristics are attributed are key questions for both developmental and neurobiologists alike. The brain is an excellent example of the cellular diversity expressed in the majority of eukaryotes. The mouse brain comprises of approximately 75 million neurons varying in morphology, electrophysiology, and preferences for synaptic partners. A powerful process in beginning to pick apart the mechanisms that specify individual characteristics of the cell, as well as their fate, is to profile gene expression patterns, chromatin states, and transcriptional networks in a cell type-specific manner, i.e., only profiling the cells of interest in a particular tissue. Depending on the organism, the questions being investigated, and the material available, certain cell type-specific profiling methods are more suitable than others. This chapter reviews the approaches presently available for selecting and isolating specific cell types and evaluates their key features.

GATA-dependent transcriptional and epigenetic control of cardiac lineage specification and differentiation

Heart progenitor cells differentiate into various cell types including pacemaker and working cardiomyocytes. Cell-type specific gene expression is achieved by combinatorial interactions between tissue-specific transcription factors (TFs), co-factors, and chromatin remodelers and DNA binding elements in regulatory regions. Dysfunction of these transcriptional networks may result in congenital heart defects. Functional analysis of the regulatory DNA sequences has contributed substantially to the identification of the transcriptional network components and combinatorial interactions regulating the tissue-specific gene programs. GATA TFs have been identified as central players in these networks. In particular, GATA binding elements have emerged as a platform to recruit broadly active histone modification enzymes and cell-type-specific co-factors to drive cell-type-specific gene programs. Here, we discuss the role of GATA factors in cell fate decisions and differentiation in the developing heart.