基于微阵列数据构建基因调控网络

随着微阵列数据的快速增长, 微阵列基因表达数据日益成为生物信息学研究的重要数据源。利用微阵列基因表达数据构建基因调控网络也成为一个研究热点。通过构建基因调控网络, 可以解读复杂的调控关系, 发现细胞内的调控模式, 并进而在系统尺度上理解生物学进程。近年来, 人们引入了多种算法来利用基因芯片数据构建基因调控网络。文章回顾了这些算法的发展历史, 尤其是其在理论和方法上的改进, 给出了一些相关的软件平台, 并预测了该领域可能的发展趋势。     

Modeling of Hysteresis in Gene Regulatory Networks

Hysteresis, observed in many gene regulatory networks, has a pivotal impact on biological systems, which enhances the robustness of cell functions. In this paper, a general model is proposed to describe the hysteretic gene regulatory network by combining the hysteresis component and the transient dynamics. The Bouc-Wen hysteresis model is modified to describe the hysteresis component in the mammalian gene regulatory networks. Rigorous mathematical analysis on the dynamical properties of the model is presented to ensure the bounded-input-bounded-output (BIBO) stability and demonstrates that the original Bouc-Wen model can only generate a clockwise hysteresis loop while the modified model can describe both clockwise and counter clockwise hysteresis loops. Simulation studies have shown that the hysteresis loops from our model are consistent with the experimental observations in three mammalian gene regulatory networks and two E.coli gene regulatory networks, which demonstrate the ability and accuracy of the mathematical model to emulate natural gene expression behavior with hysteresis. A comparison study has also been conducted to show that this model fits the experiment data significantly better than previous ones in the literature. The successful modeling of the hysteresis in all the five hysteretic gene regulatory networks suggests that the new model has the potential to be a unified framework for modeling hysteresis in gene regulatory networks and provide better understanding of the general mechanism that drives the hysteretic function.     

Efficient reverse-engineering of a developmental gene regulatory network

Understanding the complex regulatory networks underlying development and evolution of multi-cellular organisms is a major problem in biology. Computational models can be used as tools to extract the regulatory structure and dynamics of such networks from gene expression data. This approach is called reverse engineering. It has been successfully applied to many gene networks in various biological systems. However, to reconstitute the structure and non-linear dynamics of a developmental gene network in its spatial context remains a considerable challenge. Here, we address this challenge using a case study: the gap gene network involved in segment determination during early development of Drosophila melanogaster. A major problem for reverse-engineering pattern-forming networks is the significant amount of time and effort required to acquire and quantify spatial gene expression data. We have developed a simplified data processing pipeline that considerably increases the throughput of the method, but results in data of reduced accuracy compared to those previously used for gap gene network inference. We demonstrate that we can infer the correct network structure using our reduced data set, and investigate minimal data requirements for successful reverse engineering. Our results show that timing and position of expression domain boundaries are the crucial features for determining regulatory network structure from data, while it is less important to precisely measure expression levels. Based on this, we define minimal data requirements for gap gene network inference. Our results demonstrate the feasibility of reverse-engineering with much reduced experimental effort. This enables more widespread use of the method in different developmental contexts and organisms. Such systematic application of data-driven models to real-world networks has enormous potential. Only the quantitative investigation of a large number of developmental gene regulatory networks will allow us to discover whether there are rules or regularities governing development and evolution of complex multi-cellular organisms.     

基因表达时间延迟调控关系识别软件ITdGR

基因表达调控网络的深入研究有利于分子药物靶标的发现以及推新药的研发,是未来生物医学研究的重要内容。针对基因表达调控的时间延迟问题,我们初步设计开发了一套基于基因表达谱数据识别基因表达时间延迟调控关系的软件ITdGR(Identification of Time-delayed Gene Regulations)。并已经成功地将该软件应用于酿酒酵母细胞周期的基因表达谱数据中,识别出的调控关系与已有的知识相符。该软件为基因调控网络重构以及基因表达动态研究提供了一个方便和快捷的工具。     

Innexin genes and gap junction proteins in the locust frontal ganglion

Gap junctions (GJs) belong to one of the most conserved cellular structures in multicellular organisms. They probably serve similar functions in all Metazoa, providing one of the most common forms of intercellular communication. GJs are widely distributed in embryonic cells and tissues and have been attributed an important role in development, modulating cell growth and differentiation. These channels have been also implicated in mediating electrical synaptic signaling; Coupling through GJs is now accepted as a major pathway that supports network behavior and contributes to physiological rhythms. Here we focus on the physiology and molecular biology of GJs in a recently established model for the study of rhythm-generating networks and their role in behavior: the frontal ganglion (FG) of the desert locust, Schistocerca gregaria. Four novel genes of the invertebrate GJs (innexin) gene family were found to be expressed in the FG: Sg-inx1, Sg-inx2, Sg-inx3 and Sg-inx4. Immunohistochemistry revealed that some of the neurons in the FG express at least one innexin protein, INX1. We also established the presence of functional gap junction proteins in the FG and demonstrated functional electrical coupling between the neurons in the FG. This study forms the basis for further investigation of the role of GJs in network development and behavior.     

Gradient Descent Optimization in Gene Regulatory Pathways

Background

Gene Regulatory Networks (GRNs) have become a major focus of interest in recent years. Elucidating the architecture and dynamics of large scale gene regulatory networks is an important goal in systems biology. The knowledge of the gene regulatory networks further gives insights about gene regulatory pathways. This information leads to many potential applications in medicine and molecular biology, examples of which are identification of metabolic pathways, complex genetic diseases, drug discovery and toxicology analysis. High-throughput technologies allow studying various aspects of gene regulatory networks on a genome-wide scale and we will discuss recent advances as well as limitations and future challenges for gene network modeling. Novel approaches are needed to both infer the causal genes and generate hypothesis on the underlying regulatory mechanisms.

Methodology

In the present article, we introduce a new method for identifying a set of optimal gene regulatory pathways by using structural equations as a tool for modeling gene regulatory networks. The method, first of all, generates data on reaction flows in a pathway. A set of constraints is formulated incorporating weighting coefficients. Finally the gene regulatory pathways are obtained through optimization of an objective function with respect to these weighting coefficients. The effectiveness of the present method is successfully tested on ten gene regulatory networks existing in the literature. A comparative study with the existing extreme pathway analysis also forms a part of this investigation. The results compare favorably with earlier experimental results. The validated pathways point to a combination of previously documented and novel findings.

Conclusions

We show that our method can correctly identify the causal genes and effectively output experimentally verified pathways. The present method has been successful in deriving the optimal regulatory pathways for all the regulatory networks considered. The biological significance and applicability of the optimal pathways has also been discussed. Finally the usefulness of the present method on genetic engineering is depicted with an example.     

From blastocyst to gastrula: gene regulatory networks of embryonic stem cells and early mouse embryogenesis

To date, many regulatory genes and signalling events coordinating mammalian development from blastocyst to gastrulation stages have been identified by mutational analyses and reverse-genetic approaches, typically on a gene-by-gene basis. More recent studies have applied bioinformatic approaches to generate regulatory network models of gene interactions on a genome-wide scale. Such models have provided insights into the gene networks regulating pluripotency in embryonic and epiblast stem cells, as well as cell-lineage determination in vivo. Here, we review how regulatory networks constructed for different stem cell types relate to corresponding networks in vivo and provide insights into understanding the molecular regulation of the blastocyst–gastrula transition.     

Synthetic circuits reveal how mechanisms of gene regulatory networks constrain evolution

Phenotypic variation is the raw material of adaptive Darwinian evolution. The phenotypic variation found in organismal development is biased towards certain phenotypes, but the molecular mechanisms behind such biases are still poorly understood. Gene regulatory networks have been proposed as one cause of constrained phenotypic variation. However, most pertinent evidence is theoretical rather than experimental. Here, we study evolutionary biases in two synthetic gene regulatory circuits expressed in Escherichia coli that produce a gene expression stripe—a pivotal pattern in embryonic development. The two parental circuits produce the same phenotype, but create it through different regulatory mechanisms. We show that mutations cause distinct novel phenotypes in the two networks and use a combination of experimental measurements, mathematical modelling and DNA sequencing to understand why mutations bring forth only some but not other novel gene expression phenotypes. Our results reveal that the regulatory mechanisms of networks restrict the possible phenotypic variation upon mutation. Consequently, seemingly equivalent networks can indeed be distinct in how they constrain the outcome of further evolution.     

Genetic and chromosomal control of development

The last decade has been characterized by extraordinary progress in studies of the molecular mechanisms underlying the gene control of development: complexes of genes with a hierarchic principle of functioning have been identified, evolutionary-conservative systems of genes have been studied that ensure the transmembrane regulatory signaling between cells and play a key role in morphogenesis, and a concept of gene networks that coordinate gene interaction was introduced. Note that temporal and tissue-specific parameters of gene expression are correctly realized only in terms of chromosomes and are largely determined by the gene position on a chromosome and in the interphase nucleus. The epigenetic inheritance of gene status in cell generations is realized at the chromosomal level alone due to the cellular or chromosomal memory. This ontogenetic memory is an immanent property of a chromosome and cis-regulation plays a key role in its maintenance.     

Mitochondrial GFA2 is required for synergid cell death in Arabidopsis

Little is known about the molecular processes that govern female gametophyte (FG) development and function, and few FG-expressed genes have been identified. We report the identification and phenotypic analysis of 31 new FG mutants in Arabidopsis. These mutants have defects throughout development, indicating that FG-expressed genes govern essentially every step of FG development. To identify genes involved in cell death during FG development, we analyzed this mutant collection for lines with cell death defects. From this analysis, we identified one mutant, gfa2, with a defect in synergid cell death. Additionally, the gfa2 mutant has a defect in fusion of the polar nuclei. We isolated the GFA2 gene and show that it encodes a J-domain-containing protein. Of the J-domain-containing proteins in Saccharomyces cerevisiae (budding yeast), GFA2 is most similar to Mdj1p, which functions as a chaperone in the mitochondrial matrix. GFA2 is targeted to mitochondria in Arabidopsis and partially complements a yeast mdj1 mutant, suggesting that GFA2 is the Arabidopsis ortholog of yeast Mdj1p. These data suggest a role for mitochondria in cell death in plants.     

Perturbation avalanches and criticality in gene regulatory networks

Boolean networks are simplified models of gene regulatory networks. We derive an approximation of the size distribution of perturbation avalanches in Boolean networks based on known results in the theory of branching processes. We show numerically that the approximation works well for different kinds of Boolean networks. It has been suggested that gene regulatory networks may be dynamically critical. To study this, as an application of the presented theory we present a novel method for estimating an order parameter from microarray data. According to the available data and our method, we find that gene regulatory networks appear to be stable and reside near the phase transition between order and chaos.     

Reconstructing gene regulatory networks: from random to scale-free connectivity

The manipulation of organisms using combinations of gene knockout, RNAi and drug interaction experiments can be used to reveal regulatory interactions between genes. Several algorithms have been proposed that try to reconstruct the underlying regulatory networks from gene expression data sets arising from such experiments. Often these approaches assume that each gene has approximately the same number of interactions within the network, and the methods rely on prior knowledge, or the investigator's best guess, of the average network connectivity. Recent evidence points to scale-free properties in biological networks, however, where network connectivity follows a power-law distribution. For scale-free networks, the average number of regulatory interactions per gene does not satisfactorily characterise the network. With this in mind, a new reverse engineering approach is introduced that does not require prior knowledge of network connectivity and its performance is compared with other published algorithms using simulated gene expression data with biologically relevant network structures. Because this new approach does not make any assumptions about the distribution of network connections, it is suitable for application to scale-free networks.     

Evolution of robustness to damage in artificial 3-dimensional development

GReaNs is an Artificial Life platform we have built to investigate the general principles that guide evolution of multicellular development and evolution of artificial gene regulatory networks. The embryos develop in GReaNs in a continuous 3-dimensional (3D) space with simple physics. The developmental trajectories are indirectly encoded in linear genomes. The genomes are not limited in size and determine the topology of gene regulatory networks that are not limited in the number of nodes. The expression of the genes is continuous and can be modified by adding environmental noise. In this paper we evolved development of structures with a specific shape (an ellipsoid) and asymmetrical pattering (a 3D pattern inspired by the French flag problem), and investigated emergence of the robustness to damage in development and the emergence of the robustness to noise. Our results indicate that both types of robustness are related, and that including noise during evolution promotes higher robustness to damage. Interestingly, we have observed that some evolved gene regulatory networks rely on noise for proper behaviour.