Objective.Intracranial neural recordings and electrical stimulation are tools used in an increasing range of applications, including intraoperative clinical mapping and monitoring, therapeutic neuromodulation, and brain computer interface control and feedback. However, many of these applications suffer from a lack of spatial specificity and localization, both in terms of sensed neural signal and applied stimulation. This stems from limited manufacturing processes of commercial-off-the-shelf (COTS) arrays unable to accommodate increased channel density, higher channel count, and smaller contact size.Approach.Here, we describe a manufacturing and assembly approach using thin-film microfabrication for 32-channel high density subdural micro-electrocorticography (µECoG) surface arrays (contacts 1.2 mm diameter, 2 mm pitch) and intracranial electroencephalography (iEEG) depth arrays (contacts 0.5 mm × 1.5 mm, pitch 0.8 mm × 2.5 mm). Crucially, we tackle the translational hurdle and test these arrays during intraoperative studies conducted in four humans under regulatory approval.Main results.We demonstrate that the higher-density contacts provide additional unique information across the recording span compared to the density of COTS arrays which typically have electrode pitch of 8 mm or greater; 4 mm in case of specially ordered arrays. Our intracranial stimulation study results reveal that refined spatial targeting of stimulation elicits evoked potentials with differing spatial spread.Significance.Thin-film,μECoG and iEEG depth arrays offer a promising substrate for advancing a number of clinical and research applications reliant on high-resolution neural sensing and intracranial stimulation.

Early studies of deep brain stimulation (DBS) for various neurological disorders involved a temporary trial period where implanted electrodes were externalized, in which the electrical contacts exiting the patient's brain are connected to external stimulation equipment, so that stimulation efficacy could be determined before permanent implant. As the optimal brain target sites for various diseases (i.e., Parkinson's disease, essential tremor) became better established, such trial periods have fallen out of favor. However, deep brain stimulation trial periods are experiencing a modern resurgence for at least two reasons: (1) studies of newer indications such as depression or chronic pain aim to identify new targets and (2) a growing interest in adaptive DBS tools necessitates neurophysiological recordings, which are often done in the peri-surgical period. In this review, we consider the possible approaches, benefits, and risks of such inpatient trial periods with a specific focus on developing new DBS therapies for chronic pain.

Higher-order visual thalamus plays a fundamental but poorly understood role in attention-demanding tasks. To investigate how neuronal dynamics in higher-order visual thalamus are modulated by sustained attention, we performed multichannel electrophysiological recordings in the lateral posterior-pulvinar complex (LP/pulvinar) in the ferret (Mustela putorius furo). We recorded single unit activity and local field potential (LFP) during the performance of the five-choice serial reaction time task (5-CSRTT), which is used in both humans and animals as an assay of sustained attention. We found that half of the units exhibited an increasing firing rate during the delay period before stimulus onset (attention-modulated units). In contrast, the non-attention-modulated units responded to the stimulus, but not during the delay period. Spike-field coherence (SFC) of only the attention-modulated neurons significantly increased from the start of the delay period until screen touch, predominantly in the θ frequency band. In addition, θ power and θ/γ phase amplitude coupling (PAC) were elevated throughout the delay period. Our findings suggest that the θ oscillation plays a central role in orchestrating thalamic signaling during sustained attention.

The disruption of coupling between brain areas has been suggested as the mechanism underlying loss of consciousness in anesthesia. This hypothesis has been tested previously by measuring the information transfer between brain areas, and by taking reduced information transfer as a proxy for decoupling. Yet, information transfer is a function of the amount of information available in the information source-such that transfer decreases even for unchanged coupling when less source information is available. Therefore, we reconsidered past interpretations of reduced information transfer as a sign of decoupling, and asked whether impaired local information processing leads to a loss of information transfer. An important prediction of this alternative hypothesis is that changes in locally available information (signal entropy) should be at least as pronounced as changes in information transfer. We tested this prediction by recording local field potentials in two ferrets after administration of isoflurane in concentrations of 0.0%, 0.5%, and 1.0%. We found strong decreases in the source entropy under isoflurane in area V1 and the prefrontal cortex (PFC)-as predicted by our alternative hypothesis. The decrease in source entropy was stronger in PFC compared to V1. Information transfer between V1 and PFC was reduced bidirectionally, but with a stronger decrease from PFC to V1. This links the stronger decrease in information transfer to the stronger decrease in source entropy-suggesting reduced source entropy reduces information transfer. This conclusion fits the observation that the synaptic targets of isoflurane are located in local cortical circuits rather than on the synapses formed by interareal axonal projections. Thus, changes in information transfer under isoflurane seem to be a consequence of changes in local processing more than of decoupling between brain areas. We suggest that source entropy changes must be considered whenever interpreting changes in information transfer as decoupling.

Understanding the neural basis of speech perception requires that we study the human brain both at the scale of the fundamental computational unit of neurons and in their organization across the depth of cortex. Here we used high-density Neuropixels arrays1-3 to record from 685 neurons across cortical layers at nine sites in a high-level auditory region that is critical for speech, the superior temporal gyrus4,5, while participants listened to spoken sentences. Single neurons encoded a wide range of speech sound cues, including features of consonants and vowels, relative vocal pitch, onsets, amplitude envelope and sequence statistics. Neurons at each cross-laminar recording exhibited dominant tuning to a primary speech feature while also containing a substantial proportion of neurons that encoded other features contributing to heterogeneous selectivity. Spatially, neurons at similar cortical depths tended to encode similar speech features. Activity across all cortical layers was predictive of high-frequency field potentials (electrocorticography), providing a neuronal origin for macroelectrode recordings from the cortical surface. Together, these results establish single-neuron tuning across the cortical laminae as an important dimension of speech encoding in human superior temporal gyrus.

Objectives

Adult patients with epilepsy have an increased prevalence of major depressive disorder (MDD). Intracranial EEG (iEEG) captured during extended inpatient monitoring of patients with treatment-resistant epilepsy offers a particularly promising method to study MDD networks in epilepsy.

Methods

The authors used 24 hours of resting-state iEEG to examine the neural activity patterns within corticolimbic structures that reflected the presence of depressive symptoms in 13 adults with medication-refractory epilepsy. Principal component analysis was performed on the z-scored mean relative power in five standard frequency bands averaged across electrodes within a region.

Results

Principal component 3 was a statistically significant predictor of the presence of depressive symptoms (R²=0.35, p=0.014). A balanced logistic classifier model using principal component 3 alone correctly classified 78% of patients as belonging to the group with a high burden of depressive symptoms or a control group with minimal depressive symptoms (sensitivity, 75%; specificity, 80%; area under the curve=0.8, leave-one-out cross validation). Classification was dependent on beta power throughout the corticolimbic network and low-frequency cingulate power.

Conclusions

These finding suggest, for the first time, that neural features across circuits involved in epilepsy may distinguish patients who have depressive symptoms from those who do not. Larger studies are required to validate these findings and to assess their diagnostic utility in MDD.

Sustained attention requires the coordination of neural activity across multiple cortical areas in the frontoparietal network, in particular the prefrontal cortex (PFC) and posterior parietal cortex (PPC). Previous work has demonstrated that activity in these brain regions is coordinated by neuronal oscillations of the local field potential (LFP). However, the underlying coordination of activity in terms of organization of single unit (SU) spiking activity has remained poorly understood, particularly in the freely moving animal. We found that long-range functional connectivity between anatomically connected PFC and PPC was mediated by oscillations in the theta frequency band. SU activity in PFC was phase locked to theta oscillations in PPC, and spiking activity in PFC and PPC was locked to local high-gamma activity. Together, our results support a model in which frequency-specific synchronization mediates functional connectivity between and within PFC and PPC of the frontoparietal attention network in the freely moving animal.

Although oscillations during development have been characterized in a wide range of neural systems, little is known about the interaction between these network oscillations and neuronal spiking, and the interactions among different oscillation frequencies. Here we recorded the spontaneous and visual-elicited local field potential (LFP) and multi-unit activity (MUA) in the visual cortex of freely-moving juvenile ferrets before and after eye-opening. We found that both the spontaneous and visually-elicited LFP power was increased after eye-opening, especially in higher frequency bands (>30 Hz). Spike LFP phase coupling was decreased for lower frequency bands (theta and alpha) but slightly increased for higher frequencies (high-gamma band). A similar shift towards faster frequencies also occurred for phase-amplitude coupling; with maturation, the coupling of the theta/alpha/beta band amplitude to the delta phase was decreased and the high-gamma amplitude coupling to theta/alpha phase was increased. This shift towards higher frequencies was also reflected in the visual responses; the LFP oscillation became more entrained by visual stimulation with higher frequencies (>10 Hz). Taken together, these results suggest gamma oscillation as a signature of the maturation of cortical circuitry.

The hippocampus is diversely interconnected with other brain systems along its axis. Cycles of theta-frequency activity are believed to propagate from the septal to temporal pole, yet it is unclear how this one-way route supports the flexible cognitive capacities of this structure. We leveraged novel thin-film microgrid arrays conformed to the human hippocampal surface to track neural activity two-dimensionally in vivo. All oscillation frequencies identified between 1-15 Hz propagated across the tissue. Moreover, they dynamically shifted between two roughly opposite directions oblique to the long axis. This predominant propagation axis was mirrored across participants, hemispheres, and consciousness states. Directionality was modulated in a participant who performed a behavioral task, and it could be predicted by wave amplitude topography over the hippocampal surface. Our results show that propagation directions may thus represent distinct meso-scale network computations, operating along versatile spatiotemporal processing routes across the hippocampal body.

Closed-loop neurostimulation is a promising therapy being tested and clinically implemented in a growing number of neurological and psychiatric indications. This therapy is enabled by chronically implanted, bidirectional devices including the Medtronic Summit RC+S system. In order to successfully optimize therapy for patients implanted with these devices, analyses must be conducted offline on the recorded neural data, in order to inform optimal sense and stimulation parameters. The file format, volume, and complexity of raw data from these devices necessitate conversion, parsing, and time reconstruction ahead of time-frequency analyses and modeling common to standard neuroscientific analyses. Here, we provide an open-source toolbox written in Matlab which takes raw files from the Summit RC+S and transforms these data into a standardized format amenable to conventional analyses. Furthermore, we provide a plotting tool which can aid in the visualization of multiple data streams and sense, stimulation, and therapy settings. Finally, we describe an analysis module which replicates RC+S on-board power computations, a functionality which can accelerate biomarker discovery. This toolbox aims to accelerate the research and clinical advances made possible by longitudinal neural recordings and adaptive neurostimulation in people with neurological and psychiatric illnesses.