Longitudinal hair cortisol in bipolar disorder and a mechanism based on HPA dynamics.

Bipolar disorder (BD) is marked by fluctuating mood states over months to years, often with elevated cortisol levels. Elevated cortisol can also trigger mood episodes. Here, we combine longitudinal hair cortisol and mood measurements with mathematical modeling to provide a potential mechanistic link between cortisol and mood timescales in BD. Using 12 cm hair samples, representing a year of growth, we found enhanced year-scale cortisol fluctuations whose amplitude averaged 4-fold higher in BD (n = 26) participants than controls (n = 59). The proximal 2 cm of hair correlated with recent mood scores. Depression (n = 266) and mania (n = 273) scores from a longitudinal study of BD showed similar frequency spectra. These results suggest a mechanism for BD in which high emotional reactivity excites the slow timescales in the hypothalamic-pituitary-adrenal (HPA) axis to generate elevated months-scale cortisol fluctuations, triggering cortisol-induced mood episodes.

A synthetic differentiation circuit in Escherichia coli for suppressing mutant takeover.

Differentiation is crucial for multicellularity. However, it is inherently susceptible to mutant cells that fail to differentiate. These mutants outcompete normal cells by excessive self-renewal. It remains unclear what mechanisms can resist such mutant expansion. Here, we demonstrate a solution by engineering a synthetic differentiation circuit in Escherichia coli that selects against these mutants via a biphasic fitness strategy. The circuit provides tunable production of synthetic analogs of stem, progenitor, and differentiated cells. It resists mutations by coupling differentiation to the production of an essential enzyme, thereby disadvantaging non-differentiating mutants. The circuit selected for and maintained a positive differentiation rate in long-term evolution. Surprisingly, this rate remained constant across vast changes in growth conditions. We found that transit-amplifying cells (fast-growing progenitors) underlie this environmental robustness. Our results provide insight into the stability of differentiation and demonstrate a powerful method for engineering evolutionarily stable multicellular consortia.

Tradeoffs in bacterial physiology determine the efficiency of antibiotic killing.

Antibiotic effectiveness depends on a variety of factors. While many mechanistic details of antibiotic action are known, the connection between death rate and bacterial physiology is poorly understood. A common observation is that death rate in antibiotics rises linearly with growth rate; however, it remains unclear how other factors, such as environmental conditions and whole-cell physiological properties, affect bactericidal activity. To address this, we developed a high-throughput assay to precisely measure antibiotic-mediated death. We found that death rate is linear in growth rate, but the slope depends on environmental conditions. Growth under stress lowers death rate compared to nonstressed environments with similar growth rate. To understand stress's role, we developed a mathematical model of bacterial death based on resource allocation that includes a stress-response sector; we identify this sector using RNA-seq. Our model accurately predicts the minimal inhibitory concentration (MIC) with zero free parameters across a wide range of growth conditions. The model also quantitatively predicts death and MIC when sectors are experimentally modulated using cyclic adenosine monophosphate (cAMP), including protection from death at very low cAMP levels. The present study shows that different conditions with equal growth rate can have different death rates and establishes a quantitative relation between growth, death, and MIC that suggests approaches to improve antibiotic efficacy.

Timescales of Human Hair Cortisol Dynamics.

Cortisol is a major human stress hormone, secreted within minutes of acute stress. Cortisol also has slower patterns of variation: a strong circadian rhythm and a seasonal rhythm. However, longitudinal cortisol dynamics in healthy individuals over timescales of months has rarely been studied. Here, we measured longitudinal cortisol in 55 healthy participants using 12 cm of hair, which provides a retrospective measurement over one year. Individuals showed (non-seasonal) fluctuations averaging about 22% around their baseline. Fourier analysis reveals dominant slow frequencies with periods of months to a year. These frequencies can be explained by a mathematical model of the hormonal cascade that controls cortisol, the HPA axis, when including the slow timescales of tissue turnover of the glands. Measuring these dynamics is important for understanding disorders in which cortisol secretion is impaired over months, such as mood disorders, and to test models of cortisol feedback control.

A Bacterial Growth Law out of Steady State.

Bacterial growth follows simple laws in constant conditions. However, bacteria in nature often face fluctuating environments. We therefore ask whether there are growth laws that apply to changing environments. We derive a law for upshifts using an optimal resource-allocation model: the post-shift growth rate equals the geometrical mean of the pre-shift growth rate and the growth rate on saturating carbon. We test this using chemostat and batch culture experiments, as well as previous data from several species. The increase in growth rate after an upshift indicates that ribosomes have spare capacity (SC). We demonstrate theoretically that SC has the cost of slow steady-state growth but is beneficial after an upshift because it prevents large overshoots in intracellular metabolites and allows rapid response to change. We also provide predictions for downshifts. The present study quantifies the optimal degree of SC, which rises the slower the growth rate, and suggests that SC can be precisely regulated.

Dynamic Proteomics of Herpes Simplex Virus Infection.

The cellular response to viral infection is usually studied at the level of cell populations. Currently, it remains an open question whether and to what extent cell-to-cell variability impacts the course of infection. Here we address this by dynamic proteomics-imaging and tracking 400 yellow fluorescent protein (YFP)-tagged host proteins in individual cells infected by herpes simplex virus 1. By quantifying time-lapse fluorescence imaging, we analyze how cell-to-cell variability impacts gene expression from the viral genome. We identify two proteins, RFX7 and geminin, whose levels at the time of infection correlate with successful initiation of gene expression. These proteins are cell cycle markers, and we find that the position in the cell cycle at the time of infection (along with the cell motility and local cell density) can reasonably predict in which individual cells gene expression from the viral genome will commence. We find that the onset of cell division dramatically impacts the progress of infection, with 70% of dividing cells showing no additional gene expression after mitosis. Last, we identify four host proteins that are specifically modulated in infected cells, of which only one has been previously recognized. SUMO2 and RPAP3 levels are rapidly reduced, while SLTM and YTHDC1 are redistributed to form nuclear foci. These modulations are dependent on the expression of ICP0, as shown by infection with two mutant viruses that lack ICP0. Taken together, our results provide experimental validation for the long-held notion that the success of infection is dependent on the state of the host cell at the time of infection.IMPORTANCE High-throughput assays have revolutionized many fields in biology, both by allowing a more global understanding of biological processes and by deciphering rare events in subpopulations. Here we use such an assay, dynamic proteomics, to study viral infection at the single-cell level. We follow tens of thousands of individual cells infected by herpes simplex virus using fluorescence live imaging. Our results link the state of a cell at the time of virus infection with its probability to successfully initiate gene expression from the viral genome. Further, we identified three cellular proteins that were previously unknown to respond to viral infection. We conclude that dynamic proteomics provides a powerful tool to study single-cell differences during viral infection.

Optimality and sub-optimality in a bacterial growth law.

Organisms adjust their gene expression to improve fitness in diverse environments. But finding the optimal expression in each environment presents a challenge. We ask how good cells are at finding such optima by studying the control of carbon catabolism genes in Escherichia coli. Bacteria show a growth law: growth rate on different carbon sources declines linearly with the steady-state expression of carbon catabolic genes. We experimentally modulate gene expression to ask if this growth law always maximizes growth rate, as has been suggested by theory. We find that the growth law is optimal in many conditions, including a range of perturbations to lactose uptake, but provides sub-optimal growth on several other carbon sources. Combining theory and experiment, we genetically re-engineer E. coli to make sub-optimal conditions into optimal ones and vice versa. We conclude that the carbon growth law is not always optimal, but represents a practical heuristic that often works but sometimes fails.

Glucose becomes one of the worst carbon sources for E.coli on poor nitrogen sources due to suboptimal levels of cAMP.

In most conditions, glucose is the best carbon source for E. coli: it provides faster growth than other sugars, and is consumed first in sugar mixtures. Here we identify conditions in which E. coli strains grow slower on glucose than on other sugars, namely when a single amino acid (arginine, glutamate, or proline) is the sole nitrogen source. In sugar mixtures with these nitrogen sources, E. coli still consumes glucose first, but grows faster rather than slower after exhausting glucose, generating a reversed diauxic shift. We trace this counterintuitive behavior to a metabolic imbalance: levels of TCA-cycle metabolites including α-ketoglutarate are high, and levels of the key regulatory molecule cAMP are low. Growth rates were increased by experimentally increasing cAMP levels, either by adding external cAMP, by genetically perturbing the cAMP circuit or by inhibition of glucose uptake. Thus, the cAMP control circuitry seems to have a 'bug' that leads to slow growth under what may be an environmentally rare condition.

Hierarchy of non-glucose sugars in Escherichia coli.

BACKGROUND: Understanding how cells make decisions, and why they make the decisions they make, is of fundamental interest in systems biology. To address this, we study the decisions made by E. coli on which genes to express when presented with two different sugars. It is well-known that glucose, E. coli's preferred carbon source, represses the uptake of other sugars by means of global and gene-specific mechanisms. However, less is known about the utilization of glucose-free sugar mixtures which are found in the natural environment of E. coli and in biotechnology. RESULTS: Here, we combine experiment and theory to map the choices of E. coli among 6 different non-glucose carbon sources. We used robotic assays and fluorescence reporter strains to make precise measurements of promoter activity and growth rate in all pairs of these sugars. We find that the sugars can be ranked in a hierarchy: in a mixture of a higher and a lower sugar, the lower sugar system shows reduced promoter activity. The hierarchy corresponds to the growth rate supported by each sugar- the faster the growth rate, the higher the sugar on the hierarchy. The hierarchy is 'soft' in the sense that the lower sugar promoters are not completely repressed. Measurement of the activity of the master regulator CRP-cAMP shows that the hierarchy can be quantitatively explained based on differential activation of the promoters by CRP-cAMP. Comparing sugar system activation as a function of time in sugar pair mixtures at sub-saturating concentrations, we find cases of sequential activation, and also cases of simultaneous expression of both systems. Such simultaneous expression is not predicted by simple models of growth rate optimization, which predict only sequential activation. We extend these models by suggesting multi-objective optimization for both growing rapidly now and preparing the cell for future growth on the poorer sugar. CONCLUSION: We find a defined hierarchy of sugar utilization, which can be quantitatively explained by differential activation by the master regulator cAMP-CRP. The present approach can be used to understand cell decisions when presented with mixtures of conditions.

Linear superposition and prediction of bacterial promoter activity dynamics in complex conditions.

Bacteria often face complex environments. We asked how gene expression in complex conditions relates to expression in simpler conditions. To address this, we obtained accurate promoter activity dynamical measurements on 94 genes in E. coli in environments made up of all possible combinations of four nutrients and stresses. We find that the dynamics across conditions is well described by two principal component curves specific to each promoter. As a result, the promoter activity dynamics in a combination of conditions is a weighted average of the dynamics in each condition alone. The weights tend to sum up to approximately one. This weighted-average property, called linear superposition, allows predicting the promoter activity dynamics in a combination of conditions based on measurements of pairs of conditions. If these findings apply more generally, they can vastly reduce the number of experiments needed to understand how E. coli responds to the combinatorially huge space of possible environments.

Promoters maintain their relative activity levels under different growth conditions.

Most genes change expression levels across conditions, but it is unclear which of these changes represents specific regulation and what determines their quantitative degree. Here, we accurately measured activities of ~900 S. cerevisiae and ~1800 E. coli promoters using fluorescent reporters. We show that in both organisms 60-90% of promoters change their expression between conditions by a constant global scaling factor that depends only on the conditions and not on the promoter's identity. Quantifying such global effects allows precise characterization of specific regulation-promoters deviating from the global scale line. These are organized into few functionally related groups that also adhere to scale lines and preserve their relative activities across conditions. Thus, only several scaling factors suffice to accurately describe genome-wide expression profiles across conditions. We present a parameter-free passive resource allocation model that quantitatively accounts for the global scaling factors. It suggests that many changes in expression across conditions result from global effects and not specific regulation, and provides means for quantitative interpretation of expression profiles.

The last generation of bacterial growth in limiting nutrient.

BACKGROUND: Bacterial growth as a function of nutrients has been studied for decades, but is still not fully understood. In particular, the growth laws under dynamically changing environments have been difficult to explore, because of the rapidly changing conditions. Here, we address this challenge by means of a robotic assay and measure bacterial growth rate, promoter activity and substrate level at high temporal resolution across the entire growth curve in batch culture. As a model system, we study E. coli growing under nitrogen or carbon limitation, and explore the dynamics in the last generation of growth where nutrient levels can drop rapidly. RESULTS: We find that growth stops abruptly under limiting nitrogen or carbon, but slows gradually when nutrients are not limiting. By measuring growth rate at a 3 min time resolution, and inferring the instantaneous substrate level, s, we find that the reduction in growth rate µ under nutrient limitation follows Monod's law, µ=µ0(s/(k(s)+s)). By following promoter activity of different genes we found that the abrupt stop of growth under nitrogen or carbon limitation is accompanied by a pulse-like up-regulation of the expression of genes in the relevant nutrient assimilation pathways. We further find that sharp stop of growth is conditional on the presence of regulatory proteins in the assimilation pathway. CONCLUSIONS: The observed sharp stop of growth accompanied by a pulsed expression of assimilation genes allows bacteria to compensate for the drop in nutrients, suggesting a strategy used by the cells to prolong exponential growth under limiting substrate.

Promoter activity dynamics in the lag phase of Escherichia coli.

BACKGROUND: Lag phase is a period of time with no growth that occurs when stationary phase bacteria are transferred to a fresh medium. Bacteria in lag phase seem inert: their biomass does not increase. The low number of cells and low metabolic activity make it difficult to study this phase. As a consequence, it has not been studied as thoroughly as other bacterial growth phases. However, lag phase has important implications for bacterial infections and food safety. We asked which, if any, genes are expressed in the lag phase of Escherichia coli, and what is their dynamic expression pattern. RESULTS: We developed an assay based on imaging flow cytometry of fluorescent reporter cells that overcomes the challenges inherent in studying lag phase. We distinguish between lag1 phase- in which there is no biomass growth, and lag2 phase--in which there is biomass growth but no cell division. We find that in lag1 phase, most promoters are not active, except for the enzymes that utilize the specific carbon source in the medium. These genes show promoter activities that increase exponentially with time, despite the fact that the cells do not measurably increase in size. An oxidative stress promoter, katG, is also active. When cells enter lag2 and begin to grow in size, they switch to a full growth program of promoter activity including ribosomal and metabolic genes. CONCLUSIONS: The observed exponential increase in enzymes for the specific carbon source followed by an abrupt switch to production of general growth genes is a solution of an optimal control model, known as bang-bang control. The present approach contributes to the understanding of lag phase, the least studied of bacterial growth phases.

Surface growth of a motile bacterial population resembles growth in a chemostat.

The growth behavior in well-mixed bacterial cultures is relatively well understood. However, bacteria often grow in heterogeneous conditions on surfaces where their growth is dependent on spatial position, especially in the case of motile populations. For such populations, the relation between growth, motility and spatial position is unclear. We developed a microscope-based assay for quantifying in situ growth and gene expression in space and time, and we observe these parameters in populations of Escherichia coli swimming in galactose soft agar plates. We find that the bacterial density and the shape of the motile population, after an initial transient, are constant in time. By considering not only the advancing population but also the fraction that lags behind, we propose a growth model that relates spatial distribution, motility and growth rate. This model, that is similar to bacterial growth in a chemostat predicts that the fraction of the population lagging behind is inversely proportional to the velocity of the motile population. We test this prediction by modulating motility using inducible expression of the flagellar sigma factor FliA. Finally, we observe that bacteria in the chemotactic ring express higher relative levels of the chemotaxis and galactose metabolism genes fliC, fliL and galE than those that stay behind in the center of the plate.

Mode of regulation and the insulation of bacterial gene expression.

A gene can be said to be insulated from environmental variations if its expression level depends only on its cognate inducers, and not on variations in conditions. We tested the insulation of the lac promoter of E. coli and of synthetic constructs in which the transcription factor CRP acts as either an activator or a repressor, by measuring their input function-their expression as a function of inducers-in different growth conditions. We find that the promoter activities show sizable variation across conditions of 10%-100% (SD/mean). When the promoter is bound to its cognate regulator(s), variation across conditions is smaller than when it is unbound. Thus, mode of regulation affects insulation: activators seem to show better insulation at high expression levels, and repressors at low expression levels. This may explain the Savageau demand rule, in which E. coli genes needed often in the natural environment tend to be regulated by activators, and rarely needed genes by repressors. The present approach can be used to study insulation in other genes and organisms.

A genome-wide analysis of promoter-mediated phenotypic noise in Escherichia coli.

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as "phenotypic noise." In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alone.

Negative auto-regulation increases the input dynamic-range of the arabinose system of Escherichia coli.

BACKGROUND: Gene regulation networks are made of recurring regulatory patterns, called network motifs. One of the most common network motifs is negative auto-regulation, in which a transcription factor represses its own production. Negative auto-regulation has several potential functions: it can shorten the response time (time to reach halfway to steady-state), stabilize expression against noise, and linearize the gene's input-output response curve. This latter function of negative auto-regulation, which increases the range of input signals over which downstream genes respond, has been studied by theory and synthetic gene circuits. Here we ask whether negative auto-regulation preserves this function also in the context of a natural system, where it is embedded within many additional interactions. To address this, we studied the negative auto-regulation motif in the arabinose utilization system of Escherichia coli, in which negative auto-regulation is part of a complex regulatory network. RESULTS: We find that when negative auto-regulation is disrupted by placing the regulator araC under constitutive expression, the input dynamic range of the arabinose system is reduced by 10-fold. The apparent Hill coefficient of the induction curve changes from about n = 1 with negative auto-regulation, to about n = 2 when it is disrupted. We present a mathematical model that describes how negative auto-regulation can increase input dynamic-range, by coupling the transcription factor protein level to the input signal. CONCLUSIONS: Here we demonstrate that the negative auto-regulation motif in the native arabinose system of Escherichia coli increases the range of arabinose signals over which the system can respond. In this way, negative auto-regulation may help to increase the input dynamic-range while maintaining the specificity of cooperative regulatory systems. This function may contribute to explaining the common occurrence of negative auto-regulation in biological systems.

Robust control of nitrogen assimilation by a bifunctional enzyme in E. coli.

Bacteria regulate the assimilation of multiple nutrients to enable growth. How is balanced utilization achieved, despite fluctuations in the concentrations of the enzymes that make up the regulatory circuitry? Here we address this question by studying the nitrogen system of E. coli. A mechanism based on the avidity of a bifunctional enzyme, adenylyltransferase (AT/AR), to its multimeric substrate, glutamine synthetase, is proposed to maintain a robust ratio between two key metabolites, glutamine and α-ketoglutarate. This ratio is predicted to be insensitive to variations in protein levels of the core circuit and to the rate of nitrogen utilization. We find using mass spectrometry that the metabolite ratio is robust to variations in protein levels and that this robustness depends on the bifunctional enzyme. Moreover, robustness carries through to the bacteria growth rate. Interrupting avidity by adding a monofunctional AT/AR mutant to the native system abolishes robustness, as predicted by the proposed mechanism.

Invariant distribution of promoter activities in Escherichia coli.

Cells need to allocate their limited resources to express a wide range of genes. To understand how Escherichia coli partitions its transcriptional resources between its different promoters, we employ a robotic assay using a comprehensive reporter strain library for E. coli to measure promoter activity on a genomic scale at high-temporal resolution and accuracy. This allows continuous tracking of promoter activity as cells change their growth rate from exponential to stationary phase in different media. We find a heavy-tailed distribution of promoter activities, with promoter activities spanning several orders of magnitude. While the shape of the distribution is almost completely independent of the growth conditions, the identity of the promoters expressed at different levels does depend on them. Translation machinery genes, however, keep the same relative expression levels in the distribution across conditions, and their fractional promoter activity tracks growth rate tightly. We present a simple optimization model for resource allocation which suggests that the observed invariant distributions might maximize growth rate. These invariant features of the distribution of promoter activities may suggest design constraints that shape the allocation of transcriptional resources.

The incoherent feed-forward loop can generate non-monotonic input functions for genes.

Gene regulation networks contain recurring circuit patterns called network motifs. One of the most common network motif is the incoherent type 1 feed-forward loop (I1-FFL), in which an activator controls both gene and repressor of that gene. This motif was shown to act as a pulse generator and response accelerator of gene expression. Here we consider an additional function of this motif: the I1-FFL can generate a non-monotonic dependence of gene expression on the input signal. Here, we study this experimentally in the galactose system of Escherichia coli, which is regulated by an I1-FFL. The promoter activity of two of the gal operons, galETK and galP, peaks at intermediate levels of the signal cAMP. We find that mutants in which the I1-FFL is disrupted lose this non-monotonic behavior, and instead display monotonic input functions. Theoretical analysis suggests that non-monotonic input functions can be achieved for a wide range of parameters by the I1-FFL. The models also suggest regimes where a monotonic input-function can occur, as observed in the mglBAC operon regulated by the same I1-FFL. The present study thus experimentally demonstrates how upstream circuitry can affect gene input functions and how an I1-FFL functions within its natural context in the cell.