An ancient yet flexible cis-regulatory architecture allows localized Hedgehog tuning by patched/Ptch1

Atlas of nascent RNA transcripts reveals enhancer to gene linkages

Gene transcription is controlled and modulated by regulatory regions, including enhancers and promoters. These regions are abundant in unstable, non-coding bidirectional transcription. Using nascent RNA transcription data across hundreds of human samples, we identified over 800,000 regions containing bidirectional transcription. We then identify highly correlated transcription between bidirectional and gene regions. The identified correlated pairs, a bidirectional region and a gene, are enriched for disease associated SNPs and often supported by independent 3D data. We present these resources as an SQL database which serves as a resource for future studies into gene regulation, enhancer associated RNAs, and transcription factors.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal discs; in contrast, small IDE-reporter transgenes are expressed mainly in the anterior compartment. We found that En binds to the IDEs and suggest that it may directly repress IDE function and modulate En expression levels. We identified two en IDEs, O and S. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from an endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data suggest that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic memory set by the embryonic enhancers. This epigenetic regulation is similar to that of the Ultrabithorax IDEs and we suggest that the activity of late-acting enhancers in other genes may be similarly regulated.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal disks; small IDE-reporter transgenes are expressed in the anterior compartment, the opposite of what is expected. Our data show that the En protein binds to en IDEs, and we suggest that En directly represses IDE function. We identified two en IDEs, 'O' and 'S'. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from the endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data show that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic regulators. Further, the endogenous locus imparts a stability to en function not seen even in a large transgene, reflecting the importance of both positive and negative epigenetic influences that act over relatively large distances in chromatin.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

Behavioral genetics and genomics: Mendel's peas, mice, and bees

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Thermodynamic Modelling of Transcriptional Control: A Sensitivity Analysis

Modelling is a tool used to decipher the biochemical mechanisms involved in transcriptional control. Experimental evidence in genetics is usually supported by theoretical models in order to evaluate the effects of all the possible interactions that can occur in these complicated processes. Models derived from the thermodynamic method are critical in this labour because they are able to take into account multiple mechanisms operating simultaneously at the molecular micro-scale and relate them to transcriptional initiation at the tissular macro-scale. This work is devoted to adapting computational techniques to this context in order to theoretically evaluate the role played by several biochemical mechanisms. The interest of this theoretical analysis relies on the fact that it can be contrasted against those biological experiments where the response to perturbations in the transcriptional machinery environment is evaluated in terms of genetically activated/repressed regions. The theoretical reproduction of these experiments leads to a sensitivity analysis whose results are expressed in terms of the elasticity of a threshold function determining those activated/repressed regions. The study of this elasticity function in thermodynamic models already proposed in the literature reveals that certain modelling approaches can alter the balance between the biochemical mechanisms considered, and this can cause false/misleading outcomes. The reevaluation of classical thermodynamic models gives us a more accurate and complete picture of the interactions involved in gene regulation and transcriptional control, which enables more specific predictions. This sensitivity approach provides a definite advantage in the interpretation of a wide range of genetic experimental results.

Hedgehog Autoprocessing: From Structural Mechanisms to Drug Discovery

Hedgehog (Hh) signaling plays pivotal roles in embryonic development. In adults, Hh signaling is mostly turned off but its abnormal activation is involved in many types of cancer. Hh signaling is initiated by the Hh ligand, generated from the Hh precursor by a specialized autocatalytic process called Hh autoprocessing. The Hh precursor consists of an N-terminal signaling domain (HhN) and a C-terminal autoprocessing domain (HhC). During Hh autoprocessing, the precursor is cleaved between N- and C-terminal domain followed by the covalent ligation of cholesterol to the last residue of HhN, which subsequently leads to the generation of Hh ligand for Hh signaling. Hh autoprocessing is at the origin of canonical Hh signaling and precedes all downstream signaling events. Mutations in the catalytic residues in HhC can lead to congenital defects such as holoprosencephaly (HPE). The aim of this review is to provide an in-depth summary of the progresses and challenges towards an atomic level understanding of the structural mechanisms of Hh autoprocessing. We also discuss drug discovery efforts to inhibit Hh autoprocessing as a new direction in cancer therapy.

GLI transcriptional repression is inert prior to Hedgehog pathway activation

The Hedgehog (HH) pathway regulates a spectrum of developmental processes through the transcriptional mediation of GLI proteins. GLI repressors control tissue patterning by preventing sub-threshold activation of HH target genes, presumably even before HH induction, while lack of GLI repression activates most targets. Despite GLI repression being central to HH regulation, it is unknown when it first becomes established in HH-responsive tissues. Here, we investigate whether GLI3 prevents precocious gene expression during limb development. Contrary to current dogma, we find that GLI3 is inert prior to HH signaling. While GLI3 binds to most targets, loss of Gli3 does not increase target gene expression, enhancer acetylation or accessibility, as it does post-HH signaling. Furthermore, GLI repression is established independently of HH signaling, but after its onset. Collectively, these surprising results challenge current GLI pre-patterning models and demonstrate that GLI repression is not a default state for the HH pathway.

GLI repression has been presumed to be the default transcriptional state and important for pre-patterning tissues. Challenging current models, the authors show that GLI3 repression is inert in the limb bud before the onset of Hedgehog signaling.

Spatial regulation by multiple Gremlin1 enhancers provides digit development with cis-regulatory robustness and evolutionary plasticity

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

The BMP antagonist Gremlin1 balances BMP and SHH signalling, endowing limb bud development with robustness. Here, the authors identify enhancers controlling Grem1 levels in an additive, and spatial regulation in a synergistic manner, providing digit patterning with cis-regulatory robustness and evolutionary plasticity.

Basal epidermis collective migration and local Sonic hedgehog signaling promote skeletal branching morphogenesis in zebrafish fins

Teleost fish fins, like all vertebrate limbs, comprise a series of bones laid out in characteristic pattern. Each fin's distal bony rays typically branch to elaborate skeletal networks providing form and function. Zebrafish caudal fin regeneration studies suggest basal epidermal-expressed Sonic hedgehog (Shh) promotes ray branching by partitioning pools of adjacent pre-osteoblasts. This Shh role is distinct from its well-studied Zone of Polarizing Activity role establishing paired limb positional information. Therefore, we investigated if and how Shh signaling similarly functions during developmental ray branching of both paired and unpaired fins while resolving cellular dynamics of branching by live imaging. We found shha is expressed uniquely by basal epidermal cells overlying pre-osteoblast pools at the distal aspect of outgrowing juvenile fins. Lateral splitting of each shha-expressing epidermal domain followed by the pre-osteoblast pools precedes overt ray branching. We use ptch2:Kaede fish and Kaede photoconversion to identify short stretches of shha+basal epidermis and juxtaposed pre-osteoblasts as the Shh/Smoothened (Smo) active zone. Basal epidermal distal collective movements continuously replenish each shha+domain with individual cells transiently expressing and responding to Shh. In contrast, pre-osteoblasts maintain Shh/Smo activity until differentiating. The Smo inhibitor BMS-833923 prevents branching in all fins, paired and unpaired, with surprisingly minimal effects on caudal fin initial skeletal patterning, ray outgrowth or bone differentiation. Staggered BMS-833923 addition indicates Shh/Smo signaling acts throughout the branching process. We use live cell tracking to find Shh/Smo restrains the distal movement of basal epidermal cells by apparent 'tethering' to pre-osteoblasts. We propose short-range Shh/Smo signaling promotes these heterotypic associations to couple instructive basal epidermal collective movements to pre-osteoblast repositioning as a unique mode of branching morphogenesis.

Enhancer grammar in development, evolution, and disease: dependencies and interplay

Each language has standard books describing that language's grammatical rules. Biologists have searched for similar, albeit more complex, principles relating enhancer sequence to gene expression. Here, we review the literature on enhancer grammar. We introduce dependency grammar, a model where enhancers encode information based on dependencies between enhancer features shaped by mechanistic, evolutionary, and biological constraints. Classifying enhancers based on the types of dependencies may identify unifying principles relating enhancer sequence to gene expression. Such rules would allow us to read the instructions for development within genomes and pinpoint causal enhancer variants underlying disease and evolutionary changes.

Ciliopathic micrognathia is caused by aberrant skeletal differentiation and remodeling

Ciliopathies represent a growing class of diseases caused by defects in microtubule-based organelles called primary cilia. Approximately 30% of ciliopathies are characterized by craniofacial phenotypes such as craniosynostosis, cleft lip/palate and micrognathia. Patients with ciliopathic micrognathia experience a particular set of difficulties, including impaired feeding and breathing, and have extremely limited treatment options. To understand the cellular and molecular basis for ciliopathic micrognathia, we used the talpid² (ta² ), a bona fide avian model for the human ciliopathy oral-facial-digital syndrome subtype 14. Histological analyses revealed that the onset of ciliopathic micrognathia in ta² embryos occurred at the earliest stages of mandibular development. Neural crest-derived skeletal progenitor cells were particularly sensitive to a ciliopathic insult, undergoing unchecked passage through the cell cycle and subsequent increased proliferation. Furthermore, whereas neural crest-derived skeletal differentiation was initiated, osteoblast maturation failed to progress to completion. Additional molecular analyses revealed that an imbalance in the ratio of bone deposition and resorption also contributed to ciliopathic micrognathia in ta² embryos. Thus, our results suggest that ciliopathic micrognathia is a consequence of multiple aberrant cellular processes necessary for skeletal development, and provide potential avenues for future therapeutic treatments.

Robust and efficient gene regulation through localized nuclear microenvironments

Developmental enhancers drive gene expression in specific cell types during animal development. They integrate signals from many different sources mediated through the binding of transcription factors, producing specific responses in gene expression. Transcription factors often bind low-affinity sequences for only short durations. How brief, low-affinity interactions drive efficient transcription and robust gene expression is a central question in developmental biology. Localized high concentrations of transcription factors have been suggested as a possible mechanism by which to use these enhancer sites effectively. Here, we discuss the evidence for such transcriptional microenvironments, mechanisms for their formation and the biological consequences of such sub-nuclear compartmentalization for developmental decisions and evolution.

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.

GLI transcriptional repression regulates tissue-specific enhancer activity in response to Hedgehog signaling

Transcriptional repression needs to be rapidly reversible during embryonic development. This extends to the Hedgehog pathway, which primarily serves to counter GLI repression by processing GLI proteins into transcriptional activators. In investigating the mechanisms underlying GLI repression, we find that a subset of GLI binding regions, termed HH-responsive enhancers, specifically loses acetylation in the absence of HH signaling. These regions are highly enriched around HH target genes and primarily drive HH-specific transcriptional activity in the mouse limb bud. They also retain H3K27ac enrichment in limb buds devoid of GLI activator and repressor, indicating that their activity is primarily regulated by GLI repression. Furthermore, the Polycomb repression complex is not active at most of these regions, suggesting it is not a major mechanism of GLI repression. We propose a model for tissue-specific enhancer activity in which an HDAC-associated GLI repression complex regulates target genes by altering the acetylation status at enhancers.

The formation of the thumb requires direct modulation of Gli3 transcription by Hoxa13

In the tetrapod limb, the digits (fingers or toes) are the elements most subject to morphological diversification in response to functional adaptations. However, despite their functional importance, the mechanisms controlling digit morphology remain poorly understood. Here we have focused on understanding the special morphology of the thumb (digit 1), the acquisition of which was an important adaptation of the human hand. To this end, we have studied the limbs of the Hoxa13 mouse mutant that specifically fail to form digit 1. We show that, consistent with the role of Hoxa13 in Hoxd transcriptional regulation, the expression of Hoxd13 in Hoxa13 mutant limbs does not extend into the presumptive digit 1 territory, which is therefore devoid of distal Hox transcripts, a circumstance that can explain its agenesis. The loss of Hoxd13 expression, exclusively in digit 1 territory, correlates with increased Gli3 repressor activity, a Hoxd negative regulator, resulting from increased Gli3 transcription that, in turn, is due to the release from the negative modulation exerted by Hox13 paralogs on Gli3 regulatory sequences. Our results indicate that Hoxa13 acts hierarchically to initiate the formation of digit 1 by reducing Gli3 transcription and by enabling expansion of the 5'Hoxd second expression phase, thereby establishing anterior-posterior asymmetry in the handplate. Our work uncovers a mutual antagonism between Gli3 and Hox13 paralogs that has important implications for Hox and Gli3 gene regulation in the context of development and evolution.

MicroRNA-129-5p suppresses nasopharyngeal carcinoma lymphangiogenesis and lymph node metastasis by targeting ZIC2

PURPOSE

The etiology of nasopharyngeal carcinoma (NPC) is multifactorial, complex and not fully characterized yet. MicroRNAs (miRNAs or miRs) have been found to contribute to the development and progression of NPC. Here, we aimed to investigate the putative role of miR-129-5p in NPC lymphangiogenesis and lymph node metastasis (LNM), including the involvement of its target gene ZIC2 and the Hedgehog signaling pathway.

METHODS

The expression of miR-129-5p and ZIC2 in primary NPC tissues was assessed using RT-qPCR and Western blot analyses, followed by LNM and lymph vessel density (LVD) correlation analyses. A direct interaction between miR-129-5p and ZIC2 was verified using a dual-luciferase reporter assay. Gain- and loss-of-function experiments were conducted to investigate the effects of miR-129-5p and ZIC2 expression on NPC cell invasion, migration and proliferation in vitro, as well as on LDV and LNM in nude mice in vivo. Additionally, RT-qPCR and Western blot analyses were performed to determine the expression levels of Hedgehog signaling pathway-related factors.

RESULTS

We found that ZIC2 was highly expressed, and miR-129-5p was lowly expressed, in primary NPC tissues. In addition, we found that miR-129-5p can directly bind to and reduce ZIC2 expression. LVD was found to be negatively correlated with miR-129-5p and to be positively correlated with ZIC2 expression. Concomitantly, we found that miR-129-5p abrogated activation of the Hedgehog signaling pathway via ZIC2 targeting, leading to suppression of NPC cell invasion, migration and proliferation in vitro as well as suppression of LNM and LVD in vivo.

CONCLUSIONS

From our data we conclude that miR-129-5p, by decreasing ZIC2 expression, may inhibit NPC lymphangiogenesis and LNM through suppression of the Hedgehog signaling pathway.

Genomic architecture of Shh-dependent cochlear morphogenesis

The mammalian cochlea develops from a ventral outgrowth of the otic vesicle in response to Shh signaling. Mouse embryos lacking Shh or its essential signal transduction components display cochlear agenesis; however, a detailed understanding of the transcriptional network mediating this process is unclear. Here, we describe an integrated genomic approach to identify Shh-dependent genes and associated regulatory sequences that promote cochlear duct morphogenesis. A comparative transcriptome analysis of otic vesicles from mouse mutants exhibiting loss (Smo^ecko ) and gain (Shh-P1) of Shh signaling reveal a set of Shh-responsive genes partitioned into four expression categories in the ventral half of the otic vesicle. This target gene classification scheme provides novel insight into several unanticipated roles for Shh, including priming the cochlear epithelium for subsequent sensory development. We also mapped regions of open chromatin in the inner ear by ATAC-seq that, in combination with Gli2 ChIP-seq, identified inner ear enhancers in the vicinity of Shh-responsive genes. These datasets are useful entry points for deciphering Shh-dependent regulatory mechanisms involved in cochlear duct morphogenesis and establishment of its constituent cell types.

Division of labor in honey bees is associated with transcriptional regulatory plasticity in the brain

Studies in evolutionary and developmental biology show that relationships between transcription factors (TFs) and their target genes can be altered to result in novel regulatory relationships that generate phenotypic plasticity. We hypothesized that context-dependent shifts in the nervous system associated with behavior may also be linked to changes in TF-target relationships over physiological time scales. We tested this hypothesis using honey bee (Apis mellifera) division of labor as a model system by performing bioinformatic analyses of previously published brain transcriptomic profiles together with new RNAi and behavioral experiments. The bioinformatic analyses identified five TFs that exhibited strong signatures of regulatory plasticity as a function of division of labor. RNAi targeting of one of these TFs (broad complex) and a related TF that did not exhibit plasticity (fushi tarazu transcription factor 1) was administered in conjunction with automated analyses of foraging behavior in the field, laboratory assays of aggression and brood care behavior, and endocrine treatments. The results showed that changes in the regulatory relationships of these TFs were associated with behavioral state, social context and endocrine state. These findings provide the first empirical evidence that TF-target relationships in the brain are altered in conjunction with behavior and social context. They also suggest that one mechanism for this plasticity involves pleiotropic TFs high up in regulatory hierarchies producing behavior-specific transcriptional responses by activating different downstream TFs to induce discrete context-dependent transcriptional cascades. These findings provide new insights into the dynamic nature of the transcriptional regulatory architecture underlying behavior in the brain.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

Regulation of GLI1 by cis DNA elements and epigenetic marks

GLI1 is one of three transcription factors (GLI1, GLI2 and GLI3) that mediate the Hedgehog signal transduction pathway and play important roles in normal development. GLI1 and GLI2 form a positive-feedback loop and function as human oncogenes. The mouse and human GLI1 genes have untranslated 5′ exons and large introns 5′ of the translational start. Here we show that Sonic Hedgehog (SHH) stimulates occupancy in the introns by H3K27ac, H3K4me3 and the histone reader protein BRD4. H3K27ac and H3K4me3 occupancy is not significantly changed by removing BRD4 from the human intron and transcription start site (TSS) region. We identified six GLI binding sites (GBS) in the first intron of the human GLI1 gene that are in regions of high sequence conservation among mammals. GLI1 and GLI2 bind all of the GBS in vitro. Elimination of GBS1 and 4 attenuates transcriptional activation by GLI1. Elimination of GBS1, 2, and 4 attenuates transcriptional activation by GLI2. Eliminating all sites essentially eliminates reporter gene activation. Further, GLI1 binds the histone variant H2A.Z. These results suggest that GLI1 and GLI2 can regulate GLI1 expression through protein-protein interactions involving complexes of transcription factors, histone variants, and reader proteins in the regulatory intron of the GLI1 gene. GLI1 acting in trans on the GLI1 intron provides a mechanism for GLI1 positive feedback and auto-regulation. Understanding the combinatorial protein landscape in this locus will be important to interrupting the GLI positive feedback loop and providing new therapeutic approaches to cancers associated with GLI1 overexpression.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Analysis of the transcriptional logic governing differential spatial expression in Hh target genes

This work provides theoretical tools to analyse the transcriptional effects of certain biochemical mechanisms (i.e. affinity and cooperativity) that have been proposed in previous literature to explain the proper spatial expression of Hedgehog target genes involved in Drosophila development. Specifically we have focused on the expression of decapentaplegic, wingless, stripe and patched. The transcription of these genes is believed to be controlled by enhancer modules able to interpret opposing gradients of the activator and repressor forms of the transcription factor Cubitus interruptus (Ci). This study is based on a thermodynamic approach, which provides expression rates for these genes. These expression rates are controlled by transcription factors which are competing and cooperating for common binding sites. We have made mathematical representations of the different expression rates which depend on multiple factors and variables. The expressions obtained with the model have been refined to produce simpler equivalent formulae which allow for their mathematical analysis. Thanks to this, we can evaluate the correlation between the different interactions involved in transcription and the biological features observed at tissular level. These mathematical models can be applied to other morphogenes to help understand the complex transcriptional logic of opposing activator and repressor gradients.

Studying Transcriptional Enhancers: The Founder Fallacy, Validation Creep, and Other Biases

Transcriptional enhancers play a major role in regulating metazoan gene expression. Recent developments in genomics and next-generation sequencing have accelerated and revitalized the study of this important class of sequence elements. Increased interest and attention, however, has also led to troubling trends in the enhancer literature. In this Opinion, I describe some of these issues and show how they arise from shifting and nonuniform enhancer definitions, and genome-era biases. I discuss how they can lead to interpretative errors and an unduly narrow focus on certain aspects of enhancer biology to the potential exclusion of others.

Tuned polymerization of the transcription factor Yan limits off-DNA sequestration to confer context-specific repression

During development, transcriptional complexes at enhancers regulate gene expression in complex spatiotemporal patterns. To achieve robust expression without spurious activation, the affinity and specificity of transcription factor–DNA interactions must be precisely balanced. Protein–protein interactions among transcription factors are also critical, yet how their affinities impact enhancer output is not understood. The Drosophila transcription factor Yan provides a well-suited model to address this, as its function depends on the coordinated activities of two independent and essential domains: the DNA-binding ETS domain and the self-associating SAM domain. To explore how protein–protein affinity influences Yan function, we engineered mutants that increase SAM affinity over four orders of magnitude. This produced a dramatic subcellular redistribution of Yan into punctate structures, reduced repressive output and compromised survival. Cell-type specification and genetic interaction defects suggest distinct requirements for polymerization in different regulatory decisions. We conclude that tuned protein–protein interactions enable the dynamic spectrum of complexes that are required for proper regulation.

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.

Complex cis-regulatory landscape of the insulin receptor gene underlies the broad expression of a central signaling regulator

Insulin signaling plays key roles in development, growth and metabolism through dynamic control of glucose uptake, global protein translation and transcriptional regulation. Altered levels of insulin signaling are known to play key roles in development and disease, yet the molecular basis of such differential signaling remains obscure. Expression of the insulin receptor (InR) gene itself appears to play an important role, but the nature of the molecular wiring controlling InR transcription has not been elucidated. We characterized the regulatory elements driving Drosophila InR expression and found that the generally broad expression of this gene is belied by complex individual switch elements, the dynamic regulation of which reflects direct and indirect contributions of FOXO, EcR, Rbf and additional transcription factors through redundant elements dispersed throughout ∼40 kb of non-coding regions. The control of InR transcription in response to nutritional and tissue-specific inputs represents an integration of multiple cis-regulatory elements, the structure and function of which may have been sculpted by evolutionary selection to provide a highly tailored set of signaling responses on developmental and tissue-specific levels.

Polycomb and Trithorax Group Genes in Drosophila

Polycomb group (PcG) and Trithorax group (TrxG) genes encode important regulators of development and differentiation in metazoans. These two groups of genes were discovered in Drosophila by their opposing effects on homeotic gene (Hox) expression. PcG genes collectively behave as genetic repressors of Hox genes, while the TrxG genes are necessary for HOX gene expression or function. Biochemical studies showed that many PcG proteins are present in two protein complexes, Polycomb repressive complexes 1 and 2, which repress transcription via chromatin modifications. TrxG proteins activate transcription via a variety of mechanisms. Here we summarize the large body of genetic and biochemical experiments in Drosophila on these two important groups of genes.

Galectin-1 induces invasion and the epithelial-mesenchymal transition in human gastric cancer cells via non-canonical activation of the hedgehog signaling pathway

Galectin-1 (Gal-1) has been reported to be an independent prognostic indicator of poor survival in gastric cancer and overexpression of Gal-1 enhances the invasiveness of gastric cancer cells. However, the downstream mechanisms by which Gal-1 promotes invasion remains unclear. Moreover, the function of Gal-1 in the epithelial-mesenchymal transition (EMT) in gastric cancer has not yet been elucidated. In this study, we observed Gal-1 expression was upregulated and positively associated with metastasis and EMT markers in 162 human gastric cancer tissue specimens. In vitro studies showed Gal-1 induced invasion, the EMT phenotype and activated the non-canonical hedgehog (Hh) pathway in gastric cancer cell lines. Furthermore, our data revealed that Gal-1 modulated the non-canonical Hh pathway by increasing the transcription of glioma-associated oncogene-1 (Gli-1) via a Smoothened (SMO)-independent manner, and that upregulation of Gal-1 was strongly associated with gastric cancer metastasis. We conclude that Gal-1 promotes invasion and the EMT in gastric cancer cells via activation of the non-canonical Hh pathway, suggesting Gal-1 could represent a promising therapeutic target for the prevention and treatment of gastric cancer metastasis.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.