Cortical Circuit Dynamics Are Homeostatically Tuned to Criticality In Vivo

Psychedelics and disorders of consciousness: the current landscape and the path forward

Modern medicine has been shaken by the surge of psychedelic science that proposes a new approach to mitigate mental disorders, such as depression and post-traumatic stress disorder. Clinical trials to investigate whether psychedelic substances can treat psychiatric conditions are now underway, yet less discussion gravitates around their use in neurological disorders due to brain injury. One suggested implementation of brain-complexity enhancing psychedelics is to treat people with post-comatose disorders of consciousness (DoC). In this article, we discuss the rationale of this endeavour, examining possible outcomes of such experiments by postulating the existence of an optimal level of complexity. We consider the possible counterintuitive effects of both psychedelics and DoC on the functional connectivity of the default mode network and its possible impact on selfhood. We also elaborate on the role of computational modelling in providing complementary information to experimental studies, both contributing to our understanding of the treatment mechanisms and providing a path towards personalized medicine. Finally, we update the discourse surrounding the ethical considerations, encompassing clinical and scientific values.

At the onset of active whisking, the input layer of barrel cortex exhibits a 24 h window of increased excitability that depends on prior experience

The development of motor control over sensory organs is a critical milestone in sensory processing, enabling active exploration and shaping of the sensory environment. However, whether the onset of sensory organ motor control directly influences the development of corresponding sensory cortices remains unknown. Here, we exploit the late onset of whisking behavior in mice to address this question in the somatosensory system. Using ex vivo electrophysiology, we discovered a transient increase in the intrinsic excitability of excitatory neurons in layer IV of the barrel cortex, which processes whisker input, precisely coinciding with the onset of active whisking at postnatal day 14 (P14). This increase in neuronal gain was specific to layer IV, independent of changes in synaptic strength, and required prior sensory experience. Strikingly, the effect was not observed in layer II/III of the barrel cortex or in the visual cortex upon eye opening, suggesting a unique interaction between the development of active sensing and the thalamocortical input layer in the somatosensory system. Predictive modeling indicated that changes in active membrane conductances alone could reliably distinguish P14 neurons in control but not whisker-deprived hemispheres. Our findings demonstrate an experience-dependent, lamina-specific refinement of neuronal excitability tightly linked to the emergence of active whisking. This transient increase in the gain of the thalamic input layer coincides with a critical period for synaptic plasticity in downstream layers, suggesting a role in facilitating cortical maturation and sensory processing. Together, our results provide evidence for a direct interaction between the development of motor control and sensory cortex, offering new insights into the experience-dependent development and refinement of sensory systems. These findings have broad implications for understanding the interplay between motor and sensory development, and how the mechanisms of perception cooperate with behavior.

Modular Arrangement of Synaptic and Intrinsic Homeostatic Plasticity within Visual Cortical Circuits

Neocortical circuits use synaptic and intrinsic forms of homeostatic plasticity to stabilize key features of network activity, but whether these different homeostatic mechanisms act redundantly, or can be independently recruited to stabilize different network features, is unknown. Here we used pharmacological and genetic perturbations both in vitro and in vivo to determine whether synaptic scaling and intrinsic homeostatic plasticity (IHP) are arranged and recruited in a hierarchical or modular manner within L2/3 pyramidal neurons in rodent V1. Surprisingly, although the expression of synaptic scaling and IHP was dependent on overlapping trafficking pathways, they could be independently recruited by manipulating spiking activity or NMDAR signaling, respectively. Further, we found that changes in visual experience that affect NMDAR activation but not mean firing selectively trigger IHP, without recruiting synaptic scaling. These findings support a modular model in which synaptic and intrinsic homeostatic plasticity respond to and stabilize distinct aspects of network activity.

Low-dimensional criticality embedded in high-dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. We offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality-a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescales agrees with established theory of critical dynamics. Our results suggest that the cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Neural criticality from effective latent variables

Observations of power laws in neural activity data have raised the intriguing notion that brains may operate in a critical state. One example of this critical state is 'avalanche criticality', which has been observed in various systems, including cultured neurons, zebrafish, rodent cortex, and human EEG. More recently, power laws were also observed in neural populations in the mouse under an activity coarse-graining procedure, and they were explained as a consequence of the neural activity being coupled to multiple latent dynamical variables. An intriguing possibility is that avalanche criticality emerges due to a similar mechanism. Here, we determine the conditions under which latent dynamical variables give rise to avalanche criticality. We find that populations coupled to multiple latent variables produce critical behavior across a broader parameter range than those coupled to a single, quasi-static latent variable, but in both cases, avalanche criticality is observed without fine-tuning of model parameters. We identify two regimes of avalanches, both critical but differing in the amount of information carried about the latent variable. Our results suggest that avalanche criticality arises in neural systems in which activity is effectively modeled as a population driven by a few dynamical variables and these variables can be inferred from the population activity.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

The recovery of parabolic avalanches in spatially subsampled neuronal networks at criticality

Scaling relationships characterize complex systems at criticality. In the brain, these relationships are evident in scale-invariant activity cascades, so-called neuronal avalanches, quantified by power laws in avalanche size and duration. At the cellular level, neuronal avalanches are identified in spatially distributed groups of neurons that participate in cascades of coincident action potential firing. Such spatiotemporal synchronization is central to theories on brain function, yet scaling relationships in avalanche synchronization have been challenging to study when only a fraction of neurons is observed, underestimating avalanche properties. Here, we study these biases from fractional sampling in an all-to-all, balanced network of excitatory and inhibitory neurons with critical branching process dynamics. We focus on the growth of mean avalanche size with avalanche duration. For parabolic avalanches, this growth is quadratic, quantified by the scaling exponent, χ = 2 , which signifies rapid spatial expansion of coincident firing within a relatively short period of time. In contrast, χ < < 2 for fractionally sampled networks. We show that temporal coarse-graining combined with a threshold for the minimally required coincident firing in the network recovers χ = 2 , even when sampling as few as 0.1% of the neurons. In contrast, a commonly proposed 'crackling noise' approach fails to recover χ under those conditions. Our approach robustly identifies χ = 2 for ongoing neuronal activity in frontal cortex of awake mice using cellular 2-photon imaging. Our findings demonstrate how to correct scaling bias from fractional sampling and identifies rapid, scale-invariant synchronization of cell assemblies in the brain.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory. Our study based on extensive bifurcation analyses thus reveal the functional optimality and criticality of structural and physiological parameters in establishing the baseline operating state of local cortical networks.

Sleep restores an optimal computational regime in cortical networks

Sleep is assumed to subserve homeostatic processes in the brain; however, the set point around which sleep tunes circuit computations is unknown. Slow-wave activity (SWA) is commonly used to reflect the homeostatic aspect of sleep; although it can indicate sleep pressure, it does not explain why animals need sleep. This study aimed to assess whether criticality may be the computational set point of sleep. By recording cortical neuron activity continuously for 10-14 d in freely behaving rats, we show that normal waking experience progressively disrupts criticality and that sleep functions to restore critical dynamics. Criticality is perturbed in a context-dependent manner, and waking experience is causal in driving these effects. The degree of deviation from criticality predicts future sleep/wake behavior more accurately than SWA, behavioral history or other neural measures. Our results demonstrate that perturbation and recovery of criticality is a network homeostatic mechanism consistent with the core, restorative function of sleep.

Excitation-Inhibition Balance, Neural Criticality, and Activities in Neuronal Circuits

Neural activities in local circuits exhibit complex and multilevel dynamic features. Individual neurons spike irregularly, which is believed to originate from receiving balanced amounts of excitatory and inhibitory inputs, known as the excitation-inhibition balance. The spatial-temporal cascades of clustered neuronal spikes occur in variable sizes and durations, manifested as neural avalanches with scale-free features. These may be explained by the neural criticality hypothesis, which posits that neural systems operate around the transition between distinct dynamic states. Here, we summarize the experimental evidence for and the underlying theory of excitation-inhibition balance and neural criticality. Furthermore, we review recent studies of excitatory-inhibitory networks with synaptic kinetics as a simple solution to reconcile these two apparently distinct theories in a single circuit model. This provides a more unified understanding of multilevel neural activities in local circuits, from spontaneous to stimulus-response dynamics.

Neural Geometrodynamics, Complexity, and Plasticity: A Psychedelics Perspective

We explore the intersection of neural dynamics and the effects of psychedelics in light of distinct timescales in a framework integrating concepts from dynamics, complexity, and plasticity. We call this framework neural geometrodynamics for its parallels with general relativity's description of the interplay of spacetime and matter. The geometry of trajectories within the dynamical landscape of "fast time" dynamics are shaped by the structure of a differential equation and its connectivity parameters, which themselves evolve over "slow time" driven by state-dependent and state-independent plasticity mechanisms. Finally, the adjustment of plasticity processes (metaplasticity) takes place in an "ultraslow" time scale. Psychedelics flatten the neural landscape, leading to heightened entropy and complexity of neural dynamics, as observed in neuroimaging and modeling studies linking increases in complexity with a disruption of functional integration. We highlight the relationship between criticality, the complexity of fast neural dynamics, and synaptic plasticity. Pathological, rigid, or "canalized" neural dynamics result in an ultrastable confined repertoire, allowing slower plastic changes to consolidate them further. However, under the influence of psychedelics, the destabilizing emergence of complex dynamics leads to a more fluid and adaptable neural state in a process that is amplified by the plasticity-enhancing effects of psychedelics. This shift manifests as an acute systemic increase of disorder and a possibly longer-lasting increase in complexity affecting both short-term dynamics and long-term plastic processes. Our framework offers a holistic perspective on the acute effects of these substances and their potential long-term impacts on neural structure and function.

Basal forebrain cholinergic activity is necessary for upward firing rate homeostasis in the rodent visual cortex

Bidirectional homeostatic plasticity allows neurons and circuits to maintain stable firing in the face of developmental or learning-induced perturbations. In the primary visual cortex (V1), upward firing rate homeostasis (FRH) only occurs during active wake (AW) and downward during sleep, but how this behavioral state-dependent gating is accomplished is unknown. Here, we focus on how AW enables upward FRH in V1 of juvenile Long Evans rats. A major difference between quiet wake (QW), when upward FRH is absent, and AW, when it is present, is increased cholinergic (ACh) tone, and the main cholinergic projections to V1 arise from the horizontal diagonal band of the basal forebrain (HDB ACh). We therefore chemogenetically inhibited HDB ACh neurons while inducing upward homeostatic compensation using direct activity-suppression in V1. We found that synaptic scaling up and intrinsic homeostatic plasticity, two important cellular mediators of upward FRH, were both impaired when HDB ACh neurons were inhibited. Most strikingly, HDB ACh inhibition flipped the sign of intrinsic plasticity so that it became anti-homeostatic, and this effect was phenocopied by knockdown of the M1 ACh receptor in V1, indicating that this modulation of intrinsic plasticity is the result of direct actions of ACh within V1. Finally, we found that upward FRH induced by visual deprivation was completely prevented by HDB ACh inhibition. Together, our results show that HDB ACh modulation is a key enabler of upward homeostatic plasticity and FRH, and more broadly suggest that neuromodulatory inputs can segregate upward and downward homeostatic plasticity into distinct behavioral states.

Diversity-induced trivialization and resilience of neural dynamics

Heterogeneity is omnipresent across all living systems. Diversity enriches the dynamical repertoire of these systems but remains challenging to reconcile with their manifest robustness and dynamical persistence over time, a fundamental feature called resilience. To better understand the mechanism underlying resilience in neural circuits, we considered a nonlinear network model, extracting the relationship between excitability heterogeneity and resilience. To measure resilience, we quantified the number of stationary states of this network, and how they are affected by various control parameters. We analyzed both analytically and numerically gradient and non-gradient systems modeled as non-linear sparse neural networks evolving over long time scales. Our analysis shows that neuronal heterogeneity quenches the number of stationary states while decreasing the susceptibility to bifurcations: a phenomenon known as trivialization. Heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in network size and connection probability by quenching the system's dynamic volatility.

Protosequences in human cortical organoids model intrinsic states in the developing cortex

Neuronal firing sequences are thought to be the basic building blocks of neural coding and information broadcasting within the brain. However, when sequences emerge during neurodevelopment remains unknown. We demonstrate that structured firing sequences are present in spontaneous activity of human brain organoids and ex vivo neonatal brain slices from the murine somatosensory cortex. We observed a balance between temporally rigid and flexible firing patterns that are emergent phenomena in human brain organoids and early postnatal murine somatosensory cortex, but not in primary dissociated cortical cultures. Our findings suggest that temporal sequences do not arise in an experience-dependent manner, but are rather constrained by an innate preconfigured architecture established during neurogenesis. These findings highlight the potential for brain organoids to further explore how exogenous inputs can be used to refine neuronal circuits and enable new studies into the genetic mechanisms that govern assembly of functional circuitry during early human brain development.

The effect of ketamine on synaptic mistuning induced by impaired glutamate reuptake

Mistuning of synaptic transmission has been proposed to underlie many psychiatric disorders, with decreased reuptake of the excitatory neurotransmitter glutamate as one contributing factor. Synaptic tuning occurs through several diverging and converging forms of plasticity. By recording evoked field postsynaptic potentials in the CA1 area in hippocampal slices, we found that inhibiting glutamate transporters using DL-TBOA causes retuning of synaptic transmission, resulting in a new steady state with reduced synaptic strength and a lower threshold for inducing long-term synaptic potentiation (LTP). Moreover, a similar reduced threshold for LTP was observed in a rat model of depression with decreased levels of glutamate transporters. Most importantly, we found that the antidepressant ketamine counteracts the effects of increased glutamate on the various steps involved in synaptic retuning. We, therefore, propose that ketamine's mechanism of action as an antidepressant is to restore adequate synaptic tuning.

Age-related dysregulation of homeostatic control in neuronal microcircuits

Neuronal homeostasis prevents hyperactivity and hypoactivity. Age-related hyperactivity suggests homeostasis may be dysregulated in later life. However, plasticity mechanisms preventing age-related hyperactivity and their efficacy in later life are unclear. We identify the adult cortical plasticity response to elevated activity driven by sensory overstimulation, then test how plasticity changes with age. We use in vivo two-photon imaging of calcium-mediated cellular/synaptic activity, electrophysiology and c-Fos-activity tagging to show control of neuronal activity is dysregulated in the visual cortex in late adulthood. Specifically, in young adult cortex, mGluR5-dependent population-wide excitatory synaptic weakening and inhibitory synaptogenesis reduce cortical activity following overstimulation. In later life, these mechanisms are downregulated, so that overstimulation results in synaptic strengthening and elevated activity. We also find overstimulation disrupts cognition in older but not younger animals. We propose that specific plasticity mechanisms fail in later life dysregulating neuronal microcircuit homeostasis and that the age-related response to overstimulation can impact cognitive performance.

Circular and unified analysis in network neuroscience

Genuinely new discovery transcends existing knowledge. Despite this, many analyses in systems neuroscience neglect to test new speculative hypotheses against benchmark empirical facts. Some of these analyses inadvertently use circular reasoning to present existing knowledge as new discovery. Here, I discuss that this problem can confound key results and estimate that it has affected more than three thousand studies in network neuroscience over the last decade. I suggest that future studies can reduce this problem by limiting the use of speculative evidence, integrating existing knowledge into benchmark models, and rigorously testing proposed discoveries against these models. I conclude with a summary of practical challenges and recommendations.

The technology, opportunities, and challenges of Synthetic Biological Intelligence

Integrating neural cultures developed through synthetic biology methods with digital computing has enabled the early development of Synthetic Biological Intelligence (SBI). Recently, key studies have emphasized the advantages of biological neural systems in some information processing tasks. However, neither the technology behind this early development, nor the potential ethical opportunities or challenges, have been explored in detail yet. Here, we review the key aspects that facilitate the development of SBI and explore potential applications. Considering these foreseeable use cases, various ethical implications are proposed. Ultimately, this work aims to provide a robust framework to structure ethical considerations to ensure that SBI technology can be both researched and applied responsibly.

Changes in functional connectivity preserve scale-free neuronal and behavioral dynamics

Does the brain optimize itself for storage and transmission of information, and if so, how? The critical brain hypothesis is based in statistical physics and posits that the brain self-tunes its dynamics to a critical point or regime to maximize the repertoire of neuronal responses. Yet, the robustness of this regime, especially with respect to changes in the functional connectivity, remains an unsolved fundamental challenge. Here, we show that both scale-free neuronal dynamics and self-similar features of behavioral dynamics persist following significant changes in functional connectivity. Specifically, we find that the psychedelic compound ibogaine that is associated with an altered state of consciousness fundamentally alters the functional connectivity in the retrosplenial cortex of mice. Yet, the scale-free statistics of movement and of neuronal avalanches among behaviorally related neurons remain largely unaltered. This indicates that the propagation of information within biological neural networks is robust to changes in functional organization of subpopulations of neurons, opening up a new perspective on how the adaptive nature of functional networks may lead to optimality of information transmission in the brain.

Beyond pulsed inhibition: Alpha oscillations modulate attenuation and amplification of neural activity in the awake resting state

Alpha oscillations are a distinctive feature of the awake resting state of the human brain. However, their functional role in resting-state neuronal dynamics remains poorly understood. Here we show that, during resting wakefulness, alpha oscillations drive an alternation of attenuation and amplification bouts in neural activity. Our analysis indicates that inhibition is activated in pulses that last for a single alpha cycle and gradually suppress neural activity, while excitation is successively enhanced over a few alpha cycles to amplify neural activity. Furthermore, we show that long-term alpha amplitude fluctuations-the "waxing and waning" phenomenon-are an attenuation-amplification mechanism described by a power-law decay of the activity rate in the "waning" phase. Importantly, we do not observe such dynamics during non-rapid eye movement (NREM) sleep with marginal alpha oscillations. The results suggest that alpha oscillations modulate neural activity not only through pulses of inhibition (pulsed inhibition hypothesis) but also by timely enhancement of excitation (or disinhibition).

Criticality of resting-state EEG predicts perturbational complexity and level of consciousness during anesthesia

Consciousness has been proposed to be supported by electrophysiological patterns poised at criticality, a dynamical regime which exhibits adaptive computational properties, maximally complex patterns and divergent sensitivity to perturbation. Here, we investigated dynamical properties of the resting-state electroencephalogram of healthy subjects undergoing general anesthesia with propofol, xenon or ketamine. We then studied the relation of these dynamic properties with the perturbational complexity index (PCI), which has shown remarkably high sensitivity in detecting consciousness independent of behavior. All participants were unresponsive under anesthesia, while consciousness was retained only during ketamine anesthesia (in the form of vivid dreams)., enabling an experimental dissociation between unresponsiveness and unconsciousness. We estimated (i) avalanche criticality, (ii) chaoticity, and (iii) criticality-related measures, and found that states of unconsciousness were characterized by a distancing from both the edge of activity propagation and the edge of chaos. We were then able to predict individual subjects' PCI (i.e., PCImax) with a mean absolute error below 7%. Our results establish a firm link between the PCI and criticality and provide further evidence for the role of criticality in the emergence of consciousness.

Neural criticality from effective latent variables

Observations of power laws in neural activity data have raised the intriguing notion that brains may operate in a critical state. One example of this critical state is "avalanche criticality," which has been observed in various systems, including cultured neurons, zebrafish, rodent cortex, and human EEG. More recently, power laws were also observed in neural populations in the mouse under an activity coarse-graining procedure, and they were explained as a consequence of the neural activity being coupled to multiple latent dynamical variables. An intriguing possibility is that avalanche criticality emerges due to a similar mechanism. Here, we determine the conditions under which latent dynamical variables give rise to avalanche criticality. We find that populations coupled to multiple latent variables produce critical behavior across a broader parameter range than those coupled to a single, quasi-static latent variable, but in both cases, avalanche criticality is observed without fine-tuning of model parameters. We identify two regimes of avalanches, both critical but differing in the amount of information carried about the latent variable. Our results suggest that avalanche criticality arises in neural systems in which activity is effectively modeled as a population driven by a few dynamical variables and these variables can be inferred from the population activity.

Linear thinking does not reflect the newer 21st-century anesthesia concepts. A narrative review

The brain constitutes a good example of a chaotic, nonlinear biological system where large neuronal networks operate chaotically with random connectivity. This critical state is significantly affected by the anesthetic loss of consciousness induced by drugs whose pharmacological behavior has been classically based on linear kinetics and dynamics. Recent developments in pharmacology and brain monitoring during anesthesia suggest a different view that we tried to explore in this article. The concepts of effect-site for hypnotic drugs modeling a maximum effect, electroencephalographic dynamics during induction, maintenance, and recovery from anesthesia are discussed, integrated into this alternative view, and how it may be applied in daily clinical practice.

Tauopathy severely disrupts homeostatic set-points in emergent neural dynamics but not in the activity of individual neurons

The homeostatic regulation of neuronal activity is essential for robust computation; key set-points, such as firing rate, are actively stabilized to compensate for perturbations. From this perspective, the disruption of brain function central to neurodegenerative disease should reflect impairments of computationally essential set-points. Despite connecting neurodegeneration to functional outcomes, the impact of disease on set-points in neuronal activity is unknown. Here we present a comprehensive, theory-driven investigation of the effects of tau-mediated neurodegeneration on homeostatic set-points in neuronal activity. In a mouse model of tauopathy, we examine 27,000 hours of hippocampal recordings during free behavior throughout disease progression. Contrary to our initial hypothesis that tauopathy would impact set-points in spike rate and variance, we found that cell-level set-points are resilient to even the latest stages of disease. Instead, we find that tauopathy disrupts neuronal activity at the network-level, which we quantify using both pairwise measures of neuron interactions as well as measurement of the network's nearness to criticality, an ideal computational regime that is known to be a homeostatic set-point. We find that shifts in network criticality 1) track with symptoms, 2) predict underlying anatomical and molecular pathology, 3) occur in a sleep/wake dependent manner, and 4) can be used to reliably classify an animal's genotype. Our data suggest that the critical set-point is intact, but that homeostatic machinery is progressively incapable of stabilizing hippocampal networks, particularly during waking. This work illustrates how neurodegenerative processes can impact the computational capacity of neurobiological systems, and suggest an important connection between molecular pathology, circuit function, and animal behavior.

Critical dynamics arise during structured information presentation within embodied in vitro neuronal networks

Understanding how brains process information is an incredibly difficult task. Amongst the metrics characterising information processing in the brain, observations of dynamic near-critical states have generated significant interest. However, theoretical and experimental limitations associated with human and animal models have precluded a definite answer about when and why neural criticality arises with links from attention, to cognition, and even to consciousness. To explore this topic, we used an in vitro neural network of cortical neurons that was trained to play a simplified game of 'Pong' to demonstrate Synthetic Biological Intelligence (SBI). We demonstrate that critical dynamics emerge when neural networks receive task-related structured sensory input, reorganizing the system to a near-critical state. Additionally, better task performance correlated with proximity to critical dynamics. However, criticality alone is insufficient for a neuronal network to demonstrate learning in the absence of additional information regarding the consequences of previous actions. These findings offer compelling support that neural criticality arises as a base feature of incoming structured information processing without the need for higher order cognition.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Cell-autonomous diversification in bacteria arises from calcium dynamics self-organizing at a critical point

How dynamic bacterial calcium is regulated, with kinetics faster than typical mechanisms of cellular adaptation, is unknown. We discover bacterial calcium fluctuations are temporal-fractals resulting from a property known as self-organized criticality (SOC). SOC processes are poised at a phase transition separating ordered and chaotic dynamical regimes and are observed in many natural and anthropogenic systems. SOC in bacterial calcium emerges due to calcium channel coupling mediated via membrane voltage. Environmental or genetic perturbations modify calcium dynamics and the critical exponent suggesting a continuum of critical attractors. Moving along this continuum alters the collective information capacity of bacterial populations. We find that the stochastic transition from motile to sessile lifestyle is partially mediated by SOC-governed calcium fluctuations through the regulation of c-di-GMP. In summary, bacteria co-opt the physics of phase transitions to maintain dynamic calcium equilibrium, and this enables cell-autonomous population diversification during surface colonization by leveraging the stochasticity inherent at a boundary between phases.

Low dimensional criticality embedded in high dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. Here we offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality - a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescale agrees with established theory of critical dynamics. Our results suggest that cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Intrinsic neural diversity quenches the dynamic volatility of neural networks

Heterogeneity is the norm in biology. The brain is no different: Neuronal cell types are myriad, reflected through their cellular morphology, type, excitability, connectivity motifs, and ion channel distributions. While this biophysical diversity enriches neural systems' dynamical repertoire, it remains challenging to reconcile with the robustness and persistence of brain function over time (resilience). To better understand the relationship between excitability heterogeneity (variability in excitability within a population of neurons) and resilience, we analyzed both analytically and numerically a nonlinear sparse neural network with balanced excitatory and inhibitory connections evolving over long time scales. Homogeneous networks demonstrated increases in excitability, and strong firing rate correlations-signs of instability-in response to a slowly varying modulatory fluctuation. Excitability heterogeneity tuned network stability in a context-dependent way by restraining responses to modulatory challenges and limiting firing rate correlations, while enriching dynamics during states of low modulatory drive. Excitability heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in population size, connection probability, strength and variability of synaptic weights, by quenching the volatility (i.e., its susceptibility to critical transitions) of its dynamics. Together, these results highlight the fundamental role played by cell-to-cell heterogeneity in the robustness of brain function in the face of change.

Multiscale effects of excitatory-inhibitory homeostasis in lesioned cortical networks: A computational study

Stroke-related disruptions in functional connectivity (FC) often spread beyond lesioned areas and, given the localized nature of lesions, it is unclear how the recovery of FC is orchestrated on a global scale. Since recovery is accompanied by long-term changes in excitability, we propose excitatory-inhibitory (E-I) homeostasis as a driving mechanism. We present a large-scale model of the neocortex, with synaptic scaling of local inhibition, showing how E-I homeostasis can drive the post-lesion restoration of FC and linking it to changes in excitability. We show that functional networks could reorganize to recover disrupted modularity and small-worldness, but not network dynamics, suggesting the need to consider forms of plasticity beyond synaptic scaling of inhibition. On average, we observed widespread increases in excitability, with the emergence of complex lesion-dependent patterns related to biomarkers of relevant side effects of stroke, such as epilepsy, depression and chronic pain. In summary, our results show that the effects of E-I homeostasis extend beyond local E-I balance, driving the restoration of global properties of FC, and relating to post-stroke symptomatology. Therefore, we suggest the framework of E-I homeostasis as a relevant theoretical foundation for the study of stroke recovery and for understanding the emergence of meaningful features of FC from local dynamics.

Sleep and vigilance states: Embracing spatiotemporal dynamics

The classic view of sleep and vigilance states is a global stationary perspective driven by the interaction between neuromodulators and thalamocortical systems. However, recent data are challenging this view by demonstrating that vigilance states are highly dynamic and regionally complex. Spatially, sleep- and wake-like states often co-occur across distinct brain regions, as in unihemispheric sleep, local sleep in wakefulness, and during development. Temporally, dynamic switching prevails around state transitions, during extended wakefulness, and in fragmented sleep. This knowledge, together with methods monitoring brain activity across multiple regions simultaneously at millisecond resolution with cell-type specificity, is rapidly shifting how we consider vigilance states. A new perspective incorporating multiple spatial and temporal scales may have important implications for considering the governing neuromodulatory mechanisms, the functional roles of vigilance states, and their behavioral manifestations. A modular and dynamic view highlights novel avenues for finer spatiotemporal interventions to improve sleep function.

Are sleep disturbances a cause or consequence of autism spectrum disorder?

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by core symptoms such as atypical social communication, stereotyped behaviors, and restricted interests. One of the comorbid symptoms of individuals with ASD is sleep disturbance. There are two major hypotheses regarding the neural mechanism underlying ASD, i.e., the excitation/inhibition (E/I) imbalance and the altered neuroplasticity hypotheses. However, the pathology of ASD remains unclear due to inconsistent research results. This paper argues that sleep is a confounding factor, thus, must be considered when examining the pathology of ASD because sleep plays an important role in modulating the E/I balance and neuroplasticity in the human brain. Investigation of the E/I balance and neuroplasticity during sleep might enhance our understanding of the neural mechanisms of ASD. It may also lead to the development of neurobiologically informed interventions to supplement existing psychosocial therapies.

Distance to Criticality Undergoes Critical Transition Before Epileptic Seizure Attacks

Epilepsy is a common neurological disorder characterized by recurring seizures, but its underlying mechanisms remain poorly understood. Despite extensive research, there are still gaps in our knowledge about the relationship between brain dynamics and seizures. In this study, our aim is to address these gaps by proposing a novel approach to assess the role of brain network dynamics in the onset of seizures. Specifically, we investigate the relationship between brain dynamics and seizures by tracking the distance to criticality. Our hypothesis is that this distance plays a crucial role in brain state changes and that seizures may be related to critical transitions of this distance. To test this hypothesis, we develop a method to measure the evolution of the brain network's distance to the critical dynamic systems (i.e., the distance to the tipping point, DTP) using dynamic network biomarker theory and random matrix theory. The results show that the DTP of the brain decreases significantly immediately after onset of an epileptic seizure, suggesting that the brain loses its well-defined quasi-critical state during seizures. We refer to this phenomenon as the "criticality of the criticality" (COC). Furthermore, we observe that DTP exhibits a shape transition before and after the onset of the seizures. This phenomenon suggests the possibility of early warning signal (EWS) identification in the dynamic sequence of DTP, which could be utilized for seizure prediction. Our results show that the Hurst exponent, skewness, kurtosis, autocorrelation, and variance of the DTP sequence are potential EWS features. This study advances our understanding of the relationship between brain dynamics and seizures and highlights the potential for using criticality-based measures to predict and prevent seizures.

Responses of pyramidal cell somata and apical dendrites in mouse visual cortex over multiple days

The apical dendrites of pyramidal neurons in sensory cortex receive primarily top-down signals from associative and motor regions, while cell bodies and nearby dendrites are heavily targeted by locally recurrent or bottom-up inputs from the sensory periphery. Based on these differences, a number of theories in computational neuroscience postulate a unique role for apical dendrites in learning. However, due to technical challenges in data collection, little data is available for comparing the responses of apical dendrites to cell bodies over multiple days. Here we present a dataset collected through the Allen Institute Mindscope's OpenScope program that addresses this need. This dataset comprises high-quality two-photon calcium imaging from the apical dendrites and the cell bodies of visual cortical pyramidal neurons, acquired over multiple days in awake, behaving mice that were presented with visual stimuli. Many of the cell bodies and dendrite segments were tracked over days, enabling analyses of how their responses change over time. This dataset allows neuroscientists to explore the differences between apical and somatic processing and plasticity.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

Transcriptomic cell type structures in vivo neuronal activity across multiple timescales

Cell type is hypothesized to be a key determinant of a neuron's role within a circuit. Here, we examine whether a neuron's transcriptomic type influences the timing of its activity. We develop a deep-learning architecture that learns features of interevent intervals across timescales (ms to >30 min). We show that transcriptomic cell-class information is embedded in the timing of single neuron activity in the intact brain of behaving animals (calcium imaging and extracellular electrophysiology) as well as in a bio-realistic model of the visual cortex. Further, a subset of excitatory cell types are distinguishable but can be classified with higher accuracy when considering cortical layer and projection class. Finally, we show that computational fingerprints of cell types may be universalizable across structured stimuli and naturalistic movies. Our results indicate that transcriptomic class and type may be imprinted in the timing of single neuron activity across diverse stimuli.

The fractal brain: scale-invariance in structure and dynamics

The past 40 years have witnessed extensive research on fractal structure and scale-free dynamics in the brain. Although considerable progress has been made, a comprehensive picture has yet to emerge, and needs further linking to a mechanistic account of brain function. Here, we review these concepts, connecting observations across different levels of organization, from both a structural and functional perspective. We argue that, paradoxically, the level of cortical circuits is the least understood from a structural point of view and perhaps the best studied from a dynamical one. We further link observations about scale-freeness and fractality with evidence that the environment provides constraints that may explain the usefulness of fractal structure and scale-free dynamics in the brain. Moreover, we discuss evidence that behavior exhibits scale-free properties, likely emerging from similarly organized brain dynamics, enabling an organism to thrive in an environment that shares the same organizational principles. Finally, we review the sparse evidence for and try to speculate on the functional consequences of fractality and scale-freeness for brain computation. These properties may endow the brain with computational capabilities that transcend current models of neural computation and could hold the key to unraveling how the brain constructs percepts and generates behavior.

Detecting change points in neural population activity with contrastive metric learning

Finding points in time where the distribution of neural responses changes (change points) is an important step in many neural data analysis pipelines. However, in complex and free behaviors, where we see different types of shifts occurring at different rates, it can be difficult to use existing methods for change point (CP) detection because they can't necessarily handle different types of changes that may occur in the underlying neural distribution. Additionally, response changes are often sparse in high dimensional neural recordings, which can make existing methods detect spurious changes. In this work, we introduce a new approach for finding changes in neural population states across diverse activities and arousal states occurring in free behavior. Our model follows a contrastive learning approach: we learn a metric for CP detection based on maximizing the Sinkhorn divergences of neuron firing rates across two sides of a labeled CP. We apply this method to a 12-hour neural recording of a freely behaving mouse to detect changes in sleep stages and behavior. We show that when we learn a metric, we can better detect change points and also yield insights into which neurons and sub-groups are important for detecting certain types of switches that occur in the brain.

Recurrent activity in neuronal avalanches

A new statistical analysis of large neuronal avalanches observed in mouse and rat brain tissues reveals a substantial degree of recurrent activity and cyclic patterns of activation not seen in smaller avalanches. To explain these observations, we adapted a model of structural weakening in materials. In this model, dynamical weakening of neuron firing thresholds closely replicates experimental avalanche size distributions, firing number distributions, and patterns of cyclic activity. This agreement between model and data suggests that a mechanism like dynamical weakening plays a key role in recurrent activity found in large neuronal avalanches. We expect these results to illuminate the causes and dynamics of large avalanches, like those seen in seizures.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory.

Quasiuniversal scaling in mouse-brain neuronal activity stems from edge-of-instability critical dynamics

The brain is in a state of perpetual reverberant neural activity, even in the absence of specific tasks or stimuli. Shedding light on the origin and functional significance of such a dynamical state is essential to understanding how the brain transmits, processes, and stores information. An inspiring, albeit controversial, conjecture proposes that some statistical characteristics of empirically observed neuronal activity can be understood by assuming that brain networks operate in a dynamical regime with features, including the emergence of scale invariance, resembling those seen typically near phase transitions. Here, we present a data-driven analysis based on simultaneous high-throughput recordings of the activity of thousands of individual neurons in various regions of the mouse brain. To analyze these data, we construct a unified theoretical framework that synergistically combines a phenomenological renormalization group approach and techniques that infer the general dynamical state of a neural population, while designing complementary tools. This strategy allows us to uncover strong signatures of scale invariance that are "quasiuniversal" across brain regions and experiments, revealing that all the analyzed areas operate, to a greater or lesser extent, near the edge of instability.

Stochastic and deterministic dynamics in networks with excitable nodes

Networks of excitable systems provide a flexible and tractable model for various phenomena in biology, social sciences, and physics. A large class of such models undergo a continuous phase transition as the excitability of the nodes is increased. However, models of excitability that result in this continuous phase transition are based implicitly on the assumption that the probability that a node gets excited, its transfer function, is linear for small inputs. In this paper, we consider the effect of cooperative excitations, and more generally the case of a nonlinear transfer function, on the collective dynamics of networks of excitable systems. We find that the introduction of any amount of nonlinearity changes qualitatively the dynamical properties of the system, inducing a discontinuous phase transition and hysteresis. We develop a mean-field theory that allows us to understand the features of the dynamics with a one-dimensional map. We also study theoretically and numerically finite-size effects by examining the fate of initial conditions where only one node is excited in large but finite networks. Our results show that nonlinear transfer functions result in a rich effective phase diagram for finite networks, and that one should be careful when interpreting predictions of models that assume noncooperative excitations.

Scale-free behavioral dynamics directly linked with scale-free cortical dynamics

Naturally occurring body movements and collective neural activity both exhibit complex dynamics, often with scale-free, fractal spatiotemporal structure. Scale-free dynamics of both brain and behavior are important because each is associated with functional benefits to the organism. Despite their similarities, scale-free brain activity and scale-free behavior have been studied separately, without a unified explanation. Here, we show that scale-free dynamics of mouse behavior and neurons in the visual cortex are strongly related. Surprisingly, the scale-free neural activity is limited to specific subsets of neurons, and these scale-free subsets exhibit stochastic winner-take-all competition with other neural subsets. This observation is inconsistent with prevailing theories of scale-free dynamics in neural systems, which stem from the criticality hypothesis. We develop a computational model which incorporates known cell-type-specific circuit structure, explaining our findings with a new type of critical dynamics. Our results establish neural underpinnings of scale-free behavior and clear behavioral relevance of scale-free neural activity.

Neurodegenerative disease of the brain: a survey of interdisciplinary approaches

Neurodegenerative diseases of the brain pose a major and increasing global health challenge, with only limited progress made in developing effective therapies over the last decade. Interdisciplinary research is improving understanding of these diseases and this article reviews such approaches, with particular emphasis on tools and techniques drawn from physics, chemistry, artificial intelligence and psychology.

Drive competition underlies effective allostatic orchestration

Living systems ensure their fitness by self-regulating. The optimal matching of their behavior to the opportunities and demands of the ever-changing natural environment is crucial for satisfying physiological and cognitive needs. Although homeostasis has explained how organisms maintain their internal states within a desirable range, the problem of orchestrating different homeostatic systems has not been fully explained yet. In the present paper, we argue that attractor dynamics emerge from the competitive relation of internal drives, resulting in the effective regulation of adaptive behaviors. To test this hypothesis, we develop a biologically-grounded attractor model of allostatic orchestration that is embedded into a synthetic agent. Results show that the resultant neural mass model allows the agent to reproduce the navigational patterns of a rodent in an open field. Moreover, when exploring the robustness of our model in a dynamically changing environment, the synthetic agent pursues the stability of the self, being its internal states dependent on environmental opportunities to satisfy its needs. Finally, we elaborate on the benefits of resetting the model’s dynamics after drive-completion behaviors. Altogether, our studies suggest that the neural mass allostatic model adequately reproduces self-regulatory dynamics while overcoming the limitations of previous models.

Regulation of circuit organization and function through inhibitory synaptic plasticity

Diverse inhibitory neurons in the mammalian brain shape circuit connectivity and dynamics through mechanisms of synaptic plasticity. Inhibitory plasticity can establish excitation/inhibition (E/I) balance, control neuronal firing, and affect local calcium concentration, hence regulating neuronal activity at the network, single neuron, and dendritic level. Computational models can synthesize multiple experimental results and provide insight into how inhibitory plasticity controls circuit dynamics and sculpts connectivity by identifying phenomenological learning rules amenable to mathematical analysis. We highlight recent studies on the role of inhibitory plasticity in modulating excitatory plasticity, forming structured networks underlying memory formation and recall, and implementing adaptive phenomena and novelty detection. We conclude with experimental and modeling progress on the role of interneuron-specific plasticity in circuit computation and context-dependent learning.

Quasicriticality explains variability of human neural dynamics across life span

Aging impacts the brain's structural and functional organization and over time leads to various disorders, such as Alzheimer's disease and cognitive impairment. The process also impacts sensory function, bringing about a general slowing in various perceptual and cognitive functions. Here, we analyze the Cambridge Centre for Ageing and Neuroscience (Cam-CAN) resting-state magnetoencephalography (MEG) dataset-the largest aging cohort available-in light of the quasicriticality framework, a novel organizing principle for brain functionality which relates information processing and scaling properties of brain activity to brain connectivity and stimulus. Examination of the data using this framework reveals interesting correlations with age and gender of test subjects. Using simulated data as verification, our results suggest a link between changes to brain connectivity due to aging and increased dynamical fluctuations of neuronal firing rates. Our findings suggest a platform to develop biomarkers of neurological health.

Sampling effects and measurement overlap can bias the inference of neuronal avalanches

To date, it is still impossible to sample the entire mammalian brain with single-neuron precision. This forces one to either use spikes (focusing on few neurons) or to use coarse-sampled activity (averaging over many neurons, e.g. LFP). Naturally, the sampling technique impacts inference about collective properties. Here, we emulate both sampling techniques on a simple spiking model to quantify how they alter observed correlations and signatures of criticality. We describe a general effect: when the inter-electrode distance is small, electrodes sample overlapping regions in space, which increases the correlation between the signals. For coarse-sampled activity, this can produce power-law distributions even for non-critical systems. In contrast, spike recordings do not suffer this particular bias and underlying dynamics can be identified. This may resolve why coarse measures and spikes have produced contradicting results in the past.