CATaDa reveals global remodelling of chromatin accessibility during stem cell differentiation in vivo

Mucosa-like differentiation of head and neck cancer cells is inducible and drives the epigenetic loss of cell malignancy

Head and neck squamous cell carcinoma (HNSCC) is a highly malignant disease with high death rates that have remained substantially unaltered for decades. Therefore, new treatment approaches are urgently needed. Human papillomavirus-negative tumors harbor areas of terminally differentiated tissue that are characterized by cornification. Dissecting this intrinsic ability of HNSCC cells to irreversibly differentiate into non-malignant cells may have tumor-targeting potential. We modeled the cornification of HNSCC cells in a primary spheroid model and analyzed the mechanisms underlying differentiation by ATAC-seq and RNA-seq. Results were verified by immunofluorescence using human HNSCC tissue of distinct anatomical locations. HNSCC cell differentiation was accompanied by cell adhesion, proliferation stop, diminished tumor-initiating potential in immunodeficient mice, and activation of a wound-healing-associated signaling program. Small promoter accessibility increased despite overall chromatin closure. Differentiating cells upregulated KRT17 and cornification markers. Although KRT17 represents a basal stem cell marker in normal mucosa, we confirm KRT17 to represent an early differentiation marker in HNSCC tissue. Cornification was frequently found surrounding necrotic areas in human tumors, indicating an involvement of pro-inflammatory stimuli. Indeed, inflammatory mediators activated the differentiation program in primary HNSCC cells. In HNSCC tissue, distinct cell differentiation states were found to create a common tissue architecture in normal mucosa and HNSCCs. Our data demonstrate a loss of cell malignancy upon faithful HNSCC cell differentiation, indicating that targeted differentiation approaches may be therapeutically valuable. Moreover, we describe KRT17 to be a candidate biomarker for HNSCC cell differentiation and early tumor detection.

Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

Notch signaling and Bsh homeodomain activity are integrated to diversify Drosophila lamina neuron types

Notch signaling is an evolutionarily conserved pathway for specifying binary neuronal fates, yet how it specifies different fates in different contexts remains elusive. In our accompanying paper, using the Drosophila lamina neuron types (L1-L5) as a model, we show that the primary homeodomain transcription factor (HDTF) Bsh activates secondary HDTFs Ap (L4) and Pdm3 (L5) and specifies L4/L5 neuronal fates. Here we test the hypothesis that Notch signaling enables Bsh to differentially specify L4 and L5 fates. We show asymmetric Notch signaling between newborn L4 and L5 neurons, but they are not siblings; rather, Notch signaling in L4 is due to Delta expression in adjacent L1 neurons. While Notch signaling and Bsh expression are mutually independent, Notch is necessary and sufficient for Bsh to specify L4 fate over L5. The NotchON L4, compared to NotchOFF L5, has a distinct open chromatin landscape which allows Bsh to bind distinct genomic loci, leading to L4-specific identity gene transcription. We propose a novel model in which Notch signaling is integrated with the primary HDTF activity to diversify neuron types by directly or indirectly generating a distinct open chromatin landscape that constrains the pool of genes that a primary HDTF can activate.

HOPX-associated molecular programs control cardiomyocyte cell states underpinning cardiac structure and function

Genomic regulation of cardiomyocyte differentiation is central to heart development and function. This study uses genetic loss-of-function human-induced pluripotent stem cell-derived cardiomyocytes to evaluate the genomic regulatory basis of the non-DNA-binding homeodomain protein HOPX. We show that HOPX interacts with and controls cardiac genes and enhancer networks associated with diverse aspects of heart development. Using perturbation studies in vitro, we define how upstream cell growth and proliferation control HOPX transcription to regulate cardiac gene programs. We then use cell, organoid, and zebrafish regeneration models to demonstrate that HOPX-regulated gene programs control cardiomyocyte function in development and disease. Collectively, this study mechanistically links cell signaling pathways as upstream regulators of HOPX transcription to control gene programs underpinning cardiomyocyte identity and function.

Homeodomain proteins hierarchically specify neuronal diversity and synaptic connectivity

How our brain generates diverse neuron types that assemble into precise neural circuits remains unclear. Using Drosophila lamina neuron types (L1-L5), we show that the primary homeodomain transcription factor (HDTF) brain-specific homeobox (Bsh) is initiated in progenitors and maintained in L4/L5 neurons to adulthood. Bsh activates secondary HDTFs Ap (L4) and Pdm3 (L5) and specifies L4/L5 neuronal fates while repressing the HDTF Zfh1 to prevent ectopic L1/L3 fates (control: L1-L5; Bsh-knockdown: L1-L3), thereby generating lamina neuronal diversity for normal visual sensitivity. Subsequently, in L4 neurons, Bsh and Ap function in a feed-forward loop to activate the synapse recognition molecule DIP-β, thereby bridging neuronal fate decision to synaptic connectivity. Expression of a Bsh:Dam, specifically in L4, reveals Bsh binding to the DIP-β locus and additional candidate L4 functional identity genes. We propose that HDTFs function hierarchically to coordinate neuronal molecular identity, circuit formation, and function. Hierarchical HDTFs may represent a conserved mechanism for linking neuronal diversity to circuit assembly and function.

Chromatin state transitions in the Drosophila intestinal lineage identify principles of cell-type specification

Tissue homeostasis relies on rewiring of stem cell transcriptional programs into those of differentiated cells. Here, we investigate changes in chromatin occurring in a bipotent adult stem cells. Combining mapping of chromatin-associated factors with statistical modeling, we identify genome-wide transitions during differentiation in the adult Drosophila intestinal stem cell (ISC) lineage. Active, stem-cell-enriched genes transition to a repressive heterochromatin protein-1-enriched state more prominently in enteroendocrine cells (EEs) than in enterocytes (ECs), in which the histone H1-enriched Black state is preeminent. In contrast, terminal differentiation genes associated with metabolic functions follow a common path from a repressive, primed, histone H1-enriched Black state in ISCs to active chromatin states in EE and EC cells. Furthermore, we find that lineage priming has an important function in adult ISCs, and we identify histone H1 as a mediator of this process. These data define underlying principles of chromatin changes during adult multipotent stem cell differentiation.

Methods and applications for single-cell and spatial multi-omics

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Maintenance of neuronal fate and transcriptional identity

The processes that drive naive multipotent stem cells towards fully differentiated fates are increasingly well understood. However, once differentiated, the mechanisms and molecular factors involved in maintaining differentiated states and associated transcriptomes are less well studied. Neurons are a post-mitotic cell-type with highly specialised functions that largely lack the capacity for renewal. Therefore, neuronal cell identities and the transcriptional states that underpin them are locked into place by active mechanisms that prevent lineage reversion/dedifferentiation and repress cell cycling. Furthermore, individual neurons may be very long-lived, so these mechanisms must be sufficient to ensure the fidelity of neuronal transcriptomes over long time periods. This Review aims to provide an overview of recent progress in understanding how neuronal cell fate and associated gene expression are maintained and the transcriptional regulators that are involved. Maintenance of neuronal fate and subtype specification are discussed, as well as the activating and repressive mechanisms involved. The relevance of these processes to disease states, such as brain cancers and neurodegeneration is outlined. Finally, outstanding questions and hypotheses in this field are proposed.

Dam Assisted Fluorescent Tagging of Chromatin Accessibility (DAFCA) for Optical Genome Mapping in Nanochannel Arrays

Proteins and enzymes in the cell nucleus require physical access to their DNA target sites in order to perform genomic tasks such as gene activation and transcription. Hence, chromatin accessibility is a central regulator of gene expression, and its genomic profile holds essential information on the cell type and state. We utilized the E. coli Dam methyltransferase in combination with a fluorescent cofactor analogue to generate fluorescent tags in accessible DNA regions within the cell nucleus. The accessible portions of the genome are then detected by single-molecule optical genome mapping in nanochannel arrays. This method allowed us to characterize long-range structural variations and their associated chromatin structure. We show the ability to create whole-genome, allele-specific chromatin accessibility maps composed of long DNA molecules extended in silicon nanochannels.

Organoid modeling of human fetal lung alveolar development reveals mechanisms of cell fate patterning and neonatal respiratory disease

Variation in lung alveolar development is strongly linked to disease susceptibility. However, underlying cellular and molecular mechanisms are difficult to study in humans. We have identified an alveolar-fated epithelial progenitor in human fetal lungs, which we grow as self-organizing organoids that model key aspects of cell lineage commitment. Using this system, we have functionally validated cell-cell interactions in the developing human alveolar niche, showing that Wnt signaling from differentiating fibroblasts promotes alveolar-type-2 cell identity, whereas myofibroblasts secrete the Wnt inhibitor, NOTUM, providing spatial patterning. We identify a Wnt-NKX2.1 axis controlling alveolar differentiation. Moreover, we show that differential binding of NKX2.1 coordinates alveolar maturation, allowing us to model the effects of human genetic variation in NKX2.1 on alveolar differentiation. Our organoid system recapitulates key aspects of human fetal lung stem cell biology allowing mechanistic experiments to determine the cellular and molecular regulation of human development and disease.

DamID-seq: A Genome-Wide DNA Methylation Method that Captures Both Transient and Stable TF-DNA Interactions in Plant Cells

Capturing the dynamic and transient interactions of a transcription factor (TF) with its genome-wide targets whose regulation leads to plants' adaptation to their changing environment is a major technical challenge. This is a widespread problem with biochemical methods such as chromatin immunoprecipitation-sequencing (ChIP-seq) which are biased towards capturing stable TF-target gene interactions. Herein, we describe how DNA adenine methyltransferase identification and sequencing (DamID-seq) can be used to capture both transient and stable TF-target interactions by DNA methylation. The DamID technique uses a TF protein fused to a DNA adenine methyltransferase (Dam) from E. coli. When expressed in a plant cell, the Dam-TF fusion protein will methylate adenine (A) bases near the sites of TF-DNA interactions. In this way, DamID results in a permanent, stable DNA methylation mark on TF-target gene promoters, even if the target gene is only transiently "touched" by the Dam-TF fusion protein. Here we provide a step-by-step protocol to perform DamID-seq experiments in isolated plant cells for any Dam-TF fusion protein of interest. We also provide information that will enable researchers to analyze DamID-seq data to identify TF-binding sites in the genome. Our protocol includes instructions for vector cloning of the Dam-TF fusion proteins, plant cell protoplast transfections, DamID preps, library preparation, and sequencing data analysis. The protocol outlined in this chapter is performed in Arabidopsis thaliana, however, the DamID-seq workflow developed in this guide is broadly applicable to other plants and organisms.

In vivo, genome-wide profiling of endogenously tagged chromatin-binding proteins with spatial and temporal resolution using NanoDam in Drosophila

NanoDam is a technique for genome-wide profiling of the binding targets of any endogenously tagged chromatin-binding protein in vivo, without the need for antibodies, crosslinking, or immunoprecipitation. Here, we explain the procedure for NanoDam experiments in Drosophila, starting from a genetic cross, to the generation of sequencing libraries and, finally, bioinformatic analysis. This protocol can be readily adapted for use in other model systems after simple modifications. For complete details on the use and execution of this protocol, please refer to Tang et al. (2022).

From developmental to atavistic bet-hedging: How cancer cells pervert the exploitation of random single-cell phenotypic fluctuations

Stochastic gene expression plays a leading developmental role through its contribution to cell differentiation. It is also proposed to promote phenotypic diversification in malignant cells. However, it remains unclear if these two forms of cellular bet-hedging are identical or rather display distinct features. Here we argue that bet-hedging phenomena in cancer cells are more similar to those occurring in unicellular organisms than to those of normal metazoan cells. We further propose that the atavistic bet-hedging strategies in cancer originate from a hijacking of the normal developmental bet-hedging of metazoans. Finally, we discuss the constraints that may shape the atavistic bet-hedging strategies of cancer cells.

Reduced chromatin accessibility correlates with resistance to Notch activation

The Notch signalling pathway is a master regulator of cell fate transitions in development and disease. In the brain, Notch promotes neural stem cell (NSC) proliferation, regulates neuronal migration and maturation and can act as an oncogene or tumour suppressor. How NOTCH and its transcription factor RBPJ activate distinct gene regulatory networks in closely related cell types in vivo remains to be determined. Here we use Targeted DamID (TaDa), requiring only thousands of cells, to identify NOTCH and RBPJ binding in NSCs and their progeny in the mouse embryonic cerebral cortex in vivo. We find that NOTCH and RBPJ associate with a broad network of NSC genes. Repression of NSC-specific Notch target genes in intermediate progenitors and neurons correlates with decreased chromatin accessibility, suggesting that chromatin compaction may contribute to restricting NOTCH-mediated transactivation.

Genome-wide maps of nucleolus interactions reveal distinct layers of repressive chromatin domains

Eukaryotic chromosomes are folded into hierarchical domains, forming functional compartments. Nuclear periphery and nucleolus are two nuclear landmarks contributing to repressive chromosome architecture. However, while the role of nuclear lamina (NL) in genome organization has been well documented, the function of the nucleolus remains under-investigated due to the lack of methods for the identification of nucleolar associated domains (NADs). Here we have established DamID- and HiC-based methodologies to generate accurate genome-wide maps of NADs in embryonic stem cells (ESCs) and neural progenitor cells (NPCs), revealing layers of genome compartmentalization with distinct, repressive chromatin states based on the interaction with the nucleolus, NL, or both. NADs show higher H3K9me2 and lower H3K27me3 content than regions exclusively interacting with NL. Upon ESC differentiation into NPCs, chromosomes around the nucleolus acquire a more compact, rigid architecture with neural genes moving away from nucleoli and becoming unlocked for later activation. Further, histone modifications and the interaction strength within A and B compartments of NADs and LADs in ESCs set the choice to associate with NL or nucleoli upon dissociation from their respective compartments during differentiation. The methodologies here developed will make possible to include the nucleolar contribution in nuclear space and genome function in diverse biological systems.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

OUP accepted manuscript

Eukaryotic chromosomes are spatially segregated into topologically associating domains (TADs). Some TADs are attached to the nuclear lamina (NL) through lamina-associated domains (LADs). Here, we identified LADs and TADs at two stages of Drosophila spermatogenesis – in bam^Δ86 mutant testes which is the commonly used model of spermatogonia (SpG) and in larval testes mainly filled with spermatocytes (SpCs). We found that initiation of SpC-specific transcription correlates with promoters’ detachment from the NL and with local spatial insulation of adjacent regions. However, this insulation does not result in the partitioning of inactive TADs into sub-TADs. We also revealed an increased contact frequency between SpC-specific genes in SpCs implying their de novo gathering into transcription factories. In addition, we uncovered the specific X chromosome organization in the male germline. In SpG and SpCs, a single X chromosome is stronger associated with the NL than autosomes. Nevertheless, active chromatin regions in the X chromosome interact with each other more frequently than in autosomes. Moreover, despite the absence of dosage compensation complex in the male germline, randomly inserted SpG-specific reporter is expressed higher in the X chromosome than in autosomes, thus evidencing that non-canonical dosage compensation operates in SpG.

Single-Molecule Multikilobase-Scale Profiling of Chromatin Accessibility Using m6A-SMAC-Seq and m6A-CpG-GpC-SMAC-Seq

A hallmark feature of active cis-regulatory elements (CREs) in eukaryotes is their nucleosomal depletion and, accordingly, higher accessibility to enzymatic treatment. This property has been the basis of a number of sequencing-based assays for genome-wide identification and tracking the activity of CREs across different biological conditions, such as DNAse-seq, ATAC-seq , NOMeseq, and others. However, the fragmentation of DNA inherent to many of these assays and the limited read length of short-read sequencing platforms have so far not allowed the simultaneous measurement of the chromatin accessibility state of CREs located distally from each other. The combination of labeling accessible DNA with DNA modifications and nanopore sequencing has made it possible to develop such assays. Here, we provide a detailed protocol for carrying out the SMAC-seq assay (Single-Molecule long-read Accessible Chromatin mapping sequencing), in its m6A-SMAC-seq and m6A-CpG-GpC-SMAC-seq variants, together with methods for data processing and analysis, and discuss key experimental and analytical considerations for working with SMAC-seq datasets.

Profiling Protein-DNA Interactions Cell-Type-Specifically with Targeted DamID

Targeted DamID (TaDa) is a means of profiling the binding of any DNA-associated protein cell-type specifically, including transcription factors, RNA polymerase, and chromatin-modifying proteins. The technique is highly sensitive, highly reproducible, requires no mechanical disruption, cell isolation or antibody purification, and can be performed by anyone with basic molecular biology knowledge. Here, we describe the TaDa method and downstream bioinformatics data processing.

In vivo targeted DamID identifies CHD8 genomic targets in fetal mouse brain

Genetic studies of autism have revealed causal roles for chromatin remodeling gene mutations. Chromodomain helicase DNA binding protein 8 (CHD8) encodes a chromatin remodeler with significant de novo mutation rates in sporadic autism. However, relationships between CHD8 genomic function and autism-relevant biology remain poorly elucidated. Published studies utilizing ChIP-seq to map CHD8 protein-DNA interactions have high variability, consistent with technical challenges and limitations associated with this method. Thus, complementary approaches are needed to establish CHD8 genomic targets and regulatory functions in developing brain. We used in utero CHD8 Targeted DamID followed by sequencing (TaDa-seq) to characterize CHD8 binding in embryonic mouse cortex. CHD8 TaDa-seq reproduced interaction patterns observed from ChIP-seq and further highlighted CHD8 distal interactions associated with neuronal loci. This study establishes TaDa-seq as a useful alternative for mapping protein-DNA interactions in vivo and provides insights into the regulatory targets of CHD8 and autism-relevant pathophysiology associated with CHD8 mutations.

Olfactory Receptor Gene Regulation in Insects: Multiple Mechanisms for Singular Expression

The singular expression of insect olfactory receptors in specific populations of olfactory sensory neurons is fundamental to the encoding of odors in patterns of neuronal activity in the brain. How a receptor gene is selected, from among a large repertoire in the genome, to be expressed in a particular neuron is an outstanding question. Focusing on Drosophila melanogaster, where most investigations have been performed, but incorporating recent insights from other insect species, we review the multilevel regulatory mechanisms of olfactory receptor expression. We discuss how cis-regulatory elements, trans-acting factors, chromatin modifications, and feedback pathways collaborate to activate and maintain expression of the chosen receptor (and to suppress others), highlighting similarities and differences with the mechanisms underlying singular receptor expression in mammals. We also consider the plasticity of receptor regulation in response to environmental cues and internal state during the lifetime of an individual, as well as the evolution of novel expression patterns over longer timescales. Finally, we describe the mechanisms and potential significance of examples of receptor co-expression.

Tissue-specific DamID protocol using nanopore sequencing

DNA adenine methylation identification (DamID) is a powerful method to determine DNA binding profiles of proteins at a genomic scale. The method leverages the fusion between a protein of interest and the Dam methyltransferase of E. coli, which methylates proximal DNA in vivo. Here, we present an optimized procedure, which was developed for tissue-specific analyses in Caenorhabditis elegans and successfully used to footprint genes actively transcribed by RNA polymerases and to map transcription factor binding in gene regulatory regions. The present protocol details C. elegans-specific steps involved in the preparation of transgenic lines and genomic DNA samples, as well as broadly applicable steps for the DamID procedure, including the isolation of methylated DNA fragments, the preparation of multiplexed libraries, Nanopore sequencing, and data analysis. Two distinctive features of the approach are (i) the use of an efficient recombination-based strategy to selectively analyze rare cell types and (ii) the use of Nanopore sequencing, which streamlines the process. The method allows researchers to go from genomic DNA samples to sequencing results in less than a week, while being sensitive enough to report reliable DNA footprints in cell types as rare as 2 cells per animal.

Quantification, Dynamic Visualization, and Validation of Bias in ATAC-Seq Data with ataqv

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) has become the preferred method for mapping chromatin accessibility due to its time and input material efficiency. However, it can be difficult to evaluate data quality and identify sources of technical bias across samples. Here, we present ataqv, a computational toolkit for efficiently measuring, visualizing, and comparing quality control (QC) results across samples and experiments. We use ataqv to analyze 2,009 public ATAC-seq datasets; their QC metrics display a 10-fold range. Tn5 dosage experiments and statistical modeling show that technical variation in the ratio of Tn5 transposase to nuclei and sequencing flowcell density induces systematic bias in ATAC-seq data by changing the enrichment of reads across functional genomic annotations including promoters, enhancers, and transcription-factor-bound regions, with the notable exception of CTCF. ataqv can be integrated into existing computational pipelines and is freely available at https://github.com/ParkerLab/ataqv/.

A KDM5-Prospero transcriptional axis functions during early neurodevelopment to regulate mushroom body formation

Mutations in the lysine demethylase 5 (KDM5) family of transcriptional regulators are associated with intellectual disability, yet little is known regarding their spatiotemporal requirements or neurodevelopmental contributions. Utilizing the mushroom body (MB), a major learning and memory center within the Drosophila brain, we demonstrate that KDM5 is required within ganglion mother cells and immature neurons for proper axogenesis. Moreover, the mechanism by which KDM5 functions in this context is independent of its canonical histone demethylase activity. Using in vivo transcriptional and binding analyses, we identify a network of genes directly regulated by KDM5 that are critical modulators of neurodevelopment. We find that KDM5 directly regulates the expression of prospero, a transcription factor that we demonstrate is essential for MB morphogenesis. Prospero functions downstream of KDM5 and binds to approximately half of KDM5-regulated genes. Together, our data provide evidence for a KDM5-Prospero transcriptional axis that is essential for proper MB development.

Targeted molecular profiling of rare olfactory sensory neurons identifies fate, wiring, and functional determinants

Determining the molecular properties of neurons is essential to understand their development, function and evolution. Using Targeted DamID (TaDa), we characterize RNA polymerase II occupancy and chromatin accessibility in selected Ionotropic receptor (Ir)-expressing olfactory sensory neurons in Drosophila. Although individual populations represent a minute fraction of cells, TaDa is sufficiently sensitive and specific to identify the expected receptor genes. Unique Ir expression is not consistently associated with differences in chromatin accessibility, but rather to distinct transcription factor profiles. Genes that are heterogeneously expressed across populations are enriched for neurodevelopmental factors, and we identify functions for the POU-domain protein Pdm3 as a genetic switch of Ir neuron fate, and the atypical cadherin Flamingo in segregation of neurons into discrete glomeruli. Together this study reveals the effectiveness of TaDa in profiling rare neural populations, identifies new roles for a transcription factor and a neuronal guidance molecule, and provides valuable datasets for future exploration.

Mechanisms underlying the control of dynamic regulatory element activity and chromatin accessibility during metamorphosis

Cis-regulatory modules of metazoan genomes determine the when and where of gene expression during development. Here we discuss insights into the genetic and molecular mechanisms behind cis-regulatory module usage that have come from recent application of genomics assays to insect metamorphosis. Assays including FAIRE-seq, ATAC-seq, and CUT&RUN indicate that sequential changes in chromatin accessibility play a key role in mediating stage-specific cis-regulatory module activity and gene expression. We review the current understanding of what controls precisely coordinated changes in chromatin accessibility during metamorphosis and describe evidence that points to systemic hormone signaling as a primary signal to trigger genome-wide shifts in accessibility patterns and cis-regulatory module usage.

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.

Genetic approaches to revealing the principles of nuclear architecture

The spatial organization of chromosomes inside the eukaryotic nucleus is important for DNA replication, repair and gene expression. During development of multicellular organisms, different compendiums of genes are either repressed or activated to produce specific cell types. Genetic manipulation of tractable organisms is invaluable to elucidate chromosome configuration and the underlying mechanisms. Systematic inhibition of genes through RNA interference and, more recently, CRISPR/Cas9-based screens have identified new proteins with significant roles in nuclear organization. Coupling this with advances in imaging techniques, such as multiplexed DNA fluorescence in situ hybridization, and with tissue-specific genome profiling by DNA adenine methylation identification has increased our knowledge about the immense complexity and dynamics of the nucleus.

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.

Persistent epigenetic reprogramming of sweet taste by diet

Diets rich in sugar, salt, and fat alter taste perception and food preference, contributing to obesity and metabolic disorders, but the molecular mechanisms through which this occurs are unknown. Here, we show that in response to a high sugar diet, the epigenetic regulator Polycomb Repressive Complex 2.1 (PRC2.1) persistently reprograms the sensory neurons of Drosophila melanogaster flies to reduce sweet sensation and promote obesity. In animals fed high sugar, the binding of PRC2.1 to the chromatin of the sweet gustatory neurons is redistributed to repress a developmental transcriptional network that modulates the responsiveness of these cells to sweet stimuli, reducing sweet sensation. Half of these transcriptional changes persist despite returning the animals to a control diet, causing a permanent decrease in sweet taste. Our results uncover a new epigenetic mechanism that, in response to the dietary environment, regulates neural plasticity and feeding behavior to promote obesity.

μDamID: A Microfluidic Approach for Joint Imaging and Sequencing of Protein-DNA Interactions in Single Cells

DNA adenine methyltransferase identification (DamID) measures a protein's DNA-binding history by methylating adenine bases near each protein-DNA interaction site and then selectively amplifying and sequencing these methylated regions. Additionally, these interactions can be visualized using m6A-Tracer, a fluorescent protein that binds to methyladenines. Here, we combine these imaging and sequencing technologies in an integrated microfluidic platform (μDamID) that enables single-cell isolation, imaging, and sorting, followed by DamID. We use μDamID and an improved m6A-Tracer protein to generate paired imaging and sequencing data from individual human cells. We validate interactions between Lamin-B1 protein and lamina-associated domains (LADs), observe variable 3D chromatin organization and broad gene regulation patterns, and jointly measure single-cell heterogeneity in Dam expression and background methylation. μDamID provides the unique ability to compare paired imaging and sequencing data for each cell and between cells, enabling the joint analysis of the nuclear localization, sequence identity, and variability of protein-DNA interactions. A record of this paper's transparent peer review process is included in the Supplemental Information.

Tissue-Specific Transcription Footprinting Using RNA PoI DamID (RAPID) in Caenorhabditis elegans

Differential gene expression across cell types underlies development and cell physiology in multicellular organisms. Caenorhabditis elegans is a powerful, extensively used model to address these biological questions. A remaining bottleneck relates to the difficulty to obtain comprehensive tissue-specific gene transcription data, since available methods are still challenging to execute and/or require large worm populations. Here, we introduce the RNA Polymerase DamID (RAPID) approach, in which the Dam methyltransferase is fused to a ubiquitous RNA polymerase subunit to create transcriptional footprints via methyl marks on the DNA of transcribed genes. To validate the method, we determined the polymerase footprints in whole animals, in sorted embryonic blastomeres and in different tissues from intact young adults by driving tissue-specific Dam fusion expression. We obtained meaningful transcriptional footprints in line with RNA-sequencing (RNA-seq) studies in whole animals or specific tissues. To challenge the sensitivity of RAPID and demonstrate its utility to determine novel tissue-specific transcriptional profiles, we determined the transcriptional footprints of the pair of XXX neuroendocrine cells, representing 0.2% of the somatic cell content of the animals. We identified 3901 candidate genes with putatively active transcription in XXX cells, including the few previously known markers for these cells. Using transcriptional reporters for a subset of new hits, we confirmed that the majority of them were expressed in XXX cells and identified novel XXX-specific markers. Taken together, our work establishes RAPID as a valid method for the determination of RNA polymerase footprints in specific tissues of C. elegans without the need for cell sorting or RNA tagging.

Dynamic neurotransmitter specific transcription factor expression profiles during Drosophila development

The remarkable diversity of neurons in the nervous system is generated during development, when properties such as cell morphology, receptor profiles and neurotransmitter identities are specified. In order to gain a greater understanding of neurotransmitter specification we profiled the transcription state of cholinergic, GABAergic and glutamatergic neurons in vivo at three developmental time points. We identified 86 differentially expressed transcription factors that are uniquely enriched, or uniquely depleted, in a specific neurotransmitter type. Some transcription factors show a similar profile across development, others only show enrichment or depletion at specific developmental stages. Profiling of Acj6 (cholinergic enriched) and Ets65A (cholinergic depleted) binding sites in vivo reveals that they both directly bind the ChAT locus, in addition to a wide spectrum of other key neuronal differentiation genes. We also show that cholinergic enriched transcription factors are expressed in mostly non-overlapping populations in the adult brain, implying the absence of combinatorial regulation of neurotransmitter fate in this context. Furthermore, our data underlines that, similar to Caenorhabditis elegans, there are no simple transcription factor codes for neurotransmitter type specification.

This article has an associated First Person interview with the first author of the paper.

Summary: Transcriptome profiling of cholinergic, GABAergic and glutamatergic neurons in Drosophila identified multiple transcription factors as potential regulators of neurotransmitter fate.

Are Parallel Proliferation Pathways Redundant?

Are the receptor tyrosine kinase (RTK) and JAK-STAT-driven proliferation pathways 'parallel' or 'redundant'? And what about those of K-Ras4B versus N-Ras? 'Parallel' proliferation pathways accomplish a similar drug resistance outcome. Thus, are they 'redundant'? In this paper, it is argued that there is a fundamental distinction between 'parallel' and 'redundant'. Cellular proliferation pathways are influenced by the genome sequence, 3D organization and chromatin accessibility, and determined by protein availability prior to cancer emergence. In the opinion presented, if they operate the same downstream protein families, they are redundant; if evolutionary-independent, they are parallel. Thus, RTK and JAK-STAT-driven proliferation pathways are parallel; those of Ras isoforms are redundant. Our Precision Medicine Call to map cancer proliferation pathways is vastly important since it can expedite effective therapeutics.

Single-cell transcriptional diversity is a hallmark of developmental potential

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for reconstructing cellular differentiation trajectories. However, inferring both the state and direction of differentiation is challenging. Here, we demonstrate a simple, yet robust, determinant of developmental potential-the number of expressed genes per cell-and leverage this measure of transcriptional diversity to develop a computational framework (CytoTRACE) for predicting differentiation states from scRNA-seq data. When applied to diverse tissue types and organisms, CytoTRACE outperformed previous methods and nearly 19,000 annotated gene sets for resolving 52 experimentally determined developmental trajectories. Additionally, it facilitated the identification of quiescent stem cells and revealed genes that contribute to breast tumorigenesis. This study thus establishes a key RNA-based feature of developmental potential and a platform for delineation of cellular hierarchies.

Condensin I subunit Cap-G is essential for proper gene expression during the maturation of post-mitotic neurons

Condensin complexes are essential for mitotic chromosome assembly and segregation during cell divisions, however, little is known about their functions in post-mitotic cells. Here we report a role for the condensin I subunit Cap-G in Drosophila neurons. We show that, despite not requiring condensin for mitotic chromosome compaction, post-mitotic neurons express Cap-G. Knockdown of Cap-G specifically in neurons (from their birth onwards) results in developmental arrest, behavioural defects, and dramatic gene expression changes, including reduced expression of a subset of neuronal genes and aberrant expression of genes that are not normally expressed in the developing brain. Knockdown of Cap-G in mature neurons results in similar phenotypes but to a lesser degree. Furthermore, we see dynamic binding of Cap-G at distinct loci in progenitor cells and differentiated neurons. Therefore, Cap-G is essential for proper gene expression in neurons and plays an important role during the early stages of neuronal development.

Dynamics of activating and repressive histone modifications in Drosophila neural stem cell lineages and brain tumors

During central nervous system development, spatiotemporal gene expression programs mediate specific lineage decisions to generate neuronal and glial cell types from neural stem cells (NSCs). However, little is known about the epigenetic landscape underlying these highly complex developmental events. Here, we perform ChIP-seq on distinct subtypes of Drosophila FACS-purified NSCs and their differentiated progeny to dissect the epigenetic changes accompanying the major lineage decisions in vivo By analyzing active and repressive histone modifications, we show that stem cell identity genes are silenced during differentiation by loss of their activating marks and not via repressive histone modifications. Our analysis also uncovers a new set of genes specifically required for altering lineage patterns in type II neuroblasts (NBs), one of the two main Drosophila NSC identities. Finally, we demonstrate that this subtype specification in NBs, unlike NSC differentiation, requires Polycomb-group-mediated repression.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

MadID, a Versatile Approach to Map Protein-DNA Interactions, Highlights Telomere-Nuclear Envelope Contact Sites in Human Cells

Mapping the binding sites of DNA- or chromatin-interacting proteins is essential to understanding biological processes. DNA adenine methyltransferase identification (DamID) has emerged as a comprehensive method to map genome-wide occupancy of proteins of interest. A caveat of DamID is the specificity of Dam methyltransferase for GATC motifs that are not homogenously distributed in the genome. Here, we developed an optimized method named MadID, using proximity labeling of DNA by the methyltransferase M.EcoGII. M.EcoGII mediates N6-adenosine methylation in any DNA sequence context, resulting in deeper and unbiased coverage of the genome. We demonstrate, using m6A-specific immunoprecipitation and deep sequencing, that MadID is a robust method to identify protein-DNA interactions at the whole-genome level. Using MadID, we revealed contact sites between human telomeres, repetitive sequences devoid of GATC sites, and the nuclear envelope. Overall, MadID opens the way to identification of binding sites in genomic regions that were largely inaccessible.

•

MadID: mapping of protein-DNA interactions in vivo using proximity labeling

•

Deeper and unbiased genome-wide coverage using M.EcoGII, a methyltransferase

•

Identification of binding sites in previously inaccessible regions of the genome

•

Identification of telomere-nuclear envelope contact sites

Simultaneous quantification of protein-DNA contacts and transcriptomes in single cells

Protein-DNA interactions are critical to the regulation of gene expression, but it remains challenging to define how cell-to-cell heterogeneity in protein-DNA binding influences gene expression variability. Here we report a method for the simultaneous quantification of protein-DNA contacts by combining single-cell DNA adenine methyltransferase identification (DamID) with messenger RNA sequencing of the same cell (scDam&T-seq). We apply scDam&T-seq to reveal how genome-lamina contacts or chromatin accessibility correlate with gene expression in individual cells. Furthermore, we provide single-cell genome-wide interaction data on a polycomb-group protein, RING1B, and the associated transcriptome. Our results show that scDam&T-seq is sensitive enough to distinguish mouse embryonic stem cells cultured under different conditions and their different chromatin landscapes. Our method will enable the analysis of protein-mediated mechanisms that regulate cell-type-specific transcriptional programs in heterogeneous tissues.

DamID as a versatile tool for understanding gene regulation

The interaction of proteins and RNA with chromatin underlies the regulation of gene expression. The ability to profile easily these interactions is fundamental for understanding chromatin biology in vivo. DNA adenine methyltransferase identification (DamID) profiles genome-wide protein-DNA interactions without antibodies, fixation or protein pull-downs. Recently, DamID has been adapted for applications beyond simple assaying of protein-DNA interactions, such as for studying RNA-chromatin interactions, chromatin accessibility and long-range chromosome interactions. Here, we provide an overview of DamID and introduce improvements to the technology, discuss their applications and compare alternative methodologies.

Summary: This Primer provides an overview of DNA adenine methyltransferase identification (DamID), which is used to profile genome-wide chromatin interactions, and introduces recent improvements to the technology.

TaDa! Analysing cell type-specific chromatin in vivo with Targeted DamID

The emergence of neuronal diversity during development of the nervous system relies on dynamic changes in the epigenetic landscape of neural stem cells and their progeny. Targeted DamID (TaDa) is proving invaluable in identifying the genome-wide binding sites of chromatin-associated proteins in vivo, without fixation, cell isolation, or immunoprecipitation. The simplicity and efficiency of the technique have led to an ever-expanding TaDa toolbox. These tools enable profiling of gene expression and chromatin accessibility, as well as the identification of the genome-wide binding sites of chromatin complexes, transcription factors and RNAs. Here, we review these new developments, with particular emphasis on the use of TaDa in studying neuronal specification.

Neuroblast-specific open chromatin allows the temporal transcription factor, Hunchback, to bind neuroblast-specific loci

Spatial and temporal cues are required to specify neuronal diversity, but how these cues are integrated in neural progenitors remains unknown. Drosophila progenitors (neuroblasts) are a good model: they are individually identifiable with relevant spatial and temporal transcription factors known. Here we test whether spatial/temporal factors act independently or sequentially in neuroblasts. We used Targeted DamID to identify genomic binding sites of the Hunchback temporal factor in two neuroblasts (NB5-6 and NB7-4) that make different progeny. Hunchback targets were different in each neuroblast, ruling out the independent specification model. Moreover, each neuroblast had distinct open chromatin domains, which correlated with differential Hb-bound loci in each neuroblast. Importantly, the Gsb/Pax3 spatial factor, expressed in NB5-6 but not NB7-4, had genomic binding sites correlated with open chromatin in NB5-6, but not NB7-4. Our data support a model in which early-acting spatial factors like Gsb establish neuroblast-specific open chromatin domains, leading to neuroblast-specific temporal factor binding and the production of different neurons in each neuroblast lineage.

The human brain is considered to be the most complicated object in the universe, but it only takes a handful of stem cells to make one. The process depends on two types of information: signals separated across space and time. Spatial cues tell a stem cell what type of cell it is going to be, while temporal cues work as molecular clocks to generate a sequence of different neurons over time. Together, these cues generate the large array of cell types in the nervous system.

Each stem cell occupies its own space in the developing body and receives its own spatial cues, but they all follow the same timeline. For example, proteins called transcription factors act as molecular clocks and interact with specific genes, telling the cell when to turn them on or off. The same series of transcription factors operates in different stem cells, but they have different effects. So far, it has been unclear whether spatial and temporal signals work independently or sequentially to generate new cell types.

To find out, Sen et al. studied two distinct, developing stem cells in fruit flies, which receive different spatial signals. Transcription factors only work if they are able to get to their target genes. Cells can open or close access to different genes by changing the structure of the chromatin wrapping that surrounds the genes. In the experiments, a marker was used to reveal the areas of open chromatin in each of the cells. Another marker was used to track the transcription factors. The results showed that the areas of open chromatin varied between stem cells. Moreover, although both cells used the same transcription factor called Hunchback, it targeted different genes in each stem cell. This was due to changes in the chromatin wrapping: Hunchback only acted in areas where the chromatin was open. This suggests that the spatial cues first sculpt the chromatin, making some genes easier to get to than others. Then, the same transcription factors go to the accessible gene, which will differ from one stem cell to another.

These findings help us to understand how different types of brain cells develop, which may also aid us in finding a way how to engineer specific cell types. If we could turn stem cells into different types of brain cells, it might help us to treat brain diseases. This may involve giving the right spatial signal before starting the temporal cues.

Targeted DamID reveals differential binding of mammalian pluripotency factors

The precise control of gene expression by transcription factor networks is crucial to organismal development. The predominant approach for mapping transcription factor-chromatin interactions has been chromatin immunoprecipitation (ChIP). However, ChIP requires a large number of homogeneous cells and antisera with high specificity. A second approach, DamID, has the drawback that high levels of Dam methylase are toxic. Here, we modify our targeted DamID approach (TaDa) to enable cell type-specific expression in mammalian systems, generating an inducible system (mammalian TaDa or MaTaDa) to identify genome-wide protein/DNA interactions in 100 to 1000 times fewer cells than ChIP-based approaches. We mapped the binding sites of two key pluripotency factors, OCT4 and PRDM14, in mouse embryonic stem cells, epiblast-like cells and primordial germ cell-like cells (PGCLCs). PGCLCs are an important system for elucidating primordial germ cell development in mice. We monitored PRDM14 binding during the specification of PGCLCs, identifying direct targets of PRDM14 that are key to understanding its crucial role in PGCLC development. We show that MaTaDa is a sensitive and accurate method for assessing cell type-specific transcription factor binding in limited numbers of cells.