Before revision:  This work uncovers a mechanism that underlies this emergence of multiple modes under slow slow promoter kinetics, and studies it theoretically. Mathematical tools from singular perturbation theory are employed in order to analytically characterize stationary distributions of Chemical Master Equations for GRN's, in the limit of slow switching, as a mixture of Poisson distributions This approach, which may be interpreted as a finite-dimensional reduction of a countable Markov chain, offers a rigorous framework to explain phenomena such as non-genetic population heterogeneity and transcriptional bursting .
After revision:  This work uncovers a mechanism that underlies this emergence of multiple modes . A general formal approach gives insight into the influence of parameters and the prediction of how changes in GRN wiring, for example through mutations or artificial interventions, impact the possible number, location, and likelihood of alternative states. We adapt tools from the mathematical field of singular perturbation theory are employed in order to analytically characterize stationary distributions of Chemical Master Equations for GRN's, in the limit of slow switching, as a mixture of Poisson distributions This approach, which may be interpreted as a finite-dimensional reduction of a countable Markov chain, offers a rigorous framework to explain phenomena such as non-genetic population heterogeneity and transcriptional bursting .