Simultaneous recordings of gene network dynamics across large populations have revealed that cell characteristics vary considerably even in clonal lines. Inferring the variability of parameters that determine gene dynamics is key to understanding cellular behavior. However, this is complicated by the fact that the outcomes and effects of many reactions are not observable directly. Unobserved reactions can be replaced with time delays to reduce model dimensionality and simplify inference. However, the resulting models are non-Markovian, and require the development of new inference techniques. We propose a non-Markovian, hierarchical Bayesian inference framework for quantifying the variability of cellular processes within and across cells in a population. We illustrate our approach using a delayed birth-death process. In general, a distributed delay model, rather than a popular fixed delay model, is needed for inference, even if only mean reaction delays are of interest. Using in silico and experimental data we show that the proposed hierarchical framework is robust and leads to improved estimates compared to its non-hierarchical counterpart. We apply our method to data obtained using time-lapse microscopy and infer the parameters that describe the dynamics of protein production at the single cell and population level. The mean delays in protein production are larger than previously reported, have a coefficient of variation of around 0.2 across the population, and are not strongly correlated with protein production or growth rates.

The epigenetic landscape illustrates how cells differentiate into different types through the control of gene regulatory networks. Numerous studies have investigated epigenetic gene regulation but there are limited studies on how the epigenetic landscape and the presence of pathogens influence the evolution of host traits. Here we formulate a multistable decision-switch model involving many possible phenotypes with the antagonistic influence of parasitism. As expected, pathogens can drive dominant (common) phenotypes to become inferior, such as through negative frequency-dependent selection. Furthermore, novel predictions of our model show that parasitism can steer the dynamics of phenotype specification from multistable equilibrium convergence to oscillations. This oscillatory behavior could explain pathogen-mediated epimutations and excessive phenotypic plasticity. The Red Queen dynamics also occur in certain parameter space of the model, which demonstrates winnerless cyclic phenotype-switching in hosts and in pathogens. The results of our simulations elucidate how epigenetic landscape is associated with the phenotypic fitness landscape and how parasitism facilitates non-genetic phenotypic diversity.

This paper studies linear conjugacy of PL-RDK systems, which are kinetic systems with power law rate functions whose kinetic orders are identical for branching reactions, i.e. reactions with the same reactant complex. Mass action kinetics (MAK) systems are the best known examples of such systems with reactant-determined kinetic orders (RDK). We specify their kinetics with their rate vector and T matrix. The T matrix is formed from the kinetic order matrix by replacing the reactions with their reactant complexes as row indices (thus compressing identical rows of branching reactions of a reactant complex to one) and taking the transpose of the resulting matrix. The T matrix is hence the kinetic analogue of the network's matrix of complexes Y with the latter's columns of non-reactant complexes truncated away. For MAK systems, the T matrix and the truncated Y matrix are identical. We show that, on non-branching networks, a necessary condition for linear conjugacy of MAK systems and, more generally, of PL-FSK (power law factor span surjective kinetics) systems, i.e. those whose T matrix columns are pairwise different, is T = T , i.e. equality of their T matrices. This motivated our inclusion of the condition T = T in exploring extension

The epigenetic landscape illustrates how cells differentiate through the control of gene regulatory networks. Numerous studies have investigated epigenetic gene regulation but there are limited studies on how the epigenetic landscape and the presence of pathogens influence the evolution of host traits. Here, we formulate a multistable decision-switch model involving several phenotypes with the antagonistic influence of parasitism. As expected, pathogens can drive dominant (common) phenotypes to become inferior through negative frequency-dependent selection. Furthermore, novel predictions of our model show that parasitism can steer the dynamics of phenotype specification from multistable equilibrium convergence to oscillations. This oscillatory behavior could explain pathogen-mediated epimutations and excessive phenotypic plasticity. The Red Queen dynamics also occur in certain parameter space of the model, which demonstrates winnerless cyclic phenotype-switching in hosts and in pathogens. The results of our simulations elucidate the association between the epigenetic and phenotypic fitness landscapes and how parasitism facilitates non-genetic phenotypic diversity.