Transcriptional drift in aging cells: A global decontroller.

As cells age, they undergo a remarkable global change: In transcriptional drift, hundreds of genes become overexpressed while hundreds of others become underexpressed. Using archetype modeling and Gene Ontology analysis on data from aging Caenorhabditis elegans worms, we find that the up-regulated genes code for sensory proteins upstream of stress responses and down-regulated genes are growth- and metabolism-related. We observe similar trends within human fibroblasts, suggesting that this process is conserved in higher organisms. We propose a simple mechanistic model for how such global coordination of multiprotein expression levels may be achieved by the binding of a single factor that concentrates with age in C. elegans. A key implication is that a cell's own responses are part of its aging process, so unlike wear-and-tear processes, intervention might be able to modulate these effects.

Achieving Occam's razor: Deep learning for optimal model reduction.

All fields of science depend on mathematical models. Occam's razor refers to the principle that good models should exclude parameters beyond those minimally required to describe the systems they represent. This is because redundancy can lead to incorrect estimates of model parameters from data, and thus inaccurate or ambiguous conclusions. Here, we show how deep learning can be powerfully leveraged to apply Occam's razor to model parameters. Our method, FixFit, uses a feedforward deep neural network with a bottleneck layer to characterize and predict the behavior of a given model from its input parameters. FixFit has three major benefits. First, it provides a metric to quantify the original model's degree of complexity. Second, it allows for the unique fitting of data. Third, it provides an unbiased way to discriminate between experimental hypotheses that add value versus those that do not. In three use cases, we demonstrate the broad applicability of this method across scientific domains. To validate the method using a known system, we apply FixFit to recover known composite parameters for the Kepler orbit model and a dynamic model of blood glucose regulation. In the latter, we demonstrate the ability to fit the latent parameters to real data. To illustrate how the method can be applied to less well-established fields, we use it to identify parameters for a multi-scale brain model and reduce the search space for viable candidate mechanisms.

Mapping the Single-cell Differentiation Landscape of Osteosarcoma.

PURPOSE: The genetic intratumoral heterogeneity observed in human osteosarcomas (OS) poses challenges for drug development and the study of cell fate, plasticity, and differentiation, processes linked to tumor grade, cell metastasis, and survival. EXPERIMENTAL DESIGN: To pinpoint errors in OS differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directsMSCs toward adipogenic and osteoblastic fates. Incorporating pre-existing chondrocyte data, we applied trajectory analysis and non-negative matrix factorization (NMF) to generate the first human mesenchymal differentiation atlas. RESULTS: This 'roadmap' served as a reference to delineate the cellular composition of morphologically complex OS tumors and quantify each cell's lineage commitment. Projecting a bulk RNA-seq OS dataset onto this roadmap unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes. CONCLUSIONS: Our study quantifies OS differentiation and lineage, a prerequisite to better understanding lineage-specific differentiation bottlenecks that might someday be targeted therapeutically.

D-ꞵ-hydroxybutyrate stabilizes hippocampal CA3-CA1 circuit during acute insulin resistance.

The brain primarily relies on glycolysis for mitochondrial respiration but switches to alternative fuels such as ketone bodies (KBs) when less glucose is available. Neuronal KB uptake, which does not rely on glucose transporter 4 (GLUT4) or insulin, has shown promising clinical applicability in alleviating the neurological and cognitive effects of disorders with hypometabolic components. However, the specific mechanisms by which such interventions affect neuronal functions are poorly understood. In this study, we pharmacologically blocked GLUT4 to investigate the effects of exogenous KB D-ꞵ-hydroxybutyrate (D-ꞵHb) on mouse brain metabolism during acute insulin resistance (AIR). We found that both AIR and D-ꞵHb had distinct impacts across neuronal compartments: AIR decreased synaptic activity and long-term potentiation (LTP) and impaired axonal conduction, synchronization, and action potential properties, while D-ꞵHb rescued neuronal functions associated with axonal conduction, synchronization, and LTP.

Functional transcriptional signatures for tumor-type-agnostic phenotype prediction.

Cancer transcriptional patterns exhibit both shared and unique features across diverse cancer types, but whether these patterns are sufficient to characterize the full breadth of tumor phenotype heterogeneity remains an open question. We hypothesized that cancer transcriptional diversity mirrors patterns in normal tissues optimized for distinct functional tasks. Starting with normal tissue transcriptomic profiles, we use non-negative matrix factorization to derive six distinct transcriptomic phenotypes, called archetypes, which combine to describe both normal tissue patterns and variations across a broad spectrum of malignancies. We show that differential enrichment of these signatures correlates with key tumor characteristics, including overall patient survival and drug sensitivity, independent of clinically actionable DNA alterations. Additionally, we show that in HR+/HER2-breast cancers, metastatic tumors adopt transcriptomic signatures consistent with the invaded tissue. Broadly, our findings suggest that cancer often arrogates normal tissue transcriptomic characteristics as a component of both malignant progression and drug response. This quantitative framework provides a strategy for connecting the diversity of cancer phenotypes and could potentially help manage individual patients.

D-ß-hydroxybutyrate stabilizes hippocampal CA3-CA1 circuit during acute insulin resistance.

1.The brain primarily relies on glycolysis for mitochondrial respiration but switches to alternative fuels such as ketone bodies (KBs) when less glucose is available. Neuronal KB uptake, which does not rely on glucose transporter 4 (GLUT4) or insulin, has shown promising clinical applicability in alleviating the neurological and cognitive effects of disorders with hypometabolic components. However, the specific mechanisms by which such interventions affect neuronal functions are poorly understood. In this study, we pharmacologically blocked GLUT4 to investigate the effects of exogenous KB D-P-hydroxybutyrate (D-ßHb) on mouse brain metabolism during acute insulin resistance (AIR). We found that both AIR and D-ßHb had distinct impacts across neuronal compartments: AIR decreased synaptic activity and long-term potentiation (LTP) and impaired axonal conduction, synchronization, and action potential (AP) properties, while D- PHb rescued neuronal functions associated with axonal conduction, synchronization and LTP.

Transcriptional drift in aging cells: A global de-controller.

As cells age, they undergo a remarkable global change: In transcriptional drift, hundreds of genes become overexpressed while hundreds of others become underexpressed. Using archetype modeling and Gene Ontology analysis on data from aging Caenorhabditis elegans worms, we find that the upregulated genes code for sensory proteins upstream of stress responses and downregulated genes are growth- and metabolism-related. We propose a simple mechanistic model for how such global coordination of multi-protein expression levels may be achieved by the binding of a single ligand that concentrates with age. A key implication is that a cell's own responses are part of its aging process, so unlike for wear-and-tear processes, intervention might be able to modulate these effects.

Mapping the Single-cell Differentiation Landscape of Osteosarcoma.

The genetic and intratumoral heterogeneity observed in human osteosarcomas (OS) poses challenges for drug development and the study of cell fate, plasticity, and differentiation, processes linked to tumor grade, cell metastasis, and survival. To pinpoint errors in OS differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directs MSCs toward adipogenic and osteoblastic fates. Incorporating pre-existing chondrocyte data, we applied trajectory analysis and non-negative matrix factorization (NMF) to generate the first human mesenchymal differentiation atlas. This 'roadmap' served as a reference to delineate the cellular composition of morphologically complex OS tumors and quantify each cell's lineage commitment. Projecting these signatures onto a bulk RNA-seq OS dataset unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes. Our study takes the critical first step in accurately quantifying OS differentiation and lineage, a prerequisite to better understanding global differentiation bottlenecks that might someday be targeted therapeutically.

Ion gradient-driven bifurcations of a multi-scale neuronal model.

Metabolic limitations within the brain frequently arise in the context of aging and disease. As the largest consumers of energy within the brain, ion pumps that maintain the neuronal membrane potential are the most affected when energy supply becomes limited. To characterize the effects of such limitations, we analyze the ion gradients present in a conductance-based (Morris-Lecar) neural mass model. We show the existence and locations of Neimark-Sacker and period-doubling bifurcations in the sodium, calcium, and potassium reversal potentials and demonstrate that these bifurcations form physiologically relevant bounds of ion gradient variability. Within these bounds, we show how depolarization of the gradients causes decreased neural activity. We also show that the depolarization of ion gradients decreases inter-regional coherence, causing a shift in the critical point at which the coupling occurs and thereby inducing loss of synchrony between regions. In this way, we show that the Larter-Breakspear model captures ion gradient variability present at the microscale level and propagates these changes to the macroscale effects such as those observed in human neuroimaging studies.

The maximum entropy principle for compositional data.

BACKGROUND: Compositional systems, represented as parts of some whole, are ubiquitous. They encompass the abundances of proteins in a cell, the distribution of organisms in nature, and the stoichiometry of the most basic chemical reactions. Thus, a central goal is to understand how such processes emerge from the behaviors of their components and their pairwise interactions. Such a study, however, is challenging for two key reasons. Firstly, such systems are complex and depend, often stochastically, on their constituent parts. Secondly, the data lie on a simplex which influences their correlations. RESULTS: To resolve both of these issues, we provide a general and data-driven modeling tool for compositional systems called Compositional Maximum Entropy (CME). By integrating the prior geometric structure of compositions with sample-specific information, CME infers the underlying multivariate relationships between the constituent components. We provide two proofs of principle. First, we measure the relative abundances of different bacteria and infer how they interact. Second, we show that our method outperforms a common alternative for the extraction of gene-gene interactions in triple-negative breast cancer. CONCLUSIONS: CME provides novel and biologically-intuitive insights and is promising as a comprehensive quantitative framework for compositional data.

Metabolism modulates network synchrony in the aging brain.

Brain aging is associated with hypometabolism and global changes in functional connectivity. Using functional MRI (fMRI), we show that network synchrony, a collective property of brain activity, decreases with age. Applying quantitative methods from statistical physics, we provide a generative (Ising) model for these changes as a function of the average communication strength between brain regions. We find that older brains are closer to a critical point of this communication strength, in which even small changes in metabolism lead to abrupt changes in network synchrony. Finally, by experimentally modulating metabolic activity in younger adults, we show how metabolism alone-independent of other changes associated with aging-can provide a plausible candidate mechanism for marked reorganization of brain network topology.

The Refractory Period Matters: Unifying Mechanisms of Macroscopic Brain Waves.

The relationship between complex brain oscillations and the dynamics of individual neurons is poorly understood. Here we utilize maximum caliber, a dynamical inference principle, to build a minimal yet general model of the collective (mean field) dynamics of large populations of neurons. In agreement with previous experimental observations, we describe a simple, testable mechanism, involving only a single type of neuron, by which many of these complex oscillatory patterns may emerge. Our model predicts that the refractory period of neurons, which has often been neglected, is essential for these behaviors.

Inferring a network from dynamical signals at its nodes.

We give an approximate solution to the difficult inverse problem of inferring the topology of an unknown network from given time-dependent signals at the nodes. For example, we measure signals from individual neurons in the brain, and infer how they are inter-connected. We use Maximum Caliber as an inference principle. The combinatorial challenge of high-dimensional data is handled using two different approximations to the pairwise couplings. We show two proofs of principle: in a nonlinear genetic toggle switch circuit, and in a toy neural network.

Diet modulates brain network stability, a biomarker for brain aging, in young adults.

Epidemiological studies suggest that insulin resistance accelerates progression of age-based cognitive impairment, which neuroimaging has linked to brain glucose hypometabolism. As cellular inputs, ketones increase Gibbs free energy change for ATP by 27% compared to glucose. Here we test whether dietary changes are capable of modulating sustained functional communication between brain regions (network stability) by changing their predominant dietary fuel from glucose to ketones. We first established network stability as a biomarker for brain aging using two large-scale (n = 292, ages 20 to 85 y; n = 636, ages 18 to 88 y) 3 T functional MRI (fMRI) datasets. To determine whether diet can influence brain network stability, we additionally scanned 42 adults, age < 50 y, using ultrahigh-field (7 T) ultrafast (802 ms) fMRI optimized for single-participant-level detection sensitivity. One cohort was scanned under standard diet, overnight fasting, and ketogenic diet conditions. To isolate the impact of fuel type, an independent overnight fasted cohort was scanned before and after administration of a calorie-matched glucose and exogenous ketone ester (d-ß-hydroxybutyrate) bolus. Across the life span, brain network destabilization correlated with decreased brain activity and cognitive acuity. Effects emerged at 47 y, with the most rapid degeneration occurring at 60 y. Networks were destabilized by glucose and stabilized by ketones, irrespective of whether ketosis was achieved with a ketogenic diet or exogenous ketone ester. Together, our results suggest that brain network destabilization may reflect early signs of hypometabolism, associated with dementia. Dietary interventions resulting in ketone utilization increase available energy and thus may show potential in protecting the aging brain.

Perspective: Maximum caliber is a general variational principle for dynamical systems.

We review here Maximum Caliber (Max Cal), a general variational principle for inferring distributions of paths in dynamical processes and networks. Max Cal is to dynamical trajectories what the principle of maximum entropy is to equilibrium states or stationary populations. In Max Cal, you maximize a path entropy over all possible pathways, subject to dynamical constraints, in order to predict relative path weights. Many well-known relationships of non-equilibrium statistical physics-such as the Green-Kubo fluctuation-dissipation relations, Onsager's reciprocal relations, and Prigogine's minimum entropy production-are limited to near-equilibrium processes. Max Cal is more general. While it can readily derive these results under those limits, Max Cal is also applicable far from equilibrium. We give examples of Max Cal as a method of inference about trajectory distributions from limited data, finding reaction coordinates in bio-molecular simulations, and modeling the complex dynamics of non-thermal systems such as gene regulatory networks or the collective firing of neurons. We also survey its basis in principle and some limitations.

Hydrodynamic Hunters.

The Gram-negative Bdellovibrio bacteriovorus (BV) is a model bacterial predator that hunts other bacteria and may serve as a living antibiotic. Despite over 50 years since its discovery, it is suggested that BV probably collides into its prey at random. It remains unclear to what degree, if any, BV uses chemical cues to target its prey. The targeted search problem by the predator for its prey in three dimensions is a difficult problem: it requires the predator to sensitively detect prey and forecast its mobile prey's future position on the basis of previously detected signal. Here instead we find that rather than chemically detecting prey, hydrodynamics forces BV into regions high in prey density, thereby improving its odds of a chance collision with prey and ultimately reducing BV's search space for prey. We do so by showing that BV's dynamics are strongly influenced by self-generated hydrodynamic flow fields forcing BV onto surfaces and, for large enough defects on surfaces, forcing BV in orbital motion around these defects. Key experimental controls and calculations recapitulate the hydrodynamic origin of these behaviors. While BV's prey (Escherichia coli) are too small to trap BV in hydrodynamic orbit, the prey are also susceptible to their own hydrodynamic fields, substantially confining them to surfaces and defects where mobile predator and prey density is now dramatically enhanced. Colocalization, driven by hydrodynamics, ultimately reduces BV's search space for prey from three to two dimensions (on surfaces) even down to a single dimension (around defects). We conclude that BV's search for individual prey remains random, as suggested in the literature, but confined, however-by generic hydrodynamic forces-to reduced dimensionality.