Spatiotemporal dynamics of sensory responses in layer 2/3 of rat barrel cortex measured in vivo by voltage-sensitive dye imaging combined with whole-cell voltage recordings and neuron reconstructions

Induction of excitatory brain state governs plastic functional changes in visual cortical topology

Adult visual plasticity underlying local remodeling of the cortical circuitry in vivo appears to be associated with a spatiotemporal pattern of strongly increased spontaneous and evoked activity of populations of cells. Here we review and discuss pioneering work by us and others about principles of plasticity in the adult visual cortex, starting with our study which showed that a confined lesion in the cat retina causes increased excitability in the affected region in the primary visual cortex accompanied by fine-tuned restructuring of neuronal function. The underlying remodeling processes was further visualized with voltage-sensitive dye (VSD) imaging that allowed a direct tracking of retinal lesion-induced reorganization across horizontal cortical circuitries. Nowadays, application of noninvasive stimulation methods pursues the idea further of increased cortical excitability along with decreased inhibition as key factors for the induction of adult cortical plasticity. We used high-frequency transcranial magnetic stimulation (TMS), for the first time in combination with VSD optical imaging, and provided evidence that TMS-amplified excitability across large pools of neurons forms the basis for noninvasively targeting reorganization of orientation maps in the visual cortex. Our review has been compiled on the basis of these four own studies, which we discuss in the context of historical developments in the field of visual cortical plasticity and the current state of the literature. Overall, we suggest markers of LTP-like cortical changes at mesoscopic population level as a main driving force for the induction of visual plasticity in the adult. Elevations in excitability that predispose towards cortical plasticity are most likely a common property of all cortical modalities. Thus, interventions that increase cortical excitability are a promising starting point to drive perceptual and potentially motor learning in therapeutic applications.

Encoding luminance surfaces in the visual cortex of mice and monkeys: difference in responses to edge and center

Luminance and spatial contrast provide information on the surfaces and edges of objects. We investigated neural responses to black and white surfaces in the primary visual cortex (V1) of mice and monkeys. Unlike primates that use their fovea to inspect objects with high acuity, mice lack a fovea and have low visual acuity. It thus remains unclear whether monkeys and mice share similar neural mechanisms to process surfaces. The animals were presented with white or black surfaces and the population responses were measured at high spatial and temporal resolution using voltage-sensitive dye imaging. In mice, the population response to the surface was not edge-dominated with a tendency to center-dominance, whereas in monkeys the response was edge-dominated with a "hole" in the center of the surface. The population response to the surfaces in both species exhibited suppression relative to a grating stimulus. These results reveal the differences in spatial patterns to luminance surfaces in the V1 of mice and monkeys and provide evidence for a shared suppression process relative to grating.

Propagating population activity patterns during spontaneous slow waves in the thalamus of rodents

Slow waves (SWs) represent the most prominent electrophysiological events in the thalamocortical system under anesthesia and during deep sleep. Recent studies have revealed that SWs have complex spatiotemporal dynamics and propagate across neocortical regions. However, it is still unclear whether neuronal activity in the thalamus exhibits similar propagation properties during SWs. Here, we report propagating population activity in the thalamus of ketamine/xylazine-anesthetized rats and mice visualized by high-density silicon probe recordings. In both rodent species, propagation of spontaneous thalamic activity during up-states was most frequently observed in dorsal thalamic nuclei such as the higher order posterior (Po), lateral posterior (LP) or laterodorsal (LD) nuclei. The preferred direction of thalamic activity spreading was along the dorsoventral axis, with over half of the up-states exhibiting a gradual propagation in the ventral-to-dorsal direction. Furthermore, simultaneous neocortical and thalamic recordings collected under anesthesia demonstrated that there is a weak but noticeable interrelation between propagation patterns observed during cortical up-states and those displayed by thalamic population activity. In addition, using chronically implanted silicon probes, we detected propagating activity patterns in the thalamus of naturally sleeping rats during slow-wave sleep. However, in comparison to propagating up-states observed under anesthesia, these propagating patterns were characterized by a reduced rate of occurrence and a faster propagation speed. Our findings suggest that the propagation of spontaneous population activity is an intrinsic property of the thalamocortical network during synchronized brain states such as deep sleep or anesthesia. Additionally, our data implies that the neocortex may have partial control over the formation of propagation patterns within the dorsal thalamus under anesthesia.

Dopamine signaling from ganglion cells directs layer-specific angiogenesis in the retina

During central nervous system (CNS) development, a precisely patterned vasculature emerges to support CNS function. How neurons control angiogenesis is not well understood. Here, we show that the neuromodulator dopamine restricts vascular development in the retina via temporally limited production by an unexpected neuron subset. Our genetic and pharmacological experiments demonstrate that elevating dopamine levels inhibits tip-cell sprouting and vessel growth, whereas reducing dopamine production by all retina neurons increases growth. Dopamine production by canonical dopaminergic amacrine interneurons is dispensable for these events. Instead, we found that temporally restricted dopamine production by retinal ganglion cells (RGCs) modulates vascular development. RGCs produce dopamine precisely during angiogenic periods. Genetically limiting dopamine production by ganglion cells, but not amacrines, decreases angiogenesis. Conversely, elevating ganglion-cell-derived dopamine production inhibits early vessel growth. These vasculature outcomes occur downstream of vascular endothelial growth factor receptor (VEGFR) activation and Notch-Jagged1 signaling. Jagged1 is increased and subsequently inhibits Notch signaling when ganglion cell dopamine production is reduced. Our findings demonstrate that dopaminergic neural activity from a small neuron subset functions upstream of VEGFR to serve as developmental timing cue that regulates vessel growth.

Velocity of conduction between columns and layers in barrel cortex reported by parvalbumin interneurons

Inhibitory interneurons expressing parvalbumin (PV) play critical roles throughout the brain. Their rapid spiking enables them to control circuit dynamics on a millisecond time scale, and the timing of their activation by different excitatory pathways is critical to these functions. We used a genetically encoded hybrid voltage sensor to image PV interneuron voltage changes with sub-millisecond precision in primary somatosensory barrel cortex (BC) of adult mice. Electrical stimulation evoked depolarizations with a latency that increased with distance from the stimulating electrode, allowing us to determine conduction velocity. Spread of responses between cortical layers yielded an interlaminar conduction velocity and spread within layers yielded intralaminar conduction velocities in different layers. Velocities ranged from 74 to 473 μm/ms depending on trajectory; interlaminar conduction was 71% faster than intralaminar conduction. Thus, computations within columns are more rapid than between columns. The BC integrates thalamic and intracortical input for functions such as texture discrimination and sensory tuning. Timing differences between intra- and interlaminar PV interneuron activation could impact these functions. Imaging of voltage in PV interneurons reveals differences in signaling dynamics within cortical circuitry. This approach offers a unique opportunity to investigate conduction in populations of axons based on their targeting specificity.

Functional States Shape the Spatiotemporal Representation of Local and Cortex-wide Neural Activity in Mouse Sensory Cortex

The spatiotemporal representation of neural activity during rest and upon sensory stimulation in cortical areas is highly dynamic, and may be predominantly governed by cortical state. On the mesoscale level, intrinsic neuronal activity ranges from a persistent state, generally associated with a sustained depolarization of neurons, to a bimodal, slow-wave like state with bursts of neuronal activation, alternating with silent periods. These different activity states are prevalent under certain types of sedatives, or are associated with specific behavioral or vigilance conditions. Neurophysiological experiments assessing circuit activity, usually assume a constant underlying state, yet reports of variability of neuronal responses under seemingly constant conditions are common in the field. Even when a certain type of neural activity or cortical state can stably be maintained over time, the associated response properties are highly relevant for explaining experimental outcomes. Here we describe the spatiotemporal characteristics of ongoing activity and sensory evoked responses under two predominant functional states in the sensory cortices of mice: persistent activity (PA) and slow wave activity (SWA). Using electrophysiological recordings, and local and wide-field calcium recordings, we examine whether spontaneous and sensory evoked neuronal activity propagate throughout the cortex in a state dependent manner. We find that PA and SWA differ in their spatiotemporal characteristics which determine the cortical network's response to a sensory stimulus. During PA state, sensory stimulation elicits gamma-based short-latency responses which precisely follow each stimulation pulse and are prone to adaptation upon higher stimulation frequencies. Sensory responses during SWA are more variable, dependent on refractory periods following spontaneous slow waves. While spontaneous slow waves propagated in anterior-posterior direction in a majority of observations, the direction of propagation of stimulus-elicited wave depends on the sensory modality. These findings suggest that cortical state explains variance and should be considered when investigating multi-scale correlates of functional neurocircuit activity.

Spatially distributed computation in cortical circuits

The traditional view of neural computation in the cerebral cortex holds that sensory neurons are specialized, i.e., selective for certain dimensions of sensory stimuli. This view was challenged by evidence of contextual interactions between stimulus dimensions in which a neuron's response to one dimension strongly depends on other dimensions. Here, we use methods of mathematical modeling, psychophysics, and electrophysiology to address shortcomings of the traditional view. Using a model of a generic cortical circuit, we begin with the simple demonstration that cortical responses are always distributed among neurons, forming characteristic waveforms, which we call neural waves. When stimulated by patterned stimuli, circuit responses arise by interference of neural waves. Results of this process depend on interaction between stimulus dimensions. Comparison of modeled responses with responses of biological vision makes it clear that the framework of neural wave interference provides a useful alternative to the standard concept of neural computation.

Altered cortical processing of sensory input in Huntington disease mouse models

Huntington disease (HD), a hereditary neurodegenerative disorder, manifests as progressively impaired movement and cognition. Although early abnormalities of neuronal activity in striatum are well established in HD models, there are fewer in vivo studies of the cortex. Here, we record local field potentials (LFPs) in YAC128 HD model mice versus wild-type mice. In multiple cortical areas, limb sensory stimulation evokes a greater change in LFP power in YAC128 mice. Mesoscopic imaging using voltage-sensitive dyes reveals more extensive spread of evoked sensory signals across the cortical surface in YAC128 mice. YAC128 layer 2/3 sensory cortical neurons ex vivo show increased excitatory events, which could contribute to enhanced sensory responses in vivo. Cortical LFP responses to limb stimulation, visual and auditory input are also significantly increased in zQ175 HD mice. Results presented here extend knowledge of HD beyond ex vivo studies of individual neurons to the intact cortical network.

Simultaneous Two-Photon Voltage or Calcium Imaging and Multi-Channel Local Field Potential Recordings in Barrel Cortex of Awake and Anesthetized Mice

Neuronal population activity, both spontaneous and sensory-evoked, generates propagating waves in cortex. However, high spatiotemporal-resolution mapping of these waves is difficult as calcium imaging, the work horse of current imaging, does not reveal subthreshold activity. Here, we present a platform combining voltage or calcium two-photon imaging with multi-channel local field potential (LFP) recordings in different layers of the barrel cortex from anesthetized and awake head-restrained mice. A chronic cranial window with access port allows injecting a viral vector expressing GCaMP6f or the voltage-sensitive dye (VSD) ANNINE-6plus, as well as entering the brain with a multi-channel neural probe. We present both average spontaneous activity and average evoked signals in response to multi-whisker air-puff stimulations. Time domain analysis shows the dependence of the evoked responses on the cortical layer and on the state of the animal, here separated into anesthetized, awake but resting, and running. The simultaneous data acquisition allows to compare the average membrane depolarization measured with ANNINE-6plus with the amplitude and shape of the LFP recordings. The calcium imaging data connects these data sets to the large existing database of this important second messenger. Interestingly, in the calcium imaging data, we found a few cells which showed a decrease in calcium concentration in response to vibrissa stimulation in awake mice. This system offers a multimodal technique to study the spatiotemporal dynamics of neuronal signals through a 3D architecture in vivo. It will provide novel insights on sensory coding, closing the gap between electrical and optical recordings.

Optimization of Near-Infrared Fluorescence Voltage-Sensitive Dye Imaging for Neuronal Activity Monitoring in the Rodent Brain

Many currently employed clinical brain functional imaging technologies rely on indirect measures of activity such as hemodynamics resulting in low temporal and spatial resolutions. To improve upon this, optical systems were developed in conjunction with methods to deliver near-IR voltage-sensitive dye (VSD) to provide activity-dependent optical contrast to establish a clinical tool to facilitate direct monitoring of neuron depolarization through the intact skull. Following the previously developed VSD delivery protocol through the blood-brain barrier, IR-780 perchlorate VSD concentrations in the brain were varied and stimulus-evoked responses were observed. In this paper, a range of optimal VSD tissue concentrations was established that maximized fluorescence fractional change for detection of membrane potential responses to external stimuli through a series of phantom, in vitro, ex vivo, and in vivo experiments in mouse models.

Studying Synaptically Evoked Cortical Responses ex vivo With Combination of a Single Neuron Recording (Whole-Cell) and Population Voltage Imaging (Genetically Encoded Voltage Indicator)

In a typical electrophysiology experiment, synaptic stimulus is delivered in a cortical layer (1-6) and neuronal responses are recorded intracellularly in individual neurons. We recreated this standard electrophysiological paradigm in brain slices of mice expressing genetically encoded voltage indicators (GEVIs). This allowed us to monitor membrane voltages in the target pyramidal neurons (whole-cell), and population voltages in the surrounding neuropil (optical imaging), simultaneously. Pyramidal neurons have complex dendritic trees that span multiple cortical layers. GEVI imaging revealed areas of the brain slice that experienced the strongest depolarization on a specific synaptic stimulus (location and intensity), thus identifying cortical layers that contribute the most afferent activity to the recorded somatic voltage waveform. By combining whole-cell with GEVI imaging, we obtained a crude distribution of activated synaptic afferents in respect to the dendritic tree of a pyramidal cell. Synaptically evoked voltage waves propagating through the cortical neuropil (dendrites and axons) were not static but rather they changed on a millisecond scale. Voltage imaging can identify areas of brain slices in which the neuropil was in a sustained depolarization (plateau), long after the stimulus onset. Upon a barrage of synaptic inputs, a cortical pyramidal neuron experiences: (a) weak temporal summation of evoked voltage transients (EPSPs); and (b) afterhyperpolarization (intracellular recording), which are not represented in the GEVI population imaging signal (optical signal). To explain these findings [(a) and (b)], we used four voltage indicators (ArcLightD, chi-VSFP, Archon1, and di-4-ANEPPS) with different optical sensitivity, optical response speed, labeling strategy, and a target neuron type. All four imaging methods were used in an identical experimental paradigm: layer 1 (L1) synaptic stimulation, to allow direct comparisons. The population voltage signal showed paired-pulse facilitation, caused in part by additional recruitment of new neurons and dendrites. "Synaptic stimulation" delivered in L1 depolarizes almost an entire cortical column to some degree.

Nonlinear collision between propagating waves in mouse somatosensory cortex

How does cellular organization shape the spatio-temporal patterns of activity in the cortex while processing sensory information? After measuring the propagation of activity in the mouse primary somatosensory cortex (S1) in response to single whisker deflections with Voltage Sensitive Dye (VSD) imaging, we developed a two dimensional model of S1. We designed an inference method to reconstruct model parameters from VSD data, revealing that a spatially heterogeneous organization of synaptic strengths between pyramidal neurons in S1 is likely to be responsible for the heterogeneous spatio-temporal patterns of activity measured experimentally. The model shows that, for strong enough excitatory cortical interactions, whisker deflections generate a propagating wave in S1. Finally, we report that two consecutive stimuli activating different spatial locations in S1 generate two waves which collide sub-linearly, giving rise to a suppressive wave. In the inferred model, the suppressive wave is explained by a lower sensitivity to external perturbations of neural networks during activated states.

Global, Low-Amplitude Cortical State Predicts Response Outcomes in a Selective Detection Task in Mice

Spontaneous neuronal activity strongly impacts stimulus encoding and behavioral responses. We sought to determine the effects of neocortical prestimulus activity on stimulus detection. We trained mice in a selective whisker detection task, in which they learned to respond (lick) to target stimuli in one whisker field and ignore distractor stimuli in the contralateral whisker field. During expert task performance, we used widefield Ca2+ imaging to assess prestimulus and post-stimulus neuronal activity broadly across frontal and parietal cortices. We found that lower prestimulus activity correlated with enhanced stimulus detection: lower prestimulus activity predicted response versus no response outcomes and faster reaction times. The activity predictive of trial outcome was distributed through dorsal neocortex, rather than being restricted to whisker or licking regions. Using principal component analysis, we demonstrate that response trials are associated with a distinct and less variable prestimulus neuronal subspace. For single units, prestimulus choice probability was weak yet distributed broadly, with lower than chance choice probability correlating with stronger sensory and motor encoding. These findings support low amplitude and low variability as an optimal prestimulus cortical state for stimulus detection that presents globally and predicts response outcomes for both target and distractor stimuli.

Cholinergic modulation of sensory processing in awake mouse cortex

Cholinergic modulation of brain activity is fundamental for awareness and conscious sensorimotor behaviours, but deciphering the timing and significance of acetylcholine actions for these behaviours is challenging. The widespread nature of cholinergic projections to the cortex means that new insights require access to specific neuronal populations, and on a time-scale that matches behaviourally relevant cholinergic actions. Here, we use fast, voltage imaging of L2/3 cortical pyramidal neurons exclusively expressing the genetically-encoded voltage indicator Butterfly 1.2, in awake, head-fixed mice, receiving sensory stimulation, whilst manipulating the cholinergic system. Altering muscarinic acetylcholine function re-shaped sensory-evoked fast depolarisation and subsequent slow hyperpolarisation of L2/3 pyramidal neurons. A consequence of this re-shaping was disrupted adaptation of the sensory-evoked responses, suggesting a critical role for acetylcholine during sensory discrimination behaviour. Our findings provide new insights into how the cortex processes sensory information and how loss of acetylcholine, for example in Alzheimer's Disease, disrupts sensory behaviours.

An increase in dendritic plateau potentials is associated with experience-dependent cortical map reorganization

Here we describe a mechanism for cortical map plasticity. Classically, representational map changes are thought to be driven by changes within cortico-cortical circuits, e.g., Hebbian plasticity of synaptic circuits that lost vs. maintained an excitatory drive from the first-order thalamus, possibly steered by neuromodulatory forces from deep brain regions. Our work provides evidence for an additional gating mechanism, provided by plateau potentials, which are driven by higher-order thalamic feedback. Higher-order thalamic neurons are characterized by broad receptive fields, and the plateau potentials that they evoke strongly facilitate long-term potentiation and elicit spikes. We show that these features combined constitute a powerful driving force for the fusion or expansion of sensory representations within cortical maps.

The organization of sensory maps in the cerebral cortex depends on experience, which drives homeostatic and long-term synaptic plasticity of cortico-cortical circuits. In the mouse primary somatosensory cortex (S1) afferents from the higher-order, posterior medial thalamic nucleus (POm) gate synaptic plasticity in layer (L) 2/3 pyramidal neurons via disinhibition and the production of dendritic plateau potentials. Here we address whether these thalamocortically mediated responses play a role in whisker map plasticity in S1. We find that trimming all but two whiskers causes a partial fusion of the representations of the two spared whiskers, concomitantly with an increase in the occurrence of POm-driven N-methyl-D-aspartate receptor-dependent plateau potentials. Blocking the plateau potentials restores the archetypical organization of the sensory map. Our results reveal a mechanism for experience-dependent cortical map plasticity in which higher-order thalamocortically mediated plateau potentials facilitate the fusion of normally segregated cortical representations.

In silico voltage-sensitive dye imaging reveals the emergent dynamics of cortical populations

Voltage-sensitive dye imaging (VSDI) is a powerful technique for interrogating membrane potential dynamics in assemblies of cortical neurons, but with effective resolution limits that confound interpretation. To address this limitation, we developed an in silico model of VSDI in a biologically faithful digital reconstruction of rodent neocortical microcircuitry. Using this model, we extend previous experimental observations regarding the cellular origins of VSDI, finding that the signal is driven primarily by neurons in layers 2/3 and 5, and that VSDI measurements do not capture individual spikes. Furthermore, we test the capacity of VSD image sequences to discriminate between afferent thalamic inputs at various spatial locations to estimate a lower bound on the functional resolution of VSDI. Our approach underscores the power of a bottom-up computational approach for relating scales of cortical processing.

In Vivo Whole-Cell Patch-Clamp Methods: Recent Technical Progress and Future Perspectives

Brain functions are fundamental for the survival of organisms, and they are supported by neural circuits consisting of a variety of neurons. To investigate the function of neurons at the single-cell level, researchers often use whole-cell patch-clamp recording techniques. These techniques enable us to record membrane potentials (including action potentials) of individual neurons of not only anesthetized but also actively behaving animals. This whole-cell recording method enables us to reveal how neuronal activities support brain function at the single-cell level. In this review, we introduce previous studies using in vivo patch-clamp recording techniques and recent findings primarily regarding neuronal activities in the hippocampus for behavioral function. We further discuss how we can bridge the gap between electrophysiology and biochemistry.

Learning poly-synaptic paths with traveling waves

Traveling waves are commonly observed across the brain. While previous studies have suggested the role of traveling waves in learning, the mechanism remains unclear. We adopted a computational approach to investigate the effect of traveling waves on synaptic plasticity. Our results indicate that traveling waves facilitate the learning of poly-synaptic network paths when combined with a reward-dependent local synaptic plasticity rule. We also demonstrate that traveling waves expedite finding the shortest paths and learning nonlinear input/output mapping, such as exclusive or (XOR) function.

There are approximately 10¹¹ neurons with 10¹⁴ connections in the human brain. Information transmission among neurons in this large network is considered crucial for our behavior. To achieve this, multiple synaptic connections along a poly-synaptic network path must be adjusted coherently during learning. Because the previously proposed reward-dependent synaptic plasticity rule requires coactivation of presynaptic and postsynaptic neurons, learning can fail if a subset of neurons along a distant network path is inactive at the beginning of learning. We suggest that traveling waves that are initiated at an information source can mitigate this problem. We performed computer simulations of spiking neural networks with reward-dependent local synaptic plasticity rules and traveling waves. Our results show that this combination facilitates the learning and refinement of synaptic network paths. We argue that these features are a general biological strategy for maintaining and optimizing our brain function. Our research provides new insights into how complex neural networks in the brain form during learning and memory consolidation.

Optical Analysis of Acute Changes after Peripheral Nerve Injury in Spatio-Temporal Pattern of Neural Response to Forelimb Stimulation in Rat Somatosensory Cortex

Peripheral nerve injury induces functional reorganization of the central nervous system. The mechanisms underlying this reorganization have been widely studied. Our previous study involving multiple-site optical recording reported that a neural excitatory wave induced by somatic stimulation begins in a small area and propagates in the cortex. In the present study, to examine the possible role of this propagation wave in cortical reorganization, we analyzed the early changes in the spatio-temporal pattern of the sensory-evoked wave immediately, and 30 min, after nerve injury. The response to hypothenar stimulation, innervated by the ulnar nerve and adjoining the median nerve area, persisted after injury to either the ulnar or median nerve. Initially, we assessed changes in the response pattern at the focus. The latency increased after ulnar nerve injury, whereas no change was observed after median nerve injury. Similarly, no change was noted in the duration of the response signal with either nerve injury. Second, changes in the propagation wave pattern were analyzed. Ulnar nerve injury decreased the propagation velocity in the medial direction but the median nerve injury induced no changes. These results indicated that the propagation wave pattern is readily altered, even immediately after nerve injury, and suggest that this immediate change in the spatio-temporal pattern is one of the factors contributing to the cortical reorganization.

Global enhancement of cortical excitability following coactivation of large neuronal populations

Correlated activation of cortical neurons often occurs in the brain and repetitive correlated neuronal firing could cause long-term modifications of synaptic efficacy and intrinsic excitability. We found that repetitive optogenetic activation of neuronal populations in the mouse cortex caused enhancement of optogenetically evoked firing of local coactivated neurons as well as distant cortical neurons in both ipsilateral and contralateral hemispheres. This global enhancement of evoked responses required coactivation of a sufficiently large population of neurons either within one cortical area or distributed in several areas. Enhancement of neuronal firing was saturable after repeated episodes of coactivation, diminished by inhibition of N-methyl-d-aspartic acid receptors, and accompanied by elevated excitatory postsynaptic potentials, all consistent with activity-induced synaptic potentiation. Chemogenetic inhibition of neuronal activity of the thalamus decreased the enhancement effect, suggesting thalamic involvement. Thus, correlated excitation of large neuronal populations leads to global enhancement of neuronal excitability.

Spatial and Spectral Mapping and Decomposition of Neural Dynamics and Organization of the Mouse Brain with Multispectral Optoacoustic Tomography

In traditional optical imaging, limited light penetration constrains high-resolution interrogation to tissue surfaces. Optoacoustic imaging combines the superb contrast of optical imaging with deep penetration of ultrasound, enabling a range of new applications. We used multispectral optoacoustic tomography (MSOT) for functional and structural neuroimaging in mice at resolution, depth, and specificity unattainable by other neuroimaging modalities. Based on multispectral readouts, we computed hemoglobin gradient and oxygen saturation changes related to processing of somatosensory signals in different structures along the entire subcortical-cortical axis. Using temporal correlation analysis and seed-based maps, we reveal the connectivity between cortical, thalamic, and sub-thalamic formations. With the same modality, high-resolution structural tomography of intact mouse brain was achieved based on endogenous contrasts, demonstrating near-perfect matches with anatomical features revealed by histology. These results extend the limits of noninvasive observations beyond the reach of standard high-resolution neuroimaging, verifying the suitability of MSOT for small-animal studies.

Ketamine anaesthesia induces gain enhancement via recurrent excitation in granular input layers of the auditory cortex

KEY POINTS

Ketamine is a common anaesthetic agent used in research and more recently as medication in treatment of depression. It has known effects on inhibition of interneurons and cortical stimulus-locked responses, but the underlying functional network mechanisms are still elusive. Analysing population activity across all layers within the auditory cortex, we found that doses of this anaesthetic induce a stronger activation and stimulus-locked response to pure-tone stimuli. This cortical response is driven by gain enhancement of thalamocortical input processing selectively within granular layers due to an increased recurrent excitation. Time-frequency analysis indicates a higher broadband magnitude response and prolonged phase coherence in granular layers, possibly pointing to disinhibition of this recurrent excitation. These results further the understanding of ketamine's functional mechanisms, which will improve the ability to interpret physiological studies moving from anaesthetized to awake paradigms and may lead to the development of better ketamine-based depression treatments with lower side effects.

ABSTRACT

Ketamine is commonly used as an anaesthetic agent and has more recently gained attention as an antidepressant. It has been linked to increased stimulus-locked excitability, inhibition of interneurons and modulation of intrinsic neuronal oscillations. However, the functional network mechanisms are still elusive. A better understanding of these anaesthetic network effects may improve upon previous interpretations of seminal studies conducted under anaesthesia and have widespread relevance for neuroscience with awake and anaesthetized subjects as well as in medicine. Here, we investigated the effects of anaesthetic doses of ketamine (15 mg kg^-1  h^-1 i.p.) on the network activity after pure-tone stimulation within the auditory cortex of male Mongolian gerbils (Meriones unguiculatus). We used laminar current source density (CSD) analysis and subsequent layer-specific continuous wavelet analysis to investigate spatiotemporal response dynamics on cortical columnar processing in awake and ketamine-anaesthetized animals. We found thalamocortical input processing within granular layers III/IV to be significantly increased under ketamine. This layer-dependent gain enhancement under ketamine was not due to changes in cross-trial phase coherence but was rather attributed to a broadband increase in magnitude reflecting an increase in recurrent excitation. A time-frequency analysis was indicative of a prolonged period of stimulus-induced excitation possibly due to a reduced coupling of excitation and inhibition in granular input circuits - in line with the common hypothesis of cortical disinhibition via suppression of GABAergic interneurons.

Perspective on the Multiple Pathways to Changing Brain States

In this review article, we highlight several disparate ideas that are linked to changes in brain state (i.e., sleep to arousal, Down to Up, synchronized to de-synchronized). In any discussion of the brain state, we propose that the cortical pyramidal neuron has a central position. EEG recordings, which typically assess brain state, predominantly reflect the activity of cortical pyramidal neurons. This means that the dominant rhythmic activity that characterizes a particular brain state ultimately has to manifest globally across the pyramidal neuron population. During state transitions, it is the long-range connectivity of these neurons that broadcast the resultant changes in activity to many subcortical targets. Structures like the thalamus, brainstem/hypothalamic neuromodulatory systems, and respiratory systems can also strongly influence brain state, and for many decades we have been uncovering bidirectional pathways that link these structures to state changes in the cerebral cortex. More recently, movement and active behaviors have emerged as powerful drivers of state changes. Each of these systems involve different circuits distributed across the brain. Yet, for a system-wide change in brain state, there must be a collaboration between these circuits that reflects and perhaps triggers the transition between brain states. As we expand our understanding of how brain state changes, our current challenge is to understand how these diverse sets of circuits and pathways interact to produce the changes observed in cortical pyramidal neurons.

Propagating Motor Cortical Dynamics Facilitate Movement Initiation

Voluntary movement initiation involves the modulations of large groups of neurons in the primary motor cortex (M1). Yet similar modulations occur during movement planning when no movement occurs. Here, we show that a sequential spatiotemporal pattern of excitability propagates across M1 prior to the movement initiation in one of two oppositely oriented directions along the rostro-caudal axis. Using spatiotemporal patterns of intracortical microstimulation, we find that reaction time increases significantly when stimulation is delivered against, but not with, the natural propagation direction. Functional connections among M1 units emerge at movement that are oriented along the same rostro-caudal axis but not during movement planning. Finally, we show that beta amplitude profiles can more accurately decode muscle activity when they conform to the natural propagating patterns. These findings provide the first causal evidence that large-scale, propagating patterns of cortical excitability are behaviorally relevant and may be a necessary component of movement initiation.

Comparative Evaluation of Genetically Encoded Voltage Indicators

Imaging voltage using fluorescent-based sensors could be an ideal technique to probe neural circuits with high spatiotemporal resolution. However, due to insufficient signal-to-noise ratio (SNR), imaging membrane potential in mammalian preparations is still challenging. In recent years, many genetically encoded voltage indicators (GEVIs) have been developed. To compare them and guide decisions on which GEVI to use, we have characterized side by side the performance of eight GEVIs that represent different families of molecular constructs. We tested GEVIs in vitro with 1-photon imaging and in vivo with 1-photon wide-field imaging and 2-photon imaging. We find that QuasAr2 exhibited the best performance in vitro, whereas only ArcLight-MT could be used to reliably detect electrical activity in vivo with 2-photon excitation. No single GEVI was ideal for every experiment. These results provide a guide for choosing optimal GEVIs for specific applications.

Transcranial photoacoustic imaging of NMDA-evoked focal circuit dynamics in the rat hippocampus

OBJECTIVE

We report the transcranial functional photoacoustic (fPA) neuroimaging of N-methyl-D-aspartate (NMDA) evoked neural activity in the rat hippocampus. Concurrent quantitative electroencephalography (qEEG) and microdialysis were used to record real-time circuit dynamics and excitatory neurotransmitter concentrations, respectively.

APPROACH

We hypothesized that location-specific fPA voltage-sensitive dye (VSD) contrast would identify neural activity changes in the hippocampus which correlate with NMDA-evoked excitatory neurotransmission.

MAIN RESULTS

Transcranial fPA VSD imaging at the contralateral side of the microdialysis probe provided NMDA-evoked VSD responses with positive correlation to extracellular glutamate concentration changes. qEEG validated a wide range of glutamatergic excitation, which culminated in focal seizure activity after a high NMDA dose. We conclude that transcranial fPA VSD imaging can distinguish focal glutamate loads in the rat hippocampus, based on the VSD redistribution mechanism which is sensitive to the electrophysiologic membrane potential.

SIGNIFICANCE

Our results suggest the future utility of this emerging technology in both laboratory and clinical sciences as an innovative functional neuroimaging modality.

Spatiotemporal properties of whisker-evoked tactile responses in the mouse secondary somatosensory cortex

The representation of rodents' mystacial vibrissae within the primary somatosensory (S1) cortex has become a major model for studying the cortical processing of tactile sensory information. However, upon vibrissal stimulation, tactile information first reaches S1 but also, almost simultaneously, the secondary somatosensory cortex (S2). To further understand the role of S2 in the processing of whisker inputs, it is essential to characterize the spatio-temporal properties of whisker-evoked response dynamics in this area. Here we describe the topography of the whiskers representation in the mouse S2 with voltage sensitive dye imaging. Analysis of the spatial properties of the early S2 responses induced by stimulating individually 22 to 24 whiskers revealed that they are spatially ordered in a mirror symmetric map with respect to S1 responses. Evoked signals in S2 and S1 are of similar amplitude and closely correlated at the single trial level. They confirm a short delay (~3 ms) between S1 and S2 early activation. In both S1 and S2 caudo-dorsal whiskers induce stronger responses than rostro-ventral ones. Finally, analysis of early C2-evoked responses indicates a faster activation of neighboring whisker representations in S2 relative to S1, probably due to the reduced size of the whisker map in S2.

Enrichment drives emergence of functional columns and improves sensory coding in the whisker map in L2/3 of mouse S1

Sensory maps in layer (L) 2/3 of rodent cortex lack precise functional column boundaries, and instead exhibit locally heterogeneous (salt-and-pepper) tuning superimposed on smooth global topography. Could this organization be a byproduct of impoverished experience in laboratory housing? We compared whisker map somatotopy in L2/3 and L4 excitatory cells of somatosensory (S1) cortex in normally housed vs. tactile-enriched mice, using GCaMP6s imaging. Normally housed mice had a dispersed, salt-and-pepper whisker map in L2/3, but L4 was more topographically precise. Enrichment (P21 to P46-71) sharpened whisker tuning and decreased, but did not abolish, local tuning heterogeneity. In L2/3, enrichment strengthened and sharpened whisker point representations, and created functional boundaries of tuning similarity and noise correlations at column edges. Thus, enrichment drives emergence of functional columnar topography in S1, and reduces local tuning heterogeneity. These changes predict better touch detection by neural populations within each column.

Sensorimotor processing in the rodent barrel cortex

Tactile sensory information from facial whiskers provides nocturnal tunnel-dwelling rodents, including mice and rats, with important spatial and textural information about their immediate surroundings. Whiskers are moved back and forth to scan the environment (whisking), and touch signals from each whisker evoke sparse patterns of neuronal activity in whisker-related primary somatosensory cortex (wS1; barrel cortex). Whisking is accompanied by desynchronized brain states and cell-type-specific changes in spontaneous and evoked neuronal activity. Tactile information, including object texture and location, appears to be computed in wS1 through integration of motor and sensory signals. wS1 also directly controls whisker movements and contributes to learned, whisker-dependent, goal-directed behaviours. The cell-type-specific neuronal circuitry in wS1 that contributes to whisker sensory perception is beginning to be defined.

Validating Non-invasive EEG Source Imaging Using Optimal Electrode Configurations on a Representative Rat Head Model

The curtain of technical limitations impeding rat multichannel non-invasive electroencephalography (EEG) has risen. Given the importance of this preclinical model, development and validation of EEG source imaging (ESI) is essential. We investigate the validity of well-known human ESI methodologies in rats which individual tissue geometries have been approximated by those extracted from an MRI template, leading also to imprecision in electrode localizations. With the half and fifth sensitivity volumes we determine both the theoretical minimum electrode separation for non-redundant scalp EEG measurements and the electrode sensitivity resolution, which vary over the scalp because of the head geometry. According to our results, electrodes should be at least ~3 to 3.5 mm apart for an optimal configuration. The sensitivity resolution is generally worse for electrodes at the boundaries of the scalp measured region, though, by analogy with human montages, concentrates the sensitivity enough to localize sources. Cramér-Rao lower bounds of source localization errors indicate it is theoretically possible to achieve ESI accuracy at the level of anatomical structures, such as the stimulus-specific somatosensory areas, using the template. More validation for this approximation is provided through the comparison between the template and the individual lead field matrices, for several rats. Finally, using well-accepted inverse methods, we demonstrate that somatosensory ESI is not only expected but also allows exploring unknown phenomena related to global sensory integration. Inheriting the advantages and pitfalls of human ESI, rat ESI will boost the understanding of brain pathophysiological mechanisms and the evaluation of ESI methodologies, new pharmacological treatments and ESI-based biomarkers.

Repetitive transcranial magnetic stimulation recovers cortical map plasticity induced by sensory deprivation due to deafferentiation

KEY POINTS

Partial sensory deprivation (deafferentation) by removing whiskers from the rat snout resulted in a reduced responsiveness of related cortical representations. Repetitive transcranial magnetic stimulation (three blocks of intermittent theta-burst) applied for 5 days in combination with sensory exploration restored the normal responsiveness level of the deafferented barrel cortex. However, intracortical inhibition (lateral and recurrent) appeared to be reduced after repetitive transcranial magnetic stimulation, probably as the cause of improved responsiveness. Repetitive transcranial magnetic stimulation also reduced the asymmetry of the lateral spread of sensory activity.

ABSTRACT

Repetitive transcranial magnetic stimulation (rTMS) modulates human cortical excitability. It has the potential to support recovery to normal cortical function when the excitation-inhibition balance is altered (e.g. after a stroke or loss of sensory input). We tested cortical map plasticity on the basis of sensory responses (local field potentials, LFPs) and expression of neuronal activity marker proteins within the barrel cortex of rats receiving either active or sham rTMS after selective unilateral deafferentation by whiskers plucking. Rats received daily rTMS [intermittent theta-burst (iTBS), active or sham] for 5 days before exploring an enriched environment. Our previous studies indicated a disinhibitory effect of iTBS on cortical activity. Therefore, we also expected disinhibitory effects if deafferentation causes depression of sensory responses. Deafferentation resulted in an acute general reduction of sensory responsiveness and enhanced expression of inhibitory activity markers (GAD67, parvalbumin) in the deafferented hemisphere. Active but not sham-iTBS-rTMS normalized these measures. The stronger caudal-to-frontal horizontal spread of activity across barrels was reduced after deafferentation but not restored after active iTBS, despite generally increased responses. Fitting the LFP data with a computational model of different strengths and types of excitatory and inhibitory connections further revealed an iTBS-induced reduction of lateral and recurrent inhibition as the most probable scenario. Whether the disinhibitory effect of iTBS for the restoration of normal cortical function in the acute phase of depression after deafferentiation is also beneficial in humans remains to be demonstrated. As recently discussed, disinhibition appears to be required to open a window for neuronal plasticity.

Technologies to Study Action Potential Propagation With a Focus on HD-MEAs

Axons convey information in neuronal circuits via reliable conduction of action potentials (APs) from the axon initial segment (AIS) to the presynaptic terminals. Recent experimental findings increasingly evidence that the axonal function is not limited to the simple transmission of APs. Advances in subcellular-resolution recording techniques have shown that axons display activity-dependent modulation in spike shape and conduction velocity, which influence synaptic strength and latency. We briefly review here, how recent methodological developments facilitate the understanding of the axon physiology. We included the three most common methods, i.e., genetically encoded voltage imaging (GEVI), subcellular patch-clamp and high-density microelectrode arrays (HD-MEAs). We then describe the potential of using HD-MEAs in studying axonal physiology in more detail. Due to their robustness, amenability to high-throughput and high spatiotemporal resolution, HD-MEAs can provide a direct functional electrical readout of single cells and cellular ensembles at subcellular resolution. HD-MEAs can, therefore, be employed in investigating axonal pathologies, the effects of large-scale genomic interventions (e.g., with RNAi or CRISPR) or in compound screenings. A combination of extracellular microelectrode arrays (MEAs), intracellular microelectrodes and optical imaging may potentially reveal yet unexplored repertoires of axonal functions.

Focal activation of neuronal circuits induced by microstimulation in the visual cortex

OBJECTIVE

Microstimulation to the cortical tissue applied with penetrating electrodes delivers current that spreads concentrically around the electrode tip and is known to evoke focal visual sensations, i.e. phosphenes. However, to date, there is no direct evidence depicting the spatiotemporal properties of neuronal activity induced immediately after microstimulation and how such activity drives the subsequent local cortical circuits.

APPROACH

In the present study, we imaged the spatiotemporal distribution of action potentials (APs) directly induced by microstimulation and the subsequent trans-synaptic signal propagation using a voltage-sensitive dye (VSD) and a calcium-sensitive dye (CaSD) in slice preparations of the mouse primary visual cortex.

MAIN RESULTS

The directly induced APs were confined to the close vicinity of the electrode tip, and the effective distance of excitation was proportional to the square root of the current intensity. The excitation around the electrode tip in layer IV mainly propagated to layer II/III to further induce the subsequent focal activation in downstream local cortical circuits. The extent of activation in the downstream circuits was restrained by competitive interactions between excitatory and inhibitory signals. Namely, the spread of the excitation to lateral neighbor neurons along the layer II/III was confined by the delayed inhibition that also spread laterally at a faster rate.

SIGNIFICANCE

These observations indicate that dynamic interactions between excitatory and inhibitory signals play a critical role in the focal activation of a cortical circuit in response to intracortical microstimulation and, therefore, in evoking a localized phosphene.

The Three-Dimensional Signal Collection Field for Fiber Photometry in Brain Tissue

Fiber photometry is used to monitor signals from fluorescent indicators in genetically-defined neural populations in behaving animals. Recently, fiber photometry has rapidly expanded and it now provides researchers with increasingly powerful means to record neural dynamics and neuromodulatory action. However, it is not clear how to select the optimal fiber optic given the constraints and goals of a particular experiment. Here, using combined confocal/2-photon microscope, we quantitatively characterize the fluorescence collection properties of various optical fibers in brain tissue. We show that the fiber size plays a major role in defining the volume of the optically sampled brain region, whereas numerical aperture impacts the total amount of collected signal and, marginally, the shape and size of the collection volume. We show that ~80% of the effective signal arises from 105 to 106 μm3 volume extending ~200 μm from the fiber facet for 200 μm core optical fibers. Together with analytical and ray tracing collection maps, our results reveal the light collection properties of different optical fibers in brain tissue, allowing for an accurate selection of the fibers for photometry and helping for a more precise interpretation of measurements in terms of sampled volume.

Implantable Neural Probes for Brain-Machine Interfaces - Current Developments and Future Prospects

A Brain-Machine interface (BMI) allows for direct communication between the brain and machines. Neural probes for recording neural signals are among the essential components of a BMI system. In this report, we review research regarding implantable neural probes and their applications to BMIs. We first discuss conventional neural probes such as the tetrode, Utah array, Michigan probe, and electroencephalography (ECoG), following which we cover advancements in next-generation neural probes. These next-generation probes are associated with improvements in electrical properties, mechanical durability, biocompatibility, and offer a high degree of freedom in practical settings. Specifically, we focus on three key topics: (1) novel implantable neural probes that decrease the level of invasiveness without sacrificing performance, (2) multi-modal neural probes that measure both electrical and optical signals, (3) and neural probes developed using advanced materials. Because safety and precision are critical for practical applications of BMI systems, future studies should aim to enhance these properties when developing next-generation neural probes.

Multimodal Modeling of Neural Network Activity: Computing LFP, ECoG, EEG, and MEG Signals With LFPy 2.0

Recordings of extracellular electrical, and later also magnetic, brain signals have been the dominant technique for measuring brain activity for decades. The interpretation of such signals is however nontrivial, as the measured signals result from both local and distant neuronal activity. In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons. Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed. This allows for the development of computational tools implementing forward models grounded in the biophysics underlying electrical and magnetic measurement modalities. LFPy (LFPy.readthedocs.io) incorporated a well-established scheme for predicting extracellular potentials of individual neurons with arbitrary levels of biological detail. It relies on NEURON (neuron.yale.edu) to compute transmembrane currents of multicompartment neurons which is then used in combination with an electrostatic forward model. Its functionality is now extended to allow for modeling of networks of multicompartment neurons with concurrent calculations of extracellular potentials and current dipole moments. The current dipole moments are then, in combination with suitable volume-conductor head models, used to compute non-invasive measures of neuronal activity, like scalp potentials (electroencephalographic recordings; EEG) and magnetic fields outside the head (magnetoencephalographic recordings; MEG). One such built-in head model is the four-sphere head model incorporating the different electric conductivities of brain, cerebrospinal fluid, skull and scalp. We demonstrate the new functionality of the software by constructing a network of biophysically detailed multicompartment neuron models from the Neocortical Microcircuit Collaboration (NMC) Portal (bbp.epfl.ch/nmc-portal) with corresponding statistics of connections and synapses, and compute in vivo-like extracellular potentials (local field potentials, LFP; electrocorticographical signals, ECoG) and corresponding current dipole moments. From the current dipole moments we estimate corresponding EEG and MEG signals using the four-sphere head model. We also show strong scaling performance of LFPy with different numbers of message-passing interface (MPI) processes, and for different network sizes with different density of connections. The open-source software LFPy is equally suitable for execution on laptops and in parallel on high-performance computing (HPC) facilities and is publicly available on GitHub.com.

Assessing the performance of Granger-Geweke causality: Benchmark dataset and simulation framework

Nonparametric methods based on spectral factorization offer well validated tools for estimating spectral measures of causality, called Granger–Geweke Causality (GGC). In Pagnotta et al. (2018) [1] we benchmarked nonparametric GGC methods using EEG data recorded during unilateral whisker stimulations in ten rats; here, we include detailed information about the benchmark dataset. In addition, we provide codes for estimating nonparametric GGC and a simulation framework to evaluate the effects on GGC analyses of potential problems, such as the common reference problem, signal-to-noise ratio (SNR) differences between channels, and the presence of additive noise. We focus on nonparametric methods here, but these issues also affect parametric methods, which can be tested in our framework as well. Our examples allow showing that time reversal testing for GGC (tr-GGC) mitigates the detrimental effects due to SNR imbalance and presence of mixed additive noise, and illustrate that, when using a common reference, tr-GGC unambiguously detects the causal influence׳s dominant spectral component, irrespective of the characteristics of the common reference signal. Finally, one of our simulations provides an example that nonparametric methods can overcome a pitfall associated with the implementation of conditional GGC in traditional parametric methods.

Spatio-temporal characteristics of population responses evoked by microstimulation in the barrel cortex

Intra-cortical microstimulation (ICMS) is a widely used technique to artificially stimulate cortical tissue. This method revealed functional maps and provided causal links between neuronal activity and cognitive, sensory or motor functions. The effects of ICMS on neural activity depend on stimulation parameters. Past studies investigated the effects of stimulation frequency mainly at the behavioral or motor level. Therefore the direct effect of frequency stimulation on the evoked spatio-temporal patterns of cortical activity is largely unknown. To study this question we used voltage-sensitive dye imaging to measure the population response in the barrel cortex of anesthetized rats evoked by high frequency stimulation (HFS), a lower frequency stimulation (LFS) of the same duration or a single pulse stimulation. We found that single pulse and short trains of ICMS induced cortical activity extending over few mm. HFS evoked a lower population response during the sustained response and showed a smaller activation across time and space compared with LFS. Finally the evoked population response started near the electrode site and spread horizontally at a propagation velocity in accordance with horizontal connections. In summary, HFS was less effective in cortical activation compared to LFS although HFS had 5 fold more energy than LFS.

Cortical signatures of wakeful somatosensory processing

Sensory inputs carry critical information for the survival of an organism. In mice, tactile information conveyed by the whiskers is of high behavioural relevance, and is broadcasted across cortical areas beyond the primary somatosensory cortex. Mesoscopic voltage sensitive dye imaging (VSDI) of cortical population response to whisker stimulations has shown that seemingly 'simple' sensory stimuli can have extended impact on cortical circuit dynamics. Here we took advantage of genetically encoded voltage indicators (GEVIs) that allow for cell type-specific monitoring of population voltage dynamics in a chronic dual-hemisphere transcranial windowed mouse preparation to directly compare the cortex-wide broadcasting of sensory information in wakening (lightly anesthetized to sedated) and awake mice. Somatosensory-evoked cortex-wide dynamics is altered across brain states, with anatomically sequential hyperpolarising activity observed in the awake cortex. GEVI imaging revealed cortical activity maps with increased specificity, high spatial coverage, and at the timescale of cortical information processing.

Excitable neuronal assemblies with adaptation as a building block of brain circuits for velocity-controlled signal propagation

The time scale of neuronal network dynamics is determined by synaptic interactions and neuronal signal integration, both of which occur on the time scale of milliseconds. Yet many behaviors like the generation of movements or vocalizations of sounds occur on the much slower time scale of seconds. Here we ask the question of how neuronal networks of the brain can support reliable behavior on this time scale. We argue that excitable neuronal assemblies with spike-frequency adaptation may serve as building blocks that can flexibly adjust the speed of execution of neural circuit function. We show in simulations that a chain of neuronal assemblies can propagate signals reliably, similar to the well-known synfire chain, but with the crucial difference that the propagation speed is slower and tunable to the behaviorally relevant range. Moreover we study a grid of excitable neuronal assemblies as a simplified model of the somatosensory barrel cortex of the mouse and demonstrate that various patterns of experimentally observed spatial activity propagation can be explained.

Models of activity propagation in cortical networks have often been based on feedforward structures. Here we propose a model of activity propagation, called excitation chain, which does not need such a feedforward structure. The model is composed of excitable neural assemblies with spike-frequency adaptation, connected bidirectionally in a row or a grid. This prototypical neural circuit can propagate activity forwards, backwards or in both directions. Furthermore, the propagation speed is slow enough to trigger the generation of behaviors on the time scale of hundreds of milliseconds. A two-dimensional variant of the model is able to generate different activity propagation patterns, similar to spontaneous activity and stimulus-evoked responses in anesthetized mouse barrel cortex. We propose the excitation chain model as a basic component that can be employed in various ways to create spiking neural circuit models that generate signals on behavioral time scales. In contrast to abstract models of excitable media, our model makes an explicit link to the time scale of neuronal spikes.

Optical study of interactions among propagation waves of neural excitation in the rat somatosensory cortex evoked by forelimb and hindlimb stimuli

Multisite optical recording has revealed that the neural excitation wave induced by a sensory stimulation begins at a focus and propagates in the cortex. This wave is considered to be important for computation in the sensory cortex, particularly the integration of sensory information; however, the nature of this wave remains largely unknown. In the present study, we examined the interaction between two waves in the rat sensory cortex induced by hindlimb and forelimb stimuli with different interstimulus intervals. We classified the resultant patterns as follows: 1) the collision of two waves, 2) the hindlimb response being evoked while the forelimb-induced wave is passing the hindlimb focus, and 3) the hindlimb response being evoked after the forelimb-induced wave has passed the hindlimb focus. In pattern 1, the two waves fused into a single wave, but the propagation pattern differed from that predicted by the superimposition of two singly induced propagation courses. In pattern 2, the state of the interaction between the two waves varied depending on the phase of optical signals constituting the forelimb-induced wave around the hindlimb focus. Although no hindlimb-induced wave was observed in the rising phase, the propagating velocity of the forelimb-induced wave increased. At the peak, neither the hindlimb-induced response nor a modulatory effect on the forelimb-induced wave was detected. In pattern 3, the hindlimb-induced wave showed a reduced amplitude and spatial extent. These results indicate that the state of the interaction between waves was strongly influenced by the relative timing of sensory inputs.

NEW & NOTEWORTHY Sensory stimulation-induced cortical excitation propagates as a wave and spreads over a wide area of the sensory cortex. To elucidate the characteristics of this relatively unknown phenomenon, we examined the interaction between two individually induced waves in the somatosensory cortex. Either the waves collided or the preceding wave affected the emergence of the following one. Our results indicate that the state of the interaction was strongly influenced by the relative timing of sensory inputs.

Diverse Long-Range Axonal Projections of Excitatory Layer 2/3 Neurons in Mouse Barrel Cortex

Excitatory projection neurons of the neocortex are thought to play important roles in perceptual and cognitive functions of the brain by directly connecting diverse cortical and subcortical areas. However, many aspects of the anatomical organization of these inter-areal connections are unknown. Here, we studied long-range axonal projections of excitatory layer 2/3 neurons with cell bodies located in mouse primary somatosensory barrel cortex (wS1). As a population, these neurons densely projected to secondary whisker somatosensory cortex (wS2) and primary/secondary whisker motor cortex (wM1/2), with additional axon in the dysgranular zone surrounding the barrel field, perirhinal temporal association cortex and striatum. In three-dimensional reconstructions of 6 individual wS2-projecting neurons and 9 individual wM1/2-projecting neurons, we found that both classes of neurons had extensive local axon in layers 2/3 and 5 of wS1. Neurons projecting to wS2 did not send axon to wM1/2, whereas a small subset of wM1/2-projecting neurons had relatively weak projections to wS2. A small fraction of projection neurons solely targeted wS2 or wM1/2. However, axon collaterals from wS2-projecting and wM1/2-projecting neurons were typically also found in subsets of various additional areas, including the dysgranular zone, perirhinal temporal association cortex and striatum. Our data suggest extensive diversity in the axonal targets selected by individual nearby cortical long-range projection neurons with somata located in layer 2/3 of wS1.

High-performance, inexpensive setup for simultaneous multisite recording of electrophysiological signals and mesoscale voltage imaging in the mouse cortex

Simultaneous recording of optical and electrophysiological signals from multiple cortical areas may provide crucial information to expand our understanding of cortical function. However, the insertion of multiple electrodes into the brain may compromise optical imaging by both restricting the field of view and interfering with the approaches used to stabilize the specimen. Existing methods that combine electrophysiological recording and optical imaging in vivo implement either multiple surface electrodes, silicon probes, or a single electrode for deeper recordings. To address such limitation, we built a microelectrode array (hyperdrive, patent US5928143 A) compatible with wide-field imaging that allows insertion of up to 12 probes into a large brain area (8 mm diameter). The hyperdrive is comprised of a circle of individual microdrives where probes are positioned at an angle leaving a large brain area unobstructed for wide-field imaging. Multiple tetrodes and voltage-sensitive dye imaging were used for acute simultaneous registration of spontaneous and evoked cortical activity in anesthetized mice. The electrophysiological signals were used to extract local field potential (LFP) traces, multiunit, and single-unit spiking activity. To demonstrate our approach, we compared LFP and VSD signals over multiple regions of the cortex and analyzed the relationship between single-unit and global cortical population activities. The study of the interactions between cortical activity at local and global scales, such as the one presented in this work, can help to expand our knowledge of brain function.

Low-intensity focused ultrasound alters the latency and spatial patterns of sensory-evoked cortical responses in vivo

OBJECTIVE

The use of transcranial, low intensity focused ultrasound (FUS) is an emerging neuromodulation technology that shows promise for both therapeutic and research applications. Among many, one of the most exciting applications is the use of FUS to rehabilitate or augment human sensory capabilities. While there is compelling empirical evidence demonstrating this capability, basic questions regarding the spatiotemporal extent of the modulatory effects remain. Our objective was to assess the basic, yet often overlooked hypothesis that FUS in fact alters sensory-evoked neural activity within the region of the cerebral cortex at the beam's focus.

APPROACH

To address this knowledge gap, we developed an approach to optically interrogate patterns of neural activity in the cortex directly at the acoustic focus, in vivo. Implementing simultaneous wide-field optical imaging and FUS stimulation in mice, our experiments probed somatosensory-evoked electrical activity through the use of voltage sensitive dyes (VSDs) and, in transgenic mice expressing GCaMP6f, monitored associated Ca²⁺ responses.

MAIN RESULTS

Our results demonstrate that low-intensity FUS alters both the kinetics and spatial patterns of neural activity in primary somatosensory cortex at the acoustic focus. When preceded by 1 s of pulsed ultrasound at intensities below 1 W cm^-2 (I _sppa), the onset of sensory-evoked cortical responses occurred 3.0  ±  0.7 ms earlier and altered the surface spatial morphology of Ca²⁺ responses.

SIGNIFICANCE

These findings support the heretofore unconfirmed assumption that FUS-induced sensory modulation reflects, at least in part, altered reactivity in primary sensory cortex at the site of sonication. The findings are significant given the interest in using FUS to target and alter spatial aspects of sensory receptive fields on the cerebral cortex.

Somatosensory maps

Somatosensory areas containing topographic maps of the body surface are a major feature of parietal cortex. In primates, parietal cortex contains four somatosensory areas, each with its own map, with the primary cutaneous map in area 3b. Rodents have at least three parietal somatosensory areas. Maps are not isomorphic to the body surface, but magnify behaviorally important skin regions, which include the hands and face in primates, and the whiskers in rodents. Within each map, intracortical circuits process tactile information, mediate spatial integration, and support active sensation. Maps may also contain fine-scale representations of touch submodalities, or direction of tactile motion. Functional representations are more overlapping than suggested by textbook depictions of map topography. The whisker map in rodent somatosensory cortex is a canonic system for studying cortical microcircuits, sensory coding, and map plasticity. Somatosensory maps are plastic throughout life in response to altered use or injury. This chapter reviews basic principles and recent findings in primate, human, and rodent somatosensory maps.

Seizure epicenter depth and translaminar field potential synchrony underlie complex variations in tissue oxygenation during ictal initiation

Whether functional hyperemia during epileptic activity is adequate to meet the heightened metabolic demand of such events is controversial. Whereas some studies have demonstrated hyperoxia during ictal onsets, other work has reported transient hypoxic episodes that are spatially dependent on local surface microvasculature. Crucially, how laminar differences in ictal evolution can affect subsequent cerebrovascular responses has not been thus far investigated, and is likely significant in view of possible laminar-dependent neurovascular mechanisms and angioarchitecture. We addressed this open question using a novel multi-modal methodology enabling concurrent measurement of cortical tissue oxygenation, blood flow and hemoglobin concentration, alongside laminar recordings of neural activity, in a urethane anesthetized rat model of recurrent seizures induced by 4-aminopyridine. We reveal there to be a close relationship between seizure epicenter depth, translaminar local field potential (LFP) synchrony and tissue oxygenation during the early stages of recurrent seizures, whereby deep layer seizures are associated with decreased cross laminar synchrony and prolonged periods of hypoxia, and middle layer seizures are accompanied by increased cross-laminar synchrony and hyperoxia. Through comparison with functional activation by somatosensory stimulation and graded hypercapnia, we show that these seizure-related cerebrovascular responses occur in the presence of conserved neural-hemodynamic and blood flow-volume coupling. Our data provide new insights into the laminar dependency of seizure-related neurovascular responses, which may reconcile inconsistent observations of seizure-related hypoxia in the literature, and highlight a potential layer-dependent vulnerability that may contribute to the harmful effects of clinical recurrent seizures. The relevance of our findings to perfusion-related functional neuroimaging techniques in epilepsy are also discussed.

Linear distributed source modeling of local field potentials recorded with intra-cortical electrode arrays

Planar intra-cortical electrode (Utah) arrays provide a unique window into the spatial organization of cortical activity. Reconstruction of the current source density (CSD) underlying such recordings, however, requires “inverting” Poisson’s equation. For inter-laminar recordings, this is commonly done by the CSD method, which consists in taking the second-order spatial derivative of the recorded local field potentials (LFPs). Although the CSD method has been tremendously successful in mapping the current generators underlying inter-laminar LFPs, its application to planar recordings is more challenging. While for inter-laminar recordings the CSD method seems reasonably robust against violations of its assumptions, is it unclear as to what extent this holds for planar recordings. One of the objectives of this study is to characterize the conditions under which the CSD method can be successfully applied to Utah array data. Using forward modeling, we find that for spatially coherent CSDs, the CSD method yields inaccurate reconstructions due to volume-conducted contamination from currents in deeper cortical layers. An alternative approach is to “invert” a constructed forward model. The advantage of this approach is that any a priori knowledge about the geometrical and electrical properties of the tissue can be taken into account. Although several inverse methods have been proposed for LFP data, the applicability of existing electroencephalographic (EEG) and magnetoencephalographic (MEG) inverse methods to LFP data is largely unexplored. Another objective of our study therefore, is to assess the applicability of the most commonly used EEG/MEG inverse methods to Utah array data. Our main conclusion is that these inverse methods provide more accurate CSD reconstructions than the CSD method. We illustrate the inverse methods using event-related potentials recorded from primary visual cortex of a macaque monkey during a motion discrimination task.

Laminar-specific encoding of texture elements in rat barrel cortex

KEY POINTS

For rats texture discrimination is signalled by the large face whiskers by stick-slip events. Neural encoding of repetitive stick-slip events will be influenced by intrinsic properties of adaptation. We show that texture coding in the barrel cortex is laminar specific and follows a power function. Our results also show layer 2 codes for novel feature elements via robust firing rates and temporal fidelity. We conclude that texture coding relies on a subtle neural ensemble to provide important object information.

ABSTRACT

Texture discrimination by rats is exquisitely guided by fine-grain mechanical stick-slip motions of the face whiskers as they encounter, stick to and slip past successive texture-defining surface features such as bumps and grooves. Neural encoding of successive stick-slip texture events will be shaped by adaptation, common to all sensory systems, whereby receptor and neural responses to a stimulus are affected by responses to preceding stimuli, allowing resetting to signal novel information. Additionally, when a whisker is actively moved to contact and brush over surfaces, that motion itself generates neural responses that could cause adaptation of responses to subsequent stick-slip events. Nothing is known about encoding in the rat whisker system of stick-slip events defining textures of different grain or the influence of adaptation from whisker protraction or successive texture-defining stick-slip events. Here we recorded responses from halothane-anaesthetized rats in response to texture-defining stimuli applied to passive whiskers. We demonstrate that: across the columnar network of the whisker-recipient barrel cortex, adaptation in response to repetitive stick-slip events is strongest in uppermost layers and equally lower thereafter; neither whisker protraction speed nor stick-slip frequency impede encoding of stick-slip events at rates up to 34.08 Hz; and layer 2 normalizes responses to whisker protraction to resist effects on texture signalling. Thus, within laminar-specific response patterns, barrel cortex reliably encodes texture-defining elements even to high frequencies.

Age-related changes in the spatiotemporal responses to electrical stimulation in the visual cortex of rats with progressive vision loss

The Royal College of Surgeons (RCS) rat gradually loses vision due to retinal degeneration. Previous physiological studies have depicted the progressive loss of optical responses in the visual pathway, including the primary visual cortex (V1), over the course of retinal degeneration in the RCS rat. However, little is known about how the excitability of the V1 circuit changes during over the course of the gradual loss of visual signal input from the retina. We elucidated the properties of responses to electrical stimulations directly applied to V1 at different stages of vision input loss in the RCS rat in reference to those of the Long-Evans (LE) rat, using in vivo voltage-sensitive dye imaging. The V1 neuronal network of the RCS rat exhibited an excitatory response comparable to the LE rat. The excitatory response was maintained even long after total loss of the visual signal input from the retina. However, the response time-course suggested that the suppressive response was somewhat debilitated in the RCS rat. This is the first experiment demonstrating the long-term effect of retinal degeneration on cortical activities. Our findings provide the physiological fundamentals to enhance the preclinical research of cortical prostheses with the use of the RCS rat.