Boolean networks using the chi-square test for inferring large-scale gene regulatory networks
- Haseong Kim^{1},
- Jae K Lee^{2} and
- Taesung Park^{3}Email author
https://doi.org/10.1186/1471-2105-8-37
© Kim et al; licensee BioMed Central Ltd. 2007
Received: 06 February 2006
Accepted: 01 February 2007
Published: 01 February 2007
Abstract
Background
Boolean network (BN) modeling is a commonly used method for constructing gene regulatory networks from time series microarray data. However, its major drawback is that its computation time is very high or often impractical to construct large-scale gene networks. We propose a variable selection method that are not only reduces BN computation times significantly but also obtains optimal network constructions by using chi-square statistics for testing the independence in contingency tables.
Results
Both the computation time and accuracy of the network structures estimated by the proposed method are compared with those of the original BN methods on simulated and real yeast cell cycle microarray gene expression data sets. Our results reveal that the proposed chi-square testing (CST)-based BN method significantly improves the computation time, while its ability to identify all the true network mechanisms was effectively the same as that of full-search BN methods. The proposed BN algorithm is approximately 70.8 and 7.6 times faster than the original BN algorithm when the error sizes of the Best-Fit Extension problem are 0 and 1, respectively. Further, the false positive error rate of the proposed CST-based BN algorithm tends to be less than that of the original BN.
Conclusion
The CST-based BN method dramatically improves the computation time of the original BN algorithm. Therefore, it can efficiently infer large-scale gene regulatory network mechanisms.
Background
The advancement of high-throughput technologies, such as DNA chips, has enabled the study of interactions and regulations among genes on a genome-wide scale. Recently, many algorithms have been introduced to determine gene regulatory networks based on such high-throughput microarray data, including linear models [1, 2], Boolean networks [3–6], Bayesian networks [7, 8], neural networks [9], and differential equations [1, 10].
In the linear modeling of a genetic network, the expression data is fitted using a regression model, where the change in expression levels is a response for all other genes [1]. Although such standard linear modeling approaches enable the analysis of many different features of the modeled system, they are not effective in genome-wide network discovery. This is because the number of candidate parameters and models is very high and therefore it is difficult to search efficiently and reliably with tight control on many false positives.
Bayesian network algorithms have limitations with regard to the determination of an important network structure because of their complex modeling strategies (with a large number of parameters to be estimated) and a long computation time for searching all potential network structures on genome-wide expression data. These limitations of the Bayesian network may be overcome by the dynamic Bayesian network (DBN), which models the stochastic evolution of a set of random variables over time [11, 12]. Although some improvements have been proposed, the accuracy of prediction of the DBN is relatively low, and its excessive computational time is still very high [13].
Recently, studies on the hierarchical scale-free network in lower organisms [14, 15] have indicated the necessity of a network method for the simultaneous analysis of thousands of genes. The human B-cell network analysis using mutual information [16] is a type of hierarchical scale-free analysis. Although this analysis successfully constructs gene networks with thousands of genes, the method is based on mutual information between two genes; therefore, it cannot obtain the response of a target gene when more than two genes simultaneously affect the target gene.
Among these methods, the Boolean network (BN) is useful to construct gene regulatory networks observed by high-throughput microarray data because it can monitor the dynamic behavior in complex systems based on its binarization of such massive expression profiling data [17, 18]. The Boolean function of a gene in a BN can describe the behavior of the target gene according to simultaneous changes in the expression levels of the other genes.
Boolean networks
BN models were first introduced by Kauffman [3]. In these models, a gene expression is simplified with two levels: ON and OFF. A BN G(V, F) is defined by a set of node V = {x_{1}, ..., x_{ n }} and a set of Boolean functions F = {f_{1}, ..., f_{ n }}. A Boolean function f_{ i }(x_{1}, ..., x_{ k }), where i = {1, ..., n}, with k specified input nodes (indegree) is assigned to node x_{ i }. The regulation of nodes is defined by F. More specifically, given the values of the node (V) at time t - 1, the Boolean function are used to update these values at time t.
The model system has been developed into a so-called Random BN model [19]. BNs have attracted attention since the introduction of probabilistic Boolean network (PBN) models by Shmulevich et al. [6]. Many algorithms have been proposed for the inference of BNs. For example, the REVEAL algorithm has been introduced by Liang et al. [5] for causal inference by using mutual information, which is the most fundamental and general measure of correlation. Akutsu et al. [4] have constructed a BN structure based on the consistency problem, which can be used to determine the existence of a network that is consistent with the observed data. In one of the most recent studies on the BN algorithm, the Best-Fit Extension problem [20] is used for the inference of PBNs [6]. In PBNs, every node (gene) can have a chance of acquiring different Boolean functions. The probabilistic selection of Boolean functions adds flexibility in the determination of the steady state of BNs and monitoring of the dynamical network behavior for gene perturbation or intervention [18, 21].
Recently, several software packages have been developed for constructing BNs. The random BN toolbox [22] and the PBN toolbox [6] are available in Matlab. NetBuilder (version 0.94) [23] is a genetic regulatory network tool used to simulate genetic network using a BN. The BN has been widely used to describe biological processes. For example, Huang [17] has used these networks to represent cell growth, cell differentiation, and apoptosis. A transcriptional network model in yeast has been studied using a random BN [24]. Further, Johnson [25] has studied the signal transduction pathways in a B-cell ligand screen.
Advantages and disadvantages of Boolean networks
The estimation of gene regulatory networks using the BN offers several advantages. First, the BN model effectively explains the dynamic behavior of living systems [17, 18, 26]. Simplistic Boolean formalism can represent realistic complex biological phenomena such as cellular state dynamics that exhibit switch-like behavior, stability, and hysteresis [17]. It also enables the modeling of non-linear relations in complex living systems [27]. Second, Boolean algebra is an established science that provides a large set of algorithms that are already available for supervised learning in the binary domain, such as the logical analysis of data [28], and Boolean-based classification algorithms [29]. Finally, dichotomization to binary values improves the accuracy of classification and simplifies the obtained models by reducing the noise level in experimental data [30, 31].
However, the BN has some drawbacks. One of the major drawbacks is that it requires extremely high computing times to construct reliable network structures. Therefore, most BN algorithms such as REVEAL can thus be used only with a small number of genes and a low indegree value. For higher indegree values, these algorithms should be accelerated through parallelization in order to increase the search efficiency in the solution space [5]. The consistency problem [4] works in time complexity O(${2}^{{2}^{k}}$·$\left(\begin{array}{c}n\\ k\end{array}\right)$·m·n·poly(k)) (m is the number of observed time points; n, the total number of genes; and poly(k), the time required to compare a pair of examples respectively) for a fixed indegree k; this is because ${2}^{{2}^{k}}$ Boolean functions must be checked for each of the possible _{ n }C_{ k }combinations of variables and for m observations. The Best-Fit Extension problem [21] also works in time complexity O(${2}^{{2}^{k}}$·$\left(\begin{array}{c}n\\ k\end{array}\right)$·m·n·poly(k)). Although the improved consistency algorithm and Best-Fit Extension problem work in time complexity O($\left(\begin{array}{c}n\\ k\end{array}\right)$·m·n·poly(k)) [32], they still exhibit an exponential increase in the computing time for the parameters n and k. Such high computing times are a major problem in the study of large-scale gene regulatory and gene interaction systems using BNs.
Chi-square-test-based Boolean network
In order to overcome the time complexity problem of the BN method, we propose a variable selection method based on the chi-square test (CST). The proposed CST-based BN adopts the Best-Fit Extension problem, which is commonly used in the PBN to effectively determine all possible relevant Boolean functions. In our method, the maximum indegree of networks is assumed to be three. We also focus on the Boolean functions that comprise three different literals (input genes in the Boolean function). Each literal is connected by the three Boolean operators NOT(¬), AND(∧), OR(∨); for example, f = X 1 ∧ ¬X 2 ∨ X 3. Then, the time complexity of the CST-based BN reduces to O(${2}^{{2}^{k}}\cdot {\displaystyle {\sum}_{j=1}^{n}{\displaystyle {\sum}_{i=1}^{{n}_{1,j}}\cdot \left(\begin{array}{c}{n}_{2,ij}\\ k\end{array}\right)}}$·m·poly(k)), where n is the total number of genes, k is the indegree, m is the total number of time points, n_{1, j}is the number of first selected genes for the j th gene, and n_{2, ij}is the number of second selected genes when the i th gene is selected in the first step. We have found that the dichotomization of the continuous gene expression values allows us to efficiently perform the independence test for a two-way contingency table on each pairwise network mechanisms. We use the CST to identify genes that are associated with a target gene. A target gene would be expressed in accordance with a Boolean function related to the selected genes. Since the genes have only two levels, (0 and 1), we use 2 × 2 and 2 × 2 × 2 contingency tables to identify the relationship between two and three genes respectively. The proposed method is used along with the Best-Fit Extension problem. This method is described in detail in the Methods section.
Results
Simulation study
The simulation results of the original BN and CST-based BN.
Noise level | Original BN Number of Boolean functions | FPR | CST-based BN Number of Boolean functions | FPR |
---|---|---|---|---|
0 | 348(27)* | 0.922 | 286(27) | 0.899 |
0.01% | 326(22) | 0.932 | 264(22) | 0.916 |
0.02% | 269(22) | 0.918 | 207(22) | 0.893 |
0.04% | 313(20) | 0.936 | 289(20) | 0.930 |
0.06% | 264(18) | 0.931 | 244(18) | 0.926 |
0.08% | 168(12) | 0.928 | 139(12) | 0.913 |
0.1% | 271(17) | 0.937 | 216(15) | 0.930 |
0.12% | 243(16) | 0.934 | 209(16) | 0.923 |
0.14% | 614(19) | 0.969 | 514(19) | 0.963 |
0.16% | 73(9) | 0.876 | 73(9) | 0.876 |
0.18% | 269(15) | 0.944 | 217(15) | 0.930 |
0.2% | 198(9) | 0.954 | 167(9) | 0.946 |
0.22% | 54(5) | 0.907 | 54(5) | 0.907 |
0.24% | 369(3) | 0.991 | 318(3) | 0.990 |
The selected number of nodes in the first and second steps of variable selection
Node number | Number of selected nodes at the first step | Number of selected nodes at the second step |
---|---|---|
1 | 1 | 39 |
2 | 2 | 39, 38 |
3 | 1 | 39 |
4 | 1 | 39 |
5 | 1 | 39 |
6 | 0 | 0 |
7 | 8 | 39, 37, 36, 36, 32, 32, 33, 32 |
8 | 1 | 39 |
9 | 4 | 39, 37, 37, 35 |
10 | 3 | 39, 38, 37 |
11 | 0 | 0 |
12 | 9 | 39, 38, 37, 36, 32, 32, 33, 32, 31 |
13 | 0 | 0 |
14 | 0 | 0 |
15 | 0 | 0 |
16 | 4 | 39, 37, 37, 36 |
17 | 0 | 0 |
18 | 3 | 39, 38, 37 |
19 | 1 | 39 |
20 | 0 | 0 |
21 | 5 | 39, 37, 37, 36, 35 |
22 | 5 | 39, 37, 37, 36, 35 |
23 | 2 | 39, 38 |
24 | 5 | 39, 37, 37, 36, 35 |
25 | 0 | 0 |
26 | 0 | 0 |
27 | 5 | 39, 38, 37, 33, 33 |
28 | 0 | 0 |
29 | 1 | 39 |
30 | 7 | 39, 34, 35, 35, 35, 34, 33 |
31 | 10 | 39, 34, 33, 32, 31, 32, 33, 32, 31, 30 |
32 | 8 | 39, 34, 33, 32, 33, 33, 33, 32 |
33 | 1 | 39 |
34 | 0 | 0 |
35 | 1 | 39 |
36 | 2 | 39, 38 |
37 | 0 | 0 |
38 | 1 | 39 |
39 | 0 | 0 |
40 | 0 | 0 |
In summary, the CST-based BN method was approximately 6.9 times faster than the original BN method. If the network had a larger number of nodes (n), then the difference between the computing times of the two algorithms would be significantly high.
Yeast cell cycle data
In order to demonstrate the improvement in the computing times, we apply the proposed variable selection method to yeast cell cycle data [33]. The data comprises 18 time points (alpha-factor-based synchronization experiment). In this example, the computing time and accuracy of network estimation of the original BN method are compared with those of our CST-based BN method.
Comparison between the network structure estimation accuracies of the CST-based BN and the original BN
Our variable selection method significantly improves the computing time of the BNs. However, the accuracy of our method should be assessed before comparing the computing times of the two methods. The improvement in the computing times primarily depends on the cutoff statistical significance levels, α_{1} and α_{2}, for the selection of genes at time t - 1 via the CST (Methods section). Depending on the choice of appropriate values of α_{1} and α_{2}, the proposed CST-based BN method may not be able to determine some Boolean functions that can be estimated by the original BN method. This may be attributed to the missing essential variables due to the usage of extremely stringent cutoff values.
We define the error rate as the discrepancy between the Boolean functions estimated by the original BN and the proposed CST-based BN as follows:
In order to select appropriate values of α_{1} and α_{2}, the error rates were obtained for various combination of (α_{1}, α_{2}). Some of the results are shown in Figures 3 and 4. To obtain an error rate of zero, the values of α_{1} and α_{2} should be greater than 1% and 2%, respectively, when the error size is zero (Figure 3(a),(b)). However, when the error size is one, larger values of α_{1} and α_{2} are required to obtain an error rate of zero (Figure 4(a)–(d)). Based on these results, we suggest that the following values be used when the error size is 0, α_{1} = 1%, α_{2} = 2% and when the error size is 1, α_{1} = 7.5%, α_{2} = 10%.
Comparison between the computing times of the original and the proposed BN algorithms
Construction of gene networks with yeast cell cycle related 800 genes
For all 800 genes, the CST-based BN required approximately three days to construct the network structure. In order to estimate the total computation time of the original BN, we selected the first target gene and constructed the BN, which required approximately 37,011 s. Therefore, the total computation times for all 800 genes will be approximately 342 days to build the network structure. Hence, the computing time of the proposed CST-based BN method is approximately 114 times faster than that of the original BN method.
Discussions and Conclusion
Recent studies [14, 15] have emphasized that thousands of genes must be considered simultaneously in order to construct gene regulatory networks in an organism. The BN method is useful for constructing a gene regulatory network. If the gene expression data contain a considerable amount of noise, the binary transformation of these data can reduce the error [35]. The BN has been successfully used to model a nonlinear system [27] and the dynamic behavior of living systems [18, 26]. Despite these advantages, it is difficult to apply the BN method to large-scale gene regulatory network studies due to the extremely high computation times.
In order to overcome this computational drawback, we proposed the variable selection method using the CST for the two-way and three-way contingency tables of Boolean count observations. This method reduces the computation times significantly; for example, for 120 genes, the computation time is approximately 70.8 times faster than that of the original BN method. If the total number of genes and the value of k increase, the improvement in the computation time is expected to be significantly greater than the original BN method. Also the proposed method can be easily implemented with the existing BN modeling algorithms such as the PBNs by efficiently selecting only the most relevant genes for determining the Boolean functions. This method is thus demonstrated to reduce the false positive rate, which is an important problem in network studies conducted on a genome-wide scale.
In our method, the value of k is assumed to be three. However, it is possible to use a large value of k greater than 3. Since our method uses the Best-Fit Extension problem, a gene can be controlled by more than one Boolean function. Therefore, it appears that k = 3 provides a large number of Boolean functions that can model the gene regulatory network successfully.
The proposed CST-based BN used the two-step discovery of 3-indegree Boolean functions because the prediction information of addtional genes in such high-dimensional Boolean functions is mainly observed after considering, or conditional on, primary genes' effects. This strategy, in turn, resulted in much more efficient discovery of the most predictive high-dimensional Boolean functions in our BN modeling.
The result of the proposed CST-based BN may be sensitive to the sample size n. When n is small, the contingency table may contain many cells with low and zero frequencies. To ensure that the expectation value is not equal to zero, a continuity correction is used by adding a small constant 0.1 to the observed frequency in each cell [36]. This simple correction produces a successful result in the real data set [33] that contains 17 time points. However, we should be more careful while applying the CST method to data with small time points because the result of the CST can be less reliable for the sparse data set. In this case, we suggest that the CST be replaced with Fisher's exact test that provides a more reliable result for the small sample size data [36]. The small sample size problem also makes it difficult for the original BN algorithm to produce a reliable result. We think that in the near future the advancement of high throughput techniques and the cost-down of microarrays will enable us to solve the sparse data problem by producing large data sets easily.
The improvement of the computing times using the CST-based BN will significantly increase the utility and applicability of the BN to the inference of various regulatory networks, particularly those based on current large screening biological data such as microarrays. In order to apply the proposed variable selection method, we must first select the values of α_{1} and α_{2}. A small values of α cause the exclusion of essential Boolean functions and produces a high error rate. On the other hand, large values of γ cause the inclusion of most of the genes, thereby resulting in long computing times. For the yeast cell cycle data, we applied various combinations of (α_{1}, α_{2}). As shown in Figure 3, it appears that the cutoff values do not significantly affect much the accuracy of the method for the yeast cell cycle data, provided the values of α_{1} and α_{2} are greater than 1% and 2%, respectively, when the error size of the Best-Fit Extension problem is zero.
Therefore, for practical application, we suggest that α_{1} = 1% and α_{2} = 2% be used when the error size of the Best-Fit Extension problem is zero. We select the cutoff values such that the Boolean functions obtained by the CST-based BN are the same as those obtained by the original BN. However, we can use smaller cutoff values to reduce the number of false positive Boolean functions, because the original BN method tends to yield many false positive relationships, as shown in the simulation result.
In addition, a more careful dichotomization is required for a more accurate biological interpretation of the network structure. For example, since microarray data have continuous expression values with a considerably large amount of information, the dichotomization may require the selection of an appropriate threshold value depending on the biological function of each gene [35]. The performance may not be remarkable when a small number of time points and genes are available. However, we show that the proposed variable selection method is significantly more efficient for large-scale gene regulatory network studies. For example, the CST-based BN is approximately 114 times faster than the original BN for 800 genes (Figure 6).
The next step would be to perform a biological evaluation of the selected network structure. However, the main focus of our study is to improve the computation times of the BN by using the CST. Our approach allows the application of the BN to genome-wide network construction and discovery. A future study will evaluate the accuracy of the BN and compare it with other network methods such as the Bayesian network and hierarchical scale-free network.
Methods
The proposed CST-based BN consists of two steps. The first step is to determine a pair of genes that are associated with each other. The second step is to determine the third gene that is conditionally associated with the pair of genes identified in the first step.
First step for the main effect
Let n be the total number of genes. In the first step, 2 × 2 contingency tables are constructed from the dichotomized gene expression data. The p th row comprises the i th gene expression level at time t - 1 while the q th column of the table comprises the j th gene expression level at time t (i = 1, ..., n; j = 1, ..., n; p = 0,1; q = 0,1). For the i th and j th genes, a 2 × 2 contingency table is constructed with four cells: {0, 0}, {0,1}, {1, 0}, and {1,1}, where {p, q} represent the i th gene expression level at time t - 1 and the j th gene expression level at time t, respectively.
A CST statistic is then computed for testing the independence between two genes. For multinomial sampling with probabilities {π_{ pq }} in the contingency table, the null hypothesis of independence is H_{0} : π_{ pq }= π_{p+}π_{q+}(the i th gene at time t - 1 and the j th gene at time t are independent) for all p (= 0,1) and q (= 0,1). The conventional Pearson's CST can be used to test H_{0} using the observed frequency O_{ pq }and the expected frequency E_{ pq }under H_{0}. For the continuity correction, we add an arbitrary small number a to each observed frequency in order to prevent E_{ pq }from becoming zero [36]. We use a = 0.1 for the correction. Generally, {π_{p+}} and {π_{+q}} are unknown. The maximum likelihood (ML) estimates are the sample marginal proportions {$\widehat{\pi}$_{p+}= O_{p+}/O_{++}} and {$\widehat{\pi}$_{+q}= O_{+q}/O_{++}}, where O_{++} = Σ_{ p }Σ_{ q }O_{ pq }. E_{ pq }is estimated as E_{ pq }= O_{++}$\widehat{\pi}$_{p+}$\widehat{\pi}$_{+q}= O_{p+}O_{+q}/O_{++}. Therefore, the chi-square statistic is expressed as follows:
Using this CST, the significant genes are selected by an appropriate selection criterion α_{1}. A further discussion on the appropriate choice of α_{1} is provided in the Result section.
Second step for the conditional effect
Assume that the i th gene at time t - 1 is selected in the first step for the j th gene at time t. Then, a 2 × 2 × 2 contingency table can be constructed that consists of three genes – the i th and j th genes selected in the first step and an additional new gene h at time t - 1. This contingency table consists of eight cells: {0,0,0}, {0,0,1}, {0,1,0}, {0,1,1}, {1,0,0}, {1,0,1}, {1,1,0}, and {1,1,1}, where {o, p, q} represent the h th gene expression level at time t - 1, i th gene expression level at time t - 1, and j th gene expression level at time t, respectively (i = 1, ..., n; j = 1, ..., n; h = 1, ..., n; o = 0,1; p = 0,1; q = 0,1).
For the given expression value of h, there are two 2 × 2 contingency tables for the i and j genes. We focus on the conditional independence test. The null hypothesis that the i th gene at time t - 1 and the j th gene at time t are conditionally independent when the h th gene expression level at time t - 1 is given by H_{0} : π_{pq|o}= π_{p+|o}π_{q+|o}for all p (= 0, 1) and q (= 0,1), where π_{..|o}represents the conditional probability for the given o. We use the CST to test H_{0} using the observed frequency O_{ opq }and the expected frequency E_{ opq }under H_{0}. We also add 0.1 to each observed frequency for the continuity correction. The ML estimates of π_{p+|o}and π_{+q|o}are the sample conditional proportions {$\widehat{\pi}$_{p+|o}= O_{op+}/O_{o++}} and, {$\widehat{\pi}$_{+q|o}= O_{o+q}/O_{o++}}, respectively where O_{o++}= Σ_{ p }Σ_{ q }O_{ opq }. E_{pq|o}is estimated as E_{pq|o}= O_{o++}$\widehat{\pi}$_{p+|o}$\widehat{\pi}$_{+q|o}= O_{op+}O_{o+q}/O_{o++}. Then, the chi-square statistic are given by
for o = 0,1. We assume that the i th and j th genes are not independent from the first step. We select the h th gene if at least one of the two CST in the second step is significant (the p-value of the test is less than α_{2}).
The rationale for using this conditional independence test is that h affects the association between two genes i and j. The conditional test approach is very effective in the determination of a relationship for more than two genes.
Implementation of two step variable selection method
We select n_{1, j}genes at time t - 1 that are associated with the j th gene at time t in the first step where 1 ≤ n_{1, j}≤ n. If the i th gene is one of the n_{1, j}genes, we select n_{2, ij}genes in the second step for 1 ≤ n_{2, ij}≤ n - 1 (excluding the i th gene). Then, we consider all possible combinations for selected three genes (when k = 3); one gene is the i th gene selected in the first step and the other two genes are selected in the second step. These combinations are used instead of all possible combinations in the original BN algorithms. The time complexity of determining the Boolean functions for the j th gene is O(${2}^{{2}^{k}}\cdot {\displaystyle {\sum}_{i=1}^{{n}_{1,j}}\cdot \left(\begin{array}{c}{n}_{2,ij}\\ k\end{array}\right)}$·m·poly(k)) and that of variable selection using the CST is O(n^{2} + n_{1, j}·(n - 1)). Hence, the total time complexity of the proposed algorithm is expressed as follows:
As n increases, the time complexity of determining the Boolean functions dominates the time complexity of variable selection because the former increases more rapidly than the latter.
EXAMPLE
Example of dichotomized gene expression profiles
time | G1 | G2 | G3 | G4 | G5 | G6 | G7 | G8 |
---|---|---|---|---|---|---|---|---|
t1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 |
t2 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 |
t3 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 |
t4 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 |
t5 | 1 | 1 | 1 | 0 | 0 | 0 | 1 | 1 |
t6 | 0 | 1 | 0 | 0 | 1 | 1 | 1 | 1 |
t7 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 1 |
t8 | 0 | 1 | 0 | 0 | 1 | 1 | 1 | 1 |
t9 | 1 | 1 | 1 | 0 | 1 | 0 | 1 | 1 |
t10 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 |
t11 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 |
t12 | 0 | 0 | 0 | 1 | 1 | 1 | 0 | 0 |
t13 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 |
t14 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
t15 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 |
t16 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 |
t17 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 1 |
t18 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 |
Result of variable selection with 8 genes
First step | Second step | |||||
---|---|---|---|---|---|---|
j th gene(t) | i th gene(t - 1) | p-value | h th gene(t - 1) | p-value1 | p-value2 | combinations |
G1 | G1 | 0.0013 | G5 | 0.0185 | 0.0252 | _{4}C_{2} = 6 |
G6 | 0.0191 | 0.033 | ||||
G7 | 0.0059 | 0.0656 | ||||
G8 | 0.1049 | 0.0192 | ||||
G3 | 0.0253 | G2 | 0.0059 | 0.0208 | _{6}C_{2} = 15 | |
G4 | 0.0208 | 0.0059 | ||||
G5 | 0.0185 | 0.0071 | ||||
G6 | 0.0191 | 0.0058 | ||||
G7 | 0.0414 | 0.0035 | ||||
G8 | 0.1049 | 0.0021 | ||||
G3 | G1 | 0.0013 | G2 | 0.0065 | 0.0835 | _{6}C_{2} = 15 |
G4 | 0.0035 | 0.1165 | ||||
G5 | 0.0185 | 0.0252 | ||||
G6 | 0.0033 | 0.1016 | ||||
G7 | 0.033 | 0.0191 | ||||
G8 | 0.1049 | 0.0192 | ||||
G3 | 0.0078 | G5 | 0.0185 | 0.1443 | 1 | |
G7 | G4 | 0.0078 | G1 | 0.0102 | 0.2063 | _{5}C_{2} = 10 |
G2 | 0.4774 | 0.0033 | ||||
G3 | 0.0059 | 0.3535 | ||||
G5 | 0.5469 | 0.0071 | ||||
G8 | 0.0143 | 0.4864 | ||||
Total number of combinations | 47 |
Availability and requirements
Project name: Boolean networks for Large-scale gene regulatory network Project home page: http://bibs.snu.ac.kr/supplement/2006/Boolean/
Declarations
Acknowledgements
The authors would like to thank to anonymous referees whose comments were extremely helpful. This study was supported by the National Research Laboratory Program of the Korea Science and Engineering Foundation (M10500000126), the Brain Korea 21 Project of the Ministry of Education of T.P., and the US NIH grant 1R01HL081690 of J.K.L.
Authors’ Affiliations
References
- D'Haeseleer P, Wen X, Fuhrman S, Somogyi R: Linear modeling of mRNA expression levels during CNS development and injury. Pacific Symposium on Biocomputing 1999, 4: 41–52.Google Scholar
- van Someren EP, Wessels LFA, Reinders MJT: Linear Modeling of Genetic Networks from Experimental Data. Intelligent Systems for Molecular Biology 2000, 8: 355–366.Google Scholar
- Kauffman SA: Metabolic stability and epigenesis in randomly constructed genetic nets. Journal of Theoretical Biology 1969, 9: 3273–3297.Google Scholar
- Akutsu T, Miyano S, Kuhara S: Identification of genetic networks from a small number of gene expression patterns under the Boolean network model. Pacific Symposium on Biocomputing 1999, 4: 17–28.Google Scholar
- Liang S, Fuhrman S, Somogyi R: REVEAL, A general reverse engineering algorithm for inference of genetic network architectures. Pacific Symposium on Biocomputing 1998, 3: 18–29.Google Scholar
- Shmulevich I, Dougherty ER, Seungchan K, Zhang W: Probabilistic Boolean networks: A rule-based uncertainty model for gene regulatory networks. Bioinformatics 2002, 18: 261–274. 10.1093/bioinformatics/18.2.261View ArticlePubMedGoogle Scholar
- Friedman N, Goldszmidt M, Wyner A, Eds: Data analysis with baysian networks: A bootstrap approach. Proc Fifteenth Conf on Uncertainty in Artificial Intelligence (UAI). 1999.Google Scholar
- Imoto S, Goto T, Miyano S: Estimation of Genetic Networks and Functional Structures Between Genes by Using Bayesian Networks and Nonparametric Regression. Pacific Symposium on Biocomputing 2002, 7: 175–186.Google Scholar
- Weaver DC, Workman CT, Stormo GD: Modeling regulatory networks with weight matrices. Pacific Symposium on Biocomputing 1999, 4: 112–123.Google Scholar
- Chen T, He HL, Church GM: Modeling gene expression with differential equations. Pacific Symposium Biocomputing 1999, 4: 29–40.Google Scholar
- Perrin B, Ralaivola L, Mazurie A, Bottani S, Mallet J, D'Alche-Buc F: Gene networks inference using dynamic Bayesian networks. Biofincormatics 2003, 19: 138–148.Google Scholar
- Dojer N, Gambin A, Mizera A, Wilczynski B, Tiuryn J: Applying dynamic Bayesian networks to perturbed gene expression data. BMC Bioinformatics 2006, 7: 249. 10.1186/1471-2105-7-249PubMed CentralView ArticlePubMedGoogle Scholar
- Zou M, Conzen S: A new dynamic Bayesian network(DBN) approach for identifying gene regulatory networks from time course microarray data. Bioinformatics 2005, 21: 71–79. 10.1093/bioinformatics/bth463View ArticlePubMedGoogle Scholar
- Han JDJ, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, Dupuy D, Walhout AJM, Cusic ME, Roth FP, Vidal M: Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature 2004, 430: 88–93. 10.1038/nature02555View ArticlePubMedGoogle Scholar
- Jeong H, Tombor B, Albert R, Oltvai ZN, Barabasi AL: The large-scale organization of metabolic networks. Nature 2000, 407: 651–654. 10.1038/35036627View ArticlePubMedGoogle Scholar
- Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, Califano A: Reverse engineering of regulatory networks in human B cells. Nature Genetics 2005, 37: 382–390. 10.1038/ng1532View ArticlePubMedGoogle Scholar
- Huang S: Gene expression profiling, genetic networks and cellular states: An integrating concept for tumorigenesis and drug discovery. Journal of Molecular Medicine 1999, 77: 469–480. 10.1007/s001099900023View ArticlePubMedGoogle Scholar
- Shmulevich I, Gluhovsky I, Hashimoto RF, Dougherty ER, Zhang W: Steady-state analysis of genetic regulatory networks modelled by probabilistic Boolean networks. Comparative and Functional Genomics 2003, 4: 601–608. 10.1002/cfg.342PubMed CentralView ArticlePubMedGoogle Scholar
- Kauffman SA: The Origins of Order: Self-organization and Selection in Evolution. New York. Oxford University Press; 1993.Google Scholar
- Boros E, Ibaraki T, Makino K: Error-Free and Best-Fit Extensions of partially defined Boolean functions. Information and Computation 1998, 140: 254–283. 10.1006/inco.1997.2687View ArticleGoogle Scholar
- Shmulevich I, Saarinen A, Yli-Harja O, Astola J, Eds: Inference of genetic regulatory networks under the best-fit extension paradigm, in Computational and Statistical Approaches To Genomics. Boston, MA: Kluwer; 2002.Google Scholar
- Schwarzer C: Matlab Random Boolean Network Toolbox 2003.[http://www.teuscher.ch/rbntoolbox/]
- Schilstra MJ, Bolouri H: Modeling the regulation of gene expression in genetic regulatory networks.[http://strc.herts.ac.uk/bio/maria/NetBuilder]
- Kauffman SA, Peterson C, Samuelsson B, Troein C, Eds: Random Boolean network models and the yeast transcriptional network. Journal of Molecular Medicine 1999, USA. 2003, 77: 469–480.Google Scholar
- Johnson S, Ed: Boolean network inference and experiment design for the B-Cell single ligand screen. 2004. AfCS annual meeting; 2004.Google Scholar
- Shmulevich I, Dougherty ER, Zhang W: Gene perturbation and intervention in probabilistic Boolean networks. Bioinformatics 2002, 18: 1319–1331. 10.1093/bioinformatics/18.10.1319View ArticlePubMedGoogle Scholar
- Thomas R: Regulatory networks seen as asynchronous automata: a logical description. Journal of Theoretical Biology 1991, 153: 1–23. 10.1016/S0022-5193(05)80350-9View ArticleGoogle Scholar
- Boros E, Hammer PL, Ibaraki T, Kogan A: Logical analysis of numerical data. Math Program 1997, 79: 163–190. 10.1016/S0025-5610(97)00050-6Google Scholar
- Akutsu T, Miyano S, Eds: Selecting informative genes for cancer classification using gene expression data. In Proceddings of the IEEE-EURASIP Workshop on NonlinSignal and Image Processing (NSIP). Baltimore, MD; 2001.Google Scholar
- Pfahringer B, Ed: Compression-based discretization of continuous attributes. Machine Learning: Procees of the Twelfth International Conference. Edited by: Prieditis A, Russell S. San Francisco: Morgan Kaufmann; 1995.Google Scholar
- Dougherty J, Kohavi R, Sahami M, Eds: Supervised and unsupervised discretization of continuous features. In Proceedings of the Twelfth International Conference on Machine Learning. Tahoe City, CA: Morgan Kaufmann; 1995.Google Scholar
- Lahdesmaki H, Shmulevich I, Yli-Harja O: On learning gene regulatory networks under the Boolean network model. Machine Learning 2003, 52: 147–167. 10.1023/A:1023905711304View ArticleGoogle Scholar
- Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell 1998, 9: 3273–3297.PubMed CentralView ArticlePubMedGoogle Scholar
- Dougherty ER, Kim S, Chen Y: Coefficient of determination in nonlinear signal processing. Signal Process 2000, 80: 2219–2235. 10.1016/S0165-1684(00)00079-7View ArticleGoogle Scholar
- Shmulevich I, Zang W: Binary analysis and optimization-based normalization of gene expression data. Bioinformatics 2002, 18: 555–565. 10.1093/bioinformatics/18.4.555View ArticlePubMedGoogle Scholar
- Agresti A: Categorical data analysis. second edition. wiley-interscience; 2002.View ArticleGoogle Scholar
Copyright
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.