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

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.

Epigenetic pioneering by SWI/SNF family remodelers

In eukaryotic genomes, transcriptional machinery and nucleosomes compete for binding to DNA sequences; thus, a crucial aspect of gene regulatory element function is to modulate chromatin accessibility for transcription factor (TF) and RNA polymerase binding. Recent structural studies have revealed multiple modes of TF engagement with nucleosomes, but how initial "pioneering" results in steady-state DNA accessibility for further TF binding and RNA polymerase II (RNAPII) engagement has been unclear. Even less well understood is how distant sites of open chromatin interact with one another, such as when developmental enhancers activate promoters to release RNAPII for productive elongation. Here, we review evidence for the centrality of the conserved SWI/SNF family of nucleosome remodeling complexes, both in pioneering and in mediating enhancer-promoter contacts. Consideration of the nucleosome unwrapping and ATP hydrolysis activities of SWI/SNF complexes, together with their architectural features, may reconcile steady-state TF occupancy with rapid TF dynamics observed by live imaging.

Tissue-specific RNA Polymerase II promoter-proximal pause release and burst kinetics in a Drosophila embryonic patterning network

BACKGROUND

Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood.

RESULTS

Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues.

CONCLUSIONS

The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

cisDynet: An integrated platform for modeling gene-regulatory dynamics and networks

Chromatin accessibility sequencing has been widely used for uncovering genetic regulatory mechanisms and inferring gene regulatory networks. However, effectively integrating large-scale chromatin accessibility datasets has posed a significant challenge. This is due to the lack of a comprehensive end-to-end solution, as many existing tools primarily emphasize data preprocessing and overlook downstream analyses. To bridge this gap, we have introduced cisDynet, a holistic solution that combines streamlined data preprocessing using Snakemake and R functions with advanced downstream analysis capabilities. cisDynet excels in conventional data analyses, encompassing peak statistics, peak annotation, differential analysis, motif enrichment analysis, and more. Additionally, it allows to perform sophisticated data exploration, such as tissue-specific peak identification, time course data modeling, integration of RNA-seq data to establish peak-to-gene associations, constructing regulatory networks, and conducting enrichment analysis of genome-wide association study (GWAS) variants. As a proof of concept, we applied cisDynet to reanalyze comprehensive ATAC-seq datasets across various tissues from the Encyclopedia of DNA Elements (ENCODE) project. The analysis successfully delineated tissue-specific open chromatin regions (OCRs), established connections between OCRs and target genes, and effectively linked these discoveries with 1861 GWAS variants. Furthermore, cisDynet was instrumental in dissecting the time course open chromatin data of mouse embryonic development, revealing the dynamic behavior of OCRs over developmental stages and identifying key transcription factors governing differentiation trajectories. In summary, cisDynet offers researchers a user-friendly solution that minimizes the need for extensive coding, ensures the reproducibility of results, and greatly simplifies the exploration of epigenomic data.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

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.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Hybrid epigenomes reveal extensive local genetic changes to chromatin accessibility contribute to divergence in embryonic gene expression between species

Chromatin accessibility plays an important role in shaping gene expression patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system, and suggest that the regulation of chromatin is more genetically complex than was previously appreciated.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

Factors and Mechanisms That Influence Chromatin-Mediated Enhancer–Promoter Interactions and Transcriptional Regulation

Simple Summary

The physical interactions between enhancers and promoters create chromatin conformations involved in gene regulation. In cancer cells, the chromatin conformations can be altered with uncontrolled deposition of histone marks resulting in varied gene expression. Although it is not entirely comprehensive how chromatin-mediated enhancer–promoter (E–P) interactions with various histone marks can affect gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential to control gene expression. In this review, we focus on emerging views of chromatin conformations associated with the E–P interactions and factors that establish or maintain such interactions, which may regulate gene expression.

Abstract

Eukaryotic gene expression is regulated through chromatin conformation, in which enhancers and promoters physically interact (E–P interactions). How such chromatin-mediated E–P interactions affect gene expression is not yet fully understood, but the roles of histone acetylation and methylation, pioneer transcription factors, and architectural proteins such as CCCTC binding factor (CTCF) and cohesin have recently attracted attention. Moreover, accumulated data suggest that E–P interactions are mechanistically involved in biophysical events, including liquid–liquid phase separation, and in biological events, including cancers. In this review, we discuss various mechanisms that regulate eukaryotic gene expression, focusing on emerging views regarding chromatin conformations that are involved in E–P interactions and factors that establish and maintain them.

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.

Promoter architecture of Drosophila genes regulated by Myocyte enhancer factor-2

To gain understanding into the mechanisms of transcriptional activation of muscle genes, we sought to determine if genes targeted by the myogenic transcription factor Myocyte enhancer factor-2 (MEF2) were enriched for specific core promoter elements. We identified 330 known MEF2 target promoters in Drosophila, and analyzed them for for the presence and location of 17 known consensus promoter sequences. As a control, we also searched all Drosophila RNA polymerase II-dependent promoters for the same sequences. We found that promoter motifs were readily detected in the MEF2 target dataset, and that many of them were slightly enriched in frequency compared to the control dataset. A prominent sequence over-represented in the MEF2 target genes was NDM2, that appeared in over 50% of MEF2 target genes and was 2.5-fold over-represented in MEF2 targets compared to background. To test the functional significance of NDM2, we identified two promoters containing a single copy of NDM2 plus an upstream MEF2 site, and tested the activity of these promoters in vivo. Both the sticks and stones and Kahuli fragments showed strong skeletal myoblast-specific expression of a lacZ reporter in embryos. However, the timing and level of reporter expression was unaffected when the NDM2 site in either element was mutated. These studies identify variations in promoter architecture for a set of regulated genes compared to all RNA polymerase II-dependent genes, and underline the potential redundancy in the activities of some core promoter elements.

Essential role of Cp190 in physical and regulatory boundary formation

Boundaries in animal genomes delimit contact domains with enhanced internal contact frequencies and have debated functions in limiting regulatory cross-talk between domains and guiding enhancers to target promoters. Most mammalian boundaries form by stalling of chromosomal loop-extruding cohesin by CTCF, but most Drosophila boundaries form CTCF independently. However, how CTCF-independent boundaries form and function remains largely unexplored. Here, we assess genome folding and developmental gene expression in fly embryos lacking the ubiquitous boundary-associated factor Cp190. We find that sequence-specific DNA binding proteins such as CTCF and Su(Hw) directly interact with and recruit Cp190 to form most promoter-distal boundaries. Cp190 is essential for early development and prevents regulatory cross-talk between specific gene loci that pattern the embryo. Cp190 was, in contrast, dispensable for long-range enhancer-promoter communication at tested loci. Cp190 is thus currently the major player in fly boundary formation and function, revealing that diverse mechanisms evolved to partition genomes into independent regulatory domains.

ASC proneural factors are necessary for chromatin remodeling during neuroectodermal to neuroblast fate transition to ensure the timely initiation of the neural stem cell program

Background

In both Drosophila and mammals, the achaete-scute (ASC/ASCL) proneural bHLH transcription factors are expressed in the developing central and peripheral nervous systems, where they function during specification and maintenance of the neural stem cells in opposition to Notch signaling. In addition to their role in nervous system development, ASC transcription factors are oncogenic and exhibit chromatin reprogramming activity; however, the impact of ASC on chromatin dynamics during neural stem cell generation remains elusive. Here, we investigate the chromatin changes accompanying neural commitment using an integrative genetics and genomics methodology.

Results

We found that ASC factors bind equally strongly to two distinct classes of cis-regulatory elements: open regions remodeled earlier during maternal to zygotic transition by Zelda and less accessible, Zelda-independent regions. Both classes of cis-elements exhibit enhanced chromatin accessibility during neural specification and correlate with transcriptional regulation of genes involved in a variety of biological processes necessary for neuroblast function/homeostasis. We identified an ASC-Notch regulated TF network that includes likely prime regulators of neuroblast function. Using a cohort of ASC target genes, we report that ASC null neuroblasts are defectively specified, remaining initially stalled, unable to divide, and lacking expression of many proneural targets. When mutant neuroblasts eventually start proliferating, they produce compromised progeny. Reporter lines driven by proneural-bound enhancers display ASC dependency, suggesting that the partial neuroblast identity seen in the absence of ASC genes is likely driven by other, proneural-independent, cis-elements. Neuroblast impairment and the late differentiation defects of ASC mutants are corrected by ectodermal induction of individual ASC genes but not by individual members of the TF network downstream of ASC. However, in wild-type embryos, the induction of individual members of this network induces CNS hyperplasia, suggesting that they synergize with the activating function of ASC to consolidate the chromatin dynamics that promote neural specification.

Conclusions

We demonstrate that ASC proneural transcription factors are indispensable for the timely initiation of the neural stem cell program at the chromatin level by regulating a large number of enhancers in the vicinity of neural genes. This early chromatin remodeling is crucial for both neuroblast homeostasis as well as future progeny fidelity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01300-8.

Single cell enhancer activity distinguishes GABAergic and cholinergic lineages in embryonic mouse basal ganglia

During brain development, neurons are generated by spatially and temporally distinct processes that remain to be fully characterized. The ganglionic eminences (GEs) in the embryonic subpallium give rise to GABAergic and cholinergic neuron lineages that form the basal ganglia or migrate to the cerebral cortex. Beyond a limited set of canonical RNA markers, the transcriptional states of GE progenitors and immature neurons cells remain poorly defined. We combine enhancer labeling, single-cell transcriptomics using transcription factor-anchored clustering, and integration with in situ hybridization data to distinguish emerging neuronal populations in embryonic mouse basal ganglia. Our results demonstrate the specificity of enhancer-based labeling at single-cell resolution and reveal developmental origins and specification processes of critical neuronal lineages.

Enhancers integrate transcription factor signaling pathways that drive cell fate specification in the developing brain. We paired enhancer labeling and single-cell RNA-sequencing (scRNA-seq) to delineate and distinguish specification of neuronal lineages in mouse medial, lateral, and caudal ganglionic eminences (MGE, LGE, and CGE) at embryonic day (E)11.5. We show that scRNA-seq clustering using transcription factors improves resolution of regional and developmental populations, and that enhancer activities identify specific and overlapping GE-derived neuronal populations. First, we mapped the activities of seven evolutionarily conserved brain enhancers at single-cell resolution in vivo, finding that the selected enhancers had diverse activities in specific progenitor and neuronal populations across the GEs. We then applied enhancer-based labeling, scRNA-seq, and analysis of in situ hybridization data to distinguish transcriptionally distinct and spatially defined subtypes of MGE-derived GABAergic and cholinergic projection neurons and interneurons. Our results map developmental origins and specification paths underlying neurogenesis in the embryonic basal ganglia and showcase the power of scRNA-seq combined with enhancer-based labeling to resolve the complex paths of neuronal specification underlying mouse brain development.

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.

•

scATAC time course constructs regulatory trajectories of Drosophila muscle lineages

•

Combining a wild-type trajectory with mutants identifies and quantifies phenotypes de novo

•

Digital nuclear genotyping enables the processing of pooled embryos of mixed genotypes

•

This framework simultaneously uncovers cellular and molecular phenotypes in embryos

scDALI: modeling allelic heterogeneity in single cells reveals context-specific genetic regulation

While it is established that the functional impact of genetic variation can vary across cell types and states, capturing this diversity remains challenging. Current studies using bulk sequencing either ignore this heterogeneity or use sorted cell populations, reducing discovery and explanatory power. Here, we develop scDALI, a versatile computational framework that integrates information on cellular states with allelic quantifications of single-cell sequencing data to characterize cell-state-specific genetic effects. We apply scDALI to scATAC-seq profiles from developing F1 Drosophila embryos and scRNA-seq from differentiating human iPSCs, uncovering heterogeneous genetic effects in specific lineages, developmental stages, or cell types.

DamID transcriptional profiling identifies the Snail/Scratch transcription factor Kahuli as an Alk target in the Drosophila visceral mesoderm

Development of the Drosophila visceral muscle depends on Anaplastic Lymphoma Kinase (Alk) receptor tyrosine kinase (RTK) signaling, which specifies founder cells (FCs) in the circular visceral mesoderm (VM). Although Alk activation by its ligand Jelly Belly (Jeb) is well characterized, few target molecules have been identified. Here, we used targeted DamID (TaDa) to identify Alk targets in embryos overexpressing Jeb versus embryos with abrogated Alk activity, revealing differentially expressed genes, including the Snail/Scratch family transcription factor Kahuli (Kah). We confirmed Kah mRNA and protein expression in the VM, and identified midgut constriction defects in Kah mutants similar to those of pointed (pnt). ChIP and RNA-Seq data analysis defined a Kah target-binding site similar to that of Snail, and identified a set of common target genes putatively regulated by Kah and Pnt during midgut constriction. Taken together, we report a rich dataset of Alk-responsive loci in the embryonic VM and functionally characterize the role of Kah in the regulation of embryonic midgut morphogenesis.

Temporal inhibition of chromatin looping and enhancer accessibility during neuronal remodeling

During development, looping of an enhancer to a promoter is frequently observed in conjunction with temporal and tissue-specific transcriptional activation. The chromatin insulator-associated protein Alan Shepard (Shep) promotes Drosophila post-mitotic neuronal remodeling by repressing transcription of master developmental regulators, such as brain tumor (brat), specifically in maturing neurons. Since insulator proteins can promote looping, we hypothesized that Shep antagonizes brat promoter interaction with an as yet unidentified enhancer. Using chromatin conformation capture and reporter assays, we identified two enhancer regions that increase in looping frequency with the brat promoter specifically in pupal brains after Shep depletion. The brat promoters and enhancers function independently of Shep, ruling out direct repression of these elements. Moreover, ATAC-seq in isolated neurons demonstrates that Shep restricts chromatin accessibility of a key brat enhancer as well as other enhancers genome-wide in remodeling pupal but not larval neurons. These enhancers are enriched for chromatin targets of Shep and are located at Shep-inhibited genes, suggesting direct Shep inhibition of enhancer accessibility and gene expression during neuronal remodeling. Our results provide evidence for temporal regulation of chromatin looping and enhancer accessibility during neuronal maturation.

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

The hourglass model of evolutionary conservation during embryogenesis extends to developmental enhancers with signatures of positive selection

Inter-species comparisons of both morphology and gene expression within a phylum have revealed a period in the middle of embryogenesis with more similarity between species compared to earlier and later time-points. This "developmental hourglass" pattern has been observed in many phyla, yet the evolutionary constraints on gene expression, and underlying mechanisms of how this is regulated, remains elusive. Moreover, the role of positive selection on gene regulation in the more diverged earlier and later stages of embryogenesis remains unknown. Here, using DNase-seq to identify regulatory regions in two distant Drosophila species (D. melanogaster and D. virilis), we assessed the evolutionary conservation and adaptive evolution of enhancers throughout multiple stages of embryogenesis. This revealed a higher proportion of conserved enhancers at the phylotypic period, providing a regulatory basis for the hourglass expression pattern. Using an in silico mutagenesis approach, we detect signatures of positive selection on developmental enhancers at early and late stages of embryogenesis, with a depletion at the phylotypic period, suggesting positive selection as one evolutionary mechanism underlying the hourglass pattern of animal evolution.

Pioneering the developmental frontier

Coordinated changes in gene expression allow a single fertilized oocyte to develop into a complex multi-cellular organism. These changes in expression are controlled by transcription factors that gain access to discrete cis-regulatory elements in the genome, allowing them to activate gene expression. Although nucleosomes present barriers to transcription factor occupancy, pioneer transcription factors have unique properties that allow them to bind DNA in the context of nucleosomes, define cis-regulatory elements, and facilitate the subsequent binding of additional factors that determine gene expression. In this capacity, pioneer factors act at the top of gene-regulatory networks to control developmental transitions. Developmental context also influences pioneer factor binding and activity. Here we discuss the interplay between pioneer factors and development, their role in driving developmental transitions, and the influence of the cellular environment on pioneer factor binding and activity.

Cis-acting variation is common across regulatory layers but is often buffered during embryonic development

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequences, and chromatin. How DNA mutations affecting any one of these regulatory "layers" are buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eight Drosophila crosses at three embryonic stages, yielding a comprehensive data set of 240 samples spanning multiple regulatory layers. Genetic variation (allelic imbalance) impacts gene expression more frequently than chromatin features, with metabolic and environmental response genes being most often affected. Allelic imbalance in cis-regulatory elements (enhancers) is common and highly heritable, yet its functional impact does not generally propagate to gene expression. When it does, genetic variation impacts RNA levels through two alternative mechanisms involving either H3K4me3 or chromatin accessibility and H3K27ac. Changes in RNA are more predictive of variation in H3K4me3 than vice versa, suggesting a role for H3K4me3 downstream from transcription. The impact of a substantial proportion of genetic variation is consistent across embryonic stages, with 50% of allelic imbalanced features at one stage being also imbalanced at subsequent developmental stages. Crucially, buffering, as well as the magnitude and evolutionary impact of genetic variants, is influenced by regulatory complexity (i.e., number of enhancers regulating a gene), with transcription factors being most robust to cis-acting, but most influenced by trans-acting, variation.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.