Neural manifold analysis of brain circuit dynamics in health and disease.

Recent developments in experimental neuroscience make it possible to simultaneously record the activity of thousands of neurons. However, the development of analysis approaches for such large-scale neural recordings have been slower than those applicable to single-cell experiments. One approach that has gained recent popularity is neural manifold learning. This approach takes advantage of the fact that often, even though neural datasets may be very high dimensional, the dynamics of neural activity tends to traverse a much lower-dimensional space. The topological structures formed by these low-dimensional neural subspaces are referred to as "neural manifolds", and may potentially provide insight linking neural circuit dynamics with cognitive function and behavioral performance. In this paper we review a number of linear and non-linear approaches to neural manifold learning, including principal component analysis (PCA), multi-dimensional scaling (MDS), Isomap, locally linear embedding (LLE), Laplacian eigenmaps (LEM), t-SNE, and uniform manifold approximation and projection (UMAP). We outline these methods under a common mathematical nomenclature, and compare their advantages and disadvantages with respect to their use for neural data analysis. We apply them to a number of datasets from published literature, comparing the manifolds that result from their application to hippocampal place cells, motor cortical neurons during a reaching task, and prefrontal cortical neurons during a multi-behavior task. We find that in many circumstances linear algorithms produce similar results to non-linear methods, although in particular cases where the behavioral complexity is greater, non-linear methods tend to find lower-dimensional manifolds, at the possible expense of interpretability. We demonstrate that these methods are applicable to the study of neurological disorders through simulation of a mouse model of Alzheimer's Disease, and speculate that neural manifold analysis may help us to understand the circuit-level consequences of molecular and cellular neuropathology.

Neuronal signature of spatial decision-making during navigation by freely moving rats by using calcium imaging.

A challenge in spatial memory is understanding how place cell firing contributes to decision-making in navigation. A spatial recency task was created in which freely moving rats first became familiar with a spatial context over several days and thereafter were required to encode and then selectively recall one of three specific locations within it that was chosen to be rewarded that day. Calcium imaging was used to record from more than 1,000 cells in area CA1 of the hippocampus of five rats during the exploration, sample, and choice phases of the daily task. The key finding was that neural activity in the startbox rose steadily in the short period prior to entry to the arena and that this selective population cell firing was predictive of the daily changing goal on correct trials but not on trials in which the animals made errors. Single-cell and population activity measures converged on the idea that prospective coding of neural activity can be involved in navigational decision-making.

Refinements to rodent head fixation and fluid/food control for neuroscience.

The use of head fixation in mice is increasingly common in research, its use having initially been restricted to the field of sensory neuroscience. Head restraint has often been combined with fluid control, rather than food restriction, to motivate behaviour, but this too is now in use for both restrained and non-restrained animals. Despite this, there is little guidance on how best to employ these techniques to optimise both scientific outcomes and animal welfare. This article summarises current practices and provides recommendations to improve animal wellbeing and data quality, based on a survey of the community, literature reviews, and the expert opinion and practical experience of an international working group convened by the UK's National Centre for the Replacement, Refinement and Reduction of Animals in Research (NC3Rs). Topics covered include head fixation surgery and post-operative care, habituation to restraint, and the use of fluid/food control to motivate performance. We also discuss some recent developments that may offer alternative ways to collect data from large numbers of behavioural trials without the need for restraint. The aim is to provide support for researchers at all levels, animal care staff, and ethics committees to refine procedures and practices in line with the refinement principle of the 3Rs.

Ca2+ imaging of self and other in medial prefrontal cortex during social dominance interactions in a tube test.

The study of social dominance interactions between animals offers a window onto the decision-making involved in establishing dominance hierarchies and an opportunity to examine changes in social behavior observed in certain neurogenetic disorders. Competitive social interactions, such as in the widely used tube test, reflect this decision-making. Previous studies have focused on the different patterns of behavior seen in the dominant and submissive animal, neural correlates of effortful behavior believed to mediate the outcome of such encounters, and interbrain correlations of neural activity. Using a rigorous mutual information criterion, we now report that neural responses recorded with endoscopic calcium imaging in the prelimbic zone of the medial prefrontal cortex show unique correlations to specific dominance-related behaviors. Interanimal analyses revealed cell/behavior correlations that are primarily with an animal's own behavior or with the other animal's behavior, or the coincident behavior of both animals (such as pushing by one and resisting by the other). The comparison of unique and coincident cells helps to disentangle cell firing that reflects an animal's own or the other's specific behavior from situations reflecting conjoint action. These correlates point to a more cognitive rather than a solely behavioral dimension of social interactions that needs to be considered in the design of neurobiological studies of social behavior. These could prove useful in studies of disorders affecting social recognition and social engagement, and the treatment of disorders of social interaction.

A computational grid-to-place-cell transformation model indicates a synaptic driver of place cell impairment in early-stage Alzheimer's Disease.

Alzheimer's Disease (AD) is characterized by progressive neurodegeneration and cognitive impairment. Synaptic dysfunction is an established early symptom, which correlates strongly with cognitive decline, and is hypothesised to mediate the diverse neuronal network abnormalities observed in AD. However, how synaptic dysfunction contributes to network pathology and cognitive impairment in AD remains elusive. Here, we present a grid-cell-to-place-cell transformation model of long-term CA1 place cell dynamics to interrogate the effect of synaptic loss on network function and environmental representation. Synapse loss modelled after experimental observations in the APP/PS1 mouse model was found to induce firing rate alterations and place cell abnormalities that have previously been observed in AD mouse models, including enlarged place fields and lower across-session stability of place fields. Our results support the hypothesis that synaptic dysfunction underlies cognitive deficits, and demonstrate how impaired environmental representation may arise in the early stages of AD. We further propose that dysfunction of excitatory and inhibitory inputs to CA1 pyramidal cells may cause distinct impairments in place cell function, namely reduced stability and place map resolution.

Place Cells in Head-Fixed Mice Navigating a Floating Real-World Environment.

The hippocampal place cell system in rodents has provided a major paradigm for the scientific investigation of memory function and dysfunction. Place cells have been observed in area CA1 of the hippocampus of both freely moving animals, and of head-fixed animals navigating in virtual reality environments. However, spatial coding in virtual reality preparations has been observed to be impaired. Here we show that the use of a real-world environment system for head-fixed mice, consisting of an air-floating track with proximal cues, provides some advantages over virtual reality systems for the study of spatial memory. We imaged the hippocampus of head-fixed mice injected with the genetically encoded calcium indicator GCaMP6s while they navigated circularly constrained or open environments on the floating platform. We observed consistent place tuning in a substantial fraction of cells despite the absence of distal visual cues. Place fields remapped when animals entered a different environment. When animals re-entered the same environment, place fields typically remapped over a time period of multiple days, faster than in freely moving preparations, but comparable with virtual reality. Spatial information rates were within the range observed in freely moving mice. Manifold analysis indicated that spatial information could be extracted from a low-dimensional subspace of the neural population dynamics. This is the first demonstration of place cells in head-fixed mice navigating on an air-lifted real-world platform, validating its use for the study of brain circuits involved in memory and affected by neurodegenerative disorders.

Dual midbrain and forebrain origins of thalamic inhibitory interneurons.

The ubiquitous presence of inhibitory interneurons in the thalamus of primates contrasts with the sparsity of interneurons reported in mice. Here, we identify a larger than expected complexity and distribution of interneurons across the mouse thalamus, where all thalamic interneurons can be traced back to two developmental programmes: one specified in the midbrain and the other in the forebrain. Interneurons migrate to functionally distinct thalamocortical nuclei depending on their origin: the abundant, midbrain-derived class populates the first and higher order sensory thalamus while the rarer, forebrain-generated class is restricted to some higher order associative regions. We also observe that markers for the midbrain-born class are abundantly expressed throughout the thalamus of the New World monkey marmoset. These data therefore reveal that, despite the broad variability in interneuron density across mammalian species, the blueprint of the ontogenetic organisation of thalamic interneurons of larger-brained mammals exists and can be studied in mice.

Integrating new memories into the hippocampal network activity space.

By investigating the topology of neuronal co-activity, we found that mnemonic information spans multiple operational axes in the mouse hippocampus network. High-activity principal cells form the core of each memory along a first axis, segregating spatial contexts and novelty. Low-activity cells join co-activity motifs across behavioral events and enable their crosstalk along two other axes. This reveals an organizational principle for continuous integration and interaction of hippocampal memories.

Bladder pressure encoding by sacral dorsal root ganglion fibres: implications for decoding.

Objective.We aim at characterising the encoding of bladder pressure (intravesical pressure) by a population of sensory fibres. This research is motivated by the possibility to restore bladder function in elderly patients or after spinal cord injury using implanted devices, so called bioelectronic medicines. For these devices, nerve-based estimation of intravesical pressure can enable a personalized and on-demand stimulation paradigm, which has promise of being more effective and efficient. In this context, a better understanding of the encoding strategies employed by the body might in the future be exploited by informed decoding algorithms that enable a precise and robust bladder-pressure estimation.Approach.To this end, we apply information theory to microelectrode-array recordings from the cat sacral dorsal root ganglion while filling the bladder, conduct surrogate data studies to augment the data we have, and finally decode pressure in a simple informed approach.Main results.We find an encoding scheme by different main bladder neuron types that we divide into three response types (slow tonic, phasic, and derivative fibres). We show that an encoding by different bladder neuron types, each represented by multiple cells, offers reliability through within-type redundancy and high information rates through semi-independence of different types. Our subsequent decoding study shows a more robust decoding from mean responses of homogeneous cell pools.Significance.We have here, for the first time, established a link between an information theoretic analysis of the encoding of intravesical pressure by a population of sensory neurons to an informed decoding paradigm. We show that even a simple adapted decoder can exploit the redundancy in the population to be more robust against cell loss. This work thus paves the way towards principled encoding studies in the periphery and towards a new generation of informed peripheral nerve decoders for bioelectronic medicines.

Subcellular resolution three-dimensional light-field imaging with genetically encoded voltage indicators.

Significance: Light-field microscopy (LFM) enables high signal-to-noise ratio (SNR) and light efficient volume imaging at fast frame rates. Voltage imaging with genetically encoded voltage indicators (GEVIs) stands to particularly benefit from LFM's volumetric imaging capability due to high required sampling rates and limited probe brightness and functional sensitivity. Aim: We demonstrate subcellular resolution GEVI light-field imaging in acute mouse brain slices resolving dendritic voltage signals in three spatial dimensions. Approach: We imaged action potential-induced fluorescence transients in mouse brain slices sparsely expressing the GEVI VSFP-Butterfly 1.2 in wide-field microscopy (WFM) and LFM modes. We compared functional signal SNR and localization between different LFM reconstruction approaches and between LFM and WFM. Results: LFM enabled three-dimensional (3-D) localization of action potential-induced fluorescence transients in neuronal somata and dendrites. Nonregularized deconvolution decreased SNR with increased iteration number compared to synthetic refocusing but increased axial and lateral signal localization. SNR was unaffected for LFM compared to WFM. Conclusions: LFM enables 3-D localization of fluorescence transients, therefore eliminating the need for structures to lie in a single focal plane. These results demonstrate LFM's potential for studying dendritic integration and action potential propagation in three spatial dimensions.

Learning to select actions shapes recurrent dynamics in the corticostriatal system.

Learning to select appropriate actions based on their values is fundamental to adaptive behavior. This form of learning is supported by fronto-striatal systems. The dorsal-lateral prefrontal cortex (dlPFC) and the dorsal striatum (dSTR), which are strongly interconnected, are key nodes in this circuitry. Substantial experimental evidence, including neurophysiological recordings, have shown that neurons in these structures represent key aspects of learning. The computational mechanisms that shape the neurophysiological responses, however, are not clear. To examine this, we developed a recurrent neural network (RNN) model of the dlPFC-dSTR circuit and trained it on an oculomotor sequence learning task. We compared the activity generated by the model to activity recorded from monkey dlPFC and dSTR in the same task. This network consisted of a striatal component which encoded action values, and a prefrontal component which selected appropriate actions. After training, this system was able to autonomously represent and update action values and select actions, thus being able to closely approximate the representational structure in corticostriatal recordings. We found that learning to select the correct actions drove action-sequence representations further apart in activity space, both in the model and in the neural data. The model revealed that learning proceeds by increasing the distance between sequence-specific representations. This makes it more likely that the model will select the appropriate action sequence as learning develops. Our model thus supports the hypothesis that learning in networks drives the neural representations of actions further apart, increasing the probability that the network generates correct actions as learning proceeds. Altogether, this study advances our understanding of how neural circuit dynamics are involved in neural computation, revealing how dynamics in the corticostriatal system support task learning.

Neural codes - Necessary but not sufficient for understanding brain function.

Brains are information processing systems whose operational principles ultimately cannot be understood without resource to information theory. We suggest that understanding how external signals are represented in the brain is a necessary step towards employing further engineering tools (such as control theory) to understand the information processing performed by brain circuits during behaviour.

Signature Patterns for Top-Down and Bottom-Up Information Processing via Cross-Frequency Coupling in Macaque Auditory Cortex.

Predictive coding is a theoretical framework that provides a functional interpretation of top-down and bottom-up interactions in sensory processing. The theory suggests there are differences in message passing up versus down the cortical hierarchy. These differences result from the linear feedforward of prediction errors, and the nonlinear feedback of predictions. This implies that cross-frequency interactions should predominate top-down. But it remains unknown whether these differences are expressed in cross-frequency interactions in the brain. Here we examined bidirectional cross-frequency coupling across four sectors of the auditory hierarchy in the macaque. We computed two measures of cross-frequency coupling, phase-amplitude coupling (PAC) and amplitude-amplitude coupling (AAC). Our findings revealed distinct patterns for bottom-up and top-down information processing among cross-frequency interactions. Both top-down and bottom-up interactions made prominent use of low frequencies: low-to-low-frequency (theta, alpha, beta) and low-frequency-to-high- gamma couplings were predominant top-down, while low-frequency-to-low-gamma couplings were predominant bottom-up. These patterns were largely preserved across coupling types (PAC and AAC) and across stimulus types (natural and synthetic auditory stimuli), suggesting that they are a general feature of information processing in auditory cortex. Our findings suggest the modulatory effect of low frequencies on gamma-rhythms in distant regions is important for bidirectional information transfer. The finding of low-frequency-to-low-gamma interactions in the bottom-up direction suggest that nonlinearities may also play a role in feedforward message passing. Altogether, the patterns of cross-frequency interaction we observed across the auditory hierarchy are largely consistent with the predictive coding framework.

All-optical crosstalk-free manipulation and readout of Chronos-expressing neurons.

All optical neurophysiology allows manipulation and readout of neural network activity with single-cell spatial resolution and millisecond temporal resolution. Neurons can be made to express proteins that actuate transmembrane currents upon light absorption, enabling optical control of membrane potential and action potential signalling. In addition, neurons can be genetically or synthetically labelled with fluorescent reporters of changes in intracellular calcium concentration or membrane potential. Thus, to optically manipulate and readout neural activity in parallel, two spectra are involved: the action spectrum of the actuator, and the absorption spectrum of the fluorescent reporter. Due to overlap in these spectra, previous all-optical neurophysiology paradigms have been hindered by spurious activation of neuronal activity caused by the readout light. Here, we pair the blue-green absorbing optogenetic actuator, Chronos, with a deep red-emitting fluorescent calcium reporter CaSiR-1. We show that cultured Chinese hamster ovary cells transfected with Chronos do not exhibit transmembrane currents when illuminated with wavelengths and intensities suitable for exciting one-photon CaSiR-1 fluorescence. We then demonstrate crosstalk-free, high signal-to-noise ratio CaSiR-1 red fluorescence imaging at 100 frames s-1 of Chronos-mediated calcium transients evoked in neurons with blue light pulses at rates up to 20 Hz. These results indicate that the spectral separation between red light excited fluorophores, excited efficiently at or above 640 nm, with blue-green absorbing opsins such as Chronos, is sufficient to avoid spurious opsin actuation by the imaging wavelengths and therefore enable crosstalk-free all-optical neuronal manipulation and readout.

Corrigendum: Single-Neuron Level One-Photon Voltage Imaging With Sparsely Targeted Genetically Encoded Voltage Indicators.

[This corrects the article DOI: 10.3389/fncel.2019.00039.].

Correction to: PyPNS: Multiscale Simulation of a Peripheral Nerve in Python.

The original version of this article unfortunately contained a mistake. The following text: "This project has received funding from European Research Council (ERC) Synergy Grant no. 319818." is missing in the Acknowledgments.

Single-Neuron Level One-Photon Voltage Imaging With Sparsely Targeted Genetically Encoded Voltage Indicators.

Voltage imaging of many neurons simultaneously at single-cell resolution is hampered by the difficulty of detecting small voltage signals from overlapping neuronal processes in neural tissue. Recent advances in genetically encoded voltage indicator (GEVI) imaging have shown single-cell resolution optical voltage recordings in intact tissue through imaging naturally sparse cell classes, sparse viral expression, soma restricted expression, advanced optical systems, or a combination of these. Widespread sparse and strong transgenic GEVI expression would enable straightforward optical access to a densely occurring cell type, such as cortical pyramidal cells. Here we demonstrate that a recently described sparse transgenic expression strategy can enable single-cell resolution voltage imaging of cortical pyramidal cells in intact brain tissue without restricting expression to the soma. We also quantify the functional crosstalk in brain tissue and discuss optimal imaging rates to inform future GEVI experimental design.

Activatable cell-biomaterial interfacing with photo-caged peptides.

Spatio-temporally tailoring cell-material interactions is essential for developing smart delivery systems and intelligent biointerfaces. Here we report new photo-activatable cell-material interfacing systems that trigger cellular uptake of various cargoes and cell adhesion towards surfaces. To achieve this, we designed a novel photo-caged peptide which undergoes a structural transition from an antifouling ligand to a cell-penetrating peptide upon photo-irradiation. When the peptide is conjugated to ligands of interest, we demonstrate the photo-activated cellular uptake of a wide range of cargoes, including small fluorophores, proteins, inorganic (e.g., quantum dots and gold nanostars) and organic nanomaterials (e.g., polymeric particles), and liposomes. Using this system, we can remotely regulate drug administration into cancer cells by functionalizing camptothecin-loaded polymeric nanoparticles with our synthetic peptide ligands. Furthermore, we show light-controlled cell adhesion on a peptide-modified surface and 3D spatiotemporal control over cellular uptake of nanoparticles using two-photon excitation. We anticipate that the innovative approach proposed in this work will help to establish new stimuli-responsive delivery systems and biomaterials.

PyPNS: Multiscale Simulation of a Peripheral Nerve in Python.

Bioelectronic Medicines that modulate the activity patterns on peripheral nerves have promise as a new way of treating diverse medical conditions from epilepsy to rheumatism. Progress in the field builds upon time consuming and expensive experiments in living organisms. To reduce experimentation load and allow for a faster, more detailed analysis of peripheral nerve stimulation and recording, computational models incorporating experimental insights will be of great help. We present a peripheral nerve simulator that combines biophysical axon models and numerically solved and idealised extracellular space models in one environment. We modelled the extracellular space as a three-dimensional resistive continuum governed by the electro-quasistatic approximation of the Maxwell equations. Potential distributions were precomputed in finite element models for different media (homogeneous, nerve in saline, nerve in cuff) and imported into our simulator. Axons, on the other hand, were modelled more abstractly as one-dimensional chains of compartments. Unmyelinated fibres were based on the Hodgkin-Huxley model; for myelinated fibres, we adapted the model proposed by McIntyre et al. in 2002 to smaller diameters. To obtain realistic axon shapes, an iterative algorithm positioned fibres along the nerve with a variable tortuosity fit to imaged trajectories. We validated our model with data from the stimulated rat vagus nerve. Simulation results predicted that tortuosity alters recorded signal shapes and increases stimulation thresholds. The model we developed can easily be adapted to different nerves, and may be of use for Bioelectronic Medicine research in the future.

High speed functional imaging with source localized multifocal two-photon microscopy.

Multifocal two-photon microscopy (MTPM) increases imaging speed over single-focus scanning by parallelizing fluorescence excitation. The imaged fluorescence's susceptibility to crosstalk, however, severely degrades contrast in scattering tissue. Here we present a source-localized MTPM scheme optimized for high speed functional fluorescence imaging in scattering mammalian brain tissue. A rastered line array of beamlets excites fluorescence imaged with a complementary metal-oxide-semiconductor (CMOS) camera. We mitigate scattering-induced crosstalk by temporally oversampling the rastered image, generating grouped images with structured illumination, and applying Richardson-Lucy deconvolution to reassign scattered photons. Single images are then retrieved with a maximum intensity projection through the deconvolved image groups. This method increased image contrast at depths up to 112 µm in scattering brain tissue and reduced functional crosstalk between pixels during neuronal calcium imaging. Source-localization did not affect signal-to-noise ratio (SNR) in densely labeled tissue under our experimental conditions. SNR decreased at low frame rates in sparsely labeled tissue, with no effect at frame rates above 50 Hz. Our non-descanned source-localized MTPM system enables high SNR, 100 Hz capture of fluorescence transients in scattering brain, increasing the scope of MTPM to faster and smaller functional signals.