Algorithm of myogenic differentiation in higher-order organisms

A Boolean Function for Neural Induction Reveals a Critical Role of Direct Intercellular Interactions in Patterning the Ectoderm of the Ascidian Embryo

A complex system of multiple signaling molecules often produce differential gene expression patterns in animal embryos. In the ascidian embryo, four signaling ligands, Ephrin-A.d (Efna.d), Fgf9/16/20, Admp, and Gdf1/3-r, coordinately induce Otx expression in the neural lineage at the 32-cell stage. However, it has not been determined whether differential inputs of all of these signaling pathways are really necessary. It is possible that differential activation of one of these signaling pathways is sufficient and the remaining signaling pathways are activated in all cells at similar levels. To address this question, we developed a parameter-free method for determining a Boolean function for Otx expression in the present study. We treated activities of signaling pathways as Boolean values, and we also took all possible patterns of signaling gradients into consideration. We successfully determined a Boolean function that explains Otx expression in the animal hemisphere of wild-type and morphant embryos at the 32-cell stage. This Boolean function was not inconsistent with three sensing patterns, which represented whether or not individual cells received sufficient amounts of the signaling molecules. These sensing patterns all indicated that differential expression of Otx in the neural lineage is primarily determined by Efna.d, but not by differential inputs of Fgf9/16/20, Admp, and Gdf1/3-r signaling. To confirm this hypothesis experimentally, we simultaneously knocked-down Admp, Gdf1/3-r, and Fgf9/16/20, and treated this triple morphant with recombinant bFGF and BMP4 proteins, which mimic Fgf9/16/20 and Admp/Gdf1/3-r activity, respectively. Although no differential inputs of Admp, Gdf1/3-r and Fgf9/16/20 signaling were expected under this experimental condition, Otx was expressed specifically in the neural lineage. Thus, direct cell–cell interactions through Efna.d play a critical role in patterning the ectoderm of the early ascidian embryo.

It is often difficult to understand a complex system of multiple signaling molecules in animal embryos only with experimental procedures. Although theoretical analysis might solve this problem, it is often difficult to precisely determine parameters for signaling gradients and kinetics of signaling molecules. In the present study, we developed a parameter-free method for determining a Boolean function for understanding a complex signaling system using gene expression patterns of signaling molecules and geometrical configurations of individual cells within the embryo. In the ascidian embryo, four signaling ligands, Ephrin-A.d (Efna.d), Fgf9/16/20, Admp, and Gdf1/3-r, coordinately induce Otx expression in the neural lineage at the 32-cell stage. In addition to determining a Boolean function, our method determined sensing patterns, which represented whether or not individual cells received sufficient amounts of the signaling molecules. The sensing patterns predicted that differential expression of Otx in the neural lineage is primarily determined by Efna.d, but not by differential inputs of Fgf9/16/20, Admp, and Gdf1/3-r. We confirmed this prediction by an experiment. As a result, we found that only Efna.d signaling pathway is differentially activated between ectodermal cells and the remaining signaling pathways are activated in all ectodermal cells at similar levels.

SPABBATS: A pathway-discovery method based on Boolean satisfiability that facilitates the characterization of suppressor mutants

Background

Several computational methods exist to suggest rational genetic interventions that improve the productivity of industrial strains. Nonetheless, these methods are less effective to predict possible genetic responses of the strain after the intervention. This problem requires a better understanding of potential alternative metabolic and regulatory pathways able to counteract the targeted intervention.

Results

Here we present SPABBATS, an algorithm based on Boolean satisfiability (SAT) that computes alternative metabolic pathways between input and output species in a reconstructed network. The pathways can be constructed iteratively in order of increasing complexity. SPABBATS allows the accumulation of intermediates in the pathways, which permits discovering pathways missed by most traditional pathway analysis methods. In addition, we provide a proof of concept experiment for the validity of the algorithm. We deleted the genes for the glutamate dehydrogenases of the Gram-positive bacterium Bacillus subtilis and isolated suppressor mutant strains able to grow on glutamate as single carbon source. Our SAT approach proposed candidate alternative pathways which were decisive to pinpoint the exact mutation of the suppressor strain.

Conclusions

SPABBATS is the first application of SAT techniques to metabolic problems. It is particularly useful for the characterization of metabolic suppressor mutants and can be used in a synthetic biology setting to design new pathways with specific input-output requirements.