Whole-genome ChIP-chip analysis of Dorsal, Twist, and Snail suggests integration of diverse patterning processes in the Drosophila embryo

The Histone Deacetylase Gene Rpd3 Is Required for Starvation Stress Resistance

Epigenetic regulation in starvation is important but not fully understood yet. Here we identified the Rpd3 gene, a Drosophila homolog of histone deacetylase 1, as a critical epigenetic regulator for acquiring starvation stress resistance. Immunostaining analyses of Drosophila fat body revealed that the subcellular localization and levels of Rpd3 dynamically changed responding to starvation stress. In response to starvation stress, the level of Rpd3 rapidly increased, and it accumulated in the nucleolus in what appeared to be foci. These observations suggest that Rpd3 plays a role in regulation of rRNA synthesis in the nucleolus. The RT-qPCR and ChIP-qPCR analyses clarified that Rpd3 binds to the genomic region containing the rRNA promoters and activates rRNA synthesis in response to starvation stress. Polysome analyses revealed that the amount of polysomes was decreased in Rpd3 knockdown flies under starvation stress compared with the control flies. Since the autophagy-related proteins are known to be starvation stress tolerance proteins, we examined autophagy activity, and it was reduced in Rpd3 knockdown flies. Taken together, we conclude that Rpd3 accumulates in the nucleolus in the early stage of starvation, upregulates rRNA synthesis, maintains the polysome amount for translation, and finally increases stress tolerance proteins, such as autophagy-related proteins, to acquire starvation stress resistance.

Drosophila poised enhancers are generated during tissue patterning with the help of repression

Histone modifications are frequently used as markers for enhancer states, but how to interpret enhancer states in the context of embryonic development is not clear. The poised enhancer signature, involving H3K4me1 and low levels of H3K27ac, has been reported to mark inactive enhancers that are poised for future activation. However, future activation is not always observed, and alternative reasons for the widespread occurrence of this enhancer signature have not been investigated. By analyzing enhancers during dorsal-ventral (DV) axis formation in the Drosophila embryo, we find that the poised enhancer signature is specifically generated during patterning in the tissue where the enhancers are not induced, including at enhancers that are known to be repressed by a transcriptional repressor. These results suggest that, rather than serving exclusively as an intermediate step before future activation, the poised enhancer state may be a mark for spatial regulation during tissue patterning. We discuss the possibility that the poised enhancer state is more generally the result of repression by transcriptional repressors.

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.

Transcriptional Control of Developmental Cell Behaviors

Tissue-specific transcription regulators emerged as key developmental control genes, which operate in the context of complex gene regulatory networks (GRNs) to coordinate progressive cell fate specification and tissue morphogenesis. We discuss how GRNs control the individual cell behaviors underlying complex morphogenetic events. Cell behaviors classically range from mesenchymal cell motility to cell shape changes in epithelial sheets. These behaviors emerge from the tissue-specific, multiscale integration of the local activities of universal and pleiotropic effectors, which underlie modular subcellular processes including cytoskeletal dynamics, cell-cell and cell-matrix adhesion, signaling, polarity, and vesicle trafficking. Extrinsic cues and intrinsic cell competence determine the subcellular spatiotemporal patterns of effector activities. GRNs influence most subcellular activities by controlling only a fraction of the effector-coding genes, which we argue is enriched in effectors involved in reading and processing the extrinsic cues to contextualize intrinsic subcellular processes and canalize developmental cell behaviors. The properties of the transcription-cell behavior interface have profound implications for evolution and disease.

Genome-wide identification of Drosophila dorso-ventral enhancers by differential histone acetylation analysis

BACKGROUND

Drosophila dorso-ventral (DV) patterning is one of the best-understood regulatory networks to date, and illustrates the fundamental role of enhancers in controlling patterning, cell fate specification, and morphogenesis during development. Histone acetylation such as H3K27ac is an excellent marker for active enhancers, but it is challenging to obtain precise locations for enhancers as the highest levels of this modification flank the enhancer regions. How to best identify tissue-specific enhancers in a developmental system de novo with a minimal set of data is still unclear.

RESULTS

Using DV patterning as a test system, we develop a simple and effective method to identify tissue-specific enhancers de novo. We sample a broad set of candidate enhancer regions using data on CREB-binding protein co-factor binding or ATAC-seq chromatin accessibility, and then identify those regions with significant differences in histone acetylation between tissues. This method identifies hundreds of novel DV enhancers and outperforms ChIP-seq data of relevant transcription factors when benchmarked with mRNA expression data and transgenic reporter assays. These DV enhancers allow the de novo discovery of the relevant transcription factor motifs involved in DV patterning and contain additional motifs that are evolutionarily conserved and for which the corresponding transcription factors are expressed in a DV-biased fashion. Finally, we identify novel target genes of the regulatory network, implicating morphogenesis genes as early targets of DV patterning.

CONCLUSIONS

Taken together, our approach has expanded our knowledge of the DV patterning network even further and is a general method to identify enhancers in any developmental system, including mammalian development.

Genome-Wide Ultrabithorax Binding Analysis Reveals Highly Targeted Genomic Loci at Developmental Regulators and a Potential Connection to Polycomb-Mediated Regulation

Hox homeodomain transcription factors are key regulators of animal development. They specify the identity of segments along the anterior-posterior body axis in metazoans by controlling the expression of diverse downstream targets, including transcription factors and signaling pathway components. The Drosophila melanogaster Hox factor Ultrabithorax (Ubx) directs the development of thoracic and abdominal segments and appendages, and loss of Ubx function can lead for example to the transformation of third thoracic segment appendages (e.g. halters) into second thoracic segment appendages (e.g. wings), resulting in a characteristic four-wing phenotype. Here we present a Drosophila melanogaster strain with a V5-epitope tagged Ubx allele, which we employed to obtain a high quality genome-wide map of Ubx binding sites using ChIP-seq. We confirm the sensitivity of the V5 ChIP-seq by recovering 7/8 of well-studied Ubx-dependent cis-regulatory regions. Moreover, we show that Ubx binding is predictive of enhancer activity as suggested by comparison with a genome-scale resource of in vivo tested enhancer candidates. We observed densely clustered Ubx binding sites at 12 extended genomic loci that included ANTP-C, BX-C, Polycomb complex genes, and other regulators and the clustered binding sites were frequently active enhancers. Furthermore, Ubx binding was detected at known Polycomb response elements (PREs) and was associated with significant enrichments of Pc and Pho ChIP signals in contrast to binding sites of other developmental TFs. Together, our results show that Ubx targets developmental regulators via strongly clustered binding sites and allow us to hypothesize that regulation by Ubx might involve Polycomb group proteins to maintain specific regulatory states in cooperative or mutually exclusive fashion, an attractive model that combines two groups of proteins with prominent gene regulatory roles during animal development.

Global analysis of dorsoventral patterning in the wasp Nasonia reveals extensive incorporation of novelty in a regulatory network

Background

Gene regulatory networks (GRNs) underlie developmental patterning and morphogenetic processes, and changes in the interactions within the underlying GRNs are a major driver of evolutionary processes. In order to make meaningful comparisons that can provide significant insights into the evolution of regulatory networks, homologous networks from multiple taxa must be deeply characterized. One of the most thoroughly characterized GRNs is the dorsoventral (DV) patterning system of the Drosophila melanogaster embryo. We have developed the wasp Nasonia as a comparative DV patterning model because it has shown the convergent evolution of a mode of early embryonic patterning very similar to that of the fly, and it is of interest to know whether the similarity at the gross level also extends to the molecular level.

Results

We used RNAi to dorsalize and ventralize Nasonia embryos, RNAseq to quantify transcriptome-wide expression levels, and differential expression analysis to identify genes whose expression levels change in either RNAi case. This led to the identification of >100 genes differentially expressed and regulated along the DV axis. Only a handful of these genes are shared DV components in both fly and wasp. Many of those unique to Nasonia are cytoskeletal and adhesion molecules, which may be related to the divergent cell and tissue behavior observed at gastrulation. In addition, many transcription factors and signaling components are only DV regulated in Nasonia, likely reflecting the divergent upstream patterning mechanisms involved in producing the conserved pattern of cell fates observed at gastrulation. Finally, several genes that lack Drosophila orthologs show robust and distinct expression patterns. These include genes with vertebrate homologs that have been lost in the fly lineage, genes that are found only among Hymenoptera, and several genes that entered the Nasonia genome through lateral transfer from endosymbiotic bacteria.

Conclusions

Altogether, our results provide insights into how GRNs respond to new functional demands and how they can incorporate novel components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0285-y) contains supplementary material, which is available to authorized users.

Regulation of the BMP Signaling-Responsive Transcriptional Network in the Drosophila Embryo

The BMP signaling pathway has a conserved role in dorsal-ventral axis patterning during embryonic development. In Drosophila, graded BMP signaling is transduced by the Mad transcription factor and opposed by the Brinker repressor. In this study, using the Drosophila embryo as a model, we combine RNA-seq with Mad and Brinker ChIP-seq to decipher the BMP-responsive transcriptional network underpinning differentiation of the dorsal ectoderm during dorsal-ventral axis patterning. We identify multiple new BMP target genes, including positive and negative regulators of EGF signaling. Manipulation of EGF signaling levels by loss- and gain-of-function studies reveals that EGF signaling negatively regulates embryonic BMP-responsive transcription. Therefore, the BMP gene network has a self-regulating property in that it establishes a balance between its activity and that of the antagonistic EGF signaling pathway to facilitate correct patterning. In terms of BMP-dependent transcription, we identify key roles for the Zelda and Zerknüllt transcription factors in establishing the resulting expression domain, and find widespread binding of insulator proteins to the Mad and Brinker-bound genomic regions. Analysis of embryos lacking the BEAF-32 insulator protein shows reduced transcription of a peak BMP target gene and a reduction in the number of amnioserosa cells, the fate specified by peak BMP signaling. We incorporate our findings into a model for Mad-dependent activation, and discuss its relevance to BMP signal interpretation in vertebrates.

Embryogenesis involves the patterning of many different cell fates by a limited number of types of signals. One way that these signals promote a particular cell fate is through the induction of a complex, yet highly reproducible, gene expression programme that instructs changes in the cell. For example, there is a conserved role for BMP signals in specifying cell fates during dorsal-ventral axis patterning. Here, we have used genomics approaches to identify the gene expression programme implemented in response to BMP signaling during axis patterning in the Drosophila embryo. Part of the gene network downstream of BMP signaling includes members of the EGF signaling pathway, with our data highlighting reciprocal interactions between these two pathways. We have also determined genome-wide binding of BMP-responsive transcription factors to gain new insights into how the BMP gene network is activated. Our data reveal roles for specific transcription factors and insulator binding proteins, with the latter traditionally associated with the separation of transcriptional domains. Overall, our data will provide a platform for exploiting the tractability of the Drosophila embryo to determine which features of the network are critical drivers of BMP-induced cell fate changes during embryogenesis.

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.

Shadow Enhancers Mediate Dynamic Shifts of Gap Gene Expression in the Drosophila Embryo

Drosophila patterning genes often contain pairs of primary and shadow enhancers that possess overlapping activities [1-5]. It has been suggested that this regulatory "redundancy" helps ensure reliable activation of gene expression under stressful conditions such as increases in temperature [4, 5]. There is also evidence that shadow enhancers help produce sharp on/off boundaries of gene expression in response to small changes in the levels of regulatory factors, such as the maternal Bicoid gradient [6, 7]. Here, we use live-imaging methods to visualize the temporal dynamics of the gap genes Kruppel and knirps, which are essential for the patterning of the thorax and abdomen, respectively [8, 9]. Previous analyses of fixed embryos suggested anterior shifts of the Kruppel and knirps expression patterns [10]. Here, we use computational visualization methods to reveal the precise temporal dynamics of these shifts and further suggest that shadow enhancers are crucial for this process. We discuss potential mechanisms for enhancer dominance, whereby one enhancer represses the other to foster temporal dynamics.

Twist-mediated Epithelial-mesenchymal Transition Promotes Breast Tumor Cell Invasion via Inhibition of Hippo Pathway

Twist is a key transcription factor for Epithelial-mesenchymal transition (EMT), which is a cellular de-differentiation program that promotes invasion and metastasis, confers tumor cells with cancer stem cell (CSC)-like characteristics, and increases therapeutic resistance. However, the mechanisms that facilitate the functions of Twist remain unclear. Here we report that Twist overexpression increased expression of PAR1, an upstream regulator of the Hippo pathway; PAR1 promotes invasion, migration, and CSC-like properties in breast cancer by activating the transcriptional co-activator TAZ. Our study indicates that Hippo pathway inhibition is required for the increased migratory and invasiveness ability of breast cancer cells in Twist-mediated EMT.

Understanding transcriptional regulatory networks using computational models

Transcriptional regulatory networks (TRNs) encode instructions for animal development and physiological responses. Recent advances in genomic technologies and computational modeling have revolutionized our ability to construct models of TRNs. Here, we survey current computational methods for inferring TRN models using genome-scale data. We discuss their advantages and limitations. We summarize representative TRNs constructed using genome-scale data in both normal and disease development. We discuss lessons learned about the structure/function relationship of TRNs, based on examining various large-scale TRN models. Finally, we outline some open questions regarding TRNs, including how to improve model accuracy by integrating complementary data types, how to infer condition-specific TRNs, and how to compare TRNs across conditions and species in order to understand their structure/function relationship.

New methods to image transcription in living fly embryos: the insights so far, and the prospects

The regulation of transcription is a fundamental process underlying the determination of cell identity and its maintenance during development. In the last decades, most of the transcription factors, which have to be expressed at the right place and at the right time for the proper development of the fly embryo, have been identified. However, mostly because of the lack of methods to visualize transcription as the embryo develops, their coordinated spatiotemporal dynamics remains largely unexplored. Efforts have been made to decipher the transcription process with single molecule resolution at the single cell level. Recently, the fluorescent labeling of nascent RNA in developing fly embryos allowed the direct visualization of ongoing transcription at single loci within each nucleus. Together with powerful imaging and quantitative data analysis, these new methods provide unprecedented insights into the temporal dynamics of the transcription process and its intrinsic noise. Focusing on the Drosophila embryo, we discuss how the detection of single RNA molecules enhanced our comprehension of the transcription process and we outline the potential next steps made possible by these new imaging tools. In combination with genetics and theoretical analysis, these new imaging methods will aid the search for the mechanisms responsible for the robustness of development. For further resources related to this article, please visit the WIREs website.

Transcription Factor Binding Site Mapping Using ChIP-Seq

Transcription factors (TFs) play a central role in regulating gene expression in all bacteria. Yet until recently, studies of TF binding were limited to a small number of factors at a few genomic locations. Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) provides the ability to map binding sites globally for TFs, and the scalability of the technology enables the ability to map binding sites for every DNA binding protein in a prokaryotic organism. We have developed a protocol for ChIP-Seq tailored for use with mycobacteria and an analysis pipeline for processing the resulting data. The protocol and pipeline have been used to map over 100 TFs from Mycobacterium tuberculosis, as well as numerous TFs from related mycobacteria and other bacteria. The resulting data provide evidence that the long-accepted spatial relationship between TF binding site, promoter motif, and the corresponding regulated gene may be too simple a paradigm, failing to adequately capture the variety of TF binding sites found in prokaryotes. In this article we describe the protocol and analysis pipeline, the validation of these methods, and the results of applying these methods to M. tuberculosis.

Transcriptional elongation requires DNA break-induced signalling

We have previously shown that RNA polymerase II (Pol II) pause release and transcriptional elongation involve phosphorylation of the factor TRIM28 by the DNA damage response (DDR) kinases ATM and DNA-PK. Here we report a significant role for DNA breaks and DDR signalling in the mechanisms of transcriptional elongation in stimulus-inducible genes in humans. Our data show the enrichment of TRIM28 and γH2AX on serum-induced genes and the important function of DNA-PK for Pol II pause release and transcriptional activation-coupled DDR signalling on these genes. γH2AX accumulation decreases when P-TEFb is inhibited, confirming that DDR signalling results from transcriptional elongation. In addition, transcriptional elongation-coupled DDR signalling involves topoisomerase II because inhibiting this enzyme interferes with Pol II pause release and γH2AX accumulation. Our findings propose that DDR signalling is required for effective Pol II pause release and transcriptional elongation through a novel mechanism involving TRIM28, DNA-PK and topoisomerase II.

RNA polymerase II (Pol II) pause release and transcriptional elongation involve phosphorylation of TRIM28 by the DNA damage response (DDR) kinases. Here, Bunch et al. show that DDR signalling is coupled with and required for transcriptional elongation in stimulus-inducible genes and involves topoisomerase II.

Shadow Enhancers Are Pervasive Features of Developmental Regulatory Networks

Embryogenesis is remarkably robust to segregating mutations and environmental variation; under a range of conditions, embryos of a given species develop into stereotypically patterned organisms. Such robustness is thought to be conferred, in part, through elements within regulatory networks that perform similar, redundant tasks. Redundant enhancers (or "shadow" enhancers), for example, can confer precision and robustness to gene expression, at least at individual, well-studied loci. However, the extent to which enhancer redundancy exists and can thereby have a major impact on developmental robustness remains unknown. Here, we systematically assessed this, identifying over 1,000 predicted shadow enhancers during Drosophila mesoderm development. The activity of 23 elements, associated with five genes, was examined in transgenic embryos, while natural structural variation among individuals was used to assess their ability to buffer against genetic variation. Our results reveal three clear properties of enhancer redundancy within developmental systems. First, it is much more pervasive than previously anticipated, with 64% of loci examined having shadow enhancers. Their spatial redundancy is often partial in nature, while the non-overlapping function may explain why these enhancers are maintained within a population. Second, over 70% of loci do not follow the simple situation of having only two shadow enhancers-often there are three (rols), four (CadN and ade5), or five (Traf1), at least one of which can be deleted with no obvious phenotypic effects. Third, although shadow enhancers can buffer variation, patterns of segregating variation suggest that they play a more complex role in development than generally considered.

In Vivo Transcriptional Activation Using CRISPR/Cas9 in Drosophila

A number of approaches for Cas9-mediated transcriptional activation have recently been developed, allowing target genes to be overexpressed from their endogenous genomic loci. However, these approaches have thus far been limited to cell culture, and this technique has not been demonstrated in vivo in any animal. The technique involving the fewest separate components, and therefore the most amenable to in vivo applications, is the dCas9-VPR system, where a nuclease-dead Cas9 is fused to a highly active chimeric activator domain. In this study, we characterize the dCas9-VPR system in Drosophila cells and in vivo. We show that this system can be used in cell culture to upregulate a range of target genes, singly and in multiplex, and that a single guide RNA upstream of the transcription start site can activate high levels of target transcription. We observe marked heterogeneity in guide RNA efficacy for any given gene, and we confirm that transcription is inhibited by guide RNAs binding downstream of the transcription start site. To demonstrate one application of this technique in cells, we used dCas9-VPR to identify target genes for Twist and Snail, two highly conserved transcription factors that cooperate during Drosophila mesoderm development. In addition, we simultaneously activated both Twist and Snail to identify synergistic responses to this physiologically relevant combination. Finally, we show that dCas9-VPR can activate target genes and cause dominant phenotypes in vivo, providing the first demonstration of dCas9 activation in a multicellular animal. Transcriptional activation using dCas9-VPR thus offers a simple and broadly applicable technique for a variety of overexpression studies.

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 Transcriptional Coregulator LEUNIG_HOMOLOG Inhibits Light-Dependent Seed Germination in Arabidopsis

PHYTOCHROME-INTERACTING FACTOR1 (PIF1) is a basic helix-loop-helix transcription factor that inhibits light-dependent seed germination in Arabidopsis thaliana. However, it remains unclear whether PIF1 requires other factors to regulate its direct targets. Here, we demonstrate that LEUNIG_HOMOLOG (LUH), a Groucho family transcriptional corepressor, binds to PIF1 and coregulates its targets. Not only are the transcriptional profiles of the luh and pif1 mutants remarkably similar, more than 80% of the seeds of both genotypes germinate in the dark. We show by chromatin immunoprecipitation that LUH binds a subset of PIF1 targets in a partially PIF1-dependent manner. Unexpectedly, we found LUH binds and coregulates not only PIF1-activated targets but also PIF1-repressed targets. Together, our results indicate LUH functions with PIF1 as a transcriptional coregulator to inhibit seed germination.

Specification of the somatic musculature in Drosophila

The somatic muscle system formed during Drosophila embryogenesis is required for larvae to hatch, feed, and crawl. This system is replaced in the pupa by a new adult muscle set, responsible for activities such as feeding, walking, and flight. Both the larval and adult muscle systems are comprised of distinct muscle fibers to serve these specific motor functions. In this way, the Drosophila musculature is a valuable model for patterning within a single tissue: while all muscle cells share properties such as the contractile apparatus, properties such as size, position, and number of nuclei are unique for a particular muscle. In the embryo, diversification of muscle fibers relies first on signaling cascades that pattern the mesoderm. Subsequently, the combinatorial expression of specific transcription factors leads muscle fibers to adopt particular sizes, shapes, and orientations. Adult muscle precursors (AMPs), set aside during embryonic development, proliferate during the larval phases and seed the formation of the abdominal, leg, and flight muscles in the adult fly. Adult muscle fibers may either be formed de novo from the fusion of the AMPs, or are created by the binding of AMPs to an existing larval muscle. While less is known about adult muscle specification compared to the larva, expression of specific transcription factors is also important for its diversification. Increasingly, the mechanisms required for the diversification of fly muscle have found parallels in vertebrate systems and mark Drosophila as a robust model system to examine questions about how diverse cell types are generated within an organism.

Using transgenic reporter assays to functionally characterize enhancers in animals

Enhancers or cis-regulatory modules play an instructive role in regulating gene expression during animal development and in response to the environment. Despite their importance, we only have an incomplete map of enhancers in the genome and our understanding of the mechanisms governing their function is still limited. Recent advances in genomics provided powerful tools to generate genome-wide maps of potential enhancers. However, most of these methods are based on indirect measures of enhancer activity and have to be followed by functional testing. Animal transgenesis has been a valuable method to functionally test and characterize enhancers in vivo. In this review I discuss how different transgenic strategies are utilized to characterize enhancers in model organisms focusing on studies in Drosophila and mouse. I will further discuss recent large-scale transgenic efforts to systematically identify and catalog enhancers as well as highlight the challenges and future directions in the field.

From shavenbaby to the naked valley: trichome formation as a model for evolutionary developmental biology

Microtrichia or trichomes are non-sensory actin protrusions produced by the epidermal cells of many insects. Studies of trichome formation in Drosophila have over the last 30 years provided key insights towards our understanding of gene regulation, gene regulatory networks (GRNs), development, the genotype to phenotype map, and the evolution of these processes. Here we review classic studies that have used trichome formation as a model to shed light on Drosophila development as well as recent research on the architecture of the GRN underlying trichome formation. This includes the findings that both small peptides and microRNAs play important roles in the regulation and evolution of this network. In addition, we review research on the evolution of trichome patterns that has provided novel insights into the function and architecture of cis-regulatory modules, and into the genetic basis of morphological change. We conclude that further research on these apparently simple and often functionally enigmatic structures will continue to provide new and important knowledge about development and evolution.

Identification of Lineage-Specific Cis-Regulatory Modules Associated with Variation in Transcription Factor Binding and Chromatin Activity Using Ornstein-Uhlenbeck Models

Scoring the impact of noncoding variation on the function of cis-regulatory regions, on their chromatin state, and on the qualitative and quantitative expression levels of target genes is a fundamental problem in evolutionary genomics. A particular challenge is how to model the divergence of quantitative traits and to identify relationships between the changes across the different levels of the genome, the chromatin activity landscape, and the transcriptome. Here, we examine the use of the Ornstein-Uhlenbeck (OU) model to infer selection at the level of predicted cis-regulatory modules (CRMs), and link these with changes in transcription factor binding and chromatin activity. Using publicly available cross-species ChIP-Seq and STARR-Seq data we show how OU can be applied genome-wide to identify candidate transcription factors for which binding site and CRM turnover is correlated with changes in regulatory activity. Next, we profile open chromatin in the developing eye across three Drosophila species. We identify the recognition motifs of the chromatin remodelers, Trithorax-like and Grainyhead as mostly correlating with species-specific changes in open chromatin. In conclusion, we show in this study that CRM scores can be used as quantitative traits and that motif discovery approaches can be extended towards more complex models of divergence.

Islet 1 specifies the identity of hypothalamic melanocortin neurons and is critical for normal food intake and adiposity in adulthood

Food intake and body weight regulation depend on proper expression of the proopiomelanocortin gene (Pomc) in a group of neurons located in the mediobasal hypothalamus of all vertebrates. These neurons release POMC-encoded melanocortins, which are potent anorexigenic neuropeptides, and their absence from mice or humans leads to hyperphagia and severe obesity. Although the pathophysiology of hypothalamic POMC neurons is well understood, the genetic program that establishes the neuronal melanocortinergic phenotype and maintains a fully functional neuronal POMC phenotype throughout adulthood remains unknown. Here, we report that the early expression of the LIM-homeodomain transcription factor Islet 1 (ISL1) in the developing hypothalamus promotes the terminal differentiation of melanocortinergic neurons and is essential for hypothalamic Pomc expression since its initial onset and throughout the entire lifetime. We detected ISL1 in the prospective hypothalamus just before the onset of Pomc expression and, from then on, Pomc and Isl1 coexpress. ISL1 binds in vitro and in vivo to critical homeodomain binding DNA motifs present in the neuronal Pomc enhancers nPE1 and nPE2, and mutations of these sites completely disrupt the ability of these enhancers to drive reporter gene expression to hypothalamic POMC neurons in transgenic mice and zebrafish. ISL1 is necessary for hypothalamic Pomc expression during mouse and zebrafish embryogenesis. Furthermore, conditional Isl1 inactivation from POMC neurons impairs Pomc expression, leading to hyperphagia and obesity. Our results demonstrate that ISL1 specifies the identity of hypothalamic melanocortin neurons and is required for melanocortin-induced satiety and normal adiposity throughout the entire lifespan.

An evolutionarily conserved DNA architecture determines target specificity of the TWIST family bHLH transcription factors

Basic helix-loop-helix (bHLH) transcription factors recognize the canonical E-box (CANNTG) to regulate gene transcription; however, given the prevalence of E-boxes in a genome, it has been puzzling how individual bHLH proteins selectively recognize E-box sequences on their targets. TWIST is a bHLH transcription factor that promotes epithelial-mesenchymal transition (EMT) during development and tumor metastasis. High-resolution mapping of TWIST occupancy in human and Drosophila genomes reveals that TWIST, but not other bHLH proteins, recognizes a unique double E-box motif with two E-boxes spaced preferentially by 5 nucleotides. Using molecular modeling and binding kinetic analyses, we found that the strict spatial configuration in the double E-box motif aligns two TWIST-E47 dimers on the same face of DNA, thus providing a high-affinity site for a highly stable intramolecular tetramer. Biochemical analyses showed that the WR domain of TWIST dimerizes to mediate tetramer formation, which is functionally required for TWIST-induced EMT. These results uncover a novel mechanism for a bHLH transcription factor to recognize a unique spatial configuration of E-boxes to achieve target specificity. The WR-WR domain interaction uncovered here sets an example of target gene specificity of a bHLH protein being controlled allosterically by a domain outside of the bHLH region.

ChIP-nexus enables improved detection of in vivo transcription factor binding footprints

Understanding how eukaryotic enhancers are bound and regulated by specific combinations of transcription factors is still a major challenge. To better map transcription factor binding genome-wide at nucleotide resolution in vivo, we have developed a robust ChIP-exo protocol called ChIP experiments with nucleotide resolution through exonuclease, unique barcode and single ligation (ChIP-nexus), which utilizes an efficient DNA self-circularization step during library preparation. Application of ChIP-nexus to four proteins—human TBP and Drosophila NFkB, Twist and Max— demonstrates that it outperforms existing ChIP protocols in resolution and specificity, pinpoints relevant binding sites within enhancers containing multiple binding motifs and allows the analysis of in vivo binding specificities. Notably, we show that Max frequently interacts with DNA sequences next to its motif, and that this binding pattern correlates with local DNA sequence features such as DNA shape. ChIP-nexus will be broadly applicable to studying in vivo transcription factor binding specificity and its relationship to cis-regulatory changes in humans and model organisms.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

TRIM28 regulates RNA polymerase II promoter-proximal pausing and pause release

Promoter-proximal pausing of RNA polymerase II (Pol II) is a major checkpoint in transcription. An unbiased search for new human proteins that could regulate paused Pol II at the HSPA1B gene identified TRIM28. In vitro analyses indicated HSF1-dependent attenuation of Pol II pausing upon TRIM28 depletion, whereas in vivo data revealed de novo expression of HSPA1B and other known genes regulated by paused Pol II upon TRIM28 knockdown. These results were supported by genome-wide ChIP-sequencing analyses of Pol II occupancy that revealed a global role for TRIM28 in regulating Pol II pausing and pause release. Furthermore, in vivo and in vitro mechanistic studies suggest that transcription-coupled phosphorylation regulates Pol II pause release by TRIM28. Collectively, our findings identify TRIM28 as a new factor that modulates Pol II pausing and transcriptional elongation at a large number of mammalian genes.

Multiscale modeling of dorsoventral patterning in Drosophila

The role of mathematical models of signaling networks is showcased by examples from Drosophila development. Three models of consecutive stages in dorsoventral patterning are presented. We begin with a compartmental model of intracellular reactions that generates a gradient of nuclear-localized Dorsal, exhibiting constant shape and dynamic amplitude. A simple thermodynamic model of equilibrium binding explains how a spatially uniform transcription factor, Zelda, can act in combination with a graded factor, Dorsal, to cooperatively regulate gene expression borders. Finally, we formulate a dynamic and stochastic model that predicts spatiotemporal patterns of Sog expression based on known patterns of its transcription factor, Dorsal. The future of coupling multifarious models across multiple temporal and spatial scales is discussed.

Occupancy of tissue-specific cis-regulatory modules by Otx2 and TLE/Groucho for embryonic head specification

Head specification by the head-selector gene, orthodenticle (otx), is highly conserved among bilaterian lineages. However, the molecular mechanisms by which Otx and other transcription factors (TFs) interact with the genome to direct head formation are largely unknown. Here we employ ChIP-seq and RNA-seq approaches in Xenopus tropicalis gastrulae and find that occupancy of the corepressor, TLE/Groucho, is a better indicator of tissue-specific cis-regulatory modules (CRMs) than the coactivator p300, during early embryonic stages. On the basis of TLE binding and comprehensive CRM profiling, we define two distinct types of Otx2- and TLE-occupied CRMs. Using these devices, Otx2 and other head organizer TFs (for example, Lim1/Lhx1 (activator) or Goosecoid (repressor)) are able to upregulate or downregulate a large battery of target genes in the head organizer. An underlying principle is that Otx marks target genes for head specification to be regulated positively or negatively by partner TFs through specific types of CRMs.

Discrete levels of Twist activity are required to direct distinct cell functions during gastrulation and somatic myogenesis

Twist (Twi), a conserved basic helix-loop-helix transcriptional regulator, directs the epithelial-to-mesenchymal transition (EMT), and regulates changes in cell fate, cell polarity, cell division and cell migration in organisms from flies to humans. Analogous to its role in EMT, Twist has been implicated in metastasis in numerous cancer types, including breast, pancreatic and prostate. In the Drosophila embryo, Twist is essential for discrete events in gastrulation and mesodermal patterning. In this study, we derive a twi allelic series by examining the various cellular events required for gastrulation in Drosophila. By genetically manipulating the levels of Twi activity during gastrulation, we find that coordination of cell division is the most sensitive cellular event, whereas changes in cell shape are the least sensitive. Strikingly, we show that by increasing levels of Snail expression in a severe twi hypomorphic allelic background, but not a twi null background, we can reconstitute gastrulation and produce viable adult flies. Our results demonstrate that the level of Twi activity determines whether the cellular events of ventral furrow formation, EMT, cell division and mesodermal migration occur.

Zelda potentiates morphogen activity by increasing chromatin accessibility

Zygotic genome activation (ZGA) is a major genome programming event whereby the cells of the embryo begin to adopt specified fates. Experiments in Drosophila and zebrafish have revealed that ZGA depends on transcription factors that provide large-scale control of gene expression by direct and specific binding to gene regulatory sequences. Zelda (Zld) plays such a role in the Drosophila embryo, where it has been shown to control the action of patterning signals; however, the mechanisms underlying this effect remain largely unclear. A recent model proposed that Zld binding sites act as quantitative regulators of the spatiotemporal expression of genes activated by Dorsal (Dl), the morphogen that patterns the dorsoventral axis. Here we tested this model experimentally, using enhancers of brinker (brk) and short gastrulation (sog), both of which are directly activated by Dl, but at different concentration thresholds. In agreement with the model, we show that there is a clear positive correlation between the number of Zld binding sites and the spatial domain of enhancer activity. Likewise, the timing of expression could be advanced or delayed. We present evidence that Zld facilitates binding of Dl to regulatory DNA, and that this is associated with increased chromatin accessibility. Importantly, the change in chromatin accessibility is strongly correlated with the change in Zld binding, but not Dl. We propose that the ability of genome activators to facilitate readout of transcriptional input is key to widespread transcriptional induction during ZGA.

Maternal control of the Drosophila dorsal-ventral body axis

The pathway that generates the dorsal-ventral (DV) axis of the Drosophila embryo has been the subject of intense investigation over the previous three decades. The initial asymmetric signal originates during oogenesis by the movement of the oocyte nucleus to an anterior corner of the oocyte, which establishes DV polarity within the follicle through signaling between Gurken, the Drosophila Transforming Growth Factor (TGF)-α homologue secreted from the oocyte, and the Drosophila Epidermal Growth Factor Receptor (EGFR) that is expressed by the follicular epithelium cells that envelop the oocyte. Follicle cells that are not exposed to Gurken follow a ventral fate and express Pipe, a sulfotransferase that enzymatically modifies components of the inner vitelline membrane layer of the eggshell, thereby transferring DV spatial information from the follicle to the egg. These ventrally sulfated eggshell proteins comprise a localized cue that directs the ventrally restricted formation of the active Spätzle ligand within the perivitelline space between the eggshell and the embryonic membrane. Spätzle activates Toll, a transmembrane receptor in the embryonic membrane. Transmission of the Toll signal into the embryo leads to the formation of a ventral-to-dorsal gradient of the transcription factor Dorsal within the nuclei of the syncytial blastoderm stage embryo. Dorsal controls the spatially specific expression of a large constellation of zygotic target genes, the Dorsal gene regulatory network, along the embryonic DV circumference. This article reviews classic studies and integrates them with the details of more recent work that has advanced our understanding of the complex pathway that establishes Drosophila embryo DV polarity. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Enhancer diversity and the control of a simple pattern of Drosophila CNS midline cell expression

Transcriptional enhancers integrate information derived from transcription factor binding to control gene expression. One key question concerns the extent of trans- and cis-regulatory variation in how co-expressed genes are controlled. The Drosophila CNS midline cells constitute a group of neurons and glia in which expression changes can be readily characterized during specification and differentiation. Using a transgenic approach, we compare the cis-regulation of multiple genes expressed in the Drosophila CNS midline primordium cells, and show that while the expression patterns may appear alike, the target genes are not equivalent in how these common expression patterns are achieved. Some genes utilize a single enhancer that promotes expression in all midline cells, while others utilize multiple enhancers with distinct spatial, temporal, and quantitative contributions. Two regulators, Single-minded and Notch, play key roles in controlling early midline gene expression. While Single-minded is expected to control expression of most, if not all, midline primordium-expressed genes, the role of Notch in directly controlling midline transcription is unknown. Midline primordium expression of the rhomboid gene is dependent on cell signaling by the Notch signaling pathway. Mutational analysis of a rhomboid enhancer reveals at least 5 distinct types of functional cis-control elements, including a binding site for the Notch effector, Suppressor of Hairless. The results suggest a model in which Notch/Suppressor of Hairless levels are insufficient to activate rhomboid expression by itself, but does so in conjunction with additional factors, some of which, including Single-minded, provide midline specificity to Notch activation. Similarly, a midline glial enhancer from the argos gene, which is dependent on EGF/Spitz signaling, is directly regulated by contributions from both Pointed, the EGF transcriptional effector, and Single-minded. In contrast, midline primordium expression of other genes shows a strong dependence on Single-minded and varying combinations of additional transcription factors. Thus, Single-minded directly regulates midline primordium-expressed genes, but in some cases plays a primary role in directing target gene midline expression, and in others provides midline specificity to cell signaling inputs.

Impacts of the ubiquitous factor Zelda on Bicoid-dependent DNA binding and transcription in Drosophila

The Drosophila transcription factor Bicoid (Bcd) binds thousands of genomic sites during early embryogenesis, but it is unclear how many of these binding events are functionally important. Here, Small and colleagues test the role of the maternal factor Zelda (Zld) in Bcd-mediated binding and transcription. Embryos lacking Zld show enhanced Bcd binding to a subset of genomic locations, causing early activation of target genes normally silent until later stages. This study demonstrates a critical role for Zld in controlling Bcd binding and target gene activation in the early embryo.

In vivo cross-linking studies suggest that the Drosophila transcription factor Bicoid (Bcd) binds to several thousand sites during early embryogenesis, but it is not clear how many of these binding events are functionally important. In contrast, reporter gene studies have identified >60 Bcd-dependent enhancers, all of which contain clusters of the consensus binding sequence TAATCC. These studies also identified clusters of TAATCC motifs (inactive fragments) that failed to drive Bcd-dependent activation. In general, active fragments showed higher levels of Bcd binding in vivo and were enriched in predicted binding sites for the ubiquitous maternal protein Zelda (Zld). Here we tested the role of Zld in Bcd-mediated binding and transcription. Removal of Zld function and mutations in Zld sites caused significant reductions in Bcd binding to known enhancers and variable effects on the activation and spatial positioning of Bcd-dependent expression patterns. Also, insertion of Zld sites converted one of six inactive fragments into a Bcd-responsive enhancer. Genome-wide binding experiments in zld mutants showed variable effects on Bcd-binding peaks, ranging from strong reductions to significantly enhanced levels of binding. Increases in Bcd binding caused the precocious Bcd-dependent activation of genes that are normally not expressed in early embryos, suggesting that Zld controls the genome-wide binding profile of Bcd at the qualitative level and is critical for selecting target genes for activation in the early embryo. These results underscore the importance of combinatorial binding in enhancer function and provide data that will help predict regulatory activities based on DNA sequence.

A regulatory network of Drosophila germline stem cell self-renewal

Stem cells possess the capacity to generate two cells of distinct fate upon division: one cell retaining stem cell identity and the other cell destined to differentiate. These cell fates are established by cell-type-specific genetic networks. To comprehensively identify components of these networks, we performed a large-scale RNAi screen in Drosophila female germline stem cells (GSCs) covering ∼25% of the genome. The screen identified 366 genes that affect GSC maintenance, differentiation, or other processes involved in oogenesis. Comparison of GSC regulators with neural stem cell self-renewal factors identifies common and cell-type-specific self-renewal genes. Importantly, we identify the histone methyltransferase Set1 as a GSC-specific self-renewal factor. Loss of Set1 in neural stem cells does not affect cell fate decisions, suggesting a differential requirement of H3K4me3 in different stem cell lineages. Altogether, our study provides a resource that will help to further dissect the networks underlying stem cell self-renewal.

Co-activation of microRNAs by Zelda is essential for early Drosophila development

Transcription factors and microRNAs (miRNAs) are two important classes of trans-regulators in differential gene expression. Transcription factors occupy cis-regulatory motifs in DNA to activate or repress gene transcription, whereas miRNAs specifically pair with seed sites in target mRNAs to trigger mRNA decay or inhibit translation. Dynamic spatiotemporal expression patterns of transcription factors and miRNAs during development point to their stage- and tissue-specific functions. Recent studies have focused on miRNA functions during development; however, much remains to explore regarding how the expression of miRNAs is initiated and how dynamic miRNA expression patterns are achieved by transcriptional regulatory networks at different developmental stages. Here, we focused on the identification, regulation and function of miRNAs during the earliest stage of Drosophila development, when the maternal-to-zygotic transition (MZT) takes place. Eleven miRNA clusters comprise the first set of miRNAs activated in the blastoderm embryo. The transcriptional activator Zelda is required for their proper activation and regulation, and Zelda binding observed in genome-wide binding profiles is predictive of enhancer activity. In addition, other blastoderm transcription factors, comprising both activators and repressors, the activities of which are potentiated and coordinated by Zelda, contribute to the accurate temporal and spatial expression of these miRNAs, which are known to function in diverse developmental processes. Although previous genetic studies showed no early phenotypes upon loss of individual miRNAs, our analysis of the miR-1; miR-9a double mutant revealed defects in gastrulation, demonstrating the importance of co-activation of miRNAs by Zelda during the MZT.

Cell type-specific reciprocal regulation of HIF1A gene expression is dependent on 5'- and 3'-UTRs

In the present study, we demonstrated the reciprocal regulation of hypoxia-inducible factor 1 alpha (HIF1A) gene expression via untranslated region-(UTR) dependent mechanisms. A 151 nucleotide sequence found in the HIF1A 5'-UTR is sufficient for significant translational up-regulation. On the other hand, the 3'-UTR of HIF1A has been implicated in mRNA degradation. In the non-metastatic breast cancer cell line MCF7, the 3'-UTR-dependent down-regulatory machinery predominates over the 5'-UTR-dependent up-regulation of HIF1A. However, 5'-UTR-dependent up-regulation is dominant among metastatic cell lines (MDA-MB453, U87MG). It is therefore likely that the predominance of 5'-UTR-dependent translational enhancement of HIF1A is critical for the malignant phenotype of cancer cells. PTBP-1, but not HuR, is a candidate RNA binding protein for the translational control of HIF1A.

In silico evolution of the hunchback gene indicates redundancy in cis-regulatory organization and spatial gene expression

Biological development depends on the coordinated expression of genes in time and space. Developmental genes have extensive cis-regulatory regions which control their expression. These regions are organized in a modular manner, with different modules controlling expression at different times and locations. Both how modularity evolved and what function it serves are open questions. We present a computational model for the cis-regulation of the hunchback (hb) gene in the fruit fly (Drosophila). We simulate evolution (using an evolutionary computation approach from computer science) to find the optimal cis-regulatory arrangements for fitting experimental hb expression patterns. We find that the cis-regulatory region tends to readily evolve modularity. These cis-regulatory modules (CRMs) do not tend to control single spatial domains, but show a multi-CRM/multi-domain correspondence. We find that the CRM-domain correspondence seen in Drosophila evolves with a high probability in our model, supporting the biological relevance of the approach. The partial redundancy resulting from multi-CRM control may confer some biological robustness against corruption of regulatory sequences. The technique developed on hb could readily be applied to other multi-CRM developmental genes.

A graphical model method for integrating multiple sources of genome-scale data

Making effective use of multiple data sources is a major challenge in modern bioinformatics. Genome-wide data such as measures of transcription factor binding, gene expression, and sequence conservation, which are used to identify binding regions and genes that are important to major biological processes such as development and disease, can be difficult to use together due to the different biological meanings and statistical distributions of the heterogeneous data types, but each can provide valuable information for understanding the processes under study. Here we present methods for integrating multiple data sources to gain a more complete picture of gene regulation and expression. Our goal is to identify genes and cis-regulatory regions which play specific biological roles. We describe a graphical mixture model approach for data integration, examine the effect of using different model topologies, and discuss methods for evaluating the effectiveness of the models. Model fitting is computationally efficient and produces results which have clear biological and statistical interpretations. The Hedgehog and Dorsal signaling pathways in Drosophila, which are critical in embryonic development, are used as examples.

Core promoter functions in the regulation of gene expression of Drosophila dorsal target genes

Developmental processes are highly dependent on transcriptional regulation by RNA polymerase II. The RNA polymerase II core promoter is the ultimate target of a multitude of transcription factors that control transcription initiation. Core promoters consist of core promoter motifs, e.g. the initiator, TATA box, and the downstream core promoter element (DPE), which confer specific properties to the core promoter. Here, we explored the importance of core promoter functions in the dorsal-ventral developmental gene regulatory network. This network includes multiple genes that are activated by different nuclear concentrations of Dorsal, an NFκB homolog transcription factor, along the dorsal-ventral axis. We show that over two-thirds of Dorsal target genes contain DPE sequence motifs, which is significantly higher than the proportion of DPE-containing promoters in Drosophila genes. We demonstrate that multiple Dorsal target genes are evolutionarily conserved and functionally dependent on the DPE. Furthermore, we have analyzed the activation of key Dorsal target genes by Dorsal, as well as by another Rel family transcription factor, Relish, and the dependence of their activation on the DPE motif. Using hybrid enhancer-promoter constructs in Drosophila cells and embryo extracts, we have demonstrated that the core promoter composition is an important determinant of transcriptional activity of Dorsal target genes. Taken together, our results provide evidence for the importance of core promoter composition in the regulation of Dorsal target genes.

Transcriptional enhancers: from properties to genome-wide predictions

Cellular development, morphology and function are governed by precise patterns of gene expression. These are established by the coordinated action of genomic regulatory elements known as enhancers or cis-regulatory modules. More than 30 years after the initial discovery of enhancers, many of their properties have been elucidated; however, despite major efforts, we only have an incomplete picture of enhancers in animal genomes. In this Review, we discuss how properties of enhancer sequences and chromatin are used to predict enhancers in genome-wide studies. We also cover recently developed high-throughput methods that allow the direct testing and identification of enhancers on the basis of their activity. Finally, we discuss recent technological advances and current challenges in the field of regulatory genomics.

Disrupting the interaction of BRD4 with diacetylated Twist suppresses tumorigenesis in basal-like breast cancer

Twist is a key transcription activator of epithelial-mesenchymal transition (EMT). It remains unclear how Twist induces gene expression. Here we report a mechanism by which Twist recruits BRD4 to direct WNT5A expression in basal-like breast cancer (BLBC). Twist contains a "histone H4-mimic" GK-X-GK motif that is diacetylated by Tip60. The diacetylated Twist binds the second bromodomain of BRD4, whose first bromodomain interacts with acetylated H4, thereby constructing an activated Twist/BRD4/P-TEFb/RNA-Pol II complex at the WNT5A promoter and enhancer. Pharmacologic inhibition of the Twist-BRD4 association reduced WNT5A expression and suppressed invasion, cancer stem cell (CSC)-like properties, and tumorigenicity of BLBC cells. Our study indicates that the interaction with BRD4 is critical for the oncogenic function of Twist in BLBC.

A conserved role for Snail as a potentiator of active transcription

The transcription factors of the Snail family are key regulators of epithelial-mesenchymal transitions, cell morphogenesis, and tumor metastasis. Since its discovery in Drosophila ∼25 years ago, Snail has been extensively studied for its role as a transcriptional repressor. Here we demonstrate that Drosophila Snail can positively modulate transcriptional activation. By combining information on in vivo occupancy with expression profiling of hand-selected, staged snail mutant embryos, we identified 106 genes that are potentially directly regulated by Snail during mesoderm development. In addition to the expected Snail-repressed genes, almost 50% of Snail targets showed an unanticipated activation. The majority of "Snail-activated" genes have enhancer elements cobound by Twist and are expressed in the mesoderm at the stages of Snail occupancy. Snail can potentiate Twist-mediated enhancer activation in vitro and is essential for enhancer activity in vivo. Using a machine learning approach, we show that differentially enriched motifs are sufficient to predict Snail's regulatory response. In silico mutagenesis revealed a likely causative motif, which we demonstrate is essential for enhancer activation. Taken together, these data indicate that Snail can potentiate enhancer activation by collaborating with different activators, providing a new mechanism by which Snail regulates development.

Subtle changes in motif positioning cause tissue-specific effects on robustness of an enhancer's activity

Deciphering the specific contribution of individual motifs within cis-regulatory modules (CRMs) is crucial to understanding how gene expression is regulated and how this process is affected by sequence variation. But despite vast improvements in the ability to identify where transcription factors (TFs) bind throughout the genome, we are limited in our ability to relate information on motif occupancy to function from sequence alone. Here, we engineered 63 synthetic CRMs to systematically assess the relationship between variation in the content and spacing of motifs within CRMs to CRM activity during development using Drosophila transgenic embryos. In over half the cases, very simple elements containing only one or two types of TF binding motifs were capable of driving specific spatio-temporal patterns during development. Different motif organizations provide different degrees of robustness to enhancer activity, ranging from binary on-off responses to more subtle effects including embryo-to-embryo and within-embryo variation. By quantifying the effects of subtle changes in motif organization, we were able to model biophysical rules that explain CRM behavior and may contribute to the spatial positioning of CRM activity in vivo. For the same enhancer, the effects of small differences in motif positions varied in developmentally related tissues, suggesting that gene expression may be more susceptible to sequence variation in one tissue compared to another. This result has important implications for human eQTL studies in which many associated mutations are found in cis-regulatory regions, though the mechanism for how they affect tissue-specific gene expression is often not understood.

Dynamic chromatin modification sustains epithelial-mesenchymal transition following inducible expression of Snail-1

Epithelial-mesenchymal transition (EMT) is thought to contribute to cancer metastasis, but its underlying mechanisms are not well understood. To define early steps in this cellular transformation, we analyzed human mammary epithelial cells with tightly regulated expression of Snail-1, a master regulator of EMT. After Snail-1 induction, epithelial markers were repressed within 6 hr, and mesenchymal genes were induced at 24 hr. Snail-1 binding to its target promoters was transient (6–48 hr) despite continued protein expression, and it was followed by both transient and long-lasting chromatin changes. Pharmacological inhibition of selected histone acetylation and demethylation pathways suppressed the induction as well as the maintenance of Snail-1-mediated EMT. Thus, EMT involves an epigenetic switch that may be prevented or reversed with the use of small-molecule inhibitors of chromatin modifiers.

An integrative analysis reveals functional targets of GATA6 transcriptional regulation in gastric cancer

Lineage-restricted transcription factors (TFs) are frequently mutated or overexpressed in cancer and contribute toward malignant behaviors; however, the molecular bases of their oncogenic properties are largely unknown. As TF activities are difficult to inhibit directly with small molecules, the genes and pathways they regulate might represent more tractable targets for drug therapy. We studied GATA6, a TF gene that is frequently amplified or overexpressed in gastric, esophageal and pancreatic adenocarcinomas. GATA6-overexpressing gastric cancer cell lines cluster in gene expression space, separate from non-overexpressing lines. This expression clustering signifies a shared pathogenic group of genes that GATA6 may regulate through direct cis-element binding. We used chromatin immunoprecipitation and sequencing (ChIP-seq) to identify GATA6-bound genes and considered TF occupancy in relation to genes that respond to GATA6 depletion in cell lines and track with GATA6 mRNA (synexpression groups) in primary gastric cancers. Among other cellular functions, GATA6-occupied genes control apoptosis and govern the M-phase of the cell cycle. Depletion of GATA6 reduced the levels of the latter transcripts and arrested cells in G2 and M phases of the cell cycle. Synexpression in human tumor samples identified likely direct transcriptional targets substantially better than consideration only of transcripts that respond to GATA6 loss in cultured cells. Candidate target genes responded to the loss of GATA6 or its homolog GATA4 and even more to the depletion of both proteins. Many GATA6-dependent genes lacked nearby binding sites but several strongly dependent, synexpressed and GATA6-bound genes encode TFs such as MYC, HES1, RARB and CDX2. Thus, many downstream effects occur indirectly through other TFs and GATA6 activity in gastric cancer is partially redundant with GATA4. This integrative analysis of locus occupancy, gene dependency and synexpression provides a functional signature of GATA6-overexpressing gastric cancers, revealing both limits and new therapeutic directions for a challenging and frequently fatal disease.