Stable developmental patterns of gene expression without morphogen gradients.

Gene expression patterns in developing organisms are established by groups of cross-regulating target genes that are driven by morphogen gradients. As development progresses, morphogen activity is reduced, leaving the emergent pattern without stabilizing positional cues and at risk of rapid deterioration due to the inherently noisy biochemical processes at the cellular level. But remarkably, gene expression patterns remain spatially stable and reproducible over long developmental time spans in many biological systems. Here we combine spatial-stochastic simulations with an enhanced sampling method (Non-Stationary Forward Flux Sampling) and a recently developed stability theory to address how spatiotemporal integrity of a gene expression pattern is maintained in developing tissue lacking morphogen gradients. Using a minimal embryo model consisting of spatially coupled biochemical reactor volumes, we study a prototypical stripe pattern in which weak cross-repression between nearest neighbor expression domains alternates with strong repression between next-nearest neighbor domains, inspired by the gap gene system in the Drosophila embryo. We find that tuning of the weak repressive interactions to an optimal level can prolong stability of the expression patterns by orders of magnitude, enabling stable patterns over developmentally relevant times in the absence of morphogen gradients. The optimal parameter regime found in simulations of the embryo model closely agrees with the predictions of our coarse-grained stability theory. To elucidate the origin of stability, we analyze a reduced phase space defined by two measures of pattern asymmetry. We find that in the optimal regime, intact patterns are protected via restoring forces that counteract random perturbations and give rise to a metastable basin. Together, our results demonstrate that metastable attractors can emerge as a property of stochastic gene expression patterns even without system-wide positional cues, provided that the gene regulatory interactions shaping the pattern are optimally tuned.