Cis-regulatory logic in the endo16 gene: switching from a specification to a differentiation mode of control

cis-Regulatory sequences driving the expression of the Hbox12 homeobox-containing gene in the presumptive aboral ectoderm territory of the Paracentrotus lividus sea urchin embryo

Embryonic development is coordinated by networks of evolutionary conserved regulatory genes encoding transcription factors and components of cell signalling pathways. In the sea urchin embryo, a number of genes encoding transcription factors display territorial restricted expression. Among these, the zygotic Hbox12 homeobox gene is transiently transcribed in a limited number of cells of the animal-lateral half of the early Paracentrotus lividus embryo, whose descendants will constitute part of the ectoderm territory. To obtain insights on the regulation of Hbox12 expression, we have explored the cis-regulatory apparatus of the gene. In this paper, we show that the intergenic region of the tandem Hbox12 repeats drives GFP expression in the presumptive aboral ectoderm and that a 234 bp fragment, defined aboral ectoderm (AE) module, accounts for the restricted expression of the transgene. Within this module, a consensus sequence for a Sox factor and the binding of the Otx activator are both required for correct Hbox12 gene expression. Spatial restriction to the aboral ectoderm is achieved by a combination of different repressive sequence elements. Negative sequence elements necessary for repression in the endomesoderm map within the most upstream 60 bp region and nearby the Sox binding site. Strikingly, a Myb-like consensus is necessary for repression in the oral ectoderm, while down-regulation at the gastrula stage depends on a GA-rich region. These results suggest a role for Hbox12 in aboral ectoderm specification and represent our first attempt in the identification of the gene regulatory circuits involved in this process.

The sea urchin genome as a window on function

The emphasis on the sequencing of genomes seems to make this task an end in itself. However, genome sequences and the genes that are predicted from them are really an opportunity to examine the biological function of the organism constructed by that genome. This point is illustrated here by examples in which the newly annotated gene complement reveals surprises about the way Strongylocentrotus purpuratus, the purple sea urchin, goes about its business. The three topics considered here are the nature of the innate immune system; the unexpected complexity of sensory function implied by genes encoding sensory proteins; and the remarkable intricacy of the regulatory gene complement in embryogenesis.

A quantitative spatiotemporal atlas of gene expression in the Drosophila blastoderm

To fully understand animal transcription networks, it is essential to accurately measure the spatial and temporal expression patterns of transcription factors and their targets. We describe a registration technique that takes image-based data from hundreds of Drosophila blastoderm embryos, each costained for a reference gene and one of a set of genes of interest, and builds a model VirtualEmbryo. This model captures in a common framework the average expression patterns for many genes in spite of significant variation in morphology and expression between individual embryos. We establish the method's accuracy by showing that relationships between a pair of genes' expression inferred from the model are nearly identical to those measured in embryos costained for the pair. We present a VirtualEmbryo containing data for 95 genes at six time cohorts. We show that known gene-regulatory interactions can be automatically recovered from this data set and predict hundreds of new interactions.

Global regulatory logic for specification of an embryonic cell lineage

Explanation of a process of development must ultimately be couched in the terms of the genomic regulatory code. Specification of an embryonic cell lineage is driven by a network of interactions among genes encoding transcription factors. Here, we present the gene regulatory network (GRN) that directs the specification of the skeletogenic micromere lineage of the sea urchin embryo. The GRN now includes all regulatory genes expressed in this lineage up to late blastula stage, as identified in a genomewide survey. The architecture of the GRN was established by a large-scale perturbation analysis in which the expression of each gene in the GRN was cut off by use of morpholinos, and the effects on all other genes were measured quantitatively. Several cis-regulatory analyses provided additional evidence. The explanatory power of the GRN suffices to provide a causal explanation for all observable developmental functions of the micromere lineage during the specification period. These functions are: (i) initial acquisition of identity through transcriptional interpretation of localized maternal cues; (ii) activation of specific regulatory genes by use of a double negative gate; (iii) dynamic stabilization of the regulatory state by activation of a feedback subcircuit; (iv) exclusion of alternative regulatory states; (v) presentation of a signal required by the micromeres themselves and of two different signals required for development of adjacent endomesodermal lineages; and (vi) lineage-specific activation of batteries of skeletogenic genes. The GRN precisely predicts gene expression responses and provides a coherent explanation of the biology of specification.

Bio-logic: gene expression and the laws of combinatorial logic

At the heart of the development of fertilized eggs into fully formed organisms and the adaptation of cells to changed conditions are genetic regulatory networks (GRNs). In higher multicellular organisms, signal selection and multiplexing are performed at the cis-regulatory domains of genes, where combinations of transcription factors (TFs) regulate the rates at which the genes are transcribed into mRNA. To be able to act as activators or repressors of gene transcription, TFs must first bind to target sequences on the regulatory domains. Two TFs that act in concert may bind entirely independently of each other, but more often binding of the first one will alter the affinity of the other for its binding site. This article presents a systematic investigation into the effect of TF binding dependences on the predicted regulatory function of this bio-logic. Four extreme scenarios, commonly used to classify enzyme activation and inhibition patterns, for the binding of two TFs were explored: independent (the TFs bind without affecting each other's affinities), competitive (the TFs compete for the same binding site), ordered (the TFs bind in a compulsory order), and joint binding (the TFs either bind as a preformed complex, or binding of one is virtually impossible in the absence of the other). The conclusions are: (1) the laws of combinatorial logic hold only for systems with independently binding TFs; (2) systems formed according to the other scenarios can mimic the functions of their Boolean logical counterparts, but cannot be combined or decomposed in the same way; and (3) the continuously scaled output of systems consisting of competitively binding activators and repressors can be controlled more robustly than that of single TF or (quasi-)logical multi-TF systems.

TABASCO: A single molecule, base-pair resolved gene expression simulator

Background

Experimental studies of gene expression have identified some of the individual molecular components and elementary reactions that comprise and control cellular behavior. Given our current understanding of gene expression, and the goals of biotechnology research, both scientists and engineers would benefit from detailed simulators that can explicitly compute genome-wide expression levels as a function of individual molecular events, including the activities and interactions of molecules on DNA at single base pair resolution. However, for practical reasons including computational tractability, available simulators have not been able to represent genome-scale models of gene expression at this level of detail.

Results

Here we develop a simulator, TABASCO , which enables the precise representation of individual molecules and events in gene expression for genome-scale systems. We use a single molecule computational engine to track individual molecules interacting with and along nucleic acid polymers at single base resolution. Tabasco uses logical rules to automatically update and delimit the set of species and reactions that comprise a system during simulation, thereby avoiding the need for a priori specification of all possible combinations of molecules and reaction events. We confirm that single molecule, base-pair resolved simulation using TABASCO (Tabasco) can accurately compute gene expression dynamics and, moving beyond previous simulators, provide for the direct representation of intermolecular events such as polymerase collisions and promoter occlusion. We demonstrate the computational capacity of Tabasco by simulating the entirety of gene expression during bacteriophage T7 infection; for reference, the 39,937 base pair T7 genome encodes 56 genes that are transcribed by two types of RNA polymerases active across 22 promoters.

Conclusion

Tabasco enables genome-scale simulation of transcription and translation at individual molecule and single base-pair resolution. By directly representing the position and activity of individual molecules on DNA, Tabasco can directly test the effects of detailed molecular processes on system-wide gene expression. Tabasco would also be useful for studying the complex regulatory mechanisms controlling eukaryotic gene expression. The computational engine underlying Tabasco could also be adapted to represent other types of processive systems in which individual reaction events are organized across a single spatial dimension (e.g., polysaccharide synthesis).

Threshold-dominated regulation hides genetic variation in gene expression networks

Background

In dynamical models with feedback and sigmoidal response functions, some or all variables have thresholds around which they regulate themselves or other variables. A mathematical analysis has shown that when the dose-response functions approach binary or on/off responses, any variable with an equilibrium value close to one of its thresholds is very robust to parameter perturbations of a homeostatic state. We denote this threshold robustness. To check the empirical relevance of this phenomenon with response function steepnesses ranging from a near on/off response down to Michaelis-Menten conditions, we have performed a simulation study to investigate the degree of threshold robustness in models for a three-gene system with one downstream gene, using several logical input gates, but excluding models with positive feedback to avoid multistationarity. Varying parameter values representing functional genetic variation, we have analysed the coefficient of variation (CV) of the gene product concentrations in the stable state for the regulating genes in absolute terms and compared to the CV for the unregulating downstream gene. The sigmoidal or binary dose-response functions in these models can be considered as phenomenological models of the aggregated effects on protein or mRNA expression rates of all cellular reactions involved in gene expression.

Results

For all the models, threshold robustness increases with increasing response steepness. The CVs of the regulating genes are significantly smaller than for the unregulating gene, in particular for steep responses. The effect becomes less prominent as steepnesses approach Michaelis-Menten conditions. If the parameter perturbation shifts the equilibrium value too far away from threshold, the gene product is no longer an effective regulator and robustness is lost. Threshold robustness arises when a variable is an active regulator around its threshold, and this function is maintained by the feedback loop that the regulator necessarily takes part in and also is regulated by. In the present study all feedback loops are negative, and our results suggest that threshold robustness is maintained by negative feedback which necessarily exists in the homeostatic state.

Conclusion

Threshold robustness of a variable can be seen as its ability to maintain an active regulation around its threshold in a homeostatic state despite external perturbations. The feedback loop that the system necessarily possesses in this state, ensures that the robust variable is itself regulated and kept close to its threshold. Our results suggest that threshold regulation is a generic phenomenon in feedback-regulated networks with sigmoidal response functions, at least when there is no positive feedback. Threshold robustness in gene regulatory networks illustrates that hidden genetic variation can be explained by systemic properties of the genotype-phenotype map.

A cis-regulatory logic simulator

Background

A major goal of computational studies of gene regulation is to accurately predict the expression of genes based on the cis-regulatory content of their promoters. The development of computational methods to decode the interactions among cis-regulatory elements has been slow, in part, because it is difficult to know, without extensive experimental validation, whether a particular method identifies the correct cis-regulatory interactions that underlie a given set of expression data. There is an urgent need for test expression data in which the interactions among cis-regulatory sites that produce the data are known. The ability to rapidly generate such data sets would facilitate the development and comparison of computational methods that predict gene expression patterns from promoter sequence.

Results

We developed a gene expression simulator which generates expression data using user-defined interactions between cis-regulatory sites. The simulator can incorporate additive, cooperative, competitive, and synergistic interactions between regulatory elements. Constraints on the spacing, distance, and orientation of regulatory elements and their interactions may also be defined and Gaussian noise can be added to the expression values. The simulator allows for a data transformation that simulates the sigmoid shape of expression levels from real promoters. We found good agreement between sets of simulated promoters and predicted regulatory modules from real expression data. We present several data sets that may be useful for testing new methodologies for predicting gene expression from promoter sequence.

Conclusion

We developed a flexible gene expression simulator that rapidly generates large numbers of simulated promoters and their corresponding transcriptional output based on specified interactions between cis-regulatory sites. When appropriate rule sets are used, the data generated by our simulator faithfully reproduces experimentally derived data sets. We anticipate that using simulated gene expression data sets will facilitate the direct comparison of computational strategies to predict gene expression from promoter sequence. The source code is available online and as additional material. The test sets are available as additional material.

Beyond the gene

This paper is a response to the increasing difficulty biologists find in agreeing upon a definition of the gene, and indeed, the increasing disarray in which that concept finds itself. After briefly reviewing these problems, we propose an alternative to both the concept and the word gene—an alternative that, like the gene, is intended to capture the essence of inheritance, but which is both richer and more expressive. It is also clearer in its separation of what the organism statically is (what it tangibly inherits) and what it dynamically does (its functionality and behavior). Our proposal of a genetic functor, or genitor, is a sweeping extension of the classical genotype/phenotype paradigm, yet it appears to be faithful to the findings of contemporary biology, encompassing many of the recently emerging—and surprisingly complex—links between structure and functionality.

Sequential logic model deciphers dynamic transcriptional control of gene expressions

Background

Cellular signaling involves a sequence of events from ligand binding to membrane receptors through transcription factors activation and the induction of mRNA expression. The transcriptional-regulatory system plays a pivotal role in the control of gene expression. A novel computational approach to the study of gene regulation circuits is presented here.

Methodology

Based on the concept of finite state machine, which provides a discrete view of gene regulation, a novel sequential logic model (SLM) is developed to decipher control mechanisms of dynamic transcriptional regulation of gene expressions. The SLM technique is also used to systematically analyze the dynamic function of transcriptional inputs, the dependency and cooperativity, such as synergy effect, among the binding sites with respect to when, how much and how fast the gene of interest is expressed.

Principal Findings

SLM is verified by a set of well studied expression data on endo16 of Strongylocentrotus purpuratus (sea urchin) during the embryonic midgut development. A dynamic regulatory mechanism for endo16 expression controlled by three binding sites, UI, R and Otx is identified and demonstrated to be consistent with experimental findings. Furthermore, we show that during transition from specification to differentiation in wild type endo16 expression profile, SLM reveals three binary activities are not sufficient to explain the transcriptional regulation of endo16 expression and additional activities of binding sites are required. Further analyses suggest detailed mechanism of R switch activity where indirect dependency occurs in between UI activity and R switch during specification to differentiation stage.

Conclusions/Significance

The sequential logic formalism allows for a simplification of regulation network dynamics going from a continuous to a discrete representation of gene activation in time. In effect our SLM is non-parametric and model-independent, yet providing rich biological insight. The demonstration of the efficacy of this approach in endo16 is a promising step for further application of the proposed method.

The regulatory genome and the computer

The definitive feature of the many thousand cis-regulatory control modules in an animal genome is their information processing capability. These modules are "wired" together in large networks that control major processes such as development; they constitute "genomic computers." Each control module receives multiple inputs in the form of the incident transcription factors which bind to them. The functions they execute upon these inputs can be reduced to basic AND, OR and NOT logic functions, which are also the unit logic functions of electronic computers. Here we consider the operating principles of the genomic computer, the product of evolution, in comparison to those of electronic computers. For example, in the genomic computer intra-machine communication occurs by means of diffusion (of transcription factors), while in electronic computers it occurs by electron transit along pre-organized wires. There follow fundamental differences in design principle in respect to the meaning of time, speed, multiplicity of processors, memory, robustness of computation and hardware and software. The genomic computer controls spatial gene expression in the development of the body plan, and its appearance in remote evolutionary time must be considered to have been a founding requirement for animal grade life.

Gene regulation: gene control network in development

Controlling the differential expression of many thousands of genes is the most fundamental task of a developing organism. It requires an enormous computational device that has the capacity to process in parallel a vast number of regulatory inputs in the various cells of the embryo and come out with regulatory outputs that are tissue specific. The regulatory genome constitutes this computational device, comprising many thousands of processing units in the form of cis-regulatory modules. The interconnected cis-regulatory modules that control regulatory gene expression create a network that is the underlying mechanism of specification. In this review we use the gene regulatory network that governs endomesoderm specification in the sea urchin embryo to demonstrate the salient features of developmental gene regulatory networks and illustrate the information processing that is done by the regulatory sequences.

Endo16 is required for gastrulation in the sea urchin Lytechinus variegatus

The Endo16 gene encodes a large extracellular protein with several functional domains that provide some insight into the role of this protein during embryonic development. We isolated the full-length cDNA sequence from Lytechinus variegatus and utilized morpholinos to further investigate the role of Endo16 during embryonic development in this species. Endo16-deficient embryos failed to undergo gastrulation and the blastocoele became filled with dissociated cells after 24 h of incubation. Moreover, there was a delay in endoderm differentiation as assayed by staining with an antibody that recognizes Endo1. The differentiation of other cell types including oral ectoderm, primary mesenchymal cells (PMC) and secondary mesenchymal cells (SMC) appeared to be normal, although the patterns of protein expression did not resemble control embryos due to the gross morphological abnormalities elicited by the LvEndo16 morpholino. Microinjection of full-length EGFP mRNA with the LvEndo16 morpholino-targeted sequence confirmed that this phenotype can be attributed specifically to the loss of Endo16 protein. Taken together, our data suggest that Endo16 may be required for the cell-extracellular matrix (ECM) interactions that are required for endoderm differentiation in the sea urchin embryo.

High regulatory gene use in sea urchin embryogenesis: Implications for bilaterian development and evolution

A global scan of transcription factor usage in the sea urchin embryo was carried out in the context of the Strongylocentrotus purpuratus genome sequencing project, and results from six individual studies are here considered. Transcript prevalence data were obtained for over 280 regulatory genes encoding sequence-specific transcription factors of every known family, but excluding genes encoding zinc finger proteins. This is a statistically inclusive proxy for the total "regulome" of the sea urchin genome. Close to 80% of the regulome is expressed at significant levels by the late gastrula stage. Most regulatory genes must be used repeatedly for different functions as development progresses. An evolutionary implication is that animal complexity at the stage when the regulome first evolved was far simpler than even the last common bilaterian ancestor, and is thus of deep antiquity.

Analysis of gene regulatory network models with graded and binary transcriptional responses

The steep sigmoid framework developed by Plahte and Kjøglum [Plahte, E., Kjøglum, S., 2005. Analysis and generic properties of gene regulatory networks with graded response functions. Phys. D 201, 150-176, doi:10.1016/j.physd.2004.11.014] provides a uniform description of gene regulatory networks in which there may be both graded and binary transcriptional responses, as well as a method for analysing the models developed. Here we extend this framework. We show that there is a relation between the location of steady states and the feedback structure of a system, thus generalising existing results for Boolean type models. In addition, we justify underlying assumptions and generic features of the modelling framework in terms of biology and generalise the overall approach to take into account that each transcription factor only regulates one gene at a given threshold. By this assumption, the analysis of the models are greatly simplified.

Quantitative and predictive model of transcriptional control of the Drosophila melanogaster even skipped gene

Here we present a quantitative and predictive model of the transcriptional readout of the proximal 1.7 kb of the control region of the Drosophila melanogaster gene even skipped (eve). The model is based on the positions and sequence of individual binding sites on the DNA and quantitative, time-resolved expression data at cellular resolution. These data demonstrated new expression features, first reported here. The model correctly predicts the expression patterns of mutations in trans, as well as point mutations, insertions and deletions in cis. It also shows that the nonclassical expression of stripe 7 driven by this fragment is activated by the protein Caudal (Cad), and repressed by the proteins Tailless (Tll) and Giant (Gt).

Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity

Controversy surrounds the apparent rising maximums of morphological complexity during eukaryotic evolution, with organisms increasing the number and nestedness of developmental areas as evidenced by morphological elaborations reflecting area boundaries. No "predictable drive" to increase this sort of complexity has been reported. Recent genetic data and theory in the general area of gene dosage effects has engendered a robust "gene balance hypothesis," with a theoretical base that makes specific predictions as to gene content changes following different types of gene duplication. Genomic data from both chordate and angiosperm genomes fit these predictions: Each type of duplication provides a one-way injection of a biased set of genes into the gene pool. Tetraploidies and balanced segments inject bias for those genes whose products are the subunits of the most complex biological machines or cascades, like transcription factors (TFs) and proteasome core proteins. Most duplicate genes are removed after tetraploidy. Genic balance is maintained by not removing those genes that are dose-sensitive, which tends to leave duplicate "functional modules" as the indirect products (spandrels) of purifying selection. Functional modules are the likely precursors of coadapted gene complexes, a unit of natural selection. The result is a predictable drive mechanism where "drive" is used rigorously, as in "meiotic drive." Rising morphological gain is expected given a supply of duplicate functional modules. All flowering plants have survived at least three large-scale duplications/diploidizations over the last 300 million years (Myr). An equivalent period of tetraploidy and body plan evolution may have ended for animals 500 million years ago (Mya). We argue that "balanced gene drive" is a sufficient explanation for the trend that the maximums of morphological complexity have gone up, and not down, in both plant and animal eukaryotic lineages.

Comparative evaluation of reverse engineering gene regulatory networks with relevance networks, graphical gaussian models and bayesian networks

MOTIVATION

An important problem in systems biology is the inference of biochemical pathways and regulatory networks from postgenomic data. Various reverse engineering methods have been proposed in the literature, and it is important to understand their relative merits and shortcomings. In the present paper, we compare the accuracy of reconstructing gene regulatory networks with three different modelling and inference paradigms: (1) Relevance networks (RNs): pairwise association scores independent of the remaining network; (2) graphical Gaussian models (GGMs): undirected graphical models with constraint-based inference, and (3) Bayesian networks (BNs): directed graphical models with score-based inference. The evaluation is carried out on the Raf pathway, a cellular signalling network describing the interaction of 11 phosphorylated proteins and phospholipids in human immune system cells. We use both laboratory data from cytometry experiments as well as data simulated from the gold-standard network. We also compare passive observations with active interventions.

RESULTS

On Gaussian observational data, BNs and GGMs were found to outperform RNs. The difference in performance was not significant for the non-linear simulated data and the cytoflow data, though. Also, we did not observe a significant difference between BNs and GGMs on observational data in general. However, for interventional data, BNs outperform GGMs and RNs, especially when taking the edge directions rather than just the skeletons of the graphs into account. This suggests that the higher computational costs of inference with BNs over GGMs and RNs are not justified when using only passive observations, but that active interventions in the form of gene knockouts and over-expressions are required to exploit the full potential of BNs.

AVAILABILITY

Data, software and supplementary material are available from http://www.bioss.sari.ac.uk/staff/adriano/research.html

Transcription binding site prediction using Markov models

One of the main goals of analysing DNA sequences is to understand the temporal and positional information that specifies gene expression. An important step in this process is the recognition of gene expression regulatory elements. Experimental procedures for this are slow and costly. In this paper we present a computational non-supervised algorithm that facilitates the process by statistically identifying the most likely regions within a putative regulatory sequence. A probabilistic technique is presented, based on the approximation of regulatory DNA with a Markov chain, for the location of putative transcription factor binding sites in a single stretch of DNA. Hereto we developed a procedure to approximate the order of Markov model for a given DNA sequence that circumvents some of the prohibitive assumptions underlying Markov modeling. Application of the algorithm to data from 55 genes in five species shows the high sensitivity of this Markov search algorithm. Our algorithm does not require any prior knowledge in the form of description or cross-genomic comparison; it is context sensitive and takes DNA heterogeneity into account.

Confocal quantification of cis-regulatory reporter gene expression in living sea urchin

Quantification of GFP reporter gene expression at single cell level in living sea urchin embryos can now be accomplished by a new method of confocal laser scanning microscopy (CLSM). Eggs injected with a tissue-specific GFP reporter DNA construct were grown to gastrula stage and their fluorescence recorded as a series of contiguous Z-section slices that spanned the entire embryo. To measure the depth-dependent signal decay seen in the successive slices of an image stack, the eggs were coinjected with a freely diffusible internal fluorescent standard, rhodamine dextran. The measured rhodamine fluorescence was used to generate a computational correction for the depth-dependent loss of GFP fluorescence per slice. The intensity of GFP fluorescence was converted to the number of GFP molecules using a conversion constant derived from CLSM imaging of eggs injected with a measured quantity of GFP protein. The outcome is a validated method for accurately counting GFP molecules in given cells in reporter gene transfer experiments, as we demonstrate by use of an expression construct expressed exclusively in skeletogenic cells.

cis-regulatory processing of Notch signaling input to the sea urchin glial cells missing gene during mesoderm specification

The glial cells missing regulatory gene of Strongylocentrotus purpuratus (spgcm) was proposed earlier to be the genomic target of Delta/Notch (D/N) signaling required for specification of the mesodermal precursors of pigment cells. Here, we show that microinjection of a spgcm antisense morpholino oligonucleotide results in larvae without pigment cells. Microinjection of an mRNA encoding a dominant negative form of Suppressor of Hairless (dn-Su(H)) results in reduced levels of spgcm mRNA, disruption of mesodermal founder cell specification and failure to produce pigment cells. These results confirm that this gene is required for pigment cell specification. Three cis-regulatory modules of the spgcm gene were identified, which when incorporated in a GFP expression construct recapitulate the early expression pattern of this gene. Spatial expression of this GFP expression construct is severely disrupted by co-expression of dn-Su(H) mRNA, confirming that spgcm is a direct target of canonical N signaling mediated through Su(H) inputs. cis-perturbation analysis by mutation of consensus Su(H) sites identified a conserved motif paired-site and a lone site in the middle module that function both to drive expression in SMC precursors which receive the Delta signal and to repress expression in ectopic locations which lack this signal. While these Su(H) target sites provide the cis-regulatory architecture with the core of an N signaling transcriptional response switch, both the on and off outputs from this module require additional inputs.

Developmental mechanisms and cis-regulatory codes

Complex networks of transcriptional interactions control the processes of animal development. These networks begin with broad positional information that patterns the cells of the early embryo, and end with precise expression profiles that provide the functions of fully differentiated cells. At the heart of these networks are cis-regulatory modules (CRMs), which contain binding sites for regulatory proteins and control the spatial and temporal expression of genes within the network. Recent studies in several model systems have begun to decipher the 'cis-regulatory codes' of CRMs involved in various developmental processes. These studies suggest that CRMs involved in regulating co-expressed genes share sequence characteristics that can be identified by in silico approaches. They also suggest that CRMs involved in specific types of developmental events have common binding site architectures, which can be linked to their specific functions.

cis-Regulatory inputs of the wnt8 gene in the sea urchin endomesoderm network

Expression of the wnt8 gene is the key transcriptional motivator of an intercellular signaling loop which drives endomesoderm specification forward early in sea urchin embryogenesis. This gene was predicted by network perturbation analysis to be activated by inputs from the blimp1/krox gene, itself expressed zygotically in the endomesoderm during cleavage; and by a Tcf1/beta-catenin input. The implication is that zygotic expression of wnt8 is stimulated in neighboring cells by its own gene product, since reception of the Wnt8 ligand causes beta-catenin nuclearization. Here, the modular cis-regulatory system of the wnt8 gene of Strongylocentrotus purpuratus was characterized functionally, and shown to respond to blockade of both Blimp1/Krox and Tcf1/beta-catenin inputs just as does the endogenous gene. The genomic target sites for these factors were demonstrated by mutation in one of the cis-regulatory modules. The Tcf1/beta-catenin and Blimp1/Krox inputs are both necessary for normal endomesodermal expression mediated by this cis-regulatory module; thus, the genomic regulatory code underlying the predicted signaling loop thus resides in the wnt8 cis-regulatory sequence. In a second regulatory region, which initiates expression in micromere and macromere descendant cells early in cleavage, Tcf1 sites act to repress ectopic transcription in prospective ectoderm cells.

Dynamic modeling of cis-regulatory circuits and gene expression prediction via cross-gene identification

Background

Gene expression programs depend on recognition of cis elements in promoter region of target genes by transcription factors (TFs), but how TFs regulate gene expression via recognition of cis elements is still not clear. To study this issue, we define the cis-regulatory circuit of a gene as a system that consists of its cis elements and the interactions among their recognizing TFs and develop a dynamic model to study the functional architecture and dynamics of the circuit. This is in contrast to traditional approaches where a cis-regulatory circuit is constructed by a mutagenesis or motif-deletion scheme. We estimate the regulatory functions of cis-regulatory circuits using microarray data.

Results

A novel cross-gene identification scheme is proposed to infer how multiple TFs coordinate to regulate gene transcription in the yeast cell cycle and to uncover hidden regulatory functions of a cis-regulatory circuit. Some advantages of this approach over most current methods are that it is based on data obtained from intact cis-regulatory circuits and that a dynamic model can quantitatively characterize the regulatory function of each TF and the interactions among the TFs. Our method may also be applicable to other genes if their expression profiles have been examined for a sufficiently long time.

Conclusion

In this study, we have developed a dynamic model to reconstruct cis-regulatory circuits and a cross-gene identification scheme to estimate the regulatory functions of the TFs that control the regulation of the genes under study. We have applied this method to cell cycle genes because the available expression profiles for these genes are long enough. Our method not only can quantify the regulatory strengths and synergy of the TFs but also can predict the expression profile of any gene having a subset of the cis elements studied.

Computational representation of developmental genetic regulatory networks

Developmental genetic regulatory networks (GRNs) have unique architectural characteristics. They are typically large-scale, multi-layered, and organized in a nested, modular hierarchy of regulatory network kernels, function-specific building blocks, and structural gene batteries. They are also inherently multicellular and involve changing topological relationships among a growing number of cells. Reconstruction of developmental GRNs requires unique computational tools that support the above representational requirements. In addition, we argue that DNA-centered network modeling, separate descriptions of network organization and network behavior, and support for network documentation and annotation are essential requirements for computational modeling of developmental GRNs. Based on these observations, we have developed a freely available, platform-independent, open source software package (BioTapestry) which supports both the process of model construction and also model visualization, analysis, documentation, and dissemination. We provide an overview of the main features of BioTapestry. The BioTapestry software and additional documents are available from http://www.biotapestry.org. We recommend BioTapestry as the substrate for further co-development for and by the developmental biology community.

Brn1/2/4, the predicted midgut regulator of the endo16 gene of the sea urchin embryo

A specific prediction of our detailed cis-regulatory analysis of the Strongylocentrotus purpuratus (Sp) endo16 gene was that the later expression of this gene would be driven by a midgut-specific transcriptional regulator. We have now identified this factor and determined some of its functions. The cDNA sequence reveals it to be a POU domain factor related closely to the mammalian factors Brain-1, -2, and -4. The factor was termed SpBrn1/2/4 (henceforth Brn1/2/4). Quantitative measurements of transcript prevalence show that the gene is first activated in the 20-h blastula, but there remain only about 100 molecules of brn1/2/4 mRNA per embryo (only a few per endoderm cell) until an abrupt 10-fold increase occurs as gastrulation begins. Measured in the same embryos, the late rise in prevalence of endo16 transcripts follows that of brn1/2/4 transcripts. As predicted by the endo16 model, brn1/2/4 expression is confined perfectly to the midgut, coincident with the domain of endo16 expression. The kinetics of accumulation of these transcripts indicates that the switch into the late phase of endo16 expression occurs when the brn1/2/4 transcript level nears its plateau (2000 molecules mRNA per embryo), after which each endo16 gene produces about 1 mRNA every 2 min (about 380 molecules mRNA per min in the whole embryo). Arrest of Brn1/2/4 translation by MASO treatment blocks the late phase of endo16 expression and specifically abolishes expression of cis-regulatory Module B of endo16, while not affecting Module A, also as predicted. The brn1/2/4 gene lies downstream of the regulatory genes executing post-gastrular specification of the midgut, as shown by further gene expression perturbation experiments which provide an initial glimpse of the underlying network architecture.

Gastrulation in the sea urchin embryo: a model system for analyzing the morphogenesis of a monolayered epithelium

Processes of gastrulation in the sea urchin embryo have been intensively studied to reveal the mechanisms involved in the invagination of a monolayered epithelium. It is widely accepted that the invagination proceeds in two steps (primary and secondary invagination) until the archenteron reaches the apical plate, and that the constituent cells of the resulting archenteron are exclusively derived from the veg2 tier of blastomeres formed at the 60-cell stage. However, recent studies have shown that the recruitment of the archenteron cells lasts as late as the late prism stage, and some descendants of veg1 blastomeres are also recruited into the archenteron. In this review, we first illustrate the current outline of sea urchin gastrulation. Second, several factors, such as cytoskeletons, cell contact and extracellular matrix, will be discussed in relation to the cellular and mechanical basis of gastrulation. Third, differences in the manner of gastrulation among sea urchin species will be described; in some species, the archenteron does not elongate stepwise but continuously. In those embryos, bottle cells are scarcely observed, and the archenteron cells are not rearranged during invagination unlike in typical sea urchins. Attention will be also paid to some other factors, such as the turgor pressure of blastocoele and the force generated by blastocoele wall. These factors, in spite of their significance, have been neglected in the analysis of sea urchin gastrulation. Lastly, we will discuss how behavior of pigment cells defines the manner of gastrulation, because pigment cells recently turned out to be the bottle cells that trigger the initial inward bending of the vegetal plate.

R11: a cis-regulatory node of the sea urchin embryo gene network that controls early expression of SpDelta in micromeres

A gene regulatory network (GRN) controls the process by which the endomesoderm of the sea urchin embryo is specified. In this GRN, the program of gene expression unique to the skeletogenic micromere lineage is set in train by activation of the pmar1 gene. Through a double repression system, this gene is responsible for localization of expression of downstream regulatory and signaling genes to cells of this lineage. One of these genes, delta, encodes a Notch ligand, and its expression in the right place and time is crucial to the specification of the endomesoderm. Here we report a cis-regulatory element R11 that is responsible for localizing the expression of delta by means of its response to the pmar1 repression system. R11 was identified as an evolutionarily conserved genomic sequence located about 13 kb downstream of the last exon of the delta gene. We demonstrate here that this cis-regulatory element is able to drive the expression of a reporter gene in the same cells and at the same time that the endogenous delta gene is expressed, and that temporally, spatially, and quantitatively it responds to the pmar1 repression system just as predicted for the delta gene in the endomesoderm GRN. This work illustrates the application of cis-regulatory analysis to the validation of predictions of the GRN model. In addition, we introduce new methodological tools for quantitative measurement of the output of expression constructs that promise to be of general value for cis-regulatory analysis in sea urchin embryos.

cis-Regulatory activity of randomly chosen genomic fragments from the sea urchin

In order to determine the frequency and variety of cis-regulatory elements that function during embryonic development of Strongylocentrotus purpuratus, we constructed a GFP expression vector in which to test the activity of randomly chosen genomic DNA fragments that includes a promiscuous basal promoter from the endo16 gene. This vector was demonstrated to serve as a cis-regulatory element trap. We used it to carry out an initial test for the occurrence of elements that would promote GFP expression in this genome. In the screen reported here 108 different randomly chosen DNA fragments (av. 3.8 kb) were inserted in the vector, and each was injected into > 200 zygotes. Surprisingly, 13% of the fragments tested yielded detectable levels of GFP expression in the recipient embryos. Specific patterns observed included expression in endoderm, in aboral ectoderm, and in pigment cells. The majority of active constructs expressed GFP in all spatial domains of the embryo. Elements with detectable cis-regulatory activity in the embryo occur in the sample screened, on the average, about every 30 kb, and the genome must include many thousands of such elements. On further analysis one isolate was shown to contain a gut specific element as well as one that controls expression in the secondary mesenchyme cells.

cis-Regulatory control circuits in development

During development, an organism undergoes many rounds of pattern formation, generating ever-greater complexity with each ensuing round of cell division and specification. The instructions for executing this process are encoded in the cis-regulatory modules that direct the expression of developmental transcription factors and signaling molecules. Each transcription factor binding site within a cis-regulatory module contributes information about when, where, or how much a gene is turned on, and by dissecting the modules driving a given gene, all the inputs governing expression of the gene can be accurately identified. Furthermore, by mapping the output of each gene to the inputs of other genes, it is possible to reverse engineer developmental circuits and even whole networks. At this higher level of organization, common bilaterian strategies for specifying progenitor fields, locking down regulatory states, and driving development forward emerge.

Transcription of Drosophila troponin I gene is regulated by two conserved, functionally identical, synergistic elements

The Drosophila wings-up A gene encodes Troponin I. Two regions, located upstream of the transcription initiation site (upstream regulatory element) and in the first intron (intron regulatory element), regulate gene expression in specific developmental and muscle type domains. Based on LacZ reporter expression in transgenic lines, upstream regulatory element and intron regulatory element yield identical expression patterns. Both elements are required for full expression levels in vivo as indicated by quantitative reverse transcription-polymerase chain reaction assays. Three myocyte enhancer factor-2 binding sites have been functionally characterized in each regulatory element. Using exon specific probes, we show that transvection is based on transcriptional changes in the homologous chromosome and that Zeste and Suppressor of Zeste 3 gene products act as repressors for wings-up A. Critical regions for transvection and for Zeste effects are defined near the transcription initiation site. After in silico analysis in insects (Anopheles and Drosophila pseudoobscura) and vertebrates (Ratus and Coturnix), the regulatory organization of Drosophila seems to be conserved. Troponin I (TnI) is expressed before muscle progenitors begin to fuse, and sarcomere morphogenesis is affected by TnI depletion as Z discs fail to form, revealing a novel developmental role for the protein or its transcripts. Also, abnormal stoichiometry among TnI isoforms, rather than their absolute levels, seems to cause the functional muscle defects.

Gene Regulatory Network Analysis in Sea Urchin Embryos

It may safely be predicted that GRN analysis will become increasingly important. It will come to underlie the causal study of development, the major effort underway to understand the regulatory code built into animal genomes and also the evolution of these genomes. Partly by serendipity, sea urchin embryos turn out to be a superb experimental material for GRN analysis. Their natural properties have, in turn, influenced the predilections of those who work on them, and between them and us, so to speak, this is now a developmental system of which we are rapidly gaining an unusually complete understanding. The causal linkages that control development of the whole embryo will be revealed, leading all the way from the heritable genomic regulatory code to the events of embryology. The fundamental experimental operation is the perturbation analysis: Here is where causality permeates the exploration. We have in this chapter summarized in some detail the requirements for perturbation GRN analysis in sea urchin embryos. But that is not all, nor is it enough to enable the assembly of a GRN: What is required is the combined application of elegant computational methods, of gene regulation molecular biology, of genomic sequence data, and of experimental embryology. As the results crystallize together, we can begin to see how far this powerful combination of methods and ideas is going to carry us.

Using reporter genes to study cis-regulatory elements

This chapter summarizes four powerful assays for analyzing gene expression in cis-regulatory studies. The enzymatic assays (CAT, luciferase, lacZ) are currently limited by their application to embryo homogenates or fixed samples, but offer more robust analysis of gene activity than GFP. Assays based on CAT enzymatic activity or on CAT mRNA detection by WMISH are laborious but are well established for accurately quantifying gene expression and to determine spatial patterns at defined timepoints during development. LacZ assays are the current standard for spatially visualizing gene products in whole-mount fixed embryos. They are very sensitive but they provide limited temporal or quantitative information due to the perdurance of beta-galactosidase and the subtleties of the staining technique. Recently developed luciferase assays promise to be even more sensitive and accurate than the CAT and lacZ assays, and applicable to living cells and embryos. But, they have not yet been well established in invertebrate deuterostome research. GFP allows visualization of gene expression within living embryos. But because this is not an enzymatic assay, sensitivity can be a problem, particularly for weak promoters. Furthermore, imaging live embryos and quantifying gene expression in space and time (due to scattering of light by tissue, the perdurance of GFP, and other experimental details) is currently fraught with challenges. Ongoing improvements in imaging technology and the advent of multiple fluorescent proteins, as well as fluorescent and luminescent assays for vital imaging, will dramatically facilitate studies of gene expression in the coming decade.

Conservation of Endo16 expression in sea urchins despite evolutionary divergence in both cis and trans-acting components of transcriptional regulation

Evolutionary changes in transcriptional regulation undoubtedly play an important role in creating morphological diversity. However, there is little information about the evolutionary dynamics of cis-regulatory sequences. This study examines the functional consequence of evolutionary changes in the Endo16 promoter of sea urchins. The Endo16 gene encodes a large extracellular protein that is expressed in the endoderm and may play a role in cell adhesion. Its promoter has been characterized in exceptional detail in the purple sea urchin, Strongylocentrotus purpuratus. We have characterized the structure and function of the Endo16 promoter from a second sea urchin species, Lytechinus variegatus. The Endo16 promoter sequences have evolved in a strongly mosaic manner since these species diverged approximately 35 million years ago: the most proximal region (module A) is conserved, but the remaining modules (B-G) are unalignable. Despite extensive divergence in promoter sequences, the pattern of Endo16 transcription is largely conserved during embryonic and larval development. Transient expression assays demonstrate that 2.2 kb of upstream sequence in either species is sufficient to drive GFP reporter expression that correctly mimics this pattern of Endo16 transcription. Reciprocal cross-species transient expression assays imply that changes have also evolved in the set of transcription factors that interact with the Endo16 promoter. Taken together, these results suggest that stabilizing selection on the transcriptional output may have operated to maintain a similar pattern of Endo16 expression in S. purpuratus and L. variegatus, despite dramatic divergence in promoter sequence and mechanisms of transcriptional regulation.

Characterization of the upstream region that regulates the transcription of the gene for the precursor to EGF-related peptides, exogastrula-inducing peptides, of the sea urchin Anthocidaris crassispina

The EGIP gene for exogastrula-inducing peptides (EGIPs) of the sea urchin Anthocidaris crassispina, which are structurally related to the epidermal growth factor, is activated at the onset of gastrulation in subdomains of the embryonic ectoderm. We showed in our previous study that the spatial and temporal regulation of EGIP is conducted by the upstream region from -372 to +194, and that there is an enhancer element between -372 and -210. In this study, we introduced into sea urchin embryos PCR-amplified DNA containing differently truncated EGIP flanking region that was ligated to the GFP reporter gene, and examined the transient expression of the reporter gene, showing that both the -270/-238 and -249/-210 regions were essential for the enhancer activity. We further showed that there is another activating element between -65 and -21, and that even the region between -65 and +194 is sufficient for ectoderm-specific expression of the EGIP gene. The electrophoretic mobility shift assay showed that the -270/-210 enhancer region and the proximal -61/+30 region include specific binding sites for nuclear proteins of sea urchin embryos. Besides these unique sites, the presence of multiple binding sites for GCF1-like nuclear proteins have been revealed in the upstream DNA.

The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster

Expression of the Drosophila segment polarity genes is initiated by a pre-pattern of pair-rule gene products and maintained by a network of regulatory interactions throughout several stages of embryonic development. Analysis of a model of gene interactions based on differential equations showed that wild-type expression patterns of these genes can be obtained for a wide range of kinetic parameters, which suggests that the steady states are determined by the topology of the network and the type of regulatory interactions between components, not the detailed form of the rate laws. To investigate this, we propose and analyse a Boolean model of this network which is based on a binary ON/OFF representation of mRNA and protein levels, and in which the interactions are formulated as logical functions. In this model the spatial and temporal patterns of gene expression are determined by the topology of the network and whether components are present or absent, rather than the absolute levels of the mRNAs and proteins and the functional details of their interactions. The model is able to reproduce the wild-type gene expression patterns, as well as the ectopic expression patterns observed in overexpression experiments and various mutants. Furthermore, we compute explicitly all steady states of the network and identify the basin of attraction of each steady state. The model gives important insights into the functioning of the segment polarity gene network, such as the crucial role of the wingless and sloppy paired genes, and the network's ability to correct errors in the pre-pattern.

Patterning the sea urchin embryo: gene regulatory networks, signaling pathways, and cellular interactions

We discuss steps in the specification of major tissue territories of the sea urchin embryo that occur between fertilization and hatching blastula stage and the cellular interactions required to coordinate morphogenetic processes that begin after hatching. We review evidence that has led to new ideas about how this embryo is initially patterned: (1) Specification of most of the tissue territories is not direct, but proceeds gradually by progressive subdivision of broad, maternally specified domains that depend on opposing gradients in the ratios of animalizing transcription factors (ATFs) and vegetalizing (beta-catenin) transcription factors; (2) the range of maternal nuclear beta-catenin extends further than previously proposed, that is, into the animal hemisphere, where it programs many cells to adopt early aboral ectoderm characteristics; (3) cells at the extreme animal pole constitute a unique ectoderm region, lacking nuclear beta-catenin; (4) the pluripotential mesendoderm is created by the combined outputs of ATFs and nuclear beta-catenin, which initially overlap in the macromeres, and by an undefined early micromere signal; (5) later micromere signals, which activate Notch and Wnt pathways, subdivide mesendoderm into secondary mesenchyme and endoderm; and (6) oral ectoderm specification requires reprogramming early aboral ectoderm at about the hatching blastula stage. Morphogenetic processes that follow initial fate specification depend critically on continued interactions among cells in different territories. As illustrations, we discuss the regulation of (1) the ectoderm/endoderm boundary, (2) mesenchyme positioning and skeletal growth, (3) ciliated band formation, and (4) several suppressive interactions operating late in embryogenesis to limit the fates of multipotent cells.

On schemes of combinatorial transcription logic

Cells receive a wide variety of cellular and environmental signals, which are often processed combinatorially to generate specific genetic responses. Here we explore theoretically the potentials and limitations of combinatorial signal integration at the level of cis-regulatory transcription control. Our analysis suggests that many complex transcription-control functions of the type encountered in higher eukaryotes are already implementable within the much simpler bacterial transcription system. Using a quantitative model of bacterial transcription and invoking only specific protein-DNA interaction and weak glue-like interaction between regulatory proteins, we show explicit schemes to implement regulatory logic functions of increasing complexity by appropriately selecting the strengths and arranging the relative positions of the relevant protein-binding DNA sequences in the cis-regulatory region. The architectures that emerge are naturally modular and evolvable. Our results suggest that the transcription regulatory apparatus is a "programmable" computing machine, belonging formally to the class of Boltzmann machines. Crucial to our results is the ability to regulate gene expression at a distance. In bacteria, this can be achieved for isolated genes via DNA looping controlled by the dimerization of DNA-bound proteins. However, if adopted extensively in the genome, long-distance interaction can cause unintentional intergenic cross talk, a detrimental side effect difficult to overcome by the known bacterial transcription-regulation systems. This may be a key factor limiting the genome-wide adoption of complex transcription control in bacteria. Implications of our findings for combinatorial transcription control in eukaryotes are discussed.

Regulatory gene networks and the properties of the developmental process

Genomic instructions for development are encoded in arrays of regulatory DNA. These specify large networks of interactions among genes producing transcription factors and signaling components. The architecture of such networks both explains and predicts developmental phenomenology. Although network analysis is yet in its early stages, some fundamental commonalities are already emerging. Two such are the use of multigenic feedback loops to ensure the progressivity of developmental regulatory states and the prevalence of repressive regulatory interactions in spatial control processes. Gene regulatory networks make it possible to explain the process of development in causal terms and eventually will enable the redesign of developmental regulatory circuitry to achieve different outcomes.

Origins of differentiation via phenotypic plasticity

How cell types of multicellular organisms came to be differentiated is still an open issue. Here I offer a model that posits that the origins of some cell differentiation patterns were originally passive outcomes of environmental effects. As cells' contact with the external environment was diminished, their patterns of gene expression were altered, due to changes in concentrations of externally supplied substances. Later, as multicellular growth continued, the relationships of cell layers to each other shifted, producing concentration gradients of signaling molecules. These gradients emanated both from the external cell layer toward the inside and from internal cell layers to adjacent layers. In this scenario then, differentiation arose initially as a by-product of the changing patterns of gene expression and of the complex mixtures and changing concentrations of substances passing among layers. Subsequent selection would operate to stabilize the expression patterns in those cell layers whose phenotypes provide a fitness advantage to the organism.

Gene expression in fibroblasts and fibrosis: involvement in cardiac hypertrophy

Structural remodeling of the ventricular wall is a key determinant of clinical outcome in heart disease. Such remodeling involves the production and destruction of extracellular matrix proteins, cell proliferation and migration, and apoptotic and necrotic cell death. Cardiac fibroblasts are crucially involved in these processes, producing growth factors and cytokines that act as autocrine and paracrine factors, as well as extracellular matrix proteins and proteinases. Recent studies have shown that the interactions between cardiac fibroblasts and cardiomyocytes are essential for the progression of cardiac remodeling. This review addresses the functional role played by cardiac fibroblasts and the molecular mechanisms that govern their activity during cardiac hypertrophy and remodeling. A particular focus is the recent progress toward our understanding of the transcriptional regulatory mechanisms involved.

Genetic "code": representations and dynamical models of genetic components and networks

Dynamical modeling of biological systems is becoming increasingly widespread as people attempt to grasp biological phenomena in their full complexity and make sense of an accelerating stream of experimental data. We review a number of recent modeling studies that focus on systems specifically involving gene expression and regulation. These systems include bacterial metabolic operons and phase-variable piliation, bacteriophages T7 and lambda, and interacting networks of eukaryotic developmental genes. A wide range of conceptual and mathematical representations of genetic components and phenomena appears in these works. We discuss these representations in depth and give an overview of the tools currently available for creating and exploring dynamical models. We argue that for modeling to realize its full potential as a mainstream biological research technique the tools must become more general and flexible, and formal, standardized representations of biological knowledge and data must be developed.

Modeling transcriptional regulatory networks

Developmental processes in complex animals are directed by a hardwired genomic regulatory code, the ultimate function of which is to set up a progression of transcriptional regulatory states in space and time. The code specifies the gene regulatory networks (GRNs) that underlie all major developmental events. Models of GRNs are required for analysis, for experimental manipulation and, most fundamentally, for comprehension of how GRNs work. To model GRNs requires knowledge of both their overall structure, which depends upon linkage amongst regulatory genes, and the modular building blocks of which GRNs are heirarchically constructed. The building blocks consist of basic transcriptional control processes executed by one or a few functionally linked genes. We show how the functions of several such building blocks can be considered in mathematical terms, and discuss resolution of GRNs by both "top down" and "bottom up" approaches.

Splitting pairs: the diverging fates of duplicated genes

Many genes are members of large families that have arisen during evolution through gene duplication events. Our increasing understanding of gene organization at the scale of whole genomes is revealing further evidence for the extensive retention of genes that arise during duplication events of various types. Duplication is thought to be an important means of providing a substrate on which evolution can work. An understanding of gene duplication and its resolution is crucial for revealing mechanisms of genetic redundancy. Here, we consider both the theoretical framework and the experimental evidence to explain the preservation of duplicated genes.

Physical constraints and functional characteristics of transcription factor-DNA interaction

We study theoretical "design principles" for transcription factor (TF)-DNA interaction in bacteria, focusing particularly on the statistical interaction of the TFs with the genomic background (i.e., the genome without the target sites). We introduce and motivate the concept of programmability, i.e., the ability to set the threshold concentration for TF binding over a wide range merely by mutating the binding sequence of a target site. This functional demand, together with physical constraints arising from the thermodynamics and kinetics of TF-DNA interaction, leads us to a narrow range of "optimal" interaction parameters. We find that this parameter set agrees well with experimental data for the interaction parameters of a few exemplary prokaryotic TFs, which indicates that TF-DNA interaction is indeed programmable. We suggest further experiments to test whether this is a general feature for a large class of TFs.