SummaryProper lung functioning requires not only a correct structure of the conducting airway tree, but also the simultaneous development of smooth muscles and vasculature. Lung branching morphogenesis is strongly stereotyped and involves the recursive use of only three modes of branching. We have previously shown that the experimentally described interactions between Fibroblast growth factor (FGF)10, Sonic hedgehog (SHH) and Patched (Ptc) can give rise to a Turing mechanism that not only reproduces the experimentally observed wildtype branching pattern but also, in part counterintuitive, patterns in mutant mice. Here we show that, even though many proteins affect smooth muscle formation and the expression of Vegfa, an inducer of blood vessel formation, it is sufficient to add FGF9 to the FGF10/SHH/Ptc module to successfully predict simultaneously the emergence of smooth muscles in the clefts between growing lung buds, and Vegfa expression in the distal sub-epithelial mesenchyme. Our model reproduces the phenotype of both wildtype and relevant mutant mice, as well as the results of most culture conditions described in the literature.

BackgroundThe family of TGF-β ligands is large and its members are involved in many different signaling processes. These signaling processes strongly differ in type with TGF-β ligands eliciting both sustained or transient responses. Members of the TGF-β family can also act as morphogen and cellular responses would then be expected to provide a direct read-out of the extracellular ligand concentration. A number of different models have been proposed to reconcile these different behaviours. We were interested to define the set of minimal modifications that are required to change the type of signal processing in the TGF-β signaling network.ResultsTo define the key aspects for signaling plasticity we focused on the core of the TGF-β signaling network. With the help of a parameter screen we identified ranges of kinetic parameters and protein concentrations that give rise to transient, sustained, or oscillatory responses to constant stimuli, as well as those parameter ranges that enable a proportional response to time-varying ligand concentrations (as expected in the read-out of morphogens). A combination of a strong negative feedback and fast shuttling to the nucleus biases signaling to a transient rather than a sustained response, while oscillations were obtained if ligand binding to the receptor is weak and the turn-over of the I-Smad is fast. A proportional read-out required inefficient receptor activation in addition to a low affinity of receptor-ligand binding. We find that targeted modification of single parameters suffices to alter the response type. The intensity of a constant signal (i.e. the ligand concentration), on the other hand, affected only the strength but not the type of the response.ConclusionsThe architecture of the TGF-β pathway enables the observed signaling plasticity. The observed range of signaling outputs to TGF-β ligand in different cell types and under different conditions can be explained with differences in cellular protein concentrations and with changes in effective rate constants due to cross-talk with other signaling pathways. It will be interesting to uncover the exact cellular differences as well as the details of the cross-talks in future work.