LFP publications

Advancements in the field of implantable neurotechnologies have enabled the integration of hundreds of microelectrodes on ultra-thin and flexible substrates. Besides implantable components, also connectors, headstages and cables have to comply with the high-count demand, resulting in a complex and compact chain with reduced line spacing and smaller safety margins. Here, we show that epicortical recordings acquired from anesthetized rat brains with a state-of-art neural acquisition system are undoubtedly compromised by crosstalk, with signal coherence maps exhibiting a strong dependency to the routing layout. A crosstalk back-correction algorithm is developed, allowing to infer on how signals would look like under a zero-crosstalk scenario. We found that signal coherence between closely routed channels effectively drops after correction, corroborating crosstalk contamination. Our work stresses the importance of validating recorded data against the routing layout as a crucial step of data quality control, helping to come closer to ground truth data.

Closed-loop neural control is a powerful tool for both the scientific exploration of neural function and for mitigating deficiencies found in open-loop deep brain stimulation (DBS). Here, we present a protocol for artificial intelligence-guided neural control in rats using deep reinforcement learning (RL) and infrared neural stimulation (INS). We describe steps for integrating RL closed-loop control into neuroscience and neuromodulation studies. We then detail procedures for using and deploying infrared INS in chronic DBS applications. For complete details on the use and execution of this protocol, please refer to Coventry et al.1 and Coventry and Bartlett.2.

To comprehend the volumetric neural connectivity of a brain organoid, it is crucial to monitor the spatiotemporal electrophysiological signals within the organoid, known as intra-organoid signals. However, previous methods risked damaging the three-dimensional (3D) cytoarchitecture of organoids, either through sectioning or inserting rigid needle-like electrodes. Also, the limited numbers of electrodes in fixed positions with non-adjustable electrode shapes were insufficient for examining the complex neural activity throughout the organoid. Herein, we present a magnetically reshapable 3D multi-electrode array (MEA) using direct printing of liquid metals for electrophysiological analysis of brain organoids. The adaptable distribution and the softness of these printed electrodes facilitate the spatiotemporal recording of intra-organoid signals. Furthermore, the unique capability to reshape these soft electrodes within the organoid using magnetic fields allows a single electrode in the MEA to record from multiple points, effectively increasing the recording site density without the need for additional electrodes.

In this paper we aim to examine the correlation between diffusion tensor imaging parameters of anatomical connectivity and characteristics of signal transmission obtained from patient-specific transfer function models. Here, we focused on elucidating the correlation between structural and functional neural connectivity within a cohort of pediatric patients diagnosed with dystonia. Diffusion tractography images were obtained from 12 patients with dystonia prior to the deep brain stimulation surgery. For each patient, we processed the imaging data to estimate anatomical measures including fractional anisotropy, axial diffusivity, number of fibre tracts per unit area, and fibre tract length. After the implantation of temporary depth leads for each patient as part of their treatment plan, intracranial signals were recorded. Transfer function models of local field potential recordings and the corresponding measures of functional connectivity were computed for each patient. Linear mixed effect analysis was then employed to determine the relationship between transfer function measures and diffusion tractography parameters. Our results illustrate a positive correlation between fractional anisotropy, AD, and intrinsic neural transmission measures, representing amplification and spread of intrinsic neural oscillations, obtained from the transfer functions models. However, no significant correlation was found between the functional connectivity and number of fibre tracts or fibre lengths. Our findings suggest that white matter integrity, as measured by fractional anisotropy and AD, can potentially reflect the amplification and spread of intrinsic brain signals throughout the network. This study underscores the significant relationship between structural and functional connectivity, offering valuable insights into propagation of neural activity in the brain network and potential implications for optimizing non-invasive treatments and planning for neurological disorders.

How consciousness is lost in states such as sleep or anesthesia remains a mystery. To gain insight into this phenomenon, concurrent recordings of electrophysiology signals in the anterior cingulate cortex and whole‐brain functional magnetic resonance imaging (fMRI) are conducted in rats exposed to graded propofol, undergoing the transition from consciousness to unconsciousness. The results reveal that upon the loss of consciousness (LOC), there is a sharp increase in low‐frequency power of the electrophysiological signal. Additionally, fMRI signals exhibit a cascade of deactivation across a pathway including the hippocampus, thalamus, and medial prefrontal cortex (mPFC) surrounding the moment of LOC, followed by a broader increase in brain activity across the cortex during sustained unconsciousness. Furthermore, sliding window analysis demonstrates a temporary increase in synchrony of fMRI signals across the hippocampus‐thalamus‐mPFC pathway preceding LOC. These data suggest that LOC may be triggered by sequential activities in the hippocampus, thalamus, and mPFC, while wide‐spread activity increases in other cortical regions commonly observed during anesthesia‐induced unconsciousness may be a consequence, rather than a cause of LOC. Taken together, the study identifies a cascade of neural events unfolding as the brain transitions into unconsciousness, offering insight into the systems‐level neural mechanisms underpinning LOC.

Here, we characterized the dynamics of sleep spindles, focusing on their damping, which we estimated using a metric called oscillatory-Quality (o-Quality), derived by fitting an autoregressive model to electrophysiological signals, recorded from the cortex in mice. The o-Quality of sleep spindles correlates weakly with their amplitude, shows marked laminar differences and regional topography across cortical regions, reflects the level of synchrony within and between cortical networks, is strongly modulated by sleep-wake history, reflects the degree of sensory disconnection, and correlates with the strength of coupling between spindles and slow waves. As most spindle events are highly localized and not detectable with conventional low-density recording approaches, o-Quality thus emerges as a valuable metric that allows us to infer the spread and dynamics of spindle activity across the brain and directly links their spatiotemporal dynamics with local and global regulation of brain states, sleep regulation, and function.

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by ‘electrophysiology-invisible’ signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Study of the foreign body reaction to implanted electrodes in the brain is an important area of research for the future development of neuroprostheses and experimental electrophysiology. After electrode implantation in the brain, microglial activation, reactive astrogliosis, and neuronal cell death create an environment immediately surrounding the electrode that is significantly altered from its homeostatic state. Objective. To uncover physiological changes potentially affecting device function and longevity, spatial transcriptomics (ST) was implemented to identify changes in gene expression driven by electrode implantation and compare this differential gene expression to traditional metrics of glial reactivity, neuronal loss, and electrophysiological recording quality. Approach. For these experiments, rats were chronically implanted with functional Michigan-style microelectrode arrays, from which electrophysiological recordings (multi-unit activity, local field potential) were taken over a six-week time course. Brain tissue cryosections surrounding each electrode were then mounted for ST processing. The tissue was immunolabeled for neurons and astrocytes, which provided both a spatial reference for ST and a quantitative measure of glial fibrillary acidic protein and neuronal nuclei immunolabeling surrounding each implant. Main results . Results from rat motor cortex within 300 µ m of the implanted electrodes at 24 h, 1 week, and 6 weeks post-implantation showed up to 553 significantly differentially expressed (DE) genes between implanted and non-implanted tissue sections. Regression on the significant DE genes identified the 6–7 genes that had the strongest relationship to histological and electrophysiological metrics, revealing potential candidate biomarkers of recording quality and the tissue response to implanted electrodes. Significance . Our analysis has shed new light onto the potential mechanisms involved in the tissue response to implanted electrodes while generating hypotheses regarding potential biomarkers related to recorded signal quality. A new approach has been developed to understand the tissue response to electrodes implanted in the brain using genes identified through transcriptomics, and to screen those results for potential relationships with functional outcomes.

Collective studies have demonstrated that transcranial ultrasound stimulation (TUS) can elicit activation in hemodynamics, implying its potential in treating cerebral or peripheral vessel-related malfunction. The theory for hemodynamic response to TUS is neurovascular coupling (NVC) following the ultrasound-induced cellular (de)polarization. However, it was not conclusive due to the co-existence of the pathway of direct ultrasound-vessel interactions. This study thus aims to investigate and provide direct evidence for NVC pathway in a rodent model of TUS by inhibiting neural activity with sodium valproate (VPA), a GABAergic agent. Twenty Sprague-Dawley rats were randomly assigned to VPA and Saline groups. Microelectrode and optical imaging were utilized to record local field potential and relative cerebral blood flow (rCBF) during baseline, before, and after TUS periods. We found the attenuated neural activity was associated with reduced rCBF responses. These results provided direct evidence for the presence of NVC pathway in hemodynamic modulation by TUS.

Typical statistical practices in the biological sciences have been increasingly called into question due to difficulties in the replication of an increasing number of studies, many of which are confounded by the relative difficulty of null significance hypothesis testing designs and interpretation of p -values. Bayesian inference, representing a fundamentally different approach to hypothesis testing, is receiving renewed interest as a potential alternative or complement to traditional null significance hypothesis testing due to its ease of interpretation and explicit declarations of prior assumptions. Bayesian models are more mathematically complex than equivalent frequentist approaches, which have historically limited applications to simplified analysis cases. However, the advent of probability distribution sampling tools with exponential increases in computational power now allows for quick and robust inference under any distribution of data. Here we present a practical tutorial on the use of Bayesian inference in the context of neuroscientific studies in both rat electrophysiological and computational modeling data. We first start with an intuitive discussion of Bayes' rule and inference followed by the formulation of Bayesian-based regression and ANOVA models using data from a variety of neuroscientific studies. We show how Bayesian inference leads to easily interpretable analysis of data while providing an open-source toolbox to facilitate the use of Bayesian tools.

The cerebral cortex is widely considered part of the neural substrate of consciousness, but direct causal evidence is missing. Here, we tested in mice whether optogenetic activation of cortical neurons in posterior parietal cortex (PtA) or medial prefrontal cortex (mPFC) is sufficient for arousal from three behavioral states characterized by progressively deeper unresponsiveness: sleep, a coma-like state induced by muscimol injection in the midbrain, and deep sevoflurane-dexmedetomidine anesthesia. We find that cortical stimulation always awakens the mice from both NREM sleep and REM sleep, with PtA requiring weaker/shorter light pulses than mPFC. Moreover, in most cases light pulses produce both cortical activation (decrease in low frequencies) and behavioral arousal (recovery of the righting reflex) from brainstem coma, as well as cortical activation from anesthesia. These findings provide evidence that direct activation of cortical neurons is sufficient for behavioral and/or cortical arousal from sleep, brainstem coma, and anesthesia.

Objective . Therapeutic brain stimulation is conventionally delivered using constant-frequency stimulation pulses. Several recent clinical studies have explored how unconventional and irregular temporal stimulation patterns could enable better therapy. However, it is challenging to understand which irregular patterns are most effective for different therapeutic applications given the massively high-dimensional parameter space. Approach . Here we applied many irregular stimulation patterns in a single neural circuit to demonstrate how they can enable new dimensions of neural control compared to conventional stimulation, to guide future exploration of novel stimulation patterns in translational settings. We optogenetically excited the septohippocampal circuit with constant-frequency, nested pulse, sinusoidal, and randomized stimulation waveforms, systematically varying their amplitude and frequency parameters. Main results. We first found equal entrainment of hippocampal oscillations: all waveforms provided similar gamma-power increase, whereas no parameters increased theta-band power above baseline (despite the mechanistic role of the medial septum in driving hippocampal theta oscillations). We then compared each of the effects of each waveform on high-dimensional multi-band activity states using dimensionality reduction methods. Strikingly, we found that conventional stimulation drove predominantly ‘artificial’ (different from behavioral activity) effects, whereas all irregular waveforms induced activity patterns that more closely resembled behavioral activity. Significance . Our findings suggest that irregular stimulation patterns are not useful when the desired mechanism is to suppress or enhance a single frequency band. However, novel stimulation patterns may provide the greatest benefit for neural control applications where entraining a particular mixture of bands (e.g. if they are associated with different symptoms) or behaviorally-relevant activity is desired.

Objective. This study aims to reveal longitudinal changes in functional network connectivity within and across different brain structures near chronically implanted microelectrodes. While it is well established that the foreign-body response (FBR) contributes to the gradual decline of the signals recorded from brain implants over time, how the FBR affects the functional stability of neural circuits near implanted brain–computer interfaces (BCIs) remains unknown. This research aims to illuminate how the chronic FBR can alter local neural circuit function and the implications for BCI decoders. Approach. This study utilized single-shank, 16-channel,100 µ m site-spacing Michigan-style microelectrodes (3 mm length, 703 µ m2 site area) that span all cortical layers and the hippocampal CA1 region. Sex balanced C57BL6 wildtype mice (11–13 weeks old) received perpendicularly implanted microelectrode in left primary visual cortex. Electrophysiological recordings were performed during both spontaneous activity and visual sensory stimulation. Alterations in neuronal activity near the microelectrode were tested assessing cross-frequency synchronization of local field potential (LFP) and spike entrainment to LFP oscillatory activity throughout 16 weeks after microelectrode implantation. Main results . The study found that cortical layer 4, the input-receiving layer, maintained activity over the implantation time. However, layers 2/3 rapidly experienced severe impairment, leading to a loss of proper intralaminar connectivity in the downstream output layers 5/6. Furthermore, the impairment of interlaminar connectivity near the microelectrode was unidirectional, showing decreased connectivity from Layers 2/3 to Layers 5/6 but not the reverse direction. In the hippocampus, CA1 neurons gradually became unable to properly entrain to the surrounding LFP oscillations. Significance . This study provides a detailed characterization of network connectivity dysfunction over long-term microelectrode implantation periods. This new knowledge could contribute to the development of targeted therapeutic strategies aimed at improving the health of the tissue surrounding brain implants and potentially inform engineering of adaptive decoders as the FBR progresses. Our study’s understanding of the dynamic changes in the functional network over time opens the door to developing interventions for improving the long-term stability and performance of intracortical microelectrodes.

The feedback projections from cortical layer 6 (L6CT) to the sensory thalamus have long been implicated in playing a primary role in gating sensory signaling but remain poorly understood. To causally elucidate the full range of effects of these projections, we targeted silicon probe recordings to the whisker thalamocortical circuit of awake mice selectively expressing Channelrhodopsin-2 in L6CT neurons. Through optogenetic manipulation of L6CT neurons, multi-site electrophysiological recordings, and modeling of L6CT circuitry, we establish L6CT neurons as dynamic modulators of ongoing spiking in the ventral posteromedial nucleus of the thalamus (VPm), either suppressing or enhancing VPm spiking depending on L6CT neurons’ firing rate and synchrony. Differential effects across the cortical excitatory and inhibitory sub-populations point to an overall influence of L6CT feedback on cortical excitability that could have profound implications for regulating sensory signaling across a range of ethologically relevant conditions.

Sleep spindles appear to play an important role in learning new motor skills. Motor skill learning engages several brain regions with two important areas being the motor cortex (M1) and the cerebellum (CB). However, the neurophysiological processes in these areas during sleep, especially how spindle oscillations affect local and cross-region spiking, are not fully understood. We recorded an activity from the M1 and cerebellar cortex in eight rats during spontaneous activity to investigate how sleep spindles in these regions are related to local spiking as well as cross-region spiking. We found that M1 firing was significantly changed during both M1 and CB spindles, and this spiking occurred at a preferred phase of the spindle. On average, M1 and CB neurons showed most spiking at the M1 or CB spindle peaks. These neurons also developed a preferential phase locking to local or cross-area spindles with the greatest phase-locking value at spindle peaks; however, this preferential phase locking was not significant for cerebellar neurons when compared with CB spindles. Additionally, we found that the percentage of task-modulated cells in the M1 and CB that fired with nonuniform spike phase distribution during M1/CB spindle peaks were greater in the rats that learned a reach-to-grasp motor task robustly. Finally, we found that spindle band LFP coherence (for M1 and CB LFPs) showed a positive correlation with success rate in the motor task. These findings support the idea that sleep spindles in both the M1 and CB recruit neurons that participate in the awake task to support motor memory consolidation.

Temporally coordinated neural activity is central to nervous system function and purposeful behavior. Still, there is a paucity of evidence demonstrating how this coordinated activity within cortical and subcortical regions governs behavior. We investigated this between the primary motor (M1) and contralateral cerebellar cortex as rats learned a neuroprosthetic/brain-machine interface (BMI) task. In neuroprosthetic task, actuator movements are causally linked to M1 “direct” neurons that drive the decoder for successful task execution. However, it is unknown how task-related M1 activity interacts with the cerebellum. We observed a notable 3 to 6 hertz coherence that emerged between these regions’ local field potentials (LFPs) with learning that also modulated task-related spiking. We identified robust task-related indirect modulation in the cerebellum, which developed a preferential relationship with M1 task–related activity. Inhibiting cerebellar cortical and deep nuclei activity through optogenetics led to performance impairments in M1-driven neuroprosthetic control. Together, these results demonstrate that cerebellar influence is necessary for M1-driven neuroprosthetic control.

Hand movements are associated with modulations of neuronal activity across several interconnected cortical areas, including the primary motor cortex (M1) and the dorsal and ventral premotor cortices (PMd and PMv). Local field potentials (LFPs) provide a link between neuronal discharges and synaptic inputs. Our current understanding of how LFPs vary in M1, PMd, and PMv during contralateral and ipsilateral movements is incomplete. To help reveal unique features in the pattern of modulations, we simultaneously recorded LFPs in these areas in two macaque monkeys performing reach and grasp movements with either the right or left hand. The greatest effector-dependent differences were seen in M1, at low (≤13 Hz) and γ frequencies. In premotor areas, differences related to hand use were only present in low frequencies. PMv exhibited the greatest increase in low frequencies during instruction cues and the smallest effector-dependent modulation during movement execution. In PMd, δ oscillations were greater during contralateral reach and grasp, and β activity increased during contralateral grasp. In contrast, β oscillations decreased in M1 and PMv. These results suggest that while M1 primarily exhibits effector-specific LFP activity, premotor areas compute more effector-independent aspects of the task requirements, particularly during movement preparation for PMv and production for PMd. The generation of precise hand movements likely relies on the combination of complementary information contained in the unique pattern of neural modulations contained in each cortical area. Accordingly, integrating LFPs from premotor areas and M1 could enhance the performance and robustness of brain-machine interfaces.

Current soft neural probes are still operated by bulky, rigid electronics mounted to a body, which deteriorate the integrity of the device to biological systems and restrict the free behavior of a subject. We report a soft, conformable neural interface system that can monitor the single-unit activities of neurons with long-term stability. The system implements soft neural probes in the brain, and their subsidiary electronics which are directly printed on the cranial surface. The high-resolution printing of liquid metals forms soft neural probes with a cellular-scale diameter and adaptable lengths. Also, the printing of liquid metal-based circuits and interconnections along the curvature of the cranium enables the conformal integration of electronics to the body, and the cranial circuit delivers neural signals to a smartphone wirelessly. In the in-vivo studies using mice, the system demonstrates long-term recording (33 weeks) of neural activities in arbitrary brain regions. In T-maze behavioral tests, the system shows the behavior-induced activation of neurons in multiple brain regions.

Melody is a core component of music in which discrete pitches are serially arranged to convey emotion and meaning. Perception varies along several pitch-based dimensions: (i) the absolute pitch of notes, (ii) the difference in pitch between successive notes, and (iii) the statistical expectation of each note given prior context. How the brain represents these dimensions and whether their encoding is specialized for music remains unknown. We recorded high-density neurophysiological activity directly from the human auditory cortex while participants listened to Western musical phrases. Pitch, pitch-change, and expectation were selectively encoded at different cortical sites, indicating a spatial map for representing distinct melodic dimensions. The same participants listened to spoken English, and we compared responses to music and speech. Cortical sites selective for music encoded expectation, while sites that encoded pitch and pitch-change in music used the same neural code to represent equivalent properties of speech. Findings reveal how the perception of melody recruits both music-specific and general-purpose sound representations.

Stimulus-specific adaptation (SSA), the reduction of neural activity to a common stimulus that does not generalize to other, rare stimuli, is an essential property of our brain. Although well characterized in adults, it is still unknown how it develops during adolescence and what neuronal circuits are involved. Using in vivo electrophysiology and optogenetics in the lemniscal pathway of the mouse auditory system, we observed SSA to be stable from postnatal day 20 (P20) in the inferior colliculus, to develop until P30 in the auditory thalamus and even later in the primary auditory cortex (A1). We found this maturation process to be experience-dependent in A1 but not in thalamus and to be related to alterations in deep but not input layers of A1. We also identified corticothalamic projections to be implicated in thalamic SSA development. Together, our results reveal different circuits underlying the sequential SSA maturation and provide a unique perspective to understand predictive coding and surprise across sensory systems.