Transcriptional regulation by the numbers: models

The expression of genes is regularly characterized with respect to how much, how fast, when and where. Such quantitative data demands quantitative models. Thermodynamic models are based on the assumption that the level of gene expression is proportional to the equilibrium probability that RNA polymerase (RNAP) is bound to the promoter of interest. Statistical mechanics provides a framework for computing these probabilities. Within this framework, interactions of activators, repressors, helper molecules and RNAP are described by a single function, the "regulation factor". This analysis culminates in an expression for the probability of RNA polymerase binding at the promoter of interest as a function of the number of regulatory proteins in the cell.     

Mining gene expression data using a novel approach based on hidden Markov models

In this work we have developed a new framework for microarray gene expression data analysis. This framework is based on hidden Markov models. We have benchmarked the performance of this probability model-based clustering algorithm on several gene expression datasets for which external evaluation criteria were available. The results showed that this approach could produce clusters of quality comparable to two prevalent clustering algorithms, but with the major advantage of determining the number of clusters. We have also applied this algorithm to analyze published data of yeast cell cycle gene expression and found it able to successfully dig out biologically meaningful gene groups. In addition, this algorithm can also find correlation between different functional groups and distinguish between function genes and regulation genes, which is helpful to construct a network describing particular biological associations. Currently, this method is limited to time series data. Supplementary materials are available at http://www.bioinfo.tsinghua.edu.cn/~rich/hmmgep_supp/.     

Amplitude Metrics for Cellular Circadian Bioluminescence Reporters

Bioluminescence rhythms from cellular reporters have become the most common method used to quantify oscillations in circadian gene expression. These experimental systems can reveal phase and amplitude change resulting from circadian disturbances, and can be used in conjunction with mathematical models to lend further insight into the mechanistic basis of clock amplitude regulation. However, bioluminescence experiments track the mean output from thousands of noisy, uncoupled oscillators, obscuring the direct effect of a given stimulus on the genetic regulatory network. In many cases, it is unclear whether changes in amplitude are due to individual changes in gene expression level or to a change in coherence of the population. Although such systems can be modeled using explicit stochastic simulations, these models are computationally cumbersome and limit analytical insight into the mechanisms of amplitude change. We therefore develop theoretical and computational tools to approximate the mean expression level in large populations of noninteracting oscillators, and further define computationally efficient amplitude response calculations to describe phase-dependent amplitude change. At the single-cell level, a mechanistic nonlinear ordinary differential equation model is used to calculate the transient response of each cell to a perturbation, whereas population-level dynamics are captured by coupling this detailed model to a phase density function. Our analysis reveals that amplitude changes mediated at either the individual-cell or the population level can be distinguished in tissue-level bioluminescence data without the need for single-cell measurements. We demonstrate the effectiveness of the method by modeling experimental bioluminescence profiles of light-sensitive fibroblasts, reconciling the conclusions of two seemingly contradictory studies. This modeling framework allows a direct comparison between in vitro bioluminescence experiments and in silico ordinary differential equation models, and will lead to a better quantitative understanding of the factors that affect clock amplitude.     

Mathematical modelling in the post-genome era: understanding genome expression and regulation--a system theoretic approach

This paper introduces a mathematical framework for modelling genome expression and regulation. Starting with a philosophical foundation, causation is identified as the principle of explanation of change in the realm of matter. Causation is, therefore, a relationship, not between components, but between changes of states of a system. We subsequently view genome expression (formerly known as 'gene expression') as a dynamic process and model aspects of it as dynamic systems using methodologies developed within the areas of systems and control theory. We begin with the possibly most abstract but general formulation in the setting of category theory. The class of models realised are state-space models, input--output models, autoregressive models or automata. We find that a number of proposed 'gene network' models are, therefore, included in the framework presented here. The conceptual framework that integrates all of these models defines a dynamic system as a family of expression profiles. It becomes apparent that the concept of a 'gene' is less appropriate when considering mathematical models of genome expression and regulation. The main claim of this paper is that we should treat (model) the organisation and regulation of genetic pathways as what they are: dynamic systems. Microarray technology allows us to generate large sets of time series data and is, therefore, discussed with regard to its use in mathematical modelling of gene expression and regulation.     

Mathematical approaches to differentiation and gene regulation

We consider some mathematical issues raised by the modelling of gene networks. The expression of genes is governed by a complex set of regulations, which is often described symbolically by interaction graphs. These are finite oriented graphs where vertices are the genes involved in the biological system of interest and arrows describe their interactions: a positive (resp. negative) arrow from a gene to another represents an activation (resp. inhibition) of the expression of the latter gene by some product of the former. Once such an interaction graph has been established, there remains the difficult task to decide which dynamical properties of the gene network can be inferred from it, in the absence of precise quantitative data about their regulation. There mathematical tools, among others, can be of some help. In this paper we discuss a rule proposed by Thomas according to which the possibility for the network to have several stationary states implies the existence of a positive circuit in the corresponding interaction graph. We prove that, when properly formulated in rigorous terms, this rule becomes a theorem valid for several different types of formal models of gene networks. This result is already known for models of differential [C. Soulé, Graphic requirements for multistationarity, ComPlexUs 1 (2003) 123-133] or Boolean [E. Rémy, P. Ruet, D. Thieffry, Graphic requirements for multistability and attractive cycles in a boolean dynamical framework, 2005, Preprint] type. We show here that a stronger version of it holds in the differential setup when the decay of protein concentrations is taken into account. This allows us to verify also the validity of Thomas' rule in the context of piecewise-linear models. We then discuss open problems.     

Transient expression of foreign DNA during embryonic and larval development of the medaka fish (Oryzias latipes)

Summary Species of small fish are becoming useful tools for studies on vertebrate development. We have investigated the developing embryo of the Japanese medaka for its application as a transient expression system for the in vivo analysis of gene regulation and function. The temporal and spatial expression patterns of bacterial chloramphenicol acetyltransferase and galactosidase reporter genes injected in supercoiled plasmid form into the cytoplasm of one cell of the two-cell stage embryo was promoter-specific. The transient expression was found to be mosaic within the tissue and organs reflecting the unequal distribution of extrachromosomal foreign DNA and the intensive cell mixing movements that occur in fish embryogenesis. The expression data are consistent with data on DNA fate. Foreign DNA persisted during embryogenesis and was still detectable in some 3- and 9-month-old adult fish; it was found in high molecular weight form as well as in circular plasmid conformations. The DNA was replicated during early and late embryogenesis. Our data indicate that the developing medaka embryo is a powerful in vivo assay system for studies of gene regulation and function.This work contains part of the PhD thesis of C. Winkler     

Small noncoding RNAs controlling pathogenesis

Infectious diseases are a leading cause of mortality worldwide. A major challenge in achieving their eradication is a better understanding of bacterial pathogenesis processes. The recent discovery of small noncoding RNAs (sRNAs) as modulators of gene expression in response to environmental cues has brought a new insight into bacterial regulation. sRNAs coordinate complex networks of stress adaptation and virulence gene expression. sRNAs generally ensure such a regulation by pairing to mRNAs of effector and/or regulatory genes, or by binding to proteins. An updated view on bacterial models responsible for important infections illustrates the key role of sRNAs in the control of pathogenesis.     

From classical genetics to quantitative genetics to systems biology: modeling epistasis

Gene expression data has been used in lieu of phenotype in both classical and quantitative genetic settings. These two disciplines have separate approaches to measuring and interpreting epistasis, which is the interaction between alleles at different loci. We propose a framework for estimating and interpreting epistasis from a classical experiment that combines the strengths of each approach. A regression analysis step accommodates the quantitative nature of expression measurements by estimating the effect of gene deletions plus any interaction. Effects are selected by significance such that a reduced model describes each expression trait. We show how the resulting models correspond to specific hierarchical relationships between two regulator genes and a target gene. These relationships are the basic units of genetic pathways and genomic system diagrams. Our approach can be extended to analyze data from a variety of experiments, multiple loci, and multiple environments.