Explicit-duration hidden Markov model inference of UP-DOWN states from continuous signals

Spontaneous persistent activity and inactivity in vivo reveals differential cortico-entorhinal functional connectivity

Understanding the functional connectivity between brain regions and its emergent dynamics is a central challenge. Here we present a theory-experiment hybrid approach involving iteration between a minimal computational model and in vivo electrophysiological measurements. Our model not only predicted spontaneous persistent activity (SPA) during Up-Down-State oscillations, but also inactivity (SPI), which has never been reported. These were confirmed in vivo in the membrane potential of neurons, especially from layer 3 of the medial and lateral entorhinal cortices. The data was then used to constrain two free parameters, yielding a unique, experimentally determined model for each neuron. Analytic and computational analysis of the model generated a dozen quantitative predictions about network dynamics, which were all confirmed in vivo to high accuracy. Our technique predicted functional connectivity; e. g. the recurrent excitation is stronger in the medial than lateral entorhinal cortex. This too was confirmed with connectomics data. This technique uncovers how differential cortico-entorhinal dialogue generates SPA and SPI, which could form an energetically efficient working-memory substrate and influence the consolidation of memories during sleep. More broadly, our procedure can reveal the functional connectivity of large networks and a theory of their emergent dynamics.

Detection of neuronal OFF periods as low amplitude neural activity segments

BACKGROUND

During non-rapid eye movement sleep (NREM), alternating periods of synchronised high (ON period) and low (OFF period) neuronal activity are associated with high amplitude delta band (0.5-4 Hz) oscillations in neocortical electrophysiological signals termed slow waves. As this oscillation is dependent crucially on hyperpolarisation of cortical cells, there is an interest in understanding how neuronal silencing during OFF periods leads to the generation of slow waves and whether this relationship changes between cortical layers. A formal, widely adopted definition of OFF periods is absent, complicating their detection. Here, we grouped segments of high frequency neural activity containing spikes, recorded as multiunit activity from the neocortex of freely behaving mice, on the basis of amplitude and asked whether the population of low amplitude (LA) segments displayed the expected characteristics of OFF periods.

RESULTS

Average LA segment length was comparable to previous reports for OFF periods but varied considerably, from as short as 8 ms to > 1 s. LA segments were longer and occurred more frequently in NREM but shorter LA segments also occurred in half of rapid eye movement sleep (REM) epochs and occasionally during wakefulness. LA segments in all states were associated with a local field potential (LFP) slow wave that increased in amplitude with LA segment duration. We found that LA segments > 50 ms displayed a homeostatic rebound in incidence following sleep deprivation whereas short LA segments (< 50 ms) did not. The temporal organisation of LA segments was more coherent between channels located at a similar cortical depth.

CONCLUSION

We corroborate previous studies showing neural activity signals contain uniquely identifiable periods of low amplitude with distinct characteristics from the surrounding signal known as OFF periods and attribute the new characteristics of vigilance-state-dependent duration and duration-dependent homeostatic response to this phenomenon. This suggests that ON/OFF periods are currently underdefined and that their appearance is less binary than previously considered, instead representing a continuum.

UP-DOWN states and ripples differentially modulate membrane potential dynamics across DG, CA3, and CA1 in awake mice

Hippocampal ripples are transient population bursts that structure cortico-hippocampal communication and play a central role in memory processing. However, the mechanisms controlling ripple initiation in behaving animals remain poorly understood. Here we combine multisite extracellular and whole-cell recordings in awake mice to contrast the brain state and ripple modulation of subthreshold dynamics across hippocampal subfields. We find that entorhinal input to the dentate gyrus (DG) exhibits UP and DOWN dynamics with ripples occurring exclusively in UP states. While elevated cortical input in UP states generates depolarization in DG and CA1, it produces persistent hyperpolarization in CA3 neurons. Furthermore, growing inhibition is evident in CA3 throughout the course of the ripple buildup, while DG and CA1 neurons exhibit depolarization transients 100 ms before and during ripples. These observations highlight the importance of CA3 inhibition for ripple generation, while pre-ripple responses indicate a long and orchestrated ripple initiation process in the awake state.

A mechanism for inter-areal coherence through communication based on connectivity and oscillatory power

Inter-areal coherence between field potentials is a widespread phenomenon in cortex. Coherence has been hypothesized to reflect phase-synchronization between oscillators and flexibly gate communication according to behavioral and cognitive demands. We reveal an alternative mechanism where coherence is not the cause but the consequence of communication and naturally emerges because spiking activity in a sending area causes post-synaptic potentials both in the same and in other areas. Consequently, coherence depends in a lawful manner on power and phase-locking in the sender and connectivity. Changes in oscillatory power explained prominent changes in fronto-parietal and LGN-V1 coherence across behavioral conditions. Optogenetic experiments and excitatory-inhibitory network simulations identified afferent synaptic inputs rather than spiking entrainment as the principal determinant of coherence. These findings suggest that unique spectral profiles of different brain areas automatically give rise to large-scale coherence patterns that follow anatomical connectivity and continuously reconfigure as a function of behavior and cognition.

A hidden Markov model reliably characterizes ketamine-induced spectral dynamics in macaque local field potentials and human electroencephalograms

Ketamine is an NMDA receptor antagonist commonly used to maintain general anesthesia. At anesthetic doses, ketamine causes high power gamma (25-50 Hz) oscillations alternating with slow-delta (0.1-4 Hz) oscillations. These dynamics are readily observed in local field potentials (LFPs) of non-human primates (NHPs) and electroencephalogram (EEG) recordings from human subjects. However, a detailed statistical analysis of these dynamics has not been reported. We characterize ketamine's neural dynamics using a hidden Markov model (HMM). The HMM observations are sequences of spectral power in seven canonical frequency bands between 0 to 50 Hz, where power is averaged within each band and scaled between 0 and 1. We model the observations as realizations of multivariate beta probability distributions that depend on a discrete-valued latent state process whose state transitions obey Markov dynamics. Using an expectation-maximization algorithm, we fit this beta-HMM to LFP recordings from 2 NHPs, and separately, to EEG recordings from 9 human subjects who received anesthetic doses of ketamine. Our beta-HMM framework provides a useful tool for experimental data analysis. Together, the estimated beta-HMM parameters and optimal state trajectory revealed an alternating pattern of states characterized primarily by gamma and slow-delta activities. The mean duration of the gamma activity was 2.2s([1.7,2.8]s) and 1.2s([0.9,1.5]s) for the two NHPs, and 2.5s([1.7,3.6]s) for the human subjects. The mean duration of the slow-delta activity was 1.6s([1.2,2.0]s) and 1.0s([0.8,1.2]s) for the two NHPs, and 1.8s([1.3,2.4]s) for the human subjects. Our characterizations of the alternating gamma slow-delta activities revealed five sub-states that show regular sequential transitions. These quantitative insights can inform the development of rhythm-generating neuronal circuit models that give mechanistic insights into this phenomenon and how ketamine produces altered states of arousal.

Memory reactivation in rat medial prefrontal cortex occurs in a subtype of cortical UP state during slow-wave sleep

Interaction between hippocampal sharp-wave ripples (SWRs) and UP states, possibly by coordinated reactivation of memory traces, is conjectured to play an important role in memory consolidation. Recently, it was reported that SWRs were differentiated into multiple subtypes. However, whether cortical UP states can also be classified into subtypes is not known. Here, we analysed neural ensemble activity from the medial prefrontal cortex from rats trained to run a spatial sequence-memory task. Application of the hidden Markov model (HMM) with three states to epochs of UP-DOWN oscillations identified DOWN states and two subtypes of UP state (UP-1 and UP-2). The two UP subtypes were distinguished by differences in duration, with UP-1 having a longer duration than UP-2, as well as differences in the speed of population vector (PV) decorrelation, with UP-1 decorrelating more slowly than UP-2. Reactivation of recent memory sequences predominantly occurred in UP-2. Short-duration reactivating UP states were dominated by UP-2 whereas long-duration ones exhibit transitions from UP-1 to UP-2. Thus, recent memory reactivation, if it occurred within long-duration UP states, typically was preceded by a period of slow PV evolution not related to recent experience, and which we speculate may be related to previously encoded information. If that is the case, then the transition from UP-1 to UP-2 subtypes may help gradual integration of recent experience with pre-existing cortical memories by interleaving the two in the same UP state. This article is part of the Theo Murphy meeting issue 'Memory reactivation: replaying events past, present and future'.

Effect of Stimulus-Dependent Spike Timing on Population Coding of Sound Location in the Owl's Auditory Midbrain

In the auditory system, the spectrotemporal structure of acoustic signals determines the temporal pattern of spikes. Here, we investigated this effect in neurons of the barn owl's auditory midbrain (Tyto furcata) that are selective for auditory space and whether it can influence the coding of sound direction. We found that in the nucleus where neurons first become selective to combinations of sound localization cues, reproducibility of spike trains across repeated trials of identical sounds, a metric of across-trial temporal fidelity of spiking patterns evoked by a stimulus, was maximal at the sound direction that elicited the highest firing rate. We then tested the hypothesis that this stimulus-dependent patterning resulted in rate co-modulation of cells with similar frequency and spatial selectivity, driving stimulus-dependent synchrony of population responses. Tetrodes were used to simultaneously record multiple nearby units in the optic tectum (OT), where auditory space is topographically represented. While spiking of neurons in OT showed lower reproducibility across trials compared with upstream nuclei, spike-time synchrony between nearby OT neurons was highest for sounds at their preferred direction. A model of the midbrain circuit explained the relationship between stimulus-dependent reproducibility and synchrony, and demonstrated that this effect can improve the decoding of sound location from the OT output. Thus, stimulus-dependent spiking patterns in the auditory midbrain can have an effect on spatial coding. This study reports a functional connection between spike patterning elicited by spectrotemporal features of a sound and the coding of its location.

Ventral Tegmental Area Dopamine Neurons Firing Model Reveals Prenatal Nicotine Induced Alterations

The dopamine (DA) neurons found in the ventral tegmental area (VTA) are widely involved in the addiction and natural reward circuitry of the brain. Their firing patterns were shown to be important modulators of dopamine release and repetitive burst-like firing activity was highlighted as a major firing pattern of DA neurons in the VTA. In the present study we use a state space model to characterize the DA neurons firing patterns, and trace transitions of neural activity through bursting and non-bursting states. The hidden semi-Markov model (HSMM) framework, which we use, offers a statistically principled inference of bursting states and considers VTA DA firing patterns to be generated according to a Gamma process. Additionally, the explicit Gamma-based modeling of state durations allows efficient decoding of underlying neural information. Consequently, we decode and segment our single unit recordings from DA neurons in VTA according to the sequence of statistically discriminated HSMM states. The segmentation is used to study bursting state characteristics in data recorded from rats prenatally exposed to nicotine (6 mg/kg/day starting with gestational day 3) and rats from saline treated dams. Our results indicate that prenatal nicotine exposure significantly alters burst firing patterns of a subset of DA neurons in adolescent rats, suggesting nicotine exposure during gestation may induce severe effects on the neural networks involved in addiction and reward.

Synaptic dynamics and neuronal network connectivity are reflected in the distribution of times in Up states

The dynamics of neuronal networks connected by synaptic dynamics can sustain long periods of depolarization that can last for hundreds of milliseconds such as Up states recorded during sleep or anesthesia. Yet the underlying mechanism driving these periods remain unclear. We show here within a mean-field model that the residence time of the neuronal membrane potential in cortical Up states does not follow a Poissonian law, but presents several peaks. Furthermore, the present modeling approach allows extracting some information about the neuronal network connectivity from the time distribution histogram. Based on a synaptic-depression model, we find that these peaks, that can be observed in histograms of patch-clamp recordings are not artifacts of electrophysiological measurements, but rather are an inherent property of the network dynamics. Analysis of the equations reveals a stable focus located close to the unstable limit cycle, delimiting a region that defines the Up state. The model further shows that the peaks observed in the Up state time distribution are due to winding around the focus before escaping from the basin of attraction. Finally, we use in vivo recordings of intracellular membrane potential and we recover from the peak distribution, some information about the network connectivity. We conclude that it is possible to recover the network connectivity from the distribution of times that the neuronal membrane voltage spends in Up states.

Estimating latent attentional states based on simultaneous binary and continuous behavioral measures

Cognition is a complex and dynamic process. It is an essential goal to estimate latent attentional states based on behavioral measures in many sequences of behavioral tasks. Here, we propose a probabilistic modeling and inference framework for estimating the attentional state using simultaneous binary and continuous behavioral measures. The proposed model extends the standard hidden Markov model (HMM) by explicitly modeling the state duration distribution, which yields a special example of the hidden semi-Markov model (HSMM). We validate our methods using computer simulations and experimental data. In computer simulations, we systematically investigate the impacts of model mismatch and the latency distribution. For the experimental data collected from a rodent visual detection task, we validate the results with predictive log-likelihood. Our work is useful for many behavioral neuroscience experiments, where the common goal is to infer the discrete (binary or multinomial) state sequences from multiple behavioral measures.

Spontaneous persistent activity in entorhinal cortex modulates cortico-hippocampal interaction in vivo

Persistent activity is thought to mediate working memory during behavior. Can it also occur during sleep? We found that the membrane potential of medial entorhinal cortex layer III (MECIII) neurons, a gateway between neocortex and hippocampus, showed spontaneous, stochastic persistent activity in vivo in mice during Up-Down state oscillations (UDS). This persistent activity was locked to the neocortical Up states with a short delay, but persisted over several cortical UDS cycles. Lateral entorhinal neurons did not show substantial persistence, and current injections similar to those used in vitro failed to elicit persistence in vivo, implicating network mechanisms. Hippocampal CA1 neurons' spiking activity was reduced during neocortical Up states, but was increased during MECIII persistent states. These results provide, to the best of our knowledge, the first direct evidence for persistent activity in MECIII neurons in vivo and reveal its contribution to cortico-hippocampal interaction that could be involved in working memory and learning of long behavioral sequences during behavior, and memory consolidation during sleep.

Nonlinear-dynamics theory of up-down transitions in neocortical neural networks

The neurons of the neocortex show ~1-Hz synchronized transitions between an active up state and a quiescent down state. The up-down state transitions are highly coherent over large sections of the cortex, yet they are accompanied by pronounced, incoherent noise. We propose a simple model for the up-down state oscillations that allows analysis by straightforward dynamical systems theory. An essential feature is a nonuniform network geometry composed of groups of excitatory and inhibitory neurons with strong coupling inside a group and weak coupling between groups. The enhanced deterministic noise of the up state appears as the natural result of the proximity of a partial synchronization transition. The synchronization transition takes place as a function of the long-range synaptic strength linking different groups of neurons.