The log-dynamic brain: how skewed distributions affect network operations

Formation and retrieval of cell assemblies in a biologically realistic spiking neural network model of area CA3 in the mouse hippocampus

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

Distinct brain-wide presynaptic networks underlie the functional identity of individual cortical neurons

Neuronal connections provide the scaffolding for neuronal function. Revealing the connectivity of functionally identified individual neurons is necessary to understand how activity patterns emerge and support behavior. Yet, the brain-wide presynaptic wiring rules that lay the foundation for the functional selectivity of individual neurons remain largely unexplored. Cortical neurons, even in primary sensory cortex, are heterogeneous in their selectivity, not only to sensory stimuli but also to multiple aspects of behavior. Here, to investigate presynaptic connectivity rules underlying the selectivity of pyramidal neurons to behavioral state 1-12 in primary somatosensory cortex (S1), we used two-photon calcium imaging, neuropharmacology, single-cell based monosynaptic input tracing, and optogenetics. We show that behavioral state-dependent neuronal activity patterns are stable over time. These are minimally affected by neuromodulatory inputs and are instead driven by glutamatergic inputs. Analysis of brain-wide presynaptic networks of individual neurons with distinct behavioral state-dependent activity profiles revealed characteristic patterns of anatomical input. While both behavioral state-related and unrelated neurons had a similar pattern of local inputs within S1, their long-range glutamatergic inputs differed. Individual cortical neurons, irrespective of their functional properties, received converging inputs from the main S1-projecting areas. Yet, neurons that tracked behavioral state received a smaller proportion of motor cortical inputs and a larger proportion of thalamic inputs. Optogenetic suppression of thalamic inputs reduced behavioral state-dependent activity in S1, but this activity was not externally driven. Our results revealed distinct long-range glutamatergic inputs as a substrate for preconfigured network dynamics associated with behavioral state.

Sub-threshold neuronal activity and the dynamical regime of cerebral cortex

Cortical neurons exhibit temporally irregular spiking patterns and heterogeneous firing rates. These features arise in model circuits operating in a 'fluctuation-driven regime', in which fluctuations in membrane potentials emerge from the network dynamics. However, it is still debated whether the cortex operates in such a regime. We evaluated the fluctuation-driven hypothesis by analyzing spiking and sub-threshold membrane potentials of neurons in the frontal cortex of mice performing a decision-making task. We showed that while standard fluctuation-driven models successfully account for spiking statistics, they fall short in capturing the heterogeneity in sub-threshold activity. This limitation is an inevitable outcome of bombarding single-compartment neurons with a large number of pre-synaptic inputs, thereby clamping the voltage of all neurons to more or less the same average voltage. To address this, we effectively incorporated dendritic morphology into the standard models. Inclusion of dendritic morphology in the neuronal models increased neuronal selectivity and reduced error trials, suggesting a functional role for dendrites during decision-making. Our work suggests that, during decision-making, cortical neurons in high-order cortical areas operate in a fluctuation-driven regime.

A neural algorithm for computing bipartite matchings

Finding optimal bipartite matchings-e.g., matching medical students to hospitals for residency, items to buyers in an auction, or papers to reviewers for peer review-is a fundamental combinatorial optimization problem. We found a distributed algorithm for computing matchings by studying the development of the neuromuscular circuit. The neuromuscular circuit can be viewed as a bipartite graph formed between motor neurons and muscle fibers. In newborn animals, neurons and fibers are densely connected, but after development, each fiber is typically matched (i.e., connected) to exactly one neuron. We cast this synaptic pruning process as a distributed matching (or assignment) algorithm, where motor neurons "compete" with each other to "win" muscle fibers. We show that this algorithm is simple to implement, theoretically sound, and effective in practice when evaluated on real-world bipartite matching problems. Thus, insights from the development of neural circuits can inform the design of algorithms for fundamental computational problems.

Memory-efficient neurons and synapses for spike-timing-dependent-plasticity in large-scale spiking networks

This paper addresses the challenges posed by frequent memory access during simulations of large-scale spiking neural networks involving synaptic plasticity. We focus on the memory accesses performed during a common synaptic plasticity rule since this can be a significant factor limiting the efficiency of the simulations. We propose neuron models that are represented by only three state variables, which are engineered to enforce the appropriate neuronal dynamics. Additionally, memory retrieval is executed solely by fetching postsynaptic variables, promoting a contiguous memory storage and leveraging the capabilities of burst mode operations to reduce the overhead associated with each access. Different plasticity rules could be implemented despite the adopted simplifications, each leading to a distinct synaptic weight distribution (i.e., unimodal and bimodal). Moreover, our method requires fewer average memory accesses compared to a naive approach. We argue that the strategy described can speed up memory transactions and reduce latencies while maintaining a small memory footprint.

Oscillatory-Quality of sleep spindles links brain state with sleep regulation and function

Here, we characterized the dynamics of sleep spindles, focusing on their damping, which we estimated using a metric called oscillatory-Quality (o-Quality), derived by fitting an autoregressive model to electrophysiological signals, recorded from the cortex in mice. The o-Quality of sleep spindles correlates weakly with their amplitude, shows marked laminar differences and regional topography across cortical regions, reflects the level of synchrony within and between cortical networks, is strongly modulated by sleep-wake history, reflects the degree of sensory disconnection, and correlates with the strength of coupling between spindles and slow waves. As most spindle events are highly localized and not detectable with conventional low-density recording approaches, o-Quality thus emerges as a valuable metric that allows us to infer the spread and dynamics of spindle activity across the brain and directly links their spatiotemporal dynamics with local and global regulation of brain states, sleep regulation, and function.

2P-NucTag: on-demand phototagging for molecular analysis of functionally identified cortical neurons

Neural circuits are characterized by genetically and functionally diverse cell types. A mechanistic understanding of circuit function is predicated on linking the genetic and physiological properties of individual neurons. However, it remains highly challenging to map the transcriptional properties to functionally heterogeneous neuronal subtypes in mammalian cortical circuits in vivo. Here, we introduce a high-throughput two-photon nuclear phototagging (2P-NucTag) approach optimized for on-demand and indelible labeling of single neurons via a photoactivatable red fluorescent protein following in vivo functional characterization in behaving mice. We demonstrate the utility of this function-forward pipeline by selectively labeling and transcriptionally profiling previously inaccessible place and silent cells in the mouse hippocampus. Our results reveal unexpected differences in gene expression between these hippocampal pyramidal neurons with distinct spatial coding properties. Thus, 2P-NucTag opens a new way to uncover the molecular principles that govern the functional organization of neural circuits.

Linking brain activity across scales with simultaneous opto- and electrophysiology

The brain enables adaptive behavior via the dynamic coordination of diverse neuronal signals across spatial and temporal scales: from fast action potential patterns in microcircuits to slower patterns of distributed activity in brain-wide networks. Understanding principles of multiscale dynamics requires simultaneous monitoring of signals in multiple, distributed network nodes. Combining optical and electrical recordings of brain activity is promising for collecting data across multiple scales and can reveal aspects of coordinated dynamics invisible to standard, single-modality approaches. We review recent progress in combining opto- and electrophysiology, focusing on mouse studies that shed new light on the function of single neurons by embedding their activity in the context of brain-wide activity patterns. Optical and electrical readouts can be tailored to desired scales to tackle specific questions. For example, fast dynamics in single cells or local populations recorded with multi-electrode arrays can be related to simultaneously acquired optical signals that report activity in specified subpopulations of neurons, in non-neuronal cells, or in neuromodulatory pathways. Conversely, two-photon imaging can be used to densely monitor activity in local circuits while sampling electrical activity in distant brain areas at the same time. The refinement of combined approaches will continue to reveal previously inaccessible and under-appreciated aspects of coordinated brain activity.

Neuronal firing rate diversity lowers the dimension of population covariability

Populations of neurons produce activity with two central features. First, neuronal responses are very diverse - specific stimuli or behaviors prompt some neurons to emit many action potentials, while other neurons remain relatively silent. Second, the trial-to-trial fluctuations of neuronal response occupy a low dimensional space, owing to significant correlations between the activity of neurons. These two features define the quality of neuronal representation. We link these two aspects of population response using a recurrent circuit model and derive the following relation: the more diverse the firing rates of neurons in a population, the lower the effective dimension of population trial-to-trial covariability. This surprising prediction is tested and validated using simultaneously recorded neuronal populations from numerous brain areas in mice, non-human primates, and in the motor cortex of human participants. Using our relation we present a theory where a more diverse neuronal code leads to better fine discrimination performance from population activity. In line with this theory, we show that neuronal populations across the brain exhibit both more diverse mean responses and lower-dimensional fluctuations when the brain is in more heightened states of information processing. In sum, we present a key organizational principle of neuronal population response that is widely observed across the nervous system and acts to synergistically improve population representation.

A high-density multi-electrode platform examining the effects of radiation on in vitro cortical networks

Radiation therapy and stereotactic radiosurgery are common treatments for brain malignancies. However, the impact of radiation on underlying neuronal circuits is poorly understood. In the prefrontal cortex (PFC), neurons communicate via action potentials that control cognitive processes, thus it is important to understand the impact of radiation on these circuits. Here we present a novel protocol to investigate the effect of radiation on the activity and survival of PFC networks in vitro. Escalating doses of radiation were applied to PFC slices using a robotic radiosurgery platform at a standard dose rate of 10 Gy/min. High-density multielectrode array recordings of radiated slices were collected to capture extracellular activity across 4,096 channels. Radiated slices showed an increase in firing rate, functional connectivity, and complexity. Graph-theoretic measures of functional connectivity were altered following radiation. These results were compared to pharmacologically induced epileptic slices where neural complexity was markedly elevated, and functional connections were strong but remained spatially focused. Finally, propidium iodide staining revealed a dose-dependent effect of radiation on apoptosis. These findings provide a novel assay to investigate the impacts of clinically relevant doses of radiation on brain circuits and highlight the acute effects of escalating radiation doses on PFC neurons.

Failure in a population: Tauopathy disrupts homeostatic set-points in emergent dynamics despite stability in the constituent neurons

Homeostatic regulation of neuronal activity is essential for robust computation; set-points, such as firing rate, are actively stabilized to compensate for perturbations. The disruption of brain function central to neurodegenerative disease likely arises from impairments of computationally essential set-points. Here, we systematically investigated the effects of tau-mediated neurodegeneration on all known set-points in neuronal activity. We continuously tracked hippocampal neuronal activity across the lifetime of a mouse model of tauopathy. We were unable to detect effects of disease in measures of single-neuron firing activity. By contrast, as tauopathy progressed, there was disruption of network-level neuronal activity, quantified by measuring neuronal pairwise interactions and criticality, a homeostatically controlled, ideal computational regime. Deviations in criticality correlated with symptoms, predicted underlying anatomical pathology, occurred in a sleep-wake-dependent manner, and could be used to reliably classify an animal's genotype. This work illustrates how neurodegeneration may disrupt the computational capacity of neurobiological systems.

Practical parameter identifiability and handling of censored data with Bayesian inference in mathematical tumour models

Mechanistic mathematical models (MMs) are a powerful tool to help us understand and predict the dynamics of tumour growth under various conditions. In this work, we use 5 MMs with an increasing number of parameters to explore how certain (often overlooked) decisions in estimating parameters from data of experimental tumour growth affect the outcome of the analysis. In particular, we propose a framework for including tumour volume measurements that fall outside the upper and lower limits of detection, which are normally discarded. We demonstrate how excluding censored data results in an overestimation of the initial tumour volume and the MM-predicted tumour volumes prior to the first measurements, and an underestimation of the carrying capacity and the MM-predicted tumour volumes beyond the latest measurable time points. We show in which way the choice of prior for the MM parameters can impact the posterior distributions, and illustrate that reporting the most likely parameters and their 95% credible interval can lead to confusing or misleading interpretations. We hope this work will encourage others to carefully consider choices made in parameter estimation and to adopt the approaches we put forward herein.

A comprehensive data-driven model of cat primary visual cortex

Knowledge integration based on the relationship between structure and function of the neural substrate is one of the main targets of neuroinformatics and data-driven computational modeling. However, the multiplicity of data sources, the diversity of benchmarks, the mixing of observables of different natures, and the necessity of a long-term, systematic approach make such a task challenging. Here we present a first snapshot of a long-term integrative modeling program designed to address this issue in the domain of the visual system: a comprehensive spiking model of cat primary visual cortex. The presented model satisfies an extensive range of anatomical, statistical and functional constraints under a wide range of visual input statistics. In the presence of physiological levels of tonic stochastic bombardment by spontaneous thalamic activity, the modeled cortical reverberations self-generate a sparse asynchronous ongoing activity that quantitatively matches a range of experimentally measured statistics. When integrating feed-forward drive elicited by a high diversity of visual contexts, the simulated network produces a realistic, quantitatively accurate interplay between visually evoked excitatory and inhibitory conductances; contrast-invariant orientation-tuning width; center surround interactions; and stimulus-dependent changes in the precision of the neural code. This integrative model offers insights into how the studied properties interact, contributing to a better understanding of visual cortical dynamics. It provides a basis for future development towards a comprehensive model of low-level perception.

Machine learning of brain-specific biomarkers from EEG

BACKGROUND

Electroencephalography (EEG) has a long history as a clinical tool to study brain function, and its potential to derive biomarkers for various applications is far from exhausted. Machine learning (ML) can guide future innovation by harnessing the wealth of complex EEG signals to isolate relevant brain activity. Yet, ML studies in EEG tend to ignore physiological artefacts, which may cause problems for deriving biomarkers specific to the central nervous system (CNS).

METHODS

We present a framework for conceptualising machine learning from CNS versus peripheral signals measured with EEG. A signal representation based on Morlet wavelets allowed us to define traditional brain activity features (e.g. log power) and alternative inputs used by state-of-the-art ML approaches based on covariance matrices. Using more than 2600 EEG recordings from large public databases (TUAB, TDBRAIN), we studied the impact of peripheral signals and artefact removal techniques on ML models in age and sex prediction analyses.

FINDINGS

Across benchmarks, basic artefact rejection improved model performance, whereas further removal of peripheral signals using ICA decreased performance. Our analyses revealed that peripheral signals enable age and sex prediction. However, they explained only a fraction of the performance provided by brain signals.

INTERPRETATION

We show that brain signals and body signals, both present in the EEG, allow for prediction of personal characteristics. While these results may depend on specific applications, our work suggests that great care is needed to separate these signals when the goal is to develop CNS-specific biomarkers using ML.

FUNDING

All authors have been working for F. Hoffmann-La Roche Ltd.

Dynamics of spike transmission and suppression between principal cells and interneurons in the hippocampus and entorhinal cortex

Synaptic excitation and inhibition are essential for neuronal communication. However, the variables that regulate synaptic excitation and inhibition in the intact brain remain largely unknown. Here, we examined how spike transmission and suppression between principal cells (PCs) and interneurons (INTs) are modulated by activity history, brain state, cell type, and somatic distance between presynaptic and postsynaptic neurons by applying cross-correlogram analyses to datasets recorded from the dorsal hippocampus and medial entorhinal cortex (MEC) of 11 male behaving and sleeping Long Evans rats. The strength, temporal delay, and brain-state dependency of the spike transmission and suppression depended on the subregions/layers. The spike transmission probability of PC-INT excitatory pairs that showed short-term depression versus short-term facilitation was higher in CA1 and lower in CA3. Likewise, the intersomatic distance affected the proportion of PC-INT excitatory pairs that showed short-term depression and facilitation in the opposite manner in CA1 compared with CA3. The time constant of depression was longer, while that of facilitation was shorter in MEC than in CA1 and CA3. During sharp-wave ripples, spike transmission showed a larger gain in the MEC than in CA1 and CA3. The intersomatic distance affected the spike transmission gain during sharp-wave ripples differently in CA1 versus CA3. A subgroup of MEC layer 3 (EC3) INTs preferentially received excitatory inputs from and inhibited MEC layer 2 (EC2) PCs. The EC2 PC-EC3 INT excitatory pairs, most of which showed short-term depression, exhibited higher spike transmission probabilities than the EC2 PC-EC2 INT and EC3 PC-EC3 INT excitatory pairs. EC2 putative stellate cells exhibited stronger spike transmission to and received weaker spike suppression from EC3 INTs than EC2 putative pyramidal cells. This study provides detailed comparisons of monosynaptic interaction dynamics in the hippocampal-entorhinal loop, which may help to elucidate circuit operations.

Keeping Your Brain in Balance: Homeostatic Regulation of Network Function

To perform computations with the efficiency necessary for animal survival, neocortical microcircuits must be capable of reconfiguring in response to experience, while carefully regulating excitatory and inhibitory connectivity to maintain stable function. This dynamic fine-tuning is accomplished through a rich array of cellular homeostatic plasticity mechanisms that stabilize important cellular and network features such as firing rates, information flow, and sensory tuning properties. Further, these functional network properties can be stabilized by different forms of homeostatic plasticity, including mechanisms that target excitatory or inhibitory synapses, or that regulate intrinsic neuronal excitability. Here we discuss which aspects of neocortical circuit function are under homeostatic control, how this homeostasis is realized on the cellular and molecular levels, and the pathological consequences when circuit homeostasis is impaired. A remaining challenge is to elucidate how these diverse homeostatic mechanisms cooperate within complex circuits to enable them to be both flexible and stable.

Emotion brain network topology in healthy subjects following passive listening to different auditory stimuli

A large body of research establishes the efficacy of musical intervention in many aspects of physical, cognitive, communication, social, and emotional rehabilitation. However, the underlying neural mechanisms for musical therapy remain elusive. This study aimed to investigate the potential neural correlates of musical therapy, focusing on the changes in the topology of emotion brain network. To this end, a Bayesian statistical approach and a cross-over experimental design were employed together with two resting-state magnetoencephalography (MEG) as controls. MEG recordings of 30 healthy subjects were acquired while listening to five auditory stimuli in random order. Two resting-state MEG recordings of each subject were obtained, one prior to the first stimulus (pre) and one after the final stimulus (post). Time series at the level of brain regions were estimated using depth-weighted minimum norm estimation (wMNE) source reconstruction method and the functional connectivity between these regions were computed. The resultant connectivity matrices were used to derive two topological network measures: transitivity and global efficiency which are important in gauging the functional segregation and integration of brain network respectively. The differences in these measures between pre- and post-stimuli resting MEG were set as the equivalence regions. We found that the network measures under all auditory stimuli were equivalent to the resting state network measures in all frequency bands, indicating that the topology of the functional brain network associated with emotional regulation in healthy subjects remains unchanged following these auditory stimuli. This suggests that changes in the emotion network topology may not be the underlying neural mechanism of musical therapy. Nonetheless, further studies are required to explore the neural mechanisms of musical interventions especially in the populations with neuropsychiatric disorders.

Coordinated Response Modulations Enable Flexible Use of Visual Information

We use sensory information in remarkably flexible ways. We can generalize by ignoring task-irrelevant features, report different features of a stimulus, and use different actions to report a perceptual judgment. These forms of flexible behavior are associated with small modulations of the responses of sensory neurons. While the existence of these response modulations is indisputable, efforts to understand their function have been largely relegated to theory, where they have been posited to change information coding or enable downstream neurons to read out different visual and cognitive information using flexible weights. Here, we tested these ideas using a rich, flexible behavioral paradigm, multi-neuron, multi-area recordings in primary visual cortex (V1) and mid-level visual area V4. We discovered that those response modulations in V4 (but not V1) contain the ingredients necessary to enable flexible behavior, but not via those previously hypothesized mechanisms. Instead, we demonstrated that these response modulations are precisely coordinated across the population such that downstream neurons have ready access to the correct information to flexibly guide behavior without making changes to information coding or synapses. Our results suggest a novel computational role for task-dependent response modulations: they enable flexible behavior by changing the information that gets out of a sensory area, not by changing information coding within it.

Rényi entropy-complexity causality space: a novel neurocomputational tool for detecting scale-free features in EEG/iEEG data

Scale-free brain activity, linked with learning, the integration of different time scales, and the formation of mental models, is correlated with a metastable cognitive basis. The spectral slope, a key aspect of scale-free dynamics, was proposed as a potential indicator to distinguish between different sleep stages. Studies suggest that brain networks maintain a consistent scale-free structure across wakefulness, anesthesia, and recovery. Although differences in anesthetic sensitivity between the sexes are recognized, these variations are not evident in clinical electroencephalographic recordings of the cortex. Recently, changes in the slope of the power law exponent of neural activity were found to correlate with changes in Rényi entropy, an extended concept of Shannon's information entropy. These findings establish quantifiers as a promising tool for the study of scale-free dynamics in the brain. Our study presents a novel visual representation called the Rényi entropy-complexity causality space, which encapsulates complexity, permutation entropy, and the Rényi parameter q. The main goal of this study is to define this space for classical dynamical systems within theoretical bounds. In addition, the study aims to investigate how well different time series mimicking scale-free activity can be discriminated. Finally, this tool is used to detect dynamic features in intracranial electroencephalography (iEEG) signals. To achieve these goals, the study implementse the Bandt and Pompe method for ordinal patterns. In this process, each signal is associated with a probability distribution, and the causal measures of Rényi entropy and complexity are computed based on the parameter q. This method is a valuable tool for analyzing simulated time series. It effectively distinguishes elements of correlated noise and provides a straightforward means of examining differences in behaviors, characteristics, and classifications. For the iEEG experimental data, the REM state showed a greater number of significant sex-based differences, while the supramarginal gyrus region showed the most variation across different modes and analyzes. Exploring scale-free brain activity with this framework could provide valuable insights into cognition and neurological disorders. The results may have implications for understanding differences in brain function between the sexes and their possible relevance to neurological disorders.

Self-organizing recruitment of compensatory areas maximizes residual motor performance post-stroke

Whereas the orderly recruitment of compensatory motor cortical areas after stroke depends on the size of the motor cortex lesion affecting arm and hand movements, the mechanisms underlying this reorganization are unknown. Here, we hypothesized that the recruitment of compensatory areas results from the motor system's goal to optimize performance given the anatomical constraints before and after the lesion. This optimization is achieved through two complementary plastic processes: a homeostatic regulation process, which maximizes information transfer in sensory-motor networks, and a reinforcement learning process, which minimizes movement error and effort. To test this hypothesis, we developed a neuro-musculoskeletal model that controls a 7-muscle planar arm via a cortical network that includes a primary motor cortex and a premotor cortex that directly project to spinal motor neurons, and a contra-lesional primary motor cortex that projects to spinal motor neurons via the reticular formation. Synapses in the cortical areas are updated via reinforcement learning and the activity of spinal motor neurons is adjusted through homeostatic regulation. The model replicated neural, muscular, and behavioral outcomes in both non-lesioned and lesioned brains. With increasing lesion sizes, the model demonstrated systematic recruitment of the remaining primary motor cortex, premotor cortex, and contra-lesional cortex. The premotor cortex acted as a reserve area for fine motor control recovery, while the contra-lesional cortex helped avoid paralysis at the cost of poor joint control. Plasticity in spinal motor neurons enabled force generation after large cortical lesions despite weak corticospinal inputs. Compensatory activity in the premotor and contra-lesional motor cortex was more prominent in the early recovery period, gradually decreasing as the network minimized effort. Thus, the orderly recruitment of compensatory areas following strokes of varying sizes results from biologically plausible local plastic processes that maximize performance, whether the brain is intact or lesioned.

Revisiting alpha resting state dynamics underlying hallucinatory vulnerability: Insights from Hidden semi-Markov Modeling

BACKGROUND

Resting state (RS) brain activity is inherently non-stationary. Hidden semi-Markov Models (HsMM) can characterize continuous RS data as a sequence of recurring and distinct brain states along with their spatio-temporal dynamics.

NEW METHOD

Recent explorations suggest that HsMM state dynamics in the alpha frequency band link to auditory hallucination proneness (HP) in non-clinical individuals. The present study aimed to replicate these findings to elucidate robust neural correlates of hallucinatory vulnerability. Specifically, we aimed to investigate the reproducibility of HsMM states across different data sets and within-data set variants as well as the replicability of the association between alpha brain state dynamics and HP.

RESULTS

We found that most brain states are reproducible in different data sets, confirming that the HsMM characterized robust and generalizable EEG RS dynamics on a sub-second timescale. Brain state topographies and temporal dynamics of different within-data set variants showed substantial similarities and were robust against reduced data length and number of electrodes. However, the association with HP was not directly reproducible across data sets.

COMPARISON WITH EXISTING METHODS

The HsMM optimally leverages the high temporal resolution of EEG data and overcomes time-domain restrictions of other state allocation methods.

CONCLUSION

The results indicate that the sensitivity of brain state dynamics to capture individual variability in HP may depend on the data recording characteristics and individual variability in RS cognition, such as mind wandering. Future studies should consider that the order in which eyes-open and eyes-closed RS data are acquired directly influences an individual's attentional state and generation of spontaneous thoughts, and thereby might mediate the link to hallucinatory vulnerability.

L-Type Calcium Channel Modulates Low-Intensity Pulsed Ultrasound-Induced Excitation in Cultured Hippocampal Neurons

As a noninvasive technique, ultrasound stimulation is known to modulate neuronal activity both in vitro and in vivo. The latest explanation of this phenomenon is that the acoustic wave can activate the ion channels and further impact the electrophysiological properties of targeted neurons. However, the underlying mechanism of low-intensity pulsed ultrasound (LIPUS)-induced neuro-modulation effects is still unclear. Here, we characterize the excitatory effects of LIPUS on spontaneous activity and the intracellular Ca2+ homeostasis in cultured hippocampal neurons. By whole-cell patch clamp recording, we found that 15 min of 1-MHz LIPUS boosts the frequency of both spontaneous action potentials and spontaneous excitatory synaptic currents (sEPSCs) and also increases the amplitude of sEPSCs in hippocampal neurons. This phenomenon lasts for > 10 min after LIPUS exposure. Together with Ca2+ imaging, we clarified that LIPUS increases the [Ca2+]cyto level by facilitating L-type Ca2+ channels (LTCCs). In addition, due to the [Ca2+]cyto elevation by LIPUS exposure, the Ca2+-dependent CaMKII-CREB pathway can be activated within 30 min to further regulate the gene transcription and protein expression. Our work suggests that LIPUS regulates neuronal activity in a Ca2+-dependent manner via LTCCs. This may also explain the multi-activation effects of LIPUS beyond neurons. LIPUS stimulation potentiates spontaneous neuronal activity by increasing Ca2+ influx.

Maintenance of Lognormal-Like Skewed Dendritic Spine Size Distributions in Dentate Granule Cells of TNF, TNF-R1, TNF-R2, and TNF-R1/2-Deficient Mice

Dendritic spines are sites of synaptic plasticity and their head size correlates with the strength of the corresponding synapse. We recently showed that the distribution of spine head sizes follows a lognormal-like distribution even after blockage of activity or plasticity induction. As the cytokine tumor necrosis factor (TNF) influences synaptic transmission and constitutive TNF and receptor (TNF-R)-deficiencies cause changes in spine head size distributions, we tested whether these genetic alterations disrupt the lognormality of spine head sizes. Furthermore, we distinguished between spines containing the actin-modulating protein synaptopodin (SP-positive), which is present in large, strong and stable spines and those lacking it (SP-negative). Our analysis revealed that neither TNF-deficiency nor the absence of TNF-R1, TNF-R2 or TNF-R 1 and 2 (TNF-R1/R2) degrades the general lognormal-like, skewed distribution of spine head sizes (all spines, SP-positive spines, SP-negative spines). However, TNF, TNF-R1 and TNF-R2-deficiency affected the width of the lognormal distribution, and TNF-R1/2-deficiency shifted the distribution to the left. Our findings demonstrate the robustness of the lognormal-like, skewed distribution, which is maintained even in the face of genetic manipulations that alter the distribution of spine head sizes. Our observations are in line with homeostatic adaptation mechanisms of neurons regulating the distribution of spines and their head sizes.

Nondifferentiable activity in the brain

Spike raster plots of numerous neurons show vertical stripes, indicating that neurons exhibit synchronous activity in the brain. We seek to determine whether these coherent dynamics are caused by smooth brainwave activity or by something else. By analyzing biological data, we find that their cross-correlograms exhibit not only slow undulation but also a cusp at the origin, in addition to possible signs of monosynaptic connectivity. Here we show that undulation emerges if neurons are subject to smooth brainwave oscillations while a cusp results from nondifferentiable fluctuations. While modern analysis methods have achieved good connectivity estimation by adapting the models to slow undulation, they still make false inferences due to the cusp. We devise a new analysis method that may solve both problems. We also demonstrate that oscillations and nondifferentiable fluctuations may emerge in simulations of large-scale neural networks.

Probability density and information entropy of machine learning derived intracranial pressure predictions

Even with the powerful statistical parameters derived from the Extreme Gradient Boost (XGB) algorithm, it would be advantageous to define the predicted accuracy to the level of a specific case, particularly when the model output is used to guide clinical decision-making. The probability density function (PDF) of the derived intracranial pressure predictions enables the computation of a definite integral around a point estimate, representing the event's probability within a range of values. Seven hold-out test cases used for the external validation of an XGB model underwent retinal vascular pulse and intracranial pressure measurement using modified photoplethysmography and lumbar puncture, respectively. The definite integral ±1 cm water from the median (DIICP) demonstrated a negative and highly significant correlation (-0.5213±0.17, p< 0.004) with the absolute difference between the measured and predicted median intracranial pressure (DiffICPmd). The concordance between the arterial and venous probability density functions was estimated using the two-sample Kolmogorov-Smirnov statistic, extending the distribution agreement across all data points. This parameter showed a statistically significant and positive correlation (0.4942±0.18, p< 0.001) with DiffICPmd. Two cautionary subset cases (Case 8 and Case 9), where disagreement was observed between measured and predicted intracranial pressure, were compared to the seven hold-out test cases. Arterial predictions from both cautionary subset cases converged on a uniform distribution in contrast to all other cases where distributions converged on either log-normal or closely related skewed distributions (gamma, logistic, beta). The mean±standard error of the arterial DIICP from cases 8 and 9 (3.83±0.56%) was lower compared to that of the hold-out test cases (14.14±1.07%) the between group difference was statistically significant (p<0.03). Although the sample size in this analysis was limited, these results support a dual and complementary analysis approach from independently derived retinal arterial and venous non-invasive intracranial pressure predictions. Results suggest that plotting the PDF and calculating the lower order moments, arterial DIICP, and the two sample Kolmogorov-Smirnov statistic may provide individualized predictive accuracy parameters.

Preconfigured architecture of the developing mouse brain

In the adult brain, structural and functional parameters, such as synaptic sizes and neuronal firing rates, follow right-skewed and heavy-tailed distributions. While this organization is thought to have significant implications, its development is still largely unknown. Here, we address this knowledge gap by investigating a large-scale dataset recorded from the prefrontal cortex and the olfactory bulb of mice aged 4-60 postnatal days. We show that firing rates and spike train interactions have a largely stable distribution shape throughout the first 60 postnatal days and that the prefrontal cortex displays a functional small-world architecture. Moreover, early brain activity exhibits an oligarchical organization, where high-firing neurons have hub-like properties. In a neural network model, we show that analogously right-skewed and heavy-tailed synaptic parameters are instrumental to consistently recapitulate the experimental data. Thus, functional and structural parameters in the developing brain are already extremely distributed, suggesting that this organization is preconfigured and not experience dependent.

Biological direct-shortcut deep residual learning for sparse visual processing

We show, based on the following three grounds, that the primary visual cortex (V1) is a biological direct-shortcut deep residual learning neural network (ResNet) for sparse visual processing: (1) We first highlight that Gabor-like sets of basis functions, which are similar to the receptive fields of simple cells in the primary visual cortex (V1), are excellent candidates for sparse representation of natural images; i.e., images from the natural world, affirming the brain to be optimized for this. (2) We then prove that the intra-layer synaptic weight matrices of this region can be reasonably first-order approximated by identity mappings, and are thus sparse themselves. (3) Finally, we point out that intra-layer weight matrices having identity mapping as their initial approximation, irrespective of this approximation being also a reasonable first-order one or not, resemble the building blocks of direct-shortcut digital ResNets, which completes the grounds. This biological ResNet interconnects the sparsity of the final representation of the image to that of its intra-layer weights. Further exploration of this ResNet, and understanding the joint effects of its architecture and learning rules, e.g. on its inductive bias, could lead to major advancements in the area of bio-inspired digital ResNets. One immediate line of research in this context, for instance, is to study the impact of forcing the direct-shortcuts to be good first-order approximations of each building block. For this, along with the ℓ 1 -minimization posed on the basis function coefficients the ResNet finally provides at its output, another parallel one could e.g. also be posed on the weights of its residual layers.

Increasing Robustness of Intracortical Brain-Computer Interfaces for Recording Condition Changes via Data Augmentation

BACKGROUND AND OBJECTIVE

Intracortical brain-computer interfaces (iBCIs) aim to help paralyzed individuals restore their motor functions by decoding neural activity into intended movement. However, changes in neural recording conditions hinder the decoding performance of iBCIs, mainly because the neural-to-kinematic mappings shift. Conventional approaches involve either training the neural decoders using large datasets before deploying the iBCI or conducting frequent calibrations during its operation. However, collecting data for extended periods can cause user fatigue, negatively impacting the quality and consistency of neural signals. Furthermore, frequent calibration imposes a substantial computational load.

METHODS

This study proposes a novel approach to increase iBCIs' robustness against changing recording conditions. The approach uses three neural augmentation operators to generate augmented neural activity that mimics common recording conditions. Then, contrastive learning is used to learn latent factors by maximizing the similarity between the augmented neural activities. The learned factors are expected to remain stable despite varying recording conditions and maintain a consistent correlation with the intended movement.

RESULTS

Experimental results demonstrate that the proposed iBCI outperformed the state-of-the-art iBCIs and was robust to changing recording conditions across days for long-term use on one publicly available nonhuman primate dataset. It achieved satisfactory offline decoding performance, even when a large training dataset was unavailable.

CONCLUSIONS

This study paves the way for reducing the need for frequent calibration of iBCIs and collecting a large amount of annotated training data. Potential future works aim to improve offline decoding performance with an ultra-small training dataset and improve the iBCIs' robustness to severely disabled electrodes.

NeuroDots: From Single-Target to Brain-Network Modulation: Why and What Is Needed?

OBJECTIVES

Current techniques in brain stimulation are still largely based on a phrenologic approach that a single brain target can treat a brain disorder. Nevertheless, meta-analyses of brain implants indicate an overall success rate of 50% improvement in 50% of patients, irrespective of the brain-related disorder. Thus, there is still a large margin for improvement. The goal of this manuscript is to 1) develop a general theoretical framework of brain functioning that is amenable to surgical neuromodulation, and 2) describe the engineering requirements of the next generation of implantable brain stimulators that follow from this theoretic model.

MATERIALS AND METHODS

A neuroscience and engineering literature review was performed to develop a universal theoretical model of brain functioning and dysfunctioning amenable to surgical neuromodulation.

RESULTS

Even though a single target can modulate an entire network, research in network science reveals that many brain disorders are the consequence of maladaptive interactions among multiple networks rather than a single network. Consequently, targeting the main connector hubs of those multiple interacting networks involved in a brain disorder is theoretically more beneficial. We, thus, envision next-generation network implants that will rely on distributed, multisite neuromodulation targeting correlated and anticorrelated interacting brain networks, juxtaposing alternative implant configurations, and finally providing solid recommendations for the realization of such implants. In doing so, this study pinpoints the potential shortcomings of other similar efforts in the field, which somehow fall short of the requirements.

CONCLUSION

The concept of network stimulation holds great promise as a universal approach for treating neurologic and psychiatric disorders.

Advancing brain-inspired computing with hybrid neural networks

Brain-inspired computing, drawing inspiration from the fundamental structure and information-processing mechanisms of the human brain, has gained significant momentum in recent years. It has emerged as a research paradigm centered on brain-computer dual-driven and multi-network integration. One noteworthy instance of this paradigm is the hybrid neural network (HNN), which integrates computer-science-oriented artificial neural networks (ANNs) with neuroscience-oriented spiking neural networks (SNNs). HNNs exhibit distinct advantages in various intelligent tasks, including perception, cognition and learning. This paper presents a comprehensive review of HNNs with an emphasis on their origin, concepts, biological perspective, construction framework and supporting systems. Furthermore, insights and suggestions for potential research directions are provided aiming to propel the advancement of the HNN paradigm.

Synaptic architecture of a memory engram in the mouse hippocampus

Memory engrams are formed through experience-dependent remodeling of neural circuits, but their detailed architectures have remained unresolved. Using 3D electron microscopy, we performed nanoscale reconstructions of the hippocampal CA3-CA1 pathway following chemogenetic labeling of cellular ensembles with a remote history of correlated excitation during associative learning. Projection neurons involved in memory acquisition expanded their connectomes via multi-synaptic boutons without altering the numbers and spatial arrangements of individual axonal terminals and dendritic spines. This expansion was driven by presynaptic activity elicited by specific negative valence stimuli, regardless of the co-activation state of postsynaptic partners. The rewiring of initial ensembles representing an engram coincided with local, input-specific changes in the shapes and organelle composition of glutamatergic synapses, reflecting their weights and potential for further modifications. Our findings challenge the view that the connectivity among neuronal substrates of memory traces is governed by Hebbian mechanisms, and offer a structural basis for representational drifts.

Coupling of Sharp Wave Events between Zebrafish Hippocampal and Amygdala Homologs

The mammalian hippocampus exhibits spontaneous sharp wave events (1-30 Hz) with an often-present superimposed fast ripple oscillation (120-220 Hz) to form a sharp wave ripple (SWR) complex. During slow-wave sleep or quiet restfulness, SWRs result from the sequential spiking of hippocampal cell assemblies initially activated during learned or imagined experiences. Additional cortical/subcortical areas exhibit SWR events that are coupled to hippocampal SWRs, and studies in mammals suggest that coupling may be critical for the consolidation and recall of specific memories. In the present study, we have examined juvenile male and female zebrafish and show that SWR events are intrinsically generated and maintained within the telencephalon and that their hippocampal homolog, the anterodorsolateral lobe (ADL), exhibits SW events with ∼9% containing an embedded ripple (SWR). Single-cell calcium imaging coupled to local field potential recordings revealed that ∼10% of active cells in the dorsal telencephalon participate in any given SW event. Furthermore, fluctuations in cholinergic tone modulate SW events consistent with mammalian studies. Moreover, the basolateral amygdala (BLA) homolog exhibits SW events with ∼5% containing an embedded ripple. Computing the SW peak coincidence difference between the ADL and BLA showed bidirectional communication. Simultaneous coupling occurred more frequently within the same hemisphere, and in coupled events across hemispheres, the ADL more commonly preceded BLA. Together, these data suggest conserved mechanisms across species by which SW and SWR events are modulated, and memories may be transferred and consolidated through regional coupling.

Top-down modulation in canonical cortical circuits with short-term plasticity

Cortical dynamics and computations are strongly influenced by diverse GABAergic interneurons, including those expressing parvalbumin (PV), somatostatin (SST), and vasoactive intestinal peptide (VIP). Together with excitatory (E) neurons, they form a canonical microcircuit and exhibit counterintuitive nonlinear phenomena. One instance of such phenomena is response reversal, whereby SST neurons show opposite responses to top-down modulation via VIP depending on the presence of bottom-up sensory input, indicating that the network may function in different regimes under different stimulation conditions. Combining analytical and computational approaches, we demonstrate that model networks with multiple interneuron subtypes and experimentally identified short-term plasticity mechanisms can implement response reversal. Surprisingly, despite not directly affecting SST and VIP activity, PV-to-E short-term depression has a decisive impact on SST response reversal. We show how response reversal relates to inhibition stabilization and the paradoxical effect in the presence of several short-term plasticity mechanisms demonstrating that response reversal coincides with a change in the indispensability of SST for network stabilization. In summary, our work suggests a role of short-term plasticity mechanisms in generating nonlinear phenomena in networks with multiple interneuron subtypes and makes several experimentally testable predictions.

Causal connectivity measures for pulse-output network reconstruction: Analysis and applications

The causal connectivity of a network is often inferred to understand network function. It is arguably acknowledged that the inferred causal connectivity relies on the causality measure one applies, and it may differ from the network's underlying structural connectivity. However, the interpretation of causal connectivity remains to be fully clarified, in particular, how causal connectivity depends on causality measures and how causal connectivity relates to structural connectivity. Here, we focus on nonlinear networks with pulse signals as measured output, e.g., neural networks with spike output, and address the above issues based on four commonly utilized causality measures, i.e., time-delayed correlation coefficient, time-delayed mutual information, Granger causality, and transfer entropy. We theoretically show how these causality measures are related to one another when applied to pulse signals. Taking a simulated Hodgkin-Huxley network and a real mouse brain network as two illustrative examples, we further verify the quantitative relations among the four causality measures and demonstrate that the causal connectivity inferred by any of the four well coincides with the underlying network structural connectivity, therefore illustrating a direct link between the causal and structural connectivity. We stress that the structural connectivity of pulse-output networks can be reconstructed pairwise without conditioning on the global information of all other nodes in a network, thus circumventing the curse of dimensionality. Our framework provides a practical and effective approach for pulse-output network reconstruction.

A robust balancing mechanism for spiking neural networks

Dynamical balance of excitation and inhibition is usually invoked to explain the irregular low firing activity observed in the cortex. We propose a robust nonlinear balancing mechanism for a random network of spiking neurons, which works also in the absence of strong external currents. Biologically, the mechanism exploits the plasticity of excitatory-excitatory synapses induced by short-term depression. Mathematically, the nonlinear response of the synaptic activity is the key ingredient responsible for the emergence of a stable balanced regime. Our claim is supported by a simple self-consistent analysis accompanied by extensive simulations performed for increasing network sizes. The observed regime is essentially fluctuation driven and characterized by highly irregular spiking dynamics of all neurons.

Inference of network connectivity from temporally binned spike trains

BACKGROUND

Processing neural activity to reconstruct network connectivity is a central focus of neuroscience, yet the spatiotemporal requisites of biological nervous systems are challenging for current neuronal sensing modalities. Consequently, methods that leverage limited data to successfully infer synaptic connections, predict activity at single unit resolution, and decipher their effect on whole systems, can uncover critical information about neural processing. Despite the emergence of powerful methods for inferring connectivity, network reconstruction based on temporally subsampled data remains insufficiently unexplored.

NEW METHOD

We infer synaptic weights by processing firing rates within variable time bins for a heterogeneous feed-forward network of excitatory, inhibitory, and unconnected units. We assess classification and optimize model parameters for postsynaptic spike train reconstruction. We test our method on a physiological network of leaky integrate-and-fire neurons displaying bursting patterns and assess prediction of postsynaptic activity from microelectrode array data.

RESULTS

Results reveal parameters for improved prediction and performance and suggest that lower resolution data and limited access to neurons can be preferred.

COMPARISON WITH EXISTING METHOD(S)

Recent computational methods demonstrate highly improved reconstruction of connectivity from networks of parallel spike trains by considering spike lag, time-varying firing rates, and other underlying dynamics. However, these methods insufficiently explore temporal subsampling representative of novel data types.

CONCLUSIONS

We provide a framework for reverse engineering neural networks from data with limited temporal quality, describing optimal parameters for each bin size, which can be further improved using non-linear methods and applied to more complicated readouts and connectivity distributions in multiple brain circuits.

Formation and Retrieval of Cell Assemblies in a Biologically Realistic Spiking Neural Network Model of Area CA3 in the Mouse Hippocampus

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

Hierarchical organization of human physical activity

Human physical activity (HPA), a fundamental physiological signal characteristic of bodily motion is of rapidly growing interest in multidisciplinary research. Here we report the existence of hitherto unidentified hierarchical levels in the temporal organization of HPA on the ultradian scale: on the minute's scale, passive periods are followed by activity bursts of similar intensity ('quanta') that are organized into superstructures on the hours- and on the daily scale. The time course of HPA can be considered a stochastic, quasi-binary process, where quanta, assigned to task-oriented actions are organized into work packages on higher levels of hierarchy. In order to grasp the essence of this complex dynamic behaviour, we established a stochastic mathematical model which could reproduce the main statistical features of real activity time series. The results are expected to provide important data for developing novel behavioural models and advancing the diagnostics of neurological or psychiatric diseases.

A generative model of the connectome with dynamic axon growth

Connectome generative models, otherwise known as generative network models, provide insight into the wiring principles underpinning brain network organization. While these models can approximate numerous statistical properties of empirical networks, they typically fail to explicitly characterize an important contributor to brain organization - axonal growth. Emulating the chemoaffinity guided axonal growth, we provide a novel generative model in which axons dynamically steer the direction of propagation based on distance-dependent chemoattractive forces acting on their growth cones. This simple dynamic growth mechanism, despite being solely geometry-dependent, is shown to generate axonal fiber bundles with brain-like geometry and features of complex network architecture consistent with the human brain, including lognormally distributed connectivity weights, scale-free nodal degrees, small-worldness, and modularity. We demonstrate that our model parameters can be fitted to individual connectomes, enabling connectome dimensionality reduction and comparison of parameters between groups. Our work offers an opportunity to bridge studies of axon guidance and connectome development, providing new avenues for understanding neural development from a computational perspective.

Synaptic Changes in Pallidostriatal Circuits Observed in the Parkinsonian Model Triggers Abnormal Beta Synchrony with Accurate Spatio-temporal Properties across the Basal Ganglia

Excessive oscillatory activity across basal ganglia (BG) nuclei in the β frequencies (12-30 Hz) is a hallmark of Parkinson's disease (PD). While the link between oscillations and symptoms remains debated, exaggerated β oscillations constitute an important biomarker for therapeutic effectiveness in PD. The neuronal mechanisms of β-oscillation generation however remain unknown. Many existing models rely on a central role of the subthalamic nucleus (STN) or cortical inputs to BG. Contrarily, neural recordings and optogenetic manipulations in normal and parkinsonian rats recently highlighted the central role of the external pallidum (GPe) in abnormal β oscillations, while showing that the integrity of STN or motor cortex is not required. Here, we evaluate the mechanisms for the generation of abnormal β oscillations in a BG network model where neuronal and synaptic time constants, connectivity, and firing rate distributions are strongly constrained by experimental data. Guided by a mean-field approach, we show in a spiking neural network that several BG sub-circuits can drive oscillations. Strong recurrent STN-GPe connections or collateral intra-GPe connections drive γ oscillations (>40 Hz), whereas strong pallidostriatal loops drive low-β (10-15 Hz) oscillations. We show that pathophysiological strengthening of striatal and pallidal synapses following dopamine depletion leads to the emergence of synchronized oscillatory activity in the mid-β range with spike-phase relationships between BG neuronal populations in-line with experiments. Furthermore, inhibition of GPe, contrary to STN, abolishes oscillations. Our modeling study uncovers the neural mechanisms underlying PD β oscillations and may thereby guide the future development of therapeutic strategies.

Primate neocortex performs balanced sensory amplification

The sensory cortex amplifies relevant features of external stimuli. This sensitivity and selectivity arise through the transformation of inputs by cortical circuitry. We characterize the circuit mechanisms and dynamics of cortical amplification by making large-scale simultaneous measurements of single cells in awake primates and testing computational models. By comparing network activity in both driven and spontaneous states with models, we identify the circuit as operating in a regime of non-normal balanced amplification. Incoming inputs are strongly but transiently amplified by strong recurrent feedback from the disruption of excitatory-inhibitory balance in the network. Strong inhibition rapidly quenches responses, thereby permitting the tracking of time-varying stimuli.

Structure and function of the hippocampal CA3 module

The hippocampal formation is crucial for learning and memory, with submodule CA3 thought to be the substrate of pattern completion. However, the underlying synaptic and computational mechanisms of this network are not well understood. Here, we perform circuit reconstruction of a CA3 module using three dimensional (3D) electron microscopy data and combine this with functional connectivity recordings and computational simulations to determine possible CA3 network mechanisms. Direct measurements of connectivity schemes with both physiological measurements and structural 3D EM revealed a high connectivity rate, multi-fold higher than previously assumed. Mathematical modelling indicated that such CA3 networks can robustly generate pattern completion and replay memory sequences. In conclusion, our data demonstrate that the connectivity scheme of the hippocampal submodule is well suited for efficient memory storage and retrieval.

Reconstruction of phase dynamics from macroscopic observations based on linear and nonlinear response theories

We propose a method to reconstruct the phase dynamics in rhythmical interacting systems from macroscopic responses to weak inputs by developing linear and nonlinear response theories, which predict the responses in a given system. By solving an inverse problem, the method infers an unknown system: the natural frequency distribution, the coupling function, and the time delay which is inevitable in real systems. In contrast to previous methods, our method requires neither strong invasiveness nor microscopic observations. We demonstrate that the method reconstructs two phase systems from observed responses accurately. The qualitative methodological advantages demonstrated by our quantitative numerical examinations suggest its broad applicability in various fields, including brain systems, which are often observed through macroscopic signals such as electroencephalograms and functional magnetic response imaging.

Implications of variable synaptic weights for rate and temporal coding of cerebellar outputs

Purkinje cell (PC) synapses onto cerebellar nuclei (CbN) neurons allow signals from the cerebellar cortex to influence the rest of the brain. PCs are inhibitory neurons that spontaneously fire at high rates, and many PC inputs are thought to converge onto each CbN neuron to suppress its firing. It has been proposed that PCs convey information using a rate code, a synchrony and timing code, or both. The influence of PCs on CbN neuron firing was primarily examined for the combined effects of many PC inputs with comparable strengths, and the influence of individual PC inputs has not been extensively studied. Here, we find that single PC to CbN synapses are highly variable in size, and using dynamic clamp and modeling we reveal that this has important implications for PC-CbN transmission. Individual PC inputs regulate both the rate and timing of CbN firing. Large PC inputs strongly influence CbN firing rates and transiently eliminate CbN firing for several milliseconds. Remarkably, the refractory period of PCs leads to a brief elevation of CbN firing prior to suppression. Thus, individual PC-CbN synapses are suited to concurrently convey rate codes and generate precisely timed responses in CbN neurons. Either synchronous firing or synchronous pauses of PCs promote CbN neuron firing on rapid time scales for nonuniform inputs, but less effectively than for uniform inputs. This is a secondary consequence of variable input sizes elevating the baseline firing rates of CbN neurons by increasing the variability of the inhibitory conductance. These findings may generalize to other brain regions with highly variable inhibitory synapse sizes.

Simulation study on high spatio-temporal resolution acousto-electrophysiological neuroimaging

Objective.Acousto-electrophysiological neuroimaging (AENI) is a technique hypothesized to record electrophysiological activity of the brain with millimeter spatial and sub-millisecond temporal resolution. This improvement is obtained by tagging areas with focused ultrasound (fUS). Due to mechanical vibration with respect to the measuring electrodes, the electrical activity of the marked region will be modulated onto the ultrasonic frequency. The region's electrical activity can subsequently be retrieved via demodulation of the measured signal. In this study, the feasibility of this hypothesized technique is tested.Approach.This is done by calculating the forward electroencephalography response under quasi-static assumptions. The head is simplified as a set of concentric spheres. Two sizes are evaluated representing human and mouse brains. Moreover, feasibility is assessed for wet and dry transcranial, and for cortically placed electrodes. The activity sources are modeled by dipoles, with their current intensity profile drawn from a power-law power spectral density.Results.It is shown that mechanical vibration modulates the endogenous activity onto the ultrasonic frequency. The signal strength depends non-linearly on the alignment between dipole orientation, vibration direction and recording point. The strongest signal is measured when these three dependencies are perfectly aligned. The signal strengths are in the pV-range for a dipole moment of 5 nAm and ultrasonic pressures within Food and Drug Administration (FDA)-limits. The endogenous activity can then be accurately reconstructed via demodulation. Two interference types are investigated: vibrational and static. Depending on the vibrational interference, it is shown that millimeter resolution signal detection is possible also for deep brain regions. Subsequently, successful demodulation depends on the static interference, that at MHz-range has to be sub-picovolt.Significance.Our results show that mechanical vibration is a possible underlying mechanism of acousto-electrophyisological neuroimaging. This paper is a first step towards improved understanding of the conditions under which AENI is feasible.

Mapping general anesthesia states based on electro-encephalogram transition phases

Cortical electro-encephalography (EEG) served as the clinical reference for monitoring unconsciousness during general anesthesia. The existing EEG-based monitors classified general anesthesia states as underdosed, adequate, or overdosed, lacking predictive power due to the absence of transition phases among these states. In response to this limitation, we undertook an analysis of the EEG signal during isoflurane-induced general anesthesia in mice. Adopting a data-driven approach, we applied signal processing techniques to track θ- and δ-band dynamics, along with iso-electric suppressions. Combining this approach with machine learning, we successfully developed an automated algorithm. The findings of our study revealed that the dampening of the δ-band occurred several minutes before the onset of significant iso-electric suppression episodes. Furthermore, a distinct γ-frequency oscillation was observed, persisting for several minutes during the recovery phase subsequent to isoflurane-induced overdose. As a result of our research, we generated a map summarizing multiple brain states and their transitions, offering a tool for predicting and preventing overdose during general anesthesia. The transition phases identified, along with the developed algorithm, have the potential to be generalized, enabling clinicians to prevent inadequate anesthesia and, consequently, tailor anesthetic regimens to individual patients.

Exact Analysis of the Subthreshold Variability for Conductance-Based Neuronal Models with Synchronous Synaptic Inputs

The spiking activity of neocortical neurons exhibits a striking level of variability, even when these networks are driven by identical stimuli. The approximately Poisson firing of neurons has led to the hypothesis that these neural networks operate in the asynchronous state. In the asynchronous state, neurons fire independently from one another, so that the probability that a neuron experience synchronous synaptic inputs is exceedingly low. While the models of asynchronous neurons lead to observed spiking variability, it is not clear whether the asynchronous state can also account for the level of subthreshold membrane potential variability. We propose a new analytical framework to rigorously quantify the subthreshold variability of a single conductance-based neuron in response to synaptic inputs with prescribed degrees of synchrony. Technically, we leverage the theory of exchangeability to model input synchrony via jump-process-based synaptic drives; we then perform a moment analysis of the stationary response of a neuronal model with all-or-none conductances that neglects postspiking reset. As a result, we produce exact, interpretable closed forms for the first two stationary moments of the membrane voltage, with explicit dependence on the input synaptic numbers, strengths, and synchrony. For biophysically relevant parameters, we find that the asynchronous regime yields realistic subthreshold variability (voltage variance ≃4-9 mV2) only when driven by a restricted number of large synapses, compatible with strong thalamic drive. By contrast, we find that achieving realistic subthreshold variability with dense cortico-cortical inputs requires including weak but nonzero input synchrony, consistent with measured pairwise spiking correlations. We also show that, without synchrony, the neural variability averages out to zero for all scaling limits with vanishing synaptic weights, independent of any balanced state hypothesis. This result challenges the theoretical basis for mean-field theories of the asynchronous state.

Explicit mutual information for simple networks and neurons with lognormal activities

Networks with stochastic variables described by heavy-tailed lognormal distribution are ubiquitous in nature, and hence they deserve an exact information-theoretic characterization. We derive analytical formulas for mutual information between elements of different networks with correlated lognormally distributed activities. In a special case, we find an explicit expression for mutual information between neurons when neural activities and synaptic weights are lognormally distributed, as suggested by experimental data. Comparison of this expression with the case when these two variables have short tails reveals that mutual information with heavy tails for neurons and synapses is generally larger and can diverge for some finite variances in presynaptic firing rates and synaptic weights. This result suggests that evolution might prefer brains with heterogeneous dynamics to optimize information processing.

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.

Exact analysis of the subthreshold variability for conductance-based neuronal models with synchronous synaptic inputs

The spiking activity of neocortical neurons exhibits a striking level of variability, even when these networks are driven by identical stimuli. The approximately Poisson firing of neurons has led to the hypothesis that these neural networks operate in the asynchronous state. In the asynchronous state neurons fire independently from one another, so that the probability that a neuron experience synchronous synaptic inputs is exceedingly low. While the models of asynchronous neurons lead to observed spiking variability, it is not clear whether the asynchronous state can also account for the level of subthreshold membrane potential variability. We propose a new analytical framework to rigorously quantify the subthreshold variability of a single conductance-based neuron in response to synaptic inputs with prescribed degrees of synchrony. Technically we leverage the theory of exchangeability to model input synchrony via jump-process-based synaptic drives; we then perform a moment analysis of the stationary response of a neuronal model with all-or-none conductances that neglects post-spiking reset. As a result, we produce exact, interpretable closed forms for the first two stationary moments of the membrane voltage, with explicit dependence on the input synaptic numbers, strengths, and synchrony. For biophysically relevant parameters, we find that the asynchronous regime only yields realistic subthreshold variability (voltage variance ≃ 4 - 9 m V 2 ) when driven by a restricted number of large synapses, compatible with strong thalamic drive. By contrast, we find that achieving realistic subthreshold variability with dense cortico-cortical inputs requires including weak but nonzero input synchrony, consistent with measured pairwise spiking correlations. We also show that without synchrony, the neural variability averages out to zero for all scaling limits with vanishing synaptic weights, independent of any balanced state hypothesis. This result challenges the theoretical basis for mean-field theories of the asynchronous state.