Pynapple, a toolbox for data analysis in neuroscience.

Datasets collected in neuroscientific studies are of ever-growing complexity, often combining high-dimensional time series data from multiple data acquisition modalities. Handling and manipulating these various data streams in an adequate programming environment is crucial to ensure reliable analysis, and to facilitate sharing of reproducible analysis pipelines. Here, we present Pynapple, the PYthon Neural Analysis Package, a lightweight python package designed to process a broad range of time-resolved data in systems neuroscience. The core feature of this package is a small number of versatile objects that support the manipulation of any data streams and task parameters. The package includes a set of methods to read common data formats and allows users to easily write their own. The resulting code is easy to read and write, avoids low-level data processing and other error-prone steps, and is open source. Libraries for higher-level analyses are developed within the Pynapple framework but are contained within a collaborative repository of specialized and continuously updated analysis routines. This provides flexibility while ensuring long-term stability of the core package. In conclusion, Pynapple provides a common framework for data analysis in neuroscience.

State-dependent coupling of hippocampal oscillations.

Oscillations occurring simultaneously in a given area represent a physiological unit of brain states. They allow for temporal segmentation of spikes and support distinct behaviors. To establish how multiple oscillatory components co-vary simultaneously and influence neuronal firing during sleep and wakefulness in mice, we describe a multivariate analytical framework for constructing the state space of hippocampal oscillations. Examining the co-occurrence patterns of oscillations on the state space, across species, uncovered the presence of network constraints and distinct set of cross-frequency interactions during wakefulness compared to sleep. We demonstrated how the state space can be used as a canvas to map the neural firing and found that distinct neurons during navigation were tuned to different sets of simultaneously occurring oscillations during sleep. This multivariate analytical framework provides a window to move beyond classical bivariate pipelines for investigating oscillations and neuronal firing, thereby allowing to factor-in the complexity of oscillation-population interactions.

SpikeShip: A method for fast, unsupervised discovery of high-dimensional neural spiking patterns.

Neural coding and memory formation depend on temporal spiking sequences that span high-dimensional neural ensembles. The unsupervised discovery and characterization of these spiking sequences requires a suitable dissimilarity measure to spiking patterns, which can then be used for clustering and decoding. Here, we present a new dissimilarity measure based on optimal transport theory called SpikeShip, which compares multi-neuron spiking patterns based on all the relative spike-timing relationships among neurons. SpikeShip computes the optimal transport cost to make all the relative spike-timing relationships (across neurons) identical between two spiking patterns. We show that this transport cost can be decomposed into a temporal rigid translation term, which captures global latency shifts, and a vector of neuron-specific transport flows, which reflect inter-neuronal spike timing differences. SpikeShip can be effectively computed for high-dimensional neuronal ensembles, has a low (linear) computational cost that has the same order as the spike count, and is sensitive to higher-order correlations. Furthermore, SpikeShip is binless, can handle any form of spike time distributions, is not affected by firing rate fluctuations, can detect patterns with a low signal-to-noise ratio, and can be effectively combined with a sliding window approach. We compare the advantages and differences between SpikeShip and other measures like SPIKE and Victor-Purpura distance. We applied SpikeShip to large-scale Neuropixel recordings during spontaneous activity and visual encoding. We show that high-dimensional spiking sequences detected via SpikeShip reliably distinguish between different natural images and different behavioral states. These spiking sequences carried complementary information to conventional firing rate codes. SpikeShip opens new avenues for studying neural coding and memory consolidation by rapid and unsupervised detection of temporal spiking patterns in high-dimensional neural ensembles.

Heterogeneity of network and coding states in mouse CA1 place cells.

Theta sequences and phase precession shape hippocampal activity and are considered key underpinnings of memory formation. Theta sequences are sweeps of spikes from multiple cells, tracing trajectories from past to future. Phase precession is the correlation between theta firing phase and animal position. Here, we reconsider these temporal processes in CA1 and the computational principles that they are thought to obey. We find stronger heterogeneity than previously described: we identify cells that do not phase precess but reliably express theta sequences. Other cells phase precess only when medium gamma (linked to entorhinal inputs) is strongest. The same cells express more sequences, but not precession, when slow gamma (linked to CA3 inputs) dominates. Moreover, sequences occur independently in distinct cell groups. Our results challenge the view that phase precession is the mechanism underlying the emergence of theta sequences, suggesting a role for CA1 cells in multiplexing diverse computational processes.

Hippocampal gamma and sharp wave/ripples mediate bidirectional interactions with cortical networks during sleep.

Hippocampus-neocortex interactions during sleep are critical for memory processes: Hippocampally initiated replay contributes to memory consolidation in the neocortex and hippocampal sharp wave/ripples modulate cortical activity. Yet, the spatial and temporal patterns of this interaction are unknown. With voltage imaging, electrocorticography, and laminarly resolved hippocampal potentials, we characterized cortico-hippocampal signaling during anesthesia and nonrapid eye movement sleep. We observed neocortical activation transients, with statistics suggesting a quasi-critical regime, may be helpful for communication across remote brain areas. From activity transients, we identified, in a data-driven fashion, three functional networks. A network overlapping with the default mode network and centered on retrosplenial cortex was the most associated with hippocampal activity. Hippocampal slow gamma rhythms were strongly associated to neocortical transients, even more than ripples. In fact, neocortical activity predicted hippocampal slow gamma and followed ripples, suggesting that consolidation processes rely on bidirectional signaling between hippocampus and neocortex.

Replay, the default mode network and the cascaded memory systems model.

The spontaneous replay of patterns of activity related to past experiences and memories is a striking feature of brain activity, as is the coherent activation of sets of brain areas - particularly those comprising the default mode network (DMN) - during rest. We propose that these two phenomena are strongly intertwined and that their potential functions overlap. In the 'cascaded memory systems model' that we outline here, we hypothesize that the DMN forms the backbone for the propagation of replay, mediating interactions between the hippocampus and the neocortex that enable the consolidation of new memories. The DMN may also independently ignite replay cascades, which support reactivation of older memories or high-level semantic representations. We suggest that transient cortical activations, inducing long-range correlations across the neocortex, are a key mechanism supporting a hierarchy of representations that progresses from simple percepts to semantic representations of causes and, finally, to whole episodes.

The Hybrid Drive: a chronic implant device combining tetrode arrays with silicon probes for layer-resolved ensemble electrophysiology in freely moving mice.

Objective. Understanding the function of brain cortices requires simultaneous investigation at multiple spatial and temporal scales and to link neural activity to an animal's behavior. A major challenge is to measure within- and across-layer information in actively behaving animals, in particular in mice that have become a major species in neuroscience due to an extensive genetic toolkit. Here we describe the Hybrid Drive, a new chronic implant for mice that combines tetrode arrays to record within-layer information with silicon probes to simultaneously measure across-layer information.Approach. The design of our device combines up to 14 tetrodes and 2 silicon probes, that can be arranged in custom arrays to generate unique areas-specific (and multi-area) layouts.Main results. We show that large numbers of neurons and layer-resolved local field potentials can be recorded from the same brain region across weeks without loss in electrophysiological signal quality. The drive's lightweight structure (≈3.5 g) leaves animal behavior largely unchanged, compared to other tetrode drives, during a variety of experimental paradigms. We demonstrate how the data collected with the Hybrid Drive allow state-of-the-art analysis in a series of experiments linking the spiking activity of CA1 pyramidal layer neurons to the oscillatory activity across hippocampal layers.Significance. Our new device fits a gap in the existing technology and increases the range and precision of questions that can be addressed about neural computations in freely behaving mice.

Learning differentially shapes prefrontal and hippocampal activity during classical conditioning.

The ability to use sensory cues to inform goal-directed actions is a critical component of behavior. To study how sounds guide anticipatory licking during classical conditioning, we employed high-density electrophysiological recordings from the hippocampal CA1 area and the prefrontal cortex (PFC) in mice. CA1 and PFC neurons undergo distinct learning-dependent changes at the single-cell level and maintain representations of cue identity at the population level. In addition, reactivation of task-related neuronal assemblies during hippocampal awake Sharp-Wave Ripples (aSWRs) changed within individual sessions in CA1 and over the course of multiple sessions in PFC. Despite both areas being highly engaged and synchronized during the task, we found no evidence for coordinated single cell or assembly activity during conditioning trials or aSWR. Taken together, our findings support the notion that persistent firing and reactivation of task-related neural activity patterns in CA1 and PFC support learning during classical conditioning.

Reduced firing rates of pyramidal cells in the frontal cortex of APP/PS1 can be restored by acute treatment with levetiracetam.

In recent years, aberrant neural oscillations in various cortical areas have emerged as a common physiological hallmark across mouse models of amyloid pathology and patients with Alzheimer's disease. However, much less is known about the underlying effect of amyloid pathology on single cell activity. Here, we used high-density silicon probe recordings from frontal cortex area of 9-month-old APP/PS1 mice to show that local field potential power in the theta and beta band is increased in transgenic animals, whereas single-cell firing rates, specifically of putative pyramidal cells, are significantly reduced. At the same time, these sparsely firing pyramidal cells phase-lock their spiking activity more strongly to the ongoing theta and beta rhythms. Furthermore, we demonstrated that the antiepileptic drug, levetiracetam, counteracts these effects by increasing pyramidal cell firing rates in APP/PS1 mice and uncoupling pyramidal cells and interneurons. Overall, our results highlight reduced firing rates of cortical pyramidal cells as a pathophysiological phenotype in APP/PS1 mice and indicate a potentially beneficial effect of acute levetiracetam treatment.

Unsupervised clustering of temporal patterns in high-dimensional neuronal ensembles using a novel dissimilarity measure.

Temporally ordered multi-neuron patterns likely encode information in the brain. We introduce an unsupervised method, SPOTDisClust (Spike Pattern Optimal Transport Dissimilarity Clustering), for their detection from high-dimensional neural ensembles. SPOTDisClust measures similarity between two ensemble spike patterns by determining the minimum transport cost of transforming their corresponding normalized cross-correlation matrices into each other (SPOTDis). Then, it performs density-based clustering based on the resulting inter-pattern dissimilarity matrix. SPOTDisClust does not require binning and can detect complex patterns (beyond sequential activation) even when high levels of out-of-pattern "noise" spiking are present. Our method handles efficiently the additional information from increasingly large neuronal ensembles and can detect a number of patterns that far exceeds the number of recorded neurons. In an application to neural ensemble data from macaque monkey V1 cortex, SPOTDisClust can identify different moving stimulus directions on the sole basis of temporal spiking patterns.

Impaired hippocampal representation of place in the Fmr1-knockout mouse model of fragile X syndrome.

Fragile X syndrome (FXS) is an X-chromosome linked intellectual disability and the most common known inherited single gene cause of autism spectrum disorder (ASD). Building upon demonstrated deficits in neuronal plasticity and spatial memory in FXS, we investigated how spatial information processing is affected in vivo in an FXS mouse model (Fmr1-KO). Healthy hippocampal neurons (so-called place cells) exhibit place-related activity during spatial exploration, and their firing fields tend to remain stable over time. In contrast, we find impaired stability and reduced specificity of Fmr1-KO spatial representations. This is a potential biomarker for the cognitive dysfunction observed in FXS, informative on the ability to integrate sensory information into an abstract representation and successfully retain this conceptual memory. Our results provide key insight into the biological mechanisms underlying cognitive disabilities in FXS and ASD, paving the way for a targeted approach to remedy these.

Abnormal hippocampal theta and gamma hypersynchrony produces network and spike timing disturbances in the Fmr1-KO mouse model of Fragile X syndrome.

Neuronal networks can synchronize their activity through excitatory and inhibitory connections, which is conducive to synaptic plasticity. This synchronization is reflected in rhythmic fluctuations of the extracellular field. In the hippocampus, theta and gamma band LFP oscillations are a hallmark of the processing of spatial information and memory. Fragile X syndrome (FXS) is an intellectual disability and the most common genetic cause of autism spectrum disorder (Belmonte and Bourgeron, 2006). Here, we investigated how neuronal network synchronization in the mouse hippocampus is compromised by the Fmr1 mutation that causes FXS (Santos et al., 2014), relating recently observed single-cell level impairments (Arbab et al., 2017) to neuronal network aberrations. We implanted tetrodes in hippocampus of freely moving Fmr1-KO and littermate wildtype (WT) mice (Mientjes et al., 2006), to record spike trains from multiple, isolated neurons as well as LFPs in a spatial exploration paradigm. Compared to wild type mice, Fmr1-KO mice displayed greater power of hippocampal theta oscillations, and higher coherence in the slow gamma band. Additionally, spike trains of Fmr1-KO interneurons show decreased spike-count correlations and they are hypersynchronized with theta and slow gamma oscillations. The hypersynchronization of Fmr1-KO oscillations and spike timing reflects functional deficits in local networks. This network hypersynchronization pathologically decreases the heterogeneity of spike-LFP phase coupling, compromising information processing within the hippocampal circuit. These findings may reflect a pathophysiological mechanism explaining cognitive impairments in FXS and autism, in which there is anomalous processing of social and environmental cues and associated deficits in memory and cognition.

The Yin and Yang of Memory Consolidation: Hippocampal and Neocortical.

While hippocampal and cortical mechanisms of memory consolidation have long been studied, their interaction is poorly understood. We sought to investigate potential interactions with respect to trace dominance, strengthening, and interference associated with postencoding novelty or sleep. A learning procedure was scheduled in a watermaze that placed the impact of novelty and sleep in opposition. Distinct behavioural manipulations-context preexposure or interference during memory retrieval-differentially affected trace dominance and trace survival, respectively. Analysis of immediate early gene expression revealed parallel up-regulation in the hippocampus and cortex, sustained in the hippocampus in association with novelty but in the cortex in association with sleep. These findings shed light on dynamically interacting mechanisms mediating the stabilization of hippocampal and neocortical memory traces. Hippocampal memory traces followed by novelty were more dominant by default but liable to interference, whereas sleep engaged a lasting stabilization of cortical traces and consequent trace dominance after preexposure.

Functional coupling networks inferred from prefrontal cortex activity show experience-related effective plasticity.

Functional coupling networks are widely used to characterize collective patterns of activity in neural populations. Here, we ask whether functional couplings reflect the subtle changes, such as in physiological interactions, believed to take place during learning. We infer functional network models reproducing the spiking activity of simultaneously recorded neurons in prefrontal cortex (PFC) of rats, during the performance of a cross-modal rule shift task (task epoch), and during preceding and following sleep epochs. A large-scale study of the 96 recorded sessions allows us to detect, in about 20% of sessions, effective plasticity between the sleep epochs. These coupling modifications are correlated with the coupling values in the task epoch, and are supported by a small subset of the recorded neurons, which we identify by means of an automatized procedure. These potentiated groups increase their coativation frequency in the spiking data between the two sleep epochs, and, hence, participate to putative experience-related cell assemblies. Study of the reactivation dynamics of the potentiated groups suggests a possible connection with behavioral learning. Reactivation is largely driven by hippocampal ripple events when the rule is not yet learned, and may be much more autonomous, and presumably sustained by the potentiated PFC network, when learning is consolidated.

Effects of Arc/Arg3.1 gene deletion on rhythmic synchronization of hippocampal CA1 neurons during locomotor activity and sleep.

The activity-regulated cytoskeletal-associated protein/activity regulated gene (Arc/Arg3.1) is crucial for long-term synaptic plasticity and memory formation. However, the neurophysiological substrates of memory deficits occurring in the absence of Arc/Arg3.1 are unknown. We compared hippocampal CA1 single-unit and local field potential (LFP) activity in Arc/Arg3.1 knockout and wild-type mice during track running and flanking sleep periods. Locomotor activity, basic firing and spatial coding properties of CA1 cells in knockout mice were not different from wild-type mice. During active behavior, however, knockout animals showed a significantly shifted balance in LFP power, with a relative loss in high-frequency (beta-2 and gamma) bands compared to low-frequency bands. Moreover, during track-running, knockout mice showed a decrease in phase locking of spiking activity to LFP oscillations in theta, beta and gamma bands. Sleep architecture in knockout mice was not grossly abnormal. Sharp-wave ripples, which have been associated with memory consolidation and replay, showed only minor differences in dynamics and amplitude. Altogether, these findings suggest that Arc/Arg3.1 effects on memory formation are not only manifested at the level of molecular pathways regulating synaptic plasticity, but also at the systems level. The disrupted power balance in theta, beta and gamma rhythmicity and concomitant loss of spike-field phase locking may affect memory encoding during initial storage and memory consolidation stages.

A Multichannel Recording System with Optical Stimulation for Closed-Loop Optogenetic Experiments.

Selective perturbation of the activity of specific cell types in the brain tissue is essential in understanding the function of neuronal circuits involved in cognition and behavior and might also provide therapeutic neuromodulation strategies. Such selective neuronal addressing can be achieved through the optical activation of light-sensitive proteins called opsins that are expressed in specific cell populations through genetic methods-hence the name "optogenetics." In optogenetic experiments, the electrical activity of the targeted cell populations is optically triggered and monitored using arrays of microelectrodes. In closed-loop studies, the optical stimulation parameters are adjusted based on the recorded activity, ideally in real time. Here we describe the basic tools and the protocols allowing closed-loop optogenic experiments in vivo.

CA2: it's about time-and episodes.

In this issue of Neuron, Mankin et al. (2015) show that CA2, an oft-neglected hippocampal subregion, has place representations that change from one episode to the next, even as the spatial environment does not. This finding may help explain how time is encoded in episodic memories.

How to detect the Granger-causal flow direction in the presence of additive noise?

Granger-causality metrics have become increasingly popular tools to identify directed interactions between brain areas. However, it is known that additive noise can strongly affect Granger-causality metrics, which can lead to spurious conclusions about neuronal interactions. To solve this problem, previous studies have proposed the detection of Granger-causal directionality, i.e. the dominant Granger-causal flow, using either the slope of the coherency (Phase Slope Index; PSI), or by comparing Granger-causality values between original and time-reversed signals (reversed Granger testing). We show that for ensembles of vector autoregressive (VAR) models encompassing bidirectionally coupled sources, these alternative methods do not correctly measure Granger-causal directionality for a substantial fraction of VAR models, even in the absence of noise. We then demonstrate that uncorrelated noise has fundamentally different effects on directed connectivity metrics than linearly mixed noise, where the latter may result as a consequence of electric volume conduction. Uncorrelated noise only weakly affects the detection of Granger-causal directionality, whereas linearly mixed noise causes a large fraction of false positives for standard Granger-causality metrics and PSI, but not for reversed Granger testing. We further show that we can reliably identify cases where linearly mixed noise causes a large fraction of false positives by examining the magnitude of the instantaneous influence coefficient in a structural VAR model. By rejecting cases with strong instantaneous influence, we obtain an improved detection of Granger-causal flow between neuronal sources in the presence of additive noise. These techniques are applicable to real data, which we demonstrate using actual area V1 and area V4 LFP data, recorded from the awake monkey performing a visual attention task.