High resolution mapping of Twist to DNA in Drosophila embryos: Efficient functional analysis and evolutionary conservation

Sequence reliance of the Drosophila context-dependent transcription factor CLAMP

Despite binding similar cis elements in multiple locations, a single transcription factor (TF) often performs context-dependent functions at different loci. How factors integrate cis sequence and genomic context is still poorly understood and has implications for off-target effects in genetic engineering. The Drosophila context-dependent TF chromatin-linked adaptor for male-specific lethal proteins (CLAMP) targets similar GA-rich cis elements on the X-chromosome and at the histone gene locus but recruits very different, locus-specific factors. We discover that CLAMP leverages information from both cis element and local sequence to perform context-specific functions. Our observations imply the importance of other cues, including protein-protein interactions and the presence of additional cofactors.

Protein-intrinsic properties and context-dependent effects regulate pioneer factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type-specific binding and activity. The mechanisms governing pioneer factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

Transcriptional Targets of TWIST1 in Human Mesenchymal Stem/Stromal Cells Mechanistically Link Stem/Progenitor and Paracrine Functions

Despite extensive clinical testing, mesenchymal stem/stromal cell (MSC)-based therapies continue to underperform with respect to efficacy, which reflects the paucity of biomarkers that predict potency prior to patient administration. Previously, we reported that TWIST1 predicts inter-donor differences in MSC quality attributes that confer potency. To define the full spectrum of TWIST1 activity in MSCs, the present work employed integrated omics-based profiling to identify a high-confidence set of TWIST1 targets, which mapped to cellular processes related to ECM structure/organization, skeletal and circulatory system development, interferon gamma signaling, and inflammation. These targets are implicated in contributing to both stem/progenitor and paracrine activities of MSCs indicating these processes are linked mechanistically in a TWIST1-dependent manner. Targets implicated in extracellular matrix dynamics further implicate TWIST1 in modulating cellular responses to niche remodeling. Novel TWIST1-regulated genes identified herein may be prioritized for future mechanistic and functional studies.

Sequence reliance of a Drosophila context-dependent transcription factor

Despite binding similar cis elements in multiple locations, a single transcription factor often performs context-dependent functions at different loci. How factors integrate cis sequence and genomic context is still poorly understood and has implications for off-target effects in genetic engineering. The Drosophila context-dependent transcription factor CLAMP targets similar GA-rich cis elements on the X-chromosome and at the histone gene locus but recruits very different, loci-specific factors. We discover that CLAMP leverages information from both cis element and local sequence to perform context-specific functions. Our observations imply the importance of other cues, including protein-protein interactions and the presence of additional cofactors.

Protein-intrinsic properties and context-dependent effects regulate pioneer-factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type specific binding and activity. The mechanisms governing pioneer-factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head, and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin, and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

Contribution of TWIST1-EVX1 Axis in Invasiveness of Esophageal Squamous Cell Carcinoma; a Functional Study

BACKGROUND

Epithelial-mesenchymal transition (EMT) is a biological process in embryonic development and cancer progression, and different gene families, such as HOX genes, are closely related to this process.

OBJECTIVES

Our aim in this study was to investigate the correlation between TWIST1 and EVX1 mRNA expression in ESCC patients and also examine the probable regulatory function of TWIST1 on EVX1 expression in human ESCC cell line.

MATERIALS AND METHODS

TWIST1 and EVX1 gene expression patterns were assessed in ESCC patients by relative comparative Real-time PCR then correlated with their clinical characteristics. In silico analysis of the EVX1 gene was conducted. KYSE-30 cells were transduced by a retroviral system to ectopically express TWIST1, followed by qRT-PCR to reveal the correlation between TWIST1 and EVX1 gene expression.

RESULTS

The expression of TWIST1 and EVX1 was correlated to each other significantly (p=0.005) in ESCC. Of 28 patients with under/normal expression of TWIST1, 22 samples (78.57%) had over/normal expression of EVX1. TWIST1 overexpression was correlated with advanced stages of the tumor (III, IV) (P = 0.019) and lymph node metastasis. However, EVX1 under expression was associated with lymph node metastasis (p = 0.027) and invasiveness of the disease (P = 0.037) in ESCC. Furthermore, retroviral transduction enforced significant overexpression of TWIST1 in GFP-hTWIST-1 approximately 9-fold compared to GFP control cells, causing a - 8.83- fold reduction in EVX1 mRNA expression significantly.

CONCLUSIONS

Our results indicated the repressive role of TWIST1 on EVX1 gene expression in ESCC. Therefore, our findings can help dissect the molecular mechanism of ESCC tumorigenesis and discover novel therapeutic targets for ESCC invasion and metastasis.

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

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

Mechanisms of Binding Specificity among bHLH Transcription Factors

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

Dynamic patterning by morphogens illuminated by cis-regulatory studies

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

Functional Transcription Factor Target Networks Illuminate Control of Epithelial Remodelling

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

Odd-paired is a pioneer-like factor that coordinates with Zelda to control gene expression in embryos

Pioneer factors such as Zelda (Zld) help initiate zygotic transcription in Drosophila early embryos, but whether other factors support this dynamic process is unclear. Odd-paired (Opa), a zinc-finger transcription factor expressed at cellularization, controls the transition of genes from pair-rule to segmental patterns along the anterior-posterior axis. Finding that Opa also regulates expression through enhancer sog_Distal along the dorso-ventral axis, we hypothesized Opa's role is more general. Chromatin-immunoprecipitation (ChIP-seq) confirmed its in vivo binding to sog_Distal but also identified widespread binding throughout the genome, comparable to Zld. Furthermore, chromatin assays (ATAC-seq) demonstrate that Opa, like Zld, influences chromatin accessibility genome-wide at cellularization, suggesting both are pioneer factors with common as well as distinct targets. Lastly, embryos lacking opa exhibit widespread, late patterning defects spanning both axes. Collectively, these data suggest Opa is a general timing factor and likely late-acting pioneer factor that drives a secondary wave of zygotic gene expression.

Twist-dependent ratchet functioning downstream from Dorsal revealed using a light-inducible degron

Graded transcription factors are pivotal regulators of embryonic patterning, but whether their role changes over time is unclear. A light-regulated protein degradation system was used to assay temporal dependence of the transcription factor Dorsal in dorsal-ventral axis patterning of Drosophila embryos. Surprisingly, the high-threshold target gene snail only requires Dorsal input early but not late when Dorsal levels peak. Instead, late snail expression can be supported by action of the Twist transcription factor, specifically, through one enhancer, sna.distal This study demonstrates that continuous input is not required for some Dorsal targets and downstream responses, such as twist, function as molecular ratchets.

Spatiotemporal control of gene expression boundaries using a feedforward loop

BACKGROUND

A feedforward loop (FFL) is commonly observed in several biological networks. The FFL network motif has been mostly studied with respect to variation of the input signal in time, with only a few studies of FFL activity in a spatially distributed system such as morphogen-mediated tissue patterning. However, most morphogen gradients also evolve in time.

RESULTS

We studied the spatiotemporal behavior of a coherent FFL in two contexts: (a) a generic, oscillating morphogen gradient and (b) the dorsal-ventral patterning of the early Drosophila embryo by a gradient of the NF-κB homolog dorsal with its early target Twist. In both models, we found features in the dynamics of the intermediate node-phase difference and noise filtering-that were largely independent of the parameterization of the models, and thus were functions of the structure of the FFL itself. In the dorsal gradient model, we also found that proper target gene expression was not possible without including the effect of maternal pioneer factor Zelda.

CONCLUSIONS

An FFL buffers fluctuation to changes in the morphogen signal ensuring stable gene expression boundaries.

Coacting enhancers can have complementary functions within gene regulatory networks and promote canalization

Developmental genes are often regulated by multiple enhancers exhibiting similar spatiotemporal outputs, which are generally considered redundantly acting though few have been studied functionally. Using CRISPR-Cas9, we created deletions of two enhancers, brk5' and brk3', that drive similar but not identical expression of the gene brinker (brk) in early Drosophila embryos. Utilizing both in situ hybridization and quantitative mRNA analysis, we investigated the changes in the gene network state caused by the removal of one or both of the early acting enhancers. brk5' deletion generally phenocopied the gene mutant, including expansion of the BMP ligand decapentaplegic (dpp) as well as inducing variability in amnioserosa tissue cell number suggesting a loss of canalization. In contrast, brk3' deletion presented unique phenotypes including dorsal expansion of several ventrally expressed genes and a decrease in amnioserosa cell number. Similarly, deletions were made for two enhancers associated with the gene short-gastrulation (sog), sog.int and sog.dist, demonstrating that they also exhibit distinct patterning phenotypes and affect canalization. In summary, this study shows that similar gene expression driven by coacting enhancers can support distinct, and sometimes complementary, functions within gene regulatory networks and, moreover, that phenotypes associated with individual enhancer deletion mutants can provide insight into new gene functions.

Migrating cells control morphogenesis of substratum serving as track to promote directional movement of the collective

In Drosophila embryos, caudal visceral mesoderm (CVM) cells undergo bilateral migration along the trunk visceral mesoderm (TVM) in order to form midgut muscles. Mutation of FGF receptor Heartless (Htl) has been shown to cause CVM migration defects, particularly midline crossing of the bilateral groups. Here, we show that htl mutants also exhibit TVM defects including contralateral merging. Both CVM mismigration and TVM contralateral merging are attenuated by restoring FGF signaling specifically in the CVM, suggesting that migrating CVM cells influence TVM morphogenesis; however, the inverse, supplying FGF to the TVM, does not rescue CVM mismigration. In addition, we show that FGF regulates integrin expression in both tissues, but only providing a source of integrin specifically to the TVM attenuates the contralateral merging phenotype. Finally, we demonstrate that the CVM influences cell shape in the TVM, and a loss of CVM results in TVM morphological defects. In summary, this study provides insight into how a migrating collective of cells can influence their tissue substrate and supports the view that morphogenesis of tissues during development is interdependent.

Broadly expressed repressors integrate patterning across orthogonal axes in embryos

The role of spatially localized repressors in supporting embryonic patterning is well appreciated, but, alternatively, the role ubiquitously expressed repressors play in this process is not well understood. We investigated the function of two broadly expressed repressors, Runt (Run) and Suppressor of Hairless [Su(H)], in patterning the Drosophila embryo. Previous studies have shown that Run and Su(H) regulate gene expression along anterior-posterior (AP) or dorsal-ventral (DV) axes, respectively, by spatially limiting activator action, but here we characterize a different role. Our data show that broadly expressed repressors silence particular enhancers within cis-regulatory systems, blocking their expression throughout the embryo fully but transiently, and, in this manner, regulate spatiotemporal outputs along both axes. Our results suggest that Run and Su(H) regulate the temporal action of enhancers and are not dedicated regulators of one axis but, instead, act coordinately to pattern both axes, AP and DV.

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.

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.

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.

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

DOI:http://dx.doi.org/10.7554/eLife.08445.001

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

DOI:http://dx.doi.org/10.7554/eLife.08445.002

Deciphering the molecular mechanisms underlying the binding of the TWIST1/E12 complex to regulatory E-box sequences

The TWIST1 bHLH transcription factor controls embryonic development and cancer processes. Although molecular and genetic analyses have provided a wealth of data on the role of bHLH transcription factors, very little is known on the molecular mechanisms underlying their binding affinity to the E-box sequence of the promoter. Here, we used an in silico model of the TWIST1/E12 (TE) heterocomplex and performed molecular dynamics (MD) simulations of its binding to specific (TE-box) and modified E-box sequences. We focused on (i) active E-box and inactive E-box sequences, on (ii) modified active E-box sequences, as well as on (iii) two box sequences with modified adjacent bases the AT- and TA-boxes. Our in silico models were supported by functional in vitro binding assays. This exploration highlighted the predominant role of protein side-chain residues, close to the heart of the complex, at anchoring the dimer to DNA sequences, and unveiled a shift towards adjacent ((-1) and (-1*)) bases and conserved bases of modified E-box sequences. In conclusion, our study provides proof of the predictive value of these MD simulations, which may contribute to the characterization of specific inhibitors by docking approaches, and their use in pharmacological therapies by blocking the tumoral TWIST1/E12 function in cancers.

A Clinical Indications Prediction Scale Based on TWIST1 for Human Mesenchymal Stem Cells

In addition to their stem/progenitor properties, mesenchymal stem cells (MSCs) also exhibit potent effector (angiogenic, antiinflammatory, immuno-modulatory) functions that are largely paracrine in nature. It is widely believed that effector functions underlie most of the therapeutic potential of MSCs and are independent of their stem/progenitor properties. Here we demonstrate that stem/progenitor and effector functions are coordinately regulated at the cellular level by the transcription factor Twist1 and specified within populations according to a hierarchical model. We further show that manipulation of Twist1 levels by genetic approaches or by exposure to widely used culture supplements including fibroblast growth factor 2 (Ffg2) and interferon gamma (IFN-gamma) alters MSC efficacy in cell-based and in vivo assays in a predictable manner. Thus, by mechanistically linking stem/progenitor and effector functions our studies provide a unifying framework in the form of an MSC hierarchy that models the functional complexity of populations. Using this framework, we developed a CLinical Indications Prediction (CLIP) scale that predicts how donor-to-donor heterogeneity and culture conditions impact the therapeutic efficacy of MSC populations for different disease indications.

•

Mesenchymal stem cells exhibit stem/progenitor and effector (angiogenic, anti-inflammatory, immuno-modulatory) functions.

•

Twist1 coordinately regulates stem/progenitor and effector functions, which are specified hierarchically in populations.

•

Twist1 levels predict inter-population differences in therapeutic efficacy of mesenchymal stem cells for different disease indications.

Mesenchymal stem cells are being evaluated in human clinical trials for treating ischemic, inflammatory, and immunological diseases. However, most completed trials have yielded suboptimal outcomes due to the inability to predict the potency of different donor populations, which are functionally heterogeneous. We demonstrate that clinically relevant biological activities of MSCs are coordinately regulated by the transcription factor Twist1. Furthermore, we showed that TWIST1 levels reliably predict differences in the angiogenic, anti-inflammatory, and immuno-modulatory activity of populations and as such used it to develop a Clinical Indications Prediction (CLIP) scale. By predicting potency of MSC populations for different disease indications the CLIP scale is expected to dramatically improve MSC-based clinical trial outcomes.

The presence of nuclear cactus in the early Drosophila embryo may extend the dynamic range of the dorsal gradient

In a developing embryo, the spatial distribution of a signaling molecule, or a morphogen gradient, has been hypothesized to carry positional information to pattern tissues. Recent measurements of morphogen distribution have allowed us to subject this hypothesis to rigorous physical testing. In the early Drosophila embryo, measurements of the morphogen Dorsal, which is a transcription factor responsible for initiating the earliest zygotic patterns along the dorsal-ventral axis, have revealed a gradient that is too narrow to pattern the entire axis. In this study, we use a mathematical model of Dorsal dynamics, fit to experimental data, to determine the ability of the Dorsal gradient to regulate gene expression across the entire dorsal-ventral axis. We found that two assumptions are required for the model to match experimental data in both Dorsal distribution and gene expression patterns. First, we assume that Cactus, an inhibitor that binds to Dorsal and prevents it from entering the nuclei, must itself be present in the nuclei. And second, we assume that fluorescence measurements of Dorsal reflect both free Dorsal and Cactus-bound Dorsal. Our model explains the dynamic behavior of the Dorsal gradient at lateral and dorsal positions of the embryo, the ability of Dorsal to regulate gene expression across the entire dorsal-ventral axis, and the robustness of gene expression to stochastic effects. Our results have a general implication for interpreting fluorescence-based measurements of signaling molecules.

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.

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

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

Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis

We used mouse ENCODE data along with complementary data from other laboratories to study the dynamics of occupancy and the role in gene regulation of the transcription factor TAL1, a critical regulator of hematopoiesis, at multiple stages of hematopoietic differentiation. We combined ChIP-seq and RNA-seq data in six mouse cell types representing a progression from multilineage precursors to differentiated erythroblasts and megakaryocytes. We found that sites of occupancy shift dramatically during commitment to the erythroid lineage, vary further during terminal maturation, and are strongly associated with changes in gene expression. In multilineage progenitors, the likely target genes are enriched for hematopoietic growth and functions associated with the mature cells of specific daughter lineages (such as megakaryocytes). In contrast, target genes in erythroblasts are specifically enriched for red cell functions. Furthermore, shifts in TAL1 occupancy during erythroid differentiation are associated with gene repression (dissociation) and induction (co-occupancy with GATA1). Based on both enrichment for transcription factor binding site motifs and co-occupancy determined by ChIP-seq, recruitment by GATA transcription factors appears to be a stronger determinant of TAL1 binding to chromatin than the canonical E-box binding site motif. Studies of additional proteins lead to the model that TAL1 regulates expression after being directed to a distinct subset of genomic binding sites in each cell type via its association with different complexes containing master regulators such as GATA2, ERG, and RUNX1 in multilineage cells and the lineage-specific master regulator GATA1 in erythroblasts.

Autoregulatory feedback controls sequential action of cis-regulatory modules at the brinker locus

cis-regulatory modules (CRMs) act sequentially to regulate temporal expression of genes, but how the switch from one to the next is accomplished is not well understood. To provide insight, here we investigate the cis-regulatory system controlling brinker (brk) expression in the Drosophila embryo. Two distally located CRMs support expression at different times, while a promoter-proximal element (PPE) is required to support their action. In the absence of Brk protein itself or upon mutagenesis of Brk binding sites within the PPE, the late-acting CRM, specifically, is delayed. This block to late-acting CRM function appears to be removed when the early-acting CRM is also deleted. These results demonstrate that autoregulatory feedback is necessary for the early-acting CRM to disengage from the promoter so that the late-acting CRM may act. Autoregulation may be a commonly used mechanism to control sequential CRM action necessary for dynamic gene expression throughout the course of development.

•

Two CRMs support spatiotemporally distinct expression of brinker in early embryo

•

A promoter-proximal element (PPE) is required to support action of both CRMs

•

Brk protein acts at the PPE to facilitate the switch from early- to late-acting CRM

•

Autoregulation may be a common mechanism used to control sequential CRM action

Genome-wide analyses of Shavenbaby target genes reveals distinct features of enhancer organization

Background

Developmental programs are implemented by regulatory interactions between Transcription Factors (TFs) and their target genes, which remain poorly understood. While recent studies have focused on regulatory cascades of TFs that govern early development, little is known about how the ultimate effectors of cell differentiation are selected and controlled. We addressed this question during late Drosophila embryogenesis, when the finely tuned expression of the TF Ovo/Shavenbaby (Svb) triggers the morphological differentiation of epidermal trichomes.

Results

We defined a sizeable set of genes downstream of Svb and used in vivo assays to delineate 14 enhancers driving their specific expression in trichome cells. Coupling computational modeling to functional dissection, we investigated the regulatory logic of these enhancers. Extending the repertoire of epidermal effectors using genome-wide approaches showed that the regulatory models learned from this first sample are representative of the whole set of trichome enhancers. These enhancers harbor remarkable features with respect to their functional architectures, including a weak or non-existent clustering of Svb binding sites. The in vivo function of each site relies on its intimate context, notably the flanking nucleotides. Two additional cis-regulatory motifs, present in a broad diversity of composition and positioning among trichome enhancers, critically contribute to enhancer activity.

Conclusions

Our results show that Svb directly regulates a large set of terminal effectors of the remodeling of epidermal cells. Further, these data reveal that trichome formation is underpinned by unexpectedly diverse modes of regulation, providing fresh insights into the functional architecture of enhancers governing a terminal differentiation program.

Drosophila functional elements are embedded in structurally constrained sequences

Modern functional genomics uncovered numerous functional elements in metazoan genomes. Nevertheless, only a small fraction of the typical non-exonic genome contains elements that code for function directly. On the other hand, a much larger fraction of the genome is associated with significant evolutionary constraints, suggesting that much of the non-exonic genome is weakly functional. Here we show that in flies, local (30–70 bp) conserved sequence elements that are associated with multiple regulatory functions serve as focal points to a pattern of punctuated regional increase in G/C nucleotide frequencies. We show that this pattern, which covers a region tenfold larger than the conserved elements themselves, is an evolutionary consequence of a shift in the balance between gain and loss of G/C nucleotides and that it is correlated with nucleosome occupancy across multiple classes of epigenetic state. Evidence for compensatory evolution and analysis of SNP allele frequencies show that the evolutionary regime underlying this balance shift is likely to be non-neutral. These data suggest that current gaps in our understanding of genome function and evolutionary dynamics are explicable by a model of sparse sequence elements directly encoding for function, embedded into structural sequences that help to define the local and global epigenomic context of such functional elements.

A key challenge in functional genomics is to predict evolutionary dynamics from functional annotation of the genome and vice versa. Modern epigenomic studies helped assign function to numerous new sequence elements, but left most of the genome essentially uncharacterized. Evolutionary genomics, on the other hand, consistently suggests that a much larger fraction of the un-annotated genome evolves under selective pressure. We hypothesize that this function-selection gap can be attributed to sequences that facilitate the physical organization of functional elements, such as transcription factor binding sites, within chromosomes. We exemplify this by studying in detail the sequences embedding small conserved elements (CEs) in Drosophila. We show that, while CEs have typically high AT content, high GC content levels around them are maintained by a non-neutral evolutionary balance between gain and loss of GC nucleotides. This non-uniform pattern is highly correlated with nucleosome organization around CEs, potentially imposing an evolutionary constraint on as much as one quarter of the genome. We suggest this can at least partly explain the above function-selection gap. Weak evolutionary constraints on “structural” sequences (at scales ranging from one nucleosome to recently described multi-megabase topological domains) may affect genome evolution just like structural motifs shape protein evolution.

ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

Chromatin immunoprecipitation (ChIP) followed by high-throughput DNA sequencing (ChIP-seq) has become a valuable and widely used approach for mapping the genomic location of transcription-factor binding and histone modifications in living cells. Despite its widespread use, there are considerable differences in how these experiments are conducted, how the results are scored and evaluated for quality, and how the data and metadata are archived for public use. These practices affect the quality and utility of any global ChIP experiment. Through our experience in performing ChIP-seq experiments, the ENCODE and modENCODE consortia have developed a set of working standards and guidelines for ChIP experiments that are updated routinely. The current guidelines address antibody validation, experimental replication, sequencing depth, data and metadata reporting, and data quality assessment. We discuss how ChIP quality, assessed in these ways, affects different uses of ChIP-seq data. All data sets used in the analysis have been deposited for public viewing and downloading at the ENCODE (http://encodeproject.org/ENCODE/) and modENCODE (http://www.modencode.org/) portals.

Genome-wide screens for in vivo Tinman binding sites identify cardiac enhancers with diverse functional architectures

The NK homeodomain factor Tinman is a crucial regulator of early mesoderm patterning and, together with the GATA factor Pannier and the Dorsocross T-box factors, serves as one of the key cardiogenic factors during specification and differentiation of heart cells. Although the basic framework of regulatory interactions driving heart development has been worked out, only about a dozen genes involved in heart development have been designated as direct Tinman target genes to date, and detailed information about the functional architectures of their cardiac enhancers is lacking. We have used immunoprecipitation of chromatin (ChIP) from embryos at two different stages of early cardiogenesis to obtain a global overview of the sequences bound by Tinman in vivo and their linked genes. Our data from the analysis of ∼50 sequences with high Tinman occupancy show that the majority of such sequences act as enhancers in various mesodermal tissues in which Tinman is active. All of the dorsal mesodermal and cardiac enhancers, but not some of the others, require tinman function. The cardiac enhancers feature diverse arrangements of binding motifs for Tinman, Pannier, and Dorsocross. By employing these cardiac and non-cardiac enhancers in machine learning approaches, we identify a novel motif, termed CEE, as a classifier for cardiac enhancers. In vivo assays for the requirement of the binding motifs of Tinman, Pannier, and Dorsocross, as well as the CEE motifs in a set of cardiac enhancers, show that the Tinman sites are essential in all but one of the tested enhancers; although on occasion they can be functionally redundant with Dorsocross sites. The enhancers differ widely with respect to their requirement for Pannier, Dorsocross, and CEE sites, which we ascribe to their different position in the regulatory circuitry, their distinct temporal and spatial activities during cardiogenesis, and functional redundancies among different factor binding sites.

The Drosophila homeodomain protein Tinman was the first transcription factor found to control the development and differentiation of the heart in any species. In spite of that, our knowledge of the number, identities, and mode of regulation of the downstream target genes of Tinman that are necessary to exert its cardiogenic functions is still very incomplete. To address these issues, we have performed a genome-wide analysis of DNA regions associated with Tinman-binding in embryos and the genes linked to them. The combined data from our in-depth in vivo assays of sequence elements with high Tinman occupancy allow the following general conclusions: (1) The majority of such sequences are active as regulatory elements (called enhancers) in mesodermal tissues that include Tinman-expressing cells. (2) The enhancers active in the heart progenitor cells and the heart generally are dependent on tinman gene activity, whereas those active in non-cardiac mesoderm are often bound neutrally by Tinman. (3) Tinman binding motifs in most cases are essential for cardiac enhancer activity, but in some cases they can be functionally-redundant with those of other cardiogenic factors. (4) Tinman-occupied cardiac enhancers are enriched for a newly discovered binding motif for an unknown factor that is functional in vivo.

CTF: a CRF-based transcription factor binding sites finding system

Background

Identifying the location of transcription factor bindings is crucial to understand transcriptional regulation. Currently, Chromatin Immunoprecipitation followed with high-throughput Sequencing (ChIP-seq) is able to locate the transcription factor binding sites (TFBSs) accurately in high throughput and it has become the gold-standard method for TFBS finding experimentally. However, due to its high cost, it is impractical to apply the method in a very large scale. Considering the large number of transcription factors, numerous cell types and various conditions, computational methods are still very valuable to accurate TFBS identification.

Results

In this paper, we proposed a novel integrated TFBS prediction system, CTF, based on Conditional Random Fields (CRFs). Integrating information from different sources, CTF was able to capture patterns of TFBSs contained in different features (sequence, chromatin and etc) and predicted the TFBS locations with a high accuracy. We compared CTF with several existing tools as well as the PWM baseline method on a dataset generated by ChIP-seq experiments (TFBSs of 13 transcription factors in mouse genome). Results showed that CTF performed significantly better than existing methods tested.

Conclusions

CTF is a powerful tool to predict TFBSs by integrating high throughput data and different features. It can be a useful complement to ChIP-seq and other experimental methods for TFBS identification and thus improve our ability to investigate functional elements in post-genomic era.

Availability: CTF is freely available to academic users at: http://cbb.sjtu.edu.cn/~ccwei/pub/software/CTF/CTF.php

Twist1 mediates repression of chondrogenesis by β-catenin to promote cranial bone progenitor specification

The bones of the mammalian skull vault form through intramembranous ossification. Skull bones ossify directly, in a process regulated by β-catenin, instead of passing through a cartilage intermediate. We tested whether β-catenin is necessary for fate selection of intramembranous bone progenitors in the skull. Here, we show in mice that removal of β-catenin from skull bone progenitors results in the near complete transformation of the skull bones to cartilage, whereas constitutive β-catenin activation inhibits skull bone fate selection. β-catenin directly activated Twist1 expression in skull progenitors, conditional Twist1 deletion partially phenocopied the absence of β-catenin, and Twist1 deletion partially restored bone formation in the presence of constitutive β-catenin activation. Finally, Twist1 bound robustly to the 3'UTR of Sox9, the central initiator of chondrogenesis, suggesting that Twist1 might directly repress cartilage formation through Sox9. These findings provide insight into how β-catenin signaling via Twist1 actively suppresses the formation of cartilage and promotes intramembranous ossification in the skull.

Transcription factor redundancy and tissue-specific regulation: evidence from functional and physical network connectivity

Two major transcriptional regulators of Caenorhabditis elegans bodywall muscle (BWM) differentiation, hlh-1 and unc-120, are expressed in muscle where they are known to bind and regulate several well-studied muscle-specific genes. Simultaneously mutating both factors profoundly inhibits formation of contractile BWM. These observations were consistent with a simple network model in which the muscle regulatory factors drive tissue-specific transcription by binding selectively near muscle-specific targets to activate them. We tested this model by measuring the number, identity, and tissue-specificity of functional regulatory targets for each factor. Some joint regulatory targets (218) are BWM-specific and enriched for nearby HLH-1 binding. However, contrary to the simple model, the majority of genes regulated by one or both muscle factors are also expressed significantly in non-BWM tissues. We also mapped global factor occupancy by HLH-1, and created a genetic interaction map that identifies hlh-1 collaborating transcription factors. HLH-1 binding did not predict proximate regulatory action overall, despite enrichment for binding among BWM-specific positive regulatory targets of hlh-1. We conclude that these tissue-specific factors contribute much more broadly to the transcriptional output of muscle tissue than previously thought, offering a partial explanation for widespread HLH-1 occupancy. We also identify a novel regulatory connection between the BWM-specific hlh-1 network and the hlh-8/twist nonstriated muscle network. Finally, our results suggest a molecular basis for synthetic lethality in which hlh-1 and unc-120 mutant phenotypes are mutually buffered by joint additive regulation of essential target genes, with additional buffering suggested via newly identified hlh-1 interacting factors.

Characterization of enhancer function from genome-wide analyses

There has been a recent surge in the use of genome-wide methodologies to identify and annotate the transcriptional regulatory elements in the human genome. Here we review some of these methodologies and the conceptual insights about transcription regulation that have been gained from the use of genome-wide studies. It has become clear that the binding of transcription factors is itself a highly regulated process, and binding does not always appear to have functional consequences. Numerous properties have now been associated with regulatory elements that may be useful in their identification. Several aspects of enhancer function have been shown to be more widespread than was previously appreciated, including the highly combinatorial nature of transcription factor binding, the postinitiation regulation of many target genes, and the binding of enhancers at early stages to maintain their competence during development. Going forward, the integration of multiple genome-wide data sets should become a standard approach to elucidate higher-order regulatory interactions.

Lateral gene expression in Drosophila early embryos is supported by Grainyhead-mediated activation and tiers of dorsally-localized repression

The general consensus in the field is that limiting amounts of the transcription factor Dorsal establish dorsal boundaries of genes expressed along the dorsal-ventral (DV) axis of early Drosophila embryos, while repressors establish ventral boundaries. Yet recent studies have provided evidence that repressors act to specify the dorsal boundary of intermediate neuroblasts defective (ind), a gene expressed in a stripe along the DV axis in lateral regions of the embryo. Here we show that a short 12 base pair sequence ("the A-box") present twice within the ind CRM is both necessary and sufficient to support transcriptional repression in dorsal regions of embryos. To identify binding factors, we conducted affinity chromatography using the A-box element and found a number of DNA-binding proteins and chromatin-associated factors using mass spectroscopy. Only Grainyhead (Grh), a CP2 transcription factor with a unique DNA-binding domain, was found to bind the A-box sequence. Our results suggest that Grh acts as an activator to support expression of ind, which was surprising as we identified this factor using an element that mediates dorsally-localized repression. Grh and Dorsal both contribute to ind transcriptional activation. However, another recent study found that the repressor Capicua (Cic) also binds to the A-box sequence. While Cic was not identified through our A-box affinity chromatography, utilization of the same site, the A-box, by both factors Grh (activator) and Cic (repressor) may also support a "switch-like" response that helps to sharpen the ind dorsal boundary. Furthermore, our results also demonstrate that TGF-β signaling acts to refine ind CRM expression in an A-box independent manner in dorsal-most regions, suggesting that tiers of repression act in dorsal regions of the embryo.

Use of a Drosophila genome-wide conserved sequence database to identify functionally related cis-regulatory enhancers

Background: Phylogenetic footprinting has revealed that cis-regulatory enhancers consist of conserved DNA sequence clusters (CSCs). Currently, there is no systematic approach for enhancer discovery and analysis that takes full-advantage of the sequence information within enhancer CSCs. Results: We have generated a Drosophila genome-wide database of conserved DNA consisting of >100,000 CSCs derived from EvoPrints spanning over 90% of the genome. cis-Decoder database search and alignment algorithms enable the discovery of functionally related enhancers. The program first identifies conserved repeat elements within an input enhancer and then searches the database for CSCs that score highly against the input CSC. Scoring is based on shared repeats as well as uniquely shared matches, and includes measures of the balance of shared elements, a diagnostic that has proven to be useful in predicting cis-regulatory function. To demonstrate the utility of these tools, a temporally-restricted CNS neuroblast enhancer was used to identify other functionally related enhancers and analyze their structural organization. Conclusions:cis-Decoder reveals that co-regulating enhancers consist of combinations of overlapping shared sequence elements, providing insights into the mode of integration of multiple regulating transcription factors. The database and accompanying algorithms should prove useful in the discovery and analysis of enhancers involved in any developmental process. Developmental Dynamics 241:169–189, 2012. © 2011 Wiley Periodicals, Inc.

Complex interactions between cis-regulatory modules in native conformation are critical for Drosophila snail expression

It has been shown in several organisms that multiple cis-regulatory modules (CRMs) of a gene locus can be active concurrently to support similar spatiotemporal expression. To understand the functional importance of such seemingly redundant CRMs, we examined two CRMs from the Drosophila snail gene locus, which are both active in the ventral region of pre-gastrulation embryos. By performing a deletion series in a ∼25 kb DNA rescue construct using BAC recombineering and site-directed transgenesis, we demonstrate that the two CRMs are not redundant. The distal CRM is absolutely required for viability, whereas the proximal CRM is required only under extreme conditions such as high temperature. Consistent with their distinct requirements, the CRMs support distinct expression patterns: the proximal CRM exhibits an expanded expression domain relative to endogenous snail, whereas the distal CRM exhibits almost complete overlap with snail except at the anterior-most pole. We further show that the distal CRM normally limits the increased expression domain of the proximal CRM and that the proximal CRM serves as a `damper' for the expression levels driven by the distal CRM. Thus, the two CRMs interact in cis in a non-additive fashion and these interactions may be important for fine-tuning the domains and levels of gene expression.