Data-driven modelling of a gene regulatory network for cell fate decisions in the growing limb bud

Parameter optimization coupled with model selection is a convenient approach to infer gene regulatory networks from experimental gene expression data, but so far it has been limited to single cells or static tissues where growth is not significant. Here, we present a computational study in which we determine an optimal gene regulatory network from the spatiotemporal dynamics of gene expression patterns in a complex 2D growing tissue (non-isotropic and heterogeneous growth rates). We use this method to predict the regulatory mechanisms that underlie proximodistal (PD) patterning of the developing limb bud. First, we map the expression patterns of the PD markers Meis1, Hoxa11 and Hoxa13 into a dynamic description of the tissue movements that drive limb morphogenesis. Secondly, we use reverse-engineering to test how different gene regulatory networks can interpret the opposing gradients of fibroblast growth factors (FGF) and retinoic acid (RA) to pattern the PD markers. Finally, we validate and extend the best model against various previously published manipulative experiments, including exogenous application of RA, surgical removal of the FGF source and genetic ectopic expression of Meis1. Our approach identifies the most parsimonious gene regulatory network that can correctly pattern the PD markers downstream of FGF and RA. This network reveals a new model of PD regulation which we call the "crossover model", because the proximal morphogen (RA) controls the distal boundary of Hoxa11, while conversely the distal morphogens (FGFs) control the proximal boundary.