Binding affinities and cooperative interactions with bHLH activators delimit threshold responses to the dorsal gradient morphogen

Optogenetic control of YAP reveals a dynamic communication code for stem cell fate and proliferation

YAP is a transcriptional regulator that controls pluripotency, cell fate, and proliferation. How cells ensure the selective activation of YAP effector genes is unknown. This knowledge is essential to rationally control cellular decision-making. Here we leverage optogenetics, live-imaging of transcription, and cell fate analysis to understand and control gene activation and cell behavior. We reveal that cells decode the steady-state concentrations and timing of YAP activation to control proliferation, cell fate, and expression of the pluripotency regulators Oct4 and Nanog. While oscillatory YAP inputs induce Oct4 expression and proliferation optimally at frequencies that mimic native dynamics, cellular differentiation requires persistently low YAP levels. We identify the molecular logic of the Oct4 dynamic decoder, which acts through an adaptive change sensor. Our work reveals how YAP levels and dynamics enable multiplexing of information transmission for the regulation of developmental decision-making and establishes a platform for the rational control of these behaviors.

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.

Minimal synthetic enhancers reveal control of the probability of transcriptional engagement and its timing by a morphogen gradient

How enhancers interpret morphogen gradients to generate gene expression patterns is a central question in developmental biology. Recent studies have proposed that enhancers can dictate whether, when, and at what rate promoters engage in transcription, but the complexity of endogenous enhancers calls for theoretical models with too many free parameters to quantitatively dissect these regulatory strategies. To overcome this limitation, we established a minimal promoter-proximal synthetic enhancer in embryos of Drosophila melanogaster. Here, a gradient of the Dorsal activator is read by a single Dorsal DNA binding site. Using live imaging to quantify transcriptional activity, we found that a single binding site can regulate whether promoters engage in transcription in a concentration-dependent manner. By modulating the binding-site affinity, we determined that a gene's decision to transcribe and its transcriptional onset time can be explained by a simple model where the promoter traverses multiple kinetic barriers before transcription can ensue.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

Predictive modeling reveals that higher-order cooperativity drives transcriptional repression in a synthetic developmental enhancer

A challenge in quantitative biology is to predict output patterns of gene expression from knowledge of input transcription factor patterns and from the arrangement of binding sites for these transcription factors on regulatory DNA. We tested whether widespread thermodynamic models could be used to infer parameters describing simple regulatory architectures that inform parameter-free predictions of more complex enhancers in the context of transcriptional repression by Runt in the early fruit fly embryo. By modulating the number and placement of Runt binding sites within an enhancer, and quantifying the resulting transcriptional activity using live imaging, we discovered that thermodynamic models call for higher-order cooperativity between multiple molecular players. This higher-order cooperativity captures the combinatorial complexity underlying eukaryotic transcriptional regulation and cannot be determined from simpler regulatory architectures, highlighting the challenges in reaching a predictive understanding of transcriptional regulation in eukaryotes and calling for approaches that quantitatively dissect their molecular nature.

Systematic analysis of low-affinity transcription factor binding site clusters in vitro and in vivo establishes their functional relevance

Binding to binding site clusters has yet to be characterized in depth, and the functional relevance of low-affinity clusters remains uncertain. We characterized transcription factor binding to low-affinity clusters in vitro and found that transcription factors can bind concurrently to overlapping sites, challenging the notion of binding exclusivity. Furthermore, small clusters with binding sites an order of magnitude lower in affinity give rise to high mean occupancies at physiologically-relevant transcription factor concentrations. To assess whether the observed in vitro occupancies translate to transcriptional activation in vivo, we tested low-affinity binding site clusters in a synthetic and native gene regulatory network in S. cerevisiae. In both systems, clusters of low-affinity binding sites generated transcriptional output comparable to single or even multiple consensus sites. This systematic characterization demonstrates that clusters of low-affinity binding sites achieve substantial occupancies, and that this occupancy can drive expression in eukaryotic promoters.

Concentration-dependent transcriptional switching through a collective action of cis-elements

Gene expression specificity of homeobox transcription factors has remained paradoxical. WUSCHEL activates and represses CLAVATA3 transcription at lower and higher concentrations, respectively. We use computational modeling and experimental analysis to investigate the properties of the cis-regulatory module. We find that intrinsically each cis-element can only activate CLAVATA3 at a higher WUSCHEL concentration. However, together, they repress CLAVATA3 at higher WUSCHEL and activate only at lower WUSCHEL, showing that the concentration-dependent interactions among cis-elements regulate both activation and repression. Biochemical experiments show that two adjacent functional cis-elements bind WUSCHEL with higher affinity and dimerize at relatively lower levels. Moreover, increasing the distance between cis-elements prolongs WUSCHEL monomer binding window, resulting in higher CLAVATA3 activation. Our work showing a constellation of optimally spaced cis-elements of defined affinities determining activation and repression thresholds in regulating CLAVATA3 transcription provides a previously unknown mechanism of cofactor-independent regulation of transcription factor binding in mediating gene expression specificity.

FOXL2 and FOXA1 cooperatively assemble on the TP53 promoter in alternative dimer configurations

Although both the p53 and forkhead box (FOX) family proteins are key transcription factors associated with cancer progression, their direct relationship is unknown. Here, we found that FOX family proteins bind to the non-canonical homotypic cluster of the p53 promoter region (TP53). Analysis of crystal structures of FOX proteins (FOXL2 and FOXA1) bound to the p53 homotypic cluster indicated that they interact with a 2:1 stoichiometry accommodated by FOX-induced DNA allostery. In particular, FOX proteins exhibited distinct dimerization patterns in recognition of the same p53-DNA; dimer formation of FOXA1 involved protein–protein interaction, but FOXL2 did not. Biochemical and biological functional analyses confirmed the cooperative binding of FOX proteins to the TP53 promoter for the transcriptional activation of TP53. In addition, up-regulation of TP53 was necessary for FOX proteins to exhibit anti-proliferative activity in cancer cells. These analyses reveal the presence of a discrete characteristic within FOX family proteins in which FOX proteins regulate the transcription activity of the p53 tumor suppressor via cooperative binding to the TP53 promoter in alternative dimer configurations.

Generation and timing of graded responses to morphogen gradients

Morphogen gradients are known to subdivide a naive cell field into distinct zones of gene expression. Here, we examine whether morphogens can also induce a graded response within such domains. To this end, we explore the role of the Dorsal protein nuclear gradient along the dorsoventral axis in defining the graded pattern of actomyosin constriction that initiates gastrulation in early Drosophila embryos. Two complementary mechanisms for graded accumulation of mRNAs of crucial zygotic Dorsal target genes were identified. First, activation of target-gene expression expands over time from the ventral-most region of high nuclear Dorsal to lateral regions, where the levels are lower, as a result of a Dorsal-dependent activation probability of transcription sites. Thus, sites that are activated earlier will exhibit more mRNA accumulation. Second, once the sites are activated, the rate of RNA Polymerase II loading is also dependent on Dorsal levels. Morphological restrictions require that translation of the graded mRNA be delayed until completion of embryonic cell formation. Such timing is achieved by large introns, which provide a delay in production of the mature mRNAs. Spatio-temporal regulation of key zygotic genes therefore shapes the pattern of gastrulation.

Mechanisms of Binding Specificity among bHLH Transcription Factors

The transcriptome of every cell is orchestrated by the complex network of interaction between transcription factors (TFs) and their binding sites on DNA. Disruption of this network can result in many forms of organism malfunction but also can be the substrate of positive natural selection. However, understanding the specific determinants of each of these individual TF-DNA interactions is a challenging task as it requires integrating the multiple possible mechanisms by which a given TF ends up interacting with a specific genomic region. These mechanisms include DNA motif preferences, which can be determined by nucleotide sequence but also by DNA’s shape; post-translational modifications of the TF, such as phosphorylation; and dimerization partners and co-factors, which can mediate multiple forms of direct or indirect cooperative binding. Binding can also be affected by epigenetic modifications of putative target regions, including DNA methylation and nucleosome occupancy. In this review, we describe how all these mechanisms have a role and crosstalk in one specific family of TFs, the basic helix-loop-helix (bHLH), with a very conserved DNA binding domain and a similar DNA preferred motif, the E-box. Here, we compile and discuss a rich catalog of strategies used by bHLH to acquire TF-specific genome-wide landscapes of binding sites.

Bayesian Markov models improve the prediction of binding motifs beyond first order

Transcription factors (TFs) regulate gene expression by binding to specific DNA motifs. Accurate models for predicting binding affinities are crucial for quantitatively understanding of transcriptional regulation. Motifs are commonly described by position weight matrices, which assume that each position contributes independently to the binding energy. Models that can learn dependencies between positions, for instance, induced by DNA structure preferences, have yielded markedly improved predictions for most TFs on in vivo data. However, they are more prone to overfit the data and to learn patterns merely correlated with rather than directly involved in TF binding. We present an improved, faster version of our Bayesian Markov model software, BaMMmotif2. We tested it with state-of-the-art motif discovery tools on a large collection of ChIP-seq and HT-SELEX datasets. BaMMmotif2 models of fifth-order achieved a median false-discovery-rate-averaged recall 13.6% and 12.2% higher than the next best tool on 427 ChIP-seq datasets and 164 HT-SELEX datasets, respectively, while being 8 to 1000 times faster. BaMMmotif2 models showed no signs of overtraining in cross-cell line and cross-platform tests, with similar improvements on the next-best tool. These results demonstrate that dependencies beyond first order clearly improve binding models for most TFs.

Dynamic patterning by morphogens illuminated by cis-regulatory studies

Morphogen concentration changes in space as well as over time during development. However, how these dynamics are interpreted by cells to specify fate is not well understood. Here, we focus on two morphogens: the maternal transcription factors Bicoid and Dorsal, which directly regulate target genes to pattern Drosophila embryos. The actions of these factors at enhancers has been thoroughly dissected and provides a rich platform for understanding direct input by morphogens and their changing roles over time. Importantly, Bicoid and Dorsal do not work alone; we also discuss additional inputs that work with morphogens to control spatiotemporal gene expression in embryos.

A low affinity cis-regulatory BMP response element restricts target gene activation to subsets of Drosophila neurons

Retrograde BMP signaling and canonical pMad/Medea-mediated transcription regulate diverse target genes across subsets of Drosophila efferent neurons, to differentiate neuropeptidergic neurons and promote motor neuron terminal maturation. How a common BMP signal regulates diverse target genes across many neuronal subsets remains largely unresolved, although available evidence implicates subset-specific transcription factor codes rather than differences in BMP signaling. Here we examine the cis-regulatory mechanisms restricting BMP-induced FMRFa neuropeptide expression to Tv4-neurons. We find that pMad/Medea bind at an atypical, low affinity motif in the FMRFa enhancer. Converting this motif to high affinity caused ectopic enhancer activity and eliminated Tv4-neuron expression. In silico searches identified additional motif instances functional in other efferent neurons, implicating broader functions for this motif in BMP-dependent enhancer activity. Thus, differential interpretation of a common BMP signal, conferred by low affinity pMad/Medea binding motifs, can contribute to the specification of BMP target genes in efferent neuron subsets.

Gli3 utilizes Hand2 to synergistically regulate tissue-specific transcriptional networks

Despite a common understanding that Gli TFs are utilized to convey a Hh morphogen gradient, genetic analyses suggest craniofacial development does not completely fit this paradigm. Using the mouse model (Mus musculus), we demonstrated that rather than being driven by a Hh threshold, robust Gli3 transcriptional activity during skeletal and glossal development required interaction with the basic helix-loop-helix TF Hand2. Not only did genetic and expression data support a co-factorial relationship, but genomic analysis revealed that Gli3 and Hand2 were enriched at regulatory elements for genes essential for mandibular patterning and development. Interestingly, motif analysis at sites co-occupied by Gli3 and Hand2 uncovered mandibular-specific, low-affinity, 'divergent' Gli-binding motifs (dGBMs). Functional validation revealed these dGBMs conveyed synergistic activation of Gli targets essential for mandibular patterning and development. In summary, this work elucidates a novel, sequence-dependent mechanism for Gli transcriptional activity within the craniofacial complex that is independent of a graded Hh signal.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Regulation of spatiotemporal limits of developmental gene expression via enhancer grammar

The regulatory specificity of a gene is determined by the structure of its enhancers, which contain multiple transcription factor binding sites. A unique combination of transcription factor binding sites in an enhancer determines the boundary of target gene expression, and their disruption often leads to developmental defects. Despite extensive characterization of binding motifs in an enhancer, it is still unclear how each binding site contributes to overall transcriptional activity. Using live imaging, quantitative analysis, and mathematical modeling, we measured the contribution of individual binding sites in transcriptional regulation. We show that binding site arrangement within the Rho-GTPase component t48 enhancer mediates the expression boundary by mainly regulating the timing of transcriptional activation along the dorsoventral axis of Drosophila embryos. By tuning the binding affinity of the Dorsal (Dl) and Zelda (Zld) sites, we show that single site modulations are sufficient to induce significant changes in transcription. Yet, no one site seems to have a dominant role; rather, multiple sites synergistically drive increases in transcriptional activity. Interestingly, Dl and Zld demonstrate distinct roles in transcriptional regulation. Dl site modulations change spatial boundaries of t48, mostly by affecting the timing of activation and bursting frequency rather than transcriptional amplitude or bursting duration. However, modulating the binding site for the pioneer factor Zld affects both the timing of activation and amplitude, suggesting that Zld may potentiate higher Dl recruitment to target DNAs. We propose that such fine-tuning of dynamic gene control via enhancer structure may play an important role in ensuring normal development.

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

A strategy for designing allosteric modulators of transcription factor dimerization

Transcription factors in the bHLH family are potentially relevant for tumor growth. Activation requires homodimerization or heterodimerization. Thus, the dimerization step is a likely significant drug target. The oligodendrocyte transcription factor 2 (OLIG2) is overexpressed in gliomas. Here, we developed a fluorescence cross-correlation spectroscopy protocol to examine 10 compounds selected using a pharmacophore-based computational strategy targeting OLIG2 dimerization. We showed that the potency to interact with OLIG2 dimerization in live cells correlates with carcinostatic efficacy. The data indicate a promising approach toward drug development targeting transcription factor overactivity and protein–protein interaction more generally.

Transcription factors (TFs) are fundamental in the regulation of gene expression in the development and differentiation of cells. They may act as oncogenes and when overexpressed in tumors become plausible targets for the design of antitumor agents. Homodimerization or heterodimerization of TFs are required for DNA binding and the association interface between subunits, for the design of allosteric modulators, appears as a privileged structure for the pharmacophore-based computational strategy. Based on this strategy, a set of compounds were earlier identified as potential suppressors of OLIG2 dimerization and found to inhibit tumor growth in a mouse glioblastoma cell line and in a whole-animal study. To investigate whether the antitumor activity is due to the predicted mechanism of action, we undertook a study of OLIG2 dimerization using fluorescence cross-correlation spectroscopy (FCCS) of live HEK cells transfected with 2 spectrally different OLIG2 clones. The selected compounds showed an effect with potency, which correlated with the earlier observed antitumor activity. The OLIG2 proteins showed change in diffusion time under compound treatment in line with dissociation from DNA. The data suggest a general approach of drug discovery based on the design of allosteric modulators of protein–protein interaction.

Multimodal transcriptional control of pattern formation in embryonic development

Predicting how interactions between transcription factors and regulatory DNA sequence dictate rates of transcription and, ultimately, drive developmental outcomes remains an open challenge in physical biology. Using stripe 2 of the even-skipped gene in Drosophila embryos as a case study, we dissect the regulatory forces underpinning a key step along the developmental decision-making cascade: the generation of cytoplasmic mRNA patterns via the control of transcription in individual cells. Using live imaging and computational approaches, we found that the transcriptional burst frequency is modulated across the stripe to control the mRNA production rate. However, we discovered that bursting alone cannot quantitatively recapitulate the formation of the stripe and that control of the window of time over which each nucleus transcribes even-skipped plays a critical role in stripe formation. Theoretical modeling revealed that these regulatory strategies (bursting and the time window) respond in different ways to input transcription factor concentrations, suggesting that the stripe is shaped by the interplay of 2 distinct underlying molecular processes.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

Setting up for gastrulation: D. melanogaster

Drosophila melanogaster embryos develop initially as a syncytium of totipotent nuclei and subsequently, once cellularized, undergo morphogenetic movements associated with gastrulation to generate the three somatic germ layers of the embryo: mesoderm, ectoderm, and endoderm. In this chapter, we focus on the first phase of gastrulation in Drosophila involving patterning of early embryos when cells differentiate their gene expression programs. This patterning process requires coordination of multiple developmental processes including genome reprogramming at the maternal-to-zygotic transition, combinatorial action of transcription factors to support distinct gene expression, and dynamic feedback between this genetic patterning by transcription factors and changes in cell morphology. We discuss the gene regulatory programs acting during patterning to specify the three germ layers, which involve the regulation of spatiotemporal gene expression coupled to physical tissue morphogenesis.

Lighting up the central dogma for predictive developmental biology

Although the last 30years have witnessed the mapping of the wiring diagrams of the gene regulatory networks that dictate cell fate and animal body plans, specific understanding building on such network diagrams that shows how DNA regulatory regions control gene expression lags far behind. These networks have yet to yield the predictive power necessary to, for example, calculate how the concentration dynamics of input transcription factors and DNA regulatory sequence prescribes output patterns of gene expression that, in turn, determine body plans themselves. Here, we argue that reaching a predictive understanding of developmental decision-making calls for an interplay between theory and experiment aimed at revealing how the regulation of the processes of the central dogma dictate network connections and how network topology guides cells toward their ultimate developmental fate. To make this possible, it is crucial to break free from the snapshot-based understanding of embryonic development facilitated by fixed-tissue approaches and embrace new technologies that capture the dynamics of developmental decision-making at the single cell level, in living embryos.

Selective deployment of transcription factor paralogs with submaximal strength facilitates gene regulation in the immune system

In multicellular organisms, duplicated genes can diverge through tissue-specific gene expression patterns, as exemplified by highly regulated expression of Runx transcription factor paralogs with apparent functional redundancy. Here we asked what cell type-specific biologies might be supported by the selective expression of Runx paralogs during Langerhans cell and inducible regulatory T cell differentiation. We uncovered functional non-equivalence between Runx paralogs. Selective expression of native paralogs allowed integration of transcription factor activity with extrinsic signals, while non-native paralogs enforced differentiation even in the absence of exogenous inducers. DNA-binding affinity was controlled by divergent amino acids within the otherwise highly conserved RUNT domain, and evolutionary reconstruction suggested convergence of RUNT domain residues towards sub-maximal strength. Hence, the selective expression of gene duplicates in specialized cell types can synergize with the acquisition of functional differences to enable appropriate gene expression, lineage choice and differentiation in the mammalian immune system.

Mechanism of cognate sequence discrimination by the ETS-family transcription factor ETS-1

Functional evidence increasingly implicates low-affinity DNA recognition by transcription factors as a general mechanism for the spatiotemporal control of developmental genes. Although the DNA sequence requirements for affinity are well-defined, the dynamic mechanisms that execute cognate recognition are much less resolved. To address this gap, here we examined ETS1, a paradigm developmental transcription factor, as a model for which cognate discrimination remains enigmatic. Using molecular dynamics simulations, we interrogated the DNA-binding domain of murine ETS1 alone and when bound to high-and low-affinity cognate sites or to nonspecific DNA. The results of our analyses revealed collective backbone and side-chain motions that distinguished cognate versus nonspecific as well as high- versus low-affinity cognate DNA binding. Combined with binding experiments with site-directed ETS1 mutants, the molecular dynamics data disclosed a triad of residues that respond specifically to low-affinity cognate DNA. We found that a DNA-contacting residue (Gln-336) specifically recognizes low-affinity DNA and triggers the loss of a distal salt bridge (Glu-343/Arg-378) via a large side-chain motion that compromises the hydrophobic packing of two core helices. As an intact Glu-343/Arg-378 bridge is the default state in unbound ETS1 and maintained in high-affinity and nonspecific complexes, the low-affinity complex represents a unique conformational adaptation to the suboptimization of developmental enhancers.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

Insect Transcription Factors: A Landscape of Their Structures and Biological Functions in Drosophila and beyond

Transcription factors (TFs) play essential roles in the transcriptional regulation of functional genes, and are involved in diverse physiological processes in living organisms. The fruit fly Drosophila melanogaster, a simple and easily manipulated organismal model, has been extensively applied to study the biological functions of TFs and their related transcriptional regulation mechanisms. It is noteworthy that with the development of genetic tools such as CRISPR/Cas9 and the next-generation genome sequencing techniques in recent years, identification and dissection the complex genetic regulatory networks of TFs have also made great progress in other insects beyond Drosophila. However, unfortunately, there is no comprehensive review that systematically summarizes the structures and biological functions of TFs in both model and non-model insects. Here, we spend extensive effort in collecting vast related studies, and attempt to provide an impartial overview of the progress of the structure and biological functions of current documented TFs in insects, as well as the classical and emerging research methods for studying their regulatory functions. Consequently, considering the importance of versatile TFs in orchestrating diverse insect physiological processes, this review will assist a growing number of entomologists to interrogate this understudied field, and to propel the progress of their contributions to pest control and even human health.

Hunchback is counter-repressed to regulate even-skipped stripe 2 expression in Drosophila embryos

Hunchback is a bifunctional transcription factor that can activate and repress gene expression in Drosophila development. We investigated the regulatory DNA sequence features that control Hunchback function by perturbing enhancers for one of its target genes, even-skipped (eve). While Hunchback directly represses the eve stripe 3+7 enhancer, we found that in the eve stripe 2+7 enhancer, Hunchback repression is prevented by nearby sequences-this phenomenon is called counter-repression. We also found evidence that Caudal binding sites are responsible for counter-repression, and that this interaction may be a conserved feature of eve stripe 2 enhancers. Our results alter the textbook view of eve stripe 2 regulation wherein Hb is described as a direct activator. Instead, to generate stripe 2, Hunchback repression must be counteracted. We discuss how counter-repression may influence eve stripe 2 regulation and evolution.

Enhancer Logic and Mechanics in Development and Disease

Enhancers are distally located genomic cis-regulatory elements that integrate spatiotemporal cues to coordinate gene expression in a tissue-specific manner during metazoan development. Enhancer function depends on a combination of bound transcription factors and cofactors that regulate local chromatin structure, as well as on the topological interactions that are necessary for their activity. Numerous genome-wide studies concur that the vast majority of disease-associated variations occur within non-coding genomic sequences, in other words the 'cis-regulome', and this underscores their relevance for human health. Advances in DNA sequencing and genome-editing technologies have dramatically expanded our ability to identify enhancers and investigate their properties in vivo, revealing an extraordinary level of interconnectivity underlying cis-regulatory networks. We discuss here these recently developed methodologies, as well as emerging trends and remaining questions in the field of enhancer biology, and how perturbation of enhancer activities/functions results in enhanceropathies.

Capicua controls Toll/IL-1 signaling targets independently of RTK regulation

The HMG-box protein Capicua (Cic) is a conserved transcriptional repressor that functions downstream of receptor tyrosine kinase (RTK) signaling pathways in a relatively simple switch: In the absence of signaling, Cic represses RTK-responsive genes by binding to nearly invariant sites in DNA, whereas activation of RTK signaling down-regulates Cic activity, leading to derepression of its targets. This mechanism controls gene expression in both Drosophila and mammals, but whether Cic can also function via other regulatory mechanisms remains unknown. Here, we characterize an RTK-independent role of Cic in regulating spatially restricted expression of Toll/IL-1 signaling targets in Drosophila embryogenesis. We show that Cic represses those targets by binding to suboptimal DNA sites of lower affinity than its known consensus sites. This binding depends on Dorsal/NF-κB, which translocates into the nucleus upon Toll activation and binds next to the Cic sites. As a result, Cic binds to and represses Toll targets only in regions with nuclear Dorsal. These results reveal a mode of Cic regulation unrelated to the well-established RTK/Cic depression axis and implicate cooperative binding in conjunction with low-affinity binding sites as an important mechanism of enhancer regulation. Given that Cic plays a role in many developmental and pathological processes in mammals, our results raise the possibility that some of these Cic functions are independent of RTK regulation and may depend on cofactor-assisted DNA binding.

Regulatory Principles Governing Tissue Specificity of Developmental Enhancers

Transcriptional enhancers are short segments of genomic DNA (50 bp to 1 kb in length) that can work over long distances (≥1 Mb) to regulate gene expression in specific cells and tissues. Genomic assays have identified on the order of 400,000 to one million putative enhancers in the human genome (e.g., ENCODE Consortium). This suggests that a typical gene is regulated by tens of enhancers, ensuring stringent regulation of gene expression in response to a variety of intrinsic and external signals. Despite the discovery of the first transcriptional enhancer more than 30 years ago, we know surprisingly little about how enhancers regulate gene expression. In particular, the relationship between primary DNA sequence and enhancer specificity remains obscure. Here we summarize recent high-throughput studies in whole embryos aimed at the systematic identification of the sequence and organizational constraints underlying enhancer function and specificity.

Natural variation in stochastic photoreceptor specification and color preference in Drosophila

Each individual perceives the world in a unique way, but little is known about the genetic basis of variation in sensory perception. In the fly eye, the random mosaic of color-detecting R7 photoreceptor subtypes is determined by stochastic on/off expression of the transcription factor Spineless (Ss). In a genome-wide association study, we identified a naturally occurring insertion in a regulatory DNA element in ss that lowers the ratio of Ss^ON to Ss^OFF cells. This change in photoreceptor fates shifts the innate color preference of flies from green to blue. The genetic variant increases the binding affinity for Klumpfuss (Klu), a zinc finger transcriptional repressor that regulates ss expression. Klu is expressed at intermediate levels to determine the normal ratio of Ss^ON to Ss^OFF cells. Thus, binding site affinity and transcription factor levels are finely tuned to regulate stochastic expression, setting the ratio of alternative fates and ultimately determining color preference.

Nuclear microenvironments modulate transcription from low-affinity enhancers

Transcription factors bind low-affinity DNA sequences for only short durations. It is not clear how brief, low-affinity interactions can drive efficient transcription. Here, we report that the transcription factor Ultrabithorax (Ubx) utilizes low-affinity binding sites in the Drosophila melanogaster shavenbaby (svb) locus and related enhancers in nuclear microenvironments of high Ubx concentrations. Related enhancers colocalize to the same microenvironments independently of their chromosomal location, suggesting that microenvironments are highly differentiated transcription domains. Manipulating the affinity of svb enhancers revealed an inverse relationship between enhancer affinity and Ubx concentration required for transcriptional activation. The Ubx cofactor, Homothorax (Hth), was co-enriched with Ubx near enhancers that require Hth, even though Ubx and Hth did not co-localize throughout the nucleus. Thus, microenvironments of high local transcription factor and cofactor concentrations could help low-affinity sites overcome their kinetic inefficiency. Mechanisms that generate these microenvironments could be a general feature of eukaryotic transcriptional regulation.

Uncoupling evolutionary changes in DNA sequence, transcription factor occupancy and enhancer activity

Sequence variation within enhancers plays a major role in both evolution and disease, yet its functional impact on transcription factor (TF) occupancy and enhancer activity remains poorly understood. Here, we assayed the binding of five essential TFs over multiple stages of embryogenesis in two distant Drosophila species (with 1.4 substitutions per neutral site), identifying thousands of orthologous enhancers with conserved or diverged combinatorial occupancy. We used these binding signatures to dissect two properties of developmental enhancers: (1) potential TF cooperativity, using signatures of co-associations and co-divergence in TF occupancy. This revealed conserved combinatorial binding despite sequence divergence, suggesting protein-protein interactions sustain conserved collective occupancy. (2) Enhancer in-vivo activity, revealing orthologous enhancers with conserved activity despite divergence in TF occupancy. Taken together, we identify enhancers with diverged motifs yet conserved occupancy and others with diverged occupancy yet conserved activity, emphasising the need to functionally measure the effect of divergence on enhancer activity.

Systematic dissection of genomic features determining transcription factor binding and enhancer function

Enhancers regulate gene expression through the binding of sequence-specific transcription factors (TFs) to cognate motifs. Various features influence TF binding and enhancer function-including the chromatin state of the genomic locus, the affinities of the binding site, the activity of the bound TFs, and interactions among TFs. However, the precise nature and relative contributions of these features remain unclear. Here, we used massively parallel reporter assays (MPRAs) involving 32,115 natural and synthetic enhancers, together with high-throughput in vivo binding assays, to systematically dissect the contribution of each of these features to the binding and activity of genomic regulatory elements that contain motifs for PPARγ, a TF that serves as a key regulator of adipogenesis. We show that distinct sets of features govern PPARγ binding vs. enhancer activity. PPARγ binding is largely governed by the affinity of the specific motif site and higher-order features of the larger genomic locus, such as chromatin accessibility. In contrast, the enhancer activity of PPARγ binding sites depends on varying contributions from dozens of TFs in the immediate vicinity, including interactions between combinations of these TFs. Different pairs of motifs follow different interaction rules, including subadditive, additive, and superadditive interactions among specific classes of TFs, with both spatially constrained and flexible grammars. Our results provide a paradigm for the systematic characterization of the genomic features underlying regulatory elements, applicable to the design of synthetic regulatory elements or the interpretation of human genetic variation.

Stepwise Progression of Embryonic Patterning

It is long established that the graded distribution of Dorsal transcription factor influences spatial domains of gene expression along the dorsoventral (DV) axis of Drosophila melanogaster embryos. However, the more recent realization that Dorsal levels also change with time raises the question of whether these dynamics are instructive. An overview of DV axis patterning is provided, focusing on new insights identified through quantitative analysis of temporal changes in Dorsal target gene expression from one nuclear cycle to the next ('steps'). Possible roles for the stepwise progression of this patterning program are discussed including (i) tight temporal regulation of signaling pathway activation, (ii) control of gene expression cohorts, and (iii) ensuring the irreversibility of the patterning and cell fate specification process.

Syntax compensates for poor binding sites to encode tissue specificity of developmental enhancers

Transcriptional enhancers are short segments of DNA that switch genes on and off in response to a variety of intrinsic and extrinsic signals. Despite the discovery of the first enhancer more than 30 y ago, the relationship between primary DNA sequence and enhancer activity remains obscure. In particular, the importance of "syntax" (the order, orientation, and spacing of binding sites) is unclear. A high-throughput screen identified synthetic notochord enhancers that are activated by the combination of ZicL and ETS transcription factors in Ciona embryos. Manipulation of these enhancers elucidated a "regulatory code" of sequence and syntax features for notochord-specific expression. This code enabled in silico discovery of bona fide notochord enhancers, including those containing low-affinity binding sites that would be excluded by standard motif identification methods. One of the newly identified enhancers maps upstream of the known enhancer that regulates Brachyury (Ci-Bra), a key determinant of notochord specification. This newly identified Ci-Bra shadow enhancer contains binding sites with very low affinity, but optimal syntax, and therefore mediates surprisingly strong expression in the notochord. Weak binding sites are compensated by optimal syntax, whereas enhancers containing high-affinity binding affinities possess suboptimal syntax. We suggest this balance has obscured the importance of regulatory syntax, as noncanonical binding motifs are typically disregarded by enhancer detection methods. As a result, enhancers with low binding affinities but optimal syntax may be a vastly underappreciated feature of the regulatory genome.

To Be Specific or Not: The Critical Relationship Between Hox And TALE Proteins

Hox proteins are key regulatory transcription factors that act in different tissues of the embryo to provide specific spatial and temporal coordinates to each cell. These patterning functions often depend on the presence of the TALE-homeodomain class cofactors, which form cooperative DNA-binding complexes with all Hox proteins. How this family of cofactors contributes to the highly diverse and specific functions of Hox proteins in vivo remains an important unsolved question. We review here the most recent advances in understanding the molecular mechanisms underlying Hox-TALE function. In particular, we discuss the role of DNA shape, DNA-binding affinity, and protein-protein interaction flexibility in dictating Hox-TALE specificity. We propose several models to explain how these mechanisms are integrated with each other in the context of the many distinct functions that Hox and TALE factors carry out in vivo.

BDNF rs6265 methylation and genotype interact on risk for schizophrenia

Epigenetic mechanisms can mediate gene-environment interactions relevant for complex disorders. The BDNF gene is crucial for development and brain plasticity, is sensitive to environmental stressors, such as hypoxia, and harbors the functional SNP rs6265 (Val(66)Met), which creates or abolishes a CpG dinucleotide for DNA methylation. We found that methylation at the BDNF rs6265 Val allele in peripheral blood of healthy subjects is associated with hypoxia-related early life events (hOCs) and intermediate phenotypes for schizophrenia in a distinctive manner, depending on rs6265 genotype: in ValVal individuals increased methylation is associated with exposure to hOCs and impaired working memory (WM) accuracy, while the opposite is true for ValMet subjects. Also, rs6265 methylation and hOCs interact in modulating WM-related prefrontal activity, another intermediate phenotype for schizophrenia, with an analogous opposite direction in the 2 genotypes. Consistently, rs6265 methylation has a different association with schizophrenia risk in ValVals and ValMets. The relationships of methylation with BDNF levels and of genotype with BHLHB2 binding likely contribute to these opposite effects of methylation. We conclude that BDNF rs6265 methylation interacts with genotype to bridge early environmental exposures to adult phenotypes, relevant for schizophrenia. The study of epigenetic changes in regions containing genetic variation relevant for human diseases may have beneficial implications for the understanding of how genes are actually translated into phenotypes.

The Soft Touch: Low-Affinity Transcription Factor Binding Sites in Development and Evolution

Transcription factor proteins regulate gene expression by binding to specific DNA regions. Most studies of transcription factor binding sites have focused on the highest affinity sites for each factor. There is abundant evidence, however, that binding sites with a range of affinities, including very low affinities, are critical to gene regulation. Here, we present the theoretical and experimental evidence for the importance of low-affinity sites in gene regulation and development. We also discuss the implications of the widespread use of low-affinity sites in eukaryotic genomes for robustness, precision, specificity, and evolution of gene regulation.

Gene Regulatory Mechanisms Underlying the Spatial and Temporal Regulation of Target-Dependent Gene Expression in Drosophila Neurons

Neuronal differentiation often requires target-derived signals from the cells they innervate. These signals typically activate neural subtype-specific genes, but the gene regulatory mechanisms remain largely unknown. Highly restricted expression of the FMRFa neuropeptide in Drosophila Tv4 neurons requires target-derived BMP signaling and a transcription factor code that includes Apterous. Using integrase transgenesis of enhancer reporters, we functionally dissected the Tv4-enhancer of FMRFa within its native cellular context. We identified two essential but discrete cis-elements, a BMP-response element (BMP-RE) that binds BMP-activated pMad, and a homeodomain-response element (HD-RE) that binds Apterous. These cis-elements have low activity and must be combined for Tv4-enhancer activity. Such combinatorial activity is often a mechanism for restricting expression to the intersection of cis-element spatiotemporal activities. However, concatemers of the HD-RE and BMP-RE cis-elements were found to independently generate the same spatiotemporal expression as the Tv4-enhancer. Thus, the Tv4-enhancer atypically combines two low-activity cis-elements that confer the same output from distinct inputs. The activation of target-dependent genes is assumed to 'wait' for target contact. We tested this directly, and unexpectedly found that premature BMP activity could not induce early FMRFa expression; also, we show that the BMP-insensitive HD-RE cis-element is activated at the time of target contact. This led us to uncover a role for the nuclear receptor, seven up (svp), as a repressor of FMRFa induction prior to target contact. Svp is normally downregulated immediately prior to target contact, and we found that maintaining Svp expression prevents cis-element activation, whereas reducing svp gene dosage prematurely activates cis-element activity. We conclude that the target-dependent FMRFa gene is repressed prior to target contact, and that target-derived BMP signaling directly activates FMRFa gene expression through an atypical gene regulatory mechanism.

Morphogen rules: design principles of gradient-mediated embryo patterning

The Drosophila blastoderm and the vertebrate neural tube are archetypal examples of morphogen-patterned tissues that create precise spatial patterns of different cell types. In both tissues, pattern formation is dependent on molecular gradients that emanate from opposite poles. Despite distinct evolutionary origins and differences in time scales, cell biology and molecular players, both tissues exhibit striking similarities in the regulatory systems that establish gene expression patterns that foreshadow the arrangement of cell types. First, signaling gradients establish initial conditions that polarize the tissue, but there is no strict correspondence between specific morphogen thresholds and boundary positions. Second, gradients initiate transcriptional networks that integrate broadly distributed activators and localized repressors to generate patterns of gene expression. Third, the correct positioning of boundaries depends on the temporal and spatial dynamics of the transcriptional networks. These similarities reveal design principles that are likely to be broadly applicable to morphogen-patterned tissues.

Suboptimization of developmental enhancers

Transcriptional enhancers direct precise on-off patterns of gene expression during development. To explore the basis for this precision, we conducted a high-throughput analysis of the Otx-a enhancer, which mediates expression in the neural plate of Ciona embryos in response to fibroblast growth factor (FGF) signaling and a localized GATA determinant. We provide evidence that enhancer specificity depends on submaximal recognition motifs having reduced binding affinities ("suboptimization"). Native GATA and ETS (FGF) binding sites contain imperfect matches to consensus motifs. Perfect matches mediate robust but ectopic patterns of gene expression. The native sites are not arranged at optimal intervals, and subtle changes in their spacing alter enhancer activity. Multiple tiers of enhancer suboptimization produce specific, but weak, patterns of expression, and we suggest that clusters of weak enhancers, including certain "superenhancers," circumvent this trade-off in specificity and activity.

Discovery and information-theoretic characterization of transcription factor binding sites that act cooperatively

Transcription factor binding to the surface of DNA regulatory regions is one of the primary causes of regulating gene expression levels. A probabilistic approach to model protein-DNA interactions at the sequence level is through position weight matrices (PWMs) that estimate the joint probability of a DNA binding site sequence by assuming positional independence within the DNA sequence. Here we construct conditional PWMs that depend on the motif signatures in the flanking DNA sequence, by conditioning known binding site loci on the presence or absence of additional binding sites in the flanking sequence of each site's locus. Pooling known sites with similar flanking sequence patterns allows for the estimation of the conditional distribution function over the binding site sequences. We apply our model to the Dorsal transcription factor binding sites active in patterning the Dorsal-Ventral axis of Drosophila development. We find that those binding sites that cooperate with nearby Twist sites on average contain about 0.5 bits of information about the presence of Twist transcription factor binding sites in the flanking sequence. We also find that Dorsal binding site detectors conditioned on flanking sequence information make better predictions about what is a Dorsal site relative to background DNA than detection without information about flanking sequence features.

The GATA transcription factor GtaC regulates early developmental gene expression dynamics in Dictyostelium

In many systems, including the social amoeba Dictyostelium discoideum, development is often marked by dynamic morphological and transcriptional changes orchestrated by key transcription factors. However, efforts to examine sequential genome-wide changes of gene regulation in developmental processes have been fairly limited. Here we report the developmental regulatory dynamics of GtaC, a GATA-type zinc-finger transcription factor, through the analyses of serial ChIP- and RNA-sequencing data. GtaC is essential for developmental progression, decoding extracellular cAMP pulses during early development and may play a role in mediating cell-type differentiation at later stages. We find that GtaC exhibits temporally distinctive DNA-binding patterns concordant with each developmental stage. We identify direct GtaC targets and observe cotemporaneous GtaC-binding and developmental expression regulation. Our results suggest that GtaC regulates multiple physiological processes as Dictyostelium transitions from a group of unicellular amoebae to an integrated multicellular organism.

Development involves dynamic transcriptional changes. By serial ChIP- and RNA-sequencing, here, the authors show that GtaC, a GATA type transcription factor, exhibits temporally distinctive DNA binding and regulation of gene expression concordant with the development in the social amoeba Dictyostelium discoideum.

Shadow enhancers enable Hunchback bifunctionality in the Drosophila embryo

Hunchback (Hb) is a bifunctional transcription factor that activates and represses distinct enhancers. Here, we investigate the hypothesis that Hb can activate and repress the same enhancer. Computational models predicted that Hb bifunctionally regulates the even-skipped (eve) stripe 3+7 enhancer (eve3+7) in Drosophila blastoderm embryos. We measured and modeled eve expression at cellular resolution under multiple genetic perturbations and found that the eve3+7 enhancer could not explain endogenous eve stripe 7 behavior. Instead, we found that eve stripe 7 is controlled by two enhancers: the canonical eve3+7 and a sequence encompassing the minimal eve stripe 2 enhancer (eve2+7). Hb bifunctionally regulates eve stripe 7, but it executes these two activities on different pieces of regulatory DNA--it activates the eve2+7 enhancer and represses the eve3+7 enhancer. These two "shadow enhancers" use different regulatory logic to create the same pattern.

Su(H)-mediated repression positions gene boundaries along the dorsal-ventral axis of Drosophila embryos

In Drosophila embryos, a nuclear gradient of the Dorsal (Dl) transcription factor directs differential gene expression along the dorsoventral (DV) axis, translating it into distinct domains that specify future mesodermal, neural, and ectodermal territories. However, the mechanisms used to differentially position gene expression boundaries along this axis are not fully understood. Here, using a combination of approaches, including mutant phenotype analyses and chromatin immunoprecipitation, we show that the transcription factor Suppressor of Hairless, Su(H), helps define dorsal boundaries for many genes expressed along the DV axis. Synthetic reporter constructs also provide molecular evidence that Su(H) binding sites support repression and act to counterbalance activation through Dl and the ubiquitous activator Zelda. Our study highlights a role for broadly expressed repressors, like Su(H), and organization of transcription factor binding sites within cis-regulatory modules as important elements controlling spatial domains of gene expression to facilitate flexible positioning of boundaries across the entire DV axis.

Low affinity binding site clusters confer hox specificity and regulatory robustness

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.