Columnar Localization and Laminar Origin of Cortical Surface Electrical Potentials

Recording of single-unit activities with flexible micro-electrocorticographic array in rats for decoding of whole-body navigation

Objective.Micro-electrocorticographic (μECoG) arrays are able to record neural activities from the cortical surface, without the need to penetrate the brain parenchyma. Owing in part to small electrode sizes, previous studies have demonstrated that single-unit spikes could be detected from the cortical surface, and likely from Layer I neurons of the neocortex. Here we tested the ability to useμECoG arrays to decode, in rats, body position during open field navigation, through isolated single-unit activities.Approach. μECoG arrays were chronically implanted onto primary motor cortex (M1) of Wistar rats, and neural recording was performed in awake, behaving rats in an open-field enclosure. The signals were band-pass filtered between 300-3000 Hz. Threshold-crossing spikes were identified and sorted into distinct units based on defined criteria including waveform morphology and refractory period. Body positions were derived from video recordings. We used gradient-boosting machine to predict body position based on previous 100 ms of spike data, and correlation analyses to elucidate the relationship between position and spike patterns.Main results.Single-unit spikes could be extracted during chronic recording fromμECoG, and spatial position could be decoded from these spikes with a mean absolute error of prediction of 0.135 and 0.090 in the x- and y- dimensions (of a normalized range from 0 to 1), and Pearson's r of 0.607 and 0.571, respectively.Significance. μECoG can detect single-unit activities that likely arise from superficial neurons in the cortex and is a promising alternative to intracortical arrays, with the added benefit of scalability to cover large cortical surface with minimal incremental risks. More studies should be performed in human related to its use as brain-machine interface.

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets.

Large-scale single-neuron speech sound encoding across the depth of human cortex

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.

Resolving Non-identifiability Mitigates Bias in Models of Neural Tuning and Functional Coupling

In the brain, all neurons are driven by the activity of other neurons, some of which maybe simultaneously recorded, but most are not. As such, models of neuronal activity need to account for simultaneously recorded neurons and the influences of unmeasured neurons. This can be done through inclusion of model terms for observed external variables (e.g., tuning to stimuli) as well as terms for latent sources of variability. Determining the influence of groups of neurons on each other relative to other influences is important to understand brain functioning. The parameters of statistical models fit to data are commonly used to gain insight into the relative importance of those influences. Scientific interpretation of models hinge upon unbiased parameter estimates. However, evaluation of biased inference is rarely performed and sources of bias are poorly understood. Through extensive numerical study and analytic calculation, we show that common inference procedures and models are typically biased. We demonstrate that accurate parameter selection before estimation resolves model non-identifiability and mitigates bias. In diverse neurophysiology data sets, we found that contributions of coupling to other neurons are often overestimated while tuning to exogenous variables are underestimated in common methods. We explain heterogeneity in observed biases across data sets in terms of data statistics. Finally, counter to common intuition, we found that model non-identifiability contributes to bias, not variance, making it a particularly insidious form of statistical error. Together, our results identify the causes of statistical biases in common models of neural data, provide inference procedures to mitigate that bias, and reveal and explain the impact of those biases in diverse neural data sets.