Open Access

BISON: bio-interface for the semi-global analysis of network patterns

  • Christopher Besemann1,
  • Anne Denton1Email author,
  • Nathan J Carr2 and
  • Birgit M Prüβ2
Source Code for Biology and Medicine20061:8

DOI: 10.1186/1751-0473-1-8

Received: 22 July 2006

Accepted: 29 November 2006

Published: 29 November 2006



The large amount of genomics data that have accumulated over the past decade require extensive data mining. However, the global nature of data mining, which includes pattern mining, poses difficulties for users who want to study specific questions in a more local environment. This creates a need for techniques that allow a localized analysis of globally determined patterns.


We developed a tool that determines and evaluates global patterns based on protein property and network information, while providing all the benefits of a perspective that is targeted at biologist users with specific goals and interests. Our tool uses our own data mining techniques, integrated into current visualization and navigation techniques. The functionality of the tool is discussed in the context of the transcriptional network of regulation in the enteric bacterium Escherichia coli. Two biological questions were asked: (i) Which functional categories of proteins (identified by hidden Markov models) are regulated by a regulator with a specific domain? (ii) Which regulators are involved in the regulation of proteins that contain a common hidden Markov model? Using these examples, we explain the gene-centered and pattern-centered analysis that the tool permits.


In summary, we have a tool that can be used for a wide variety of applications in biology, medicine, or agriculture. The pattern mining engine is global in the way that patterns are determined across the entire network. The tool still permits a localized analysis for users who want to analyze a subportion of the total network. We have named the tool BISON (B io-I nterface for the S emi-global analysis O f N etwork patterns).


Research on biological networks is a well-established part of bioinformatics [1]. Examples of biological networks include regulatory networks [2], protein-protein interactions [3, 4], and domain-fusion networks [5, 6], among others. Typical objectives are to gain information about the over-all structure and evolution of the network in question [7, 8]. Protein function and other annotations are rarely included in network studies and, if so, the results are normally limited to the statistics of similarity or dissimilarity between neighbors [1] or a correlation of function with traditional subgraph statistics [9]. Note that in the following, we will use the terms network and graph interchangeably. We will sometimes refer to proteins in a network as nodes of a graph, and to regulatory or other interactions as edges.

The large amount of annotation and network data that has accumulated over the past decades requires the use of data mining techniques. Pattern mining is a subset of data mining that has the goal of identifying frequently occurring combinations of items of information. We will refer to pieces of information, such as domain and functional information, as properties. Initial work considered simple types of item information [10, 11]. Pattern mining techniques have also been used to find frequent subgraphs of larger graphs [12, 13]. The most general case of pattern mining considers any combination of relational tables [14]. Recently, the specific problem of finding patterns that involve networks and item data has gained importance [1517].

BISON integrates our own pattern mining techniques with modern graph visualization and navigation techniques. Combinations of visualization and navigation techniques have been used previously [1825]. Graph visualization techniques address complexity and size of networks [26].

We demonstrate the usefulness of BISON through two examples within the E. coli network of transcriptional regulation. The first example uses FlhD/FlhC, a transcriptional regulator that was originally described as an activator of more than 50 flagellar genes [27] and later recognized as a global regulator of metabolism [28]. Expression of the flhD operon is a target point for many global regulators and global signals [2932]. The portion of the E. coli transcriptional network that centers around FlhD/FlhC was summarized [33]. We will use this system to demonstrate how diverse data such as microarray data can be integrated with existing data and analyzed by BISON in the context of the entire regulatory network.

The second example focuses on ABC transporters, protein complexes that form continuous channels through both cellular membranes that are specific for certain substrates and require the hydrolysis of ATP to provide energy for the transport process (for a review, please, see [34]). Different regulators have been described for the regulation of the many ABC transporters [35]. To our knowledge, this study is the first attempt to summarize their regulations.

In summary, we take the study of biological networks beyond its traditional focus on network structure and move it towards a more function-oriented view that looks at meaningful patterns in a localized context and provides targeted information to biologists working on a limited number of genes.

Results and discussion

This study presents an application [BISON; see Additional file 1] that combines our global pattern mining engine (an extension of [17]) with modern navigation and network visualization techniques [36, 37]. Fig. 1 is a schematic of BISON. The underlying pattern mining engine is shown in the top portion of the Figure. It operates on the full network in a global fashion. The bottom part describes the BISON interface including a network visualization unit that uses graph navigation capability and navigation capability using a modern graphical user interface.
Figure 1

Schematic of the BISON data flow and design. Top part, underlying pattern mining engine; botton part, BISON interface.

BISON connects regulatory network information with properties of the involved proteins. Properties can be any information that is associated with proteins, including experimental information and sequence information, such as annotations. The majority of the property data that are currently included in BISON are hidden Markov model domains (HMMs). Fig. 2A presents as small portion of such a regulatory network with the respective HMMs of the involved proteins.
Figure 2

Concept of a graph-relational pattern. Panel A: Portion of a regulatory network with property data. Panel B, pattern from that network.

We refer to the global patterns as graph-relational patterns. Fig. 2B illustrates the concept of a graph-relational pattern. Patterns are defined as frequently occurring combinations of properties in regulator proteins and the proteins that are encoded by their regulated genes. Patterns express structural attributes of the graph as well as specific properties of the relational records. Consequently, they form a bridge between regulation and function. Throughout this manuscript, properties for regulators are indicated (0) and properties for regulated proteins are indicated (1). The pattern in Fig. 2B, {(0).(hmm.response_reg),(1).(hmm.abc_tran)}, occurs at least five times within the entire network, a portion of which is displayed as Panel A.

It is important to note that BISON determines the frequency of occurrence of a certain pattern across the entire network and not only over the portion of the network that is displayed. Before a combination is considered a pattern, the combination has to occur at least a user specified number of times. The pattern repository currently provided with BISON is based on a cutoff of five times. It is also noteworthy that the pattern search is not limited to single properties on each site. All combinations of properties that satisfy the cutoff are being tested. While the current sub-graph type is set at single interactions (two proteins), the application also has the ability to find patterns that involve more than two proteins. Examples of these larger patterns are the hierarchical regulation of a final target gene, the simultaneous regulation of several genes with distinct properties, and co-regulation of a target gene by multiple regulators. Instructions on how to change these settings are provided with the User Manual that accompanies this manuscript [see Additional file 2].

Fig. 3 further details the concept of a graph-relational pattern. The combination of the network data, the property data, and the shape of the pattern (Panel A) will yield the combinations of genes that eventually form the pattern. Panel B includes a sampling of property combinations that have the potential of leading to patterns. The two columns labeled 'Descriptors' contain the property information. The two columns labeled 'Gene' contain the network information. Each gene encoding a regulator in the 'Gene 0' column is listed with the regulated gene in the 'Gene 1' column. The first line is a combination that contributes to the pattern. In this example, the pattern has a single property on each side. The second and third lines show how a pattern with multiple properties may be formed. The fourth line shows how gene ontology (GO) information can contribute to potential patterns. Of all the potential property combinations that are listed in Panel B, only the combination of the hmm.response_reg domain with the hmm.abc_tran domain (printed in bold) actually led to the pattern that is displayed in Panel C. The Pattern Information Page in BISON only lists such combinations of properties that actually form the pattern.
Figure 3

Details of a graph-relational pattern. Panel A: Network data, property data and the shape of the network define the combination of properties that can form a pattern. Panel B: Sampling of property combinations that have the potential to form patterns. Panel C: Pattern that was formed by the combinations (bold print) in Panel B.

We now provide two examples of the type of biological questions that the pattern search engine of BISON can answer. Our examples will use the transcriptional network of the Escherichia coli K-12 strain. The first question asked is what functional categories of proteins are regulated by the global regulator FlhD/FlhC. This is an important question when analyzing microarray data, avoiding the publication of long lists of genes. For this purpose, the gene list of a previous microarray experiment [28] was integrated into BISON. Instructions on how users can integrate their own microarray data into BISON are included in the User Manual.

Using the pattern-centered analysis function of BISON (Fig. 4), we were able to determine all patterns that contain the hmm.flhd domain in the regulator protein (left portion of the navigation page), the genes that are associated with each pattern (right portion of the navigation page), and the full pattern information for each pattern (pattern information page).
Figure 4

Pattern-centered screen shot of BISON. Panel A, Pattern Information page; Panel B, Network Visualization page; Panel C, Navigation page.

Navigation page (Fig. 4C): The part of the 'Descriptors' column, designated (0), indicates properties found in the regulator, while the part, designated (1), indicates properties found in the regulated proteins. The 'Links' column indicates whether regulation is positive (+) or negative (-). Selecting a pattern on the left side will yield the associated genes on the right side. For example, the hmm.flhd domain frequently occurred in combination with the flg__bb_rod domain. The FlhD/FlhC regulated proteins that contain this domain are FlgB, FlgC, FlgE, FlgF, and FlgG.

Network visualization page (Fig. 4B): Understanding of the patterns is supported by a visual display of the network environment. This display is limited to those proteins that have an immediate regulatory interaction with this gene to keep the computational and visual complexity under control. Selecting a gene on the right side of the navigation page will form the network around this gene in the network visualization page. This will also perform a switch from pattern-centered analysis to gene-centered analysis by converting the pattern information page into the gene information page (Fig. 5). This function is designed for users who want to further investigate properties of a gene or protein. A detailed list of genes and protein functions is now provided in Panel A. A link is provided to an external information source which is currently the Kyoto Encyclopedia of Genes and Genomes [38, 39].
Figure 5

Gene-centered screen shot of BISON. Panel A, Gene Information page; Panel B, Network Visualization page; Panel C, Navigation page.

Pattern information page (Fig. 4A): This page contains a summary of all the information that is required for the pattern analysis. The columns are designated as described for Fig. 3B. The two 'Descriptors' columns contain the properties, the two 'Gene' columns the network information. In our example, FlhD/FlhC with the hmm.flhd domain regulates all the genes whose encoded proteins contain the hmm.flg_bb_rod domain (flgC, flgD, flgE, flgF, and flgG). This information can be obtained for all the patterns in the left portion of the navigation page.

Table 1 summarizes a selection of the results from the FlhD/FlhC analysis. Starting from Fig. 4, HMMs from the regulated proteins were taken from the 'Descriptors' column (labeled (1)) and regulated genes whose proteins contain this HMM were taken from the 'Gene 1' column. This was done for a selection of the patterns. For the sake of an easier discussion, proteins containing related HMMs are grouped into functional categories. For example, proteins that contain the hmm.flg_bb_rod domain and proteins that contain the hmm.mcpsignal domain are discussed together as flagellar proteins. Proteins that contain the hmm.pts_eiic domain and proteins that contain the hmm.abc_tran domain are discussed together as transport proteins. The terminal reductases of anaerobic metabolism contain up to three specific domains. These, too, are discussed together. Functional categories are indicated as bold printed subheadings in Table 1.
Table 1

Functional categories of proteins whose genes are regulated by FlhD/FlhC

Regulator gene1

HMM of the regulator2

Regulated genes3

HMM of regulated proteins4

Flagellar proteins




flgB, flgC, flgE, fleF, flgG




tar, tsr, tap, trg, aer


Anaerobic metabolism




dmsA, ynfE, ynfF, napA, narG, narZ

molydopterin_oxidoreductase, molydop_binding, molybdopterin

Transport proteins




nagE, frvB, treB, murP, mngA, chbC, malX, fruA, ptsG




rbsA, oppD, nikD, fepC, proV, nikE, btuD, uup, metN, thiQ, phnK, msbA, fhuC, mglA, ybhF, dppD


Transcriptional regulators




cynR, dsdC, xapR, tdcA, ydcI, ybbS, nhaR, argP, yeeV

lysr_substrate, hth_1

1The gene 'flhD'.

2HMM for FlhD.

3Genes that are regulated by FlhD/FlhC are taken from the 'Gene 1' column in BISON that is associated with a selection of the property combinations in the 'Descriptors' column.

4HMMs for the regulated proteins are taken from the 'Descriptors' column in BISON.

Two functional categories were chosen for discussion as control groups of genes whose regulation by FlhD/FlhC was known. The remaining two functional categories were chosen as examples for functional categories that were newly recognized as regulated by FlhD/FlhC.

The two functional categories of proteins that were known to be regulated by FlhD/FlhC are flagellar proteins and anaerobic metabolism [28]. The group of flagella proteins that are covered with this analysis, includes the above mentioned flagella basal body rod proteins (hmm.flg_bb_rod) and the methyl-accepting chemotaxis proteins (hmm.mcpsignal). Since FlhD/FlhC was initially described as a transcriptional regulator of all flagellar genes [27], this confirms previous results and will serve as our primary positive control. Please, note that FlhD/FlhC regulates more flagellar genes than just the ten genes that are listed in Table 1. A complete list of FlhD/FlhC regulates genes can be obtained with the gene-centered analysis (gene information page). Pattern-centered analysis can lead to further hypotheses on domains that co-regulated genes share. Lowering the cutoff from five to three (number of combinations that are considered a pattern) yields the hmm.duf domain within FlgK, FlgG, and FlgE. This is so far a domain of unknown function. The observation that it is common among FlhD-regulated proteins may assist in determining its associated function.

As a secondary positive control, we confirmed the effect of FlhD/FlhC upon enzymes of anaerobic metabolism. As published before [28], several proteins that contain the HMMs characteristic for terminal reductases were regulated by FlhD/FlhC. Interestingly, BISON identified two proteins of this category (YnfE, YnfF: [40]) as regulated by FlhD/FlhC that were missed by the previous analysis [28].

In addition to the controls mentioned above, we were able to establish new information using BISON. Two functional categories of genes regulated by the regulator of interest, FlhD/FlhC, were better defined or newly established, transporters and transcriptional regulators. While it was previously recognized [28] that FlhD/FlhC regulated many genes encoding transporters, no attempt had been undertaken to classify these transporters. Our current analysis revealed two HMMs within the functional category of transport proteins that were regulated by FlhD/FlhC. These are pts_eiic (phosphotransferases) and abc_tran (ABC transporters). Each HMM is found in many proteins whose genes are regulated by FlhD/FlhC (Table 1). The HMM lysr_substrate in combination with hth_1 (helix-turn-helix) is an indication for DNA binding and transcriptional regulation. Nine proteins contain these two HMMs, all of these are regulated by FlhD/FlhC. This functional group of transcriptional regulators is a very interesting new finding, supporting our idea of FlhD/FlhC being part of a larger network of transcriptional regulation [33].

In summary, BISON identified two functional categories that were previously described as regulated by FlhD/FlhC (flagellar proteins and anaerobic metabolism), further detailed the category of transport proteins, and identified a new functional category, transcriptional regulators. In addition, new proteins were recognized as regulated by FlhD/FlhC in the category that relates to anaerobic respiration. We interpret this as an indication that pattern mining is a useful tool for the analysis of complex microarray and network data.

The second biological question that was asked was what functional categories of proteins (indicated by their HMM) are involved in the regulation of ABC transporter genes. A search was performed for patterns containing 'ABC'. After selecting a pattern involving 'ABC', the hmm.abc_tran domain appears in the 'Descriptor' column labeled (1). HMMs of proteins that regulate ABC transporter genes were taken from the 'Descriptor' column, labeled (0). The genes encoding these regulators were taken from the 'Gene 0' column. The specific ABC transporter genes that are regulated by this regulator were taken from the 'Gene 1' column. These regulations are summarized in Table 2 and correspond to the Pattern Information page.
Table 2

Regulators that affect the expression levels of ABC transporter genes

Regulator gene1

HMM of the regulator2

Regulated genes3

HMM of regulated proteins4

Two-component systems








ugpC, pstB, phnL, phnK, phnC




cydC, cydD, ccmA








tauB, fepC, ycjV, nikD, nikE, malK







hiska, response_reg

tauB, fepC, ycjV, nikD, nikE, malK




yehX, cysA, proV, ugpC, dppD




yehX, cysA, proV, ugpC, dppD




glnQ, potG, dppF, dppD, hisP, yhdZ


DNA binding proteins




modC, ccmA




tauB, cysA










1Genes encoding regulators that affect the expression levels of ABC transporter genes are taken from the 'Gene 0' column on the navigation page.

2HMMs for the regulators are taken from the 'Descriptors' column on the navigation page. HMMs are indicative of the functional categories of regulators (bold subheadings) that affect the expression levels of ABC transporters.

3Genes encoding ABC transporters that are affected by the regulators in column 2 are taken from the 'Gene 1' column on the navigation page.

4HMMs for the regulated proteins are taken from the 'Descriptors' column on the navigation page. HMMs are indicative of ABC transporters.

Proteins regulating ABC transporters were grouped into functional categories again, based upon their properties (HMM). Two-component systems each consist of a histidine kinase (hmm.hiska) and a response regulator (hmm.response_reg) (for a review, please, see [41]). RcsCDB is a rare case of a three-component system, where the first component (RcsC) contains both functional domains [42]. In the cases of RcsCDB, EnvZ/OmpR, and TorS/TorR, histidine kinases and response regulators have been identified as regulators of ABC transporter genes (Table 2). In the cases of PhoR/PhoB, NarX/NarL, NarQ/NarP, and NtrB/NtrC (synonym GlnG that was identified as a regulator of ABC transporter genes by its hth_8 domain), only the response regulator was found to regulate ABC transporter genes. Overall, the contribution of two-component systems to the regulation of ABC transporters seems to be large. Considering the small degree of overlap between the regulated genes, it seems like many two-component systems are specific for a certain set of ABC transporter genes.

An example of how BISON can be used to create hypotheses that can be further examined experimentally is given as an extension of the above study and involves the second functional categories of proteins that regulate ABC transporter genes, DNA binding proteins: ModE (hmm.hth_1) is a known repressor of the modA operon that encodes a molybdate specific transporter [43] and the ccmA operon that encodes a haem transport system [44]. While early studies with the ccmA operon showed that transcription was induced during anaerobic growth, regulation by known regulators of anaerobic respiration (FNR, ArcB/ArcA) could not be detected [44]. With this study, we found a regulation of ccmA by NarX/NarL and NarQ/NarP (Table 2). Both these two-component systems are global regulators during anaerobic growth in the presence of nitrate [45]. This leads to the hypothesis that regulation of ccmA by these systems might explain the regulation by anaerobiosis. Further experimentation might confirm this hypothesis.

In summary, regulation of ABC transporters seems to be as global as any regulation of metabolism. Interactions between regulators and their target genes are summarized in Fig. 6.
Figure 6

Schematic of the regulation of ABC transporters. Regulators are boxed, regulated ABC transporter genes are in italics.


We have developed a tool for the analysis of networks and global patterns. BISON is global in its design, yet allows the biology user to ask specific questions. We provide two examples for the kind of questions that can be answered with BISON and examine results from a microarray experiment in the context of the existing regulatory network. We also provide an example of how BISON can be used to create hypotheses. BISON provides a function oriented analysis for complex data. It bridges a gap between studies of single interactions and global pattern mining techniques. For this reason, we call our analysis 'semi-global'. Access to BISON is provided through BioMed Central [Additional file 1] and the NDSU Computer Science website [46].


Components of BISON

BISON is a tool designed to support exploration and graphical visualization of network patterns. It was developed in the Java 5.0 programming language. The schematic in Fig. 1 illustrates the structure of BISON. BISON takes raw network and property input as well as results provided by our own external Perl library for frequent (graph-relational) pattern mining (SGR Perl Library). We consider differential patterns that express differences in the properties of connected nodes as discussed [17]. In addition, patterns that represent similarity across edges are also discovered. Similarity and dissimilarity patterns are treated separately. This allows us to maintain the expression power of patterns we developed earlier [17] and yet to provide a larger set of patterns to the user. BISON is divided into a graphical visualization part and a text database interface for analysis and search. The network visualization component makes use of the Java Universal Network/Graph (JUNG) Framework [37]. The source code for BISON, documentation, and the pattern library are contained in the supplemental files.

Sources for network data

Network data were obtained from public databases, as well as experiments of our own and other researchers. A summary is given in Table 3. Four sources are described as follows: RegulonDB [2, 47] provides 2,537 interactions defined by 142 regulators and 1,059 targets. Microarray data for all of the E. coli two-component systems [48] contain 1,028 interactions between 40 regulators and 372 targets. A previous compilation of data [33] contains 1,969 interactions over 26 regulators and 856 targets. Finally, a previous microarray experiment from our own laboratory [28] contributed 886 interactions over 2 regulators and 444 targets. The total number of interactions is 6,227 between 186 regulators and 1,934 targets. These sources were formatted for use in BISON. Since we wanted to link network data to node property information, we resolved all gene names to a common identifier (Blattner IDs [49]). Gene names that could not be linked to Blattner IDs were excluded from the network since they could not contribute to useful patterns.
Table 3

Network data




Regulated genes


File name


























Sources for property data

Property data were obtained from public sources, supplemented by our own set of properties. A summary is given in Table 4. The initial property data were obtained from the E. coli genome project [49, 50] that includes 4,285 proteins and contains Gene Ontology (GO) classifications. The project yielded a total of 62 GO annotations (prefix "GO" in BISON). To extend these data, we searched the protein sequences for protein family domains. First, we gathered default data from the Pfam site [51, 52] from which we extracted 1,032 annotations covering 2,271 proteins (prefix "PF" in BISON). In addition, we used the HMMER profile hidden Markov model software [53, 54] that is used by Pfam to identify potential domain annotations using different criteria. Sequences from the E. coli Genome Project [49, 50] were tested with the Pfam A domain models. Domains were detected at a cut-off e-value of 1e-10. From this software setup, we extracted 1,747 annotations covering 3,124 proteins (prefix "hmm"). Note that some proteins received multiple HMM annotations and some annotations overlap with the default Pfam annotations.
Table 4

Annotation data






E. coli Genome Project



















Data input files

All input data are collected in a single data directory (default_data). The directory includes a configuration file (bison.config) read by BISON that specifies names of data files. Files required by BISON: An entity file (ecoli_entity.txt) lists the nodes of the network and the set of properties available for each node. Specific node IDs (Blattner IDs) are gathered from this file. An alias file (ecoli_alias.txt) specifies the default gene names for the nodes. A synonym file (ecoli_syn.txt) lists additional names for the nodes. Finally, a pattern file (patterns.out) stores the patterns of entities and annotations discovered in the network. This file is generated by the pattern mining engine. BISON can accept results from other pattern generation libraries, provided the format requirements of the pattern output file are satisfied. The configuration file also contains a list of external web links that are used in combination with protein IDs to construct hyperlinks in the object information page.

Four network files (*.net) contain data from the four different sources we used (Table 3). Additional network files can be added to allow the user to integrate their own data and understand them in the context of the existing network. A line will have to be added in the bison.config file that lists the name of the new data file. Similarly, protein property information can be supplemented. If additional data are to be included in the pattern discovery, then the pattern mining script has to be re-run with the new data. This will result in a new pattern file (see above). Instructions for the addition of user-specific files are contained in the User Manual.

JUNG Java extensions

We extended the JUNG Java library [37] as a basis for the graph visualization component. By default, the Fruchterman-Reingold node layout [36] is used for graph visualization if the displayed portion of the graph is smaller than 100 nodes, and a circle layout is used otherwise. BISON also resizes nodes to fit into the visualization space. Since this can cause problems with reading labels on the graph, text information is also provided in the application.

Pattern mining

Graph-Relational pattern analysis (Fig. 2, top portion, indicated as SGR) is implemented as a Java library within BISON. A summary of the mining process is outlined below:
  1. 1.

    Input network data and property data

  2. 2.

    Join the input data according to the requested pattern shape

  3. 3.

    Find frequent patterns in the joined table of data

  4. 4.

    Compute statistical significance measure

  5. 5.

    Output patterns for use in BISON


Step 1: Input files for the pattern mining contain the network data and the property data. The goal of pattern mining is to characterize interaction structures given as pattern shapes or subgraphs and the properties associated with the nodes involved in those structures. Pattern shapes can be pairs of regulators and regulated genes, or regulatory structures of three or more proteins.

Steps 2 and 3: The algorithm is an extension of our previous work [17]. In addition to analyzing differences between nodes, this algorithm allows the study of similarities. Patterns are separated into groups that represent combinations of similar and dissimilar annotations between different nodes in the network structure.

Step 4: The p-value for frequent patterns from step 3 is calculated, using a Chi-Squared test based on the contingency table of all items in the pattern. The objective of the test is to determine whether the pattern could have occurred randomly. We consider independence of all items in the pattern as null hypothesis. Note that both absence and presence of properties defines a pattern (see [17] for a discussion). That means that for a pattern (0).A (1).B in a 1-edge shape a four-dimensional contingency table is constructed based on variables 0.A, 1.A, 0.B, 1.B. Note also that a pattern may appear as significant because of co-occurrence of items within individual proteins. A p-value of 0.001 was chosen as cutoff. Based on this cutoff, all patterns that were discussed in the manuscript are significant. We are currently working on algorithms to specifically determine the significance of patterns with respect to the regulatory interactions and take the sparseness of the data into account.

Step 5: Writes the results to an output file (patterns.out) that is compatible with BISON.

Publication of the BISON application as part of an Open Access journal and making the code available in an Open Source format helps to facilitate the accessibility of the techniques. We will deliver updates and improvements as the BISON application progresses. We appreciate citation of this publication for users who want to include results obtained with BISON in their manuscripts.



Bio-Interface for the Semi-global analysis Of Network patterns


Hidden Markov Model


Java Universal Network/Graph Framework


K yoto E ncyclopedia of Genes and G enomes


S ignificant G raph-R elational pattern mining library.



The authors thank Shelley Horne, Megan Townsend, and Preeti Sule (all North Dakota State University, Fargo ND) for testing BISON and making useful suggestions. We thank both of our anonymous reviewers for their help improving the manuscript and the BISON interface.

AD and CB were supported by grant IDM-0415190 from the National Science Foundation. BMP was supported by the EPSCoR (through grant #EPS-0132289 from the National Science Foundation) program of North Dakota and a 'Biosecurity, Disease Surveillance, and Food Safety' earmark grant (through USDA-APHIS). CB and NC were supported by the 'Small Grants Program to foster Research Collaboration with the ND Agricultural Experiment Station' from North Dakota State University. CB was also supported in part by the NDSU Presidential Fellowship program.

Authors’ Affiliations

Department of Computer Sciences, North Dakota State University
Department of Veterinary and Microbiological Sciences, North Dakota State University


  1. Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, Chen Y, Cheng X, Chua G, Friesen H, Goldberg DS, Haynes J, Humphries C, He G, Hussein S, Ke L, Krogan N, Li Z, Levinson JN, Lu H, Menard P, Munyana C, Parsons AB, Ryan O, Tonikian R, Roberts T, Sdicu AM, Shapiro J, Sheikh B, Suter B, Wong SL, Zhang LV, Zhu H, Burd CG, Munro S, Sander C, Rine J, Greenblatt J, Peter M, Bretscher A, Bell G, Roth FP, Brown GW, Andrews B, Bussey H, Boone C: Global mapping of the yeast genetic interaction network. Science. 2004, 303: 808-813. 10.1126/science.1091317.View ArticlePubMed
  2. Salgado H, Gama-Castro S, Peralta-Gil M, az-Peredo E, Sanchez-Solano F, Santos-Zavaleta A, Martinez-Flores I, Jimenez-Jacinto V, Bonavides-Martinez C, Segura-Salazar J, Martinez-Antonio A, Collado-Vides J: RegulonDB (version 5.0): Escherichia coli K-12 transcriptional regulatory network, operon organization, and growth conditions. Nucl Acids Res. 2006, 34: D394-D397. 10.1093/nar/gkj156.PubMed CentralView ArticlePubMed
  3. Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y: A comprehensive two-hybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci U S A. 2001, 98: 4569-4574. 10.1073/pnas.061034498.PubMed CentralView ArticlePubMed
  4. Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, Qureshi-Emili A, Li Y, Godwin B, Conover D, Kalbfleisch T, Vijayadamodar G, Yang M, Johnston M, Fields S, Rothberg JM: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000, 403: 623-627. 10.1038/35001009.View ArticlePubMed
  5. Truong K, Ikura M: Domain fusion analysis by applying relational algebra to protein sequence and domain databases. BMC Bioinformatics. 2003, 4: 16-10.1186/1471-2105-4-16.PubMed CentralView ArticlePubMed
  6. Marcotte EM, Pellegrini M, Ng HL, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285: 751-753. 10.1126/science.285.5428.751.View ArticlePubMed
  7. Wagner A, Wu W: Compactness and cycles in signal transduction and transcriptional regulation networks: a signature of natural selection?. Adv Complex Syst. 2005, 7: 419-432. 10.1142/S0219525904000251.View Article
  8. Przulj N, Corneil DG, Jurisica I: Modeling interactome: scale-free or geometric?. Bioinformatics. 2004, 20: 3508-3515. 10.1093/bioinformatics/bth436.View ArticlePubMed
  9. Yu H, Zhu X, Greenbaum D, Karro J, Gerstein M: TopNet: a tool for comparing biological sub-networks, correlating protein properties with topological statistics. Nucleic Acids Res. 2004, 32: 328-337. 10.1093/nar/gkh164.PubMed CentralView ArticlePubMed
  10. Agrawal R, Imielinski T, Swami AN: Mining association rules between sets of items in large databases: 1993/5/26. Washington, D.C.; 1993:207-216.
  11. Agrawal R, Srikant R: Fast algorithms for mining association rules in large databases: 1994/9/12. Santiago, Chile; 1994:487-499
  12. Kuramochi M, Karypis G: Frequent subgraph discovery: 2001/11/29. San Jose, CA; 2001:313-320.
  13. Yan X, Han J: gSpan: graph-based substructure pattern mining: 2002. Maebashi City, Japan; 2002.
  14. Goethals B, Van den Bussche J: Relational association rules: getting WARMeR: 2002/9/16; London, UK. London, UK; 2002:125-139.
  15. Oyama T, Kitano K, Satou K, Ito T: Extraction of knowledge on protein-protein interaction by association rule discovery. Bioinformatics. 2002, 18: 705-714. 10.1093/bioinformatics/18.5.705.View ArticlePubMed
  16. Besemann C, Denton A, Yekkirala A, Hutchison R, Anderson M: Differential association rules: understanding annotations in protein interaction networks. Advanced Data Mining Technologies in Bioinformatics. Edited by: Hsu HH. 2006, Tamkang University, Taiwan, Idea Group Publishing
  17. Besemann C, Denton A, Yekkirala A, Hutchison R, Anderson M: Differential association rule mining for the study of protein-protein interaction networks. Edited by: Zaki MJ, Morishita S and Rigoutsos I. Seattle, Washington; 2004:72-80.
  18. Osprey. 2006, [http://​biodata.​mshri.​on.​ca/​osprey/​servlet/​Index]
  19. VisANT. 2006, [http://​visant.​bu.​edu/​]
  20. Cytoscape. 2006, [http://​www.​cytoscape.​org/​]
  21. Aragues R, Jaeggi D, Oliva B: PIANA: protein interactions and network analysis. Bioinformatics. 2006, 22: 1015-1017. 10.1093/bioinformatics/btl072.View ArticlePubMed
  22. Enright AJ, Ouzounis CA: BioLayout - an automatic graph layout algorithm for similarity visualization. Bioinformatics. 2001, 17: 853-854. 10.1093/bioinformatics/17.9.853.View ArticlePubMed
  23. Hanisch D, Sohler F, Zimmer R: ToPNet--an application for interactive analysis of expression data and biological networks. Bioinformatics. 2004, 20: 1470-1471. 10.1093/bioinformatics/bth096.View ArticlePubMed
  24. Iragne F, Nikolski M, Mathieu B, Auber D, Sherman D: ProViz: protein interaction visualization and exploration. Bioinformatics. 2005, 21: 272-274. 10.1093/bioinformatics/bth494.View ArticlePubMed
  25. Kolpakov FA, Ananko EA, Kolesov GB, Kolchanov NA: GeneNet: a gene network database and its automated visualization. Bioinformatics. 1998, 14: 529-537. 10.1093/bioinformatics/14.6.529.View ArticlePubMed
  26. Ju BH, Han K: Complexity management in visualizing protein interaction networks. Bioinformatics. 2003, 19 Suppl 1: i177-i179. 10.1093/bioinformatics/btg1022.View ArticlePubMed
  27. Bartlett DH, Frantz BB, Matsumura P: Flagellar transcriptional activators FlbB and FlaI: gene sequences and 5' consensus sequences of operons under FlbB and FlaI control. J Bacteriol. 1988, 170: 1575-1581.PubMed CentralPubMed
  28. Prüß BM, Campbell JW, Van Dyk TK, Zhu C, Kogan Y, Matsumura P: FlhD/FlhC is a regulator of anaerobic respiration and the Entner-Doudoroff pathway through induction of the methyl-accepting chemotaxis protein Aer. J Bacteriol. 2003, 185: 534-543. 10.1128/JB.185.2.534-543.2003.View ArticlePubMed
  29. Bertin P, Terao E, Lee EH, Lejeune P, Colson C, Danchin A, Collatz E: The H-NS protein is involved in the biogenesis of flagella in Escherichia coli. J Bacteriol. 1994, 176: 5537-5540.PubMed CentralPubMed
  30. Francez-Charlot A, Laugel B, Van GA, Dubarry N, Wiorowski F, Castanie-Cornet MP, Gutierrez C, Cam K: RcsCDB His-Asp phosphorelay system negatively regulates the flhDC operon in Escherichia coli. Mol Microbiol. 2003, 49: 823-832. 10.1046/j.1365-2958.2003.03601.x.View ArticlePubMed
  31. Lehnen D, Blumer C, Polen T, Wackwitz B, Wendisch VF, Unden G: LrhA as a new transcriptional key regulator of flagella, motility and chemotaxis genes in Escherichia coli. Mol Microbiol. 2002, 45: 521-532. 10.1046/j.1365-2958.2002.03032.x.View ArticlePubMed
  32. Shin S, Park C: Modulation of flagellar expression in Escherichia coli by acetyl phosphate and the osmoregulator OmpR. J Bacteriol. 1995, 177: 4696-4702.PubMed CentralPubMed
  33. Prüß BM, Besemann C, Denton A, Wolfe AJ: A complex transcription network controls the early stages of biofilm development by Escherichia coli. J Bacteriol. 2006, 188: 3731-3739. 10.1128/JB.01780-05.PubMed CentralView ArticlePubMed
  34. Holland IB, Schmitt L, Young J: Type 1 protein secretion in bacteria, the ABC-transporter dependent pathway (review). Mol Membr Biol. 2005, 22: 29-39.View ArticlePubMed
  35. Rodionov DA, Hebbeln P, Gelfand MS, Eitinger T: Comparative and functional genomic analysis of prokaryotic nickel and cobalt uptake transporters: evidence for a novel group of ATP-binding cassette transporters. J Bacteriol. 2006, 188: 317-327. 10.1128/JB.188.1.317-327.2006.PubMed CentralView ArticlePubMed
  36. Fruchterman TMJ, Reingold EM: Graph drawing by force-directed placement. Software - Practice and Experience. 1991, 21: 1129-1164. 10.1002/spe.4380211102.View Article
  37. JUNG: Java Universal Network/Graph Framework. 2006, [http://​jung.​sourceforge.​net]
  38. Kanehisa M, Goto S, Hattori M, oki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Research. 2006, 34: D354-D357. 10.1093/nar/gkj102.PubMed CentralView ArticlePubMed
  39. KEGG: Kyoto Encyclopedia of Genes and Genomes. 2006, [http://​www.​genome.​jp/​kegg]
  40. Lubitz SP, Weiner JH: The Escherichia coliynfEFGHI operon encodes polypeptides which are paralogues of dimethyl sulfoxide reductase (DmsABC). Arch Biochem Biophys. 2003, 418: 205-216. 10.1016/ ArticlePubMed
  41. Galperin MY: Bacterial signal transduction network in a genomic perspective. Environ Microbiol. 2004, 6: 552-567. 10.1111/j.1462-2920.2004.00633.x.PubMed CentralView ArticlePubMed
  42. Appleby JL, Parkinson JS, Bourret RB: Signal transduction via the multi-step phosphorelay: not necessarily a road less traveled. Cell. 1996, 86: 845-848. 10.1016/S0092-8674(00)80158-0.View ArticlePubMed
  43. Grunden AM, Self WT, Villain M, Blalock JE, Shanmugam KT: An analysis of the binding of repressor protein ModE to modABCD (molybdate transport) operator/promoter DNA of Escherichia coli. J Biol Chem. 1999, 274: 24308-24315. 10.1074/jbc.274.34.24308.View ArticlePubMed
  44. Tanapongpipat S, Reid E, Cole JA, Crooke H: Transcriptional control and essential roles of the Escherichia coli ccm gene products in formate-dependent nitrite reduction and cytochrome c synthesis. Biochem J. 1998, 334 ( Pt 2): 355-365.View Article
  45. Overton TW, Griffiths L, Patel MD, Hobman JL, Penn CW, Cole JA, Constantinidou C: Microarray analysis of gene regulation by oxygen, nitrate, nitrite, FNR, NarL and NarP during anaerobic growth of Escherichia coli: new insights into microbial physiology. Biochem Soc Trans. 2006, 34: 104-107. 10.1042/BST0340104.View ArticlePubMed
  46. BISON. 2006, [http://​denton.​cs.​ndsu.​nodak.​edu/​bison/​]
  47. RegulonDB. 2006, [http://​regulondb.​ccg.​unam.​mx/​index.​html]
  48. Oshima T, Aiba H, Masuda Y, Kanaya S, Sugiura M, Wanner BL, Mori H, Mizuno T: Transcriptome analysis of all two-component regulatory system mutants of Escherichia coli K-12. Mol Microbiol. 2002, 46: 281-291. 10.1046/j.1365-2958.2002.03170.x.View ArticlePubMed
  49. E. coli genome project. 2006, [http://​www.​genome.​wisc.​edu/​]
  50. Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, Gregor J, Davis NW, Kirkpatrick HA, Goeden MA, Rose DJ, Mau B, Shao Y: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1474. 10.1126/science.277.5331.1453.View ArticlePubMed
  51. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-D141. 10.1093/nar/gkh121.PubMed CentralView ArticlePubMed
  52. Pfam. 2006, [http://​www.​sanger.​ac.​uk/​Software/​Pfam/​]
  53. Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.View ArticlePubMed
  54. HMMER. 2004, [http://​hmmer.​janelia.​org/​]


© Besemann et al; licensee BioMed Central Ltd. 2006

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.