Volume Conduction Coupling of Whisker-Evoked Cortical LFP in the Mouse Olfactory Bulb

Local contribution to the somatosensory evoked potentials in rat's thalamus

Local Field Potential (LFP), despite its name, often reflects remote activity. Depending on the orientation and synchrony of their sources, both oscillations and more complex waves may passively spread in brain tissue over long distances and be falsely interpreted as local activity at such distant recording sites. Here we show that the whisker-evoked potentials in the thalamic nuclei are of local origin up to around 6 ms post stimulus, but the later (7-15 ms) wave is overshadowed by a negative component reaching from cortex. This component can be analytically removed and local thalamic LFP can be recovered reliably using Current Source Density analysis. We used model-based kernel CSD (kCSD) method which allowed us to study the contribution of local and distant currents to LFP from rat thalamic nuclei and barrel cortex recorded with multiple, non-linear and non-regular multichannel probes. Importantly, we verified that concurrent recordings from the cortex are not essential for reliable thalamic CSD estimation. The proposed framework can be used to analyze LFP from other brain areas and has consequences for general LFP interpretation and analysis.

Respiratory rhythm modulates membrane potential and spiking of nonolfactory neurons

In recent years, several studies have shown a respiratory drive of the local field potential (LFP) in numerous brain areas so that the respiratory rhythm could be considered as a master clock promoting communication between distant brain locations. However, outside of the olfactory system, it remains unknown whether the respiratory rhythm could shape membrane potential (MP) oscillations. To fill this gap, we co-recorded MP and LFP activities in different nonolfactory brain areas, medial prefrontal cortex (mPFC), primary somatosensory cortex (S1), primary visual cortex (V1), and hippocampus (HPC), in urethane-anesthetized rats. Using respiratory cycle-by-cycle analysis, we observed that respiration could modulate both MP and spiking discharges in all recorded areas during episodes that we called respiration-related oscillations (RRo). Further quantifications revealed that RRo episodes were transient in most neurons (5 consecutive respiratory cycles in average). RRo development in MP was largely correlated with the presence of respiratory modulation in the LFP. By showing that the respiratory rhythm influenced brain activities deep to the MP of nonolfactory neurons, our data support the idea that respiratory rhythm could mediate long-range communication between brain areas.NEW & NOTEWORTHY In this study, we evidenced strong respiratory-driven oscillations of neuronal membrane potential and spiking discharge in various nonolfactory areas of the mammal brain. These oscillations were found in the medial prefrontal cortex, primary somatosensory cortex, primary visual cortex, and hippocampus. These findings support the idea that respiratory rhythm could be used as a common clock to set the dynamics of large-scale neuronal networks on the same slow rhythm.

Extracting electromyographic signals from multi-channel LFPs using independent component analysis without direct muscular recording

Electromyography (EMG) has been commonly used for the precise identification of animal behavior. However, it is often not recorded together with in vivo electrophysiology due to the need for additional surgeries and setups and the high risk of mechanical wire disconnection. While independent component analysis (ICA) has been used to reduce noise from field potential data, there has been no attempt to proactively use the removed "noise," of which EMG signals are thought to be one of the major sources. Here, we demonstrate that EMG signals can be reconstructed without direct EMG recording using the "noise" ICA component from local field potentials. The extracted component is highly correlated with directly measured EMG, termed IC-EMG. IC-EMG is useful for measuring an animal's sleep/wake, freezing response, and non-rapid eye movement (NREM)/REM sleep states consistently with actual EMG. Our method has advantages in precise and long-term behavioral measurement in wide-ranging in vivo electrophysiology experiments.

Theoretical considerations and supporting evidence for the primary role of source geometry on field potential amplitude and spatial extent

Field potential (FP) recording is an accessible means to capture the shifts in the activity of neuron populations. However, the spatial and composite nature of these signals has largely been ignored, at least until it became technically possible to separate activities from co-activated sources in different structures or those that overlap in a volume. The pathway-specificity of mesoscopic sources has provided an anatomical reference that facilitates transcending from theoretical analysis to the exploration of real brain structures. We review computational and experimental findings that indicate how prioritizing the spatial geometry and density of sources, as opposed to the distance to the recording site, better defines the amplitudes and spatial reach of FPs. The role of geometry is enhanced by considering that zones of the active populations that act as sources or sinks of current may arrange differently with respect to each other, and have different geometry and densities. Thus, observations that seem counterintuitive in the scheme of distance-based logic alone can now be explained. For example, geometric factors explain why some structures produce FPs and others do not, why different FP motifs generated in the same structure extend far while others remain local, why factors like the size of an active population or the strong synchronicity of its neurons may fail to affect FPs, or why the rate of FP decay varies in different directions. These considerations are exemplified in large structures like the cortex and hippocampus, in which the role of geometrical elements and regional activation in shaping well-known FP oscillations generally go unnoticed. Discovering the geometry of the sources in play will decrease the risk of population or pathway misassignments based solely on the FP amplitude or temporal pattern.

Localized chemogenetic silencing of inhibitory neurons: a novel mouse model of focal cortical epileptic activity

Focal cortical epilepsies are frequently refractory to available anticonvulsant drug therapies. One key factor contributing to this state is the limited availability of animal models that allow to reliably study focal cortical seizures and how they recruit surrounding brain areas in vivo. In this study, we selectively expressed the inhibitory chemogenetic receptor, hM4D, in GABAergic neurons in focal cortical areas using viral gene transfer. GABAergic silencing using Clozapine-N-Oxide (CNO) demonstrated reliable induction of local epileptiform events in the electroencephalogram signal of awake freely moving mice. Anesthetized mice experiments showed consistent induction of focal epileptiform-events in both the barrel cortex (BC) and the medial prefrontal cortex (mPFC), accompanied by high-frequency oscillations, a known characteristic of human seizures. Epileptiform-events showed propagation indication with favored propagation pathways: from the BC on 1 hemisphere to its counterpart and from the BC to the mPFC, but not vice-versa. Lastly, sensory whisker-pad stimulation evoked BC epileptiform events post-CNO, highlighting the potential use of this model in studying sensory-evoked seizures. Combined, our results show that targeted chemogenetic inhibition of GABAergic neurons using hM4D can serve as a novel, versatile, and reliable model of focal cortical epileptic activity suitable for systematically studying cortical ictogenesis in different cortical areas.

High-frequency neuronal signal better explains multi-phase BOLD response

Visual stimulation-evoked blood-oxygen-level dependent (BOLD) responses can exhibit more complex temporal dynamics than a simple monophasic response. For instance, BOLD responses sometimes include a phase of positive response followed by a phase of post-stimulus undershoot. Whether the BOLD response during these phases reflects the underlying neuronal signal fluctuations or is contributed by non-neuronal physiological factors remains elusive. When presenting blocks of sustained (i.e. DC) light ON-OFF stimulations to unanesthetized rats, we observed that the response following a decrease in illumination (i.e. OFF stimulation-evoked BOLD response) in the visual cortices displayed reproducible multiple phases, including an initial positive BOLD response, followed by an undershoot and then an overshoot before the next ON trial. This multi-phase BOLD response did not result from the entrainment of the periodic stimulation structure. When we measured the neural correlates of these responses, we found that the high-frequency band from the LFP power (300 - 3000 Hz, multi-unit activity (MUA)), but not the power in the gamma band (30 - 100 Hz) exhibited the same multiphasic dynamics as the BOLD signal. This study suggests that the post-stimulus phases of the BOLD response can be better explained by the high-frequency neuronal signal.

Neural underpinning of a respiration-associated resting-state fMRI network

Respiration can induce motion and CO2 fluctuation during resting-state fMRI (rsfMRI) scans, which will lead to non-neural artifacts in the rsfMRI signal. In the meantime, as a crucial physiologic process, respiration can directly drive neural activity change in the brain, and may thereby modulate the rsfMRI signal. Nonetheless, this potential neural component in the respiration-fMRI relationship is largely unexplored. To elucidate this issue, here we simultaneously recorded the electrophysiology, rsfMRI, and respiration signals in rats. Our data show that respiration is indeed associated with neural activity changes, evidenced by a phase-locking relationship between slow respiration variations and the gamma-band power of the electrophysiological signal recorded in the anterior cingulate cortex. Intriguingly, slow respiration variations are also linked to a characteristic rsfMRI network, which is mediated by gamma-band neural activity. In addition, this respiration-related brain network disappears when brain-wide neural activity is silenced at an isoelectrical state, while the respiration is maintained, further confirming the necessary role of neural activity in this network. Taken together, this study identifies a respiration-related brain network underpinned by neural activity, which represents a novel component in the respiration-rsfMRI relationship that is distinct from respiration-related rsfMRI artifacts. It opens a new avenue for investigating the interactions between respiration, neural activity, and resting-state brain networks in both healthy and diseased conditions.

In vivo direct imaging of neuronal activity at high temporospatial resolution

There has been a long-standing demand for noninvasive neuroimaging methods that can detect neuronal activity at both high temporal and high spatial resolution. We present a two-dimensional fast line-scan approach that enables direct imaging of neuronal activity with millisecond precision while retaining the high spatial resolution of magnetic resonance imaging (MRI). This approach was demonstrated through in vivo mouse brain imaging at 9.4 tesla during electrical whisker-pad stimulation. In vivo spike recording and optogenetics confirmed the high correlation of the observed MRI signal with neural activity. It also captured the sequential and laminar-specific propagation of neuronal activity along the thalamocortical pathway. This high-resolution, direct imaging of neuronal activity will open up new avenues in brain science by providing a deeper understanding of the brain's functional organization, including the temporospatial dynamics of neural networks.

Site-dependent shaping of field potential waveforms

The activity of neuron populations gives rise to field potentials (FPs) that extend beyond the sources. Their mixing in the volume dilutes the original temporal motifs in a site-dependent manner, a fact that has received little attention. And yet, it potentially rids of physiological significance the time-frequency parameters of individual waves (amplitude, phase, duration). This is most likely to happen when a single source or a local origin is erroneously assumed. Recent studies using spatial treatment of these signals and anatomically realistic modeling of neuron aggregates provide convincing evidence for the multisource origin and site-dependent blend of FPs. Thus, FPs generated in primary structures like the neocortex and hippocampus reach far and cross-contaminate each other but also, they add and even impose their temporal traits on distant regions. Furthermore, both structures house neurons that act as spatially distinct (but overlapped) FP sources whose activation is state, region, and time dependent, making the composition of so-called local FPs highly volatile and strongly site dependent. Since the spatial reach cannot be predicted without source geometry, it is important to assess whether waveforms and temporal motifs arise from a single source; otherwise, those from each of the co-active sources should be sought.

Cross Laminar Traveling Components of Field Potentials due to Volume Conduction of Non-Traveling Neuronal Activity in Macaque Sensory Cortices

Field potentials (FPs) reflect neuronal activities in the brain, and often exhibit traveling peaks across recording sites. While traveling FPs are interpreted as propagation of neuronal activity, not all studies directly reveal such propagating patterns of neuronal activation. Neuronal activity is associated with transmembrane currents that form dipoles and produce negative and positive fields. Thereby, FP components reverse polarity between those fields and have minimal amplitudes at the center of dipoles. Although their amplitudes could be smaller, FPs are never flat even around these reversals. What occurs around the reversal has not been addressed explicitly, although those are rationally in the middle of active neurons. We show that sensory FPs around the reversal appeared with peaks traveling across cortical laminae in macaque sensory cortices. Interestingly, analyses of current source density did not depict traveling patterns but lamina-delimited current sinks and sources. We simulated FPs produced by volume conduction of a simplified 2 dipoles' model mimicking sensory cortical laminar current source density components. While FPs generated by single dipoles followed the temporal patterns of the dipole moments without traveling peaks, FPs generated by concurrently active dipole moments appeared with traveling components in the vicinity of dipoles by superimposition of individually non-traveling FPs generated by single dipoles. These results indicate that not all traveling FP are generated by traveling neuronal activity, and that recording positions need to be taken into account to describe FP peak components around active neuronal populations.SIGNIFICANCE STATEMENT Field potentials (FPs) generated by neuronal activity in the brain occur with fields of opposite polarity. Likewise, in the cerebral cortices, they have mirror-imaged waveforms in upper and lower layers. We show that FPs appear like traveling across the cortical layers. Interestingly, the traveling FPs occur without traveling components of current source density, which represents transmembrane currents associated with neuronal activity. These seemingly odd findings are explained using current source density models of multiple dipoles. Concurrently active, non-traveling dipoles produce FPs as mixtures of FPs produced by individual dipoles, and result in traveling FP waveforms as the mixing ratio depends on the distances from those dipoles. The results suggest that not all traveling FP components are associated with propagating neuronal activity.

Decision Signals in the Local Field Potentials of Early and Mid-Level Macaque Visual Cortex

The neural basis of perceptual decision making has typically been studied using measurements of single neuron activity, though decisions are likely based on the activity of large neuronal ensembles. Local field potentials (LFPs) may, in some cases, serve as a useful proxy for population activity and thus be useful for understanding the neural basis of perceptual decision making. However, little is known about whether LFPs in sensory areas include decision-related signals. We therefore analyzed LFPs recorded using two 48-electrode arrays implanted in primary visual cortex (V1) and area V4 of macaque monkeys trained to perform a fine orientation discrimination task. We found significant choice information in low (0-30 Hz) and higher (70-500 Hz) frequency components of the LFP, but little information in gamma frequencies (30-70 Hz). Choice information was more robust in V4 than V1 and stronger in LFPs than in simultaneously measured spiking activity. LFP-based choice information included a global component, common across electrodes within an area. Our findings reveal the presence of robust choice-related signals in the LFPs recorded in V1 and V4 and suggest that LFPs may be a useful complement to spike-based analyses of decision making.

Nasal respiration is necessary for ketamine-dependent high frequency network oscillations and behavioral hyperactivity in rats

Changes in oscillatory activity are widely reported after subanesthetic ketamine, however their mechanisms of generation are unclear. Here, we tested the hypothesis that nasal respiration underlies the emergence of high-frequency oscillations (130–180 Hz, HFO) and behavioral activation after ketamine in freely moving rats. We found ketamine 20 mg/kg provoked “fast” theta sniffing in rodents which correlated with increased locomotor activity and HFO power in the OB. Bursts of ketamine-dependent HFO were coupled to “fast” theta frequency sniffing. Theta coupling of HFO bursts were also found in the prefrontal cortex and ventral striatum which, although of smaller amplitude, were coherent with OB activity. Haloperidol 1 mg/kg pretreatment prevented ketamine-dependent increases in fast sniffing and instead HFO coupling to slower basal respiration. Consistent with ketamine-dependent HFO being driven by nasal respiration, unilateral naris blockade led to an ipsilateral reduction in ketamine-dependent HFO power compared to the control side. Bilateral nares blockade reduced ketamine-induced hyperactivity and HFO power and frequency. These findings suggest that nasal airflow entrains ketamine-dependent HFO in diverse brain regions, and that the OB plays an important role in the broadcast of this rhythm.

Local and Volume-Conducted Contributions to Cortical Field Potentials

Brain field potentials (FPs) can reach far from their sources, making difficult to know which waves come from where. We show that modern algorithms efficiently segregate the local and remote contributions to cortical FPs by recovering the generator-specific spatial voltage profiles. We investigated experimentally and numerically the local and remote origin of FPs in different cortical areas in anesthetized rats. All cortices examined show significant state, layer, and region dependent contribution of remote activity, while the voltage profiles help identify their subcortical or remote cortical origin. Co-activation of different cortical modules can be discriminated by the distinctive spatial features of the corresponding profiles. All frequency bands contain remote activity, thus influencing the FP time course, in cases drastically. The reach of different FP patterns is boosted by spatial coherence and curved geometry of the sources. For instance, slow cortical oscillations reached the entire brain, while hippocampal theta reached only some portions of the cortex. In anterior cortices, most alpha oscillations have a remote origin, while in the visual cortex the remote theta and gamma even surpass the local contribution. The quantitative approach to local and distant FP contributions helps to refine functional connectivity among cortical regions, and their relation to behavior.

nNOS-expressing interneurons control basal and behaviorally evoked arterial dilation in somatosensory cortex of mice

Cortical neural activity is coupled to local arterial diameter and blood flow. However, which neurons control the dynamics of cerebral arteries is not well understood. We dissected the cellular mechanisms controlling the basal diameter and evoked dilation in cortical arteries in awake, head-fixed mice. Locomotion drove robust arterial dilation, increases in gamma band power in the local field potential (LFP), and increases calcium signals in pyramidal and neuronal nitric oxide synthase (nNOS)-expressing neurons. Chemogenetic or pharmocological modulation of overall neural activity up or down caused corresponding increases or decreases in basal arterial diameter. Modulation of pyramidal neuron activity alone had little effect on basal or evoked arterial dilation, despite pronounced changes in the LFP. Modulation of the activity of nNOS-expressing neurons drove changes in the basal and evoked arterial diameter without corresponding changes in population neural activity.

Layer 4 organization and respiration locking in the rodent nose somatosensory cortex

Rodents and other mammals acquire sensory information by precisely orchestrated head, whisker, and respiratory movements. We have, however, only limited information about integration of these signals. In the somatosensory domain, the integration of somatosensory information with other modalities is particularly pertinent for body parts such as eyes, ears, and nose, which serve another modality. Here we analyzed the nose/nostril representation in the rodent somatosensory cortex. We identified the representation of the nose/nostril in the rat somatosensory cortex by receptive field mapping and subsequent histological reconstruction. In tangential somatosensory cortical sections, the rat nostril cortex was evident as a prominent stripe-like recess of layer 4 revealed by cytochrome-c oxidase reactivity or by antibodies against the vesicular glutamate-transporter-2 (identifying thalamic afferents). We compared flattened somatosensory cortices of various rodents including rats, mice, gerbils, chinchillas, and chipmunks. We found that such a nose/nostril module was evident as a region with thinned or absent layer 4 at the expected somatotopic position of the nostril. Extracellular spike activity was strongly modulated by respiration in the rat somatosensory cortex, and field potential recordings revealed a stronger locking of nostril recording sites to respiration than for whisker/barrel cortex recoding sites. We conclude that the rodent nose/nostril representation has a conserved architecture and specifically interfaces with respiration signals.NEW & NOTEWORTHY We characterized the rodent nose somatosensory cortex. The nostril representation appeared as a kind of "hole" (i.e., as a stripe-like recess of layer 4) in tangential cortical sections. Neural activity in nose somatosensory cortex was locked to respiration, and simultaneous field recordings indicate that this locking was specific to this region. Our results reveal previously unknown cytoarchitectonic and physiological properties of the rodent nose somatosensory cortex, potentially enabling it to integrate multiple sensory modalities.

Volume-Conducted Origin of the Field Potential at the Lateral Habenula

Field potentials (FPs) are easily reached signals that provide information about the brain's processing. However, FP should be interpreted cautiously since their biophysical bases are complex. The lateral habenula (LHb) is a brain structure involved in the encoding of aversive motivational values. Previous work indicates that the activity of the LHb is relevant for hippocampal-dependent learning. Moreover, it has been proposed that the interaction of the LHb with the hippocampal network is evidenced by the synchronization of LHb and hippocampal FPs during theta rhythm. However, the origin of the habenular FP has not been analyzed. Hence, its validity as a measurement of LHb activity has not been proven. In this work, we used electrophysiological recordings in anesthetized rats and feed-forward modeling to investigate biophysical basis of the FP recorded in the LHb. Our results indicate that the FP in the LHb during theta rhythm is a volume-conducted signal from the hippocampus. This result highlight that FPs must be thoroughly analyzed before its biological interpretation and argues against the use of the habenular FP signal as a readout of the activity of the LHb.

Twitches, Blinks, and Fidgets: Important Generators of Ongoing Neural Activity

Animals and humans continuously engage in small, spontaneous motor actions, such as blinking, whisking, and postural adjustments ("fidgeting"). These movements are accompanied by changes in neural activity in sensory and motor regions of the brain. The frequency of these motions varies in time, is affected by sensory stimuli, arousal levels, and pathology. These fidgeting behaviors can be entrained by sensory stimuli. Fidgeting behaviors will cause distributed, bilateral functional activation in the 0.01 to 0.1 Hz frequency range that will show up in functional magnetic resonance imaging and wide-field calcium neuroimaging studies, and will contribute to the observed functional connectivity among brain regions. However, despite the large potential of these behaviors to drive brain-wide activity, these fidget-like behaviors are rarely monitored. We argue that studies of spontaneous and evoked brain dynamics in awake animals and humans should closely monitor these fidgeting behaviors. Differences in these fidgeting behaviors due to arousal or pathology will "contaminate" ongoing neural activity, and lead to apparent differences in functional connectivity. Monitoring and accounting for the brain-wide activations by these behaviors is essential during experiments to differentiate fidget-driven activity from internally driven neural dynamics.

Vascular and neural basis of the BOLD signal

Neural activity in the brain is usually coupled to increases in local cerebral blood flow, leading to the increase in oxygenation that generates the BOLD fMRI signal. Recent work has begun to elucidate the vascular and neural mechanisms underlying the BOLD signal. The dilatory response is distributed throughout the vascular network. Arteries actively dilate within a second following neural activity increases, while venous distensions are passive and have a time course that last tens of seconds. Vasodilation, and thus local blood flow, is controlled by the activity of both neurons and astrocytes via multiple different pathways. The relationship between sensory-driven neural activity and the vascular dynamics in sensory areas are well-captured with a linear convolution model. However, depending on the behavioral state or brain region, the coupling between neural activity and hemodynamic signals can be weak or even inverted.

Correlation Structure in Micro-ECoG Recordings is Described by Spatially Coherent Components

Electrocorticography (ECoG) is becoming more prevalent due to improvements in fabrication and recording technology as well as its ease of implantation compared to intracortical electrophysiology, larger cortical coverage, and potential advantages for use in long term chronic implantation. Given the flexibility in the design of ECoG grids, which is only increasing, it remains an open question what geometry of the electrodes is optimal for an application. Conductive polymer, PEDOT:PSS, coated microelectrodes have an advantage that they can be made very small without losing low impedance. This makes them suitable for evaluating the required granularity of ECoG recording in humans and experimental animals. We used two-dimensional (2D) micro-ECoG grids to record intra-operatively in humans and during acute implantations in mouse with separation distance between neighboring electrodes (i.e., pitch) of 0.4 mm and 0.2/0.25 mm respectively. To assess the spatial properties of the signals, we used the average correlation between electrodes as a function of the pitch. In agreement with prior studies, we find a strong frequency dependence in the spatial scale of correlation. By applying independent component analysis (ICA), we find that the spatial pattern of correlation is largely due to contributions from multiple spatially extended, time-locked sources present at any given time. Our analysis indicates the presence of spatially structured activity down to the sub-millimeter spatial scale in ECoG despite the effects of volume conduction, justifying the use of dense micro-ECoG grids.

Electrocorticography (ECoG) is a type of electrophysiological monitoring that uses electrodes placed directly on the exposed surface of the brain. ECoG is a promising technique for studying the brain, and EcoG signals can be used to control brain-computer interfaces. Advances have made it possible to record simultaneously with an increasing number of smaller, and more closely spaced electrodes. However, a property of electrical recording from outside the brain is that common signals appear on different electrodes at different locations, and this affects decisions about how to best distribute a limited number of electrodes to maximize the information that can be gathered. Large spacing of electrodes around one centimeter apart on the brain’s surface has proven useful for clinical and research use, but how much benefit there is to recording from more locations in a smaller area remains to be answered. We found that we can explain the commonality between the different locations as the combination of different patterns of brain activity that are present at multiple electrode locations, and that signals recorded from very closely spaced electrodes, around a millimeter or less apart, are able to identify patterns that are at this small scale.

Differential recordings of local field potential: A genuine tool to quantify functional connectivity

Local field potential (LFP) recording is a very useful electrophysiological method to study brain processes. However, this method is criticized for recording low frequency activity in a large area of extracellular space potentially contaminated by distal activity. Here, we theoretically and experimentally compare ground-referenced (RR) with differential recordings (DR). We analyze electrical activity in the rat cortex with these two methods. Compared with RR, DR reveals the importance of local phasic oscillatory activities and their coherence between cortical areas. Finally, we show that DR provides a more faithful assessment of functional connectivity caused by an increase in the signal to noise ratio, and of the delay in the propagation of information between two cortical structures.