Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions

Bone morphogenetic protein signaling: the pathway and its regulation

In the mid-1960s, bone morphogenetic proteins (BMPs) were first identified in the extracts of bone to have the remarkable ability to induce heterotopic bone. When the Drosophila gene decapentaplegic (dpp) was first identified to share sequence similarity with mammalian BMP2/BMP4 in the late-1980s, it became clear that secreted BMP ligands can mediate processes other than bone formation. Following this discovery, collaborative efforts between Drosophila geneticists and mammalian biochemists made use of the strengths of their respective model systems to identify BMP signaling components and delineate the pathway. The ability to conduct genetic modifier screens in Drosophila with relative ease was critical in identifying the intracellular signal transducers for BMP signaling and the related transforming growth factor-beta/activin signaling pathway. Such screens also revealed a host of genes that encode other core signaling components and regulators of the pathway. In this review, we provide a historical account of this exciting time of gene discovery and discuss how the field has advanced over the past 30 years. We have learned that while the core BMP pathway is quite simple, composed of 3 components (ligand, receptor, and signal transducer), behind the versatility of this pathway lies multiple layers of regulation that ensures precise tissue-specific signaling output. We provide a sampling of these discoveries and highlight many questions that remain to be answered to fully understand the complexity of BMP signaling.

Antennapedia: The complexity of a master developmental transcription factor

Hox genes encode transcription factors that play an important role in establishing the basic body plan of animals. In Drosophila, Antennapedia is one of the five genes that make up the Antennapedia complex (ANT-C). Antennapedia determines the identity of the second thoracic segment, known as the mesothorax. Misexpression of Antennapedia at different developmental stages changes the identity of the mesothorax, including the muscles, nervous system, and cuticle. In Drosophila, Antennapedia has two distinct promoters highly regulated throughout development by several transcription factors. Antennapedia proteins are found with other transcription factors in different ANTENNAPEDIA transcriptional complexes to regulate multiple subsets of target genes. In this review, we describe the different mechanisms that regulate the expression and function of Antennapedia and the role of this Hox gene in the development of Drosophila.

Global repression by tailless during segmentation

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes. With ectopic misexpression of tll we detected unexpected repression of almost all pair-rule stripes of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz). Examining Tll embryonic ChIP-chip data with regions mapped as Cis-Regulatory Modules (CRMs) of pair-rule stripes we verified Tll interactions to these regions. With the ChIP-chip data we also verified Tll interactions to the CRMs of gap domains and in the misexpression assay, Tll-mediated repression on Kruppel (Kr), kni (kni) and giant (gt) according to their differential sensitivity to Tll. These results with gap genes confirmed previous data from the literature and argue against indirect repression roles of Tll in the striped pattern. Moreover, the prediction of Tll binding sites in the CRMs of eve stripes and the mathematical modeling of their removal using an experimentally validated theoretical framework shows effects on eve stripes compatible with the absence of a repressor binding to the CRMs. In addition, modeling increased tll levels in the embryo results in the differential repression of eve stripes, agreeing well with the results of the misexpression assay. In genetic assays we investigated eve 5, that is strongly repressed by the ectopic domain and representative of more central stripes not previously implied to be under direct regulation of tll. While this stripe is little affected in tll-, its posterior border is expanded in gt- but detected with even greater expansion in gt-;tll-. We end up by discussing tll with key roles in combinatorial repression mechanisms to contain the expression of medial patterns of the segmentation cascade in the extremities of the embryo.

Poly(ADP-ribosyl)ating enzymes coordinate changes in the expression of metabolic genes with developmental progression

Metabolism, known to be temporally regulated to meet evolving energy demands, plays a crucial role in shaping developmental pace. Recent studies have demonstrated that two key proteins PARP1 and PARG play a regulatory role in the transcription of both morphogenic and metabolic genes. Intriguingly, in Drosophila, the depletion of PARP1 or PARG proteins causes a developmental arrest before pupation, resulting in individuals unable to complete their development. This phenotype highlights the critical involvement of poly(ADP-ribosyl)ating enzymes in regulating the metamorphic process. In this study, we provide compelling evidence that these enzymes intricately coordinate transcriptional changes in both developmental and metabolic pathways during metamorphosis. Specifically, they promote the expression of genes crucial for pupation, while simultaneously negatively regulating the expression of metabolic genes before the transition to the pupal stage. Additionally, these enzymes suppress the expression of genes that are no longer required during this transformative period. Our findings shed light on the intricate interplay between poly(ADP-ribosyl)ating enzymes, developmental processes, and metabolic regulation before metamorphosis and highlight a new role of poly(ADP-ribosyl)ating enzymes in the global regulation of transcription.

Gcm counteracts Toll-induced inflammation and impacts hemocyte number through cholinergic signaling

Hemocytes, the myeloid-like immune cells of Drosophila, fulfill a variety of functions that are not completely understood, ranging from phagocytosis to transduction of inflammatory signals. We here show that downregulating the hemocyte-specific Glial cell deficient/Glial cell missing (Glide/Gcm) transcription factor enhances the inflammatory response to the constitutive activation of the Toll pathway. This correlates with lower levels of glutathione S-transferase, suggesting an implication of Glide/Gcm in reactive oxygen species (ROS) signaling and calling for a widespread anti-inflammatory potential of Glide/Gcm. In addition, our data reveal the expression of acetylcholine receptors in hemocytes and that Toll activation affects their expressions, disclosing a novel aspect of the inflammatory response mediated by neurotransmitters. Finally, we provide evidence for acetylcholine receptor nicotinic acetylcholine receptor alpha 6 (nAchRalpha6) regulating hemocyte proliferation in a cell autonomous fashion and for non-cell autonomous cholinergic signaling regulating the number of hemocytes. Altogether, this study provides new insights on the molecular pathways involved in the inflammatory response.

Identification of transcription factor high accumulation DNA zones

BACKGROUND

Transcription factors (TF) play a crucial role in the regulation of gene transcription; alterations of their activity and binding to DNA areas are strongly involved in cancer and other disease onset and development. For proper biomedical investigation, it is hence essential to correctly trace TF dense DNA areas, having multiple bindings of distinct factors, and select DNA high occupancy target (HOT) zones, showing the highest accumulation of such bindings. Indeed, systematic and replicable analysis of HOT zones in a large variety of cells and tissues would allow further understanding of their characteristics and could clarify their functional role.

RESULTS

Here, we propose, thoroughly explain and discuss a full computational procedure to study in-depth DNA dense areas of transcription factor accumulation and identify HOT zones. This methodology, developed as a computationally efficient parametric algorithm implemented in an R/Bioconductor package, uses a systematic approach with two alternative methods to examine transcription factor bindings and provide comparative and fully-reproducible assessments. It offers different resolutions by introducing three distinct types of accumulation, which can analyze DNA from single-base to region-oriented levels, and a moving window, which can estimate the influence of the neighborhood for each DNA base under exam.

CONCLUSIONS

We quantitatively assessed the full procedure by using our implemented software package, named TFHAZ, in two example applications of biological interest, proving its full reliability and relevance.

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

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

Single-cell dissection of aggression in honeybee colonies

Understanding how genotypic variation results in phenotypic variation is especially difficult for collective behaviour because group phenotypes arise from complex interactions among group members. A genome-wide association study identified hundreds of genes associated with colony-level variation in honeybee aggression, many of which also showed strong signals of positive selection, but the influence of these 'colony aggression genes' on brain function was unknown. Here we use single-cell (sc) transcriptomics and gene regulatory network (GRN) analyses to test the hypothesis that genetic variation for colony aggression influences individual differences in brain gene expression and/or gene regulation. We compared soldiers, which respond to territorial intrusion with stinging attacks, and foragers, which do not. Colony environment showed stronger influences on soldier-forager differences in brain gene regulation compared with brain gene expression. GRN plasticity was strongly associated with colony aggression, with larger differences in GRN dynamics detected between soldiers and foragers from more aggressive relative to less aggressive colonies. The regulatory dynamics of subnetworks composed of genes associated with colony aggression genes were more strongly correlated with each other across different cell types and brain regions relative to other genes, especially in brain regions involved with olfaction and vision and multimodal sensory integration, which are known to mediate bee aggression. These results show how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.

Finding meaning in chaos: a selection signature for functional interactions and its use in molecular biology

A primary goal of biomedical research is to elucidate molecular mechanisms, particularly those responsible for human traits, either normal or pathological. Yet achieving this goal is difficult if not impossible when the traits of interest lack tractable models and so cannot be dissected through time-honoured approaches like forward genetics or reconstitution. Arguably, no biological problem has hindered scientific progress more than this: the inability to dissect a trait's mechanism without a tractable likeness of the trait. At root, forward genetics and reconstitution are powerful approaches because they assay for specific molecular functions. Here, we discuss an alternative way to uncover important mechanistic interactions, namely, to assay for positive natural selection. If an interaction has been selected for, then it must perform an important function, a function that significantly promotes reproductive success. Accordingly, selection is a consequence and indicator of function, and uncovering multimolecular selection will reveal important functional interactions. We propose a selection signature for interactions and review recent selection-based approaches through which to dissect traits that are not inherently tractable. The review includes proof-of-principle studies in which important interactions were uncovered by screening for selection. In sum, screens for selection appear feasible when screens for specific functions are not. Selection screens thus constitute a novel tool through which to reveal the mechanisms that shape the fates of organisms.

A bistable autoregulatory module in the developing embryo commits cells to binary expression fates

Bistable autoactivation has been proposed as a mechanism for cells to adopt binary fates during embryonic development. However, it is unclear whether the autoactivating modules found within developmental gene regulatory networks are bistable, unless their parameters are quantitatively determined. Here, we combine in vivo live imaging with mathematical modeling to dissect the binary cell fate dynamics of the fruit fly pair-rule gene fushi tarazu (ftz), which is regulated by two known enhancers: the early (non-autoregulating) element and the autoregulatory element. Live imaging of transcription and protein concentration in the blastoderm revealed that binary Ftz fates are achieved as Ftz expression rapidly transitions from being dictated by the early element to the autoregulatory element. Moreover, we discovered that Ftz concentration alone is insufficient to activate the autoregulatory element, and that this element only becomes responsive to Ftz at a prescribed developmental time. Based on these observations, we developed a dynamical systems model and quantitated its kinetic parameters directly from experimental measurements. Our model demonstrated that the ftz autoregulatory module is indeed bistable and that the early element transiently establishes the content of the binary cell fate decision to which the autoregulatory module then commits. Further in silico analysis revealed that the autoregulatory element locks the Ftz fate quickly, within 35 min of exposure to the transient signal of the early element. Overall, our work confirms the widely held hypothesis that autoregulation can establish developmental fates through bistability and, most importantly, provides a framework for the quantitative dissection of cellular decision-making.

Precisely timed regulation of enhancer activity defines the binary expression pattern of Fushi tarazu in the Drosophila embryo

The genes that drive development each typically have many different enhancers. Properly coordinating the action of these different enhancers is crucial to correctly specifying cell-fate decisions, yet it remains poorly understood how their activity is choregraphed in time. To shed light on this question, we used recently developed single-cell live imaging tools to dissect the regulation of Fushi tarazu (Ftz) in Drosophila melanogaster embryos. Ftz is a transcription factor that is expressed in asymmetric stripes by two distinct enhancers: autoregulatory and zebra. The anterior edge of each stripe needs to be sharply defined to specify essential lineage boundaries. Here, we tracked how boundary cells commit to either a high-Ftz or low-Ftz fate by measuring Ftz protein traces in real time and simultaneously quantifying transcription from the endogenous locus and individual enhancers. This revealed that the autoregulatory enhancer does not establish this fate choice. Instead, it perpetuates the decision defined by zebra. This is contrary to the prevailing view that autoregulation drives the fate decision by causing bi-stable Ftz expression. Furthermore, we showed that the autoregulatory enhancer is not activated based on a Ftz-concentration threshold but through a timing-based mechanism. We hypothesize that this is regulated by several ubiquitously expressed pioneer-like transcription factors, which have recently been shown to act as timers in the embryo. Our work provides new insight into how precisely timed enhancer activity can directly regulate the dynamics of gene regulatory networks, which may be a general mechanism for ensuring that embryogenesis runs like clockwork.

Transcriptional activators in the early Drosophila embryo perform different kinetic roles

Combinatorial regulation of gene expression by transcription factors (TFs) may in part arise from kinetic synergy-wherein TFs regulate different steps in the transcription cycle. Kinetic synergy requires that TFs play distinguishable kinetic roles. Here, we used live imaging to determine the kinetic roles of three TFs that activate transcription in the Drosophila embryo-Zelda, Bicoid, and Stat92E-by introducing their binding sites into the even-skipped stripe 2 enhancer. These TFs influence different sets of kinetic parameters, and their influence can change over time. All three TFs increased the fraction of transcriptionally active nuclei; Zelda also shortened the first-passage time into transcription and regulated the interval between transcription events. Stat92E also increased the lifetimes of active transcription. Different TFs can therefore play distinct kinetic roles in activating the transcription. This has consequences for understanding the composition and flexibility of regulatory DNA sequences and the biochemical function of TFs. A record of this paper's transparent peer review process is included in the supplemental information.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

An additive repression mechanism sets the anterior limits of anterior pair-rule stripes 1

Segments are repeated anatomical units forming the body of insects. In Drosophila, the specification of the body takes place during the blastoderm through the segmentation cascade. Pair-rule genes such as hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) are of the intermediate level of the cascade and each pair-rule gene is expressed in seven transversal stripes along the antero-posterior axis of the embryo. Stripes are formed by independent cis-regulatory modules (CRMs) under the regulation of transcription factors of maternal source and of gap proteins of the first level of the cascade. The initial blastoderm of Drosophila is a syncytium and it also coincides with the mid-blastula transition when thousands of zygotic genes are transcribed and their products are able to diffuse in the cytoplasm. Thus, we anticipated a complex regulation of the CRMs of the pair-rule stripes. The CRMs of h 1, eve 1, run 1, ftz 1 are able to be activated by bicoid (bcd) throughout the anterior blastoderm and several lines of evidence indicate that they are repressed by the anterior gap genes slp1 (sloppy-paired 1), tll (tailless) and hkb (huckebein). The modest activity of these repressors led to the premise of a combinatorial mechanism regulating the expression of the CRMs of h 1, eve 1, run 1, ftz 1 in more anterior regions of the embryo. We tested this possibility by progressively removing the repression activities of slp1, tll and hkb. In doing so, we were able to expose a mechanism of additive repression limiting the anterior borders of stripes 1. Stripes 1 respond depending on their distance from the anterior end and repressors operating at different levels.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Genome organization controls transcriptional dynamics during development

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

A dual role for DNA binding by Runt in activation and repression of sloppy paired transcription

This work investigates the role of DNA binding by Runt in regulating the sloppy paired 1 (slp1) gene and in particular two distinct cis-regulatory elements that mediate regulation by Runt and other pair-rule transcription factors during Drosophila segmentation. We find that a DNA-binding-defective form of Runt is ineffective at repressing both the distal (DESE) and proximal (PESE) early stripe elements of slp1 and is also compromised for DESE-dependent activation. The function of Runt-binding sites in DESE is further investigated using site-specific transgenesis and quantitative imaging techniques. When DESE is tested as an autonomous enhancer, mutagenesis of the Runt sites results in a clear loss of Runt-dependent repression but has little to no effect on Runt-dependent activation. Notably, mutagenesis of these same sites in the context of a reporter gene construct that also contains the PESE enhancer results in a significant reduction of DESE-dependent activation as well as the loss of repression observed for the autonomous mutant DESE enhancer. These results provide strong evidence that DNA binding by Runt directly contributes to the regulatory interplay of interactions between these two enhancers in the early embryo.

Fitting thermodynamic-based models: Incorporating parameter sensitivity improves the performance of an evolutionary algorithm

A detailed comprehension of transcriptional regulation is critical to understanding the genetic control of development and disease across many different organisms. To more fully investigate the complex molecular interactions controlling the precise expression of genes, many groups have constructed mathematical models to complement their experimental approaches. A critical step in such studies is choosing the most appropriate parameter estimation algorithm to enable detailed analysis of the parameters that contribute to the models. In this study, we develop a novel set of evolutionary algorithms that use a pseudo-random Sobol Set to construct the initial population and incorporate parameter sensitivities into the adaptation of mutation rates, using local, global, and hybrid strategies. Comparison of the performance of these new algorithms to a number of current state-of-the-art global parameter estimation algorithms on a range of continuous test functions, as well as synthetic biological data representing models of gene regulatory systems, reveals improved performance of the new algorithms in terms of runtime, error and reproducibility. In addition, by analyzing the ability of these algorithms to fit datasets of varying quality, we provide the experimentalist with a guide to how the algorithms perform across a range of noisy data. These results demonstrate the improved performance of the new set of parameter estimation algorithms and facilitate meaningful integration of model parameters and predictions in our understanding of the molecular mechanisms of gene regulation.

Extracellular matrix protein N-glycosylation mediates immune self-tolerance in Drosophila melanogaster

In order to respond to infection, hosts must distinguish pathogens from their own tissues. This allows for the precise targeting of immune responses against pathogens and also ensures self-tolerance, the ability of the host to protect self tissues from immune damage. One way to maintain self-tolerance is to evolve a self signal and suppress any immune response directed at tissues that carry this signal. Here, we characterize the Drosophila tuSz1 mutant strain, which mounts an aberrant immune response against its own fat body. We demonstrate that this autoimmunity is the result of two mutations: 1) a mutation in the GCS1 gene that disrupts N-glycosylation of extracellular matrix proteins covering the fat body, and 2) a mutation in the Drosophila Janus Kinase ortholog that causes precocious activation of hemocytes. Our data indicate that N-glycans attached to extracellular matrix proteins serve as a self signal and that activated hemocytes attack tissues lacking this signal. The simplicity of this invertebrate self-recognition system and the ubiquity of its constituent parts suggests it may have functional homologs across animals.

Annotating the Insect Regulatory Genome

An ever-growing number of insect genomes is being sequenced across the evolutionary spectrum. Comprehensive annotation of not only genes but also regulatory regions is critical for reaping the full benefits of this sequencing. Driven by developments in sequencing technologies and in both empirical and computational discovery strategies, the past few decades have witnessed dramatic progress in our ability to identify cis-regulatory modules (CRMs), sequences such as enhancers that play a major role in regulating transcription. Nevertheless, providing a timely and comprehensive regulatory annotation of newly sequenced insect genomes is an ongoing challenge. We review here the methods being used to identify CRMs in both model and non-model insect species, and focus on two tools that we have developed, REDfly and SCRMshaw. These resources can be paired together in a powerful combination to facilitate insect regulatory annotation over a broad range of species, with an accuracy equal to or better than that of other state-of-the-art methods.

A model of active transcription hubs that unifies the roles of active promoters and enhancers

An essential questions of gene regulation is how large number of enhancers and promoters organize into gene regulatory loops. Using transcription-factor binding enrichment as an indicator of enhancer strength, we identified a portion of H3K27ac peaks as potentially strong enhancers and found a universal pattern of promoter and enhancer distribution: At actively transcribed regions of length of ∼200-300 kb, the numbers of active promoters and enhancers are inversely related. Enhancer clusters are associated with isolated active promoters, regardless of the gene's cell-type specificity. As the number of nearby active promoters increases, the number of enhancers decreases. At regions where multiple active genes are closely located, there are few distant enhancers. With Hi-C analysis, we demonstrate that the interactions among the regulatory elements (active promoters and enhancers) occur predominantly in clusters and multiway among linearly close elements and the distance between adjacent elements shows a preference of ∼30 kb. We propose a simple rule of spatial organization of active promoters and enhancers: Gene transcriptions and regulations mainly occur at local active transcription hubs contributed dynamically by multiple elements from linearly close enhancers and/or active promoters. The hub model can be represented with a flower-shaped structure and implies an enhancer-like role of active promoters.

Transcriptional Regulation by (Super)Enhancers: From Discovery to Mechanisms

Accurate control of gene expression in the right cell at the right moment is of fundamental importance to animal development and homeostasis. At the heart of gene regulation lie the enhancers, a class of gene regulatory elements that ensures precise spatiotemporal activation of gene transcription. Mammalian genomes are littered with enhancers, which are frequently organized in cooperative clusters such as locus control regions and superenhancers. Here, we discuss our current knowledge of enhancer biology, including an overview of the discovery of the various enhancer subsets and the mechanistic models used to explain their gene regulatory function.

Inference of dynamic spatial GRN models with multi-GPU evolutionary computation

Reverse engineering mechanistic gene regulatory network (GRN) models with a specific dynamic spatial behavior is an inverse problem without analytical solutions in general. Instead, heuristic machine learning algorithms have been proposed to infer the structure and parameters of a system of equations able to recapitulate a given gene expression pattern. However, these algorithms are computationally intensive as they need to simulate millions of candidate models, which limits their applicability and requires high computational resources. Graphics processing unit (GPU) computing is an affordable alternative for accelerating large-scale scientific computation, yet no method is currently available to exploit GPU technology for the reverse engineering of mechanistic GRNs from spatial phenotypes. Here we present an efficient methodology to parallelize evolutionary algorithms using GPU computing for the inference of mechanistic GRNs that can develop a given gene expression pattern in a multicellular tissue area or cell culture. The proposed approach is based on multi-CPU threads running the lightweight crossover, mutation and selection operators and launching GPU kernels asynchronously. Kernels can run in parallel in a single or multiple GPUs and each kernel simulates and scores the error of a model using the thread parallelism of the GPU. We tested this methodology for the inference of spatiotemporal mechanistic gene regulatory networks (GRNs)-including topology and parameters-that can develop a given 2D gene expression pattern. The results show a 700-fold speedup with respect to a single CPU implementation. This approach can streamline the extraction of knowledge from biological and medical datasets and accelerate the automatic design of GRNs for synthetic biology applications.

Natural variation in the regulation of neurodevelopmental genes modifies flight performance in Drosophila

The winged insects of the order Diptera are colloquially named for their most recognizable phenotype: flight. These insects rely on flight for a number of important life history traits, such as dispersal, foraging, and courtship. Despite the importance of flight, relatively little is known about the genetic architecture of flight performance. Accordingly, we sought to uncover the genetic modifiers of flight using a measure of flies’ reaction and response to an abrupt drop in a vertical flight column. We conducted a genome wide association study (GWAS) using 197 of the Drosophila Genetic Reference Panel (DGRP) lines, and identified a combination of additive and marginal variants, epistatic interactions, whole genes, and enrichment across interaction networks. Egfr, a highly pleiotropic developmental gene, was among the most significant additive variants identified. We functionally validated 13 of the additive candidate genes’ (Adgf-A/Adgf-A2/CG32181, bru1, CadN, flapper (CG11073), CG15236, flippy (CG9766), CREG, Dscam4, form3, fry, Lasp/CG9692, Pde6, Snoo), and introduce a novel approach to whole gene significance screens: PEGASUS_flies. Additionally, we identified ppk23, an Acid Sensing Ion Channel (ASIC) homolog, as an important hub for epistatic interactions. We propose a model that suggests genetic modifiers of wing and muscle morphology, nervous system development and function, BMP signaling, sexually dimorphic neural wiring, and gene regulation are all important for the observed differences flight performance in a natural population. Additionally, these results represent a snapshot of the genetic modifiers affecting drop-response flight performance in Drosophila, with implications for other insects.

Insect flight is a widely recognizable phenotype of many winged insects, hence the name: flies. While fruit flies, or Drosophila melanogaster, are a genetically tractable model, flight performance is a highly integrative phenotype, and therefore challenging to identify comprehensively which genetic modifiers contribute to its genetic architecture. Accordingly, we screened 197 Drosophila Genetic Reference Panel lines for their ability to react and respond to an abrupt drop. Using several computational approaches, we identified additive, marginal, and epistatic variants, as well as whole genes and altered sub-networks of gene-gene and protein-protein interaction networks that contribute to variation in flight performance. More generally, we demonstrate the benefits of employing multiple methodologies to elucidate the genetic architecture of complex traits. Many variants and genes mapped to regions of the genome that affect neurodevelopment, wing and muscle development, and regulation of gene expression. We also introduce PEGASUS_flies, a Drosophila-adapted version of the PEGASUS platform first used in human studies, to infer gene-level significance of association based on the gene’s distribution of individual variant P-values. Our results contribute to the debate over the relative importance of individual, additive factors and epistatic, or higher order, interactions, in the mapping of genotype to phenotype.

AP-1 subunits converge promiscuously at enhancers to potentiate transcription

The AP-1 transcription factor (TF) dimer contributes to many biological processes and environmental responses. AP-1 can be composed of many interchangeable subunits. Unambiguously determining the binding locations of these subunits in the human genome is challenging because of variable antibody specificity and affinity. Here, we definitively establish the genome-wide binding patterns of five AP-1 subunits by using CRISPR to introduce a common antibody tag on each subunit. We find limited evidence for strong dimerization preferences between subunits at steady state and find that, under a stimulus, dimerization patterns reflect changes in the transcriptome. Further, our analysis suggests that canonical AP-1 motifs indiscriminately recruit all AP-1 subunits to genomic sites, which we term AP-1 hotspots. We find that AP-1 hotspots are predictive of cell type–specific gene expression and of genomic responses to glucocorticoid signaling (more so than super-enhancers) and are significantly enriched in disease-associated genetic variants. Together, these results support a model where promiscuous binding of many AP-1 subunits to the same genomic location play a key role in regulating cell type–specific gene expression and environmental responses.

Comparison between instrumental variable and mediation-based methods for reconstructing causal gene networks in yeast

Causal gene networks model the flow of information within a cell. Reconstructing causal networks from omics data is challenging because correlation does not imply causation. When genomics and transcriptomics data from a segregating population are combined, genomic variants can be used to orient the direction of causality between gene expression traits. Instrumental variable methods use a local expression quantitative trait locus (eQTL) as a randomized instrument for a gene's expression level, and assign target genes based on distal eQTL associations. Mediation-based methods additionally require that distal eQTL associations are mediated by the source gene. A detailed comparison between these methods has not yet been conducted, due to the lack of a standardized implementation of different methods, the limited sample size of most multi-omics datasets, and the absence of ground-truth networks for most organisms. Here we used Findr, a software package providing uniform implementations of instrumental variable, mediation, and coexpression-based methods, a recent dataset of 1012 segregants from a cross between two budding yeast strains, and the Yeastract database of known transcriptional interactions to compare causal gene network inference methods. We found that causal inference methods result in a significant overlap with the ground-truth, whereas coexpression did not perform better than random. A subsampling analysis revealed that the performance of mediation saturates at large sample sizes, due to a loss of sensitivity when residual correlations become significant. Instrumental variable methods on the other hand contain false positive predictions, due to genomic linkage between eQTL instruments. Instrumental variable and mediation-based methods also have complementary roles for identifying causal genes underlying transcriptional hotspots. Instrumental variable methods correctly predicted STB5 targets for a hotspot centred on the transcription factor STB5, whereas mediation failed due to Stb5p auto-regulating its own expression. Mediation suggests a new candidate gene, DNM1, for a hotspot on Chr XII, whereas instrumental variable methods could not distinguish between multiple genes located within the hotspot. In conclusion, causal inference from genomics and transcriptomics data is a powerful approach for reconstructing causal gene networks, which could be further improved by the development of methods to control for residual correlations in mediation analyses, and for genomic linkage and pleiotropic effects from transcriptional hotspots in instrumental variable analyses.

Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data

Single-cell RNA-sequencing (scRNAseq) technologies are rapidly evolving. Although very informative, in standard scRNAseq experiments, the spatial organization of the cells in the tissue of origin is lost. Conversely, spatial RNA-seq technologies designed to maintain cell localization have limited throughput and gene coverage. Mapping scRNAseq to genes with spatial information increases coverage while providing spatial location. However, methods to perform such mapping have not yet been benchmarked. To fill this gap, we organized the DREAM Single-Cell Transcriptomics challenge focused on the spatial reconstruction of cells from the Drosophila embryo from scRNAseq data, leveraging as silver standard, genes with in situ hybridization data from the Berkeley Drosophila Transcription Network Project reference atlas. The 34 participating teams used diverse algorithms for gene selection and location prediction, while being able to correctly localize clusters of cells. Selection of predictor genes was essential for this task. Predictor genes showed a relatively high expression entropy, high spatial clustering and included prominent developmental genes such as gap and pair-rule genes and tissue markers. Application of the top 10 methods to a zebra fish embryo dataset yielded similar performance and statistical properties of the selected genes than in the Drosophila data. This suggests that methods developed in this challenge are able to extract generalizable properties of genes that are useful to accurately reconstruct the spatial arrangement of cells in tissues.

Functional Transcription Factor Target Networks Illuminate Control of Epithelial Remodelling

Cell identity is governed by gene expression, regulated by transcription factor (TF) binding at cis-regulatory modules. Decoding the relationship between TF binding patterns and gene regulation is nontrivial, remaining a fundamental limitation in understanding cell decision-making. We developed the NetNC software to predict functionally active regulation of TF targets; demonstrated on nine datasets for the TFs Snail, Twist, and modENCODE Highly Occupied Target (HOT) regions. Snail and Twist are canonical drivers of epithelial to mesenchymal transition (EMT), a cell programme important in development, tumour progression and fibrosis. Predicted "neutral" (non-functional) TF binding always accounted for the majority (50% to 95%) of candidate target genes from statistically significant peaks and HOT regions had higher functional binding than most of the Snail and Twist datasets examined. Our results illuminated conserved gene networks that control epithelial plasticity in development and disease. We identified new gene functions and network modules including crosstalk with notch signalling and regulation of chromatin organisation, evidencing networks that reshape Waddington's epigenetic landscape during epithelial remodelling. Expression of orthologous functional TF targets discriminated breast cancer molecular subtypes and predicted novel tumour biology, with implications for precision medicine. Predicted invasion roles were validated using a tractable cell model, supporting our approach.

Distinct Roles of Broadly Expressed Repressors Support Dynamic Enhancer Action and Change in Time

How broadly expressed repressors regulate gene expression is incompletely understood. To gain insight, we investigated how Suppressor of Hairless-Su(H)-and Runt regulate expression of bone morphogenetic protein (BMP) antagonist short-gastrulation via the sog_Distal enhancer. A live imaging protocol was optimized to capture this enhancer's spatiotemporal output throughout the early Drosophila embryo, finding in this context that Runt regulates transcription initiation, Su(H) regulates transcription rate, and both factors control spatial expression. Furthermore, whereas Su(H) functions as a dedicated repressor, Runt temporally switches from repressor to activator. Our results demonstrate that broad repressors play temporally distinct roles and contribute to dynamic gene expression. Both Run and Su(H)'s ability to influence the spatiotemporal domains of gene expression may serve to counterbalance activators and function in this manner as important regulators of the maternal-to-zygotic transition in early embryos.

Daughterless, the Drosophila orthologue of TCF4, is required for associative learning and maintenance of the synaptic proteome

Mammalian transcription factor 4 (TCF4) has been linked to schizophrenia and intellectual disabilities, such as Pitt-Hopkins syndrome (PTHS). Here, we show that similarly to mammalian TCF4, fruit fly orthologue Daughterless (Da) is expressed widely in the Drosophila brain. Furthermore, silencing of da, using several central nervous system-specific Gal4 driver lines, impairs appetitive associative learning of the larvae and leads to decreased levels of the synaptic proteins Synapsin (Syn) and Discs large 1 (Dlg1), suggesting the involvement of Da in memory formation. Here, we demonstrate that Syn and dlg1 are direct target genes of Da in adult Drosophila heads, as Da binds to the regulatory regions of these genes and the modulation of Da levels alter the levels of Syn and dlg1 mRNA. Silencing of da also affects negative geotaxis of the adult flies, suggesting the impairment of locomotor function. Overall, our findings suggest that Da regulates Drosophila larval memory and adult negative geotaxis, possibly via its synaptic target genes Syn and dlg1 These behavioural phenotypes can be further used as a PTHS model to screen for therapeutics.This article has an associated First Person interview with the first author of the paper.

A Comprehensive Drosophila melanogaster Transcription Factor Interactome

Combinatorial interactions among transcription factors (TFs) play essential roles in generating gene expression specificity and diversity in metazoans. Using yeast 2-hybrid (Y2H) assays on nearly all sequence-specific Drosophila TFs, we identified 1,983 protein-protein interactions (PPIs), more than doubling the number of currently known PPIs among Drosophila TFs. For quality assessment, we validated a subset of our interactions using MITOMI and bimolecular fluorescence complementation assays. We combined our interactome with prior PPI data to generate an integrated Drosophila TF-TF binary interaction network. Our analysis of ChIP-seq data, integrating PPI and gene expression information, uncovered different modes by which interacting TFs are recruited to DNA. We further demonstrate the utility of our Drosophila interactome in shedding light on human TF-TF interactions. This study reveals how TFs interact to bind regulatory elements in vivo and serves as a resource of Drosophila TF-TF binary PPIs for understanding tissue-specific gene regulation.

Combinatorial regulation by transcription factors (TFs) is one mechanism for achieving condition and tissue-specific gene regulation. Shokri et al. mapped TF-TF interactions between most Drosophila TFs, reporting a comprehensive TF-TF network integrated with previously known interactions. They used this network to discern distinct TF-DNA binding modes.

REDfly: the transcriptional regulatory element database for Drosophila

The REDfly database provides a comprehensive curation of experimentally-validated Drosophila transcriptional cis-regulatory elements and includes information on DNA sequence, experimental evidence, patterns of regulated gene expression, and more. Now in its thirteenth year, REDfly has grown to over 23 000 records of tested reporter gene constructs and 2200 tested transcription factor binding sites. Recent developments include the start of curation of predicted cis-regulatory modules in addition to experimentally-verified ones, improved search and filtering, and increased interaction with the authors of curated papers. An expanded data model that will capture information on temporal aspects of gene regulation, regulation in response to environmental and other non-developmental cues, sexually dimorphic gene regulation, and non-endogenous (ectopic) aspects of reporter gene expression is under development and expected to be in place within the coming year. REDfly is freely accessible at http://redfly.ccr.buffalo.edu, and news about database updates and new features can be followed on Twitter at @REDfly_database.

Zygotic pioneer factor activity of Odd-paired/Zic is necessary for late function of the Drosophila segmentation network

Because chromatin determines whether information encoded in DNA is accessible to transcription factors, dynamic chromatin states in development may constrain how gene regulatory networks impart embryonic pattern. To determine the interplay between chromatin states and regulatory network function, we performed ATAC-seq on Drosophila embryos during the establishment of the segmentation network, comparing wild-type and mutant embryos in which all graded maternal patterning inputs are eliminated. While during the period between zygotic genome activation and gastrulation many regions maintain stable accessibility, cis-regulatory modules (CRMs) within the network undergo extensive patterning-dependent changes in accessibility. A component of the network, Odd-paired (opa), is necessary for pioneering accessibility of late segmentation network CRMs. opa-driven changes in accessibility are accompanied by equivalent changes in gene expression. Interfering with the timing of opa activity impacts the proper patterning of expression. These results indicate that dynamic systems for chromatin regulation directly impact the reading of embryonic patterning information.

Eukaryotic transcription factors can track and control their target genes using DNA antennas

Eukaryotic transcription factors (TF) function by binding to short 6-10 bp DNA recognition sites located near their target genes, which are scattered through vast genomes. Such process surmounts enormous specificity, efficiency and celerity challenges using a molecular mechanism that remains poorly understood. Combining biophysical experiments, theory and bioinformatics, we dissect the interplay between the DNA-binding domain of Engrailed, a Drosophila TF, and the regulatory regions of its target genes. We find that Engrailed binding affinity is strongly amplified by the DNA regions flanking the recognition site, which contain long tracts of degenerate recognition-site repeats. Such DNA organization operates as an antenna that attracts TF molecules in a promiscuous exchange among myriads of intermediate affinity binding sites. The antenna ensures a local TF supply, enables gene tracking and fine control of the target site's basal occupancy. This mechanism illuminates puzzling gene expression data and suggests novel engineering strategies to control gene expression.

Two pair-rule responsive enhancers regulate wingless transcription in the Drosophila blastoderm embryo

BACKGROUND

While many developmentally relevant enhancers act in a modular fashion, there is growing evidence for nonadditive interactions between distinct cis-regulatory enhancers. We investigated if nonautonomous enhancer interactions underlie transcription regulation of the Drosophila segment polarity gene, wingless.

RESULTS

We identified two wg enhancers active at the blastoderm stage: wg 3613u, located from -3.6 to -1.3 kb upstream of the wg transcription start site (TSS) and 3046d, located in intron two of the wg gene, from 3.0 to 4.6 kb downstream of the TSS. Genetic experiments confirm that Even Skipped (Eve), Fushi-tarazu (Ftz), Runt, Odd-paired (Opa), Odd-skipped (Odd), and Paired (Prd) contribute to spatially regulated wg expression. Interestingly, there are enhancer specific differences in response to the gain or loss of function of pair-rule gene activity. Although each element recapitulates aspects of wg expression, a composite reporter containing both enhancers more faithfully recapitulates wg regulation than would be predicted from the sum of their individual responses.

CONCLUSION

These results suggest that the regulation of wg by pair-rule genes involves nonadditive interactions between distinct cis-regulatory enhancers.

Exploring a Drosophila Transcription Factor Interaction Network to Identify Cis-Regulatory Modules

Multiple transcription factors (TFs) bind to specific sites in the genome and interact among themselves to form the cis-regulatory modules (CRMs). They are essential in modulating the expression of genes, and it is important to study this interplay to understand gene regulation. In the present study, we integrated experimentally identified TF binding sites collected from published studies with computationally predicted TF binding sites to identify Drosophila CRMs. Along with the detection of the previously known CRMs, this approach identified novel protein combinations. We determined high-occupancy target sites, where a large number of TFs bind. Investigating these sites revealed that Giant, Dichaete, and Knirp are highly enriched in these locations. A common TAG team motif was observed at these sites, which might play a role in recruiting other TFs. While comparing the binding sites at distal and proximal promoters, we found that certain regulatory TFs, such as Zelda, were highly enriched in enhancers. Our study has shown that, from the information available concerning the TF binding sites, the real CRMs could be predicted accurately and efficiently. Although we only may claim co-occurrence of these proteins in this study, it may actually point to their interaction (as known interaction proteins typically co-occur together). Such an integrative approach can, therefore, help us to provide a better understanding of the interplay among the factors, even though further experimental verification is required.

The genome-wide multi-layered architecture of chromosome pairing in early Drosophila embryos

Genome organization involves cis and trans chromosomal interactions, both implicated in gene regulation, development, and disease. Here, we focus on trans interactions in Drosophila, where homologous chromosomes are paired in somatic cells from embryogenesis through adulthood. We first address long-standing questions regarding the structure of embryonic homolog pairing and, to this end, develop a haplotype-resolved Hi-C approach to minimize homolog misassignment and thus robustly distinguish trans-homolog from cis contacts. This computational approach, which we call Ohm, reveals pairing to be surprisingly structured genome-wide, with trans-homolog domains, compartments, and interaction peaks, many coinciding with analogous cis features. We also find a significant genome-wide correlation between pairing, transcription during zygotic genome activation, and binding of the pioneer factor Zelda. Our findings reveal a complex, highly structured organization underlying homolog pairing, first discovered a century ago in Drosophila. Finally, we demonstrate the versatility of our haplotype-resolved approach by applying it to mammalian embryos.

Identification of evolutionarily conserved downstream core promoter elements required for the transcriptional regulation of Fushi tarazu target genes

The regulation of transcription initiation is critical for developmental and cellular processes. RNA polymerase II (Pol II) is recruited by the basal transcription machinery to the core promoter where Pol II initiates transcription. The core promoter encompasses the region from -40 to +40 bp relative to the +1 transcription start site (TSS). Core promoters may contain one or more core promoter motifs that confer specific properties to the core promoter, such as the TATA box, initiator (Inr) and motifs that are located downstream of the TSS, namely, motif 10 element (MTE), the downstream core promoter element (DPE) and the Bridge, a bipartite core promoter element. We had previously shown that Caudal, an enhancer-binding homeodomain transcription factor and a key regulator of the Hox gene network, is a DPE-specific activator. Interestingly, pair-rule proteins have been implicated in enhancer-promoter communication at the engrailed locus. Fushi tarazu (Ftz) is an enhancer-binding homeodomain transcription factor encoded by the ftz pair-rule gene. Ftz works in concert with its co-factor, Ftz-F1, to activate transcription. Here, we examined whether Ftz and Ftz-F1 activate transcription with a preference for a specific core promoter motif. Our analysis revealed that similarly to Caudal, Ftz and Ftz-F1 activate the promoter containing a TATA box mutation to significantly higher levels than the promoter containing a DPE mutation, thus demonstrating a preference for the DPE motif. We further discovered that Ftz target genes are enriched for a combination of functional downstream core promoter elements that are conserved among Drosophila species. Thus, the unique combination (Inr, Bridge and DPE) of functional downstream core promoter elements within Ftz target genes highlights the complexity of transcriptional regulation via the core promoter in the transcription of different developmental gene regulatory networks.

ATAC-seq reveals regional differences in enhancer accessibility during the establishment of spatial coordinates in the Drosophila blastoderm

Establishment of spatial coordinates during Drosophila embryogenesis relies on differential regulatory activity of axis patterning enhancers. Concentration gradients of activator and repressor transcription factors (TFs) provide positional information to each enhancer, which in turn promotes transcription of a target gene in a specific spatial pattern. However, the interplay between an enhancer regulatory activity and its accessibility as determined by local chromatin organization is not well understood. We profiled chromatin accessibility with ATAC-seq in narrow, genetically tagged domains along the antero-posterior axis in the Drosophila blastoderm. We demonstrate that one-quarter of the accessible genome displays significant regional variation in its ATAC-seq signal immediately after zygotic genome activation. Axis patterning enhancers are enriched among the most variable intervals, and their accessibility changes correlate with their regulatory activity. In an embryonic domain where an enhancer receives a net activating TF input and promotes transcription, it displays elevated accessibility in comparison to a domain where it receives a net repressive input. We propose that differential accessibility is a signature of patterning cis-regulatory elements in the Drosophila blastoderm and discuss potential mechanisms by which accessibility of enhancers may be modulated by activator and repressor TFs.

Transcriptional and post-transcriptional regulation of extra macrochaetae during Drosophila adult peripheral neurogenesis

Regulation of the Drosophila ID protein Extra macrochaetae (Emc) is important because reduced Emc levels have been proposed to favor proneural gene activity and thereby define a prepattern for neurogenesis. Recent studies suggest a major role for post-translational control of Emc levels. To further define the mechanisms of Emc regulation, we identified two redundant cis-regulatory regions by germline transformation-rescue experiments that make use of new molecularly-defined emc mutants. We distinguished the mechanisms by which Daughterless (Da) regulated Emc expression, finding post-translational regulation in most tissues, and additional transcriptional regulation in the eye imaginal disc posterior to the morphogenetic furrow. Dpp and Hh signaling pathways repressed Emc transcriptionally and post-translationally within the morphogenetic furrow of the eye disc, whereas Wg signaling repressed Emc expression at the anterior margin of the wing imaginal disc. Although the emc 3' UTR is potentially regulatory, no effect of miRNA pathways on Emc protein levels was discernible. Our work supports recent evidence that post-transcriptional mechanisms contribute more to regulation of Emc protein levels than transcriptional mechanisms do.

Inference of Transcription Factor Regulation Patterns Using Gene Expression Covariation in Natural Populations of Drosophila melanogaster

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Exploiting regulatory heterogeneity to systematically identify enhancers with high accuracy

Identifying functional enhancer elements in metazoan systems is a major challenge. Large-scale validation of enhancers predicted by ENCODE reveal false-positive rates of at least 70%. We used the pregrastrula-patterning network of Drosophila melanogaster to demonstrate that loss in accuracy in held-out data results from heterogeneity of functional signatures in enhancer elements. We show that at least two classes of enhancers are active during early Drosophila embryogenesis and that by focusing on a single, relatively homogeneous class of elements, greater than 98% prediction accuracy can be achieved in a balanced, completely held-out test set. The class of well-predicted elements is composed predominantly of enhancers driving multistage segmentation patterns, which we designate segmentation driving enhancers (SDE). Prediction is driven by the DNA occupancy of early developmental transcription factors, with almost no additional power derived from histone modifications. We further show that improved accuracy is not a property of a particular prediction method: after conditioning on the SDE set, naïve Bayes and logistic regression perform as well as more sophisticated tools. Applying this method to a genome-wide scan, we predict 1,640 SDEs that cover 1.6% of the genome. An analysis of 32 SDEs using whole-mount embryonic imaging of stably integrated reporter constructs chosen throughout our prediction rank-list showed >90% drove expression patterns. We achieved 86.7% precision on a genome-wide scan, with an estimated recall of at least 98%, indicating high accuracy and completeness in annotating this class of functional elements.

Glucocorticoid receptor recruits to enhancers and drives activation by motif-directed binding

Glucocorticoids are potent steroid hormones that regulate immunity and metabolism by activating the transcription factor (TF) activity of glucocorticoid receptor (GR). Previous models have proposed that DNA binding motifs and sites of chromatin accessibility predetermine GR binding and activity. However, there are vast excesses of both features relative to the number of GR binding sites. Thus, these features alone are unlikely to account for the specificity of GR binding and activity. To identify genomic and epigenetic contributions to GR binding specificity and the downstream changes resultant from GR binding, we performed hundreds of genome-wide measurements of TF binding, epigenetic state, and gene expression across a 12-h time course of glucocorticoid exposure. We found that glucocorticoid treatment induces GR to bind to nearly all pre-established enhancers within minutes. However, GR binds to only a small fraction of the set of accessible sites that lack enhancer marks. Once GR is bound to enhancers, a combination of enhancer motif composition and interactions between enhancers then determines the strength and persistence of GR binding, which consequently correlates with dramatic shifts in enhancer activation. Over the course of several hours, highly coordinated changes in TF binding and histone modification occupancy occur specifically within enhancers, and these changes correlate with changes in the expression of nearby genes. Following GR binding, changes in the binding of other TFs precede changes in chromatin accessibility, suggesting that other TFs are also sensitive to genomic features beyond that of accessibility.

Evidence for the temporal regulation of insect segmentation by a conserved sequence of transcription factors

Long-germ insects, such as the fruit fly Drosophila melanogaster, pattern their segments simultaneously, whereas short-germ insects, such as the beetle Tribolium castaneum, pattern their segments sequentially, from anterior to posterior. While the two modes of segmentation at first appear quite distinct, much of this difference might simply reflect developmental heterochrony. We now show here that, in both Drosophila and Tribolium, segment patterning occurs within a common framework of sequential Caudal, Dichaete, and Odd-paired expression. In Drosophila these transcription factors are expressed like simple timers within the blastoderm, while in Tribolium they form wavefronts that sweep from anterior to posterior across the germband. In Drosophila, all three are known to regulate pair-rule gene expression and influence the temporal progression of segmentation. We propose that these regulatory roles are conserved in short-germ embryos, and that therefore the changing expression profiles of these genes across insects provide a mechanistic explanation for observed differences in the timing of segmentation. In support of this hypothesis we demonstrate that Odd-paired is essential for segmentation in Tribolium, contrary to previous reports.

Patterns of chromatin accessibility along the anterior-posterior axis in the early Drosophila embryo

As the Drosophila embryo transitions from the use of maternal RNAs to zygotic transcription, domains of open chromatin, with relatively low nucleosome density and specific histone marks, are established at promoters and enhancers involved in patterned embryonic transcription. However it remains unclear how regions of activity are established during early embryogenesis, and if they are the product of spatially restricted or ubiquitous processes. To shed light on this question, we probed chromatin accessibility across the anterior-posterior axis (A-P) of early Drosophila melanogaster embryos by applying a transposon based assay for chromatin accessibility (ATAC-seq) to anterior and posterior halves of hand-dissected, cellular blastoderm embryos. We find that genome-wide chromatin accessibility is highly similar between the two halves, with regions that manifest significant accessibility in one half of the embryo almost always accessible in the other half, even for promoters that are active in exclusively one half of the embryo. These data support previous studies that show that chromatin accessibility is not a direct result of activity, and point to a role for ubiquitous factors or processes in establishing chromatin accessibility at promoters in the early embryo. However, in concordance with similar works, we find that at enhancers active exclusively in one half of the embryo, we observe a significant skew towards greater accessibility in the region of their activity, highlighting the role of patterning factors such as Bicoid in this process.

VIPER: Visualization Pipeline for RNA-seq, a Snakemake workflow for efficient and complete RNA-seq analysis

Background

RNA sequencing has become a ubiquitous technology used throughout life sciences as an effective method of measuring RNA abundance quantitatively in tissues and cells. The increase in use of RNA-seq technology has led to the continuous development of new tools for every step of analysis from alignment to downstream pathway analysis. However, effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts.

Results

Using the workflow management system Snakemake we have developed a user friendly, fast, efficient, and comprehensive pipeline for RNA-seq analysis. VIPER (Visualization Pipeline for RNA-seq analysis) is an analysis workflow that combines some of the most popular tools to take RNA-seq analysis from raw sequencing data, through alignment and quality control, into downstream differential expression and pathway analysis. VIPER has been created in a modular fashion to allow for the rapid incorporation of new tools to expand the capabilities. This capacity has already been exploited to include very recently developed tools that explore immune infiltrate and T-cell CDR (Complementarity-Determining Regions) reconstruction abilities. The pipeline has been conveniently packaged such that minimal computational skills are required to download and install the dozens of software packages that VIPER uses.

Conclusions

VIPER is a comprehensive solution that performs most standard RNA-seq analyses quickly and effectively with a built-in capacity for customization and expansion.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2139-9) contains supplementary material, which is available to authorized users.

Dynamic genome wide expression profiling of Drosophila head development reveals a novel role of Hunchback in retinal glia cell development and blood-brain barrier integrity

Drosophila melanogaster head development represents a valuable process to study the developmental control of various organs, such as the antennae, the dorsal ocelli and the compound eyes from a common precursor, the eye-antennal imaginal disc. While the gene regulatory network underlying compound eye development has been extensively studied, the key transcription factors regulating the formation of other head structures from the same imaginal disc are largely unknown. We obtained the developmental transcriptome of the eye-antennal discs covering late patterning processes at the late 2^nd larval instar stage to the onset and progression of differentiation at the end of larval development. We revealed the expression profiles of all genes expressed during eye-antennal disc development and we determined temporally co-expressed genes by hierarchical clustering. Since co-expressed genes may be regulated by common transcriptional regulators, we combined our transcriptome dataset with publicly available ChIP-seq data to identify central transcription factors that co-regulate genes during head development. Besides the identification of already known and well-described transcription factors, we show that the transcription factor Hunchback (Hb) regulates a significant number of genes that are expressed during late differentiation stages. We confirm that hb is expressed in two polyploid subperineurial glia cells (carpet cells) and a thorough functional analysis shows that loss of Hb function results in a loss of carpet cells in the eye-antennal disc. Additionally, we provide for the first time functional data indicating that carpet cells are an integral part of the blood-brain barrier. Eventually, we combined our expression data with a de novo Hb motif search to reveal stage specific putative target genes of which we find a significant number indeed expressed in carpet cells.

The development of different cell types must be tightly coordinated, and the eye-antennal imaginal discs of Drosophila melanogaster represent an excellent model to study the molecular mechanisms underlying this coordination. These imaginal discs contain the anlagen of nearly all adult head structures, such as the antennae, the head cuticle, the ocelli and the compound eyes. While large scale screens have been performed to unravel the gene regulatory network underlying compound eye development, a comprehensive understanding of genome wide expression dynamics throughout head development is still missing to date. We studied the genome wide gene expression dynamics during eye-antennal disc development in D. melanogaster to identify new central regulators of the underlying gene regulatory network. Expression based gene clustering and transcription factor motif enrichment analyses revealed a central regulatory role of the transcription factor Hunchback (Hb). We confirmed that hb is expressed in two polyploid retinal subperineurial glia cells (carpet cells). Our functional analysis shows that Hb is necessary for carpet cell development and we show for the first time that the carpet cells are an integral part of the blood-brain barrier.

Iterative random forests to discover predictive and stable high-order interactions

We developed a predictive, stable, and interpretable tool: the iterative random forest algorithm (iRF). iRF discovers high-order interactions among biomolecules with the same order of computational cost as random forests. We demonstrate the efficacy of iRF by finding known and promising interactions among biomolecules, of up to fifth and sixth order, in two data examples in transcriptional regulation and alternative splicing.

Genomics has revolutionized biology, enabling the interrogation of whole transcriptomes, genome-wide binding sites for proteins, and many other molecular processes. However, individual genomic assays measure elements that interact in vivo as components of larger molecular machines. Understanding how these high-order interactions drive gene expression presents a substantial statistical challenge. Building on random forests (RFs) and random intersection trees (RITs) and through extensive, biologically inspired simulations, we developed the iterative random forest algorithm (iRF). iRF trains a feature-weighted ensemble of decision trees to detect stable, high-order interactions with the same order of computational cost as the RF. We demonstrate the utility of iRF for high-order interaction discovery in two prediction problems: enhancer activity in the early Drosophila embryo and alternative splicing of primary transcripts in human-derived cell lines. In Drosophila, among the 20 pairwise transcription factor interactions iRF identifies as stable (returned in more than half of bootstrap replicates), 80% have been previously reported as physical interactions. Moreover, third-order interactions, e.g., between Zelda (Zld), Giant (Gt), and Twist (Twi), suggest high-order relationships that are candidates for follow-up experiments. In human-derived cells, iRF rediscovered a central role of H3K36me3 in chromatin-mediated splicing regulation and identified interesting fifth- and sixth-order interactions, indicative of multivalent nucleosomes with specific roles in splicing regulation. By decoupling the order of interactions from the computational cost of identification, iRF opens additional avenues of inquiry into the molecular mechanisms underlying genome biology.

Discovery of cell-type specific DNA motif grammar in cis-regulatory elements using random Forest

Background

It has been observed that many transcription factors (TFs) can bind to different genomic loci depending on the cell type in which a TF is expressed in, even though the individual TF usually binds to the same core motif in different cell types. How a TF can bind to the genome in such a highly cell-type specific manner, is a critical research question. One hypothesis is that a TF requires co-binding of different TFs in different cell types. If this is the case, it may be possible to observe different combinations of TF motifs – a motif grammar – located at the TF binding sites in different cell types. In this study, we develop a bioinformatics method to systematically identify DNA motifs in TF binding sites across multiple cell types based on published ChIP-seq data, and address two questions: (1) can we build a machine learning classifier to predict cell-type specificity based on motif combinations alone, and (2) can we extract meaningful cell-type specific motif grammars from this classifier model.

Results

We present a Random Forest (RF) based approach to build a multi-class classifier to predict the cell-type specificity of a TF binding site given its motif content. We applied this RF classifier to two published ChIP-seq datasets of TF (TCF7L2 and MAX) across multiple cell types. Using cross-validation, we show that motif combinations alone are indeed predictive of cell types. Furthermore, we present a rule mining approach to extract the most discriminatory rules in the RF classifier, thus allowing us to discover the underlying cell-type specific motif grammar.

Conclusions

Our bioinformatics analysis supports the hypothesis that combinatorial TF motif patterns are cell-type specific.

Electronic supplementary material

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