Spatiotemporal network motif reveals the biological traits of developmental gene regulatory networks in Drosophila melanogaster
© Kim et al.; licensee BioMed Central Ltd. 2012
Received: 16 December 2011
Accepted: 1 May 2012
Published: 1 May 2012
Network motifs provided a “conceptual tool” for understanding the functional principles of biological networks, but such motifs have primarily been used to consider static network structures. Static networks, however, cannot be used to reveal time- and region-specific traits of biological systems. To overcome this limitation, we proposed the concept of a “spatiotemporal network motif,” a spatiotemporal sequence of network motifs of sub-networks which are active only at specific time points and body parts.
On the basis of this concept, we analyzed the developmental gene regulatory network of the Drosophila melanogaster embryo. We identified spatiotemporal network motifs and investigated their distribution pattern in time and space. As a result, we found how key developmental processes are temporally and spatially regulated by the gene network. In particular, we found that nested feedback loops appeared frequently throughout the entire developmental process. From mathematical simulations, we found that mutual inhibition in the nested feedback loops contributes to the formation of spatial expression patterns.
Taken together, the proposed concept and the simulations can be used to unravel the design principle of developmental gene regulatory networks.
To uncover the governing principles underlying complex biological processes, it is important to understand the relationship between topological structures and the dynamical characteristics of gene regulatory networks [1–4]. One promising method of investigation is to disassemble the large regulatory network into its more basic, constituent building blocks called network motifs, which recur within a network much more often than expected in random networks. Network motifs are considered to have been evolutionarily selected because of their functional advantages .
Most previous studies have identified network motifs of biological networks by implicitly assuming that all the links in a network can be active or working at the same time. However, such approaches may not be applicable to developmental networks where all genes and interactions do not operate simultaneously due to spatial and temporal variations. Some network motif approaches have partially considered spatial or temporal information on biological networks [6, 7]. Papatsenko analyzed the dynamics of network motifs for a spatial stripe pattern formation, only in early embryogenesis , while Kim et al. explored the dynamics for temporal network motifs . Nevertheless, patterns of spatiotemporal variations in gene regulatory networks have not yet been explored.
In this paper, we propose a novel concept called the “spatiotemporal network motif,” which is a sequence of network motifs in sub-networks that are spatiotemporally active. These network motifs are constructed by re-organizing the regulations between spatiotemporally expressed genes.
We applied this approach to the developmental gene regulatory network of D. melanogaster. First, we identified a spatio-temporal sequence of network motifs which change according to developmental stages and regions. Then, we analyzed the pattern of spatio-temporal network motifs and their dynamics (we only considered three-node network motifs for simplicity and for a concise illustration of the method). As a result, we found that the most frequently observed structure in the spatiotemporal network motif pattern is the feed-forward loop structure. This result implies that signal-processing via feed-forward loops is required throughout all of the development stages [6, 8]. Another important network motif that we identified was nested feedback loops, where one feedback loop is nested inside another feedback or feed-forward loop. Such nested feedback loops were considered necessary for the development of a central nervous system, as they should be stable and robust against both noise and small perturbations . This result suggests that nested feedback loops might play an important role in the elaborate regulation of developmental processes. Interestingly, we found that most nested feedback loops had mutual inhibitory structures. Through mathematical simulations, we showed that such mutual inhibitory structures can enable exclusive spatial expression of gap genes.
Taken together, the proposed concept and the simulations can reveal time- and region-specific biological traits in dynamic processes such as the developmental gene regulatory network, and can be widely used to investigate the relationship between dynamic network structures and their regulatory functions.
Identification of spatio-temporal network motifs
Let us consider an example network with gene expression data measurements taken in three regions (R1, R2 and R3) at three time points (T1, T2 and T3), as shown in Figure 1b. In this case, we should consider nine active sub-networks. Figure 1b shows the network motif pattern of the nine active sub-networks. The network motifs of the R3 sub-network do not vary as time evolves, whereas the patterns of network motifs in the R1 or R2 sub-networks vary. In addition, at each fixed point in time, the network motifs of the active sub-networks were different in all three regions. For instance, Motif 3 was found in R1, but not in other regions. This suggested that Motif 3 is necessary for some particular biological function only found in R1. Likewise, we can use the proposed approach to identify network structures that are necessary at some specific region and time.
Unravelling region- and time-specific biological traits of developmental processes
We have summarized body part(s), network motifs, triple genes, and GO terms for three periods (early, middle, and late stages)
GO term (p-value)
bcd, hb, kni
dors ectoderm AISN
Kr, hb, gt
vent ectoderm AISN
procephalic ectoderm AISN
dors ectoderm AISN
eve, en, zen
Post endoderm PR
tll, Kr, salm
tll, hb, salm
vent nerve cord PR P3
cas, Kr, nub
vent nerve cord PR P3
Kr, hb, nub
vent ectoderm PR
hb, en, eve
odd, prd, en
vent epidermis PR
vent epidermis PR
nkd, en, wg
slp1, en, wg
vent epidermis PR
dors epidermis PR
Kr, hb, nub
vent nerve cord
gt, Kr, hkb
hb, Kr, Nos
In the middle period of embryogenesis (Stages 7–10), one of the key events is the completion of gastrulation. During gastrulation, the morphology of an embryo is rearranged to form the three germ layers: ectoderm, mesoderm, and endoderm . Another triplet of genes—tll Kr, and salm—were expressed in the endoderm tissue. The gene tll is one of the gastrulation genes and salm is well-known to influence the development of gut. Another triplet of genes—cas hb Kr, and nub—were related to determining the fate of the ganglion mother cell (GO: 0007402), and were expressed in the ventral nerve cord tissue. This suggested that the identified genes might influence the development of the central nervous system (CNS). The midline precursors undergo a synchronous cell division to give rise to 16 midline progenitor cells per segment. Then, until stage 12, these 16 midline cells go through cell shape changes, cell division, and differentiation to form the midline primordium (PR) .
Illustration of the GO analysis for late embryogenesis
GO:0048513 organ development
GO:0007423 sensory organ development
GO:0007507 heart development
GO:0061061 muscle structure development
GO:0007498 mesoderm development
GO:0014016 neuroblast differentiation
GO:0007400 neuroblast fate determination
GO:0007419 ventral cord development
GO:0008544 epidermis development
Dynamical and structural analysis of the gap gene network
We found that nested feedback loops (Motif 6, 10, 11, 12 and 13) were frequently observed in the developmental network, especially in the gap gene network. The combination of feed-forward loops and feedback loops can induce an emergent property . For example, in stages 4–16, Motifs 4–8 and 10–13 are commonly observed in the ectoderm. These stages include the sequential development of the central nervous system and epidermis in the ectoderm tissue. These processes are both important in ectoderm-specific development, and are considered to require stable and robust regulation against both noises and small perturbations . The network of Drosophila neurogenesis also features nested feedback loops . Nakajima et al. suggested that nested feedback loops composed of hb Kr nub, and cas are precisely regulated by three different kinds of links. The authors also found a minimum network which can reproduce the sequential expression pattern of the four genes. However, the minimum network is less robust against parameter variations than the original network with nested feedback loops. Hence, we conclude that the nested feedback loops induce an emergent property for the elaborate and robust regulations of developmental processes.
Why are the gap genes connected with each other through such a complex structure? In order to understand the dynamics of the complex network structure, we performed Boolean simulations based on the gap gene network (Figure 3b). Then, we identified attractors, where an attractor means a set towards which a dynamical system evolves over time (Figure 3c). Each attractor can be mapped into the lineage-associated transcription factor activities [21–23] and can represent a developmental cell fate. For example, in the model of myeloid progenitor cell differentiation, there are two attractors where GATA-1 and PU.1 are exclusively expressed, which correspond to erythroid/megakaryocyte and myeloid-monocytic fates, respectively . Here, an attractor is expressed as a combination of the digits 1 and 0, where 1 denotes an expressed state of a gene and 0 denotes an unexpressed state of a gene in a specific region. As a result of the Boolean simulation for five gap genes (hb, gt, Kr, kni, and tll), we identified eight attractors within the gap gene network (Additional file 1: Figure S3 and Table S1). The gap genes are activated or repressed by the maternal effect genes (bcd and cad) (Figure 3B). The concentration of the protein bcd is high in the anterior region , while the concentration of the protein cad is high in the posterior region. The spatial-specific concentration of maternal effect genes affects the spatial-specific concentration of gap genes, which in turn results in distinct steady states which correspond to attractors. The attractors can be classified into two classes: attractors corresponding to the case where only one out of five gap genes is expressed (attractors 01000, 00100, 00010, 00001 and 10000) and attractors corresponding to the case where more than one gene are expressed (attractors 10001 and 00101). For example, an attractor 00100 means that only Kr is expressed and the rest of genes are unexpressed. It is well known that among gap genes, only Kr is expressed in the central region [25–27] . Hence, the attractor 00100 corresponds to a state of gap gene expression at a special region (i.e., the body part of 50% egg length (EL). Likewise, another attractor 10001 corresponds to the spatial expression pattern in the posterior pole region (Figure 3d) where only two genes, hb and tll, are expressed [28, 29], and the rest of genes are unexpressed. Similarily, two other attractors 00010 and 01000 also correspond to a set of spatial expression patterns of gap genes. Hence, we conclude that the complex structure of the gap gene network with nested feedback loops contributes to the formation of spatial expression patterns.
The development of multi-cellular organisms relies on the coordinated spatiotemporal regulation of gene expressions. To unravel the organizing principles of developmental gene regulatory networks, it is crucial to understand the relationship between the structure and function of spatiotemporal subnetworks. However, thus far, such spatiotemporal networks have not been investigated. Hence, we proposed a new concept called a “spatiotemporal network motif,” which is a sequence of network motifs in time and space, and we applied this concept to analyze the developmental gene regulatory network of D. melanogaster. We found that the results of our approach coincide with spatially specific processes in early, middle, and late embryogenesis.
We also identified patterns of spatiotemporal network motifs and analyzed the relationship between network structures and their biological functions (Figure 2). We found that the most frequently observed structure in the spatiotemporal network motif pattern was the feed-forward loop structure (Motif 4) (79% for 19 sub-networks). This finding is also well supported by recent studies: Motif 4 is found ubiquitously in the D. melanogaster gene regulatory network ; Motif 4 is the core structure of the D. melanogaster gene regulatory network ; Motif 4 plays an essential role in the D. melanogaster central nervous system . This structure, with various regulation types, has several important functions in a biological network, such as detecting persistent signals, generating pulse, and accelerating response . For example, in the feed-forward loop consisting of the three genes eve, en, and hb, the genes are related to specifying cell fate (GO: 0001708) and commitment (GO: 0045165). These triple genes form a coherent type-1 of Motif 4, such that eve and en are activated by hb, and en is activated by eve in ectoderm tissue at stages 4–6. In this developmental process, the dynamics of the coherent type-1  of Motif 4 can be used as a persistent signal detector, which enables the system to respond only to persistent signals while neglecting short-term signals. This means that developmental processes related to cell fate must be robust in relation to noises, and explains how the developmental network deals with noises via the structure of the coherent type-1 of Motif 4.
Interestingly, we found that nested feedback loops were frequently observed in the gap gene network and most of the nested feedback loops contain mutual inhibition structures. Based on Boolean simulations, we showed that the gap gene network has a significantly large number of attractors (eight attractors) and such many attractors in the network are attributable to mutual inhibition. Hence we infer that the gap gene network might have evolved to induce a large number of attractors (by increasing the number of mutual inhibitions) which correspond to various developmental states.
The interlinked incoherent feed-forward loop structure is a key regulatory structure for stripe formation at 4–6 stages in the maternal region  and we identified this network motif (ID 6) at the same spatiotemporal developmental stages. The triple genes (gt, Kr and eve) of the network motif ID 6 are also consistent with the previous study . In addition, it is well-known that the feed-forward loop is a crucial structure to DV (Droso-Ventral) axis formation at 4–6 stages in the maternal region  and we also identified this network motif (ID 4) at the same spatiotemporal stages. From these, we conclude we could infer the design principles of Drosophila development in a holistic manner using our approach.
Network motifs cannot uniquely determine the whole dynamical properties of a regulatory network. In general, the dynamics of a regulatory network depends on multiple factors such as initial conditions, cellular environments, and randomness [34, 35]. However, some particular dynamical properties can be determined by certain network structures . For instance, bistable switching cannot be realized without a positive feedback loop in the regulatory network. So, understanding the relationship between network structures and dynamics may still be useful as we can infer some possible dynamical characteristics of a network from its structure. The proposed approach guides us to find specific network motifs (e.g., positive feedback) at a specific spatiotemporal stage and therefore we can estimate possible dynamical properties (e.g., bistable switch) related to the identified network motifs. Taken together, our approach is useful to infer developmental functions of spatiotemporally varying cells based on identification of network motifs.
Most of the previous studies identified network motifs of the whole regulatory networks integrated from various literatures without considering particular biological contexts (e.g., environmental conditions, developmental stages, etc.). The key difference of our study from the previous ones is the identification of network motifs depending on active sub-network assuming that only part of genes may express under some particular spatiotemporal condition. Such a concept has not been proposed so far. This concept provides us (time- and space-) varying patterns of network motifs in terms of time and region simultaneously. For instance, we can find out many types of network motifs at the 4–8 stages, while there is no network motif at the 9–16 stages in the maternal region (Figure 2). From this, we can infer that the maternal region requires more complex regulation through several types of network motifs at the former stages compared to the latter stages.
The topological structures of the network motifs that we discovered in this study are not new by themselves. However, the sequence of time- and space-varying network motifs is new. Furthermore, we could associate the dynamical properties of identified network motifs and spatiotemporal developmental processes of Drosophila. For instance, it is well known that the major developmental process at 4–6 stages is differentiation, but there is no differentiation at 1–3 stages in the maternal region. Interestingly, we found the mutual inhibition network motif (ID 8), a key network motif for the differentiation process, at 4–6 stages but not at 1–3 stages. Together, we can infer a specific developmental process at a specific developmental stage from the dynamical properties of the identified network motifs. Furthermore, the presented approach provides us with a useful and single framework in which we can investigate the whole developmental process in a comprehensive view.
We proposed a novel concept called the “spatiotemporal network motif,” which is a sequence of network motifs in sub-networks that are spatiotemporally active. Since the proposed approach is based on the network motif framework, many important issues must still be addressed, including the reliability of the constructed network and the justification of DEG selection. Nonetheless, the proposed approach can provide a good framework for improving our understanding of developmental processes and identifying key regulatory processes. Also, by applying the proposed approach to a developmental network, we can gain new insights into the organizing principles of a developmental network whose structures change spatially and temporally, and can be widely used to investigate the relationship between dynamic network structures and their regulatory functions.
Identification of active sub-networks
In order to reconstruct active sub-networks from a whole gene regulatory network, we need two kinds of information: an integrated gene regulatory network and information on the differentially expressed genes (DEGs) in each region and stage. To acquire the DEGs, we incorporated spatiotemporal information from the BDGP database (see the following subsection—“The BDGP database for spatio-temporal information” for details) . The gene regulatory network data of D. melanogaster were retrieved from the TRANSFAC database , which is a manually curated transcriptional network database. The network consists of 155 nodes and 377 links (see Additional file 2). To identify which sub-network is active at each time point and region, we used the following information: if two genes that share an edge in the gene regulatory network are also included in the DEG set, then those two genes and the connecting link constitute the active sub-network. By repeating this procedure, we obtain a set of spatio-temporally active sub-networks (Additional file 1: Figure S2).
Microarray data and BDGP database for spatiotemporal information
In this dataset, the expression levels were measured at 22 sequential time points for 16 h (stages 1–16) during embryogenesis . All the data were measured with the time interval of 1 h and the sampling time points during the first 6 h were selected with 1-h overlapping period (from 1.5-2.5 h to 6.5-7.5 h). All the data were verified through the morphological trait of each developmental stage . Therefore, the time points we used were not arbitrarily chosen, but carefully selected to represent molecular profiles of each developmental stage. From the given time-series microarray data, we transformed the time window into the corresponding developmental stage using FlyMove .
The BDGP database contains 97,732 digital photographs and expression data for 7,152 genes. The images representing gene expression patterns were classified according to pre-defined stage range annotations. For temporal information, the BDGP data for the first 16 stages of embryogenesis were divided into six stage ranges (stages 1–3, 4–6, 7–8, 9–10, 11–12, and 13–16). For annotation information (http://bigbrain.lbl.gov/cgi-bin/benb/vocab_selector.pl?noreturn=1http://www.fruitfly-.org/ex/Annotation.htm), we obtained a controlled vocabulary to annotate gene expression patterns during embryogenesis. For spatial information, we used the annotations to group developmental structures into six organ systems : Maternal (Ma), Endoderm (En), Mesoderm (Me), Ectoderm (Ec), Central Nervous System (CNS), and Epidermis (Epi). Each body part was also subdivided into four types of anatomical structures (with a suffix): anlage in statu nascendi (AISN); anlage; primordium (PR) (primordium usually develops from an anlage, and can give rise to one or more differentiated organs); differentiated organs. Additional file 1: Figure S2 summarizes the procedure for selecting the DEG sets, and the results are provided in Additional file 2. The procedure for reconstructing spatiotemporally active sub-networks (for details, see ) is illustrated in Additional file 1: Figure S3.
Motif enrichment analysis
In order to identify and visualize network motifs using MAVisto, we generated 1,000 randomized networks by randomly reshuffling links while preserving the in- and out-degree of each node in the original network . Next, we identified two-node and three-node motifs for the gene regulatory network of D. melanogaster (using 0.01 as the cut-off threshold for network motifs).
GO (gene ontology) analysis
GO annotation is used to indicate gene traits. GO analysis is a useful tool for both small- and large-scale analysis. In this study, GO functional annotations were obtained from the GO database . To evaluate the statistical significance of the overlap between selected genes, we used AmiGo  to perform GO terms enrichment . FlyBase  was used as a database filter in the GO analysis.
Boolean network modeling and simulation analysis
We are grateful to A. Fletcher, O. Cominetti, and J.-K. Won for their valuable comments on the manuscript. This work was supported by the National Research Foundation of Korea (NRF) grants funded by the Korea Government, the Ministry of Education, Science & Technology (MEST) (2009–0086964, 2010–0017662, and 2011–0006314). It was also supported by the WCU (World Class University) program (R32-2008-000-10218-0) through the NRF funded by MEST and by GIST Systems Biology Infrastructure Establishment Grant (2012).
- Adryan B, Teichmann SA: The developmental expression dynamics of Drosophila melanogaster transcription factors. Genome Biol. 2010, 11: R40-10.1186/gb-2010-11-4-r40.View ArticleGoogle Scholar
- Damle S, Davidson EH: Precise cis-regulatory control of spatial and temporal expression of the alx-1 gene in the skeletogenic lineage of s. purpuratus. Dev Biol. 2011, 357: 505-517. 10.1016/j.ydbio.2011.06.016.View ArticleGoogle Scholar
- Hinman VF, Nguyen AT, Cameron RA, Davidson EH: Developmental gene regulatory network architecture across 500 million years of echinoderm evolution. Proc Natl Acad Sci U S A. 2003, 100: 13356-13361. 10.1073/pnas.2235868100.View ArticleGoogle Scholar
- Smith J, Theodoris C, Davidson EH: A gene regulatory network subcircuit drives a dynamic pattern of gene expression. Science. 2007, 318: 794-797. 10.1126/science.1146524.View ArticleGoogle Scholar
- Alon U: An Introduction to Systems Biology. 2006, CRC Press, London, UKGoogle Scholar
- Kim MS, Kim JR, Cho KH: Dynamic network rewiring determines temporal regulatory functions in Drosophila melanogaster development processes. BioEssays. 2010, 32: 505-513. 10.1002/bies.200900169.View ArticleGoogle Scholar
- Papatsenko D: Stripe formation in the early fly embryo: principles, models, and networks. BioEssays. 2009, 31: 1172-1180. 10.1002/bies.200900096.View ArticleGoogle Scholar
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298: 824-827. 10.1126/science.298.5594.824.View ArticleGoogle Scholar
- Prill RJ, Iglesias PA, Levchenko A: Dynamic properties of network motifs contribute to biological network organization. PLoS Biol. 2005, 3: e343-10.1371/journal.pbio.0030343.View ArticleGoogle Scholar
- Tomancak P, Beaton A, Weiszmann R, Kwan E, Shu S, Lewis SE, Richards S, Ashburner M, Hartenstein V, Celniker SE, Rubin GM: Systematic determination of patterns of gene expression during Drosophila embryogenesis. Genome Biol. 2002, 3: RESEARCH0088-View ArticleGoogle Scholar
- Leptin M, Grunewald B: Cell shape changes during gastrulation in Drosophila. Development. 1990, 110: 73-84.Google Scholar
- Kearney JB, Wheeler SR, Estes P, Parente B, Crews ST: Gene expression profiling of the developing Drosophila CNS midline cells. Dev Biol. 2004, 275: 473-492. 10.1016/j.ydbio.2004.08.047.View ArticleGoogle Scholar
- Nakajima A, Isshiki T, Kaneko K, Ishihara S: Robustness under functional constraint: the genetic network for temporal expression in Drosophila neurogenesis. PLoS Comput Biol. 2010, 6: e1000760-10.1371/journal.pcbi.1000760.View ArticleGoogle Scholar
- Brandman O, Ferrell JE: Li R, Meyer T: Interlinked fast and slow positive feedback loops drive reliable cell decisions. Science. 2005, 310: 496-498. 10.1126/science.1113834.View ArticleGoogle Scholar
- Sanchez L, Thieffry D: A logical analysis of the Drosophila gap-gene system. J Theor Biol. 2001, 211: 115-141. 10.1006/jtbi.2001.2335.View ArticleGoogle Scholar
- Sanchez L, Thieffry D: Segmenting the fly embryo: a logical analysis of the pair-rule cross-regulatory module. J Theor Biol. 2003, 224: 517-537. 10.1016/S0022-5193(03)00201-7.View ArticleGoogle Scholar
- Seo CH, Kim JR, Kim MS, Cho KH: Hub genes with positive feedbacks function as master switches in developmental gene regulatory networks. Bioinformatics. 2009, 25: 1898-1904. 10.1093/bioinformatics/btp316.View ArticleGoogle Scholar
- Davidson EH, Levine MS: Properties of developmental gene regulatory networks. Proc Natl Acad Sci U S A. 2008, 105: 20063-20066. 10.1073/pnas.0806007105.View ArticleGoogle Scholar
- Gilbert SF: Developmental Biology. 2006, Sinauer Associates, Sunderland MAGoogle Scholar
- Jaeger J, Blagov M, Kosman D, Kozlov KN, Manu , Myasnikova E, Surkova S, Vanario-Alonso CE, Samsonova M, Sharp DH, Reinitz J: Dynamical analysis of regulatory interactions in the gap gene system of Drosophila melanogaster. Genetics. 2004, 167: 1721-1737. 10.1534/genetics.104.027334.View ArticleGoogle Scholar
- Enver T, Pera M, Peterson C, Andrews PW: Stem cell states, fates, and the rules of attraction. Cell Stem Cell. 2009, 4: 387-397. 10.1016/j.stem.2009.04.011.View ArticleGoogle Scholar
- Huang S: Reprogramming cell fates: reconciling rarity with robustness. BioEssays. 2009, 31: 546-560. 10.1002/bies.200800189.View ArticleGoogle Scholar
- Rolls ET, Webb TJ: Cortical attractor network dynamics with diluted connectivity. Brain Res. 2012, 1434: 212-225.View ArticleGoogle Scholar
- Manu , Surkova S, Spirov AV, Gursky VV, Janssens H, Kim AR, Radulescu O, Vanario-Alonso CE, Sharp DH, Samsonova M, Reinitz J: Canalization of gene expression and domain shifts in the Drosophila blastoderm by dynamical attractors. PLoS Comput Biol. 2009, 5: e1000303-10.1371/journal.pcbi.1000303.View ArticleGoogle Scholar
- Gaul U, Weigel D: Regulation of Kruppel expression in the anlage of the Malpighian tubules in the Drosophila embryo. Mech Dev. 1990, 33: 57-67. 10.1016/0925-4773(90)90135-9.View ArticleGoogle Scholar
- Hoch M, Schroder C, Seifert E, Jackle H: cis-acting control elements for Kruppel expression in the Drosophila embryo. EMBO J. 1990, 9: 2587-2595.Google Scholar
- Manu , Surkova S, Spirov AV, Gursky VV, Janssens H, Kim AR, Radulescu O, Vanario-Alonso CE, Sharp DH, Samsonova M, Reinitz J: Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation. PLoS Biol. 2009, 7: e1000049-View ArticleGoogle Scholar
- Bronner G, Jackle H: Control and function of terminal gap gene activity in the posterior pole region of the Drosophila embryo. Mech Dev. 1991, 35: 205-211. 10.1016/0925-4773(91)90019-3.View ArticleGoogle Scholar
- Margolis JS, Borowsky ML, Steingrimsson E, Shim CW, Lengyel JA, Posakony JW: Posterior stripe expression of hunchback is driven from two promoters by a common enhancer element. Development. 1995, 121: 3067-3077.Google Scholar
- Cotterell J, Sharpe J: An atlas of gene regulatory networks reveals multiple three-gene mechanisms for interpreting morphogen gradients. Mol Syst Biol. 2010, 6: 425-View ArticleGoogle Scholar
- Karlsson D, Baumgardt M, Thor S: Segment-specific neuronal subtype specification by the integration of anteroposterior and temporal cues. PLoS Biol. 2010, 8: e1000368-10.1371/journal.pbio.1000368.View ArticleGoogle Scholar
- Mangan S, Alon U: Structure and function of the feed-forward loop network motif. Proc Natl Acad Sci U S A. 2003, 100: 11980-11985. 10.1073/pnas.2133841100.View ArticleGoogle Scholar
- Reeves GT, Trisnadi N, Truong TV, Nahmad M: Katz S. 2012, Dorsal-Ventral Gene Expression in the Drosophila Embryo Reflects the Dynamics and Precision of the Dorsal Nuclear Gradient. Dev Cell, Stathopoulos AGoogle Scholar
- Ingram PJ, Stumpf MP, Stark J: Network motifs: structure does not determine function. BMC Genomics. 2006, 7: 108-10.1186/1471-2164-7-108.View ArticleGoogle Scholar
- Zhu XM, Yin L, Hood L, Ao P: Calculating biological behaviors of epigenetic states in the phage lambda life cycle. Funct Integr Genomics. 2004, 4: 188-195.View ArticleGoogle Scholar
- Tyson JJ, Chen KC, Novak B: Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr Opin Cell Biol. 2003, 15: 221-231. 10.1016/S0955-0674(03)00017-6.View ArticleGoogle Scholar
- Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, et al: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, 34: D108-110. 10.1093/nar/gkj143.View ArticleGoogle Scholar
- Arbeitman MN, Furlong EE, Imam F, Johnson E, Null BH, Baker BS, Krasnow MA, Scott MP, Davis RW, White KP: Gene expression during the life cycle of Drosophila melanogaster. Science. 2002, 297: 2270-2275. 10.1126/science.1072152.View ArticleGoogle Scholar
- Weigmann K, Klapper R, Strasser T, Rickert C, Technau G, Jackle H, Janning W, Klambt C: FlyMove–a new way to look at development of Drosophila. Trends Genet. 2003, 19: 310-311. 10.1016/S0168-9525(03)00050-7.View ArticleGoogle Scholar
- Tomancak P, Berman BP, Beaton A, Weiszmann R, Kwan E, Hartenstein V, Celniker SE, Rubin GM: Global analysis of patterns of gene expression during Drosophila embryogenesis. Genome Biol. 2007, 8: R145-10.1186/gb-2007-8-7-r145.View ArticleGoogle Scholar
- Schreiber F, Schwobbermeyer H: MAVisto: a tool for the exploration of network motifs. Bioinformatics. 2005, 21: 3572-3574. 10.1093/bioinformatics/bti556.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29.Google Scholar
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25: 288-289. 10.1093/bioinformatics/btn615.View ArticleGoogle Scholar
- Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R, Zhang H: FlyBase: enhancing Drosophila Gene Ontology annotations. Nucleic Acids Res. 2009, 37: D555-559. 10.1093/nar/gkn788.View ArticleGoogle Scholar
- Li F, Long T, Lu Y, Ouyang Q, Tang C: The yeast cell-cycle network is robustly designed. Proc Natl Acad Sci U S A. 2004, 101: 4781-4786. 10.1073/pnas.0305937101.View ArticleGoogle Scholar