Emergent spatiotemporal population dynamics with cell-length control of synthetic microbial consortia

The increased complexity of synthetic microbial biocircuits highlights the need for distributed cell functionality due to concomitant increases in metabolic and regulatory burdens imposed on single-strain topologies. Distributed systems, however, introduce additional challenges since consortium composition and spatiotemporal dynamics of constituent strains must be robustly controlled to achieve desired circuit behaviors. Here, we address these challenges with a modeling-based investigation of emergent spatiotemporal population dynamics using cell-length control in monolayer, two-strain bacterial consortia. We demonstrate that with dynamic control of a strain’s division length, nematic cell alignment in close-packed monolayers can be destabilized. We find that this destabilization confers an emergent, competitive advantage to smaller-length strains—but by mechanisms that differ depending on the spatial patterns of the population. We used complementary modeling approaches to elucidate underlying mechanisms: an agent-based model to simulate detailed mechanical and signaling interactions between the competing strains, and a reductive, stochastic lattice model to represent cell-cell interactions with a single rotational parameter. Our modeling suggests that spatial strain-fraction oscillations can be generated when cell-length control is coupled to quorum-sensing signaling in negative feedback topologies. Our research employs novel methods of population control and points the way to programming strain fraction dynamics in consortial synthetic biology.     

Adaptive molecular networks controlling chemotactic migration: dynamic inputs and selection of the network architecture

Eukaryotic signalling networks underlying the cell''s ability to sense the gradient of chemotactic cues frequently have the dual property of perfect adaptation to spatially homogeneous inputs, and persistent activation by inputs that are spatially graded. This property is also shared by bacterial chemotaxis networks, raising the question of whether these two types of chemotactic processes also have similar organization of the underlying biomolecular processes. Interestingly, perfect adaptation can only be achieved robustly by a handful of mechanisms, and while eukaryotic chemotactic networks appear to rely on one of these—the incoherent feed-forward loop, bacterial chemotaxis depends on another—the negative feedback loop. In this review, we discuss how this conclusion can be reached even if the details of the molecular networks are incompletely understood. Furthermore, we argue that the use of distinct network architectures is not accidental and may be a consequence of the nature of the signalling inputs and the limitations of the sensory properties of different cell types.     

Learning the rules of collective cell migration using deep attention networks

Collective, coordinated cellular motions underpin key processes in all multicellular organisms, yet it has been difficult to simultaneously express the ‘rules’ behind these motions in clear, interpretable forms that effectively capture high-dimensional cell-cell interaction dynamics in a manner that is intuitive to the researcher. Here we apply deep attention networks to analyze several canonical living tissues systems and present the underlying collective migration rules for each tissue type using only cell migration trajectory data. We use these networks to learn the behaviors of key tissue types with distinct collective behaviors—epithelial, endothelial, and metastatic breast cancer cells—and show how the results complement traditional biophysical approaches. In particular, we present attention maps indicating the relative influence of neighboring cells to the learned turning decisions of a ‘focal cell’–the primary cell of interest in a collective setting. Colloquially, we refer to this learned relative influence as ‘attention’, as it serves as a proxy for the physical parameters modifying the focal cell’s future motion as a function of each neighbor cell. These attention networks reveal distinct patterns of influence and attention unique to each model tissue. Endothelial cells exhibit tightly focused attention on their immediate forward-most neighbors, while cells in more expansile epithelial tissues are more broadly influenced by neighbors in a relatively large forward sector. Attention maps of ensembles of more mesenchymal, metastatic cells reveal completely symmetric attention patterns, indicating the lack of any particular coordination or direction of interest. Moreover, we show how attention networks are capable of detecting and learning how these rules change based on biophysical context, such as location within the tissue and cellular crowding. That these results require only cellular trajectories and no modeling assumptions highlights the potential of attention networks for providing further biological insights into complex cellular systems.     

Adjusting phenotypes by noise control

Genetically identical cells can show phenotypic variability. This is often caused by stochastic events that originate from randomness in biochemical processes involving in gene expression and other extrinsic cellular processes. From an engineering perspective, there have been efforts focused on theory and experiments to control noise levels by perturbing and replacing gene network components. However, systematic methods for noise control are lacking mainly due to the intractable mathematical structure of noise propagation through reaction networks. Here, we provide a numerical analysis method by quantifying the parametric sensitivity of noise characteristics at the level of the linear noise approximation. Our analysis is readily applicable to various types of noise control and to different types of system; for example, we can orthogonally control the mean and noise levels and can control system dynamics such as noisy oscillations. As an illustration we applied our method to HIV and yeast gene expression systems and metabolic networks. The oscillatory signal control was applied to p53 oscillations from DNA damage. Furthermore, we showed that the efficiency of orthogonal control can be enhanced by applying extrinsic noise and feedback. Our noise control analysis can be applied to any stochastic model belonging to continuous time Markovian systems such as biological and chemical reaction systems, and even computer and social networks. We anticipate the proposed analysis to be a useful tool for designing and controlling synthetic gene networks.     

Complexity and diversity in sparse code priors improve receptive field characterization of Macaque V1 neurons

System identification techniques—projection pursuit regression models (PPRs) and convolutional neural networks (CNNs)—provide state-of-the-art performance in predicting visual cortical neurons’ responses to arbitrary input stimuli. However, the constituent kernels recovered by these methods are often noisy and lack coherent structure, making it difficult to understand the underlying component features of a neuron’s receptive field. In this paper, we show that using a dictionary of diverse kernels with complex shapes learned from natural scenes based on efficient coding theory, as the front-end for PPRs and CNNs can improve their performance in neuronal response prediction as well as algorithmic data efficiency and convergence speed. Extensive experimental results also indicate that these sparse-code kernels provide important information on the component features of a neuron’s receptive field. In addition, we find that models with the complex-shaped sparse code front-end are significantly better than models with a standard orientation-selective Gabor filter front-end for modeling V1 neurons that have been found to exhibit complex pattern selectivity. We show that the relative performance difference due to these two front-ends can be used to produce a sensitive metric for detecting complex selectivity in V1 neurons.     

Modulating RhoA effectors induces transitions to oscillatory and more wavelike RhoA dynamics in Caenorhabditis elegans zygotes

Pulsatile RhoA dynamics underlie a wide range of cell and tissue behaviors. The circuits that produce these dynamics in different cells share common architectures based on fast positive and delayed negative feedback through F-actin, but they can produce very different spatiotemporal patterns of RhoA activity. However, the underlying causes of this variation remain poorly understood. Here we asked how this variation could arise through modulation of actin network dynamics downstream of active RhoA in early Caenorhabditis elegans embryos. We find that perturbing two RhoA effectors—formin and anillin—induce transitions from nonrecurrent focal pulses to either large noisy oscillatory pulses (formin depletion) or noisy oscillatory waves (anillin depletion). In both cases these transitions could be explained by changes in local F-actin levels and depletion dynamics, leading to changes in spatial and temporal patterns of RhoA inhibition. However, the underlying mechanisms for F-actin depletion are distinct, with different dependencies on myosin II activity. Thus, modulating actomyosin network dynamics could shape the spatiotemporal dynamics of RhoA activity for different physiological or morphogenetic functions.     

Model guided trait-specific co-expression network estimation as a new perspective for identifying molecular interactions and pathways

A wide variety of 1) parametric regression models and 2) co-expression networks have been developed for finding gene-by-gene interactions underlying complex traits from expression data. While both methodological schemes have their own well-known benefits, little is known about their synergistic potential. Our study introduces their methodological fusion that cross-exploits the strengths of individual approaches via a built-in information-sharing mechanism. This fusion is theoretically based on certain trait-conditioned dependency patterns between two genes depending on their role in the underlying parametric model. Resulting trait-specific co-expression network estimation method 1) serves to enhance the interpretation of biological networks in a parametric sense, and 2) exploits the underlying parametric model itself in the estimation process. To also account for the substantial amount of intrinsic noise and collinearities, often entailed by expression data, a tailored co-expression measure is introduced along with this framework to alleviate related computational problems. A remarkable advance over the reference methods in simulated scenarios substantiate the method’s high-efficiency. As proof-of-concept, this synergistic approach is successfully applied in survival analysis, with acute myeloid leukemia data, further highlighting the framework’s versatility and broad practical relevance.     

Interplay between Graph Topology and Correlations of Third Order in Spiking Neuronal Networks

The study of processes evolving on networks has recently become a very popular research field, not only because of the rich mathematical theory that underpins it, but also because of its many possible applications, a number of them in the field of biology. Indeed, molecular signaling pathways, gene regulation, predator-prey interactions and the communication between neurons in the brain can be seen as examples of networks with complex dynamics. The properties of such dynamics depend largely on the topology of the underlying network graph. In this work, we want to answer the following question: Knowing network connectivity, what can be said about the level of third-order correlations that will characterize the network dynamics? We consider a linear point process as a model for pulse-coded, or spiking activity in a neuronal network. Using recent results from theory of such processes, we study third-order correlations between spike trains in such a system and explain which features of the network graph (i.e. which topological motifs) are responsible for their emergence. Comparing two different models of network topology—random networks of Erdős-Rényi type and networks with highly interconnected hubs—we find that, in random networks, the average measure of third-order correlations does not depend on the local connectivity properties, but rather on global parameters, such as the connection probability. This, however, ceases to be the case in networks with a geometric out-degree distribution, where topological specificities have a strong impact on average correlations.     

Failure of Delayed Feedback Deep Brain Stimulation for Intermittent Pathological Synchronization in Parkinson’s Disease

Suppression of excessively synchronous beta-band oscillatory activity in the brain is believed to suppress hypokinetic motor symptoms of Parkinson’s disease. Recently, a lot of interest has been devoted to desynchronizing delayed feedback deep brain stimulation (DBS). This type of synchrony control was shown to destabilize the synchronized state in networks of simple model oscillators as well as in networks of coupled model neurons. However, the dynamics of the neural activity in Parkinson’s disease exhibits complex intermittent synchronous patterns, far from the idealized synchronous dynamics used to study the delayed feedback stimulation. This study explores the action of delayed feedback stimulation on partially synchronized oscillatory dynamics, similar to what one observes experimentally in parkinsonian patients. We employ a computational model of the basal ganglia networks which reproduces experimentally observed fine temporal structure of the synchronous dynamics. When the parameters of our model are such that the synchrony is unphysiologically strong, the feedback exerts a desynchronizing action. However, when the network is tuned to reproduce the highly variable temporal patterns observed experimentally, the same kind of delayed feedback may actually increase the synchrony. As network parameters are changed from the range which produces complete synchrony to those favoring less synchronous dynamics, desynchronizing delayed feedback may gradually turn into synchronizing stimulation. This suggests that delayed feedback DBS in Parkinson’s disease may boost rather than suppress synchronization and is unlikely to be clinically successful. The study also indicates that delayed feedback stimulation may not necessarily exhibit a desynchronization effect when acting on a physiologically realistic partially synchronous dynamics, and provides an example of how to estimate the stimulation effect.     

Unintentional Interpersonal Synchronization Represented as a Reciprocal Visuo-Postural Feedback System: A Multivariate Autoregressive Modeling Approach

People’s behaviors synchronize. It is difficult, however, to determine whether synchronized behaviors occur in a mutual direction—two individuals influencing one another—or in one direction—one individual leading the other, and what the underlying mechanism for synchronization is. To answer these questions, we hypothesized a non-leader-follower postural sway synchronization, caused by a reciprocal visuo-postural feedback system operating on pairs of individuals, and tested that hypothesis both experimentally and via simulation. In the behavioral experiment, 22 participant pairs stood face to face either 20 or 70 cm away from each other wearing glasses with or without vision blocking lenses. The existence and direction of visual information exchanged between pairs of participants were systematically manipulated. The time series data for the postural sway of these pairs were recorded and analyzed with cross correlation and causality. Results of cross correlation showed that postural sway of paired participants was synchronized, with a shorter time lag when participant pairs could see one another’s head motion than when one of the participants was blindfolded. In addition, there was less of a time lag in the observed synchronization when the distance between participant pairs was smaller. As for the causality analysis, noise contribution ratio (NCR), the measure of influence using a multivariate autoregressive model, was also computed to identify the degree to which one’s postural sway is explained by that of the other’s and how visual information (sighted vs. blindfolded) interacts with paired participants’ postural sway. It was found that for synchronization to take place, it is crucial that paired participants be sighted and exert equal influence on one another by simultaneously exchanging visual information. Furthermore, a simulation for the proposed system with a wider range of visual input showed a pattern of results similar to the behavioral results.     

Social structure of primate interaction networks facilitates the emergence of cooperation

Animal cooperation has puzzled biologists for a long time as its existence seems to contravene the basic notion of evolutionary biology that natural selection favours ‘selfish’ genes that promote only their own well-being. Evolutionary game theory has shown that cooperators can prosper in populations of selfish individuals if they occur in clusters, interacting more frequently with each other than with the selfish. Here we show that social networks of primates possess the necessary social structure to promote the emergence of cooperation. By simulating evolutionary dynamics of cooperative behaviour on interaction networks of 70 primate groups, we found that for most groups network reciprocity augmented the fixation probability for cooperation. The variation in the strength of this effect can be partly explained by the groups’ community modularity—a network measure for the groups’ heterogeneity. Thus, given selective update and partner choice mechanisms, network reciprocity has the potential to explain socially learned forms of cooperation in primate societies.     

Systems biology of cellular rhythms

Rhythms abound in biological systems, particularly at the cellular level where they originate from the feedback loops present in regulatory networks. Cellular rhythms can be investigated both by experimental and modeling approaches, and thus represent a prototypic field of research for systems biology. They have also become a major topic in synthetic biology. We review advances in the study of cellular rhythms of biochemical rather than electrical origin by considering a variety of oscillatory processes such as Ca++ oscillations, circadian rhythms, the segmentation clock, oscillations in p53 and NF-κB, synthetic oscillators, and the oscillatory dynamics of cyclin-dependent kinases driving the cell cycle. Finally we discuss the coupling between cellular rhythms and their robustness with respect to molecular noise.     

Network Events on Multiple Space and Time Scales in Cultured Neural Networks and in a Stochastic Rate Model

Cortical networks, in-vitro as well as in-vivo, can spontaneously generate a variety of collective dynamical events such as network spikes, UP and DOWN states, global oscillations, and avalanches. Though each of them has been variously recognized in previous works as expression of the excitability of the cortical tissue and the associated nonlinear dynamics, a unified picture of the determinant factors (dynamical and architectural) is desirable and not yet available. Progress has also been partially hindered by the use of a variety of statistical measures to define the network events of interest. We propose here a common probabilistic definition of network events that, applied to the firing activity of cultured neural networks, highlights the co-occurrence of network spikes, power-law distributed avalanches, and exponentially distributed ‘quasi-orbits’, which offer a third type of collective behavior. A rate model, including synaptic excitation and inhibition with no imposed topology, synaptic short-term depression, and finite-size noise, accounts for all these different, coexisting phenomena. We find that their emergence is largely regulated by the proximity to an oscillatory instability of the dynamics, where the non-linear excitable behavior leads to a self-amplification of activity fluctuations over a wide range of scales in space and time. In this sense, the cultured network dynamics is compatible with an excitation-inhibition balance corresponding to a slightly sub-critical regime. Finally, we propose and test a method to infer the characteristic time of the fatigue process, from the observed time course of the network’s firing rate. Unlike the model, possessing a single fatigue mechanism, the cultured network appears to show multiple time scales, signalling the possible coexistence of different fatigue mechanisms.     

Cell Lineages and the Logic of Proliferative Control

It is widely accepted that the growth and regeneration of tissues and organs is tightly controlled. Although experimental studies are beginning to reveal molecular mechanisms underlying such control, there is still very little known about the control strategies themselves. Here, we consider how secreted negative feedback factors (“chalones”) may be used to control the output of multistage cell lineages, as exemplified by the actions of GDF11 and activin in a self-renewing neural tissue, the mammalian olfactory epithelium (OE). We begin by specifying performance objectives—what, precisely, is being controlled, and to what degree—and go on to calculate how well different types of feedback configurations, feedback sensitivities, and tissue architectures achieve control. Ultimately, we show that many features of the OE—the number of feedback loops, the cellular processes targeted by feedback, even the location of progenitor cells within the tissue—fit with expectations for the best possible control. In so doing, we also show that certain distinctions that are commonly drawn among cells and molecules—such as whether a cell is a stem cell or transit-amplifying cell, or whether a molecule is a growth inhibitor or stimulator—may be the consequences of control, and not a reflection of intrinsic differences in cellular or molecular character.