Gene Regulatory Interactions Research Articles

Beyond the French Flag Model: Exploiting Spatial and Gene Regulatory Interactions for Positional Information

A crucial step in the early development of multicellular organisms involves the establishment of spatial patterns of gene expression which later direct proliferating cells to take on different cell fates. These patterns enable the cells to infer their global position within a tissue or an organism by reading out local gene expression levels. The patterning system is thus said to encode positional information, a concept that was formalized recently in the framework of information theory. Here we introduce a toy model of patterning in one spatial dimension, which can be seen as an extension of Wolpert’s paradigmatic “French Flag” model, to patterning by several interacting, spatially coupled genes subject to intrinsic and extrinsic noise. Our model, a variant of an Ising spin system, allows us to systematically explore expression patterns that optimally encode positional information. We find that optimal patterning systems use positional cues, as in the French Flag model, together with gene-gene interactions to generate combinatorial codes for position which we call “Counter” patterns. Counter patterns can also be stabilized against noise and variations in system size or morphogen dosage by longer-range spatial interactions of the type invoked in the Turing model. The simple setup proposed here qualitatively captures many of the experimentally observed properties of biological patterning systems and allows them to be studied in a single, theoretically consistent framework.

Studying Dynamic Features in Myocardial Infarction Progression by Integrating miRNA-Transcription Factor Co-Regulatory Networks and Time-Series RNA Expression Data from Peripheral Blood Mononuclear Cells

Myocardial infarction (MI) is a serious heart disease and a leading cause of mortality and morbidity worldwide. Although some molecules (genes, miRNAs and transcription factors (TFs)) associated with MI have been studied in a specific pathological context, their dynamic characteristics in gene expressions, biological functions and regulatory interactions in MI progression have not been fully elucidated to date. In the current study, we analyzed time-series RNA expression data from peripheral blood mononuclear cells. We observed that significantly differentially expressed genes were sharply up- or down-regulated in the acute phase of MI, and then changed slowly until the chronic phase. Biological functions involved at each stage of MI were identified. Additionally, dynamic miRNA–TF co-regulatory networks were constructed based on the significantly differentially expressed genes and miRNA–TF co-regulatory motifs, and the dynamic interplay of miRNAs, TFs and target genes were investigated. Finally, a new panel of candidate diagnostic biomarkers (STAT3 and ICAM1) was identified to have discriminatory capability for patients with or without MI, especially the patients with or without recurrent events. The results of the present study not only shed new light on the understanding underlying regulatory mechanisms involved in MI progression, but also contribute to the discovery of true diagnostic biomarkers for MI.

TRaCE+: Ensemble inference of gene regulatory networks from transcriptional expression profiles of gene knock-out experiments.

BackgroundThe inference of gene regulatory networks (GRNs) from transcriptional expression profiles is challenging, predominantly due to its underdetermined nature. One important consequence of underdetermination is the existence of many possible solutions to this inference. Our previously proposed ensemble inference algorithm TRaCE addressed this issue by inferring an ensemble of network directed graphs (digraphs) using differential gene expressions from gene knock-out (KO) experiments. However, TRaCE could not deal with the mode of the transcriptional regulations (activation or repression), an important feature of GRNs.ResultsIn this work, we developed a new algorithm called TRaCE+ for the inference of an ensemble of signed GRN digraphs from transcriptional expression data of gene KO experiments. The sign of the edges indicates whether the regulation is an activation (positive) or a repression (negative). TRaCE+ generates the upper and lower bounds of the ensemble, which define uncertain regulatory interactions that could not be verified by the data. As demonstrated in the case studies using Escherichia coli GRN and 100-gene gold-standard GRNs from DREAM 4 network inference challenge, by accounting for regulatory signs, TRaCE+ could extract more information from the KO data than TRaCE, leading to fewer uncertain edges. Importantly, iterating TRaCE+ with an optimal design of gene KOs could resolve the underdetermined issue of GRN inference in much fewer KO experiments than using TRaCE.ConclusionsTRaCE+ expands the applications of ensemble GRN inference strategy by accounting for the mode of the gene regulatory interactions. In comparison to TRaCE, TRaCE+ enables a better utilization of gene KO data, thereby reducing the cost of tackling underdetermined GRN inference. TRaCE+ subroutines for MATLAB are freely available at the following website: http://www.cabsel.ethz.ch/tools/trace.html.

A pancreatic exocrine-like cell regulatory circuit operating in the upper stomach of the sea urchin Strongylocentrotus purpuratus larva.

BackgroundDigestive cells are present in all metazoans and provide the energy necessary for the whole organism. Pancreatic exocrine cells are a unique vertebrate cell type involved in extracellular digestion of a wide range of nutrients. Although the organization and regulation of this cell type is intensively studied in vertebrates, its evolutionary history is still unknown. In order to understand which are the elements that define the pancreatic exocrine phenotype, we have analyzed the expression of genes that contribute to specification and function of this cell-type in an early branching deuterostome, the sea urchin Strongylocentrotus purpuratus.ResultsWe defined the spatial and temporal expression of sea urchin orthologs of pancreatic exocrine genes and described a unique population of cells clustered in the upper stomach of the sea urchin embryo where exocrine markers are co-expressed. We used a combination of perturbation analysis, drug and feeding experiments and found that in these cells of the sea urchin embryo gene expression and gene regulatory interactions resemble that of bona fide pancreatic exocrine cells. We show that the sea urchin Ptf1a, a key transcriptional activator of digestive enzymes in pancreatic exocrine cells, can substitute for its vertebrate ortholog in activating downstream genes.ConclusionsCollectively, our study is the first to show with molecular tools that defining features of a vertebrate cell-type, the pancreatic exocrine cell, are shared by a non-vertebrate deuterostome. Our results indicate that the functional cell-type unit of the vertebrate pancreas may evolutionarily predate the emergence of the pancreas as a discrete organ. From an evolutionary perspective, these results encourage to further explore the homologs of other vertebrate cell-types in traditional or newly emerging deuterostome systems.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-016-0686-0) contains supplementary material, which is available to authorized users.

Modular reorganization of the global network of gene regulatory interactions during perinatal human brain development.

BackgroundDuring early development of the nervous system, gene expression patterns are known to vary widely depending on the specific developmental trajectories of different structures. Observable changes in gene expression profiles throughout development are determined by an underlying network of precise regulatory interactions between individual genes. Elucidating the organizing principles that shape this gene regulatory network is one of the central goals of developmental biology. Whether the developmental programme is the result of a dynamic driven by a fixed architecture of regulatory interactions, or alternatively, the result of waves of regulatory reorganization is not known.ResultsHere we contrast these two alternative models by examining existing expression data derived from the developing human brain in prenatal and postnatal stages. We reveal a sharp change in gene expression profiles at birth across brain areas. This sharp division between foetal and postnatal profiles is not the result of pronounced changes in level of expression of existing gene networks. Instead we demonstrate that the perinatal transition is marked by the widespread regulatory rearrangement within and across existing gene clusters, leading to the emergence of new functional groups. This rearrangement is itself organized into discrete blocks of genes, each targeted by a distinct set of transcriptional regulators and associated to specific biological functions.ConclusionsOur results provide evidence of an acute modular reorganization of the regulatory architecture of the brain transcriptome occurring at birth, reflecting the reassembly of new functional associations required for the normal transition from prenatal to postnatal brain development.Electronic supplementary materialThe online version of this article (doi:10.1186/s12861-016-0111-3) contains supplementary material, which is available to authorized users.

Data-based Reconstruction of Gene Regulatory Networks of Fungal Pathogens.

In the emerging field of systems biology of fungal infection, one of the central roles belongs to the modeling of gene regulatory networks (GRNs). Utilizing omics-data, GRNs can be predicted by mathematical modeling. Here, we review current advances of data-based reconstruction of both small-scale and large-scale GRNs for human pathogenic fungi. The advantage of large-scale genome-wide modeling is the possibility to predict central (hub) genes and thereby indicate potential biomarkers and drug targets. In contrast, small-scale GRN models provide hypotheses on the mode of gene regulatory interactions, which have to be validated experimentally. Due to the lack of sufficient quantity and quality of both experimental data and prior knowledge about regulator–target gene relations, the genome-wide modeling still remains problematic for fungal pathogens. While a first genome-wide GRN model has already been published for Candida albicans, the feasibility of such modeling for Aspergillus fumigatus is evaluated in the present article. Based on this evaluation, opinions are drawn on future directions of GRN modeling of fungal pathogens. The crucial point of genome-wide GRN modeling is the experimental evidence, both used for inferring the networks (omics ‘first-hand’ data as well as literature data used as prior knowledge) and for validation and evaluation of the inferred network models.

How to Predict Molecular Interactions between Species?

Organisms constantly interact with other species through physical contact which leads to changes on the molecular level, for example the transcriptome. These changes can be monitored for all genes, with the help of high-throughput experiments such as RNA-seq or microarrays. The adaptation of the gene expression to environmental changes within cells is mediated through complex gene regulatory networks. Often, our knowledge of these networks is incomplete. Network inference predicts gene regulatory interactions based on transcriptome data. An emerging application of high-throughput transcriptome studies are dual transcriptomics experiments. Here, the transcriptome of two or more interacting species is measured simultaneously. Based on a dual RNA-seq data set of murine dendritic cells infected with the fungal pathogen Candida albicans, the software tool NetGenerator was applied to predict an inter-species gene regulatory network. To promote further investigations of molecular inter-species interactions, we recently discussed dual RNA-seq experiments for host-pathogen interactions and extended the applied tool NetGenerator (Schulze et al., 2015). The updated version of NetGenerator makes use of measurement variances in the algorithmic procedure and accepts gene expression time series data with missing values. Additionally, we tested multiple modeling scenarios regarding the stimuli functions of the gene regulatory network. Here, we summarize the work by Schulze et al. (2015) and put it into a broader context. We review various studies making use of the dual transcriptomics approach to investigate the molecular basis of interacting species. Besides the application to host-pathogen interactions, dual transcriptomics data are also utilized to study mutualistic and commensalistic interactions. Furthermore, we give a short introduction into additional approaches for the prediction of gene regulatory networks and discuss their application to dual transcriptomics data. We conclude that the application of network inference on dual-transcriptomics data is a promising approach to predict molecular inter-species interactions.

Gene Regulatory Network Inference of Immunoresponsive Gene 1 (IRG1) Identifies Interferon Regulatory Factor 1 (IRF1) as Its Transcriptional Regulator in Mammalian Macrophages.

Immunoresponsive gene 1 (IRG1) is one of the highest induced genes in macrophages under pro-inflammatory conditions. Its function has been recently described: it codes for immune-responsive gene 1 protein/cis-aconitic acid decarboxylase (IRG1/CAD), an enzyme catalysing the production of itaconic acid from cis-aconitic acid, a tricarboxylic acid (TCA) cycle intermediate. Itaconic acid possesses specific antimicrobial properties inhibiting isocitrate lyase, the first enzyme of the glyoxylate shunt, an anaplerotic pathway that bypasses the TCA cycle and enables bacteria to survive on limited carbon conditions. To elucidate the mechanisms underlying itaconic acid production through IRG1 induction in macrophages, we examined the transcriptional regulation of IRG1. To this end, we studied IRG1 expression in human immune cells under different inflammatory stimuli, such as TNFα and IFNγ, in addition to lipopolysaccharides. Under these conditions, as previously shown in mouse macrophages, IRG1/CAD accumulates in mitochondria. Furthermore, using literature information and transcription factor prediction models, we re-constructed raw gene regulatory networks (GRNs) for IRG1 in mouse and human macrophages. We further implemented a contextualization algorithm that relies on genome-wide gene expression data to infer putative cell type-specific gene regulatory interactions in mouse and human macrophages, which allowed us to predict potential transcriptional regulators of IRG1. Among the computationally identified regulators, siRNA-mediated gene silencing of interferon regulatory factor 1 (IRF1) in macrophages significantly decreased the expression of IRG1/CAD at the gene and protein level, which correlated with a reduced production of itaconic acid. Using a synergistic approach of both computational and experimental methods, we here shed more light on the transcriptional machinery of IRG1 expression and could pave the way to therapeutic approaches targeting itaconic acid levels.

GeneNetFinder2: Improved Inference of Dynamic Gene Regulatory Relations with Multiple Regulators.

A gene involved in complex regulatory interactions may have multiple regulators since gene expression in such interactions is often controlled by more than one gene. Another thing that makes gene regulatory interactions complicated is that regulatory interactions are not static, but change over time during the cell cycle. Most research so far has focused on identifying gene regulatory relations between individual genes in a particular stage of the cell cycle. In this study we developed a method for identifying dynamic gene regulations of several types from the time-series gene expression data. The method can find gene regulations with multiple regulators that work in combination or individually as well as those with single regulators. The method has been implemented as the second version of GeneNetFinder (hereafter called GeneNetFinder2) and tested on several gene expression datasets. Experimental results with gene expression data revealed the existence of genes that are not regulated by individual genes but rather by a combination of several genes. Such gene regulatory relations cannot be found by conventional methods. Our method finds such regulatory relations as well as those with multiple, independent regulators or single regulators, and represents gene regulatory relations as a dynamic network in which different gene regulatory relations are shown in different stages of the cell cycle. GeneNetFinder2 is available at http://bclab.inha.ac.kr/GeneNetFinder and will be useful for modeling dynamic gene regulations with multiple regulators.

Forward Individualized Medicine from Personal Genomes to Interactomes.

When considering the variation in the genome, transcriptome, proteome and metabolome, and their interaction with the environment, every individual can be rightfully considered as a unique biological entity. Individualized medicine promises to take this uniqueness into account to optimize disease treatment and thereby improve health benefits for every patient. The success of individualized medicine relies on a precise understanding of the genotype-phenotype relationship. Although omics technologies advance rapidly, there are several challenges that need to be overcome: Next generation sequencing can efficiently decipher genomic sequences, epigenetic changes, and transcriptomic variation in patients, but it does not automatically indicate how or whether the identified variation will cause pathological changes. This is likely due to the inability to account for (1) the consequences of gene-gene and gene-environment interactions, and (2) (post)transcriptional as well as (post)translational processes that eventually determine the concentration of key metabolites. The technologies to accurately measure changes in these latter layers are still under development, and such measurements in humans are also mainly restricted to blood and circulating cells. Despite these challenges, it is already possible to track dynamic changes in the human interactome in healthy and diseased states by using the integration of multi-omics data. In this review, we evaluate the potential value of current major bioinformatics and systems biology-based approaches, including genome wide association studies, epigenetics, gene regulatory and protein-protein interaction networks, and genome-scale metabolic modeling. Moreover, we address the question whether integrative analysis of personal multi-omics data will help understanding of personal genotype-phenotype relationships.

Transcription Factor Activity Mapping of a Tissue-Specific In Vivo Gene Regulatory Network

A wealth of physical interaction data between transcription factors (TFs) and DNA has been generated, but these interactions often do not have apparent regulatory consequences. Thus, equating physical interaction data with gene regulatory networks (GRNs) is problematic. Here, we comprehensively assay TF activity, rather than binding, to construct a network of gene regulatory interactions in the C. elegans intestine. By manually observing the in vivo tissue-specific knockdown of 921 TFs on a panel of 19 fluorescent transcriptional reporters, we identified a GRN of 411 interactions between 19 promoters and 177 TFs. This GRN shows only modest overlap with physical interactions, indicating that many regulatory interactions are indirect. We applied nested effects modeling to uncover information flow between TFs in the intestine that converges on a small set of physical TF-promoter interactions. We found numerous cell nonautonomous regulatory interactions, illustrating tissue-to-tissue communication. Altogether, our study illuminates the complexity of gene regulation in the context of a living animal.

Character trees from transcriptome data: Origin and individuation of morphological characters and the so-called "species signal".

We elaborate a framework for investigating the evolutionary history of morphological characters. We argue that morphological character trees generated by phylogenetic analysis of transcriptomes provide a useful tool for identifying causal gene expression differences underlying the development and evolution of morphological characters. They also enable rigorous testing of different models of morphological character evolution and origination, including the hypothesis that characters originate via divergence of repeated ancestral characters. Finally, morphological character trees provide evidence that character transcriptomes undergo concerted evolution. We argue that concerted evolution of transcriptomes can explain the so-called "species signal" found in several recent comparative transcriptome studies. The species signal is the phenomenon that transcriptomes cluster by species rather than character type, even though the characters are older than the respective species. We suggest the species signal is a natural consequence of concerted gene expression evolution resulting from mutations that alter gene regulatory network interactions shared by the characters under comparison. Thus, character trees generated from transcriptomes allow us to investigate the variational independence, or individuation, of morphological characters at the level of genetic programs.

Inferring interaction type in gene regulatory networks using co-expression data.

BackgroundKnowledge of interaction types in biological networks is important for understanding the functional organization of the cell. Currently information-based approaches are widely used for inferring gene regulatory interactions from genomics data, such as gene expression profiles; however, these approaches do not provide evidence about the regulation type (positive or negative sign) of the interaction.ResultsThis paper describes a novel algorithm, “Signing of Regulatory Networks” (SIREN), which can infer the regulatory type of interactions in a known gene regulatory network (GRN) given corresponding genome-wide gene expression data. To assess our new approach, we applied it to three different benchmark gene regulatory networks, including Escherichia coli, prostate cancer, and an in silico constructed network. Our new method has approximately 68, 70, and 100 percent accuracy, respectively, for these networks. To showcase the utility of SIREN algorithm, we used it to predict previously unknown regulation types for 454 interactions related to the prostate cancer GRN.ConclusionsSIREN is an efficient algorithm with low computational complexity; hence, it is applicable to large biological networks. It can serve as a complementary approach for a wide range of network reconstruction methods that do not provide information about the interaction type.Electronic supplementary materialThe online version of this article (doi:10.1186/s13015-015-0054-4) contains supplementary material, which is available to authorized users.

Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes as a Model for Heart Development and Congenital Heart Disease

Congenital heart disease (CHD) remains a significant health problem, with a growing population of survivors with chronic disease. Despite intense efforts to understand the genetic basis of CHD in humans, the etiology of most CHD is unknown. Furthermore, new models of CHD are required to better understand the development of CHD and to explore novel therapies for this patient population. In this review, we highlight the role that human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes can serve to enhance our understanding of the development, pathophysiology and potential therapeutic targets for CHD. We highlight the use of hiPSC-derived cardiomyocytes to model gene regulatory interactions, cell-cell interactions and tissue interactions contributing to CHD. We further emphasize the importance of using hiPSC-derived cardiomyocytes as personalized research models. The use of hiPSCs presents an unprecedented opportunity to generate disease-specific cellular models, investigate the underlying molecular mechanisms of disease and uncover new therapeutic targets for CHD.

Expression noise facilitates the evolution of gene regulation.

Although it is often tacitly assumed that gene regulatory interactions are finely tuned, how accurate gene regulation could evolve from a state without regulation is unclear. Moreover, gene expression noise would seem to impede the evolution of accurate gene regulation, and previous investigations have provided circumstantial evidence that natural selection has acted to lower noise levels. By evolving synthetic Escherichia coli promoters de novo, we here show that, contrary to expectations, promoters exhibit low noise by default. Instead, selection must have acted to increase the noise levels of highly regulated E. coli promoters. We present a general theory of the interplay between gene expression noise and gene regulation that explains these observations. The theory shows that propagation of expression noise from regulators to their targets is not an unwanted side-effect of regulation, but rather acts as a rudimentary form of regulation that facilitates the evolution of more accurate regulation.

Yeast protein-protein interaction network model based on biological experimental data

Duplication and divergence have been widely recognized as the two dominant evolutionary forces in shaping biological networks, e.g., gene regulatory networks and protein-protein interaction (PPI) networks. It has been shown that the network growth models constructed on the principle of duplication and divergence can recapture the topological properties of real PPI networks. However, such network models only consider the evolution processes. How to select the model parameters with the real biological experimental data has not been presented. Therefore, based on the real PPI network statistical data, a yeast PPI network model is constructed. The simulation results indicate that the topological characteristics of the constructed network model are well consistent with those of real PPI networks, especially on sparseness, scale-free, small-world, hierarchical modularity, and disassortativity.

Analysis of gene expression patterns and levels in maize hybrids and their parents.

Heterosis has greatly contributed to conventional plant breeding and is widely used to increase crop plant productivity. However, although some studies have explored the mechanisms of heterosis at the genomic and transcriptome level, these mechanisms still remain unclear. The growth and development of maize seedlings and immature embryos have an important impact on subsequent production. This study investigated differentially expressed genes (DEGs) between parents and reciprocal hybrids in the seedling leaves, roots, and immature embryo 15 days after pollination using amplified fragment length polymorphism (AFLP)-based transcript profiling (cDNA-AFLP). We isolated 180, 170, and 108 genes from the leaves, roots, and immature embryos, respectively, that were differentially expressed between hybrids and parents. Sequencing and functional analysis revealed that 107 transcript-derived fragments in the roots and leaves and 90 in the immature embryos were involved in known functions, whereas many DEGs had roles in plant growth and development, photosynthesis, signal transduction, and seed germination. Quantitative reverse-transcription polymerase chain reaction analysis of relative expression levels between reciprocal hybrids and both parental genotypes of selected genes produced results that were consistent with cDNA-AFLP. We validated the expression patterns of 15 selected genes related to heterosis formation and revealed that most showed non-additive expression in one or both hybrids, including dominant, underdominant, and overdominant expression. This indicates that gene-regulatory interactions among parental alleles play an important role in heterosis during the early developmental stages of maize.

K-core decomposition of a protein domain co-occurrence network reveals lower cancer mutation rates for interior cores.

BackgroundNetwork biology currently focuses primarily on metabolic pathways, gene regulatory, and protein-protein interaction networks. While these approaches have yielded critical information, alternative methods to network analysis will offer new perspectives on biological information. A little explored area is the interactions between domains that can be captured using domain co-occurrence networks (DCN). A DCN can be used to study the function and interaction of proteins by representing protein domains and their co-existence in genes and by mapping cancer mutations to the individual protein domains to identify signals.ResultsThe domain co-occurrence network was constructed for the human proteome based on PFAM domains in proteins. Highly connected domains in the central cores were identified using the k-core decomposition technique. Here we show that these domains were found to be more evolutionarily conserved than the peripheral domains. The somatic mutations for ovarian, breast and prostate cancer diseases were obtained from the TCGA database. We mapped the somatic mutations to the individual protein domains and the local false discovery rate was used to identify significantly mutated domains in each cancer type. Significantly mutated domains were found to be enriched in cancer disease pathways. However, we found that the inner cores of the DCN did not contain any of the significantly mutated domains. We observed that the inner core protein domains are highly conserved and these domains co-exist in large numbers with other protein domains.ConclusionMutations and domain co-occurrence networks provide a framework for understanding hierarchal designs in protein function from a network perspective. This study provides evidence that a majority of protein domains in the inner core of the DCN have a lower mutation frequency and that protein domains present in the peripheral regions of the k-core contribute more heavily to the disease. These findings may contribute further to drug development.Electronic supplementary materialThe online version of this article (doi:10.1186/s13336-015-0016-6) contains supplementary material, which is available to authorized users.

Network representations of immune system complexity.

The mammalian immune system is a dynamic multiscale system composed of a hierarchically organized set of molecular, cellular, and organismal networks that act in concert to promote effective host defense. These networks range from those involving gene regulatory and protein–protein interactions underlying intracellular signaling pathways and single‐cell responses to increasingly complex networks of in vivo cellular interaction, positioning, and migration that determine the overall immune response of an organism. Immunity is thus not the product of simple signaling events but rather nonlinear behaviors arising from dynamic, feedback‐regulated interactions among many components. One of the major goals of systems immunology is to quantitatively measure these complex multiscale spatial and temporal interactions, permitting development of computational models that can be used to predict responses to perturbation. Recent technological advances permit collection of comprehensive datasets at multiple molecular and cellular levels, while advances in network biology support representation of the relationships of components at each level as physical or functional interaction networks. The latter facilitate effective visualization of patterns and recognition of emergent properties arising from the many interactions of genes, molecules, and cells of the immune system. We illustrate the power of integrating ‘omics’ and network modeling approaches for unbiased reconstruction of signaling and transcriptional networks with a focus on applications involving the innate immune system. We further discuss future possibilities for reconstruction of increasingly complex cellular‐ and organism‐level networks and development of sophisticated computational tools for prediction of emergent immune behavior arising from the concerted action of these networks. WIREs Syst Biol Med 2015, 7:13–38. doi: 10.1002/wsbm.1288This article is categorized under: 1Analytical and Computational Methods > Computational Methods2Laboratory Methods and Technologies > Macromolecular Interactions, Methods

Gene regulatory systems that control gene expression in the Ciona embryo.

Transcriptional control of gene expression is one of the most important regulatory systems in animal development. Specific gene expression is basically determined by combinatorial regulation mediated by multiple sequence-specific transcription factors. The decoding of animal genomes has provided an opportunity for us to systematically examine gene regulatory networks consisting of successive layers of control of gene expression. It remains to be determined to what extent combinatorial regulation encoded in gene regulatory networks can explain spatial and temporal gene-expression patterns. The ascidian Ciona intestinalis is one of the animals in which the gene regulatory network has been most extensively studied. In this species, most specific gene expression patterns in the embryo can be explained by combinations of upstream regulatory genes encoding transcription factors and signaling molecules. Systematic scrutiny of gene expression patterns and regulatory interactions at the cellular resolution have revealed incomplete parts of the network elucidated so far, and have identified novel regulatory genes and novel regulatory mechanisms.