Transcriptional enhancers in animal development and evolution

The chromatin "landscape" of a murine adult β-globin gene is unaffected by deletion of either the gene promoter or a downstream enhancer

In mammals, the complex tissue- and developmental-specific expression of genes within the β-globin cluster is known to be subject to control by the gene promoters, by a locus control region (LCR) located upstream of the cluster, and by sequence elements located across the intergenic regions. Despite extensive investigation, however, the complement of sequences that is required for normal regulation of chromatin structure and gene expression within the cluster is not fully defined. To further elucidate regulation of the adult β-globin genes, we investigate the effects of two deletions engineered within the endogenous murine β-globin locus. First, we find that deletion of the β2-globin gene promoter, while eliminating β2-globin gene expression, results in no additional effects on chromatin structure or gene expression within the cluster. Notably, our observations are not consistent with competition among the β-globin genes for LCR activity. Second, we characterize a novel enhancer located 3′ of the β2-globin gene, but find that deletion of this sequence has no effect whatsoever on gene expression or chromatin structure. This observation highlights the difficulty in assigning function to enhancer sequences identified by the chromatin “landscape” or even by functional assays.

Transcription factors: specific DNA binding and specific gene regulation

Specific recognition of cis-regulatory regions is essential for correct gene regulation in response to developmental and environmental signals. Such DNA sequences are recognized by transcription factors (TFs) that recruit the transcriptional machinery. Achievement of specific sequence recognition is not a trivial problem; many TFs recognize similar consensus DNA-binding sites and a genome can harbor thousands of consensus or near-consensus sequences, both functional and nonfunctional. Although genomic technologies have provided large-scale snapshots of TF binding, a full understanding of the mechanistic and quantitative details of specific recognition in the context of gene regulation is lacking. Here, we explore the various ways in which TFs recognizing similar consensus sites distinguish their own targets from a large number of other sequences to ensure specific cellular responses.

The role of self-organization in developmental evolution

In developmental and evolutionary biology, particular emphasis has been given to the relationship between transcription factors and the cognate cis-regulatory elements of their target genes. These constitute the gene regulatory networks that control expression and are assumed to causally determine the formation of structures and body plans. Comparative analysis has, however, established a broad sequence homology among species that nonetheless display quite different anatomies. Transgenic experiments have also confirmed that many developmentally important elements are, in fact, functionally interchangeable. Although dependent upon the appropriate degree of gene expression, the actual construction of specific structures appears not directly linked to the functions of gene products alone. Instead, the self-formation of complex patterns, due in large part to epigenetic and non-genetic determinants, remains a persisting theme in the study of ontogeny and regenerative medicine. Recent evidence indeed points to the existence of a self-organizing process, operating through a set of intrinsic rules and forces, which imposes coordination and a holistic order upon cells and tissue. This has been repeatedly demonstrated in experiments on regeneration as well as in the autonomous formation of structures in vitro. The process cannot be wholly attributed to the functional outcome of protein-protein interactions or to concentration gradients of diffusible chemicals. This phenomenon is examined here along with some of the methodological and theoretical approaches that are now used in understanding the causal basis for self-organization in development and its evolution.

Dissection of thousands of cell type-specific enhancers identifies dinucleotide repeat motifs as general enhancer features

Gene expression is determined by genomic elements called enhancers, which contain short motifs bound by different transcription factors (TFs). However, how enhancer sequences and TF motifs relate to enhancer activity is unknown, and general sequence requirements for enhancers or comprehensive sets of important enhancer sequence elements have remained elusive. Here, we computationally dissect thousands of functional enhancer sequences from three different Drosophila cell lines. We find that the enhancers display distinct cis-regulatory sequence signatures, which are predictive of the enhancers’ cell type-specific or broad activities. These signatures contain transcription factor motifs and a novel class of enhancer sequence elements, dinucleotide repeat motifs (DRMs). DRMs are highly enriched in enhancers, particularly in enhancers that are broadly active across different cell types. We experimentally validate the importance of the identified TF motifs and DRMs for enhancer function and show that they can be sufficient to create an active enhancer de novo from a nonfunctional sequence. The function of DRMs as a novel class of general enhancer features that are also enriched in human regulatory regions might explain their implication in several diseases and provides important insights into gene regulation.

Nucleosomes are stably evicted from enhancers but not promoters upon induction of certain pro-inflammatory genes in mouse macrophages

Chromatin is thought to act as a barrier for binding of cis-regulatory transcription factors (TFs) to their sites on DNA and recruitment of the transcriptional machinery. Here we have analyzed changes in nucleosome occupancy at the enhancers as well as at the promoters of three pro-inflammatory genes when they are induced by bacterial lipopolysaccharides (LPS) in primary mouse macrophages. We find that nucleosomes are removed from the distal enhancers of IL12B and IL1A, as well as from the distal and proximal enhancers of IFNB1, and that clearance of enhancers correlates with binding of various cis-regulatory TFs. We further show that for IFNB1 the degree of nucleosome removal correlates well with the level of induction of the gene under different conditions. Surprisingly, we find that nucleosome occupancy at the promoters of IL12B and IL1A does not change significantly when the genes are induced, and that a considerably fraction of the cells is occupied by nucleosomes at any given time. We hypothesize that competing nucleosomes at the promoters of IL12B and IL1A may play a role in limiting the size of transcriptional bursts in individual cells, which may be important for controlling cytokine production in a population of immune cells.

Hormone-responsive enhancer-activity maps reveal predictive motifs, indirect repression, and targeting of closed chromatin

Steroid hormones act as important developmental switches, and their nuclear receptors regulate many genes. However, few hormone-dependent enhancers have been characterized, and important aspects of their sequence architecture, cell-type-specific activating and repressing functions, or the regulatory roles of their chromatin structure have remained unclear. We used STARR-seq, a recently developed enhancer-screening assay, and ecdysone signaling in two different Drosophila cell types to derive genome-wide hormone-dependent enhancer-activity maps. We demonstrate that enhancer activation depends on cis-regulatory motif combinations that differ between cell types and can predict cell-type-specific ecdysone targeting. Activated enhancers are often not accessible prior to induction. Enhancer repression following hormone treatment seems independent of receptor motifs and receptor binding to the enhancer, as we show using ChIP-seq, but appears to rely on motifs for other factors, including Eip74. Our strategy is applicable to study signal-dependent enhancers for different pathways and across organisms.

Imogene: identification of motifs and cis-regulatory modules underlying gene co-regulation

Cis-regulatory modules (CRMs) and motifs play a central role in tissue and condition-specific gene expression. Here we present Imogene, an ensemble of statistical tools that we have developed to facilitate their identification and implemented in a publicly available software. Starting from a small training set of mammalian or fly CRMs that drive similar gene expression profiles, Imogene determines de novocis-regulatory motifs that underlie this co-expression. It can then predict on a genome-wide scale other CRMs with a regulatory potential similar to the training set. Imogene bypasses the need of large datasets for statistical analyses by making central use of the information provided by the sequenced genomes of multiple species, based on the developed statistical tools and explicit models for transcription factor binding site evolution. We test Imogene on characterized tissue-specific mouse developmental CRMs. Its ability to identify CRMs with the same specificity based on its de novo created motifs is comparable to that of previously evaluated ‘motif-blind’ methods. We further show, both in flies and in mammals, that Imogene de novo generated motifs are sufficient to discriminate CRMs related to different developmental programs. Notably, purely relying on sequence data, Imogene performs as well in this discrimination task as a previously reported learning algorithm based on Chromatin Immunoprecipitation (ChIP) data for multiple transcription factors at multiple developmental stages.

Evolutionary conservation of the eumetazoan gene regulatory landscape

Despite considerable differences in morphology and complexity of body plans among animals, a great part of the gene set is shared among Bilateria and their basally branching sister group, the Cnidaria. This suggests that the common ancestor of eumetazoans already had a highly complex gene repertoire. At present it is therefore unclear how morphological diversification is encoded in the genome. Here we address the possibility that differences in gene regulation could contribute to the large morphological divergence between cnidarians and bilaterians. To this end, we generated the first genome-wide map of gene regulatory elements in a nonbilaterian animal, the sea anemone Nematostella vectensis. Using chromatin immunoprecipitation followed by deep sequencing of five chromatin modifications and a transcriptional cofactor, we identified over 5000 enhancers in the Nematostella genome and could validate 75% of the tested enhancers in vivo. We found that in Nematostella, but not in yeast, enhancers are characterized by the same combination of histone modifications as in bilaterians, and these enhancers preferentially target developmental regulatory genes. Surprisingly, the distribution and abundance of gene regulatory elements relative to these genes are shared between Nematostella and bilaterian model organisms. Our results suggest that complex gene regulation originated at least 600 million yr ago, predating the common ancestor of eumetazoans.

Rapid and pervasive changes in genome-wide enhancer usage during mammalian development

Enhancers are distal regulatory elements that can activate tissue-specific gene expression and are abundant throughout mammalian genomes. Although substantial progress has been made toward genome-wide annotation of mammalian enhancers, their temporal activity patterns and global contributions in the context of developmental in vivo processes remain poorly explored. Here we used epigenomic profiling for H3K27ac, a mark of active enhancers, coupled to transgenic mouse assays to examine the genome-wide utilization of enhancers in three different mouse tissues across seven developmental stages. The majority of the ∼90,000 enhancers identified exhibited tightly temporally restricted predicted activity windows and were associated with stage-specific biological functions and regulatory pathways in individual tissues. Comparative genomic analysis revealed that evolutionary conservation of enhancers decreases following midgestation across all tissues examined. The dynamic enhancer activities uncovered in this study illuminate rapid and pervasive temporal in vivo changes in enhancer usage that underlie processes central to development and disease.

Cooperative activation of Xenopus rhodopsin transcription by paired-like transcription factors

Background

In vertebrates, rod photoreceptor-specific gene expression is regulated by the large Maf and Pax-like transcription factors, Nrl/LNrl and Crx/Otx5. The ubiquitous occurrence of their target DNA binding sites throughout rod-specific gene promoters suggests that multiple transcription factor interactions within the promoter are functionally important. Cooperative action by these transcription factors activates rod-specific genes such as rhodopsin. However, a quantitative mechanistic explanation of transcriptional rate determinants is lacking.

Results

We investigated the contributions of various paired-like transcription factors and their cognate cis-elements to rhodopsin gene activation using cultured cells to quantify activity. The Xenopus rhodopsin promoter (XOP) has a bipartite structure, with ~200 bp proximal to the start site (RPP) coordinating cooperative activation by Nrl/LNrl-Crx/Otx5 and the adjacent 5300 bp upstream sequence increasing the overall expression level. The synergistic activation by Nrl/LNrl-Crx/Otx5 also occurred when XOP was stably integrated into the genome. We determined that Crx/Otx5 synergistically activated transcription independently and additively through the two Pax-like cis-elements, BAT1 and Ret4, but not through Ret1. Other Pax-like family members, Rax1 and Rax2, do not synergistically activate XOP transcription with Nrl/LNrl and/or Crx/Otx5; rather they act as co-activators via the Ret1 cis-element.

Conclusions

We have provided a quantitative model of cooperative transcriptional activation of the rhodopsin promoter through interaction of Crx/Otx5 with Nrl/LNrl at two paired-like cis-elements proximal to the NRE and TATA binding site. Further, we have shown that Rax genes act in cooperation with Crx/Otx5 with Nrl/LNrl as co-activators of rhodopsin transcription.

A comparison of midline and tracheal gene regulation during Drosophila development

Within the Drosophila embryo, two related bHLH-PAS proteins, Single-minded and Trachealess, control development of the central nervous system midline and the trachea, respectively. These two proteins are bHLH-PAS transcription factors and independently form heterodimers with another bHLH-PAS protein, Tango. During early embryogenesis, expression of Single-minded is restricted to the midline and Trachealess to the trachea and salivary glands, whereas Tango is ubiquitously expressed. Both Single-minded/Tango and Trachealess/Tango heterodimers bind to the same DNA sequence, called the CNS midline element (CME) within cis-regulatory sequences of downstream target genes. While Single-minded/Tango and Trachealess/Tango activate some of the same genes in their respective tissues during embryogenesis, they also activate a number of different genes restricted to only certain tissues. The goal of this research is to understand how these two related heterodimers bind different enhancers to activate different genes, thereby regulating the development of functionally diverse tissues. Existing data indicates that Single-minded and Trachealess may bind to different co-factors restricted to various tissues, causing them to interact with the CME only within certain sequence contexts. This would lead to the activation of different target genes in different cell types. To understand how the context surrounding the CME is recognized by different bHLH-PAS heterodimers and their co-factors, we identified and analyzed novel enhancers that drive midline and/or tracheal expression and compared them to previously characterized enhancers. In addition, we tested expression of synthetic reporter genes containing the CME flanked by different sequences. Taken together, these experiments identify elements overrepresented within midline and tracheal enhancers and suggest that sequences immediately surrounding a CME help dictate whether a gene is expressed in the midline or trachea.

Assessing constraints on the path of regulatory sequence evolution

Structural and functional constraints are known to play a major role in restricting the path of evolution of protein activities. However, constraints acting on evolving transcriptional regulatory sequences, e.g. enhancers, are largely unknown. Recently, we elucidated how a novel expression pattern of the Neprilysin-1 (Nep1) gene in the optic lobe of Drosophila santomea evolved via co-option of existing enhancer activities. Drosophila santomea, which has diverged from Drosophila yakuba by approximately 400 000 years has accumulated four fixed mutations that each contribute to the full activity of this enhancer. Recreating and testing the optic lobe enhancer of the ancestor of D. santomea and D. yakuba revealed that the strong D. santomea enhancer activity evolved from a weak ancestral activity. Because each mutation on the path from the D. yakuba/santomea ancestor to modern-day D. santomea contributes to the newly derived optic lobe enhancer activity, we sought here to use this system to study the path of evolution of enhancer sequences. We inferred likely paths of evolution of this enhancer by observing the transcriptional output of all possible intermediate steps between the ancestral D. yakuba/santomea enhancer and the modern D. santomea enhancer. Many possible paths had epistatic and cooperative effects. Furthermore, we found that several paths significantly increased ectopic transcriptional activity or affected existing enhancer activities from which the novel activity was co-opted. We suggest that these attributes highlight constraints that guide the path of evolution of enhancers.

Topology of mammalian developmental enhancers and their regulatory landscapes

How a complex animal can arise from a fertilized egg is one of the oldest and most fascinating questions of biology, the answer to which is encoded in the genome. Body shape and organ development, and their integration into a functional organism all depend on the precise expression of genes in space and time. The orchestration of transcription relies mostly on surrounding control sequences such as enhancers, millions of which form complex regulatory landscapes in the non-coding genome. Recent research shows that high-order chromosome structures make an important contribution to enhancer functionality by triggering their physical interactions with target genes.

Expression-based segmentation of the Drosophila genome

Background

It is generally accepted that gene order in eukaryotes is nonrandom, with adjacent genes often sharing expression patterns across tissues, and that this organization may be important for gene regulation. Here we describe a novel method, based on an explicit probability model instead of correlation analysis, for identifying coordinately expressed gene clusters (‘coexpression segments’), apply it to Drosophila melanogaster, and look for epigenetic associations using publicly available data.

Results

We find that two-thirds of Drosophila genes fall into multigenic coexpression segments, and that such segments are of two main types, housekeeping and tissue-restricted. Consistent with correlation-based studies, we find that adjacent genes within the same segment tend to be physically closer to each other than to the adjacent genes in different segments, and that tissue-restricted segments are enriched for testis-expressed genes. Our segmentation pattern correlates with Hi-C based physical interaction domains, but segments are generally much smaller than domains. Intersegment regions (including those which do not correspond to physical domain boundaries) are enriched for insulator binding sites.

Conclusions

We describe a novel approach for identifying coexpression clusters that does not require arbitrary cutoff values or heuristics, and find that coexpression of adjacent genes is widespread in the Drosophila genome. Coexpression segments appear to reflect a level of regulatory organization related to, but below that of physical interaction domains, and depending in part on insulator binding.

Low-affinity transcription factor binding sites shape morphogen responses and enhancer evolution

In the era of functional genomics, the role of transcription factor (TF)-DNA binding affinity is of increasing interest: for example, it has recently been proposed that low-affinity genomic binding events, though frequent, are functionally irrelevant. Here, we investigate the role of binding site affinity in the transcriptional interpretation of Hedgehog (Hh) morphogen gradients. We noted that enhancers of several Hh-responsive Drosophila genes have low predicted affinity for Ci, the Gli family TF that transduces Hh signalling in the fly. Contrary to our initial hypothesis, improving the affinity of Ci/Gli sites in enhancers of dpp, wingless and stripe, by transplanting optimal sites from the patched gene, did not result in ectopic responses to Hh signalling. Instead, we found that these enhancers require low-affinity binding sites for normal activation in regions of relatively low signalling. When Ci/Gli sites in these enhancers were altered to improve their binding affinity, we observed patterning defects in the transcriptional response that are consistent with a switch from Ci-mediated activation to Ci-mediated repression. Synthetic transgenic reporters containing isolated Ci/Gli sites confirmed this finding in imaginal discs. We propose that the requirement for gene activation by Ci in the regions of low-to-moderate Hh signalling results in evolutionary pressure favouring weak binding sites in enhancers of certain Hh target genes.

Discover context-specific combinatorial transcription factor interactions by integrating diverse ChIP-Seq data sets

Combinatorial interactions among transcription factors (TFs) are critical for integrating diverse intrinsic and extrinsic signals, fine-tuning regulatory output and increasing the robustness and plasticity of regulatory systems. Current knowledge about combinatorial regulation is rather limited due to the lack of suitable experimental technologies and bioinformatics tools. The rapid accumulation of ChIP-Seq data has provided genome-wide occupancy maps for a large number of TFs and chromatin modification marks for identifying enhancers without knowing individual TF binding sites. Integration of the two data types has not been researched extensively, resulting in underused data and missed opportunities. We describe a novel method for discovering frequent combinatorial occupancy patterns by multiple TFs at enhancers. Our method is based on probabilistic item set mining and takes into account uncertainty in both types of ChIP-Seq data. By joint analysis of 108 TFs in four human cell types, we found that cell–type-specific interactions among TFs are abundant and that the majority of enhancers have flexible architecture. We show that several families of transposable elements disproportionally overlap with enhancers with combinatorial patterns, suggesting that these transposable element families play an important role in the evolution of combinatorial regulation.

Development: lights, camera, action--the Drosophila embryo goes live!

Live imaging of developmental gene expression in Drosophila embryos opens up exciting new prospects for understanding gene regulation during development.

Confirmation of genetic homogeneity of syndactyly type IV and triphalangeal thumb-polysyndactyly syndrome in a Chinese family and review of the literature

Syndactyly type IV (SD4) is inherited in an autosomal dominant fashion and characterized by complete cutaneous syndactyly of all fingers accompanied with polydactyly. Triphalangeal thumb-polysyndactyly syndrome (TPTPS) consists of a triphalangeal thumb, polydactyly, and syndactyly and is transmitted in an autosomal dominant manner with variable expression. Genomic duplications of the long-range limb-specific cis-regulator (ZRS) cause SD4 and TPTPS. Here, we report two individuals from a Chinese family with syndactyly. One individual had overlapping clinical symptoms of TPTPS and SD4, while the other had a typical SD4 with postaxial polydactyly of the toe. Results of quantitative PCR suggested that the duplication of ZRS involved all affected individuals, and array comparative genomic hybridization detected its size as 115.3 kb.

CONCLUSION

This work confirms the genetic homogeneity of SD4 and TPTPS. Our result expands the spectrum of ZRS duplications. TPTPS and SD4 should be considered as a continuum of phenotypes.

Dynamic model for the coordination of two enhancers of broad by EGFR signaling

Although it is widely appreciated that a typical developmental control gene is regulated by multiple enhancers, coordination of enhancer activities remains poorly understood. We propose a mechanism for such coordination in Drosophila oogenesis, when the expression of the transcription factor Broad (BR) evolves from a uniform to a two-domain pattern that prefigures the formation of two respiratory eggshell appendages. This change reflects sequential activities of two enhancers of the br gene, early and late, both of which are controlled by the epidermal growth factor receptor (EGFR) pathway. The late enhancer controls br in the appendage-producing cells, but the function of the early enhancer remained unclear. We found that the early enhancer is essential for the activity of the late enhancer and induction of eggshell appendages. This requirement can be explained by a mechanism whereby the BR protein produced by the early enhancer protects the late enhancer from EGFR-dependent repression. We illustrate this complex mechanism using a computational model that correctly predicts the wild-type dynamics of BR expression and its response to genetic perturbations.

Distinct regulatory mechanisms act to establish and maintain Pax3 expression in the developing neural tube

Pattern formation in developing tissues is driven by the interaction of extrinsic signals with intrinsic transcriptional networks that together establish spatially and temporally restricted profiles of gene expression. How this process is orchestrated at the molecular level by genomic cis-regulatory modules is one of the central questions in developmental biology. Here we have addressed this by analysing the regulation of Pax3 expression in the context of the developing spinal cord. Pax3 is induced early during neural development in progenitors of the dorsal spinal cord and is maintained as pattern is subsequently elaborated, resulting in the segregation of the tissue into dorsal and ventral subdivisions. We used a combination of comparative genomics and transgenic assays to define and dissect several functional cis-regulatory modules associated with the Pax3 locus. We provide evidence that the coordinated activity of two modules establishes and refines Pax3 expression during neural tube development. Mutational analyses of the initiating element revealed that in addition to Wnt signaling, Nkx family homeodomain repressors restrict Pax3 transcription to the presumptive dorsal neural tube. Subsequently, a second module mediates direct positive autoregulation and feedback to maintain Pax3 expression. Together, these data indicate a mechanism by which transient external signals are converted into a sustained expression domain by the activities of distinct regulatory elements. This transcriptional logic differs from the cross-repression that is responsible for the spatiotemporal patterns of gene expression in the ventral neural tube, suggesting that a variety of circuits are deployed within the neural tube regulatory network to establish and elaborate pattern formation.

The complex organization of tissues is established precisely and reproducibly during development. In the vertebrate neural tube, as in many other tissues, the interplay between extrinsic morphogens and intrinsic transcription factors produces spatial patterns of gene expression that delineate precursors for specific cell types. One such transcription factor, Pax3, defines the precursors of all sensory neuron subtypes and distinguishes them from precursors fated to give rise to the motor circuits. To gain insight into the molecular mechanisms by which the spinal cord is segregated into these two functional domains, we analysed the genomic regulatory sequences responsible for controlling Pax3 activity. We identified two regions of the genome, the coordinated activity of which establishes and refines Pax3 activity. We showed that the combination of activating signals from secreted Wnt factors together with Nkx family homeodomain repressors restrict Pax3 activity to the presumptive sensory region of the neural tissue. Subsequently, Pax3 acts to directly potentiate its own transcription and this autoregulation sustains Pax3 expression at later developmental stages. Together, our study reveals the way in which intrinsic and extrinsic signals are integrated by cells and converted into a sustained pattern of gene activity in the developing nervous system.

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.

Transcription factor binding in human cells occurs in dense clusters formed around cohesin anchor sites

During cell division, transcription factors (TFs) are removed from chromatin twice, during DNA synthesis and during condensation of chromosomes. How TFs can efficiently find their sites following these stages has been unclear. Here, we have analyzed the binding pattern of expressed TFs in human colorectal cancer cells. We find that binding of TFs is highly clustered and that the clusters are enriched in binding motifs for several major TF classes. Strikingly, almost all clusters are formed around cohesin, and loss of cohesin decreases both DNA accessibility and binding of TFs to clusters. We show that cohesin remains bound in S phase, holding the nascent sister chromatids together at the TF cluster sites. Furthermore, cohesin remains bound to the cluster sites when TFs are evicted in early M phase. These results suggest that cohesin-binding functions as a cellular memory that promotes re-establishment of TF clusters after DNA replication and chromatin condensation.

Cellular resolution models for even skipped regulation in the entire Drosophila embryo

Transcriptional control ensures genes are expressed in the right amounts at the correct times and locations. Understanding quantitatively how regulatory systems convert input signals to appropriate outputs remains a challenge. For the first time, we successfully model even skipped (eve) stripes 2 and 3+7 across the entire fly embryo at cellular resolution. A straightforward statistical relationship explains how transcription factor (TF) concentrations define eve's complex spatial expression, without the need for pairwise interactions or cross-regulatory dynamics. Simulating thousands of TF combinations, we recover known regulators and suggest new candidates. Finally, we accurately predict the intricate effects of perturbations including TF mutations and misexpression. Our approach imposes minimal assumptions about regulatory function; instead we infer underlying mechanisms from models that best fit the data, like the lack of TF-specific thresholds and the positional value of homotypic interactions. Our study provides a general and quantitative method for elucidating the regulation of diverse biological systems. DOI:http://dx.doi.org/10.7554/eLife.00522.001.

Reprogramming to pluripotency using designer TALE transcription factors targeting enhancers

The modular DNA recognition code of the transcription-activator-like effectors (TALEs) from plant pathogenic bacterial genus Xanthomonas provides a powerful genetic tool to create designer transcription factors (dTFs) targeting specific DNA sequences for manipulating gene expression. Previous studies have suggested critical roles of enhancers in gene regulation and reprogramming. Here, we report dTF activator targeting the distal enhancer of the Pou5f1 (Oct4) locus induces epigenetic changes, reactivates its expression, and substitutes exogenous OCT4 in reprogramming mouse embryonic fibroblast cells (MEFs) to induced pluripotent stem cells (iPSCs). Similarly, dTF activator targeting a Nanog enhancer activates Nanog expression and reprograms epiblast stem cells (EpiSCs) to iPSCs. Conversely, dTF repressors targeting the same genetic elements inhibit expression of these loci, and effectively block reprogramming. This study indicates that dTFs targeting specific enhancers can be used to study other biological processes such as transdifferentiation or directed differentiation of stem cells.

•

The dTF for the Oct4 distal enhancer replaces OCT4 in reprogramming MEFs to iPSCs

•

The repressor dTF induces ESC differentiation and blocks reprogramming

•

The dTF for the Nanog reprograms EpiSCs to iPSCs, whereas the repressor blocks it

•

Targeting enhancers by dTFs can be an effective way to control gene expression

Drosophila distal-less and Rotund bind a single enhancer ensuring reliable and robust bric-a-brac2 expression in distinct limb morphogenetic fields

Most identified Drosophila appendage-patterning genes encode DNA-binding proteins, whose cross-regulatory interactions remain to be better characterized at the molecular level, notably by studying their direct binding to tissue-specific transcriptional enhancers. A fine-tuned spatio-temporal expression of bric-a-brac2 (bab2) along concentric rings is essential for proper proximo-distal (P-D) differentiation of legs and antennae. However, within the genetic interaction landscape governing limb development, no transcription factor directly controlling bab2 expression has been identified to date. Using site-targeted GFP reporter assay and BAC recombineering, we show here that restricted bab2 expression in leg and antennal imaginal discs relies on a single 567-bp-long cis-regulatory module (CRM), termed LAE (for leg and antennal enhancer). We show that this CRM (i) is necessary and sufficient to ensure normal bab2 activity in developing leg and antenna, and (ii) is structurally and functionally conserved among Drosophilidae. Through deletion and site-directed mutagenesis approaches, we identified within the LAE essential sequence motifs required in both leg and antennal tissues. Using genetic and biochemical tests, we establish that in the LAE (i) a key TAAT-rich activator motif interacts with the homeodomain P-D protein Distal-less (Dll) and (ii) a single T-rich activator motif binds the C2H2 zinc-finger P-D protein Rotund (Rn), leading to bab2 up-regulation respectively in all or specifically in the proximal-most ring(s), both in leg and antenna. Joint ectopic expression of Dll and Rn is sufficient to cell-autonomously activate endogenous bab2 and LAE-driven reporter expression in wing and haltere cells. Our findings indicate that accuracy, reliability and robustness of developmental gene expression do not necessarily require cis-regulatory information redundancy.

In insects, leg and antenna are homologous limbs, though derive from a single ancestral appendage. In Drosophila, leg and antennal development along the proximo-distal (P-D) axis relies on relatively-well known genetic cascades, in which most appendage-patterning genes encode transcription factors (TF). However, their cross-regulatory interactions remain to be better characterized at the molecular level. A fine-tuned expression of the bric-a-brac2 (bab2) gene is essential for normal leg and antennal segmentation. However, within the genetic cascades governing P-D limb development, no TF directly controlling bab2 expression has been identified to date. We show here that restricted bab2 expression in developing leg and antenna is governed by a single enhancer, termed LAE, which is necessary and sufficient in-vivo to ensure bab2 functions there. We show that leg and antennal cis-regulatory elements are closely associated and that essential LAE sites interact with Distal-less (Dll) and Rotund (Rn) TFs, leading to bab2 activation in all or specifically in the proximal-most expressing cells, respectively. Finally, joint ectopic expression of Dll and Rn is sufficient to instruct wing and haltere cells to up-regulate bab2. Taken together, our work indicates that a single enhancer is necessary and sufficient to reliably govern bab2 expression in distinct morphogenetic fields.

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.

Enhancers: five essential questions

It is estimated that the human genome contains hundreds of thousands of enhancers, so understanding these gene-regulatory elements is a crucial goal. Several fundamental questions need to be addressed about enhancers, such as how do we identify them all, how do they work, and how do they contribute to disease and evolution? Five prominent researchers in this field look at how much we know already and what needs to be done to answer these questions.

Human genetic variation within neural crest enhancers: molecular and phenotypic implications

Developmental gene expression programmes are coordinated by the specialized distal cis-regulatory elements called enhancers, which integrate lineage- and signalling-dependent inputs to guide morphogenesis. In previous work, we characterized the genome-wide repertoire of active enhancers in human neural crest cells (hNCC), an embryonic cell population with critical roles in craniofacial development. We showed that in hNCC, co-occupancy of a master regulator TFAP2A with nuclear receptors NR2F1 and NR2F2 correlates with the presence of permissive enhancer chromatin states. Here, we take advantage of pre-existing human genetic variation to further explore potential cooperation between TFAP2A and NR2F1/F2. We demonstrate that isolated single nucleotide polymorphisms affecting NR2F1/F2-binding sites within hNCC enhancers can alter TFAP2A occupancy and overall chromatin features at the same enhancer allele. We propose that a similar strategy can be used to elucidate other cooperative relationships between transcription factors involved in developmental transitions. Using the neural crest and its major contribution to human craniofacial phenotypes as a paradigm, we discuss how genetic variation might modulate the molecular properties and activity of enhancers, and ultimately impact human phenotypic diversity.

High-throughput identification of long-range regulatory elements and their target promoters in the human genome

Enhancer elements are essential for tissue-specific gene regulation during mammalian development. Although these regulatory elements are often distant from their target genes, they affect gene expression by recruiting transcription factors to specific promoter regions. Because of this long-range action, the annotation of enhancer element-target promoter pairs remains elusive. Here, we developed a novel analysis methodology that takes advantage of Hi-C data to comprehensively identify these interactions throughout the human genome. To do this, we used a geometric distribution-based model to identify DNA-DNA interaction hotspots that contact gene promoters with high confidence. We observed that these promoter-interacting hotspots significantly overlap with known enhancer-associated histone modifications and DNase I hypersensitive sites. Thus, we defined thousands of candidate enhancer elements by incorporating these features, and found that they have a significant propensity to be bound by p300, an enhancer binding transcription factor. Furthermore, we revealed that their target genes are significantly bound by RNA Polymerase II and demonstrate tissue-specific expression. Finally, we uncovered that these elements are generally found within 1 Mb of their targets, and often regulate multiple genes. In total, our study presents a novel high-throughput workflow for confident, genome-wide discovery of enhancer-target promoter pairs, which will significantly improve our understanding of these regulatory interactions.

RFECS: a random-forest based algorithm for enhancer identification from chromatin state

Transcriptional enhancers play critical roles in regulation of gene expression, but their identification in the eukaryotic genome has been challenging. Recently, it was shown that enhancers in the mammalian genome are associated with characteristic histone modification patterns, which have been increasingly exploited for enhancer identification. However, only a limited number of cell types or chromatin marks have previously been investigated for this purpose, leaving the question unanswered whether there exists an optimal set of histone modifications for enhancer prediction in different cell types. Here, we address this issue by exploring genome-wide profiles of 24 histone modifications in two distinct human cell types, embryonic stem cells and lung fibroblasts. We developed a Random-Forest based algorithm, RFECS (Random Forest based Enhancer identification from Chromatin States) to integrate histone modification profiles for identification of enhancers, and used it to identify enhancers in a number of cell-types. We show that RFECS not only leads to more accurate and precise prediction of enhancers than previous methods, but also helps identify the most informative and robust set of three chromatin marks for enhancer prediction.

Enhancers are regions in the genome that can activate the expression of a gene irrespective of their location with respect to the gene. Identifying these elements is critical in understanding regulatory differences between different cell-types. Since enhancers lack characteristic sequence features and can be far away from the gene they regulate, their identification is not trivial. Experimentally determining the genome-wide binding sites of transcriptional co-activator p300 is one way of finding enhancers but it can only identify a subset of enhancers. A few years ago, it was observed that the binding sites of p300 are marked by distinctive, post-translational histone modifications. Several groups have exploited this discovery to predict genome-wide enhancers based on their similarity to the histone modification profiles of p300 binding sites. We here report a novel algorithm for this purpose and show that it has much greater accuracy than existing methods. Another unique feature of our algorithm is the ability to automatically deduce the most informative set of histone modifications required for enhancer prediction. We expect that this method will become increasingly useful with the expanding number of known histone modifications and rapid accumulation of epigenomic datasets for various cell types and species.

PRISM offers a comprehensive genomic approach to transcription factor function prediction

The human genome encodes 1500–2000 different transcription factors (TFs). ChIP-seq is revealing the global binding profiles of a fraction of TFs in a fraction of their biological contexts. These data show that the majority of TFs bind directly next to a large number of context-relevant target genes, that most binding is distal, and that binding is context specific. Because of the effort and cost involved, ChIP-seq is seldom used in search of novel TF function. Such exploration is instead done using expression perturbation and genetic screens. Here we propose a comprehensive computational framework for transcription factor function prediction. We curate 332 high-quality nonredundant TF binding motifs that represent all major DNA binding domains, and improve cross-species conserved binding site prediction to obtain 3.3 million conserved, mostly distal, binding site predictions. We combine these with 2.4 million facts about all human and mouse gene functions, in a novel statistical framework, in search of enrichments of particular motifs next to groups of target genes of particular functions. Rigorous parameter tuning and a harsh null are used to minimize false positives. Our novel PRISM (predicting regulatory information from single motifs) approach obtains 2543 TF function predictions in a large variety of contexts, at a false discovery rate of 16%. The predictions are highly enriched for validated TF roles, and 45 of 67 (67%) tested binding site regions in five different contexts act as enhancers in functionally matched cells.

Genome-wide quantitative enhancer activity maps identified by STARR-seq

Genomic enhancers are important regulators of gene expression, but their identification is a challenge, and methods depend on indirect measures of activity. We developed a method termed STARR-seq to directly and quantitatively assess enhancer activity for millions of candidates from arbitrary sources of DNA, which enables screens across entire genomes. When applied to the Drosophila genome, STARR-seq identifies thousands of cell type-specific enhancers across a broad continuum of strengths, links differential gene expression to differences in enhancer activity, and creates a genome-wide quantitative enhancer map. This map reveals the highly complex regulation of transcription, with several independent enhancers for both developmental regulators and ubiquitously expressed genes. STARR-seq can be used to identify and quantify enhancer activity in other eukaryotes, including humans.

Enhancer-associated H3K4 monomethylation by Trithorax-related, the Drosophila homolog of mammalian Mll3/Mll4

Monomethylation of histone H3 on Lys 4 (H3K4me1) and acetylation of histone H3 on Lys 27 (H3K27ac) are histone modifications that are highly enriched over the body of actively transcribed genes and on enhancers. Although in yeast all H3K4 methylation patterns, including H3K4me1, are implemented by Set1/COMPASS (complex of proteins associated with Set1), there are three classes of COMPASS-like complexes in Drosophila that could carry out H3K4me1 on enhancers: dSet1, Trithorax, and Trithorax-related (Trr). Here, we report that Trr, the Drosophila homolog of the mammalian Mll3/4 COMPASS-like complexes, can function as a major H3K4 monomethyltransferase on enhancers in vivo. Loss of Trr results in a global decrease of H3K4me1 and H3K27ac levels in various tissues. Assays with the cut wing margin enhancer implied a functional role for Trr in enhancer-mediated processes. A genome-wide analysis demonstrated that Trr is required to maintain the H3K4me1 and H3K27ac chromatin signature that resembles the histone modification patterns described for enhancers. Furthermore, studies in the mammalian system suggested a role for the Trr homolog Mll3 in similar processes. Since Trr and mammalian Mll3/4 complexes are distinguished by bearing a unique subunit, the H3K27 demethylase UTX, we propose a model in which the H3K4 monomethyltransferases Trr/Mll3/Mll4 and the H3K27 demethylase UTX cooperate to regulate the transition from inactive/poised to active enhancers.

A dynamic and intricate regulatory network determines Pseudomonas aeruginosa virulence

Pseudomonas aeruginosa is a metabolically versatile bacterium that is found in a wide range of biotic and abiotic habitats. It is a major human opportunistic pathogen causing numerous acute and chronic infections. The critical traits contributing to the pathogenic potential of P. aeruginosa are the production of a myriad of virulence factors, formation of biofilms and antibiotic resistance. Expression of these traits is under stringent regulation, and it responds to largely unidentified environmental signals. This review is focused on providing a global picture of virulence gene regulation in P. aeruginosa. In addition to key regulatory pathways that control the transition from acute to chronic infection phenotypes, some regulators have been identified that modulate multiple virulence mechanisms. Despite of a propensity for chaotic behaviour, no chaotic motifs were readily observed in the P. aeruginosa virulence regulatory network. Having a ‘birds-eye’ view of the regulatory cascades provides the forum opportunities to pose questions, formulate hypotheses and evaluate theories in elucidating P. aeruginosa pathogenesis. Understanding the mechanisms involved in making P. aeruginosa a successful pathogen is essential in helping devise control strategies.

A misplaced lncRNA causes brachydactyly in humans

Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis- and in trans-acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.

Long Noncoding RNAs in Cardiac Development and Pathophysiology

Heart function requires sophisticated regulatory networks to orchestrate organ development, physiological responses, and environmental adaptation. Until recently, it was thought that these regulatory networks are composed solely of protein-mediated transcriptional control and signaling systems; consequently, it was thought that cardiac disease involves perturbation of these systems. However, it is becoming evident that RNA, long considered to function primarily as the platform for protein production, may in fact play a major role in most, if not all, aspects of gene regulation, especially the epigenetic processes that underpin organogenesis. These include not only well-validated classes of regulatory RNAs, such as microRNAs, but also tens of thousands of long noncoding RNAs that are differentially expressed across the entire genome of humans and other animals. Here, we review this emerging landscape, summarizing what is known about their functions and their role in cardiac biology, and provide a toolkit to assist in exploring this previously hidden layer of gene regulation that may underpin heart adaptation and complex heart diseases.

A novel mammal-specific three partite enhancer element regulates node and notochord-specific Noto expression

The vertebrate organizer and notochord have conserved, essential functions for embryonic development and patterning. The restricted expression of developmental regulators in these tissues is directed by specific cis-regulatory modules (CRMs) whose sequence conservation varies considerably. Some CRMs have been conserved throughout vertebrates and likely represent ancestral regulatory networks, while others have diverged beyond recognition but still function over a wide evolutionary range. Here we identify and characterize a mammalian-specific CRM required for node and notochord specific (NNC) expression of NOTO, a transcription factor essential for node morphogenesis, nodal cilia movement and establishment of laterality in mouse. A 523 bp enhancer region (NOCE) upstream the Noto promoter was necessary and sufficient for NNC expression from the endogenous Noto locus. Three subregions in NOCE together mediated full activity in vivo. Binding sites for known transcription factors in NOCE were functional in vitro but dispensable for NOCE activity in vivo. A FOXA2 site in combination with a novel motif was necessary for NOCE activity in vivo. Strikingly, syntenic regions in non-mammalian vertebrates showed no recognizable sequence similarities. In contrast to its activity in mouse NOCE did not drive NNC expression in transgenic fish. NOCE represents a novel, mammal-specific CRM required for the highly restricted Noto expression in the node and nascent notochord and thus regulates normal node development and function.

A tissue-specific chromatin loop activates the erythroid ankyrin-1 promoter

The human ankyrin-1 gene (ANK1) contains 3 tissue-specific alternative promoters. We have shown previously that the erythroid-specific ankyrin 1 (ANK1E) core promoter contains a 5' DNase I hypersensitive site (HS) with barrier insulator function that prevents gene silencing in vitro and in vivo. Mutations in the ANK1E barrier region lead to decreased ANK1 mRNA levels and hereditary spherocytosis. In this report, we demonstrate a second ANK1E regulatory element located in an adjacent pair of DNase I HS located 5.6 kb 3' of the ANK1E promoter at the 3' boundary of an erythroid-specific DNase I-sensitive chromatin domain. The 3' regulatory element exhibits enhancer activity in vitro and in transgenic mice, and it has the histone modifications associated with an enhancer element. One of the ANK1E 3'HS contains an NF-E2 binding site that is required for enhancer function. We show that a chromatin loop brings the 3' enhancer and NF-E2 into proximity with the 5' barrier region including the ANK1E core promoter. These observations demonstrate a model for the tissue-specific activation of alternative promoters that may be applicable to the ∼ 30% of mammalian genes with alternative promoters that exhibit distinct expression patterns.

Theria-specific homeodomain and cis-regulatory element evolution of the Dlx3-4 bigene cluster in 12 different mammalian species

The mammalian Dlx3 and Dlx4 genes are configured as a bigene cluster, and their respective expression patterns are controlled temporally and spatially by cis-elements that largely reside within the intergenic region of the cluster. Previous work revealed that there are conspicuously conserved elements within the intergenic region of the Dlx3-4 bigene clusters of mouse and human. In this paper we have extended these analyses to include 12 additional mammalian taxa (including a marsupial and a monotreme) in order to better define the nature and molecular evolutionary trends of the coding and non-coding functional elements among morphologically divergent mammals. Dlx3-4 regions were fully sequenced from 12 divergent taxa of interest. We identified three theria-specific amino acid replacements in homeodomain of Dlx4 gene that functions in placenta. Sequence analyses of constrained nucleotide sites in the intergenic non-coding region showed that many of the intergenic conserved elements are highly conserved and have evolved slowly within the mammals. In contrast, a branchial arch/craniofacial enhancer I37-2 exhibited accelerated evolution at the branch between the monotreme and therian common ancestor despite being highly conserved among therian species. Functional analysis of I37-2 in transgenic mice has shown that the equivalent region of the platypus fails to drive transcriptional activity in branchial arches. These observations, taken together with our molecular evolutionary data, suggest that theria-specific episodic changes in the I37-2 element may have contributed to craniofacial innovation at the base of the mammalian lineage.

Small and long non-coding RNAs in cardiac homeostasis and regeneration

Cardiovascular diseases and in particular heart failure are major causes of morbidity and mortality in the Western world. Recently, the notion of promoting cardiac regeneration as a means to replace lost cardiomyocytes in the damaged heart has engendered considerable research interest. These studies envisage the utilization of both endogenous and exogenous cellular populations, which undergo highly specialized cell fate transitions to promote cardiomyocyte replenishment. Such transitions are under the control of regenerative gene regulatory networks, which are enacted by the integrated execution of specific transcriptional programs. In this context, it is emerging that the non-coding portion of the genome is dynamically transcribed generating thousands of regulatory small and long non-coding RNAs, which are central orchestrators of these networks. In this review, we discuss more particularly the biological roles of two classes of regulatory non-coding RNAs, i.e. microRNAs and long non-coding RNAs, with a particular emphasis on their known and putative roles in cardiac homeostasis and regeneration. Indeed, manipulating non-coding RNA-mediated regulatory networks could provide keys to unlock the dormant potential of the mammalian heart to regenerate. This should ultimately improve the effectiveness of current regenerative strategies and discover new avenues for repair. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Cardiac Pathways of Differentiation, Metabolism and Contraction.

Noncoding transcription at enhancers: general principles and functional models

Mammalian genomes are extensively transcribed outside the borders of protein-coding genes. Genome-wide studies recently demonstrated that cis-regulatory genomic elements implicated in transcriptional control, such as enhancers and locus-control regions, represent major sites of extragenic noncoding transcription. Enhancer-templated transcripts provide a quantitatively small contribution to the total amount of cellular nonribosomal RNA; nevertheless, the possibility that enhancer transcription and the resulting enhancer RNAs may, in some cases, have functional roles, rather than represent mere transcriptional noise at accessible genomic regions, is supported by an increasing amount of experimental data. In this article we review the current knowledge on enhancer transcription and its functional implications.

Transcription factors: from enhancer binding to developmental control

Developmental progression is driven by specific spatiotemporal domains of gene expression, which give rise to stereotypically patterned embryos even in the presence of environmental and genetic variation. Views of how transcription factors regulate gene expression are changing owing to recent genome-wide studies of transcription factor binding and RNA expression. Such studies reveal patterns that, at first glance, seem to contrast with the robustness of the developmental processes they encode. Here, we review our current knowledge of transcription factor function from genomic and genetic studies and discuss how different strategies, including extensive cooperative regulation (both direct and indirect), progressive priming of regulatory elements, and the integration of activities from multiple enhancers, confer specificity and robustness to transcriptional regulation during development.

Mechanisms of transcriptional precision in animal development

We review recently identified mechanisms of transcriptional control that ensure reliable and reproducible patterns of gene expression in natural populations of developing embryos, despite inherent fluctuations in gene regulatory processes, variations in genetic backgrounds and exposure to diverse environmental conditions. These mechanisms are not responsible for switching genes on and off. Instead, they control the fine-tuning of gene expression and ensure regulatory precision. Several such mechanisms are discussed, including redundant binding sites within transcriptional enhancers, shadow enhancers, and 'poised' enhancers and promoters, as well as the role of 'redundant' gene interactions within regulatory networks. We propose that such regulatory mechanisms provide population fitness and 'fine-tune' the spatial and temporal control of gene expression.

HOT DNAs: a novel class of developmental enhancers

Enhancers mediate localized patterns of gene expression during development. A common feature of "traditional" enhancers is the presence of clustered binding motifs for sequence-specific transcription factors (TFs). In this issue of Genes & Development, Kvon and colleagues (pp. 908-913) present new evidence that HOT (highly occupied transcription) DNAs direct specific patterns of gene expression, despite being depleted for TF-binding motifs.

Transcriptional integration of Wnt and Nodal pathways in establishment of the Spemann organizer

Signaling inputs from multiple pathways are essential for the establishment of distinct cell and tissue types in the embryo. Therefore, multiple signals must be integrated to activate gene expression and confer cell fate, but little is known about how this occurs at the level of target gene promoters. During early embryogenesis, Wnt and Nodal signals are required for formation of the Spemann organizer, which is essential for germ layer patterning and axis formation. Signaling by both Wnt and Nodal pathways is required for the expression of multiple organizer genes, suggesting that integration of these signals is required for organizer formation. Here, we demonstrate transcriptional cooperation between the Wnt and Nodal pathways in the activation of the organizer genes Goosecoid (Gsc), Cerberus (Cer), and Chordin (Chd). Combined Wnt and Nodal signaling synergistically activates transcription of these organizer genes. Effectors of both pathways occupy the Gsc, Cer and Chd promoters and effector occupancy is enhanced with active Wnt and Nodal signaling. This suggests that, at organizer gene promoters, a stable transcriptional complex containing effectors of both pathways forms in response to combined Wnt and Nodal signaling. Consistent with this idea, the histone acetyltransferase p300 is recruited to organizer promoters in a Wnt and Nodal effector-dependent manner. Taken together, these results offer a mechanism for spatial and temporal restriction of organizer gene transcription by the integration of two major signaling pathways, thus establishing the Spemann organizer domain.

In silico identification and in vivo validation of a set of evolutionary conserved plant root-specific cis-regulatory elements

Marker genes are specifically expressed in a tissue, organ or time of development. Here we used a computational screen to identify marker genes of the root in Arabidopsis thaliana. We mined the existing transcriptome datasets for genes having high expression in roots while being low in all other organs under a wide range of growth conditions. We show that the root-specificity of these genes is conserved in the sister species Arabidopsis lyrata, indicating that their expression pattern is under selective pressure. We delineated the cis-regulatory elements responsible for root-specific expression and validated two third of those in planta as bona fide root-specific regulatory sequences. We identified three motifs over-represented in these sequences, which mutation resulted in alteration of root-specific expression, demonstrating that these motifs are functionally relevant. In addition, the three motifs are also over-represented in the cis-regulatory regions of the A. lyrata orthologs of our root-specific genes, and this despite an overall low degree of sequence conservation of these regions. Our results provide a resource to assess root-identity in the model genus Arabidopsis and shed light on the evolutionary history of gene regulation in plants.

Versatility of cooperative transcriptional activation: a thermodynamical modeling analysis for greater-than-additive and less-than-additive effects

We derive a statistical model of transcriptional activation using equilibrium thermodynamics of chemical reactions. We examine to what extent this statistical model predicts synergy effects of cooperative activation of gene expression. We determine parameter domains in which greater-than-additive and less-than-additive effects are predicted for cooperative regulation by two activators. We show that the statistical approach can be used to identify different causes of synergistic greater-than-additive effects: nonlinearities of the thermostatistical transcriptional machinery and three-body interactions between RNA polymerase and two activators. In particular, our model-based analysis suggests that at low transcription factor concentrations cooperative activation cannot yield synergistic greater-than-additive effects, i.e., DNA transcription can only exhibit less-than-additive effects. Accordingly, transcriptional activity turns from synergistic greater-than-additive responses at relatively high transcription factor concentrations into less-than-additive responses at relatively low concentrations. In addition, two types of re-entrant phenomena are predicted. First, our analysis predicts that under particular circumstances transcriptional activity will feature a sequence of less-than-additive, greater-than-additive, and eventually less-than-additive effects when for fixed activator concentrations the regulatory impact of activators on the binding of RNA polymerase to the promoter increases from weak, to moderate, to strong. Second, for appropriate promoter conditions when activator concentrations are increased then the aforementioned re-entrant sequence of less-than-additive, greater-than-additive, and less-than-additive effects is predicted as well. Finally, our model-based analysis suggests that even for weak activators that individually induce only negligible increases in promoter activity, promoter activity can exhibit greater-than-additive responses when transcription factors and RNA polymerase interact by means of three-body interactions. Overall, we show that versatility of transcriptional activation is brought about by nonlinearities of transcriptional response functions and interactions between transcription factors, RNA polymerase and DNA.