  Gene expression patterns are established by cross-regulating target genes
that interpret morphogen gradients. However, as development progresses,
morphogen activity is reduced, leaving the emergent pattern without stabilizing
positional cues. The pattern then can be deteriorated by the intrinsically
noisy biochemical processes acting at the cellular level. But remarkably, the
established gene expression patterns remain spatially and temporally stable in
many biological systems. Here we combine spatial-stochastic simulations with an
enhanced sampling method 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 stripe
pattern in which weak cross-repression between nearest neighbor domians
alternates with strong repression between next-nearest neighbor domains,
inspired by the gap gene system in the Drosophila embryo. We find that
fine-tuning of the weak repressive interactions to an optimal level can
significantly increase temporal stability of the expression patterns, allowing
for stable patterns over developmentally relevant times in the absence of
morphogen gradients. The numerically determined optimal parameters closely
agree with the predictions of the stability theory. By analizing the dynamics
of factors characterizing pattern integrity, we trace back the stability
enhancement to the emergence of restoring forces, maintaining the pattern in a
meta-stable basin. Altogether 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.
