Thalamic epileptic spikes disrupt sleep spindles in patients with epileptic encephalopathy.

In severe epileptic encephalopathies, epileptic activity contributes to progressive cognitive dysfunction. Epileptic encephalopathies share the trait of spike-wave activation during non-REM sleep (EE-SWAS), a sleep stage dominated by sleep spindles, which are brain oscillations known to coordinate offline memory consolidation. Epileptic activity has been proposed to hijack the circuits driving these thalamocortical oscillations, thereby contributing to cognitive impairment. Using a unique dataset of simultaneous human thalamic and cortical recordings in subjects with and without EE-SWAS, we provide evidence for epileptic spike interference of thalamic sleep spindle production in patients with EE-SWAS. First, we show that epileptic spikes and sleep spindles are both predicted by slow oscillations during stage two sleep (N2), but at different phases of the slow oscillation. Next, we demonstrate that sleep-activated cortical epileptic spikes propagate to the thalamus (thalamic spike rate increases after a cortical spike, P ≈ 0). We then show that epileptic spikes in the thalamus increase the thalamic spindle refractory period (P ≈ 0). Finally, we show that in three patients with EE-SWAS, there is a downregulation of sleep spindles for 30 s after each thalamic spike (P < 0.01). These direct human thalamocortical observations support a proposed mechanism for epileptiform activity to impact cognitive function, wherein epileptic spikes inhibit thalamic sleep spindles in epileptic encephalopathy with spike and wave activation during sleep.

Spike ripples localize the epileptogenic zone best: an international intracranial study.

We evaluated whether spike ripples, the combination of epileptiform spikes and ripples, provide a reliable and improved biomarker for the epileptogenic zone compared with other leading interictal biomarkers in a multicentre, international study. We first validated an automated spike ripple detector on intracranial EEG recordings. We then applied this detector to subjects from four centres who subsequently underwent surgical resection with known 1-year outcomes. We evaluated the spike ripple rate in subjects cured after resection [International League Against Epilepsy Class 1 outcome (ILAE 1)] and those with persistent seizures (ILAE 2-6) across sites and recording types. We also evaluated available interictal biomarkers: spike, spike-gamma, wideband high frequency oscillation (HFO, 80-500 Hz), ripple (80-250 Hz) and fast ripple (250-500 Hz) rates using previously validated automated detectors. The proportion of resected events was computed and compared across subject outcomes and biomarkers. Overall, 109 subjects were included. Most spike ripples were removed in subjects with ILAE 1 outcome (P < 0.001), and this was qualitatively observed across all sites and for depth and subdural electrodes (P < 0.001 and P < 0.001, respectively). Among ILAE 1 subjects, the mean spike ripple rate was higher in the resected volume (0.66/min) than in the non-removed tissue (0.08/min, P < 0.001). A higher proportion of spike ripples were removed in subjects with ILAE 1 outcomes compared with ILAE 2-6 outcomes (P = 0.06). Among ILAE 1 subjects, the proportion of spike ripples removed was higher than the proportion of spikes (P < 0.001), spike-gamma (P < 0.001), wideband HFOs (P < 0.001), ripples (P = 0.009) and fast ripples (P = 0.009) removed. At the individual level, more subjects with ILAE 1 outcomes had the majority of spike ripples removed (79%, 38/48) than spikes (69%, P = 0.12), spike-gamma (69%, P = 0.12), wideband HFOs (63%, P = 0.03), ripples (45%, P = 0.01) or fast ripples (36%, P < 0.001) removed. Thus, in this large, multicentre cohort, when surgical resection was successful, the majority of spike ripples were removed. Furthermore, automatically detected spike ripples localize the epileptogenic tissue better than spikes, spike-gamma, wideband HFOs, ripples and fast ripples.

Multiple sources of fast traveling waves during human seizures: resolving a controversy.

During human seizures organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N=11 patients, 2 females, N=31 seizures) from previous human studies to analyze the traveling wave source. We find - inconsistent with both existing theories - a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo We conclude that combining both existing theories can generate the diversity of ictal traveling waves.Significance StatementThe source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Focal Sleep Spindle Deficits Reveal Focal Thalamocortical Dysfunction and Predict Cognitive Deficits in Sleep Activated Developmental Epilepsy.

Childhood epilepsy with centrotemporal spikes (CECTS) is the most common focal epilepsy syndrome, yet the cause of this disease remains unknown. Now recognized as a mild epileptic encephalopathy, children exhibit sleep-activated focal epileptiform discharges and cognitive difficulties during the active phase of the disease. The association between the abnormal electrophysiology and sleep suggests disruption to thalamocortical circuits. Thalamocortical circuit dysfunction resulting in pathologic epileptiform activity could hinder the production of sleep spindles, a brain rhythm essential for memory processes. Despite this pathophysiologic connection, the relationship between spindles and cognitive symptoms in epileptic encephalopathies has not been previously evaluated. A significant challenge limiting such work has been the poor performance of available automated spindle detection methods in the setting of sharp activities, such as epileptic spikes. Here, we validate a robust new method to accurately measure sleep spindles in patients with epilepsy. We then apply this detector to a prospective cohort of male and female children with CECTS with combined high-density EEGs during sleep and cognitive testing at varying time points of disease. We show that: (1) children have a transient, focal deficit in spindles during the symptomatic phase of disease; (2) spindle rate anticorrelates with spike rate; and (3) spindle rate, but not spike rate, predicts performance on cognitive tasks. These findings demonstrate focal thalamocortical circuit dysfunction and provide a pathophysiological explanation for the shared seizures and cognitive symptoms in CECTS. Further, this work identifies sleep spindles as a potential treatment target of cognitive dysfunction in this common epileptic encephalopathy.SIGNIFICANCE STATEMENT Childhood epilepsy with centrotemporal spikes is the most common idiopathic focal epilepsy syndrome, characterized by self-limited focal seizures and cognitive symptoms. Here, we provide the first evidence that focal thalamocortical circuit dysfunction underlies the shared seizures and cognitive dysfunction observed. In doing so, we identify sleep spindles as a mechanistic biomarker, and potential treatment target, of cognitive dysfunction in this common developmental epilepsy and provide a novel method to reliably quantify spindles in brain recordings from patients with epilepsy.

Direct Discriminative Decoder Models for Analysis of High-Dimensional Dynamical Neural Data.

With the accelerated development of neural recording technology over the past few decades, research in integrative neuroscience has become increasingly reliant on data analysis methods that are scalable to high-dimensional recordings and computationally tractable. Latent process models have shown promising results in estimating the dynamics of cognitive processes using individual models for each neuron's receptive field. However, scaling these models to work on high-dimensional neural recordings remains challenging. Not only is it impractical to build receptive field models for individual neurons of a large neural population, but most neural data analyses based on individual receptive field models discard the local history of neural activity, which has been shown to be critical in the accurate inference of the underlying cognitive processes. Here, we propose a novel, scalable latent process model that can directly estimate cognitive process dynamics without requiring precise receptive field models of individual neurons or brain nodes. We call this the direct discriminative decoder (DDD) model. The DDD model consists of (1) a discriminative process that characterizes the conditional distribution of the signal to be estimated, or state, as a function of both the current neural activity and its local history, and (2) a state transition model that characterizes the evolution of the state over a longer time period. While this modeling framework inherits advantages of existing latent process modeling methods, its computational cost is tractable. More important, the solution can incorporate any information from the history of neural activity at any timescale in computing the estimate of the state process. There are many choices in building the discriminative process, including deep neural networks or gaussian processes, which adds to the flexibility of the framework. We argue that these attributes of the proposed methodology, along with its applicability to different modalities of neural data, make it a powerful tool for high-dimensional neural data analysis. We also introduce an extension of these methods, called the discriminative-generative decoder (DGD). The DGD includes both discriminative and generative processes in characterizing observed data. As a result, we can combine physiological correlates like behavior with neural data to better estimate underlying cognitive processes. We illustrate the methods, including steps for inference and model identification, and demonstrate applications to multiple data analysis problems with high-dimensional neural recordings. The modeling results demonstrate the computational and modeling advantages of the DDD and DGD methods.

Risk-taking bias in human decision-making is encoded via a right-left brain push-pull system.

A person's decisions vary even when options stay the same, like when a gambler changes bets despite constant odds of winning. Internal bias (e.g., emotion) contributes to this variability and is shaped by past outcomes, yet its neurobiology during decision-making is not well understood. To map neural circuits encoding bias, we administered a gambling task to 10 participants implanted with intracerebral depth electrodes in cortical and subcortical structures. We predicted the variability in betting behavior within and across patients by individual bias, which is estimated through a dynamical model of choice. Our analysis further revealed that high-frequency activity increased in the right hemisphere when participants were biased toward risky bets, while it increased in the left hemisphere when participants were biased away from risky bets. Our findings provide electrophysiological evidence that risk-taking bias is a lateralized push-pull neural system governing counterintuitive and highly variable decision-making in humans.

Assessing Goodness-of-Fit in Marked Point Process Models of Neural Population Coding via Time and Rate Rescaling.

Marked point process models have recently been used to capture the coding properties of neural populations from multiunit electrophysiological recordings without spike sorting. These clusterless models have been shown in some instances to better describe the firing properties of neural populations than collections of receptive field models for sorted neurons and to lead to better decoding results. To assess their quality, we previously proposed a goodness-of-fit technique for marked point process models based on time rescaling, which for a correct model produces a set of uniform samples over a random region of space. However, assessing uniformity over such a region can be challenging, especially in high dimensions. Here, we propose a set of new transformations in both time and the space of spike waveform features, which generate events that are uniformly distributed in the new mark and time spaces. These transformations are scalable to multidimensional mark spaces and provide uniformly distributed samples in hypercubes, which are well suited for uniformity tests. We discuss the properties of these transformations and demonstrate aspects of model fit captured by each transformation. We also compare multiple uniformity tests to determine their power to identify lack-of-fit in the rescaled data. We demonstrate an application of these transformations and uniformity tests in a simulation study. Proofs for each transformation are provided in the appendix.

Scalp recorded spike ripples predict seizure risk in childhood epilepsy better than spikes.

In the past decade, brief bursts of fast oscillations in the ripple range have been identified in the scalp EEG as a promising non-invasive biomarker for epilepsy. However, investigation and clinical application of this biomarker have been limited because standard approaches to identify these brief, low amplitude events are difficult, time consuming, and subjective. Recent studies have demonstrated that ripples co-occurring with epileptiform discharges ('spike ripple events') are easier to detect than ripples alone and have greater pathological significance. Here, we used objective techniques to quantify spike ripples and test whether this biomarker predicts seizure risk in childhood epilepsy. We evaluated spike ripples in scalp EEG recordings from a prospective cohort of children with a self-limited epilepsy syndrome, benign epilepsy with centrotemporal spikes, and healthy control children. We compared the rate of spike ripples between children with epilepsy and healthy controls, and between children with epilepsy during periods of active disease (active, within 1 year of seizure) and after a period of sustained seizure-freedom (seizure-free, >1 year without seizure), using semi-automated and automated detection techniques. Spike ripple rate was higher in subjects with active epilepsy compared to healthy controls (P = 0.0018) or subjects with epilepsy who were seizure-free ON or OFF medication (P = 0.0018). Among epilepsy subjects with spike ripples, each month seizure-free decreased the odds of a spike ripple by a factor of 0.66 [95% confidence interval (0.47, 0.91), P = 0.021]. Comparing the diagnostic accuracy of the presence of at least one spike ripple versus a classic spike event to identify group, we found comparable sensitivity and negative predictive value, but greater specificity and positive predictive value of spike ripples compared to spikes (P = 0.016 and P = 0.006, respectively). We found qualitatively consistent results using a fully automated spike ripple detector, including comparison with an automated spike detector. We conclude that scalp spike ripple events identify disease and track with seizure risk in this epilepsy population, using both semi-automated and fully automated detection methods, and that this biomarker outperforms analysis of spikes alone in categorizing seizure risk. These data provide evidence that spike ripples are a specific non-invasive biomarker for seizure risk in benign epilepsy with centrotemporal spikes and support future work to evaluate the utility of this biomarker to guide medication trials and tapers in these children and predict seizure risk in other at-risk populations.

Dynamic connectivity modulates local activity in the core regions of the default-mode network.

Segregation and integration are distinctive features of large-scale brain activity. Although neuroimaging studies have been unraveling their neural correlates, how integration takes place over segregated modules remains elusive. Central to this problem is the mechanism by which a brain region adjusts its activity according to the influence it receives from other regions. In this study, we explore how dynamic connectivity between two regions affects the neural activity within a participating region. Combining functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) in the same group of subjects, we analyzed resting-state data from the core of the default-mode network. We observed directed influence from the posterior cingulate cortex (PCC) to the anterior cingulate cortex (ACC) in the 10-Hz range. This time-varying influence was associated with the power alteration in the ACC: strong influence corresponded with a decrease of power around 13-16 Hz and an increase of power in the lower (1-7 Hz) and higher (30-55 Hz) ends of the spectrum. We also found that the amplitude of the 30- to 55-Hz activity was coupled to the phase of the 3- to 4-Hz activity in the ACC. These results characterized the local spectral changes associated with network interactions. The specific spectral information both highlights the functional roles of PCC-ACC connectivity in the resting state and provides insights into the dynamic relationship between local activity and coupling dynamics of a network.

Decoding Hidden Cognitive States From Behavior and Physiology Using a Bayesian Approach.

Cognitive processes, such as learning and cognitive flexibility, are both difficult to measure and to sample continuously using objective tools because cognitive processes arise from distributed, high-dimensional neural activity. For both research and clinical applications, that dimensionality must be reduced. To reduce dimensionality and measure underlying cognitive processes, we propose a modeling framework in which a cognitive process is defined as a low-dimensional dynamical latent variable-called a cognitive state, which links high-dimensional neural recordings and multidimensional behavioral readouts. This framework allows us to decompose the hard problem of modeling the relationship between neural and behavioral data into separable encoding-decoding approaches. We first use a state-space modeling framework, the behavioral decoder, to articulate the relationship between an objective behavioral readout (e.g., response times) and cognitive state. The second step, the neural encoder, involves using a generalized linear model (GLM) to identify the relationship between the cognitive state and neural signals, such as local field potential (LFP). We then use the neural encoder model and a Bayesian filter to estimate cognitive state using neural data (LFP power) to generate the neural decoder. We provide goodness-of-fit analysis and model selection criteria in support of the encoding-decoding result. We apply this framework to estimate an underlying cognitive state from neural data in human participants (N=8) performing a cognitive conflict task. We successfully estimated the cognitive state within the 95% confidence intervals of that estimated using behavior readout for an average of 90% of task trials across participants. In contrast to previous encoder-decoder models, our proposed modeling framework incorporates LFP spectral power to encode and decode a cognitive state. The framework allowed us to capture the temporal evolution of the underlying cognitive processes, which could be key to the development of closed-loop experiments and treatments.

A common goodness-of-fit framework for neural population models using marked point process time-rescaling.

A critical component of any statistical modeling procedure is the ability to assess the goodness-of-fit between a model and observed data. For spike train models of individual neurons, many goodness-of-fit measures rely on the time-rescaling theorem and assess model quality using rescaled spike times. Recently, there has been increasing interest in statistical models that describe the simultaneous spiking activity of neuron populations, either in a single brain region or across brain regions. Classically, such models have used spike sorted data to describe relationships between the identified neurons, but more recently clusterless modeling methods have been used to describe population activity using a single model. Here we develop a generalization of the time-rescaling theorem that enables comprehensive goodness-of-fit analysis for either of these classes of population models. We use the theory of marked point processes to model population spiking activity, and show that under the correct model, each spike can be rescaled individually to generate a uniformly distributed set of events in time and the space of spike marks. After rescaling, multiple well-established goodness-of-fit procedures and statistical tests are available. We demonstrate the application of these methods both to simulated data and real population spiking in rat hippocampus. We have made the MATLAB and Python code used for the analyses in this paper publicly available through our Github repository at https://github.com/Eden-Kramer-Lab/popTRT .

Capturing Spike Variability in Noisy Izhikevich Neurons Using Point Process Generalized Linear Models.

To understand neural activity, two broad categories of models exist: statistical and dynamical. While statistical models possess rigorous methods for parameter estimation and goodness-of-fit assessment, dynamical models provide mechanistic insight. In general, these two categories of models are separately applied; understanding the relationships between these modeling approaches remains an area of active research. In this letter, we examine this relationship using simulation. To do so, we first generate spike train data from a well-known dynamical model, the Izhikevich neuron, with a noisy input current. We then fit these spike train data with a statistical model (a generalized linear model, GLM, with multiplicative influences of past spiking). For different levels of noise, we show how the GLM captures both the deterministic features of the Izhikevich neuron and the variability driven by the noise. We conclude that the GLM captures essential features of the simulated spike trains, but for near-deterministic spike trains, goodness-of-fit analyses reveal that the model does not fit very well in a statistical sense; the essential random part of the GLM is not captured.

Measuring the signal-to-noise ratio of a neuron.

The signal-to-noise ratio (SNR), a commonly used measure of fidelity in physical systems, is defined as the ratio of the squared amplitude or variance of a signal relative to the variance of the noise. This definition is not appropriate for neural systems in which spiking activity is more accurately represented as point processes. We show that the SNR estimates a ratio of expected prediction errors and extend the standard definition to one appropriate for single neurons by representing neural spiking activity using point process generalized linear models (PP-GLM). We estimate the prediction errors using the residual deviances from the PP-GLM fits. Because the deviance is an approximate χ(2) random variable, we compute a bias-corrected SNR estimate appropriate for single-neuron analysis and use the bootstrap to assess its uncertainty. In the analyses of four systems neuroscience experiments, we show that the SNRs are -10 dB to -3 dB for guinea pig auditory cortex neurons, -18 dB to -7 dB for rat thalamic neurons, -28 dB to -14 dB for monkey hippocampal neurons, and -29 dB to -20 dB for human subthalamic neurons. The new SNR definition makes explicit in the measure commonly used for physical systems the often-quoted observation that single neurons have low SNRs. The neuron's spiking history is frequently a more informative covariate for predicting spiking propensity than the applied stimulus. Our new SNR definition extends to any GLM system in which the factors modulating the response can be expressed as separate components of a likelihood function.

Rapid classification of hippocampal replay content for real-time applications.

Sharp-wave ripple (SWR) events in the hippocampus replay millisecond-timescale patterns of place cell activity related to the past experience of an animal. Interrupting SWR events leads to learning and memory impairments, but how the specific patterns of place cell spiking seen during SWRs contribute to learning and memory remains unclear. A deeper understanding of this issue will require the ability to manipulate SWR events based on their content. Accurate real-time decoding of SWR replay events requires new algorithms that are able to estimate replay content and the associated uncertainty, along with software and hardware that can execute these algorithms for biological interventions on a millisecond timescale. Here we develop an efficient estimation algorithm to categorize the content of replay from multiunit spiking activity. Specifically, we apply real-time decoding methods to each SWR event and then compute the posterior probability of the replay feature. We illustrate this approach by classifying SWR events from data recorded in the hippocampus of a rat performing a spatial memory task into four categories: whether they represent outbound or inbound trajectories and whether the activity is replayed forward or backward in time. We show that our algorithm can classify the majority of SWR events in a recording epoch within 20 ms of the replay onset with high certainty, which makes the algorithm suitable for a real-time implementation with short latencies to incorporate into content-based feedback experiments.

Examination of rhythmicity of extracellularly recorded neurons in the entorhinal cortex.

A number of studies have examined the theta-rhythmic modulation of neuronal firing in the hippocampal circuit. For extracellular recordings, this is often done by examining spectral properties of the spike-time autocorrelogram, most significantly, for validating the presence or absence of theta modulation across species. These techniques can show significant rhythmicity for high firing rate, highly rhythmic neurons; however, they are substantially biased by several factors including the peak firing rate of the neuron, the amount of time spent in the neuron's receptive field, and other temporal properties of the rhythmicity such as cycle-skipping. These limitations make it difficult to examine rhythmic modulation in neurons with low firing rates or when an animal has short dwell times within the firing field and difficult to compare rhythmicity under disparate experimental conditions when these factors frequently differ. Here, we describe in detail the challenges that researchers face when using these techniques and apply our findings to recent recordings from bat entorhinal grid cells, suggesting that they may have lacked enough data to examine theta rhythmicity robustly. We describe a more sensitive and statistically rigorous method using maximum likelihood estimation (MLE) of a parametric model of the lags within the autocorrelation window, which helps to alleviate some of the problems of traditional methods and was also unable to detect rhythmicity in bat grid cells. Using large batteries of simulated data, we explored the boundaries for which the MLE technique and the theta index can detect rhythmicity. The MLE technique is less sensitive to many features of the autocorrelogram and provides a framework for statistical testing to detect rhythmicity as well as changes in rhythmicity in individual sessions providing a substantial improvement over previous methods.

Clusterless Decoding of Position from Multiunit Activity Using a Marked Point Process Filter.

Point process filters have been applied successfully to decode neural signals and track neural dynamics. Traditionally these methods assume that multiunit spiking activity has already been correctly spike-sorted. As a result, these methods are not appropriate for situations where sorting cannot be performed with high precision, such as real-time decoding for brain-computer interfaces. Because the unsupervised spike-sorting problem remains unsolved, we took an alternative approach that takes advantage of recent insights into clusterless decoding. Here we present a new point process decoding algorithm that does not require multiunit signals to be sorted into individual units. We use the theory of marked point processes to construct a function that characterizes the relationship between a covariate of interest (in this case, the location of a rat on a track) and features of the spike waveforms. In our example, we use tetrode recordings, and the marks represent a four-dimensional vector of the maximum amplitudes of the spike waveform on each of the four electrodes. In general, the marks may represent any features of the spike waveform. We then use Bayes's rule to estimate spatial location from hippocampal neural activity. We validate our approach with a simulation study and experimental data recorded in the hippocampus of a rat moving through a linear environment. Our decoding algorithm accurately reconstructs the rat's position from unsorted multiunit spiking activity. We then compare the quality of our decoding algorithm to that of a traditional spike-sorting and decoding algorithm. Our analyses show that the proposed decoding algorithm performs equivalent to or better than algorithms based on sorted single-unit activity. These results provide a path toward accurate real-time decoding of spiking patterns that could be used to carry out content-specific manipulations of population activity in hippocampus or elsewhere in the brain.

Human seizures self-terminate across spatial scales via a critical transition.

Why seizures spontaneously terminate remains an unanswered fundamental question of epileptology. Here we present evidence that seizures self-terminate via a discontinuous critical transition or bifurcation. We show that human brain electrical activity at various spatial scales exhibits common dynamical signatures of an impending critical transition--slowing, increased correlation, and flickering--in the approach to seizure termination. In contrast, prolonged seizures (status epilepticus) repeatedly approach, but do not cross, the critical transition. To support these results, we implement a computational model that demonstrates that alternative stable attractors, representing the ictal and postictal states, emulate the observed dynamics. These results suggest that self-terminating seizures end through a common dynamical mechanism. This description constrains the specific biophysical mechanisms underlying seizure termination, suggests a dynamical understanding of status epilepticus, and demonstrates an accessible system for studying critical transitions in nature.

Impact of brain tissue filtering on neurostimulation fields: a modeling study.

Electrical neurostimulation techniques, such as deep brain stimulation (DBS) and transcranial magnetic stimulation (TMS), are increasingly used in the neurosciences, e.g., for studying brain function, and for neurotherapeutics, e.g., for treating depression, epilepsy, and Parkinson's disease. The characterization of electrical properties of brain tissue has guided our fundamental understanding and application of these methods, from electrophysiologic theory to clinical dosing-metrics. Nonetheless, prior computational models have primarily relied on ex-vivo impedance measurements. We recorded the in-vivo impedances of brain tissues during neurosurgical procedures and used these results to construct MRI guided computational models of TMS and DBS neurostimulatory fields and conductance-based models of neurons exposed to stimulation. We demonstrated that tissues carry neurostimulation currents through frequency dependent resistive and capacitive properties not typically accounted for by past neurostimulation modeling work. We show that these fundamental brain tissue properties can have significant effects on the neurostimulatory-fields (capacitive and resistive current composition and spatial/temporal dynamics) and neural responses (stimulation threshold, ionic currents, and membrane dynamics). These findings highlight the importance of tissue impedance properties on neurostimulation and impact our understanding of the biological mechanisms and technological potential of neurostimulatory methods.

Characterizing context-dependent differential firing activity in the hippocampus and entorhinal cortex.

The rat hippocampus and entorhinal cortex have been shown to possess neurons with place fields that modulate their firing properties under different behavioral contexts. Such context-dependent changes in neural activity are commonly studied through electrophysiological experiments in which a rat performs a continuous spatial alternation task on a T-maze. Previous research has analyzed context-based differential firing during this task by describing differences in the mean firing activity between left-turn and right-turn experimental trials. In this article, we develop qualitative and quantitative methods to characterize and compare changes in trial-to-trial firing rate variability for sets of experimental contexts. We apply these methods to cells in the CA1 region of hippocampus and in the dorsocaudal medial entorhinal cortex (dcMEC), characterizing the context-dependent differences in spiking activity during spatial alternation. We identify a subset of cells with context-dependent changes in firing rate variability. Additionally, we show that dcMEC populations encode turn direction uniformly throughout the T-maze stem, whereas CA1 populations encode context at major waypoints in the spatial trajectory. Our results suggest scenarios in which individual cells that sparsely provide information on turn direction might combine in the aggregate to produce a robust population encoding.