Inferring models of gene expression dynamics

We study the problem of identifying genetic networks in which expression dynamics are modeled by a differential equation that uses logical rules to specify time derivatives. We make three main contributions. First, we describe computationally efficient procedures for identifying the structure and dynamics of such networks from expression time series. Second, we derive predictions for the expected amount of data needed to identify randomly generated networks. Third, if expression values are available for only some of the genes, we show that the structure of the network for these "visible" genes can be identified and that the size and overall complexity of the network can be estimated. We validate these procedures and predictions using simulation experiments based on randomly generated networks with up to 30,000 genes and 17 distinct regulators per gene and on a network that models floral morphogenesis in Arabidopsis thaliana.     

State space modeling of yeast gene expression dynamics

Combined interaction of all the genes forms a central part of the functional system of a cell. Thus, especially the data-based modeling of the gene expression network is currently one of the main challenges in the field of systems biology. However, the problem is an extremely high-dimensional and complex one, so that normal identification methods are usually not applicable specially if aiming at dynamic models. We propose in this paper a subspace identification approach, which is well suited for high-dimensional system modeling and the presented modified version can also handle the underdetermined case with less data samples than variables (genes). The algorithm is applied to two public stress-response data sets collected from yeast Saccharomyces cerevisiae. The obtained dynamic state space model is tested by comparing the simulation results with the measured data. It is shown that the identified model can relatively well describe the dynamics of the general stress-related changes in the expression of the complete yeast genome. However, it seems inevitable that more precise modeling of the dynamics of the whole genome would require experiments especially designed for systemic modeling.     

Collective dynamics of gene expression in cell populations

The phenotypic state of the cell is commonly thought to be determined by the set of expressed genes. However, given the apparent complexity of genetic networks, it remains open what processes stabilize a particular phenotypic state. Moreover, it is not clear how unique is the mapping between the vector of expressed genes and the cell's phenotypic state. To gain insight on these issues, we study here the expression dynamics of metabolically essential genes in twin cell populations. We show that two yeast cell populations derived from a single steady-state mother population and exhibiting a similar growth phenotype in response to an environmental challenge, displayed diverse expression patterns of essential genes. The observed diversity in the mean expression between populations could not result from stochastic cell-to-cell variability, which would be averaged out in our large cell populations. Remarkably, within a population, sets of expressed genes exhibited coherent dynamics over many generations. Thus, the emerging gene expression patterns resulted from collective population dynamics. It suggests that in a wide range of biological contexts, gene expression reflects a self-organization process coupled to population-environment dynamics.     

Inferring gene expression dynamics via functional regression analysis

Background

The recent emergence of high-throughput automated image acquisition technologies has forever changed how cell biologists collect and analyze data. Historically, the interpretation of cellular phenotypes in different experimental conditions has been dependent upon the expert opinions of well-trained biologists. Such qualitative analysis is particularly effective in detecting subtle, but important, deviations in phenotypes. However, while the rapid and continuing development of automated microscope-based technologies now facilitates the acquisition of trillions of cells in thousands of diverse experimental conditions, such as in the context of RNA interference (RNAi) or small-molecule screens, the massive size of these datasets precludes human analysis. Thus, the development of automated methods which aim to identify novel and biological relevant phenotypes online is one of the major challenges in high-throughput image-based screening. Ideally, phenotype discovery methods should be designed to utilize prior/existing information and tackle three challenging tasks, i.e. restoring pre-defined biological meaningful phenotypes, differentiating novel phenotypes from known ones and clarifying novel phenotypes from each other. Arbitrarily extracted information causes biased analysis, while combining the complete existing datasets with each new image is intractable in high-throughput screens.

Results

Here we present the design and implementation of a novel and robust online phenotype discovery method with broad applicability that can be used in diverse experimental contexts, especially high-throughput RNAi screens. This method features phenotype modelling and iterative cluster merging using improved gap statistics. A Gaussian Mixture Model (GMM) is employed to estimate the distribution of each existing phenotype, and then used as reference distribution in gap statistics. This method is broadly applicable to a number of different types of image-based datasets derived from a wide spectrum of experimental conditions and is suitable to adaptively process new images which are continuously added to existing datasets. Validations were carried out on different dataset, including published RNAi screening using Drosophila embryos [Additional files 1, 2], dataset for cell cycle phase identification using HeLa cells [Additional files 1, 3, 4] and synthetic dataset using polygons, our methods tackled three aforementioned tasks effectively with an accuracy range of 85%–90%. When our method is implemented in the context of a Drosophila genome-scale RNAi image-based screening of cultured cells aimed to identifying the contribution of individual genes towards the regulation of cell-shape, it efficiently discovers meaningful new phenotypes and provides novel biological insight. We also propose a two-step procedure to modify the novelty detection method based on one-class SVM, so that it can be used to online phenotype discovery. In different conditions, we compared the SVM based method with our method using various datasets and our methods consistently outperformed SVM based method in at least two of three tasks by 2% to 5%. These results demonstrate that our methods can be used to better identify novel phenotypes in image-based datasets from a wide range of conditions and organisms.

Conclusion

We demonstrate that our method can detect various novel phenotypes effectively in complex datasets. Experiment results also validate that our method performs consistently under different order of image input, variation of starting conditions including the number and composition of existing phenotypes, and dataset from different screens. In our findings, the proposed method is suitable for online phenotype discovery in diverse high-throughput image-based genetic and chemical screens.     

Evaluating gene expression dynamics using pairwise RNA FISH data

Recently, a novel approach has been developed to study gene expression in single cells with high time resolution using RNA Fluorescent In Situ Hybridization (FISH). The technique allows individual mRNAs to be counted with high accuracy in wild-type cells, but requires cells to be fixed; thus, each cell provides only a "snapshot" of gene expression. Here we show how and when RNA FISH data on pairs of genes can be used to reconstruct real-time dynamics from a collection of such snapshots. Using maximum-likelihood parameter estimation on synthetically generated, noisy FISH data, we show that dynamical programs of gene expression, such as cycles (e.g., the cell cycle) or switches between discrete states, can be accurately reconstructed. In the limit that mRNAs are produced in short-lived bursts, binary thresholding of the FISH data provides a robust way of reconstructing dynamics. In this regime, prior knowledge of the type of dynamics--cycle versus switch--is generally required and additional constraints, e.g., from triplet FISH measurements, may also be needed to fully constrain all parameters. As a demonstration, we apply the thresholding method to RNA FISH data obtained from single, unsynchronized cells of Saccharomyces cerevisiae. Our results support the existence of metabolic cycles and provide an estimate of global gene-expression noise. The approach to FISH data presented here can be applied in general to reconstruct dynamics from snapshots of pairs of correlated quantities including, for example, protein concentrations obtained from immunofluorescence assays.     

The complexity of gene expression dynamics revealed by permutation entropy

Background

High complexity is considered a hallmark of living systems. Here we investigate the complexity of temporal gene expression patterns using the concept of Permutation Entropy (PE) first introduced in dynamical systems theory. The analysis of gene expression data has so far focused primarily on the identification of differentially expressed genes, or on the elucidation of pathway and regulatory relationships. We aim to study gene expression time series data from the viewpoint of complexity.

Results

Applying the PE complexity metric to abiotic stress response time series data in Arabidopsis thaliana, genes involved in stress response and signaling were found to be associated with the highest complexity not only under stress, but surprisingly, also under reference, non-stress conditions. Genes with house-keeping functions exhibited lower PE complexity. Compared to reference conditions, the PE of temporal gene expression patterns generally increased upon stress exposure. High-complexity genes were found to have longer upstream intergenic regions and more cis-regulatory motifs in their promoter regions indicative of a more complex regulatory apparatus needed to orchestrate their expression, and to be associated with higher correlation network connectivity degree. Arabidopsis genes also present in other plant species were observed to exhibit decreased PE complexity compared to Arabidopsis specific genes.

Conclusions

We show that Permutation Entropy is a simple yet robust and powerful approach to identify temporal gene expression profiles of varying complexity that is equally applicable to other types of molecular profile data.     

Geometry and dynamics of activity-dependent homeostatic regulation in neurons

To maintain activity in a functional range, neurons constantly adjust membrane excitability to changing intra- and extracellular conditions. Such activity-dependent homeostatic regulation (ADHR) is critical for normal processing of the nervous system and avoiding pathological conditions. Here, we posed a homeostatic regulation problem for the classical Morris-Lecar (ML) model. The problem was motivated by the phenomenon of the functional recovery of stomatogastric neurons in crustaceans in the absence of neuromodulation. In our study, the regulation of the ionic conductances in the ML model depended on the calcium current or the intracellular calcium concentration. We found an asymptotic solution to the problem under the assumption of slow regulation. The solution provides a full account of the regulation in the case of correlated or anticorrelated changes of the maximal conductances of the calcium and potassium currents. In particular, the solution shows how the target and parameters of the regulation determine which perturbations of the conductances can be compensated by the ADHR. In some cases, the sets of compensated initial perturbations are not convex. On the basis of our analysis we formulated specific questions for subsequent experimental and theoretical studies of ADHR.     

Brain corticosteroid receptor gene expression and neuroendocrine dynamics during aging

The present study examined the stress responsiveness of the hypothalamic-pituitary-adrenal axis in relation to the properties of corticosteroid receptors in the brain and pituitary in old (30 months) and young (3 months) male Brown Norway rats. The data demonstrate that circulating ACTH rather than the corticosteroid plasma level was elevated under basal conditions and following stress. Furthermore, a reduction of mineralocorticoid receptor (MR) number in the hippocampus and of glucocorticoid receptor (GR) number in the hypothalamus and the pituitary correspond to increased neuroendocrine responsiveness and negative feedback following stress. The changes in receptor binding do not parallel the changes in the amount of MR and GR mRNA measured with in situ hybridization. This suggests that the processing rather than the receptor gene expression is affected in senescence.     

Visualization of the dynamics of gene expression in the living mouse

Reporter genes can monitor the status and activity of recombinant genomes in a diverse array of organisms, from bacteria and yeast to plants and animals. We have combined luciferase reporter genes with a conditional gene expression system based on regulatory elements from the lac operon of Escherichia coli to visualize the dynamics of gene expression in realtime in the living mouse. Using this technology, we have determined the rate of gene induction and repression, the level of target gene activity in response to different doses of inducer, and the schedule of induction during early embryogenesis of both the endogenous and the experimentally manipulated programs of mammalian gene expression associated with the HD/Hdh locus. The combination of in vivo imaging and lac regulation is a powerful tool for generating conditional transgenic mice that can be screened rapidly for optimal regulation and expression patterns, and for monitoring the induction and repression of regulated genes noninvasively in the living animal.     

Evolution of robustness to noise and mutation in gene expression dynamics

Phenotype of biological systems needs to be robust against mutation in order to sustain themselves between generations. On the other hand, phenotype of an individual also needs to be robust against fluctuations of both internal and external origins that are encountered during growth and development. Is there a relationship between these two types of robustness, one during a single generation and the other during evolution? Could stochasticity in gene expression have any relevance to the evolution of these types of robustness? Robustness can be defined by the sharpness of the distribution of phenotype; the variance of phenotype distribution due to genetic variation gives a measure of 'genetic robustness', while that of isogenic individuals gives a measure of 'developmental robustness'. Through simulations of a simple stochastic gene expression network that undergoes mutation and selection, we show that in order for the network to acquire both types of robustness, the phenotypic variance induced by mutations must be smaller than that observed in an isogenic population. As the latter originates from noise in gene expression, this signifies that the genetic robustness evolves only when the noise strength in gene expression is larger than some threshold. In such a case, the two variances decrease throughout the evolutionary time course, indicating increase in robustness. The results reveal how noise that cells encounter during growth and development shapes networks' robustness to stochasticity in gene expression, which in turn shapes networks' robustness to mutation. The necessary condition for evolution of robustness, as well as the relationship between genetic and developmental robustness, is derived quantitatively through the variance of phenotypic fluctuations, which are directly measurable experimentally.     

How to cluster gene expression dynamics in response to environmental signals

Organisms usually cope with change in the environment by altering the dynamic trajectory of gene expression to adjust the complement of active proteins. The identification of particular sets of genes whose expression is adaptive in response to environmental changes helps to understand the mechanistic base of gene-environment interactions essential for organismic development. We describe a computational framework for clustering the dynamics of gene expression in distinct environments through Gaussian mixture fitting to the expression data measured at a set of discrete time points. We outline a number of quantitative testable hypotheses about the patterns of dynamic gene expression in changing environments and gene-environment interactions causing developmental differentiation. The future directions of gene clustering in terms of incorporations of the latest biological discoveries and statistical innovations are discussed. We provide a set of computational tools that are applicable to modeling and analysis of dynamic gene expression data measured in multiple environments.     

Cell dynamics and gene expression control in tissue homeostasis and development

During tissue and organ development and maintenance, the dynamic regulation of cellular proliferation and differentiation allows cells to build highly elaborate structures. The development of the vertebrate retina or the maintenance of adult intestinal crypts, for instance, involves the arrangement of newly created cells with different phenotypes, the proportions of which need to be tightly controlled. While some of the basic principles underlying these processes developing and maintaining these organs are known, much remains to be learnt from how cells encode the necessary information and use it to attain those complex but reproducible arrangements. Here, we review the current knowledge on the principles underlying cell population dynamics during tissue development and homeostasis. In particular, we discuss how stochastic fate assignment, cell division, feedback control and cellular transition states interact during organ and tissue development and maintenance in multicellular organisms. We propose a framework, involving the existence of a transition state in which cells are more susceptible to signals that can affect their gene expression state and influence their cell fate decisions. This framework, which also applies to systems much more amenable to quantitative analysis like differentiating embryonic stem cells, links gene expression programmes with cell population dynamics.     

Measuring single-cell gene expression dynamics in bacteria using fluorescence time-lapse microscopy

Quantitative single-cell time-lapse microscopy is a powerful method for analyzing gene circuit dynamics and heterogeneous cell behavior. We describe the application of this method to imaging bacteria by using an automated microscopy system. This protocol has been used to analyze sporulation and competence differentiation in Bacillus subtilis, and to quantify gene regulation and its fluctuations in individual Escherichia coli cells. The protocol involves seeding and growing bacteria on small agarose pads and imaging the resulting microcolonies. Images are then reviewed and analyzed using our laboratory's custom MATLAB analysis code, which segments and tracks cells in a frame-to-frame method. This process yields quantitative expression data on cell lineages, which can illustrate dynamic expression profiles and facilitate mathematical models of gene circuits. With fast-growing bacteria, such as E. coli or B. subtilis, image acquisition can be completed in 1 d, with an additional 1-2 d for progressing through the analysis procedure.