Entrainment of heterogeneous glycolytic oscillations in single cells

Lysine Demethylase 4A is a Centrosome Associated Protein Required for Centrosome Integrity and Genomic Stability

Centrosomes play a fundamental role in nucleating and organizing microtubules in the cell and are vital for faithful chromosome segregation and maintenance of genomic stability. Loss of structural or functional integrity of centrosomes causes genomic instability and is a driver of oncogenesis. The lysine demethylase 4A (KDM4A) is an epigenetic 'eraser' of chromatin methyl marks, which we show also localizes to the centrosome with single molecule resolution. We additionally discovered KDM4A demethylase enzymatic activity is required to maintain centrosome homeostasis, and is required for centrosome integrity, a new functionality unlinked to altered expression of genes regulating centrosome number. We find rather, that KDM4A interacts with both mother and daughter centriolar proteins to localize to the centrosome in all stages of mitosis. Loss of KDM4A results in supernumerary centrosomes and accrual of chromosome segregation errors including chromatin bridges and micronuclei, markers of genomic instability. In summary, these data highlight a novel role for an epigenetic 'eraser' regulating centrosome integrity, mitotic fidelity, and genomic stability at the centrosome.

Quantifying microbial robustness in dynamic environments using microfluidic single-cell cultivation

BACKGROUND

Microorganisms must respond to changes in their environment. Analysing the robustness of functions (i.e. performance stability) to such dynamic perturbations is of great interest in both laboratory and industrial settings. Recently, a quantification method capable of assessing the robustness of various functions, such as specific growth rate or product yield, across different conditions, time frames, and populations has been developed for microorganisms grown in a 96-well plate. In micro-titer-plates, environmental change is slow and undefined. Dynamic microfluidic single-cell cultivation (dMSCC) enables the precise maintenance and manipulation of microenvironments, while tracking single cells over time using live-cell imaging. Here, we combined dMSCC and a robustness quantification method to a pipeline for assessing performance stability to changes occurring within seconds or minutes.

RESULTS

Saccharomyces cerevisiae CEN.PK113-7D, harbouring a biosensor for intracellular ATP levels, was exposed to glucose feast-starvation cycles, with each condition lasting from 1.5 to 48 min over a 20 h period. A semi-automated image and data analysis pipeline was developed and applied to assess the performance and robustness of various functions at population, subpopulation, and single-cell resolution. We observed a decrease in specific growth rate but an increase in intracellular ATP levels with longer oscillation intervals. Cells subjected to 48 min oscillations exhibited the highest average ATP content, but the lowest stability over time and the highest heterogeneity within the population.

CONCLUSION

The proposed pipeline enabled the investigation of function stability in dynamic environments, both over time and within populations. The strategy allows for parallelisation and automation, and is easily adaptable to new organisms, biosensors, cultivation conditions, and oscillation frequencies. Insights on the microbial response to changing environments will guide strain development and bioprocess optimisation.

Dynamic fluctuations in a bacterial metabolic network

The operation of the central metabolism is typically assumed to be deterministic, but dynamics and high connectivity of the metabolic network make it potentially prone to generating fluctuations. However, time-resolved measurements of metabolite levels in individual cells that are required to characterize such fluctuations remained a challenge, particularly in small bacterial cells. Here we use single-cell metabolite measurements based on Förster resonance energy transfer, combined with computer simulations, to explore the real-time dynamics of the metabolic network of Escherichia coli. We observe that steplike exposure of starved E. coli to glycolytic carbon sources elicits large periodic fluctuations in the intracellular concentration of pyruvate in individual cells. These fluctuations are consistent with predicted oscillatory dynamics of E. coli metabolic network, and they are primarily controlled by biochemical reactions around the pyruvate node. Our results further indicate that fluctuations in glycolysis propagate to other cellular processes, possibly leading to temporal heterogeneity of cellular states within a population.

Principles, mechanisms and functions of entrainment in biological oscillators

Entrainment is a phenomenon in which two oscillators interact with each other, typically through physical or chemical means, to synchronize their oscillations. This phenomenon occurs in biology to coordinate processes from the molecular to organismal scale. Biological oscillators can be entrained within a single cell, between cells or to an external input. Using six illustrative examples of entrainable biological oscillators, we discuss the distinctions between entrainment and synchrony and explore features that contribute to a system's propensity to entrain. Entrainment can either enhance or reduce the heterogeneity of oscillations within a cell population, and we provide examples and mechanisms of each case. Finally, we discuss the known functions of entrainment and discuss potential functions from an evolutionary perspective.

Origins of oscillatory dynamics in the model of reactive oxygen species in the rhizosphere

Oscillatory processes are essential for normal functioning and survival of biological systems, and reactive oxygen species have a prominent role in many of them. A mechanism representing the dynamics of these species in the rhizosphere is analyzed using stoichiometric network analysis with the aim to determine its capabilities to simulate various dynamical states, including oscillations. A detailed analysis has shown that unstable steady states result from four destabilizing feedback cycles, among which the cycle involving hydroquinone, an electron acceptor, and its semi-reduced form is the dominant one responsible for the existence of saddle-node and Andronov-Hopf bifurcations. This requires a higher steady-state concentration for the reduced electron acceptor compared to that of the remaining species, where the level of oxygen steady-state concentration determines whether the Andronov-Hopf or saddle-node bifurcation will occur.

Synchronisation of glycolytic activity in yeast cells

Glycolysis is the central metabolic pathway of almost every cell and organism. Under appropriate conditions, glycolytic oscillations may occur in individual cells as well as in entire cell populations or tissues. In many biological systems, glycolytic oscillations drive coherent oscillations of other metabolites, for instance in cardiomyocytes near anorexia, or in pancreas where they lead to a pulsatile release of insulin. Oscillations at the population or tissue level require the cells to synchronize their metabolism. We review the progress achieved in studying a model organism for glycolytic oscillations, namely yeast. Oscillations may occur on the level of individual cells as well as on the level of the cell population. In yeast, the cell-to-cell interaction is realized by diffusion-mediated intercellular communication via a messenger molecule. The present mini-review focuses on the synchronisation of glycolytic oscillations in yeast. Synchronisation is a quorum-sensing phenomenon because the collective oscillatory behaviour of a yeast cell population ceases when the cell density falls below a threshold. We review the question, under which conditions individual cells in a sparse population continue or cease to oscillate. Furthermore, we provide an overview of the pathway leading to the onset of synchronized oscillations. We also address the effects of spatial inhomogeneities (e.g., the formation of spatial clusters) on the collective dynamics, and also review the emergence of travelling waves of glycolytic activity. Finally, we briefly review the approaches used in numerical modelling of synchronized cell populations.

Intercellular communication induces glycolytic synchronization waves between individually oscillating cells

Many organs have internal structures with spatially differentiated and sometimes temporally synchronized groups of cells. The mechanisms leading to such differentiation and coordination are not well understood. Here we design a diffusion-limited microfluidic system to mimic a multicellular organ structure with peripheral blood flow and test whether a group of individually oscillating yeast cells could form subpopulations of spatially differentiated and temporally synchronized cells. Upon substrate addition, the dynamic response at single-cell level shows glycolytic oscillations, leading to wave fronts traveling through the monolayered population and to synchronized communities at well-defined positions in the cell chamber. A detailed mechanistic model with the architectural structure of the flow chamber incorporated successfully predicts the spatial-temporal experimental data, and allows for a molecular understanding of the observed phenomena. The intricate interplay of intracellular biochemical reaction networks leading to the oscillations, combined with intercellular communication via metabolic intermediates and fluid dynamics of the reaction chamber, is responsible for the generation of the subpopulations of synchronized cells. This mechanism, as analyzed from the model simulations, is experimentally tested using different concentrations of cyanide stress solutions. The results are reproducible and stable, despite cellular heterogeneity, and the spontaneous community development is reminiscent of a zoned cell differentiation often observed in multicellular organs.

Partial synchronisation of glycolytic oscillations in yeast cell populations

The transition between synchronized and asynchronous behaviour of immobilized yeast cells of the strain Saccharomyces carlsbergensis was investigated by monitoring the autofluorescence of the coenzyme NADH. In populations of intermediate cell densities the individual cells remained oscillatory, whereas on the level of the cell population both a partially synchronized and an asynchronous state were accessible for experimental studies. In the partially synchronized state, the mean oscillatory frequency was larger than that of the cells in the asynchronous state. This suggests that synchronisation occurred due to entrainment by the cells that oscillated more rapidly. This is typical for synchronisation due to phase advancement. Furthermore, the synchronisation of the frequency of the glycolytic oscillations preceded the synchronisation of their phases. However, the cells did not synchronize completely, as the distribution of the oscillatory frequencies only narrowed but did not collapse to a unique frequency. Cells belonging to spatially denser clusters showed a slightly enhanced local synchronisation during the episode of partial synchronisation. Neither the clusters nor a transition from partially synchronized glycolytic oscillations to travelling glycolytic waves did substantially affect the degree of partial synchronisation. Chimera states, i.e., the coexistence of a synchronized and an asynchronous part of the population, could not be found.

Emergence of collective oscillations in adaptive cells

Collective oscillations of cells in a population appear under diverse biological contexts. Here, we establish a set of common principles by categorising the response of individual cells against a time-varying signal. A positive intracellular signal relay of sufficient gain from participating cells is required to sustain the oscillations, together with phase matching. The two conditions yield quantitative predictions for the onset cell density and frequency in terms of measured single-cell and signal response functions. Through mathematical constructions, we show that cells that adapt to a constant stimulus fulfil the phase requirement by developing a leading phase in an active frequency window that enables cell-to-signal energy flow. Analysis of dynamical quorum sensing in several cellular systems with increasing biological complexity reaffirms the pivotal role of adaptation in powering oscillations in an otherwise dissipative cell-to-cell communication channel. The physical conditions identified also apply to synthetic oscillatory systems.

There are many examples of cell populations exhibiting density-dependent collective oscillatory behaviour. Here, the authors show that sustained collective oscillations emerge when cells anticipate variation in signal and attempt to amplify it, a property that can be linked to adaptation.

Phosphofructokinase controls the acetaldehyde-induced phase shift in isolated yeast glycolytic oscillators

The response of oscillatory systems to external perturbations is crucial for emergent properties such as synchronisation and phase locking and can be quantified in a phase response curve (PRC). In individual, oscillating yeast cells, we characterised experimentally the phase response of glycolytic oscillations for external acetaldehyde pulses and followed the transduction of the perturbation through the system. Subsequently, we analysed the control of the relevant system components in a detailed mechanistic model. The observed responses are interpreted in terms of the functional coupling and regulation in the reaction network. We find that our model quantitatively predicts the phase-dependent phase shift observed in the experimental data. The phase shift is in agreement with an adaptation leading to synchronisation with an external signal. Our model analysis establishes that phosphofructokinase plays a key role in the phase shift dynamics as shown in the PRC and adaptation time to external perturbations. Specific mechanism-based interventions, made possible through such analyses of detailed models, can improve upon standard trial and error methods, e.g. melatonin supplementation to overcome jet-lag, which are error-prone, specifically, since the effects are phase dependent and dose dependent. The models by Gustavsson and Goldbeter discussed in the text can be obtained from the JWS Online simulation database: (https://jjj.bio.vu.nl/models/gustavsson5 and https://jjj.bio.vu.nl/models/goldbeter1).

Biological Oscillators in Nanonetworks-Opportunities and Challenges

One of the major issues in molecular communication-based nanonetworks is the provision and maintenance of a common time knowledge. To stay true to the definition of molecular communication, biological oscillators are the potential solutions to achieve that goal as they generate oscillations through periodic fluctuations in the concentrations of molecules. Through the lens of a communication systems engineer, the scope of this survey is to explicitly classify, for the first time, existing biological oscillators based on whether they are found in nature or not, to discuss, in a tutorial fashion, the main principles that govern the oscillations in each oscillator, and to analyze oscillator parameters that are most relevant to communication engineer researchers. In addition, the survey highlights and addresses the key open research issues pertaining to several physical aspects of the oscillators and the adoption and implementation of the oscillators to nanonetworks. Moreover, key research directions are discussed.

Sensitivity analysis of the noise-induced oscillatory multistability in Higgins model of glycolysis

A phenomenon of the noise-induced oscillatory multistability in glycolysis is studied. As a basic deterministic skeleton, we consider the two-dimensional Higgins model. The noise-induced generation of mixed-mode stochastic oscillations is studied in various parametric zones. Probabilistic mechanisms of the stochastic excitability of equilibria and noise-induced splitting of randomly forced cycles are analysed by the stochastic sensitivity function technique. A parametric zone of supersensitive Canard-type cycles is localized and studied in detail. It is shown that the generation of mixed-mode stochastic oscillations is accompanied by the noise-induced transitions from order to chaos.

Effect of the periodic properties of toxic stress on the oscillatory behaviour of glycolysis in yeast-evidence of a toxic effect frequency

Starving and nondividing yeast cells induce changes in the electron donor nicotinamide adenine dinucleotide (NADH) levels in a cyclic and wave-like manner for over 90min. Yeast suspensions were used to examine the toxic effects of contaminants on the cyclic behaviour of metabolite changes during anaerobic glycolysis. The cyclic behaviour NADH levels in yeast cell suspensions starved for 2 to 5h was studied after the addition of 10mM glucose for 5min followed by 10mM KCN to block aerobic glycolysis. The effects of three toxic elements (CuSO₄, silver nanoparticles-nAg, and GdCl₃), known for their potential to alter glycolsysis, on NADH levels over time were examined during the 3-h starvation step. The data were analyzed using spectral analysis (Fourier transformation) to characterize the cyclic behaviour of NADH levels during anaerobic glycolysis. Increasing the starvation time by 3h increased the amplitude of changes in NADH levels with characteristic periods of 3 to 8min. Longer starvation times decreased the amplitude of oscillations during these periods, with the appearance of NADH changes at higher frequencies. Moreover, the amplitude changes in NADH were proportional to the starvation time. Exposure to the above chemicals during the 3-h starvation time led to the formation of higher frequencies with concentration-dependent amplitude changes. In comparison with nAg and Gd³⁺, Cu²⁺ was the most toxic (decreased viability the most) and produce changes at higher frequencies as well. It is noteworthy that each element produced a characteristic change in the frequency profiles, which suggests different mechanisms of action in which the severity of toxicity shifted NADH changes at higher frequencies. In conclusion, the appearance of synchronized oscillations in dense yeast populations following synchronization stress could be induced by starvation and exposure to chemicals. However, synchronicity could be abolished when cells desynchronize as a result of loss of cell viability, which contributes to heterogeneity in yeast populations, translating into NADH changes at higher frequencies. This is the first report on the influence of environmental contaminants on the cyclic or wave-like behaviour of biochemical changes in cells.

Entrainment of a Bacterial Synthetic Gene Oscillator through Proteolytic Queueing

Internal chemical oscillators (chemical clocks) direct the behavior of numerous biological systems, and maintenance of a given period and phase among many such oscillators may be important for their proper function. However, both environmental variability and fundamental molecular noise can cause biochemical oscillators to lose coherence. One solution to maintaining coherence is entrainment, where an external signal provides a cue that resets the phase of the oscillators. In this work, we study the entrainment of gene networks by a queueing interaction established by competition between proteins for a common proteolytic pathway. Principles of queueing entrainment are investigated for an established synthetic oscillator in Escherichia coli. We first explore this theoretically using a standard chemical reaction network model and a map-based model, both of which suggest that queueing entrainment can be achieved through pulsatile production of an additional protein competing for a common degradation pathway with the oscillator proteins. We then use a combination of microfluidics and fluorescence microscopy to verify that pulse trains modulating the production rate of a fluorescent protein targeted to the same protease (ClpXP) as the synthetic oscillator can entrain the oscillator.

Modelling chronotaxicity of cellular energy metabolism to facilitate the identification of altered metabolic states

Altered cellular energy metabolism is a hallmark of many diseases, one notable example being cancer. Here, we focus on the identification of the transition from healthy to abnormal metabolic states. To do this, we study the dynamics of energy production in a cell. Due to the thermodynamic openness of a living cell, the inability to instantaneously match fluctuating supply and demand in energy metabolism results in nonautonomous time-varying oscillatory dynamics. However, such oscillatory dynamics is often neglected and treated as stochastic. Based on experimental evidence of metabolic oscillations, we show that changes in metabolic state can be described robustly by alterations in the chronotaxicity of the corresponding metabolic oscillations, i.e. the ability of an oscillator to resist external perturbations. We also present a method for the identification of chronotaxicity, applicable to general oscillatory signals and, importantly, apply this to real experimental data. Evidence of chronotaxicity was found in glycolytic oscillations in real yeast cells, verifying that chronotaxicity could be used to study transitions between metabolic states.

The 6th International Conference on Analysis of Microbial Cells at the Single Cell Level, Retz, Austria, 19-22 July 2015

The 6th International Conference on Analysis of Microbial Cells at the Single Cell Level, held in Retz, Austria from 19 to 22 July 2015, brought together experts from different areas working with bacterial, yeast and mammalian cell systems. The conference highlighted the importance of dissecting cell behaviour down to the single cell level, as analysis of mixed populations can obscure crucial cell-to-cell variations. The sessions covered advances in the fields of image analysis and microscopy, flow cytometry and cell sorting as well as bioinformatics, including recent developments and new applications of existing tools. In addition, a high speed poster slam session contributed to the lively discussions and exchange of expertise among academic and industrial researchers.