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.

Geroprotective interventions converge on gene expression programs of reduced inflammation and restored fatty acid metabolism.

Understanding the mechanisms of geroprotective interventions is central to aging research. We compare four prominent interventions: senolysis, caloric restriction, in vivo partial reprogramming, and heterochronic parabiosis. Using published mice transcriptomic data, we juxtapose these interventions against normal aging. We find a gene expression program common to all four interventions, in which inflammation is reduced and several metabolic processes, especially fatty acid metabolism, are increased. Normal aging exhibits the inverse of this signature across multiple organs and tissues. A similar inverse signature arises in three chronic inflammation disease models in a non-aging context, suggesting that the shift in metabolism occurs downstream of inflammation. Chronic inflammation is also shown to accelerate transcriptomic age. We conclude that a core mechanism of geroprotective interventions acts through the reduction of inflammation with downstream effects that restore fatty acid metabolism. This supports the notion of directly targeting genes associated with these pathways to mitigate age-related deterioration.

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.

Major depressive disorder and bistability in an HPA-CNS toggle switch.

Major depressive disorder (MDD) is the most common psychiatric disorder. It has a complex and heterogeneous etiology. Most treatments take weeks to show effects and work well only for a fraction of the patients. Thus, new concepts are needed to understand MDD and its dynamics. One of the strong correlates of MDD is increased activity and dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis which produces the stress hormone cortisol. Existing mathematical models of the HPA axis describe its operation on the scale of hours, and thus are unable to explore the dynamic on the scale of weeks that characterizes many aspects of MDD. Here, we propose a mathematical model of MDD on the scale of weeks, a timescale provided by the growth of the HPA hormone glands under control of HPA hormones. We add to this the mutual inhibition of the HPA axis and the hippocampus and other regions of the central nervous system (CNS) that forms a toggle switch. The model shows bistability between euthymic and depressed states, with a slow timescale of weeks in its dynamics. It explains why prolonged but not acute stress can trigger a self-sustaining depressive episode that persists even after the stress is removed. The model explains the weeks timescale for drugs to take effect, as well as the dysregulation of the HPA axis in MDD, based on gland mass changes. This understanding of MDD dynamics may help to guide strategies for treatment.

Excitable dynamics of flares and relapses in autoimmune diseases.

Many autoimmune disorders exhibit flares in which symptoms erupt and then decline, as exemplified by multiple sclerosis (MS) in its relapsing-remitting form. Existing mathematical models of autoimmune flares often assume regular oscillations, failing to capture the stochastic and non-periodic nature of flare-ups. We suggest that autoimmune flares are driven by excitable dynamics triggered by stochastic events auch as stress, infection and other factors. Our minimal model, involving autoreactive and regulatory T-cells, demonstrates this concept. Autoimmune response initiates antigen-induced expansion through positive feedback, while regulatory cells counter the autoreactive cells through negative feedback. The model explains the decrease in MS relapses during pregnancy and the subsequent surge postpartum, based on lymphocyte dynamics. Additionally, it identifies potential therapeutic targets, predicting significant reduction in relapse rate from mild adjustments of regulatory T cell activity or production. These findings indicate that excitable dynamics may underlie flare-ups across various autoimmune disorders, potentially informing treatment strategies.

The tumor microenvironment shows a hierarchy of cell-cell interactions dominated by fibroblasts.

The tumor microenvironment (TME) is comprised of non-malignant cells that interact with each other and with cancer cells, critically impacting cancer biology. The TME is complex, and understanding it requires simplifying approaches. Here we provide an experimental-mathematical approach to decompose the TME into small circuits of interacting cell types. We find, using female breast cancer single-cell-RNA-sequencing data, a hierarchical network of interactions, with cancer-associated fibroblasts (CAFs) at the top secreting factors primarily to tumor-associated macrophages (TAMs). This network is composed of repeating circuit motifs. We isolate the strongest two-cell circuit motif by culturing fibroblasts and macrophages in-vitro, and analyze their dynamics and transcriptomes. This isolated circuit recapitulates the hierarchy of in-vivo interactions, and enables testing the effect of ligand-receptor interactions on cell dynamics and function, as we demonstrate by identifying a mediator of CAF-TAM interactions - RARRES2, and its receptor CMKLR1. Thus, the complexity of the TME may be simplified by identifying small circuits, facilitating the development of strategies to modulate the TME.

Damage dynamics and the role of chance in the timing of E. coli cell death.

Genetically identical cells in the same stressful condition die at different times. The origin of this stochasticity is unclear; it may arise from different initial conditions that affect the time of demise, or from a stochastic damage accumulation mechanism that erases the initial conditions and instead amplifies noise to generate different lifespans. To address this requires measuring damage dynamics in individual cells over the lifespan, but this has rarely been achieved. Here, we used a microfluidic device to measure membrane damage in 635 carbon-starved Escherichia coli cells at high temporal resolution. We find that initial conditions of damage, size or cell-cycle phase do not explain most of the lifespan variation. Instead, the data points to a stochastic mechanism in which noise is amplified by a rising production of damage that saturates its own removal. Surprisingly, the relative variation in damage drops with age: cells become more similar to each other in terms of relative damage, indicating increasing determinism with age. Thus, chance erases initial conditions and then gives way to increasingly deterministic dynamics that dominate the lifespan distribution.

Autoimmune thyroid diseases as a cost of physiological autoimmune surveillance.

Graves' disease (GD) and Hashimoto's thyroiditis (HT) are common autoimmune diseases of the thyroid gland, causing hyperthyroidism and hypothyroidism, respectively. Despite their opposing clinical manifestation, they have several enigmatic links. Here, we propose that GD and HT have the same fundamental origin: both diseases are the cost of a beneficial physiological process called autoimmune surveillance of hypersecreting mutants. Autoreactive T cells selectively eliminate mutant cells that hypersecrete the hormones and threaten to become toxic nodules. These T cells can trigger a humoral response in susceptible individuals, leading to the production of antibodies against thyroid antigens. This shared origin can explain similarities in incidence and risk factors between HT and GD, despite their opposite clinical phenotypes.

Sex-dependent control of pheromones on social organization within groups of wild house mice.

Dominance hierarchy is a fundamental social phenomenon in a wide range of mammalian species, critically affecting fitness and health. Here, we investigate the role of pheromone signals in the control of social hierarchies and individual personalities within groups of wild mice. For this purpose, we combine high-throughput behavioral phenotyping with computational tools in freely interacting groups of wild house mice, males and females, in an automated, semi-natural system. We show that wild mice form dominance hierarchies in both sexes but use sex-specific strategies, displaying distinct male-typical and female-typical behavioral personalities that were also associated with social ranking. Genetic disabling of VNO-mediated pheromone detection generated opposite behavioral effects within groups, enhancing social interactions in males and reducing them in females. Behavioral personalities in the mutated mice displayed mixtures of male-typical and female-typical behaviors, thus blurring sex differences. In addition, rank-associated personalities were abolished despite the fact that both sexes of mutant mice formed stable hierarchies. These findings suggest that group organization is governed by pheromone-mediated sex-specific neural circuits and pave the way to investigate the mechanisms underlying sexual dimorphism in dominance hierarchies under naturalistic settings.

Rules for body fat interventions based on an operating point mechanism.

Interventions to reduce fat are important for human health. However, they can have opposing effects such as exercise that decreases fat but increases food intake, or coherent effects such as leptin resistance which raises both. Furthermore, some interventions show an overshoot in food intake, such as recovery from a diet, whereas others do not. To explain these properties we present a graphical framework called the operating point model, based on leptin control of feeding behavior. Steady-state fat and food intake is given by the intersection of two experimental curves - steady-state fat at a given food intake and ad libitum food intake at a given fat level. Depending on which curve an intervention shifts, it has opposing or coherent effects with or without overshoot, in excellent agreement with rodent data. The model also explains the quadratic relation between leptin and fat in humans. These concepts may guide the understanding of fat regulation disorders.

An autocrine signaling circuit in hepatic stellate cells underlies advanced fibrosis in nonalcoholic steatohepatitis.

Advanced hepatic fibrosis, driven by the activation of hepatic stellate cells (HSCs), affects millions worldwide and is the strongest predictor of mortality in nonalcoholic steatohepatitis (NASH); however, there are no approved antifibrotic therapies. To identify antifibrotic drug targets, we integrated progressive transcriptomic and morphological responses that accompany HSC activation in advanced disease using single-nucleus RNA sequencing and tissue clearing in a robust murine NASH model. In advanced fibrosis, we found that an autocrine HSC signaling circuit emerged that was composed of 68 receptor-ligand interactions conserved between murine and human NASH. These predicted interactions were supported by the parallel appearance of markedly increased direct stellate cell-cell contacts in murine NASH. As proof of principle, pharmacological inhibition of one such autocrine interaction, neurotrophic receptor tyrosine kinase 3-neurotrophin 3, inhibited human HSC activation in culture and reversed advanced murine NASH fibrosis. In summary, we uncovered a repertoire of antifibrotic drug targets underlying advanced fibrosis in vivo. The findings suggest a therapeutic paradigm in which stage-specific therapies could yield enhanced antifibrotic efficacy in patients with advanced hepatic fibrosis.

Prostanoid receptors in hyperfiltration-mediated glomerular injury: Novel agonists and antagonists reveal opposing roles for EP2 and EP4 receptors.

Increased fluid-flow shear stress (FFSS) contributes to hyperfiltration-induced podocyte and glomerular injury resulting in progression of chronic kidney disease (CKD). We reported that increased FFSS in vitro and in vivo upregulates PGE2 receptor EP2 (but not EP4 expression), COX2-PGE2 -EP2 axis, and EP2-linked Akt-GSK3ß-ß-catenin signaling pathway in podocytes. To understand and use the disparities between PGE2 receptors, specific agonists, and antagonists of EP2 and EP4 were used to assess phosphorylation of Akt, GSK3ß and ß-catenin in podocytes using Western blotting, glomerular filtration barrier function using in vitro albumin permeability (Palb ) assay, and mitigation of hyperfiltration-induced injury in unilaterally nephrectomized (UNX) mice at 1 and 6 months. Results show an increase in Palb by PGE2 , EP2 agonist (EP2AGO ) and EP4 antagonist (EP4ANT ), but not by EP2 antagonist (EP2ANT ) or EP4 agonist (EP4AGO ). Pretreatment with EP2ANT blocked the effect of PGE2 or EP2AGO on Palb . Modulation of EP2 and EP4 also induced opposite effects on phosphorylation of Akt and ß-Catenin. Individual agonists or antagonists of EP2 or EP4 did not induce significant improvement in albuminuria in UNX mice. However, treatment with a combination EP2ANT + EP4AGO for 1 or 6 months caused a robust decrease in albuminuria. EP2ANT + EP4AGO combination did not impact adaptive hypertrophy or increased serum creatinine. Observed differences between expression of EP2 and EP4 on the glomerular barrier highlight these receptors as potential targets for intervention. Safe and effective mitigating effect of EP2ANT + EP4AGO presents a novel opportunity to delay the progression of hyperfiltration-associated CKD as seen in transplant donors.

Dynamics of thyroid diseases and thyroid-axis gland masses.

Thyroid disorders are common and often require lifelong hormone replacement. Treating thyroid disorders involves a fascinating and troublesome delay, in which it takes many weeks for serum thyroid-stimulating hormone (TSH) concentration to normalize after thyroid hormones return to normal. This delay challenges attempts to stabilize thyroid hormones in millions of patients. Despite its importance, the physiological mechanism for the delay is unclear. Here, we present data on hormone delays from Israeli medical records spanning 46 million life-years and develop a mathematical model for dynamic compensation in the thyroid axis, which explains the delays. The delays are due to a feedback mechanism in which peripheral thyroid hormones and TSH control the growth of the thyroid and pituitary glands; enlarged or atrophied glands take many weeks to recover upon treatment due to the slow turnover of the tissues. The model explains why thyroid disorders such as Hashimoto's thyroiditis and Graves' disease have both subclinical and clinical states and explains the complex inverse relation between TSH and thyroid hormones. The present model may guide approaches to dynamically adjust the treatment of thyroid disorders.

A Mechanism for Ovulation Number Control.

Every menstrual cycle, many follicles begin to develop but only a specific number ovulate. This ovulation number determines how many offspring are produced per litter, and differs between species. The physiological mechanism that controls ovulation number is unknown; a class of mathematical models can explain it, but these models have no physiological basis. Here, we suggest a physiological mechanism for ovulation number control, which enables selection of a specific number of follicles out of many, and analyze it in a mathematical model of follicular growth. The mechanism is based on a signal, intra-follicular androgen concentration, that measures follicle size relative to the other follicles. This signal has a biphasic effect, suppressing follicles that are too large or too small compared to others. The ovulation number is determined by the androgen inhibitory thresholds. The model has a scaling symmetry that explains why the dominant follicles grow linearly with time, as observed in human ultrasound data. This approach also explains how chronic hyperandrogenism disrupts ovulation in polycystic ovary syndrome (PCOS), a leading cause of infertility. We propose specific experiments for testing the proposed mechanism.

The dopamine circuit as a reward-taxis navigation system.

Studying the brain circuits that control behavior is challenging, since in addition to their structural complexity there are continuous feedback interactions between actions and sensed inputs from the environment. It is therefore important to identify mathematical principles that can be used to develop testable hypotheses. In this study, we use ideas and concepts from systems biology to study the dopamine system, which controls learning, motivation, and movement. Using data from neuronal recordings in behavioral experiments, we developed a mathematical model for dopamine responses and the effect of dopamine on movement. We show that the dopamine system shares core functional analogies with bacterial chemotaxis. Just as chemotaxis robustly climbs chemical attractant gradients, the dopamine circuit performs 'reward-taxis' where the attractant is the expected value of reward. The reward-taxis mechanism provides a simple explanation for scale-invariant dopaminergic responses and for matching in free operant settings, and makes testable quantitative predictions. We propose that reward-taxis is a simple and robust navigation strategy that complements other, more goal-directed navigation mechanisms.

Strategies and foundations for scientific discovery in longitudinal studies of bipolar disorder.

Bipolar disorder (BD) is a complex and dynamic condition with a typical onset in late adolescence or early adulthood followed by an episodic course with intervening periods of subthreshold symptoms or euthymia. It is complicated by the accumulation of comorbid medical and psychiatric disorders. The etiology of BD remains unknown and no reliable biological markers have yet been identified. This is likely due to lack of comprehensive ontological framework and, most importantly, the fact that most studies have been based on small nonrepresentative clinical samples with cross-sectional designs. We propose to establish large, global longitudinal cohorts of BD studied consistently in a multidimensional and multidisciplinary manner to determine etiology and help improve treatment. Herein we propose collection of a broad range of data that reflect the heterogenic phenotypic manifestations of BD that include dimensional and categorical measures of mood, neurocognitive, personality, behavior, sleep and circadian, life-story, and outcomes domains. In combination with genetic and biological information such an approach promotes the integrating and harmonizing of data within and across current ontology systems while supporting a paradigm shift that will facilitate discovery and become the basis for novel hypotheses.

Controls for Phylogeny and Robust Analysis in Pareto Task Inference.

Understanding the tradeoffs faced by organisms is a major goal of evolutionary biology. One of the main approaches for identifying these tradeoffs is Pareto task inference (ParTI). Two recent papers claim that results obtained in ParTI studies are spurious due to phylogenetic dependence (Mikami T, Iwasaki W. 2021. The flipping t-ratio test: phylogenetically informed assessment of the Pareto theory for phenotypic evolution. Methods Ecol Evol. 12(4):696-706) or hypothetical p-hacking and population-structure concerns (Sun M, Zhang J. 2021. Rampant false detection of adaptive phenotypic optimization by ParTI-based Pareto front inference. Mol Biol Evol. 38(4):1653-1664). Here, we show that these claims are baseless. We present a new method to control for phylogenetic dependence, called SibSwap, and show that published ParTI inference is robust to phylogenetic dependence. We show how researchers avoided p-hacking by testing for the robustness of preprocessing choices. We also provide new methods to control for population structure and detail the experimental tests of ParTI in systems ranging from ammonites to cancer gene expression. The methods presented here may help to improve future ParTI studies.

Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Temporal fluctuations in chemotaxis gain implement a simulated-tempering strategy for efficient navigation in complex environments.

Bacterial chemotaxis is a major testing ground for systems biology, including the role of fluctuations and individual variation. Individual bacteria vary in their tumbling frequency and adaptation time. Recently, large cell-cell variation was also discovered in chemotaxis gain, which determines the sensitivity of the tumbling rate to attractant gradients. Variation in gain is puzzling, because low gain impairs chemotactic velocity. Here, we provide a functional explanation for gain variation by establishing a formal analogy between chemotaxis and algorithms for sampling probability distributions. We show that temporal fluctuations in gain implement simulated tempering, which allows sampling of attractant distributions with many local peaks. Periods of high gain allow bacteria to detect and climb gradients quickly, and periods of low gain allow them to move to new peaks. Gain fluctuations thus allow bacteria to thrive in complex environments, and more generally they may play an important functional role for organism navigation.