Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Biochim Biophys Acta. Author manuscript; available in PMC 2010 April 1.
Published in final edited form as:
PMCID: PMC2762351

Visualization, documentation, analysis, and communication of large scale gene regulatory networks


Genetic regulatory networks (GRNs) are complex, large-scale, and spatially and temporally distributed. These characteristics impose challenging demands on computational GRN modeling tools, and there is a need for custom modeling tools. In this paper, we report on our ongoing development of BioTapestry, an open source, freely available computational tool designed specifically for GRN modeling. We also outline our future development plans, and give some examples of current applications of BioTapestry.


As our understanding of genetic regulatory networks (GRNs) increases, ever more complex networks are studied. Ad-hoc ways of describing such networks (e.g. using generic drawing tools) are inefficient and inadequate, and there is an increasing need for specialized software.

In this paper, we present a software tool we have developed, called BioTapestry ( [1], which has been designed from the ground-up to model GRNs. BioTapestry is free, open source, and runs on all popular computer platforms. Below, we illustrate some of the ways in which BioTapestry facilitates GRN modeling. Since BioTapestry is an active and ongoing project, we will also outline the requirements that will guide our future development efforts. We follow this with a roadmap on how to get started using BioTapestry, and finally give some examples of current applications of the software.

The need for a specialized GRN modeling tool

The architecture of a GRN arises directly from the DNA sequence of the genome, and a GRN model is directly testable by DNA manipulations. Thus, the representation of GRNs must be genome oriented, with specific emphasis placed on the predicted DNA inputs that form the basis of the model. Furthermore, the GRN needs to be viewable at a number of different levels, from the whole, to the subcircuits, to the cis-regulatory DNA, and to the nucleotide sequence.

General-purpose network layout and presentation tools do not provide an appropriate level and style of abstraction for modeling GRNs. Many pathway modeling tools represent molecular interaction networks at the level of biochemical reactions. Because of the large number of reactions involved, representing GRNs as a set of biochemical reactions can result in overwhelmingly complex diagrams and obscure the regulatory architecture of GRNs. Moreover, the necessary biochemical data is rarely available to characterize such detailed views. On the other hand, overly abstract representations, such as those used in graph visualization software, lead to ambiguous network diagrams that convey little information.

Figure 1A shows a common graph-layout style GRN diagram. Compare this to the BioTapestry diagram in Figure 1B. The BioTapestry view immediately coveys a number of key concepts absent from the graph view. Firstly, cis-regulatory relationships are very easy to see and decipher in the BioTapestry view. Second, all nodes are not equal in BioTapestry layouts. BioTapestry uses automated layout templates to highlight regulatory relationships among the genes. For example in this figure, upstream regulators are placed near the top and to the left, while downstream genes are cascaded towards the right and bottom.

Figure 1
Part A: A regulatory network presented using a common graph-layout style GRN diagram. Part B: The same network rendered using BioTapestry. Also shown (inset, upper right) is a screenshot displaying a subnetwork view of the GRN at a specific time in a ...

A key feature of GRNs is that a single gene will typically perform different regulatory interactions in different cells and at different times. A single static view of a GRN cannot convey the way a gene becomes part of different processes and functional modules in different cells and times. Figure 1B also shows how BioTapestry provides a hierarchical representation of GRNs which allow a user to track a GRN within a given group of cells over time, or to compare GRN state different between different cells at any given time.

Finally, BioTapestry is designed to facilitate the process of GRN model building and provides extensive support for network annotation and curation. We discuss these features in depth below.

GRN-specific representation

BioTapestry supports a symbolic representation of genes, their products, and their interactions, which emphasizes regulatory and experimentally-derived network features.

Representing GRNs at the cis-regulatory level

The most important concept to be communicated in GRN visualization is how the transcription of a gene is regulated by other genes in the network. This crucial information must be instantly recognizable from a cursory inspection of the network diagram. The representation of a gene, and in particular the cis-regulatory region of the gene, must be unique, structured, and organized in a fashion that stands out and is quickly understandable.

BioTapestry depicts a gene with the commonly used shorthand representation shown in Figure 2A. The key feature is an explicit schematic representation of the cis-regulatory modules of the gene. As shown, multiple binding sites for the same transcription factor, and multiple cis-regulatory modules can be depicted. Potentially important regulatory features, such as the spatial ordering of transcription factor binding sites on DNA, are also preserved. Furthermore, each regulatory input can be provided with a colored annotation tag for documentation purposes. For example, we often use colored diamond symbols to indicate the type of experimental evidence available for each binding site.

Figure 2
BioTapestry's symbols and their use. Part A: The collection of symbols for genes and for representing off-DNA interactions and chemical species. Genes are shown in (1). The thick horizontal line is an abstract representation of the cis-regulatory region ...

The above “cartoon” representation facilitates whole-network visualization of GRNs. To allow more in-depth evaluation, BioTapestry provides a data page for each network element (e.g. a gene or an interaction). These data pages can be customized by the user to display tables, figures and other data, or illustrations of the internal cis-regulatory logic of a gene. We plan to enhance BioTapestry in the future to permit multiple drill-down displays, and to make it easier to customize the data page.

Compact representation of off-DNA interactions

BioTapestry depicts off-DNA interactions as simple, compact and distinct symbols that quickly provide a sense of the general nature of the process, and its regulatory inputs and outputs, while hiding the details. In this way, complex processes such as signal transduction are modeled in terms of their regulatory role within the GRN of interest. This makes it much easier to understand the GRN at a glance, and allows uncluttered visualization of large-scale GRNs.

BioTapestry's collection of symbols for off-DNA actions and interactions are shown in Figure 2A. They are designed to provide enough information so that a viewer can “mentally fill in” the details from general knowledge. For example the canonical Wnt pathway may be communicated by a single labeled symbol. In general, pathways that do not include multiple regulatory inputs can be summarized by a single input-output symbol to avoid clutter. In cases where the details do need be accessible to the user, right-clicking on the symbol allows the user to view the same type of (user-customizable) pop-up data page containing tables, figures and other data, as is provided for genes.

Simplified representation of post-transcriptional processes

Post translation processes that are not differentially regulated within a GRN of interest are not represented explicitly in BioTapestry. The outputs of transcription factor genes are typically shown as direct inputs into the regulated gene targets, with an implicit understanding of what that simplification represents. Sometimes, there are post-transcriptional steps that are critical components of the regulatory behavior of the network, e.g. translation inhibition. In these cases, an explicit series of one or more linked off-DNA symbols can be inserted into, and replace, the simple direct input link (see Figure 2B). Similarly, if the gene creates multiple distinct products that are relevant to the regulatory function, the single gene output can be split into several tagged links, one for each product (Figure 2C). The same approach can be used to model the regulation of gene outputs via interactions with micro-RNAs, such as destruction of mRNA or interference with translation. Transcription of miRNAs is regulated in the same fashion as any other genes and can be modeled in BioTapestry in the same way as for any other gene. Regulation by miRNAs is modeled as post-transcriptional interaction, as illustrated in Figure 2D.

Maximizing GRN readability

BioTapestry uses a variety of strategies to facilitate the visualization of large numbers of genetic linkages:

  • Links are not each drawn separately, but are bundled together and drawn as a group, as shown in Figure 3A. This is an extremely efficient approach that significantly cuts down on visual clutter that results when there are many genes in a network.
    Figure 3
    Features of links in BioTapestry. Part A: By bundling together all the separate links from the same source, BioTapestry can efficiently represent links and cur down visual clutter. Part B: Color is an essential component of link appearance in BioTapestry. ...
  • Coloring is used to distinguish between adjacent and overlapping lines. BioTapestry automatically assigns a color from a palette of visually distinct colors to each link source, and then uses that same color for all of the source's outbound links. This makes it much easier to visually tie a distant link source to its corresponding link traces while inspecting the network (Figure 3B).
  • Unique layout algorithms try to take advantage of the bundled link style to improve readability.
  • Interactive tools to find link sources and targets, and pop-up tool tips that identify link source names, help for quick identification of network connectivity.
  • “Branch bubbles” (as illustrated in Figure 3B) can be optionally activated to mark true link intersections, which can significantly improve readability and eliminate crossing ambiguities when the network grows large.

Although the outputs from a single source are typically rendered in a uniform fashion, BioTapestry allows the user to optionally highlight particular links by properties such as thickness, color, or line style. In addition, certain presentation properties of links, genes, and other nodes can be assigned to specific model properties. For example, the thickness of a link can be tied to the type of experimental evidence available. Such links are rendered accordingly even if the model is laid out differently (see Figure 3C).

Representing different network subsets and states

The same underlying GRN behaves differently in different cell types, spatial domains, and environmental conditions, and at different times. BioTapestry is designed to help the user to organize these varying views of the network state in a coherent fashion, while helping the user to understand how these views are derived from the single underlying GRN. As illustrated in Figure 4 (see also [1]), BioTapestry uses a three-level hierarchy to describe a GRN:

Figure 4
Hierarchical models in BioTapestry. Clockwise from upper left: View from the Genome (VfG), View from all Nuclei (VfA), View from the Nucleus (VfN). See text for further discussion of each model level.
  1. The View from the Genome (VfG) provides a summary of all inputs into each gene, regardless of when and where those inputs are relevant. Only one copy of each network element is shown.
  2. The View from All nuclei (VfA) is derived from the VfG, and contains the interactions present in different regions over the entire time period of interest. Each region in a VfA is a subset of the VfG, and sub-networks may be duplicated in different regions.
  3. Views from the Nucleus (VfN): Each VfN describes a specific state of the network at a particular time and place. Inactive portions of the network are indicated in gray, while the active elements are shown colored.

Each of these hierarchical views provides a different perspective on the GRN. A researcher can start their exploration of the network at any level, depending on the data available and the researcher's interests. For example, VfGs offer a natural perspective on each gene's full regulatory program within a GRN. However, to study functional motifs in the network (e.g. mutual exclusion), which are highly dependent on specific temporal and spatial conditions, VfN diagrams would be the most appropriate view.

BioTapestry makes it easy to create and organize GRN models using this approach. The hierarchical framework is general enough to be useful in many different ways. For example, one can use the lower levels of the hierarchy to depict variations in network behavior due to different experimental conditions. Alternately, submodels can be used to highlight network components discovered with particular experimental methods, or highlight subsets of genes or interactions that meet some significant selection criteria.

In BioTapestry, all the network models are automatically kept consistent across additions and deletions of network elements. For example, when a user inserts a node into an existing link in the root BioTapestry model, the program propagates the new node to all submodels that include an instance of the link. Since links may be laid out differently, and there may be multiple link copies spanning multiple regions in submodels, this can be a complex process and its automation greatly enhances model consistency and integrity. Along these same lines, we plan to add even more features that will help the user to easily propagate newly added network features to targeted regions and submodels.

Using the GRN model as an information portal

As we discussed above, network elements are simplified, abstract representations of complex subsystems, be they detailed cis-regulatory logic networks controlling a gene or complex off-DNA interactions such as lengthy signal transduction pathways. Often, it is useful to have access to more detailed information.

In BioTapestry, right-clicking on any gene or symbol in a network give the user an option to pop up a data display page. For genes, this page is typically configured to display raw perturbation data, the generic expression and interaction data tables that are used to drive the dynamic models, or arbitrary user-specified text. However, the page is customizable using small code plug-ins. In particular, it is easy to write a plug-in that displays data from a web server.

The plug-in approach allows for great flexibility in what data can be displayed, which is appropriate given the wide variety of information that may be appropriate to show. For example, what evidence is available to support the conclusion that a regulatory input is direct? What data have been used to determine that a gene expresses in a particular cell type and time? What are the exact cell type, genetic background and experimental conditions that a model is based on? Which particular member of a family does a gene or gene-product refer to? Potentially, all these types of information can be essential to fully understanding a model and its limitations.

In future versions of BioTapestry, we plan to simplify the data page customization process so that users can install commonly used options without needing to create special code modules.

Visualizing network dynamics

It is much easier to understand a sequence of GRN state changes through animations and interactive manipulations. A key feature of the BioTapestry Editor is that the user can click on genes and linkages to query their properties (e.g. all target genes, experimental evidence, or alternative paths between two nodes). These interactive features make it much easier to see the underlying organization within a large and busy static view.

BioTapestry also provides strong support for representing network dynamics:

  • Interactive time slider controls allow visualizations of sequences of GRN state changes. The visualizations are driven automatically from data tables, which can be built from either experimental data or simulation runs.
  • The user can create visualization sequence paths, each of which is a series of network views that can be easily traversed in the specified order. This allows the user to smoothly and quickly follow the state of a GRN through a sequence that may not be contiguous in the model hierarchy (see Figure 5).
    Figure 5
    User-defined model paths and intermediate expression levels are two new features of BioTapestry. The three screenshots show how a user-defined model viewing path can improve navigation. In many cases, following the temporal progress of a given region ...
  • Individual static GRN views can show intermediate levels of gene expression and link activity (not just binary on/off, fully colored or completely grey) in a variety of easy to understand ways. For example, variable expression levels can be depicted using pie glyphs, color saturation, and line thickness. The user can pick and choose between the desired representations (see Figure 5).

Support for continuously variable expression levels in dynamic time-slider models is a planned future enhancement.

Making interactive GRN models available on the web

BioTapestry is written in Java, which is a freely available, cross-platform web technology. The program can be run as a web-based application using the Java Web Start facility. This means that a user with Java installed on their computer can click on a link on a web page, and the program will be downloaded to their computer and start to run automatically. Using this framework, a stripped-down “read-only” version of the program, called the BioTapestry Viewer, allows web-based viewing and interactive exploration of published GRN models. This feature enables a GRN model to act as an interactive and dynamic information portal for dissemination of research results and classroom teaching. The facility also encourages the development of community consensus models for widely-studied GRNs, shared over the web and built up by long-range collaboration. Since it is possible to tie the annotation data displayed for each network feature to a web page (using the plug-in facility), collaborative web technologies such as wikis can be used to enable community discussion and feedback.

Presenting the GRN model within its spatial and temporal context

BioTapestry allows images to be displayed alongside each network view in the model hierarchy. This facility can be used to indicate the cells of interest in the embryo at the appropriate developmental stage, or to indicate other spatial and contextual information using illustrative cartoons. In this way, images can be used to provide biological context for GRN models (Figure 6). This is particularly useful in systems involving multiple signaling events and changing cellular neighborhoods (e.g. development). When combined with the viewing path feature described above, the smooth progression of GRN states can be tracked using the images as well as the network model.

Figure 6
Images can be displayed alongside each model, including each time point of a dynamic model, as shown here. The pictures provide the physical context for the abstract network depicted in the main window.

In future BioTapestry development, we plan to introduce 3D pictures and to integrate the image representation with the navigation functions.

Interacting with other computational tools

BioTapestry top-level networks can be exported to other software tools using the Systems Biology Markup Language (SBML) ( [2]. There are currently over 110 SBML-compatible software packages. BioTapestry also can import and export Cytoscape ( [3] interaction files. Furthermore, ongoing BioTapestry development aims to support the Gaggle framework ( [4]. Integration with the Gaggle will allow BioTapestry to interactively exchange network models and data with other tools.

A particularly popular feature of BioTapestry is that it can read entire hierarchical network descriptions from spreadsheets. Such spreadsheets could be generated manually, or by ad-hoc processing of experimental data, or automatically by computational data analysis pipelines. For example, Figure 7 shows a single subnetwork view of a very large dataset, from the Halobacterium EGRIN project ( [4], that was built using the BioTapestry spreadsheet import feature.

Figure 7
The improved comma-separated value (CSV) file import and automatic layout capabilities of BioTapestry were used to generate a large network, with an extensive set of detailed submodels such as the one shown here, for the EGRIN (Environment and Gene Regulatory ...

Data analysis and network construction

A common challenge in GRN modeling is distinguishing between direct and indirect linkages. To help meet this challenge, BioTapestry provides a display of alternative paths between a source and target gene (Figure 8A). The user can then inspect the supporting data to determine whether each linkage is direct. Tools like this help the researcher to build the model up from raw data and explore possible network architectures consistent with the data. We are currently working to enhance BioTapestry to provide more evidence visualization tools and a supporting analysis pipeline that can build a plausible network model from the perturbation and expression data.

Figure 8
BioTapestry network exploration tools. Part A: BioTapestry provides a display to show alternative paths between a source and target gene. Tools such as these help the researcher to investigate possible alternative explanations of perturbation data while ...

Network exploration

BioTapestry now includes an enhanced search tool (Figure 8B) that allows the user to search for all nodes matching or partially matching a given name. It also allows the user to find all sources or targets of a given node, optionally selecting the relevant network segments. Another tool allows the user to select the sources or targets of a given link segment. Even more search enhancements are anticipated. For example, we plan to modify the search tool to work across multiple models as well as the current model. Another planned feature will allow the user to select all nodes and links within n hops of a selected subnetwork, thus making it possible to find the local neighborhoods of that subnetwork.

Making it easy to update and maintain the network

In general, a GRN model will never be definitive and complete. Over time, the model will be expanded, pruned, revised, and refined; interactions that were once thought to be direct turn out to be indirect, and vice versa. Modifying and maintaining increasingly complex models over an extended period of time is challenging, and BioTapestry provides features to handle these issues. As one example, it provides an incremental layout feature so that new network elements can be added while retaining most of an existing layout. Since this is an important long-term goal, we are continuing to add improvements that simplify creation and modification of the network. For example, recent BioTapestry enhancements allow model hierarchy trees to be duplicated; regions in a model can also be duplicated. Nodes can be inserted into links and automatically propagated to submodels.

Providing extensive layout support customized to GRNs

Whether the layout of complex GRN models is meaningful and understandable is a subjective judgment dependent on factors such as the user's points of interest, and focus. The commonplace algorithms for laying out network diagrams typically do not take biological considerations into account in their layout. BioTapestry layout algorithms address this challenge in the following ways:

  • BioTapestry layout algorithms take into account the meaning of different symbol types, and use this information to organize nodes in an understandable way. For example, since genes are centrally important to understanding GRNs, they receive priority treatment in the layout process, with other symbols serving a supporting role (see Figure 3C). Furthermore, these small gene-centric blocks can be grouped into distinct units, and then these mid-scale units are laid out according to a high-level organizational strategy.
  • BioTapestry layout algorithms now use different layout templates to organize GRNs in pre-configured ways. For example, one template specifies that regulatory genes are located in a different region of the diagram from downstream target genes. Another template also separates regulatory genes from downstream target genes, and also specifies that genes are positioned such that spatial location is indicative of time of expression onset.
  • As was mentioned previously, the support of model hierarchies can make layout issues much more complex, since it is highly desirable to be able to maintain layout consistency among the various highly related but still unique views of the model. It is frequently important that the user be able to move between different representations of pieces of the network while still being able to tie that back to a mental model of the network as a whole. BioTapestry currently uses a network layout subsetting strategy for propagating an overall layout to submodels such that strict geometric ordering is maintained across views while compressing out unused space.
  • BioTapestry works to provide a variety of layout support tools to help a user who is drawing networks by hand, such as per-link auto layout, or global color assignment tools that can eliminate link crossing ambiguities.

Supporting multiple levels of abstraction within network views

The level of abstraction that we have used in BioTapestry so far has worked well for up to medium-sized networks, but as networks grow in size and complexity, new ways of organizing and thinking about network elements are needed.

The simplest way in which BioTapestry aids the understanding of GRNs is interactivity. A typical GRN presentation can be hard to understand when first viewed; it is only after interactively interrogating it and studying the various levels of the hierarchy that the network organization and functional features become apparent.

One of the ways in we hope we to facilitate understanding of large-scale GRNs is by introducing process diagrams as an additional level of representation within BioTapestry (see Figure 9A for an example). Another form of higher-level grouping, but typically dealing with smaller network chunks, is the identification of functional blocks [5] (for an example, see Figure 9B). Identification of functional blocks (e.g. feedback loops) also allows the user to view the GRN as a smaller set of interacting units, each with a clearly understood function.

Figure 9
Proposed methods for introducing higher level abstractions into BioTapestry. Part A: A way to integrate process diagrams into BioTapestry, applied to the PMC domain of the sea urchin endomesoderm network. From left to right: a process diagram (see [6 ...

These additional views, together with the network representations we have already discussed, imply a natural ordering of four levels of abstraction that are appropriate for looking at GRNs:

  1. High-level process diagrams
  2. Medium-level functional blocks
  3. Fine-grained gene-centric view
  4. Detailed network descriptions that model the logic embedded in cis-regulatory modules, and the specifics steps of off-DNA interactions

We plan to enhance BioTapestry to support these different representations in a manner that allows the user to view the network at the chosen level of detail, and to switch between these views as needed.

Getting started with BioTapestry

The BioTapestry Editor and associated tutorials are available freely at The only prerequisite is that the freely-available Java Runtime Environment (JRE), from Sun Microsystems, is installed on your computer. This is commonly the case. If not, JRE can be installed easily by following instructions on the BioTapestry home page. With Java installed, clicking on the BioTapestry start link will cause the Java Web Start system to download the software to your computer and run it. This system allows users to be kept up-to-date with the latest version. When a user is not connected to the internet, they can still run the program from a desktop icon using the associated files saved automatically when the BioTapestry Java Web Start link was last accessed. Note that although the software is downloaded and maintained via the web, your data stays on your machine as locally saved files and is never uploaded to the server.

Also on the BioTapestry home page are a set of online tutorials that are designed to help you learn how to use the program by taking you step-by-step through simple examples that highlight important and common program operations. For example, BioTapestry supports several different ways of creating networks. Firstly, GRN models can be created manually. In this case, BioTapestry provides the user with an easy way to draw the network by hand, with fine-grained manual control over network placement and link layout. This method is taught in a Quick Start tutorial. At the other end of the spectrum, for large-scale regulatory networks based upon large amounts of data (e.g. from high-throughput experiments or computational analyses), BioTapestry is able to import large network descriptions from spreadsheet files, and automatically layout the large number of network elements in a coherent fashion. This method is described in a separate tutorial on Building Networks from Comma-Separated Value Files. Between these two extremes, BioTapestry allows users to specify networks interactively through a set of dialog boxes that guide the user. BioTapestry then automatically generates the network layout, thereby avoiding mundane layout tasks; this method is covered in a tutorial on Building Networks from Interaction Tables

In addition to these and other tutorials, there is an extensive online Frequently Asked Questions (FAQ) list that covers many topics in considerable depth.

Recent applications of BioTapestry

BioTapestry was initially developed to model the sea urchin endomesoderm specification networks [8]. A regularly updated BioTapestry viewer for this network is available at:

It has since been used to model GRNs in a wide variety of organisms, including drosophila [9], arabidopsis [10], ciona [11], and yeast [12].

Also, an increasing number of projects are taking advantage of BioTapestry's ability to share interactive GRN models over the web. For example, an interactive model of mouse ventral neural tube specification [13] is available at:

This model demonstrates how multiple VfA models can be created in BioTapestry to track a continuously varying set of developmental domains as development progresses.

The EGRIN (Environment and Gene Regulatory Influence Network) model recently developed for Halobacterium salinarum [4] provides an example of using BioTapestry's CSV input facilities and new auto layout algorithms to visualize large networks generated from computational analysis of high-throughput data. See (Figure 7), and:

In this case, BioTapestry's model hierarchy provides an excellent solution for showing the many different states of a genome-scale model.

The zebrafish developmental GRN presented in this issue [14] is an example of a BioTapestry model developed by merging data from public databases with additional local and 3rd party experimental observations. An interactive web model is available at:

Finally, an interactive web model of the mammalian T-cell developmental GRN [18], which is available at: is an example of a BioTapestry model developed by merging data from local and 3rd party experimental observations with additional public sources.

Looking ahead

More and more GRNs are being characterized every day. The resulting models are increasingly complex, and integrate very large volumes of experimental data. As the size and complexity of GRN models grows, four inherent capabilities of BioTapestry will prove increasingly essential.

Firstly, BioTapestry network diagrams present an integrated view of (i) the high-level architecture of the network, (ii) the cis-regulatory features of individual genes, and (iii) the supporting experimental evidence.

Second, BioTapestry's hierarchical views of GRNs highlight regulatory differences among cells (visualized in VfNs), as well as regulatory changes over time (visualized with ‘slide shows’), while at the same time emphasizing the relationship of dynamic GRN modules to genomic organization (visualized in VfAs and VfGs).

Third, as larger and larger networks are studied, the size and complexity of datasets is increasingly making their ad-hoc interpretation difficult and error-prone. BioTapestry supports a structured process for curating and translating experimental data into GRN models.

Fourth, our ongoing work to integrate process diagrams, and functional block representations into BioTapestry models facilitates the comparison of GRNs from different organisms (see for example [15,16,17]), allowing new insights into the functional and logical architecture of GRNs, and the fundamental principles underlying genetic control of cellular function.


BioTapestry development is supported by NIGMS grant GM061005.


Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.


1. Longabaugh WJR, Davidson EH, Bolouri H. Computational representation of developmental genetic regulatory networks. Dev Biol. 2005;283:1–16. [PubMed]
2. Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novere N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J. The Systems Biology Markup Language (SBML): A Medium for Representation and Exchange of Biochemical Network Models. Bioinformatics. 2003;19(4):524–531. [PubMed]
3. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003;13:2498–2504. [PubMed]
4. Bonneau R, Facciotti MT, Reiss DJ, Schmid AK, Pan M, Kaur A, Thorsson V, Shannon P, Johnson MH, Bare JC, Longabaugh W, Vuthoori M, Whitehead K, Madar A, Suzuki L, Mori T, Chang D, DiRuggiero J, Johnson CH, Hood L, Baliga NS. A Predictive Model for Transcriptional Control of Physiology in a Free Living Cell. Cell. 2007;131:1354–1365. [PubMed]
5. Longabaugh W, Bolouri H. Understanding the Dynamic Behavior of Genetic Regulatory Networks by Functional Decomposition. Current Genomics. 2006;7(6):333–341. [PMC free article] [PubMed]
6. Oliveri P, Tu Q, Davidson EH. Global Regulatory Logic for Specification of an Embryonic Cell Lineage. Proc Natl Acad Sci USA. in press. [PubMed]
7. Oliveri P, Walton KD, Davidson EH, McClay DR. Repression of mesodermal fate by foxa, a key endoderm regulator of the sea urchin embryo. Development. 2006;133:4173–4181. [PubMed]
8. Davidson EH, Rast JP, Oliveri P, Ransick A, Calestani C, Yuh C, Minokawa T, Amore G, Hinman V, Arenas-Mena C, Otim O, Brown CT, Livi CB, Lee PY, Revilla R, Rust AG, Pan Z, Schilstra MJ, Clarke PJC, Arnone MI, Rowen L, Cameron RA, McClay DR, Hood L, Bolouri H. A Genomic Regulatory Network for Development. Science. 2002;295(5560):1669–1678. [PubMed]
9. Levine M, Davidson EH. Gene regulatory networks for development. Proc Natl Acad Sci USA. 2005;102:4936–4942. [PubMed]
10. Wellmer F, Alves-Ferreira M, Dubois A, Riechmann JL, Meyerowitz EM. Genome-Wide Analysis of Gene Expression during Early Arabidopsis Flower Development. PLoS Genetics. 2006;2(7):e117. [PMC free article] [PubMed]
11. Stathopoulos A, Levine M. Genomic regulatory networks and animal development. Dev Cell. 2005;9(4):449–62. [PubMed]
12. Smith JJ, Ramsey SA, Marelli M, Marzolf B, Hwang D, Saleem RA, Rachubinski RA, Aitchison JD. Transcriptional responses to fatty acid are coordinated by combinatorial control. Molecular Systems Biology. 2007;3:115. [PMC free article] [PubMed]
13. Vokes SA, Ji H, McCuine S, Tenzen T, Giles S, Zhong S, Longabaugh WJR, Davidson EH, Wong WH, McMahon AP. Genomic characterization of Gli-activator targets in sonic hedgehog-mediated neural patterning. Development. 2007;134:1977–1989. [PubMed]
14. Chan TM, Longabaugh W, Bolouri H, Chen HL, Tseng WF, Chao CH, Jang TH, Lin YI, Hsu HI, Huang CB, Hung SC, Yuh CH. Developmental Gene Regulatory Networks of the Zebrafish Embryo. this issue. [PubMed]
15. Oliveri P, Davidson EH. Built to Run, Not Fail. Science. 2007;315(5818):1510–1511. [PubMed]
16. Hinman VF, Davidson EH. Evolutionary plasticity of developmental gene regulatory network architecture. Proc Natl Acad Sci USA. 2007;104(49):19404–19409. [PubMed]
17. Davidson EH, Erwin DH. Gene Regulatory Networks and the Evolution of Animal Body Plans. Science. 2006;311(5762):796–800. [PubMed]
18. Georgescu C, Longabaugh W, Scripture-Adams D, David-Fung E, Yui M, Zarnegar M, Bolouri H, Rothenberg EV. A gene regulatory network armature for T lymphocyte specification. PNAS. 2008;105(51):20100–20105. [PubMed]