Master Equation: Biological Applications and Thermodynamic Description

It is well known that many realistic mathematical models of biological systems, such as cell growth, cellular development and differentiation, gene expression, gene regulatory networks, enzyme cascades, synaptic plasticity, aging and population growth need to include stochasticity. These systems are not isolated, but rather subject to intrinsic and extrinsic fluctuations, which leads to a quasi equilibrium state (homeostasis). The natural framework is provided by Markov processes and the Master equation (ME) describes the temporal evolution of the probability of each state, specified by the number of units of each species. The ME is a relevant tool for modeling realistic biological systems and allow also to explore the behavior of open systems. These systems may exhibit not only the classical thermodynamic equilibrium states but also the nonequilibrium steady states (NESS). This thesis deals with biological problems that can be treat with the Master equation and also with its thermodynamic consequences. It is organized into six chapters with four new scientific works, which are grouped in two parts: (1) Biological applications of the Master equation: deals with the stochastic properties of a toggle switch, involving a protein compound and a miRNA cluster, known to control the eukaryotic cell cycle and possibly involved in oncogenesis and with the propose of a one parameter family of master equations for the evolution of a population having the logistic equation as mean field limit. (2) Nonequilibrium thermodynamics in terms of the Master equation: where we study the dynamical role of chemical fluxes that characterize the NESS of a chemical network and we propose a one parameter parametrization of BCM learning, that was originally proposed to describe plasticity processes, to study the differences between systems in DB and NESS.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Differential expression of hypoxia pathway genes in honey bee (Apis mellifera L.) caste development

Diphenism in social bees is essentially contingent on nutrient-induced cellular and systemic physiological responses resulting in divergent gene expression patterns. Analyses of juvenile hormone (JH) titers and functional genomics assays of the insulin-insulin-like signaling (IIS) pathway and its associated branch, target-of-rapamycin (TOR), revealed systemic responses underlying honey bee (Apis mellifera) caste development. Nevertheless, little attention has been paid to cellular metabolic responses. Following up earlier investigations showing major caste differences in oxidative metabolism and mitochondrial physiology, we herein identified honey bee homologs of hypoxia signaling factors, HIF alpha/Sima, HIF beta/Tango and PHD/Fatiga and we investigated their transcript levels throughout critical stages of larval development. Amsima, Amtango and Amfatiga showed correlated transcriptional activity, with two peaks of occurring in both queens and workers, the first one shortly after the last larval molt and the second during the cocoon-spinning phase. Transcript levels for the three genes were consistently higher in workers. As there is no evidence for major microenvironmental differences in oxygen levels within the brood nest area, this appears to be an inherent caste character. Quantitative PCR analyses on worker brain, ovary, and leg imaginal discs showed that these tissues differ in transcript levels. Being a highly conserved pathway and linked to IIS/TOR, the hypoxia gene expression pattern seen in honey bee larvae denotes that the hypoxia pathway has undergone a transformation, at least during larval development, from a response to environmental oxygen concentrations to an endogenous regulatory factor in the diphenic development of honey bee larvae. (C) 2010 Elsevier Ltd. All rights reserved.

Transcriptional changes in the nuc-2A mutant strain of Neurospora crassa cultivated under conditions of phosphate shortage

The molecular mechanism that controls the response to phosphate shortage in Neurospora crassa involves four regulatory genes - nuc-2, preg, pgov, and nuc-1. Phosphate shortage is sensed by the nuc-2 gene, the product of which inhibits the functioning of the PREG-PGOV complex. This allows the translocation of the transcriptional factor NUC-1 into the nucleus, which activates the transcription of phosphate-repressible phosphatases. The nuc-2A mutant strain of N. crassa carries a loss-of-function mutation in the nuc-2 gene, which encodes an ankyrin-like repeat protein. In this study, we identified transcripts that are downregutated in the nuc-2A mutant strain. Functional grouping of these expressed sequence tags allowed the identification of genes that play essential roles in different cellular processes such as transport, transcriptional regulation, signal transduction, metabolism, protein synthesis, protein fate, and development. These results reveal novel aspects of the phosphorus-sensing network in N. crassa. (C) 2009 Elsevier GmbH. All rights reserved.

Mouse interleukin-2 receptor alpha gene expression. Delimitation of cis-acting regulatory elements in transgenic mice and by mapping of DNase-I hypersensitive sites.

The alpha chain of the interleukin-2 receptor (IL-2R alpha) is a key regulator of lymphocyte proliferation. To analyze the mechanisms controlling its expression in normal cells, we used the 5'-flanking region (base pairs -2539/+93) of the mouse gene to drive chloramphenicol acetyltransferase expression in four transgenic mouse lines. Constitutive transgene activity was restricted to lymphoid organs. In mature T lymphocytes, transgene and endogenous IL-2R alpha gene expression was stimulated by concanavalin A and up-regulated by IL-2 with very similar kinetics. In thymic T cell precursors, IL-1 and IL-2 cooperatively induced transgene and IL-2R alpha gene expression. These results show that regulation of the endogenous IL-2R alpha gene occurs mainly at the transcriptional level. They demonstrate that cis-acting elements in the 5'-flanking region present in the transgene confer correct tissue specificity and inducible expression in mature T cells and their precursors in response to antigen, IL-1, and IL-2. In a complementary approach, we screened the 5' end of the endogenous IL-2R alpha gene for DNase-I hypersensitive sites. We found three lymphocyte specific DNase-I hypersensitive sites. Two, at -0.05 and -5.3 kilobase pairs, are present in resting T cells. A third site appears at -1.35 kilobase pairs in activated T cells. It co-localizes with IL-2-responsive elements identified by transient transfection experiments.

Characterization of the murine high Km glucose transporter GLUT2 gene and its transcriptional regulation by glucose in a differentiated insulin-secreting cell line.

In pancreatic beta-cells, the high Km glucose transporter GLUT2 catalyzes the first step in glucose-induced insulin secretion by glucose uptake. Expression of the transporter has been reported to be modulated by glucose either at the protein or mRNA levels. In this study we used the differentiated insulinoma cell line INS-1 which expresses high levels of GLUT2 and show that the expression of GLUT2 is regulated by glucose at the transcriptional level. By run-on transcription assays we showed that glucose induced GLUT2 gene transcription 3-4-fold in INS-1 cells which was paralleled by a 1.7-2.3-fold increase in cytoplasmic GLUT2 mRNA levels. To determine whether glucose regulatory sequences were present in the promoter region of GLUT2, we cloned and characterized a 1.4-kilobase region of mouse genomic DNA located 5' of the translation initiation site. By RNase protection assays and primer extension, we determined that multiple transcription initiation sites were present at positions -55, -64, and -115 from the first coding ATG and which were identified in liver, intestine, kidney, and beta-cells mRNAs. Plasmids were constructed with the mouse promoter region linked to the reporter gene chloramphenicol acetyltransferase (CAT), and transiently and stably transfected in the INS-1 cells. Glucose induced a concentration-dependent increase in CAT activity which reached a maximum of 3.6-fold at 20 mM glucose. Similar CAT constructs made of the human GLUT2 promoter region and the CAT gene displayed the same glucose-dependent increase in transcriptional activity when transfected into INS-1 cells. Comparison of the mouse and human promoter regions revealed sequence identity restricted to a few stretches of sequences which suggests that the glucose responsive element(s) may be conserved in these common sequences.

Comparative study and identification of potent eukaryotic transcriptional repressors in gene switch systems.

In mammalian cells, proper gene regulation is achieved by the complex interplay of transcription factors that activate or repress gene expression by binding to the regulatory regions of target promoters. While transcriptional activators have been extensively characterised and classified into functional groups, relatively little is known about the comparative strength and cell type-specificity of transcriptional repressors. Here, we have compared the ability of a series of eukaryotic repression domains to silence basal and activated transcription. A series of the most potent repression domains was further tested in the context of a gene therapy gene-switch system in various cell types. The results indicate that the analysed repression domains exert varying silencing activities in different promoter contexts. Furthermore, their potential for gene silencing varies also depending on the cellular context. When multimerised within one chimeric repressor protein, particular combinations of repressor domains were found to display synergistic repressing effects and efficient repression in a panel of cell lines. This approach thus allowed the identification of transcriptional repressors that are both potent and versatile in terms of cellular specificity as a basis for gene switch systems.

Implication of DNA methylation and PAX5 factor in the transcriptional regulation of hTERT, the human telomerase reverse transcriptase gene

RESUME La télomérase est une enzyme dite "d'immortalité" qui permet aux cellules de maintenir la longueur de leurs télomères, ce qui confère une capacité de réplication illimitée aux cellules reproductrices et cancéreuses. A l'inverse, les cellules somatiques normales, qui n'expriment pas la télomérase, ont une capacité de réplication limitée. La sous-unité catalytique de la télomérase, hTERT, est définie comme le facteur limitant l'activité télomérasique. Entre activateurs et répresseurs, le rôle de la méthylation de l'ADN et de l'acétylation des histones, de nombreux modèles ont été suggérés. La découverte de l'implication de CTCF dans la régulation transcriptionnelle de hTERT explique en partie le mécanisme de répression de la télomérase dans la plupart des cellules somatiques et sa réactivation dans les cellules tumorales. Dans les cellules télomérase-positives, l'activité inhibitrice de CTCF est bloquée par un mécanisme dépendent ou non de la méthylation. Dans la plupart des carcinomes, une hyperméthylation de la région 5' de hTERT bloque l'effet inhibiteur de CTCF, alors qu'une petite région hypométhylée permet un faible niveau de transcription du gène. Nous avons démontré que la protéine MBD2 se lie spécifiquement sur la région 5' méthylée de hTERT dans différentes lignées cellulaires et qu'elle est impliquée dans la répression partielle de la transcription de hTERT dans les cellules tumorales méthylées. Par contre, nous avons montré que dans les lymphocytes B normaux et néoplasiques, la régulation de hTERT est indépendante de la méthylation. Dans ces cellules, le facteur PAX5 se lie sur la région 5' de hTERT en aval du site d'initiation de la traduction (ATG). L'expression exogène de PAX5 dans les cellules télomérase-négatives active la transcription de hTERT, alors que la répression de PAX5 dans les cellules lymphomateuses inhibe la transcription du gène. PAX5 est donc directement impliqué dans l'activation de l'expression de hTERT dans les lymphocytes B exprimant la télomérase. Ces résultats révèlent des différences entre les niveaux de méthylation de hTERT dans les cellules de carcinomes et les lymphocytes B exprimant la télomérase. La méthylation de hTERT en tant que biomarqueur de cancer a été évaluée, puis appliquée à la détection de métastases. Nous avons ainsi montré que la méthylation de hTERT est positivement corrélée au diagnostic cytologique dans les liquides céphalorachidiens. Nos résultats conduisent à un modèle de régulation de hTERT, qui aide à comprendre comment la transcription de ce gène est régulée par CTCF, avec un mécanisme lié ou non à la méthylation du gène hTERT. La méthylation de hTERT s'est aussi révélée être un nouveau et prometteur biomarqueur de cancer. SUMMARY Human telomerase is an "immortalizing" enzyme that enables cells to maintain telomere length, allowing unlimited replicative capacity to reproductive and cancer cells. Conversely, normal somatic cells that do not express telomerase have a finite replicative capacity. The catalytic subunit of telomerase, hTERT, is defined as the limiting factor for telomerase activity. Between activators and repressors, and the role of DNA methylation and histone acetylation, an abundance of hTERT regulatory models have been suggested. The discovery of the implication of CTCF in the transcriptional regulation of hTERT in part explained the mechanism of silencing of telomerase in most somatic cells and its reactivation in neoplastic cells. In telomerase-positive cells, the inhibitory activity of CTCF is blocked by methylation-dependent and -independent mechanisms. In most carcinoma cells, hypermethylation of the hTERT 5' region has been shown to block the inhibitory effect of CTCF, while a short hypomethylated region allows a low transcription level of the gene. We have demonstrated that MBD2 protein specifically binds the methylated 5' region of hTERT in different cell lines and is therefore involved in the partial repression of hTERT transcription in methylated tumor cells. In contrast, we have shown that in normal and neoplastic B cells, hTERT regulation is methylation-independent. The PAX5 factor has been shown to bind to the hTERT 5'region downstream of the ATG translational start site. Ectopic expression of PAX5 in telomerase-negative cells or repression of PAX5 expression in B lymphoma cells respectively activated and repressed hTERT transcription. Thus, PAX5 is strongly implicated in hTERT expression activation in telomerase-positive B cells. These results reveal differences between the hTERT methylation patterns in telomerase-positive carcinoma cells and telomerase-positive normal B cells. The potential of hTERT methylation as a cancer biomarker was evaluated and applied to the detection of metastasis. We have shown that hTERT methylation correlates with the cytological diagnosis in cerebrospinal fluids. Our results suggest a model of hTERT gene regulation, which helps us to better understand how hTERT transcription is regulated by CTCF in methylation-dependant and independent mechanisms. Our data also indicate that hTERT methylation is a promising new cancer biomarker.

Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

Role of Sam68 in post-transcriptional gene regulation.

The STAR family of proteins links signaling pathways to various aspects of post-transcriptional regulation and processing of RNAs. Sam68 belongs to this class of heteronuclear ribonucleoprotein particle K (hnRNP K) homology (KH) single domain-containing family of RNA-binding proteins that also contains some domains predicted to bind critical components in signal transduction pathways. In response to phosphorylation and other post-transcriptional modifications, Sam68 has been shown to have the ability to link signal transduction pathways to downstream effects regulating RNA metabolism, including transcription, alternative splicing or RNA transport. In addition to its function as a docking protein in some signaling pathways, this prototypic STAR protein has been identified to have a nuclear localization and to take part in the formation of both nuclear and cytosolic multi-molecular complexes such as Sam68 nuclear bodies and stress granules. Coupling with other proteins and RNA targets, Sam68 may play a role in the regulation of differential expression and mRNA processing and translation according to internal and external signals, thus mediating important physiological functions, such as cell death, proliferation or cell differentiation.

The transcriptional repressor REST determines the cell-specific expression of the human MAPK8IP1 gene encoding IB1 (JIP-1).

Islet-brain 1 (IB1) is the human and rat homologue of JIP-1, a scaffold protein interacting with the c-Jun amino-terminal kinase (JNK). IB1 expression is mostly restricted to the endocrine pancreas and to the central nervous system. Herein, we explored the transcriptional mechanism responsible for this preferential islet and neuronal expression of IB1. A 731-bp fragment of the 5' regulatory region of the human MAPK8IP1 gene was isolated from a human BAC library and cloned upstream of a luciferase reporter gene. This construct drove high transcriptional activity in both insulin-secreting and neuron-like cells but not in unrelated cell lines. Sequence analysis of this promoter region revealed the presence of a neuron-restrictive silencer element (NRSE) known to bind repressor zinc finger protein REST. This factor is not expressed in insulin-secreting and neuron-like cells. By mobility shift assay, we confirmed that REST binds to the NRSE present in the IB1 promoter. Once transiently transfected in beta-cell lines, the expression vector encoding REST repressed IB1 transcriptional activity. The introduction of a mutated NRSE in the 5' regulating region of the IB1 gene abolished the repression activity driven by REST in insulin-secreting beta cells and relieved the low transcriptional activity of IB1 observed in unrelated cells. Moreover, transfection in non-beta and nonneuronal cell lines of an expression vector encoding REST lacking its transcriptional repression domain relieved IB1 promoter activity. Last, the REST-mediated repression of IB1 could be abolished by trichostatin A, indicating that deacetylase activity is required to allow REST repression. Taken together, these data establish a critical role for REST in the control of the tissue-specific expression of the human IB1 gene.