Cell signaling networks regulate a variety of developmental and physiological processes, and changes in their response to external stimuli are often implicated in disease initiation and progression. To elucidate how different responses can arise from conserved signaling networks, we have developed a mathematical model of the well-characterized Caenorhabditis vulval development network involving EGF, Wnt and Notch signaling that recapitulates biologically observed behaviors. We experimentally block a specific element of the EGF pathway (MEK), and find different behaviors in vulval development in two Caenorhabditis species, C. elegans and C. briggsae. When we separate our parameters into subsets that correspond to these two responses, they yield model behaviors that are consistent with observed experimental results, despite the initial parameter grouping based on perturbation in a single node of the EGF pathway. Finally, our analysis predicts specific parameters that may be critical for the theoretically and experimentally observed differences, suggesting modifications that might allow intentional switching between the two species' responses. Our results indicate that all manipulations within a signal transduction pathway do not yield the same outcome, and provide a framework to identify the specific genetic perturbations within a conserved network that will confer unique behaviors on the network.

Studies of vulval development in the nematode C. elegans have identified many genes that are involved in cell division and differentiation processes. Some of these encode components of conserved signal transduction pathways mediated by EGF, Notch, and Wnt. To understand how developmental mechanisms change during evolution, we are doing a comparative analysis of vulva formation in C. briggsae, a species that is closely related to C. elegans. Here, we report 14 mutations in 7 Multivulva (Muv) genes in C. briggsae that inhibit inappropriate division of vulval precursors. We have developed a new efficient and cost-effective gene mapping method to localize Muv mutations to small genetic intervals on chromosomes, thus facilitating cloning and functional studies. We demonstrate the utility of our method by determining molecular identities of three of the Muv genes that include orthologs of Cel-lin-1 (ETS) and Cel-lin-31 (Winged-Helix) of the EGF-Ras pathway and Cel-pry-1 (Axin), of the Wnt pathway. The remaining four genes reside in regions that lack orthologs of known C. elegans Muv genes. Inhibitor studies demonstrate that the Muv phenotype of all four new genes is dependent on the activity of the EGF pathway kinase, MEK. One of these, Cbr-lin(gu167), shows modest increase in the expression of Cbr-lin-3/EGF compared to wild type. These results argue that while Cbr-lin(gu167) may act upstream of Cbr-lin-3/EGF, the other three genes influence the EGF pathway downstream or in parallel to Cbr-lin-3. Overall, our findings demonstrate that the genetic program underlying a conserved developmental process includes both conserved and divergent functional contributions.

No abstract