Dual topologies of myotomal collagen XV and Tenascin C act in concert to guide and shape developing motor axons.

During development, motor axons are guided toward muscle target by various extrinsic cues including extracellular matrix (ECM) proteins whose identities and cellular source remain poorly characterized. Here, using single-cell RNAseq of sorted GFP+ cells from smyhc1:gfp-injected zebrafish embryos, we unravel the slow muscle progenitors (SMP) pseudotemporal trajectory at the single-cell level and show that differentiating SMPs are a major source of ECM proteins. The SMP core-matrisome was characterized and computationally predicted to form a basement membrane-like structure tailored for motor axon guidance, including basement membrane-associated ECM proteins, as collagen XV-B, one of the earliest core-matrisome gene transcribed in differentiating SMPs and the glycoprotein Tenascin C. To investigate how contact-mediated guidance cues are organized along the motor path to exert their function in vivo, we used microscopy-based methods to analyze and quantify motor axon navigation in tnc and col15a1b knock-out fish. We show that motor axon shape and growth rely on the timely expression of the attractive cue Collagen XV-B that locally provides axons with a permissive soft microenvironment and separately organizes the repulsive cue Tenascin C into a unique functional dual topology. Importantly, bioprinted micropatterns that mimic this in vivo ECM topology were sufficient to drive directional motor axon growth. Our study offers evidence that not only the composition of ECM cues but their topology critically influences motor axon navigation in vertebrates with potential applications in regenerative medicine for peripheral nerve injury as regenerating nerves follow their original path.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries.

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

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.

Sequential actions of EOMES and T-BET promote stepwise maturation of natural killer cells.

EOMES and T-BET are related T-box transcription factors that control natural killer (NK) cell development. Here we demonstrate that EOMES and T-BET regulate largely distinct gene sets during this process. EOMES is dominantly expressed in immature NK cells and drives early lineage specification by inducing hallmark receptors and functions. By contrast, T-BET is dominant in mature NK cells, where it induces responsiveness to IL-12 and represses the cell cycle, likely through transcriptional repressors. Regardless, many genes with distinct functions are co-regulated by the two transcription factors. By generating two gene-modified mice facilitating chromatin immunoprecipitation of endogenous EOMES and T-BET, we show a strong overlap in their DNA binding targets, as well as extensive epigenetic changes during NK cell differentiation. Our data thus suggest that EOMES and T-BET may distinctly govern, via differential expression and co-factors recruitment, NK cell maturation by inserting partially overlapping epigenetic regulations.

m6A RNA methylation regulates promoter- proximal pausing of RNA polymerase II.

RNA polymerase II (RNAP II) pausing is essential to precisely control gene expression and is critical for development of metazoans. Here, we show that the m6A RNA modification regulates promoter-proximal RNAP II pausing in Drosophila cells. The m6A methyltransferase complex (MTC) and the nuclear reader Ythdc1 are recruited to gene promoters. Depleting the m6A MTC leads to a decrease in RNAP II pause release and in Ser2P occupancy on the gene body and affects nascent RNA transcription. Tethering Mettl3 to a heterologous gene promoter is sufficient to increase RNAP II pause release, an effect that relies on its m6A catalytic domain. Collectively, our data reveal an important link between RNAP II pausing and the m6A RNA modification, thus adding another layer to m6A-mediated gene regulation.

Zeb1 represses TCR signaling, promotes the proliferation of T cell progenitors and is essential for NK1.1+ T cell development.

T cell development proceeds under the influence of a network of transcription factors (TFs). The precise role of Zeb1, a member of this network, remains unclear. Here, we report that Zeb1 expression is induced early during T cell development in CD4-CD8- double-negative (DN) stage 2 (DN2). Zeb1 expression was further increased in the CD4+CD8+ double-positive (DP) stage before decreasing in more mature T cell subsets. We performed an exhaustive characterization of T cells in Cellophane mice that bear Zeb1 hypomorphic mutations. The Zeb1 mutation profoundly affected all thymic subsets, especially DN2 and DP cells. Zeb1 promoted the survival and proliferation of both cell populations in a cell-intrinsic manner. In the periphery of Cellophane mice, the number of conventional T cells was near normal, but invariant NKT cells, NK1.1+ γδ T cells and Ly49+ CD8 T cells were virtually absent. This suggested that Zeb1 regulates the development of unconventional T cell types from DP progenitors. A transcriptomic analysis of WT and Cellophane DP cells revealed that Zeb1 regulated the expression of multiple genes involved in the cell cycle and TCR signaling, which possibly occurred in cooperation with Tcf1 and Heb. Indeed, Cellophane DP cells displayed stronger signaling than WT DP cells upon TCR engagement in terms of the calcium response, phosphorylation events, and expression of early genes. Thus, Zeb1 is a key regulator of the cell cycle and TCR signaling during thymic T cell development. We propose that thymocyte selection is perturbed in Zeb1-mutated mice in a way that does not allow the survival of unconventional T cell subsets.

3D genome organisation in Drosophila.

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

Functional Consequences of Chromosomal Rearrangements on Gene Expression: Not So Deleterious After All?

Chromosomes are folded and organized into topologically associating domains (TADs) which provide a framework for the interaction of enhancers with the promoter of their target gene(s). Structural rearrangements observed during evolution or in disease contexts suggest that changes in genome organization strongly affect gene expression and can have drastic phenotypic effects. In this review, I will discuss how recent genomic engineering experiments reveal a more contrasted picture, suggesting that TADs are important but not always essential for gene expression regulation.

Highly rearranged chromosomes reveal uncoupling between genome topology and gene expression.

Chromatin topology is intricately linked to gene expression, yet its functional requirement remains unclear. Here, we comprehensively assessed the interplay between genome topology and gene expression using highly rearranged chromosomes (balancers) spanning ~75% of the Drosophila genome. Using transheterozyte (balancer/wild-type) embryos, we measured allele-specific changes in topology and gene expression in cis, while minimizing trans effects. Through genome sequencing, we resolved eight large nested inversions, smaller inversions, duplications and thousands of deletions. These extensive rearrangements caused many changes to chromatin topology, disrupting long-range loops, topologically associating domains (TADs) and promoter interactions, yet these are not predictive of changes in expression. Gene expression is generally not altered around inversion breakpoints, indicating that mis-appropriate enhancer-promoter activation is a rare event. Similarly, shuffling or fusing TADs, changing intra-TAD connections and disrupting long-range inter-TAD loops does not alter expression for the majority of genes. Our results suggest that properties other than chromatin topology ensure productive enhancer-promoter interactions.

Chromosome topology guides the Drosophila Dosage Compensation Complex for target gene activation.

X chromosome dosage compensation in Drosophila requires chromosome-wide coordination of gene activation. The male-specific lethal dosage compensation complex (DCC) identifies and binds to X-chromosomal high-affinity sites (HAS) from which it boosts transcription. A sub-class of HAS, PionX sites, represent first contacts on the X. Here, we explored the chromosomal interactions of representative PionX sites by high-resolution 4C and determined the global chromosome conformation by Hi-C in sex-sorted embryos. Male and female X chromosomes display similar nuclear architecture, concordant with clustered, constitutively active genes. PionX sites, like HAS, are evenly distributed in the active compartment and engage in short- and long-range interactions beyond compartment boundaries. Long-range, inter-domain interactions between DCC binding sites are stronger in males, suggesting that the complex refines chromatin organization. By de novo induction of DCC in female cells, we monitored the extent of activation surrounding PionX sites. This revealed a remarkable range of DCC action not only in linear proximity, but also at megabase distance if close in space, suggesting that DCC profits from pre-existing chromosome folding to activate genes.

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

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

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

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

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

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

FourCSeq: analysis of 4C sequencing data.

MOTIVATION: Circularized Chromosome Conformation Capture (4C) is a powerful technique for studying the spatial interactions of a specific genomic region called the 'viewpoint' with the rest of the genome, both in a single condition or comparing different experimental conditions or cell types. Observed ligation frequencies typically show a strong, regular dependence on genomic distance from the viewpoint, on top of which specific interaction peaks are superimposed. Here, we address the computational task to find these specific peaks and to detect changes between different biological conditions. RESULTS: We model the overall trend of decreasing interaction frequency with genomic distance by fitting a smooth monotonically decreasing function to suitably transformed count data. Based on the fit, z-scores are calculated from the residuals, and high z-scores are interpreted as peaks providing evidence for specific interactions. To compare different conditions, we normalize fragment counts between samples, and call for differential contact frequencies using the statistical method DESEQ2: adapted from RNA-Seq analysis. AVAILABILITY AND IMPLEMENTATION: A full end-to-end analysis pipeline is implemented in the R package FourCSeq available at www.bioconductor.org. CONTACT: felix.klein@embl.de or whuber@embl.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enhancer loops appear stable during development and are associated with paused polymerase.

Developmental enhancers initiate transcription and are fundamental to our understanding of developmental networks, evolution and disease. Despite their importance, the properties governing enhancer-promoter interactions and their dynamics during embryogenesis remain unclear. At the ß-globin locus, enhancer-promoter interactions appear dynamic and cell-type specific, whereas at the HoxD locus they are stable and ubiquitous, being present in tissues where the target genes are not expressed. The extent to which preformed enhancer-promoter conformations exist at other, more typical, loci and how transcription is eventually triggered is unclear. Here we generated a high-resolution map of enhancer three-dimensional contacts during Drosophila embryogenesis, covering two developmental stages and tissue contexts, at unprecedented resolution. Although local regulatory interactions are common, long-range interactions are highly prevalent within the compact Drosophila genome. Each enhancer contacts multiple enhancers, and promoters with similar expression, suggesting a role in their co-regulation. Notably, most interactions appear unchanged between tissue context and across development, arising before gene activation, and are frequently associated with paused RNA polymerase. Our results indicate that the general topology governing enhancer contacts is conserved from flies to humans and suggest that transcription initiates from preformed enhancer-promoter loops through release of paused polymerase.

PHD and TFIIS-Like domains of the Bye1 transcription factor determine its multivalent genomic distribution.

The BYpass of Ess1 (Bye1) protein is a putative S. cerevisiae transcription factor homologous to the human cancer-associated PHF3/DIDO family of proteins. Bye1 contains a Plant Homeodomain (PHD) and a TFIIS-like domain. The Bye1 PHD finger interacts with tri-methylated lysine 4 of histone H3 (H3K4me3) while the TFIIS-like domain binds to RNA polymerase (Pol) II. Here, we investigated the contribution of these structural features to Bye1 recruitment to chromatin as well as its function in transcriptional regulation. Genome-wide analysis of Bye1 distribution revealed at least two distinct modes of association with actively transcribed genes: within the core of Pol II- and Pol III-transcribed genes concomitant with the presence of the TFIIS transcription factor and, additionally, with promoters of a subset of Pol II-transcribed genes. Specific loss of H3K4me3 abolishes Bye1 association to gene promoters, but doesn't affect its binding within gene bodies. Genetic interactions suggested an essential role of Bye1 in cell fitness under stress conditions compensating the absence of TFIIS. Furthermore, BYE1 deletion resulted in the attenuation of GAL genes expression upon galactose-mediated induction indicating its positive role in transcription regulation. Together, these findings point to a bimodal role of Bye1 in regulation of Pol II transcription. It is recruited via its PHD domain to H3K4 tri-methylated promoters at early steps of transcription. Once Pol II is engaged into elongation, Bye1 binds directly to the transcriptional machinery, modulating its progression along the gene.

Tissue-specific analysis of chromatin state identifies temporal signatures of enhancer activity during embryonic development.

Chromatin modifications are associated with many aspects of gene expression, yet their role in cellular transitions during development remains elusive. Here, we use a new approach to obtain cell type-specific information on chromatin state and RNA polymerase II (Pol II) occupancy within the multicellular Drosophila melanogaster embryo. We directly assessed the relationship between chromatin modifications and the spatio-temporal activity of enhancers. Rather than having a unique chromatin state, active developmental enhancers show heterogeneous histone modifications and Pol II occupancy. Despite this complexity, combined chromatin signatures and Pol II presence are sufficient to predict enhancer activity de novo. Pol II recruitment is highly predictive of the timing of enhancer activity and seems dependent on the timing and location of transcription factor binding. Chromatin modifications typically demarcate large regulatory regions encompassing multiple enhancers, whereas local changes in nucleosome positioning and Pol II occupancy delineate single active enhancers. This cell type-specific view identifies dynamic enhancer usage, an essential step in deciphering developmental networks.

Genomic binding of Pol III transcription machinery and relationship with TFIIS transcription factor distribution in mouse embryonic stem cells.

RNA polymerase (Pol) III synthesizes the tRNAs, the 5S ribosomal RNA and a small number of untranslated RNAs. In vitro, it also transcribes short interspersed nuclear elements (SINEs). We investigated the distribution of Pol III and its associated transcription factors on the genome of mouse embryonic stem cells using a highly specific tandem ChIP-Seq method. Only a subset of the annotated class III genes was bound and thus transcribed. A few hundred SINEs were associated with the Pol III transcription machinery. We observed that Pol III and its transcription factors were present at 30 unannotated sites on the mouse genome, only one of which was conserved in human. An RNA was associated with >80% of these regions. More than 2200 regions bound by TFIIIC transcription factor were devoid of Pol III. These sites were associated with cohesins and often located close to CTCF-binding sites, suggesting that TFIIIC might cooperate with these factors to organize the chromatin. We also investigated the genome-wide distribution of the ubiquitous TFIIS variant, TCEA1. We found that, as in Saccharomyces cerevisiae, TFIIS is associated with class III genes and also with SINEs suggesting that TFIIS is a Pol III transcription factor in mammals.

Analyzing transcription factor occupancy during embryo development using ChIP-seq.

Accurately assessing the binding of transcription factors to cis-regulatory elements in vivo is an essential step toward understanding the mechanisms that govern embryonic development. Genome-wide transcription factor location analysis has been facilitated by the development of high-density tiling arrays (ChIP-on-chip), and more recently by next-generation sequencing technologies, which are used to sequence the DNA fragments obtained from chromatin immunoprecipitation experiments (ChIP-seq). This chapter provides a detailed protocol of the different steps required to generate a successful ChIP-seq library, starting from embryo collection and fixation to chromatin preparation, immunoprecipitation, and finally library preparation. The protocol is optimized for Drosophila embryos, but can be adapted to any organism. The obtained library is suitable for sequencing on an Illumina GAIIx platform.