A modular and adaptable analysis pipeline to compare slow cerebral rhythms across heterogeneous datasets.

Neuroscience is moving toward a more integrative discipline where understanding brain function requires consolidating the accumulated evidence seen across experiments, species, and measurement techniques. A remaining challenge on that path is integrating such heterogeneous data into analysis workflows such that consistent and comparable conclusions can be distilled as an experimental basis for models and theories. Here, we propose a solution in the context of slow-wave activity (<1 Hz), which occurs during unconscious brain states like sleep and general anesthesia and is observed across diverse experimental approaches. We address the issue of integrating and comparing heterogeneous data by conceptualizing a general pipeline design that is adaptable to a variety of inputs and applications. Furthermore, we present the Collaborative Brain Wave Analysis Pipeline (Cobrawap) as a concrete, reusable software implementation to perform broad, detailed, and rigorous comparisons of slow-wave characteristics across multiple, openly available electrocorticography (ECoG) and calcium imaging datasets.

Adaptation Shapes Local Cortical Reactivity: From Bifurcation Diagram and Simulations to Human Physiological and Pathological Responses.

Human studies employing intracerebral and transcranial perturbations suggest that the input-output properties of cortical circuits are dramatically affected during sleep in healthy subjects as well as in awake patients with multifocal and focal brain injury. In all these conditions, cortical circuits react to direct stimulation with an initial activation followed by suppression of activity (Off-period) that disrupts the build-up of sustained causal interactions typically observed in healthy wakefulness. The transition to this stereotypical response has important clinical implications, being associated with loss of consciousness or loss of functions. Here, we provide a mechanistic explanation of these findings by means of simulations of a cortical-like module endowed with activity-dependent adaptation and mean-field theory. First, we show that fundamental aspects of the local responses elicited in humans by direct cortical stimulation can be replicated by systematically varying the relationships between adaptation strength and excitation level in the network. Then, we reveal a region in the adaptation-excitation parameter space of crucial relevance for both physiological and pathologic conditions, where spontaneous activity and responses to perturbation diverge in their ability to reveal Off-periods. Finally, we substantiate through simulations of connected cortical-like modules the role of adaptation mechanisms in preventing cortical neurons from engaging in reciprocal causal interactions, as suggested by empirical studies. These modeling results provide a general theoretical framework and a mechanistic interpretation for a body of neurophysiological measurements that bears critical relevance for physiological states as well as for the assessment and rehabilitation of brain-injured patients.

Self-Consistent Stochastic Dynamics for Finite-Size Networks of Spiking Neurons.

Despite the huge number of neurons composing a brain network, ongoing activity of local cell assemblies is intrinsically stochastic. Fluctuations in their instantaneous rate of spike firing ν(t) scale with the size of the assembly and persist in isolated networks, i.e., in the absence of external sources of noise. Although deterministic chaos due to the quenched disorder of the synaptic couplings underlies this seemingly stochastic dynamics, an effective theory for the network dynamics of a finite assembly of spiking neurons is lacking. Here, we fill this gap by extending the so-called population density approach including an activity- and size-dependent stochastic source in the Fokker-Planck equation for the membrane potential density. The finite-size noise embedded in this stochastic partial derivative equation is analytically characterized leading to a self-consistent and nonperturbative description of ν(t) valid for a wide class of spiking neuron networks. Power spectra of ν(t) are found in excellent agreement with those from detailed simulations both in the linear regime and across a synchronization phase transition, when a size-dependent smearing of the critical dynamics emerges.

A novel network approach to multiscale biological regulation.

Modeling systems at multiple interacting scales is probably the most relevant task for pursuing a physically motivated explanation of biological regulation. In a new study, Smart and Zilman develop a convincing, albeit preliminary, model of the interplay between the cell microscale and the macroscopic tissue organization in biological systems.

Deterministic and Stochastic Components of Cortical Down States: Dynamics and Modulation.

Slow oscillations are an emergent activity of the cerebral cortex network consisting of alternating periods of activity (Up states) and silence (Down states). Up states are periods of persistent cortical activity that share properties with that of underlying wakefulness. However, the occurrence of Down states is almost invariably associated with unconsciousness, both in animal models and clinical studies. Down states have been attributed relevant functions, such as being a resetting mechanism or breaking causal interactions between cortical areas. But what do Down states consist of? Here, we explored in detail the network dynamics (e.g., synchronization and phase) during these silent periods in vivo (male mice), in vitro (ferrets, either sex), and in silico, investigating various experimental conditions that modulate them: anesthesia levels, excitability (electric fields), and excitation/inhibition balance. We identified metastability as two complementary phases composing such quiescence states: a highly synchronized "deterministic" period followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamical properties of the resulting rhythm, as well as the responsiveness to incoming inputs or refractoriness. We propose detailed Up and Down state cycle dynamics that bridge cortical properties emerging at the mesoscale with their underlying mechanisms at the microscale, providing a key to understanding unconscious states.SIGNIFICANCE STATEMENT The cerebral cortex expresses slow oscillations consisting of Up (active) and Down (silent) states. Such activity emerges not only in slow wave sleep, but also under anesthesia and in brain lesions. Down states functionally disconnect the network, and are associated with unconsciousness. Based on a large collection of data, novel data analysis approaches and computational modeling, we thoroughly investigate the nature of Down states. We identify two phases: a highly synchronized "deterministic" period, followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamic properties of the resulting rhythm and responsiveness to incoming inputs. This finding reconciles different theories of slow rhythm generation and provides clues about how the brain switches from conscious to unconscious brain states.

Neuronal population dynamics during motor plan cancellation in nonhuman primates.

To understand the cortical neuronal dynamics behind movement generation and control, most studies have focused on tasks where actions were planned and then executed using different instances of visuomotor transformations. However, to fully understand the dynamics related to movement control, one must also study how movements are actively inhibited. Inhibition, indeed, represents the first level of control both when different alternatives are available and only one solution could be adopted and when it is necessary to maintain the current position. We recorded neuronal activity from a multielectrode array in the dorsal premotor cortex (PMd) of monkeys performing a countermanding reaching task that requires, in a subset of trials, them to cancel a planned movement before its onset. In the analysis of the neuronal state space of PMd, we found a subspace in which activities conveying temporal information were confined during active inhibition and position holding. Movement execution required activities to escape from this subspace toward an orthogonal subspace and, furthermore, surpass a threshold associated with the maturation of the motor plan. These results revealed further details in the neuronal dynamics underlying movement control, extending the hypothesis that neuronal computation confined in an "output-null" subspace does not produce movements.

Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia.

In the arousal process, the brain restores its integrative activity from the synchronized state of slow wave activity (SWA). The mechanisms underpinning this state transition remain, however, to be elucidated. Here we simultaneously probed neuronal assemblies throughout the whole cortex with micro-electrocorticographic recordings in mice. We investigated the progressive shaping of propagating SWA at different levels of isoflurane. We found a form of memory of the wavefront shapes at deep anesthesia, tightly alternating posterior-anterior-posterior patterns. At low isoflurane, metastable patterns propagated in more directions, reflecting an increased complexity. The wandering across these mesostates progressively increased its randomness, as predicted by simulations of a network of spiking neurons, and confirmed in our experimental data. The complexity increase is explained by the elevated excitability of local assemblies with no modifications of the network connectivity. These results shed new light on the functional reorganization of the cortical network as anesthesia fades out.

Energy-Based Hierarchical Clustering of Cortical Slow Waves in Multi-Electrode Recordings.

The recent development of novel multi-electrode recording technologies has revealed the existence of traveling patterns of cortical activity in many species and under different states of awareness. Among these, slow activation waves occurring under sleep and anesthesia have been widely investigated as they provide unique insights into network features such as excitability, connectivity, structure, and dynamics of the cerebral cortex. Such characterization is usually based on clustering methods which are constrained by a priori assumptions as to the number of clusters to be used or rely on wave-by-wave pattern reconstruction. Here, we introduce a new computational tool based on modal analysis of fluid flows which is robustly applied to multivariate electrophysiological data from cortical networks, namely the Energy-based Hierarchical Waves Clustering method (EHWC). EHWC is composed of three main steps: (1) detecting the occurrence of global waves; (2) reducing the data dimensionality via singular value decomposition; (3) clustering hierarchically the singled-out waves. The analysis does not require the single-channel contribution to the waves, which is a typical bottleneck in this kind of analysis due to the unavoidable intrinsic variability of locally recorded activity. For testing and validation, here we used in vivo extracellular recordings from mice cortex under three different levels of anesthesia. As a result, we found slow waves with an increasing number of propagation modes as the anesthesia level decreases, giving an estimate of the increasing complexity of network dynamics. This and other wave's features replicate and extend the findings from previous literature, paving the way to extend the same approach to non-invasive electrophysiological recordings like EEG and fMRI used clinically for the characterization of brain dynamics and clinical stratification in brain lesions.

Binocular rivalry reveals an out-of-equilibrium neural dynamics suited for decision-making.

In ambiguous or conflicting sensory situations, perception is often 'multistable' in that it perpetually changes at irregular intervals, shifting abruptly between distinct alternatives. The interval statistics of these alternations exhibits quasi-universal characteristics, suggesting a general mechanism. Using binocular rivalry, we show that many aspects of this perceptual dynamics are reproduced by a hierarchical model operating out of equilibrium. The constitutive elements of this model idealize the metastability of cortical networks. Independent elements accumulate visual evidence at one level, while groups of coupled elements compete for dominance at another level. As soon as one group dominates perception, feedback inhibition suppresses supporting evidence. Previously unreported features in the serial dependencies of perceptual alternations compellingly corroborate this mechanism. Moreover, the proposed out-of-equilibrium dynamics satisfies normative constraints of continuous decision-making. Thus, multistable perception may reflect decision-making in a volatile world: integrating evidence over space and time, choosing categorically between hypotheses, while concurrently evaluating alternatives.

Attractor competition enriches cortical dynamics during awakening from anesthesia.

Slow oscillations (≲ 1 Hz), a hallmark of slow-wave sleep and deep anesthesia across species, arise from spatiotemporal patterns of activity whose complexity increases as wakefulness is approached and cognitive functions emerge. The arousal process constitutes an open window to the unknown mechanisms underlying the emergence of such dynamical richness in awake cortical networks. Here, we investigate the changes in network dynamics as anesthesia fades out in the rat visual cortex. Starting from deep anesthesia, slow oscillations gradually increase their frequency, eventually expressing maximum regularity. This stage is followed by the abrupt onset of an infra-slow (~0.2 Hz) alternation between sleep-like oscillations and activated states. A population rate model reproduces this transition driven by an increased excitability that brings it to periodically cross a critical point. Based on our model, dynamical richness emerges as a competition between two metastable attractor states, a conclusion strongly supported by the data.

Up and Down States During Slow Oscillations in Slow-Wave Sleep and Different Levels of Anesthesia.

Slow oscillations are a pattern of synchronized network activity generated by the cerebral cortex. They consist of Up and Down states, which are periods of activity interspersed with periods of silence, respectively. However, even when this is a unique dynamic regime of transitions between Up and Down states, this pattern is not constant: there is a range of oscillatory frequencies (0.1-4 Hz), and the duration of Up vs. Down states during the cycles is variable. This opens many questions. Is there a constant relationship between the duration of Up and Down states? How much do they vary across conditions and oscillatory frequencies? Are there different sub regimes within the slow oscillations? To answer these questions, we aimed to explore a concrete aspect of slow oscillations, Up and Down state durations, across three conditions: deep anesthesia, light anesthesia, and slow-wave sleep (SWS), in the same chronically implanted rats. We found that light anesthesia and SWS have rather similar properties, occupying a small area of the Up and Down state duration space. Deeper levels of anesthesia occupy a larger region of this space, revealing that a large variety of Up and Down state durations can emerge within the slow oscillatory regime. In a network model, we investigated the network parameters that can explain the different points within our bifurcation diagram in which slow oscillations are expressed.

Modulation of cortical slow oscillations and complexity across anesthesia levels.

The ability of different groups of cortical neurons to engage in causal interactions that are at once differentiated and integrated results in complex dynamic patterns. Complexity is low during periods of unconsciousness (deep sleep, anesthesia, unresponsive wakefulness syndrome) in which the brain tends to generate a stereotypical pattern consisting of alternating active and silent periods of neural activity-slow oscillations- and is high during wakefulness. But how is cortical complexity built up? Is it a continuum? An open question is whether cortical complexity can vary within the same brain state. Here we recorded with 32-channel multielectrode arrays from the cortical surface of the mouse and used both spontaneous dynamics (wave propagation entropy and functional complexity) and a perturbational approach (a variation of the perturbation complexity index) to measure complexity at different anesthesia levels. Variations in anesthesia level within the bistable regime of slow oscillations (0.1-1.5 Hz) resulted in a modulation of the slow oscillation frequency. Both perturbational and spontaneous complexity increased with decreasing anesthesia levels, in correlation with the decrease in coherence of the underlying network. Changes in complexity level are related to, but not dependent on, changes in excitability. We conclude that cortical complexity can vary within a single brain state dominated by slow oscillations, building up to the higher complexity associated with consciousness.

Analysis Pipeline for Extracting Features of Cortical Slow Oscillations.

Cortical slow oscillations (≲1 Hz) are an emergent property of the cortical network that integrate connectivity and physiological features. This rhythm, highly revealing of the characteristics of the underlying dynamics, is a hallmark of low complexity brain states like sleep, and represents a default activity pattern. Here, we present a methodological approach for quantifying the spatial and temporal properties of this emergent activity. We improved and enriched a robust analysis procedure that has already been successfully applied to both in vitro and in vivo data acquisitions. We tested the new tools of the methodology by analyzing the electrocorticography (ECoG) traces recorded from a custom 32-channel multi-electrode array in wild-type isoflurane-anesthetized mice. The enhanced analysis pipeline, named SWAP (Slow Wave Analysis Pipeline), detects Up and Down states, enables the characterization of the spatial dependency of their statistical properties, and supports the comparison of different subjects. The SWAP is implemented in a data-independent way, allowing its application to other data sets (acquired from different subjects, or with different recording tools), as well as to the outcome of numerical simulations. By using the SWAP, we report statistically significant differences in the observed slow oscillations (SO) across cortical areas and cortical sites. Computing cortical maps by interpolating the features of SO acquired at the electrode positions, we give evidence of gradients at the global scale along an oblique axis directed from fronto-lateral toward occipito-medial regions, further highlighting some heterogeneity within cortical areas. The results obtained using the SWAP will be essential for producing data-driven brain simulations. A spatial characterization of slow oscillations will also trigger a discussion on the role of, and the interplay between, the different regions in the cortex, improving our understanding of the mechanisms of generation and propagation of delta rhythms and, more generally, of cortical properties.

Dimensional reduction in networks of non-Markovian spiking neurons: Equivalence of synaptic filtering and heterogeneous propagation delays.

Message passing between components of a distributed physical system is non-instantaneous and contributes to determine the time scales of the emerging collective dynamics. In biological neuron networks this is due in part to local synaptic filtering of exchanged spikes, and in part to the distribution of the axonal transmission delays. How differently these two kinds of communication protocols affect the network dynamics is still an open issue due to the difficulties in dealing with the non-Markovian nature of synaptic transmission. Here, we develop a mean-field dimensional reduction yielding to an effective Markovian dynamics of the population density of the neuronal membrane potential, valid under the hypothesis of small fluctuations of the synaptic current. Within this limit, the resulting theory allows us to prove the formal equivalence between the two transmission mechanisms, holding for any synaptic time scale, integrate-and-fire neuron model, spike emission regimes and for different network states even when the neuron number is finite. The equivalence holds even for larger fluctuations of the synaptic input, if white noise currents are incorporated to model other possible biological features such as ionic channel stochasticity.

Scaling of a Large-Scale Simulation of Synchronous Slow-Wave and Asynchronous Awake-Like Activity of a Cortical Model With Long-Range Interconnections.

Cortical synapse organization supports a range of dynamic states on multiple spatial and temporal scales, from synchronous slow wave activity (SWA), characteristic of deep sleep or anesthesia, to fluctuating, asynchronous activity during wakefulness (AW). Such dynamic diversity poses a challenge for producing efficient large-scale simulations that embody realistic metaphors of short- and long-range synaptic connectivity. In fact, during SWA and AW different spatial extents of the cortical tissue are active in a given timespan and at different firing rates, which implies a wide variety of loads of local computation and communication. A balanced evaluation of simulation performance and robustness should therefore include tests of a variety of cortical dynamic states. Here, we demonstrate performance scaling of our proprietary Distributed and Plastic Spiking Neural Networks (DPSNN) simulation engine in both SWA and AW for bidimensional grids of neural populations, which reflects the modular organization of the cortex. We explored networks up to 192 × 192 modules, each composed of 1,250 integrate-and-fire neurons with spike-frequency adaptation, and exponentially decaying inter-modular synaptic connectivity with varying spatial decay constant. For the largest networks the total number of synapses was over 70 billion. The execution platform included up to 64 dual-socket nodes, each socket mounting 8 Intel Xeon Haswell processor cores @ 2.40 GHz clock rate. Network initialization time, memory usage, and execution time showed good scaling performances from 1 to 1,024 processes, implemented using the standard Message Passing Interface (MPI) protocol. We achieved simulation speeds of between 2.3 × 109 and 4.1 × 109 synaptic events per second for both cortical states in the explored range of inter-modular interconnections.

A Network Activity Reconfiguration Underlies the Transition from Goal to Action.

Neurons in prefrontal cortex (PF) represent mnemonic information about current goals until the action can be selected and executed. However, the neuronal dynamics underlying the transition from goal into specific actions are poorly understood. Here, we show that the goal-coding PF network is dynamically reconfigured from mnemonic to action selection states and that such reconfiguration is mediated by cell assemblies with heterogeneous excitability. We recorded neuronal activity from PF while monkeys selected their actions on the basis of memorized goals. Many PF neurons encoded the goal, but only a minority of them did so across both memory retention and action selection stages. Interestingly, about half of this minority of neurons switched their goal preference across the goal-action transition. Our computational model led us to propose a PF network composed of heterogeneous cell assemblies with single-state and bistable local dynamics able to produce both dynamical stability and input susceptibility simultaneously.

Slow Waves in Cortical Slices: How Spontaneous Activity is Shaped by Laminar Structure.

Cortical slow oscillations (SO) of neural activity spontaneously emerge and propagate during deep sleep and anesthesia and are also expressed in isolated brain slices and cortical slabs. We lack full understanding of how SO integrate the different structural levels underlying local excitability of cell assemblies and their mutual interaction. Here, we focus on ongoing slow waves (SWs) in cortical slices reconstructed from a 16-electrode array designed to probe the neuronal activity at multiple spatial scales. In spite of the variable propagation patterns observed, we reproducibly found a smooth strip of loci leading the SW fronts, overlapping cortical layers 4 and 5, along which Up states were the longest and displayed the highest firing rate. Propagation modes were uncorrelated in time, signaling a memoryless generation of SWs. All these features could be modeled by a multimodular large-scale network of spiking neurons with a specific balance between local and intermodular connectivity. Modules work as relaxation oscillators with a weakly stable Down state and a peak of local excitability to model layers 4 and 5. These conditions allow for both optimal sensitivity to the network structure and richness of propagation modes, both of which are potential substrates for dynamic flexibility in more general contexts.

State-dependent mean-field formalism to model different activity states in conductance-based networks of spiking neurons.

More interest has been shown in recent years to large-scale spiking simulations of cerebral neuronal networks, coming both from the presence of high-performance computers and increasing details in experimental observations. In this context it is important to understand how population dynamics are generated by the designed parameters of the networks, which is the question addressed by mean-field theories. Despite analytic solutions for the mean-field dynamics already being proposed for current-based neurons (CUBA), a complete analytic description has not been achieved yet for more realistic neural properties, such as conductance-based (COBA) network of adaptive exponential neurons (AdEx). Here, we propose a principled approach to map a COBA on a CUBA. Such an approach provides a state-dependent approximation capable of reliably predicting the firing-rate properties of an AdEx neuron with noninstantaneous COBA integration. We also applied our theory to population dynamics, predicting the dynamical properties of the network in very different regimes, such as asynchronous irregular and synchronous irregular (slow oscillations). This result shows that a state-dependent approximation can be successfully introduced to take into account the subtle effects of COBA integration and to deal with a theory capable of correctly predicting the activity in regimes of alternating states like slow oscillations.