Estimated time series of mRNA expressions.

<p>These time series are examples at <i>α</i> = 0.25 and <i>C</i> = 0.1, and we use the first optimal solution. Compared with the experimental data, each expression is classified into True Positive (TP), True Negative (TN), False Positive (FP), or False Negative (FN). The other parameters are fixed at <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], <i>ϵ</i> = 0.5, and <i>θ</i> = 1.</p

Estimating optimal sparseness of developmental gene networks using a semi-quantitative model

<div><p>To estimate gene regulatory networks, it is important that we know the number of connections, or sparseness of the networks. It can be expected that the robustness to perturbations is one of the factors determining the sparseness. We reconstruct a semi-quantitative model of gene networks from gene expression data in embryonic development and detect the optimal sparseness against perturbations. The dense networks are robust to connection-removal perturbation, whereas the sparse networks are robust to misexpression perturbation. We show that there is an optimal sparseness that serves as a trade-off between these perturbations, in agreement with the optimal result of validation for testing data. These results suggest that the robustness to the two types of perturbations determines the sparseness of gene networks.</p></div

Number of estimated connections as a function of sparseness parameter <i>α</i>.

<p>We show the average number for each soft-margin parameter <i>C</i> in 1,000 optimal solutions. The other parameters are fixed at <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], and <i>θ</i> = 1.</p

Model structure.

<p>Boxes show variables: (A) binary protein, (B) mRNA, (C) binary mRNA, and (D) protein. Arrows imply biological processes: (A)→(B) transcription, (B)→(C) mRNA export, (C)→(D) translation, and (D)→(A) protein import.</p

Robustness of perturbed networks as a function of sparseness parameter <i>α</i> at various noise intensities (A) <i>λ</i> = 10, (B) <i>λ</i> = 20, and (C) <i>λ</i> = 30.

<p>Both connection-removal and misexpression perturbations are applied to the estimated networks. As the robustness, we show the average F-measure for each number of connection removals <i>n</i>_<i>c</i> in 1,000 optimal solutions. The filled mark indicates the maximum value in each profile. The other parameters are fixed at <i>μ</i> = 0.01, <i>C</i> = 0.1, <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], <i>ϵ</i> = 0.5, and <i>θ</i> = 1.</p

Time course of cell differentiation in the model.

<p>A signal can be transferred between only neighbouring domains.</p

Accuracy for gene perturbation data as a function of sparseness parameter <i>α</i>.

<p>We show the average F-measure for each soft-margin parameter <i>C</i> in 1,000 optimal solutions. The optimal F-measure is observed at (<i>α</i>, <i>C</i>) = (0.25, 0.1) (indicated by an arrow). The other parameters are fixed at <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], <i>ϵ</i> = 0.5, and <i>θ</i> = 1.</p

Accuracy for time-series data as a function of sparseness parameter <i>α</i>.

<p>We show the average F-measure for each soft-margin parameter <i>C</i> in 1,000 optimal solutions. F-measure is calculated in the experimental data after 6 h. The other parameters are fixed at <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], <i>ϵ</i> = 0.5, and <i>θ</i> = 1.</p

Robustness of perturbed networks as a function of sparseness parameter <i>α</i>.

<p>As the robustness, we show the average F-measure in 1,000 optimal solutions. The types of perturbations are (A) connection-removal perturbation (<i>n</i>_<i>c</i> = 10), (B) misexpression perturbation (<i>λ</i> = 10, <i>μ</i> = 0.01), and (C) the addition of both perturbations. Optimal robustness in (C) is observed at <i>α</i> = 0.25 (indicated by an arrow). The other parameters are fixed at <i>C</i> = 0.1, <i>d</i> = 2 [h], <i>d</i>_<i>min</i> = 0 [h], <i>d</i>_<i>max</i> = 2 [h], <i>ϵ</i> = 0.5, and <i>θ</i> = 1.</p

MOESM3 of The role of retinoic acid signaling in starfish metamorphosis

Additional file 3: Table S1. Accession numbers of the gene used for construction of phylogenic tree. Amino acid sequences to construct tree were obtained from Uniprot or Echinobase ( http://www.echinobase.org/Echinobase/ ), Genbank