Neural correlates of tactile detection: a combined magnetoencephalography and biophysically based computational modeling study

Neocortical inhibitory imbalance predicts successful sensory detection

Perceptual success depends on fast-spiking, parvalbumin-positive interneurons (FS/PVs). However, competing theories of optimal rate and correlation in pyramidal (PYR) firing make opposing predictions regarding the underlying FS/PV dynamics. We addressed this with population calcium imaging of FS/PVs and putative PYR neurons during threshold detection. In primary somatosensory and visual neocortex, a distinct PYR subset shows increased rate and spike-count correlations on detected trials ("hits"), while most show no rate change and decreased correlations. A larger fraction of FS/PVs predicts hits with either rate increases or decreases. Using computational modeling, we found that inhibitory imbalance, created by excitatory "feedback" and interactions between FS/PV pools, can account for the data. Rate-decreasing FS/PVs increase rate and correlation in a PYR subset, while rate-increasing FS/PVs reduce correlations and offset enhanced excitation in PYR neurons. These findings indicate that selection of informative PYR ensembles, through transient inhibitory imbalance, is a common motif of optimal neocortical processing.

Biophysical Modeling of Frontocentral ERP Generation Links Circuit-Level Mechanisms of Action-Stopping to a Behavioral Race Model

Human frontocentral event-related potentials (FC-ERPs) are ubiquitous neural correlates of cognition and control, but their generating multiscale mechanisms remain mostly unknown. We used the Human Neocortical Neurosolver's biophysical model of a canonical neocortical circuit under exogenous thalamic and cortical drive to simulate the cell and circuit mechanisms underpinning the P2, N2, and P3 features of the FC-ERP observed after Stop-Signals in the Stop-Signal task (SST; N = 234 humans, 137 female). We demonstrate that a sequence of simulated external thalamocortical and corticocortical drives can produce the FC-ERP, similar to what has been shown for primary sensory cortices. We used this model of the FC-ERP to examine likely circuit-mechanisms underlying FC-ERP features that distinguish between successful and failed action-stopping. We also tested their adherence to the predictions of the horse-race model of the SST, with specific hypotheses motivated by theoretical links between the P3 and Stop process. These simulations revealed that a difference in P3 onset between successful and failed Stops is most likely due to a later arrival of thalamocortical drive in failed Stops, rather than, for example, a difference in the effective strength of the input. In contrast, the same model predicted that early thalamocortical drives underpinning the P2 and N2 differed in both strength and timing across stopping accuracy conditions. Overall, this model generates novel testable predictions of the thalamocortical dynamics underlying FC-ERP generation during action-stopping. Moreover, it provides a detailed cellular and circuit-level interpretation that supports links between these macroscale signatures and predictions of the behavioral race model.

Neuronal Modeling of Cross-Sensory Visual Evoked Magnetoencephalography Responses in the Auditory Cortex

Previous studies have demonstrated that auditory cortex activity can be influenced by cross-sensory visual inputs. Intracortical laminar recordings in nonhuman primates have suggested a feedforward (FF) type profile for auditory evoked but feedback (FB) type for visual evoked activity in the auditory cortex. To test whether cross-sensory visual evoked activity in the auditory cortex is associated with FB inputs also in humans, we analyzed magnetoencephalography (MEG) responses from eight human subjects (six females) evoked by simple auditory or visual stimuli. In the estimated MEG source waveforms for auditory cortex regions of interest, auditory evoked response showed peaks at 37 and 90 ms and visual evoked response at 125 ms. The inputs to the auditory cortex were modeled through FF- and FB-type connections targeting different cortical layers using the Human Neocortical Neurosolver (HNN), which links cellular- and circuit-level mechanisms to MEG signals. HNN modeling suggested that the experimentally observed auditory response could be explained by an FF input followed by an FB input, whereas the cross-sensory visual response could be adequately explained by just an FB input. Thus, the combined MEG and HNN results support the hypothesis that cross-sensory visual input in the auditory cortex is of FB type. The results also illustrate how the dynamic patterns of the estimated MEG source activity can provide information about the characteristics of the input into a cortical area in terms of the hierarchical organization among areas.

Agranular frontal cortical microcircuit underlying cognitive control in macaques

The error-related negativity and an N2-component recorded over medial frontal cortex index core functions of cognitive control. While they are known to originate from agranular frontal areas, the underlying microcircuit mechanisms remain elusive. Most insights about microcircuit function have been derived from variations of the so-called canonical microcircuit model. These microcircuit architectures are based extensively on studies from granular sensory cortical areas in monkeys, cats, and rodents. However, evidence has shown striking cytoarchitectonic differences across species and differences in the functional relationships across cortical layers in agranular compared to granular sensory areas. In this minireview, we outline a tentative microcircuit model underlying cognitive control in the agranular frontal cortex of primates. The model incorporates the main GABAergic interneuron subclasses with specific laminar arrangements and target regions on pyramidal cells. We emphasize the role of layer 5 pyramidal cells in error and conflict detection. We offer several specific questions necessary for creating a specific intrinsic microcircuit model of the agranular frontal cortex.

Isolating Neural Signatures of Conscious Speech Perception with a No-Report Sine-Wave Speech Paradigm

Identifying neural correlates of conscious perception is a fundamental endeavor of cognitive neuroscience. Most studies so far have focused on visual awareness along with trial-by-trial reports of task-relevant stimuli, which can confound neural measures of perceptual awareness with postperceptual processing. Here, we used a three-phase sine-wave speech paradigm that dissociated between conscious speech perception and task relevance while recording EEG in humans of both sexes. Compared with tokens perceived as noise, physically identical sine-wave speech tokens that were perceived as speech elicited a left-lateralized, near-vertex negativity, which we interpret as a phonological version of a perceptual awareness negativity. This response appeared between 200 and 300 ms after token onset and was not present for frequency-flipped control tokens that were never perceived as speech. In contrast, the P3b elicited by task-irrelevant tokens did not significantly differ when the tokens were perceived as speech versus noise and was only enhanced for tokens that were both perceived as speech and relevant to the task. Our results extend the findings from previous studies on visual awareness and speech perception and suggest that correlates of conscious perception, across types of conscious content, are most likely to be found in midlatency negative-going brain responses in content-specific sensory areas.

Neuronal modeling of magnetoencephalography responses in auditory cortex to auditory and visual stimuli

Previous studies have demonstrated that auditory cortex activity can be influenced by crosssensory visual inputs. Intracortical recordings in non-human primates (NHP) have suggested a bottom-up feedforward (FF) type laminar profile for auditory evoked but top-down feedback (FB) type for cross-sensory visual evoked activity in the auditory cortex. To test whether this principle applies also to humans, we analyzed magnetoencephalography (MEG) responses from eight human subjects (six females) evoked by simple auditory or visual stimuli. In the estimated MEG source waveforms for auditory cortex region of interest, auditory evoked responses showed peaks at 37 and 90 ms and cross-sensory visual responses at 125 ms. The inputs to the auditory cortex were then modeled through FF and FB type connections targeting different cortical layers using the Human Neocortical Neurosolver (HNN), which consists of a neocortical circuit model linking the cellular- and circuit-level mechanisms to MEG. The HNN models suggested that the measured auditory response could be explained by an FF input followed by an FB input, and the crosssensory visual response by an FB input. Thus, the combined MEG and HNN results support the hypothesis that cross-sensory visual input in the auditory cortex is of FB type. The results also illustrate how the dynamic patterns of the estimated MEG/EEG source activity can provide information about the characteristics of the input into a cortical area in terms of the hierarchical organization among areas.

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.

Stimulus awareness is associated with secondary somatosensory cortex activation in an inattentional numbness paradigm

While inattentional blindness and deafness studies have revealed neural correlates of consciousness (NCC) without the confound of task relevance in the visual and auditory modality, comparable studies for the somatosensory modality are lacking. Here, we investigated NCC using functional magnetic resonance imaging (fMRI) in an inattentional numbness paradigm. Participants (N = 44) received weak electrical stimulation on the left hand while solving a demanding visual task. Half of the participants were informed that task-irrelevant weak tactile stimuli above the detection threshold would be applied during the experiment, while the other half expected stimuli below the detection threshold. Unexpected awareness assessments after the experiment revealed that altogether 10 participants did not consciously perceive the somatosensory stimuli during the visual task. Awareness was not significantly modulated by prior information. The fMRI data show that awareness of stimuli led to increased activation in the contralateral secondary somatosensory cortex. We found no significant effects of stimulus awareness in the primary somatosensory cortex or frontoparietal areas. Thus, our results support the hypothesis that somatosensory stimulus awareness is mainly based on activation in higher areas of the somatosensory cortex and does not require strong activation in extended anterior or posterior networks, which is usually seen when perceived stimuli are task-relevant.

Cortical origin of theta error signals

A multi-scale approach elucidated the origin of the error-related-negativity (ERN), with its associated theta-rhythm, and the post-error-positivity (Pe) in macaque supplementary eye field (SEF). Using biophysical modeling, synaptic inputs to a subpopulation of layer-3 (L3) and layer-5 (L5) pyramidal cells (PCs) were optimized to reproduce error-related spiking modulation and inter-spike intervals. The intrinsic dynamics of dendrites in L5 but not L3 error PCs generate theta rhythmicity with random phases. Saccades synchronized the phases of the theta-rhythm, which was magnified on errors. Contributions from error PCs to the laminar current source density (CSD) observed in SEF were negligible and could not explain the observed association between error-related spiking modulation in L3 PCs and scalp-EEG. CSD from recorded laminar field potentials in SEF was comprised of multipolar components, with monopoles indicating strong electro-diffusion, dendritic/axonal electrotonic current leakage outside SEF, or violations of the model assumptions. Our results also demonstrate the involvement of secondary cortical regions, in addition to SEF, particularly for the later Pe component. The dipolar component from the observed CSD paralleled the ERN dynamics, while the quadrupolar component paralleled the Pe. These results provide the most advanced explanation to date of the cellular mechanisms generating the ERN.

Biophysical modeling of frontocentral ERP generation links circuit-level mechanisms of action-stopping to a behavioral race model

Human frontocentral event-related potentials (FC-ERPs) are ubiquitous neural correlates of cognition and control, but their generating multiscale mechanisms remain mostly unknown. We used the Human Neocortical Neurosolver(HNN)'s biophysical model of a canonical neocortical circuit under exogenous thalamic and cortical drive to simulate the cell and circuit mechanisms underpinning the P2, N2, and P3 features of the FC-ERP observed after Stop-Signals in the Stop-Signal task (SST). We demonstrate that a sequence of simulated external thalamocortical and cortico-cortical drives can produce the FC-ERP, similar to what has been shown for primary sensory cortices. We used this model of the FC-ERP to examine likely circuit-mechanisms underlying FC-ERP features that distinguish between successful and failed action-stopping. We also tested their adherence to the predictions of the horse-race model of the SST, with specific hypotheses motivated by theoretical links between the P3 and Stop process. These simulations revealed that a difference in P3 onset between successful and failed Stops is most likely due to a later arrival of thalamocortical drive in failed Stops, rather than, for example, a difference in effective strength of the input. In contrast, the same model predicted that early thalamocortical drives underpinning the P2 and N2 differed in both strength and timing across stopping accuracy conditions. Overall, this model generates novel testable predictions of the thalamocortical dynamics underlying FC-ERP generation during action-stopping. Moreover, it provides a detailed cellular and circuit-level interpretation that supports links between these macroscale signatures and predictions of the behavioral race model.

Laminar specificity of the auditory perceptual awareness negativity: A biophysical modeling study

How perception of sensory stimuli emerges from brain activity is a fundamental question of neuroscience. To date, two disparate lines of research have examined this question. On one hand, human neuroimaging studies have helped us understand the large-scale brain dynamics of perception. On the other hand, work in animal models (mice, typically) has led to fundamental insight into the micro-scale neural circuits underlying perception. However, translating such fundamental insight from animal models to humans has been challenging. Here, using biophysical modeling, we show that the auditory awareness negativity (AAN), an evoked response associated with perception of target sounds in noise, can be accounted for by synaptic input to the supragranular layers of auditory cortex (AC) that is present when target sounds are heard but absent when they are missed. This additional input likely arises from cortico-cortical feedback and/or non-lemniscal thalamic projections and targets the apical dendrites of layer-5 (L5) pyramidal neurons. In turn, this leads to increased local field potential activity, increased spiking activity in L5 pyramidal neurons, and the AAN. The results are consistent with current cellular models of conscious processing and help bridge the gap between the macro and micro levels of perception-related brain activity.

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.

Laminar Specificity of the Auditory Perceptual Awareness Negativity: A Biophysical Modeling Study

How perception of sensory stimuli emerges from brain activity is a fundamental question of neuroscience. To date, two disparate lines of research have examined this question. On one hand, human neuroimaging studies have helped us understand the large-scale brain dynamics of perception. On the other hand, work in animal models (mice, typically) has led to fundamental insight into the micro-scale neural circuits underlying perception. However, translating such fundamental insight from animal models to humans has been challenging. Here, using biophysical modeling, we show that the auditory awareness negativity (AAN), an evoked response associated with perception of target sounds in noise, can be accounted for by synaptic input to the supragranular layers of auditory cortex (AC) that is present when target sounds are heard but absent when they are missed. This additional input likely arises from cortico-cortical feedback and/or non-lemniscal thalamic projections and targets the apical dendrites of layer-V pyramidal neurons (PNs). In turn, this leads to increased local field potential activity, increased spiking activity in layer-V PNs, and the AAN. The results are consistent with current cellular models of conscious processing and help bridge the gap between the macro and micro levels of perception-related brain activity.

Author Summary

To date, our understanding of the brain basis of conscious perception has mostly been restricted to large-scale, network-level activity that can be measured non-invasively in human subjects. However, we lack understanding of how such network-level activity is supported by individual neurons and neural circuits. This is at least partially because conscious perception is difficult to study in experimental animals, where such detailed characterization of neural activity is possible. To address this gap, we used biophysical modeling to gain circuit-level insight into an auditory brain response known as the auditory awareness negativity (AAN). This response can be recorded non-invasively in humans and is associated with perceptual awareness of sounds of interest. Our model shows that the AAN likely arises from specific cortical layers and cell types. These data help bridge the gap between circuit- and network-level theories of consciousness, and could lead to new, targeted treatments for perceptual dysfunction and disorders of consciousness.

Stereo-EEG recordings extend known distributions of canonical movement-related oscillations

Objective.Previous electrophysiological research has characterized canonical oscillatory patterns associated with movement mostly from recordings of primary sensorimotor cortex. Less work has attempted to decode movement based on electrophysiological recordings from a broader array of brain areas such as those sampled by stereoelectroencephalography (sEEG), especially in humans. We aimed to identify and characterize different movement-related oscillations across a relatively broad sampling of brain areas in humans and if they extended beyond brain areas previously associated with movement.Approach.We used a linear support vector machine to decode time-frequency spectrograms time-locked to movement, and we validated our results with cluster permutation testing and common spatial pattern decoding.Main results.We were able to accurately classify sEEG spectrograms during a keypress movement task versus the inter-trial interval. Specifically, we found these previously-described patterns: beta (13-30 Hz) desynchronization, beta synchronization (rebound), pre-movement alpha (8-15 Hz) modulation, a post-movement broadband gamma (60-90 Hz) increase and an event-related potential. These oscillatory patterns were newly observed in a wide range of brain areas accessible with sEEG that are not accessible with other electrophysiology recording methods. For example, the presence of beta desynchronization in the frontal lobe was more widespread than previously described, extending outside primary and secondary motor cortices.Significance.Our classification revealed prominent time-frequency patterns which were also observed in previous studies that used non-invasive electroencephalography and electrocorticography, but here we identified these patterns in brain regions that had not yet been associated with movement. This provides new evidence for the anatomical extent of the system of putative motor networks that exhibit each of these oscillatory patterns.

Resolving the mesoscopic missing link: Biophysical modeling of EEG from cortical columns in primates

Event-related potentials (ERP) are among the most widely measured indices for studying human cognition. While their timing and magnitude provide valuable insights, their usefulness is limited by our understanding of their neural generators at the circuit level. Inverse source localization offers insights into such generators, but their solutions are not unique. To address this problem, scientists have assumed the source space generating such signals comprises a set of discrete equivalent current dipoles, representing the activity of small cortical regions. Based on this notion, theoretical studies have employed forward modeling of scalp potentials to understand how changes in circuit-level dynamics translate into macroscopic ERPs. However, experimental validation is lacking because it requires in vivo measurements of intracranial brain sources. Laminar local field potentials (LFP) offer a mechanism for estimating intracranial current sources. Yet, a theoretical link between LFPs and intracranial brain sources is missing. Here, we present a forward modeling approach for estimating mesoscopic intracranial brain sources from LFPs and predict their contribution to macroscopic ERPs. We evaluate the accuracy of this LFP-based representation of brain sources utilizing synthetic laminar neurophysiological measurements and then demonstrate the power of the approach in vivo to clarify the source of a representative cognitive ERP component. To that end, LFP was measured across the cortical layers of visual area V4 in macaque monkeys performing an attention demanding task. We show that area V4 generates dipoles through layer-specific transsynaptic currents that biophysically recapitulate the ERP component through the detailed forward modeling. The constraints imposed on EEG production by this method also revealed an important dissociation between computational and biophysical contributors. As such, this approach represents an important bridge between laminar microcircuitry, through the mesoscopic activity of cortical columns to the patterns of EEG we measure at the scalp.

Non-zero mean alpha oscillations revealed with computational model and empirical data

Ongoing oscillations and evoked responses are two main types of neuronal activity obtained with diverse electrophysiological recordings (EEG/MEG/iEEG/LFP). Although typically studied separately, they might in fact be closely related. One possibility to unite them is to demonstrate that neuronal oscillations have non-zero mean which predicts that stimulus- or task-triggered amplitude modulation of oscillations can contribute to the generation of evoked responses. We validated this mechanism using computational modelling and analysis of a large EEG data set. With a biophysical model, we indeed demonstrated that intracellular currents in the neuron are asymmetric and, consequently, the mean of alpha oscillations is non-zero. To understand the effect that neuronal currents exert on oscillatory mean, we varied several biophysical and morphological properties of neurons in the network, such as voltage-gated channel densities, length of dendrites, and intensity of incoming stimuli. For a very large range of model parameters, we observed evidence for non-zero mean of oscillations. Complimentary, we analysed empirical rest EEG recordings of 90 participants (50 young, 40 elderly) and, with spatio-spectral decomposition, detected at least one spatially-filtred oscillatory component of non-zero mean alpha oscillations in 93% of participants. In order to explain a complex relationship between the dynamics of amplitude-envelope and corresponding baseline shifts, we performed additional simulations with simple oscillators coupled with different time delays. We demonstrated that the extent of spatial synchronisation may obscure macroscopic estimation of alpha rhythm modulation while leaving baseline shifts unchanged. Overall, our results predict that amplitude modulation of neural oscillations should at least partially explain the generation of evoked responses. Therefore, inference about changes in evoked responses with respect to cognitive conditions, age or neuropathologies should be constructed while taking into account oscillatory neuronal dynamics.

Thalamocortical Mechanisms Regulating the Relationship between Transient Beta Events and Human Tactile Perception

Transient neocortical events with high spectral power in the 15-29 Hz beta band are among the most reliable predictors of sensory perception. Prestimulus beta event rates in primary somatosensory cortex correlate with sensory suppression, most effectively 100-300 ms before stimulus onset. However, the neural mechanisms underlying this perceptual association are unknown. We combined human magnetoencephalography (MEG) measurements with biophysical neural modeling to test potential cellular and circuit mechanisms that underlie observed correlations between prestimulus beta events and tactile detection. Extending prior studies, we found that simulated bursts from higher-order, nonlemniscal thalamus were sufficient to drive beta event generation and to recruit slow supragranular inhibition acting on a 300 ms timescale to suppress sensory information. Further analysis showed that the same beta-generating mechanism can lead to facilitated perception for a brief period when beta events occur simultaneously with tactile stimulation before inhibition is recruited. These findings were supported by close agreement between model-derived predictions and empirical MEG data. The postevent suppressive mechanism explains an array of studies that associate beta with decreased processing, whereas the during-event facilitatory mechanism may demand a reinterpretation of the role of beta events in the context of coincident timing.

Somatosensory Plasticity in Hemiplegic Cerebral Palsy Following Constraint Induced Movement Therapy

BACKGROUND

Children with hemiplegic cerebral palsy (HCP) experience upper limb somatosensory and motor deficits. Although constraint-induced movement therapy (CIMT) improves motor function, its impact on somatosensory function remains underinvestigated.

OBJECTIVE

The objective of this study was to evaluate somatosensory perception and related brain responses in children with HCP, before and after a somatosensory enhanced CIMT protocol, as measured using clinical sensory and motor assessments and magnetoencephalography.

METHODS

Children with HCP attended a somatosensory enhanced CIMT camp. Clinical somatosensory (tactile registration, 2-point discrimination, stereognosis, proprioception, kinesthesia) and motor outcomes (Quality of Upper Extremity Skills [QUEST] Total/Grasp, Jebsen-Taylor Hand Function Test, grip strength, Assisting Hand Assessment), as well as latency and amplitude of magnetoencephalography somatosensory evoked fields (SEF), were assessed before and after the CIMT camp with paired sample t-tests or Wilcoxon signed-rank tests.

RESULTS

Twelve children with HCP (mean age: 7.5 years, standard deviation: 2.4) participated. Significant improvements in tactile registration for the affected (hemiplegic) hand (Z = 2.39, P = 0.02) were observed in addition to statistically and clinically significant improvements in QUEST total (t = 3.24, P = 0.007), QUEST grasp (t = 3.24, P = 0.007), Assisting Hand Assessment (Z = 2.25, P = 0.03), and Jebsen-Taylor Hand Function Test (t = -2.62, P = 0.03). A significant increase in the SEF peak amplitude was also found in the affected hand 100 ms after stimulus onset (t = -2.22, P = 0.04).

CONCLUSIONS

Improvements in somatosensory clinical function and neural processing in the affected primary somatosensory cortex in children with HCP were observed after a somatosensory enhanced CIMT program. Further investigation is warranted to continue to evaluate the effectiveness of a sensory enhanced CIMT program in larger samples and controlled study designs.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

Neural Basis of Somatosensory Spatial and Temporal Discrimination in Humans: The Role of Sensory Detection

While detecting somatic stimuli from the external environment, an accurate determination of their spatial and temporal properties is essential for human behavior. Whether and how detection relates to human capacity for somatosensory spatial discrimination (SD) and temporal discrimination (TD) remains unclear. Here, participants underwent functional magnetic resonance imaging scanning when simply detecting vibrotactile stimuli of the leg, judging their location (SD), or deciding their number in time (TD). By conceptualizing tactile discrimination as consisting of detection and determination processes, we found that tactile detection elicited activation specifically involved in SD within the right inferior and superior parietal lobules, 2 regions previously implicated in the control of spatial attention. These 2 regions remained activated in the determination process, during which functional connectivity between these 2 regions predicted individual SD ability. In contrast, tactile detection produced little activation specifically related to TD. Participants' TD ability was implemented in brain regions implicated in coding temporal structures of somatic stimuli (primary somatosensory cortex) and time estimation (anterior cingulate, pre-supplementary motor area, and putamen). Together, our findings indicate a close link between somatosensory detection and SD (but not TD) at the neural level, which aids in explaining why we can promptly respond toward detected somatic stimuli.

Pre-stimulus alpha oscillation and post-stimulus cortical activity differ in localization between consciously perceived and missed near-threshold somatosensory stimuli

Conscious perception of a near‐threshold (NT) stimulus is characterized by the pre‐ and post‐stimulus brain state. However, the power of pre‐stimulus neural oscillations and strength of post‐stimulus cortical activity that lead to conscious perception have rarely been examined in individual cortical areas. This is because most previous electro‐ and magnetoencephalography (EEG and MEG, respectively) studies involved scalp‐ and sensor‐level analyses. Therefore, we recorded MEG during a continuous NT somatosensory stimulus detection task and applied the reconstructed source data in order to identify cortical areas where the post‐stimulus cortical activity and pre‐stimulus alpha oscillation predict the conscious perception of NT somatosensory stimuli. We found that the somatosensory hierarchical processing areas, prefrontal areas and cortical areas belonging to the default mode network showed stronger cortical activity for consciously perceived trials in the post‐stimulus period, but the cortical activity in primary somatosensory area (SI) is independent of conscious perception during the early stage of NT stimulus processing. In addition, we revealed that the pre‐stimulus alpha oscillation only in SI is predictive of conscious perception. These findings suggest that the bottom‐up stream of somatosensory information flow following SI and pre‐stimulus alpha activity fluctuation in SI as a top‐down modulation are crucial constituents of conscious perception.

Action-locked Neural Responses in Auditory Cortex to Self-generated Sounds

Sensory perception is a product of interactions between the internal state of an organism and the physical attributes of a stimulus. It has been shown across the animal kingdom that perception and sensory-evoked physiological responses are modulated depending on whether or not the stimulus is the consequence of voluntary actions. These phenomena are often attributed to motor signals sent to relevant sensory regions that convey information about upcoming sensory consequences. However, the neurophysiological signature of action-locked modulations in sensory cortex, and their relationship with perception, is still unclear. In the current study, we recorded neurophysiological (using Magnetoencephalography) and behavioral responses from 16 healthy subjects performing an auditory detection task of faint tones. Tones were either generated by subjects' voluntary button presses or occurred predictably following a visual cue. By introducing a constant temporal delay between button press/cue and tone delivery, and applying source-level analysis, we decoupled action-locked and auditory-locked activity in auditory cortex. We show action-locked evoked-responses in auditory cortex following sound-triggering actions and preceding sound onset. Such evoked-responses were not found for button-presses that were not coupled with sounds, or sounds delivered following a predictive visual cue. Our results provide evidence for efferent signals in human auditory cortex that are locked to voluntary actions coupled with future auditory consequences.

Perceptual awareness negativity: a physiological correlate of sensory consciousness

Much research on the neural correlates of consciousness (NCC) has focused on two evoked potentials, the P3b and the visual or auditory awareness negativity (VAN, AAN). Surveying a broad range of recent experimental evidence, we find that repeated failures to observe the P3b during conscious perception eliminate it as a putative NCC. Neither the VAN nor the AAN have been dissociated from consciousness; furthermore, a similar neural signal correlates with tactile consciousness. These awareness negativities can be maximal contralateral to the evoking stimulus, are likely generated in underlying sensory cortices, and point to the existence of a generalized perceptual awareness negativity (PAN) reflecting the onset of sensory consciousness.

Dissociating Perceptual Awareness and Postperceptual Processing: The P300 Is Not a Reliable Marker of Somatosensory Target Detection

A central challenge in the study of conscious perception lies in dissociating the neural correlates of perceptual awareness from those reflecting its precursors and consequences. No-report paradigms have been instrumental in this endeavor, demonstrating that the event-related potential P300, recorded from the human scalp, reflects reports rather than awareness. However, these paradigms cannot probe the degree to which stimuli are consciously processed from trial to trial and, thus, leave open the possibility that the P300 is a genuine correlate of conscious access enabling reports. Here, instead of removing report requirements, we took the opposite approach and equated postperceptual task demands across conscious and unconscious trials by orthogonalizing target detection and overt reports in a somatosensory detection task. We used Bayesian model selection to track the transformation from physical to perceptual processing stages in the EEG data of 24 male and female participants and show that the early P50 component scaled with physical stimulus intensity, whereas the N140 component was the first correlate of target detection. The late P300 component was elicited for both perceived and unperceived stimuli and was not substantially modulated by target detection. This was in stark contrast to a control experiment using a classical direct report task, which replicated the P50 and N140 effects but additionally showed a strong effect of target detection in the P300 time range. Our results demonstrate the task dependence of the P300 in the somatosensory modality and show that late cortical potentials dissociate from perceptual awareness even when stimuli are always reported.SIGNIFICANCE STATEMENT The time it takes for sensory information to enter our conscious experience can be an indicator of the neural processing stages that lead to perceptual awareness. However, because many cognitive processes routinely correlate with perception, isolating those signals that uniquely reflect perceptual awareness is not a trivial task. Here, we show that late electroencephalography signals cease to correlate with somatosensory awareness when common task confounds are controlled. Importantly, by balancing report requirements instead of abolishing them, we show that the lack of late effects cannot be explained by a lack of conscious access. Instead, we propose that conscious access occurs earlier, at ∼150 ms, supporting the view that early activity in sensory cortices is a neural correlate of conscious perception.

Neural Mechanisms Underlying Human Auditory Evoked Responses Revealed By Human Neocortical Neurosolver

Auditory evoked fields (AEFs) are commonly studied, yet their underlying neural mechanisms remain poorly understood. Here, we used the biophysical modelling software Human Neocortical Neurosolver (HNN) whose foundation is a canonical neocortical circuit model to interpret the cell and network mechanisms contributing to macroscale AEFs elicited by a simple tone, measured with magnetoencephalography. We found that AEFs can be reproduced by activating the neocortical circuit through a layer specific sequence of feedforward and feedback excitatory synaptic drives, similar to prior simulation of somatosensory evoked responses, supporting the notion that basic structures and activation patterns are preserved across sensory regions. We also applied the modeling framework to develop and test predictions on neural mechanisms underlying AEF differences in the left and right hemispheres, as well as in hemispheres contralateral and ipsilateral to the presentation of the auditory stimulus. We found that increasing the strength of the excitatory synaptic cortical feedback inputs to supragranular layers simulates the commonly observed right hemisphere dominance, while decreasing the input latencies and simultaneously increasing the number of cells contributing to the signal accounted for the contralateral dominance. These results provide a direct link between human data and prior animal studies and lay the foundation for future translational research examining the mechanisms underlying alteration in this fundamental biomarker of auditory processing in healthy cognition and neuropathology.

Neural surprise in somatosensory Bayesian learning

Tracking statistical regularities of the environment is important for shaping human behavior and perception. Evidence suggests that the brain learns environmental dependencies using Bayesian principles. However, much remains unknown about the employed algorithms, for somesthesis in particular. Here, we describe the cortical dynamics of the somatosensory learning system to investigate both the form of the generative model as well as its neural surprise signatures. Specifically, we recorded EEG data from 40 participants subjected to a somatosensory roving-stimulus paradigm and performed single-trial modeling across peri-stimulus time in both sensor and source space. Our Bayesian model selection procedure indicates that evoked potentials are best described by a non-hierarchical learning model that tracks transitions between observations using leaky integration. From around 70ms post-stimulus onset, secondary somatosensory cortices are found to represent confidence-corrected surprise as a measure of model inadequacy. Indications of Bayesian surprise encoding, reflecting model updating, are found in primary somatosensory cortex from around 140ms. This dissociation is compatible with the idea that early surprise signals may control subsequent model update rates. In sum, our findings support the hypothesis that early somatosensory processing reflects Bayesian perceptual learning and contribute to an understanding of its underlying mechanisms.

Our environment features statistical regularities, such as a drop of rain predicting imminent rainfall. Despite the importance for behavior and survival, much remains unknown about how these dependencies are learned, particularly for somatosensation. As surprise signalling about novel observations indicates a mismatch between one’s beliefs and the world, it has been hypothesized that surprise computation plays an important role in perceptual learning. By analyzing EEG data from human participants receiving sequences of tactile stimulation, we compare different formulations of surprise and investigate the employed underlying learning model. Our results indicate that the brain estimates transitions between observations. Furthermore, we identified different signatures of surprise computation and thereby provide a dissociation of the neural correlates of belief inadequacy and belief updating. Specifically, early surprise responses from around 70ms were found to signal the need for changes to the model, with encoding of its subsequent updating occurring from around 140ms. These results provide insights into how somatosensory surprise signals may contribute to the learning of environmental statistics.

Vibration Alert to the Brain: Evoked and Induced MEG Responses to High-Frequency Vibrotactile Stimuli on the Index Finger of Dominant and Non-dominant Hand

Background: In recent years, vibrotactile haptic feedback technology has been widely used for user interfaces in the mobile devices. Although functional neuroimaging studies have investigated human brain responses to different types of tactile inputs, the neural mechanisms underlying high-frequency vibrotactile perception are still relatively unknown. Our aim was to investigate neuromagnetic brain responses to high-frequency vibrotactile stimulation, using magnetoencephalography (MEG).

Methods: We measured 152-channel whole-head MEG in 30 healthy, right-handed volunteers (aged 20–28 years, 15 females). A total of 300 vibrotactile stimuli were presented at the tip of either the left index finger or the right index finger in two separate sessions. Sinusoidal vibrations at 150 Hz for 200 ms were generated with random inter-stimulus intervals between 1.6 and 2.4 s. Both time-locked analysis and time-frequency analysis were performed to identify peak responses and oscillatory modulations elicited by high-frequency vibrations. The significance of the evoked and induced responses for dominant and non-dominant hand stimulation conditions was statistically tested, respectively. The difference in responses between stimulation conditions was also statistically evaluated.

Results: Prominent peak responses were observed at 56 ms (M50) and at 100 ms (M100) for both stimulation conditions. The M50 response revealed clear dipolar field patterns in the contralateral side with significant cortical activations in the contralateral primary sensorimotor area, whereas the M100 response was not as prominent as the M50. Vibrotactile stimulation induced significant suppression of both alpha (8–12 Hz) and beta (20–30 Hz) band activity during the mid-latency period (0.2–0.4 s), primarily in sensorimotor areas contralateral to the stimulation side. In addition, a significant alpha enhancement effect in posterior regions was accompanied with alpha suppressions in sensorimotor regions. The alpha suppression was observed in a broader distribution of cortical areas for the non-dominant hand stimulation.

Conclusion: Our data demonstrate that high-frequency tactile vibrations, which is known to primarily activate Pacinian corpuscles, elicit somatosensory M50 and M100 responses in the evoked fields and induce modulations of alpha and beta band oscillations during mid-latency periods. Our study is also consistent with that the primary sensorimotor area is significantly involved in the processing of high-frequency vibrotactile information with contralateral dominance.

Differences in visually induced MEG oscillations reflect differences in deep cortical layer activity

Neural activity is organized at multiple scales, ranging from the cellular to the whole brain level. Connecting neural dynamics at different scales is important for understanding brain pathology. Neurological diseases and disorders arise from interactions between factors that are expressed in multiple scales. Here, we suggest a new way to link microscopic and macroscopic dynamics through combinations of computational models. This exploits results from statistical decision theory and Bayesian inference. To validate our approach, we used two independent MEG datasets. In both, we found that variability in visually induced oscillations recorded from different people in simple visual perception tasks resulted from differences in the level of inhibition specific to deep cortical layers. This suggests differences in feedback to sensory areas and each subject's hypotheses about sensations due to differences in their prior experience. Our approach provides a new link between non-invasive brain imaging data, laminar dynamics and top-down control.

Biophysically detailed forward modeling of the neural origin of EEG and MEG signals

Electroencephalography (EEG) and magnetoencephalography (MEG) are among the most important techniques for non-invasively studying cognition and disease in the human brain. These signals are known to originate from cortical neural activity, typically described in terms of current dipoles. While the link between cortical current dipoles and EEG/MEG signals is relatively well understood, surprisingly little is known about the link between different kinds of neural activity and the current dipoles themselves. Detailed biophysical modeling has played an important role in exploring the neural origin of intracranial electric signals, like extracellular spikes and local field potentials. However, this approach has not yet been taken full advantage of in the context of exploring the neural origin of the cortical current dipoles that are causing EEG/MEG signals. Here, we present a method for reducing arbitrary simulated neural activity to single current dipoles. We find that the method is applicable for calculating extracranial signals, but less suited for calculating intracranial electrocorticography (ECoG) signals. We demonstrate that this approach can serve as a powerful tool for investigating the neural origin of EEG/MEG signals. This is done through example studies of the single-neuron EEG contribution, the putative EEG contribution from calcium spikes, and from calculating EEG signals from large-scale neural network simulations. We also demonstrate how the simulated current dipoles can be used directly in combination with detailed head models, allowing for simulated EEG signals with an unprecedented level of biophysical details. In conclusion, this paper presents a framework for biophysically detailed modeling of EEG and MEG signals, which can be used to better our understanding of non-inasively measured neural activity in humans.

Assessment of haptic memory using somatosensory change-related cortical responses

Haptic memory briefly retains somatosensory information for later use; however, how and which cortical areas are affected by haptic memory remain unclear. We used change-related cortical responses to investigate the relationship between the somatosensory cortex and haptic memory objectively. Electrical pulses, at 50 Hz with a duration of 500 ms, were randomly applied to the second, third, and fourth fingers of the right and left hands at an even probability every 800 ms. Each stimulus was labeled as D (preceded by a different side) or S (preceded by the same side). The D stimuli were further classified into 1D, 2D, and 3D, according to the number of different preceding stimuli. The S stimuli were similarly divided into 1S and 2S. The somatosensory-evoked magnetic fields obtained were divided into four components via a dipole analysis, and each component's amplitudes were measured using the source strength waveform. The results showed that the preceding event did not affect the amplitude of the earliest 20-30 ms response in the primary somatosensory cortex. However, in the subsequent three components, the cortical activity amplitude was largest in 3D, followed by 2D, 1D, and S. These results indicate that such modulatory effects occurred somewhere in the somatosensory processing pathway higher than Brodmann's area 3b. To the best of our knowledge, this is the first study to demonstrate the existence of haptic memory for somatosensory laterality and its impact on the somatosensory cortex using change-related cortical responses without contamination from peripheral effects.

Statistical decision theory and multiscale analyses of human brain data

BACKGROUND

In the era of Big Data, large scale electrophysiological data from animal and human studies are abundant. These data contain information at multiple spatiotemporal scales. However, current approaches for the analysis of electrophysiological data often focus on a single spatiotemporal scale only.

NEW METHOD

We discuss a multiscale approach for the analysis of electrophysiological data. This is based on combining neural models that describe brain data at different scales. It allows us to make laminar-specific inferences about neurobiological properties of cortical sources using non invasive human electrophysiology data.

RESULTS

We provide a mathematical proof of this approach using statistical decision theory. We also consider its extensions to brain imaging studies including data from the same subjects performing different tasks. As an illustration, we show that changes in gamma oscillations between different people might originate from differences in recurrent connection strengths of inhibitory interneurons in layers 5/6.

COMPARISON WITH EXISTING METHODS

This is a new approach that follows up on our recent work. It is different from other approaches where the scale of spatiotemporal dynamics is fixed.

CONCLUSIONS

We discuss a multiscale approach for the analysis of human MEG data. This uses a neural mass model that includes constraints informed by a compartmental model. This has two advantages. First, it allows us to find differences in cortical laminar dynamics and understand neurobiological properties like neuromodulation, excitation to inhibition balance etc. using non invasive data. Second, it allows us to validate macroscale models by exploiting animal data.

Older Age Increases the Amplitude of Muscle Stretch-Induced Cortical Beta-Band Suppression But Does not Affect Rebound Strength

Healthy aging is associated with deterioration of the sensorimotor system, which impairs balance and somatosensation. However, the exact age-related changes in the cortical processing of sensorimotor integration are unclear. This study investigated primary sensorimotor cortex (SM1) oscillations in the 15-30 Hz beta band at rest and following (involuntary) rapid stretches to the triceps surae muscles (i.e., proprioceptive stimulation) of young and older adults. A custom-built, magnetoencephalography (MEG)-compatible device was used to deliver rapid (190°·s^-1) ankle rotations as subjects sat passively in a magnetically-shielded room while MEG recorded their cortical signals. Eleven young (age 25 ± 3 years) and 12 older (age 70 ± 3 years) adults matched for physical activity level demonstrated clear 15-30 Hz beta band suppression and rebound in response to the stretches. A sub-sample (10 young and nine older) were tested for dynamic balance control on a sliding platform. Older adults had greater cortical beta power pre-stretch (e.g., right leg: 4.0 ± 1.6 fT vs. 5.6 ± 1.7 fT, P = 0.044) and, subsequently, greater normalized movement-related cortical beta suppression post-proprioceptive stimulation (e.g., right leg: -5.8 ± 1.3 vs. -7.6 ± 1.7, P = 0.01) than young adults. Furthermore, poorer balance was associated with stronger cortical beta suppression following proprioceptive stimulation (r = -0.478, P = 0.038, n = 19). These results provide further support that cortical processing of proprioception is hindered in older adults, potentially (adversely) influencing sensorimotor integration. This was demonstrated by the impairment of prompt motor action control, i.e., regaining perturbed balance. Finally, SM1 cortex beta suppression to a proprioceptive stimulus seems to indicate poorer sensorimotor functioning in older adults.

Human Neocortical Neurosolver (HNN), a new software tool for interpreting the cellular and network origin of human MEG/EEG data

Magneto- and electro-encephalography (MEG/EEG) non-invasively record human brain activity with millisecond resolution providing reliable markers of healthy and disease states. Relating these macroscopic signals to underlying cellular- and circuit-level generators is a limitation that constrains using MEG/EEG to reveal novel principles of information processing or to translate findings into new therapies for neuropathology. To address this problem, we built Human Neocortical Neurosolver (HNN, https://hnn.brown.edu) software. HNN has a graphical user interface designed to help researchers and clinicians interpret the neural origins of MEG/EEG. HNN's core is a neocortical circuit model that accounts for biophysical origins of electrical currents generating MEG/EEG. Data can be directly compared to simulated signals and parameters easily manipulated to develop/test hypotheses on a signal's origin. Tutorials teach users to simulate commonly measured signals, including event related potentials and brain rhythms. HNN's ability to associate signals across scales makes it a unique tool for translational neuroscience research.

Neuronal correlates of full and partial visual conscious perception

Stimuli may induce only partial consciousness-an intermediate between null and full consciousness-where the presence but not identity of an object can be reported. The differences in the neuronal basis of full and partial consciousness are poorly understood. We investigated if evoked and oscillatory activity could dissociate full from partial conscious perception. We recorded human cortical activity with magnetoencephalography (MEG) during a visual perception task in which stimulus could be either partially or fully perceived. Partial consciousness was associated with an early increase in evoked activity and theta/low-alpha-band oscillations while full consciousness was also associated with late evoked activity and beta-band oscillations. Full from partial consciousness was dissociated by stronger evoked activity and late increase in theta oscillations that were localized to higher-order visual regions and posterior parietal and prefrontal cortices. Our results reveal both evoked activity and theta oscillations dissociate partial and full consciousness.

Somatosensory Cortex Efficiently Processes Touch Located Beyond the Body

The extent to which a tool is an extension of its user is a question that has fascinated writers and philosophers for centuries [1]. Despite two decades of research [2-7], it remains unknown how this could be instantiated at the neural level. To this aim, the present study combined behavior, electrophysiology and neuronal modeling to characterize how the human brain could treat a tool like an extended sensory "organ." As with the body, participants localize touches on a hand-held tool with near-perfect accuracy [7]. This behavior is owed to the ability of the somatosensory system to rapidly and efficiently use the tool as a tactile extension of the body. Using electroencephalography (EEG), we found that where a hand-held tool was touched was immediately coded in the neural dynamics of primary somatosensory and posterior parietal cortices of healthy participants. We found similar neural responses in a proprioceptively deafferented patient with spared touch perception, suggesting that location information is extracted from the rod's vibrational patterns. Simulations of mechanoreceptor responses [8] suggested that the speed at which these patterns are processed is highly efficient. A second EEG experiment showed that touches on the tool and arm surfaces were localized by similar stages of cortical processing. Multivariate decoding algorithms and cortical source reconstruction provided further evidence that early limb-based processes were repurposed to map touch on a tool. We propose that an elementary strategy the human brain uses to sense with tools is to recruit primary somatosensory dynamics otherwise devoted to the body.

Persistent Gamma Spiking in SI Nonsensory Fast Spiking Cells Predicts Perceptual Success

Gamma oscillations (30-55 Hz) are hypothesized to temporally coordinate sensory encoding, enabling perception. However, fast spiking interneurons (FS), key gamma generators, can be highly sensory responsive, as is the gamma band local field potential (LFP). How can FS-mediated gamma act as an impartial temporal reference for sensory encoding, when the sensory drive itself presumably perturbs the pre-established rhythm? Combining tetrode recording in SI barrel cortex with controlled psychophysics, we found a unique FS subtype that was not sensory responsive and spiked regularly at gamma range intervals (gamma regular nonsensory FS [grnsFS]). Successful detection was predicted by a further increase in gamma regular spiking of grnsFS, persisting from before to after sensory onset. In contrast, broadband LFP power, including gamma, negatively predicted detection and did not cohere with gamma band spiking by grnsFS. These results suggest that a distinct FS subtype mediates perceptually relevant oscillations, independent of the LFP and sensory drive.

Bayesian Modelling of Induced Responses and Neuronal Rhythms

Neural rhythms or oscillations are ubiquitous in neuroimaging data. These spectral responses have been linked to several cognitive processes; including working memory, attention, perceptual binding and neuronal coordination. In this paper, we show how Bayesian methods can be used to finesse the ill-posed problem of reconstructing-and explaining-oscillatory responses. We offer an overview of recent developments in this field, focusing on (i) the use of MEG data and Empirical Bayes to build hierarchical models for group analyses-and the identification of important sources of inter-subject variability and (ii) the construction of novel dynamic causal models of intralaminar recordings to explain layer-specific activity. We hope to show that electrophysiological measurements contain much more spatial information than is often thought: on the one hand, the dynamic causal modelling of non-invasive (low spatial resolution) electrophysiology can afford sub-millimetre (hyper-acute) resolution that is limited only by the (spatial) complexity of the underlying (dynamic causal) forward model. On the other hand, invasive microelectrode recordings (that penetrate different cortical layers) can reveal laminar-specific responses and elucidate hierarchical message passing and information processing within and between cortical regions at a macroscopic scale. In short, the careful and biophysically grounded modelling of sparse data enables one to characterise the neuronal architectures generating oscillations in a remarkable detail.

Neural basis of somatosensory target detection independent of uncertainty, relevance, and reports

Research on somatosensory awareness has yielded highly diverse findings with putative neural correlates ranging from activity within somatosensory cortex to activation of widely distributed frontoparietal networks. Divergent results from previous studies may reside in cognitive processes that often coincide with stimulus awareness in experimental settings. To scrutinise the specific relevance of regions implied in the target detection network, we used functional magnetic resonance imaging (n = 27) on a novel somatosensory detection task that explicitly controls for stimulus uncertainty, behavioural relevance, overt reports, and motor responses. Using Bayesian Model Selection, we show that responses reflecting target detection are restricted to secondary somatosensory cortex, whereas activity in insular, cingulate, and motor regions is best explained in terms of stimulus uncertainty and overt reports. Our results emphasise the role of sensory-specific cortex for the emergence of perceptual awareness and dissect the contribution of the frontoparietal network to classical detection tasks.

Induced and Evoked Properties of Vibrotactile Adaptation in the Primary Somatosensory Cortex

Prolonged exposure to afferent stimulation (“adaptation”) can cause profound short-term changes in the responsiveness of cortical sensory neurons. While several models have been proposed that link adaptation to single-neuron dynamics, including GABAergic inhibition, the process is currently imperfectly understood at the whole-brain level in humans. Here, we used magnetoencephalography (MEG) to examine the neurophysiological correlates of adaptation within SI in humans. In one condition, a 25 Hz adapting stimulus (5 s) was followed by a 1 s 25 Hz probe (“same”), and in a second condition, the adapting stimulus was followed by a 1 s 180 Hz probe (“different”). We hypothesized that changes in the mu-beta activity band (reflecting GABAergic processing) would be modulated differently between the “same” and “different” probe stimuli. We show that the primary somatosensory (SI) mu-beta response to the “same” probe is significantly reduced (p = 0.014) compared to the adapting stimulus, whereas the mu-beta response to the “different” probe is not (p = n.s.). This reduction may reflect sharpening of the spatiotemporal pattern of activity after adaptation. The stimulus onset mu-beta response did not differ between a 25 Hz adapting stimulus and a 180 Hz probe, suggesting that the mu-beta response is independent of stimulus frequency. Furthermore, we show a sustained evoked and induced desynchronization for the duration of the adapting stimulus, consistent with invasive studies. Our findings are important in understanding the neurophysiology underlying short-term and stimulus-induced plasticity in the human brain and shows that the brain response to tactile stimulation is altered after only brief stimulation.

A Finite Element Solution of the Forward Problem in EEG for Multipolar Sources

Multipolar source models have been presented in the context of electro/magnetoencephalography (E/MEG) to compensate for the limitations of the classical equivalent current dipole to represent realistic generators of brain activity. Although there exist several reports accounting for the advantages of multipolar components over single dipoles, there is still no available numerical implementation in fully personalized scenarios. In this paper, we present, for the first time, a finite element framework for simulating EEG signals generated by multipolar current sources in individualized, heterogeneous, and anisotropic head models. This formulation is based on the subtraction approach, guaranteeing the existence and uniqueness of the solution. In particular, we analyze the cases of monopolar, dipolar, and quadrupolar source components, for which we study their performance in idealized and realistic head models. Numerical solutions are compared with analytical formulas in multi-layered spherical models. Such formulas are available in the case of monopolar and dipolar sources, and here derived for the quadrupolar components. We finally illustrate their advantages in the description of extended current generators using a realistic head model. The framework presented here enables further analysis towards the estimation of biophysically principled source parameters from standard E/MEG experiments.

A Prospective Study of the Impact of Transcranial Alternating Current Stimulation on EEG Correlates of Somatosensory Perception

The (8-12 Hz) neocortical alpha rhythm is associated with shifts in attention across sensory systems, and is thought to represent a sensory gating mechanism for the inhibitory control of cortical processing. The present preliminary study sought to explore whether alpha frequency transcranial alternating current stimulation (tACS) could modulate endogenous alpha power in the somatosensory system, and whether the hypothesized modulation would causally impact perception of tactile stimuli at perceptual threshold. We combined electroencephalography (EEG) with simultaneous brief and intermittent tACS applied over primary somatosensory cortex at individuals' endogenous alpha frequency during a tactile detection task (n = 12 for EEG, n = 20 for behavior). EEG-measured pre-stimulus alpha power was higher on non-perceived than perceived trials, and analogous perceptual correlates emerged in early components of the tactile evoked response. Further, baseline normalized tactile detection performance was significantly lower during alpha than sham tACS, but the effect did not last into the post-tACS time period. Pre- to post-tACS changes in alpha power were linearly dependent upon baseline state, such that alpha power tended to increase when pre-tACS alpha power was low, and decrease when it was high. However, these observations were comparable in both groups, and not associated with evidence of tACS-induced alpha power modulation. Nevertheless, the tactile stimulus evoked response potential (ERP) revealed a potentially lasting impact of alpha tACS on circuit dynamics. The post-tACS ERP was marked by the emergence of a prominent peak ∼70 ms post-stimulus, which was not discernible post-sham, or in either pre-stimulation condition. Computational neural modeling designed to simulate macroscale EEG signals supported the hypothesis that the emergence of this peak could reflect synaptic plasticity mechanisms induced by tACS. The primary lesson learned in this study, which commanded a small sample size, was that while our experimental paradigm provided some evidence of an influence of tACS on behavior and circuit dynamics, it was not sufficient to induce observable causal effects of tACS on EEG-measured alpha oscillations. We discuss limitations and suggest improvements that may help further delineate a causal influence of tACS on cortical dynamics and perception in future studies.

Dynamic large-scale network synchronization from perception to action

Sensory-guided actions entail the processing of sensory information, generation of perceptual decisions, and the generation of appropriate actions. Neuronal activity underlying these processes is distributed into sensory, fronto-parietal, and motor brain areas, respectively. How the neuronal processing is coordinated across these brain areas to support functions from perception to action remains unknown. We investigated whether phase synchronization in large-scale networks coordinate these processes. We recorded human cortical activity with magnetoencephalography (MEG) during a task in which weak somatosensory stimuli remained unperceived or were perceived. We then assessed dynamic evolution of phase synchronization in large-scale networks from source-reconstructed MEG data by using advanced analysis approaches combined with graph theory. Here we show that perceiving and reporting of weak somatosensory stimuli is correlated with sustained strengthening of large-scale synchrony concurrently in delta/theta (3-7 Hz) and gamma (40-60 Hz) frequency bands. In a data-driven network localization, we found this synchronization to dynamically connect the task-relevant, that is, the fronto-parietal, sensory, and motor systems. The strength and temporal pattern of interareal synchronization were also correlated with the response times. These data thus show that key brain areas underlying perception, decision-making, and actions are transiently connected by large-scale dynamic phase synchronization in the delta/theta and gamma bands.

A Blind Module Identification Approach for Predicting Effective Connectivity Within Brain Dynamical Subnetworks

Model-based network discovery measures, such as the brain effective connectivity, require fitting of generative process models to measurements obtained from key areas across the network. For distributed dynamic phenomena, such as generalized seizures and slow-wave sleep, studying effective connectivity from real-time recordings is significantly complicated since (i) outputs from only a subnetwork can be practically measured, and (ii) exogenous subnetwork inputs are unobservable. Model fitting, therefore, constitutes a challenging blind module identification or model inversion problem for finding both the parameters and the many unknown inputs of the subnetwork. We herein propose a novel estimation framework for identifying nonlinear dynamic subnetworks in the case of slowly-varying, otherwise unknown local inputs. Starting with approximate predictions obtained using Cubature Kalman filtering, residuals of local output predictions are utilized to improve upon local input estimates. The algorithm performance is tested on both simulated and clinical EEG of induced seizures under electroconvulsive therapy (ECT). For the simulated network, the algorithm significantly boosted the estimation accuracy for inputs and connections from noisy EEG. For the clinical data, the algorithm predicted increased subnetwork inputs during the pre-stimulus anesthesia condition. Importantly, it predicted an increased frontocentral connectivity during the generalized seizure that is commensurate with electrode placement and that corroborates the clinical hypothesis of increased frontal focality of therapeutic ECT seizures. The proposed framework can be extended to account for several input configurations and can in principle be applied to study effective connectivity within brain subnetworks defined at the microscale (cortical lamina interaction) or at the macroscale (sensory integration).

Adaptive optimal input design and parametric estimation of nonlinear dynamical systems: application to neuronal modeling

OBJECTIVE

Many physical models of biological processes including neural systems are characterized by parametric nonlinear dynamical relations between driving inputs, internal states, and measured outputs of the process. Fitting such models using experimental data (data assimilation) is a challenging task since the physical process often operates in a noisy, possibly non-stationary environment; moreover, conducting multiple experiments under controlled and repeatable conditions can be impractical, time consuming or costly. The accuracy of model identification, therefore, is dictated principally by the quality and dynamic richness of collected data over single or few experimental sessions. Accordingly, it is highly desirable to design efficient experiments that, by exciting the physical process with smart inputs, yields fast convergence and increased accuracy of the model.

APPROACH

We herein introduce an adaptive framework in which optimal input design is integrated with square root cubature Kalman filters (OID-SCKF) to develop an online estimation procedure that first, converges significantly quicker, thereby permitting model fitting over shorter time windows, and second, enhances model accuracy when only few process outputs are accessible. The methodology is demonstrated on common nonlinear models and on a four-area neural mass model with noisy and limited measurements. Estimation quality (speed and accuracy) is benchmarked against high-performance SCKF-based methods that commonly employ dynamically rich informed inputs for accurate model identification.

MAIN RESULTS

For all the tested models, simulated single-trial and ensemble averages showed that OID-SCKF exhibited (i) faster convergence of parameter estimates and (ii) lower dependence on inter-trial noise variability with gains up to around 1000 ms in speed and 81% increase in variability for the neural mass models. In terms of accuracy, OID-SCKF estimation was superior, and exhibited considerably less variability across experiments, in identifying model parameters of (a) systems with challenging model inversion dynamics and (b) systems with fewer measurable outputs that directly relate to the underlying processes.

SIGNIFICANCE

Fast and accurate identification therefore carries particular promise for modeling of transient (short-lived) neuronal network dynamics using a spatially under-sampled set of noisy measurements, as is commonly encountered in neural engineering applications.

Non-invasive laminar inference with MEG: Comparison of methods and source inversion algorithms

Magnetoencephalography (MEG) is a direct measure of neuronal current flow; its anatomical resolution is therefore not constrained by physiology but rather by data quality and the models used to explain these data. Recent simulation work has shown that it is possible to distinguish between signals arising in the deep and superficial cortical laminae given accurate knowledge of these surfaces with respect to the MEG sensors. This previous work has focused around a single inversion scheme (multiple sparse priors) and a single global parametric fit metric (free energy). In this paper we use several different source inversion algorithms and both local and global, as well as parametric and non-parametric fit metrics in order to demonstrate the robustness of the discrimination between layers. We find that only algorithms with some sparsity constraint can successfully be used to make laminar discrimination. Importantly, local t-statistics, global cross-validation and free energy all provide robust and mutually corroborating metrics of fit. We show that discrimination accuracy is affected by patch size estimates, cortical surface features, and lead field strength, which suggests several possible future improvements to this technique. This study demonstrates the possibility of determining the laminar origin of MEG sensor activity, and thus directly testing theories of human cognition that involve laminar- and frequency-specific mechanisms. This possibility can now be achieved using recent developments in high precision MEG, most notably the use of subject-specific head-casts, which allow for significant increases in data quality and therefore anatomically precise MEG recordings.

Section

Analysis methods.

Classifications

Source localization: inverse problem; Source localization: other.

•

Laminar inferences can be made with MEG using both local and global fit metrics.

•

Source inversion algorithms with sparsity constraints performed best.

•

Classification is affected by patch size estimates, anatomy, and lead field strength.

The rate of transient beta frequency events predicts behavior across tasks and species

Beta oscillations (15-29Hz) are among the most prominent signatures of brain activity. Beta power is predictive of healthy and abnormal behaviors, including perception, attention and motor action. In non-averaged signals, beta can emerge as transient high-power 'events'. As such, functionally relevant differences in averaged power across time and trials can reflect changes in event number, power, duration, and/or frequency span. We show that functionally relevant differences in averaged beta power in primary somatosensory neocortex reflect a difference in the number of high-power beta events per trial, i.e. event rate. Further, beta events occurring close to the stimulus were more likely to impair perception. These results are consistent across detection and attention tasks in human magnetoencephalography, and in local field potentials from mice performing a detection task. These results imply that an increased propensity of beta events predicts the failure to effectively transmit information through specific neocortical representations.

Beyond the Status Quo: A Role for Beta Oscillations in Endogenous Content (Re)Activation

Among the rhythms of the brain, oscillations in the beta frequency range (∼13-30 Hz) have been considered the most enigmatic. Traditionally associated with sensorimotor functions, beta oscillations have recently become more broadly implicated in top-down processing, long-range communication, and preservation of the current brain state. Here, we extend and refine these views based on accumulating new findings of content-specific beta-synchronization during endogenous information processing in working memory (WM) and decision making. We characterize such content-specific beta activity as short-lived, flexible network dynamics supporting the endogenous (re)activation of cortical representations. Specifically, we suggest that beta-mediated ensemble formation within and between cortical areas may awake, rather than merely preserve, an endogenous cognitive set in the service of current task demands. This proposal accommodates key aspects of content-specific beta modulations in monkeys and humans, integrates with timely computational models, and outlines a functional role for beta that fits its transient temporal characteristics.

Open Ephys electroencephalography (Open Ephys + EEG): a modular, low-cost, open-source solution to human neural recording

OBJECTIVE

Electroencephalography (EEG) offers a unique opportunity to study human neural activity non-invasively with millisecond resolution using minimal equipment in or outside of a lab setting. EEG can be combined with a number of techniques for closed-loop experiments, where external devices are driven by specific neural signals. However, reliable, commercially available EEG systems are expensive, often making them impractical for individual use and research development. Moreover, by design, a majority of these systems cannot be easily altered to the specification needed by the end user. We focused on mitigating these issues by implementing open-source tools to develop a new EEG platform to drive down research costs and promote collaboration and innovation.

APPROACH

Here, we present methods to expand the open-source electrophysiology system, Open Ephys (www.openephys.org), to include human EEG recordings. We describe the equipment and protocol necessary to interface various EEG caps with the Open Ephys acquisition board, and detail methods for processing data. We present applications of Open Ephys  +  EEG as a research tool and discuss how this innovative EEG technology lays a framework for improved closed-loop paradigms and novel brain-computer interface experiments.

MAIN RESULTS

The Open Ephys  +  EEG system can record reliable human EEG data, as well as human EMG data. A side-by-side comparison of eyes closed 8-14 Hz activity between the Open Ephys  +  EEG system and the Brainvision ActiCHamp EEG system showed similar average power and signal to noise.

SIGNIFICANCE

Open Ephys  +  EEG enables users to acquire high-quality human EEG data comparable to that of commercially available systems, while maintaining the price point and extensibility inherent to open-source systems.

Linking canonical microcircuits and neuronal activity: Dynamic causal modelling of laminar recordings

Neural models describe brain activity at different scales, ranging from single cells to whole brain networks. Here, we attempt to reconcile models operating at the microscopic (compartmental) and mesoscopic (neural mass) scales to analyse data from microelectrode recordings of intralaminar neural activity. Although these two classes of models operate at different scales, it is relatively straightforward to create neural mass models of ensemble activity that are equipped with priors obtained after fitting data generated by detailed microscopic models. This provides generative (forward) models of measured neuronal responses that retain construct validity in relation to compartmental models. We illustrate our approach using cross spectral responses obtained from V1 during a visual perception paradigm that involved optogenetic manipulation of the basal forebrain. We find that the resulting neural mass model can distinguish between activity in distinct cortical layers - both with and without optogenetic activation - and that cholinergic input appears to enhance (disinhibit) superficial layer activity relative to deep layers. This is particularly interesting from the perspective of predictive coding, where neuromodulators are thought to boost prediction errors that ascend the cortical hierarchy.

Overview of MEG

Magnetoencephalography (MEG) is a method to study electrical activity in the human brain by recording the neuromagnetic field outside the head. MEG, like electroencephalography (EEG), provides an excellent, millisecond-scale time resolution, and allows the estimation of the spatial distribution of the underlying activity, in favorable cases with a localization accuracy of a few millimeters. To detect the weak neuromagnetic signals, superconducting sensors, magnetically shielded rooms, and advanced signal processing techniques are used. The analysis and interpretation of MEG data typically involves comparisons between subject groups and experimental conditions using various spatial, temporal, and spectral measures of cortical activity and connectivity. The application of MEG to cognitive neuroscience studies is illustrated with studies of spoken language processing in subjects with normal and impaired reading ability. The mapping of spatiotemporal patterns of activity within networks of cortical areas can provide useful information about the functional architecture of the brain related to sensory and cognitive processing, including language, memory, attention, and perception.