Ecdysone-Induced 3D Chromatin Reorganization Involves Active Enhancers Bound by Pipsqueak and Polycomb

Cis element length variability does not confer differential transcription factor occupancy at the D. melanogaster histone locus

Histone genes require precise regulation to maintain histone homeostasis and ensure nucleosome stoichiometry. Animal histone genes often have unique clustered genomic organization. However, there is variability of histone gene number and organization as well as differential regulation of the histone genes across species. The Drosophila melanogaster histone locus has unique organizational characteristics as it exists as a series of ∼100 highly regular, tandemly repeated arrays of the 5 replication-dependent histone genes at a single locus. Yet D. melanogaster are viable with only 12 transgenic histone gene arrays. We hypothesized that the histone genes across the locus are differentially regulated. We discovered that the GA-repeat within the H3/H4 promoter is the only variable sequence across the histone gene arrays. The H3/H4 promoter GA-repeat is targeted by CLAMP to promote histone gene expression. We also show two additional GA-binding transcription factors, GAGA Factor and Pipsqueak, target the GA-repeat. When we further examined CLAMP and GAF targeting, we determined that neither CLAMP nor GAF show bias for any GA-repeat lengths. Furthermore, we found that the distribution of GA-repeats targeted by both CLAMP and GAF do not change throughout early development. Together our results suggest that the transcription factors targeting the H3/H4 GA-repeat do not impact differential regulation of the histone genes, but indicate that future studies should interrogate additional cis elements or factors that impact histone gene regulation.

Dan forms condensates in neuroblasts and regulates nuclear architecture and progenitor competence in vivo

Genome organization is thought to underlie cell type specific gene expression, yet how it is regulated in progenitors to produce cellular diversity is unknown. In Drosophila, a developmentally-timed genome reorganization in neural progenitors terminates competence to produce early-born neurons. These events require downregulation of Distal antenna (Dan), part of the conserved pipsqueak DNA-binding superfamily. Here we find that Dan forms liquid-like condensates with high protein mobility, and whose size and subnuclear distribution are balanced with its DNA-binding. Further, we identify a LARKS domain, a structural motif associated with condensate-forming proteins. Deleting just 13 amino acids from LARKS abrogates Dan's ability to retain the early-born neural fate gene, hunchback, in the neuroblast nuclear interior and maintain competence in vivo. Conversely, domain-swapping with LARKS from known phase-separating proteins rescues Dan's effects on competence. Together, we provide in vivo evidence for condensate formation and the regulation of progenitor nuclear architecture underlying neuronal diversification.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

Mechanistic drivers of chromatin organization into compartments

The human genome is not just a simple string of DNA, it is a complex and dynamic entity intricately folded within the cell's nucleus. This three-dimensional organization of chromatin, the combination of DNA and proteins in the nucleus, is crucial for many biological processes and has been prominently studied for its intricate relationship to gene expression. Indeed, the transcriptional machinery does not operate in isolation but interacts intimately with the folded chromatin structure. Techniques for chromatin conformation capture, including genome-wide sequencing approaches, have revealed key organizational features of chromatin, such as the formation of loops by CCCTC-binding factor (CTCF) and the division of loci into chromatin compartments. While much of the recent research and reviews have focused on CTCF loops, we discuss several new revelations that have emerged concerning chromatin compartments, with a particular focus on what is known about mechanistic drivers of compartmentalization. These insights challenge the traditional views of chromatin organization and reveal the complexity behind the formation and maintenance of chromatin compartments.

Drosophila TET acts with PRC1 to activate gene expression independently of its catalytic activity

Enzymes of the ten-eleven translocation (TET) family play a key role in the regulation of gene expression by oxidizing 5-methylcytosine (5mC), a prominent epigenetic mark in many species. Yet, TET proteins also have less characterized noncanonical modes of action, notably in Drosophila, whose genome is devoid of 5mC. Here, we show that Drosophila TET activates the expression of genes required for larval central nervous system (CNS) development mainly in a catalytic-independent manner. Genome-wide profiling shows that TET is recruited to enhancer and promoter regions bound by Polycomb group complex (PcG) proteins. We found that TET interacts and colocalizes on chromatin preferentially with Polycomb repressor complex 1 (PRC1) rather than PRC2. Furthermore, PRC1 but not PRC2 is required for the activation of TET target genes. Last, our results suggest that TET and PRC1 binding to activated genes is interdependent. These data highlight the importance of TET noncatalytic function and the role of PRC1 for gene activation in the Drosophila larval CNS.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

The CLAMP GA-binding transcription factor regulates heat stress-induced transcriptional repression by associating with 3D loop anchors

In order to survive when exposed to heat stress (HS), organisms activate stress response genes and repress constitutive gene expression to prevent the accumulation of potentially toxic RNA and protein products. Although many studies have elucidated the mechanisms that drive HS-induced activation of stress response genes across species, little is known about repression mechanisms or how genes are targeted for activation versus repression context-specifically. The mechanisms of heat stress-regulated activation have been well-studied in Drosophila, in which the GA-binding transcription factor GAF is important for activating genes upon heat stress. Here, we show that a functionally distinct GA-binding transcription factor (TF) protein, CLAMP (Chromatin-linked adaptor for MSL complex proteins), is essential for repressing constitutive genes upon heat stress but not activation of the canonical heat stress pathway. HS induces loss of CLAMP-associated 3D chromatin loop anchors associated with different combinations of GA-binding TFs prior to HS if a gene becomes repressed versus activated. Overall, we demonstrate that CLAMP promotes repression of constitutive genes upon HS, and repression and activation are associated with the loss of CLAMP-associated 3D chromatin loops bound by different combinations of GA-binding TFs.

Sex-specific splicing occurs genome-wide during early Drosophila embryogenesis

Sex-specific splicing is an essential process that regulates sex determination and drives sexual dimorphism. Yet, how early in development widespread sex-specific transcript diversity occurs was unknown because it had yet to be studied at the genome-wide level. We use the powerful Drosophila model to show that widespread sex-specific transcript diversity occurs early in development, concurrent with zygotic genome activation. We also present a new pipeline called time2Splice to quantify changes in alternative splicing over time. Furthermore, we determine that one of the consequences of losing an essential maternally deposited pioneer factor called CLAMP (chromatin-linked adapter for MSL proteins) is altered sex-specific splicing of genes involved in diverse biological processes that drive development. Overall, we show that sex-specific differences in transcript diversity exist even at the earliest stages of development..

Polycomb Recruiters Inside and Outside of the Repressed Domains

The establishment and stable inheritance of individual patterns of gene expression in different cell types are required for the development of multicellular organisms. The important epigenetic regulators are the Polycomb group (PcG) and Trithorax group (TrxG) proteins, which control the silenced and active states of genes, respectively. In Drosophila, the PcG/TrxG group proteins are recruited to the DNA regulatory sequences termed the Polycomb response elements (PREs). The PREs are composed of the binding sites for different DNA-binding proteins, the so-called PcG recruiters. Currently, the role of the PcG recruiters in the targeting of the PcG proteins to PREs is well documented. However, there are examples where the PcG recruiters are also implicated in the active transcription and in the TrxG function. In addition, there is increasing evidence that the genome-wide PcG recruiters interact with the chromatin outside of the PREs and overlap with the proteins of differing regulatory classes. Recent studies of the interactomes of the PcG recruiters significantly expanded our understanding that they have numerous interactors besides the PcG proteins and that their functions extend beyond the regulation of the PRE repressive activity. Here, we summarize current data about the functions of the PcG recruiters.

Identification of ecdysone receptor target genes in the worker honey bee brains during foraging behavior

Ecdysone signaling plays central roles in morphogenesis and female ovarian development in holometabolous insects. In the European honey bee (Apis mellifera L.), however, ecdysone receptor (EcR) is expressed in the brains of adult workers, which have already undergone metamorphosis and are sterile with shrunken ovaries, during foraging behavior. Aiming at unveiling the significance of EcR signaling in the worker brain, we performed chromatin-immunoprecipitation sequencing of EcR to search for its target genes using the brains of nurse bees and foragers. The majority of the EcR targets were common between the nurse bee and forager brains and some of them were known ecdysone signaling-related genes. RNA-sequencing analysis revealed that some EcR target genes were upregulated in forager brains during foraging behavior and some were implicated in the repression of metabolic processes. Single-cell RNA-sequencing analysis revealed that EcR and its target genes were expressed mostly in neurons and partly in glial cells in the optic lobes of the forager brain. These findings suggest that in addition to its role during development, EcR transcriptionally represses metabolic processes during foraging behavior in the adult worker honey bee brain.

GAGA-associated factor fosters loop formation in the Drosophila genome

The impact of genome organization on the control of gene expression persists as a major challenge in regulatory biology. Most efforts have focused on the role of CTCF-enriched boundary elements and TADs, which enable long-range DNA-DNA associations via loop extrusion processes. However, there is increasing evidence for long-range chromatin loops between promoters and distal enhancers formed through specific DNA sequences, including tethering elements, which bind the GAGA-associated factor (GAF). Previous studies showed that GAF possesses amyloid properties in vitro, bridging separate DNA molecules. In this study, we investigated whether GAF functions as a looping factor in Drosophila development. We employed Micro-C assays to examine the impact of defined GAF mutants on genome topology. These studies suggest that the N-terminal POZ/BTB oligomerization domain is important for long-range associations of distant GAGA-rich tethering elements, particularly those responsible for promoter-promoter interactions that coordinate the activities of distant paralogous genes.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

Behavioral Model Deployment for the Transportation Projects within a Smart City Ecosystem: Cases of Germany and South Korea

This research focused on a behavioral model as a significant tangible enabler for smart city plans and initiatives across Asian and EU regions as per transportation projects. This study aimed to create a behavioral model to serve as a planning tool for policymakers, planners, and implementers of transportation initiatives in smart cities. The paper discusses the validity of the proposed model framework for fostering the diffusion of a successful smart city project transformation in a general smart city ecosystem and particularly within the transportation industry. The framework was verified using three different methods: literature review to give a speculative understanding of current smart city approaches; case studies from Germany and South Korea smart city ecosystems that were selected and applied against the behavioral model; and finally, desktop research (behavioral model) performed for smart city project development. As a result, the authors recognized key variables for deriving a possible successful behavioral model as a suggested efficient framework for further smart city strategic projects. Researchers developed and tested, with two validated examples, the suggested behavioral model for smart city projects with a focus on the transportation industry. Results of this survey could help stakeholders in different countries analyze factors influencing decision-making processes and adopt smart city projects to local business environments.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Kinetic principles underlying pioneer function of GAGA transcription factor in live cells

How pioneer factors interface with chromatin to promote accessibility for transcription control is poorly understood in vivo. Here, we directly visualize chromatin association by the prototypical GAGA pioneer factor (GAF) in live Drosophila hemocytes. Single-particle tracking reveals that most GAF is chromatin bound, with a stable-binding fraction showing nucleosome-like confinement residing on chromatin for more than 2 min, far longer than the dynamic range of most transcription factors. These kinetic properties require the full complement of GAF's DNA-binding, multimerization and intrinsically disordered domains, and are autonomous from recruited chromatin remodelers NURF and PBAP, whose activities primarily benefit GAF's neighbors such as Heat Shock Factor. Evaluation of GAF kinetics together with its endogenous abundance indicates that, despite on-off dynamics, GAF constitutively and fully occupies major chromatin targets, thereby providing a temporal mechanism that sustains open chromatin for transcriptional responses to homeostatic, environmental and developmental signals.

Comparative interactome analysis of the PRE DNA-binding factors: purification of the Combgap-, Zeste-, Psq-, and Adf1-associated proteins

The Polycomb group (PcG) and Trithorax group (TrxG) proteins are key epigenetic regulators controlling the silenced and active states of genes in multicellular organisms, respectively. In Drosophila, PcG/TrxG proteins are recruited to the chromatin via binding to specific DNA sequences termed polycomb response elements (PREs). While precise mechanisms of the PcG/TrxG protein recruitment remain unknown, the important role is suggested to belong to sequence-specific DNA-binding factors. At the same time, it was demonstrated that the PRE DNA-binding proteins are not exclusively localized to PREs but can bind other DNA regulatory elements, including enhancers, promoters, and boundaries. To gain an insight into the PRE DNA-binding protein regulatory network, here, using ChIP-seq and immuno-affinity purification coupled to the high-throughput mass spectrometry, we searched for differences in abundance of the Combgap, Zeste, Psq, and Adf1 PRE DNA-binding proteins. While there were no conspicuous differences in co-localization of these proteins with other functional transcription factors, we show that Combgap and Zeste are more tightly associated with the Polycomb repressive complex 1 (PRC1), while Psq interacts strongly with the TrxG proteins, including the BAP SWI/SNF complex. The Adf1 interactome contained Mediator subunits as the top interactors. In addition, Combgap efficiently interacted with AGO2, NELF, and TFIID. Combgap, Psq, and Adf1 have architectural proteins in their networks. We further investigated the existence of direct interactions between different PRE DNA-binding proteins and demonstrated that Combgap-Adf1, Psq-Dsp1, and Pho-Spps can interact in the yeast two-hybrid assay. Overall, our data suggest that Combgap, Psq, Zeste, and Adf1 are associated with the protein complexes implicated in different regulatory activities and indicate their potential multifunctional role in the regulation of transcription.

Pri smORF Peptides Are Wide Mediators of Ecdysone Signaling, Contributing to Shape Spatiotemporal Responses

There is growing evidence that peptides encoded by small open-reading frames (sORF or smORF) can fulfill various cellular functions and define a novel class regulatory molecules. To which extend transcripts encoding only smORF peptides compare with canonical protein-coding genes, yet remain poorly understood. In particular, little is known on whether and how smORF-encoding RNAs might need tightly regulated expression within a given tissue, at a given time during development. We addressed these questions through the analysis of Drosophila polished rice (pri, a.k.a. tarsal less or mille pattes), which encodes four smORF peptides (11-32 amino acids in length) required at several stages of development. Previous work has shown that the expression of pri during epidermal development is regulated in the response to ecdysone, the major steroid hormone in insects. Here, we show that pri transcription is strongly upregulated by ecdysone across a large panel of cell types, suggesting that pri is a core component of ecdysone response. Although pri is produced as an intron-less short transcript (1.5 kb), genetic assays reveal that the developmental functions of pri require an unexpectedly large array of enhancers (spanning over 50 kb), driving a variety of spatiotemporal patterns of pri expression across developing tissues. Furthermore, we found that separate pri enhancers are directly activated by the ecdysone nuclear receptor (EcR) and display distinct regulatory modes between developmental tissues and/or stages. Alike major developmental genes, the expression of pri in a given tissue often involves several enhancers driving apparently redundant (or shadow) expression, while individual pri enhancers can harbor pleiotropic functions across tissues. Taken together, these data reveal the broad role of Pri smORF peptides in ecdysone signaling and show that the cis-regulatory architecture of the pri gene contributes to shape distinct spatial and temporal patterns of ecdysone response throughout development.

CLAMP and Zelda function together to promote Drosophila zygotic genome activation

During the essential and conserved process of zygotic genome activation (ZGA), chromatin accessibility must increase to promote transcription. Drosophila is a well-established model for defining mechanisms that drive ZGA. Zelda (ZLD) is a key pioneer transcription factor (TF) that promotes ZGA in the Drosophila embryo. However, many genomic loci that contain GA-rich motifs become accessible during ZGA independent of ZLD. Therefore, we hypothesized that other early TFs that function with ZLD have not yet been identified, especially those that are capable of binding to GA-rich motifs such as chromatin-linked adaptor for male-specific lethal (MSL) proteins (CLAMP). Here, we demonstrate that Drosophila embryonic development requires maternal CLAMP to (1) activate zygotic transcription; (2) increase chromatin accessibility at promoters of specific genes that often encode other essential TFs; and (3) enhance chromatin accessibility and facilitate ZLD occupancy at a subset of key embryonic promoters. Thus, CLAMP functions as a pioneer factor that plays a targeted yet essential role in ZGA.

GAF is essential for zygotic genome activation and chromatin accessibility in the early Drosophila embryo

Following fertilization, the genomes of the germ cells are reprogrammed to form the totipotent embryo. Pioneer transcription factors are essential for remodeling the chromatin and driving the initial wave of zygotic gene expression. In Drosophila melanogaster, the pioneer factor Zelda is essential for development through this dramatic period of reprogramming, known as the maternal-to-zygotic transition (MZT). However, it was unknown whether additional pioneer factors were required for this transition. We identified an additional maternally encoded factor required for development through the MZT, GAGA Factor (GAF). GAF is necessary to activate widespread zygotic transcription and to remodel the chromatin accessibility landscape. We demonstrated that Zelda preferentially controls expression of the earliest transcribed genes, while genes expressed during widespread activation are predominantly dependent on GAF. Thus, progression through the MZT requires coordination of multiple pioneer-like factors, and we propose that as development proceeds control is gradually transferred from Zelda to GAF.

Proximity-dependent biotin labelling reveals CP190 as an EcR/Usp molecular partner

Proximity-dependent biotin labelling revealed undescribed participants of the ecdysone response in Drosophila. Two labelling enzymes (BioID2 and APEX2) were fused to EcR or Usp to biotin label the surrounding proteins. The EcR/Usp heterodimer was found to collaborate with nuclear pore subunits, chromatin remodelers, and architectural proteins. Many proteins identified through proximity-dependent labelling with EcR/Usp were described previously as functional components of an ecdysone response, corroborating the potency of this labelling method. A link to ecdysone response was confirmed for some newly discovered regulators by immunoprecipitation of prepupal nuclear extract with anti-EcR antibodies and functional experiments in Drosophila S2 cells. A more in-depth study was conducted to clarify the association of EcR/Usp with one of the detected proteins, CP190, a well-described cofactor of Drosophila insulators. CP190 was found to co-immunoprecipitate with the EcR subunit of EcR/Usp in a 20E-independent manner. ChIP-Seq experiments revealed only partial overlapping between CP190 and EcR bound sites in the Drosophila genome and complete absence of CP190 binding at 20E-dependent enhancers. Analysis of Hi-C data demonstrated an existence of remote interactions between 20E-dependent enhancers and CP190 sites which suggests formation of a protein complex between EcR/Usp and CP190 through the space. Our results support the previous concept that CP190 has a role in stabilization of specific chromatin loops for proper activation of transcription of genes regulated by 20E hormone.

Analysis of Hi-C data using SIP effectively identifies loops in organisms from C. elegans to mammals

Chromatin loops are a major component of 3D nuclear organization, visually apparent as intense point-to-point interactions in Hi-C maps. Identification of these loops is a critical part of most Hi-C analyses. However, current methods often miss visually evident CTCF loops in Hi-C data sets from mammals, and they completely fail to identify high intensity loops in other organisms. We present SIP, Significant Interaction Peak caller, and SIPMeta, which are platform independent programs to identify and characterize these loops in a time- and memory-efficient manner. We show that SIP is resistant to noise and sequencing depth, and can be used to detect loops that were previously missed in human cells as well as loops in other organisms. SIPMeta corrects for a common visualization artifact by accounting for Manhattan distance to create average plots of Hi-C and HiChIP data. We then demonstrate that the use of SIP and SIPMeta can lead to biological insights by characterizing the contribution of several transcription factors to CTCF loop stability in human cells. We also annotate loops associated with the SMC component of the dosage compensation complex (DCC) in Caenorhabditis elegans and demonstrate that loop anchors represent bidirectional blocks for symmetrical loop extrusion. This is in contrast to the asymmetrical extrusion until unidirectional blockage by CTCF that is presumed to occur in mammals. Using HiChIP and multiway ligation events, we then show that DCC loops form a network of strong interactions that may contribute to X Chromosome-wide condensation in C. elegans hermaphrodites.