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

Beta: bursts of cognition

Beta oscillations are linked to the control of goal-directed processing of sensory information and the timing of motor output. Recent evidence demonstrates they are not sustained but organized into intermittent high-power bursts mediating timely functional inhibition. This implies there is a considerable moment-to-moment variation in the neural dynamics supporting cognition. Beta bursts thus offer new opportunities for studying how sensory inputs are selectively processed, reshaped by inhibitory cognitive operations and ultimately result in motor actions. Recent method advances reveal diversity in beta bursts that provide deeper insights into their function and the underlying neural circuit activity motifs. We propose that brain-wide, spatiotemporal patterns of beta bursting reflect various cognitive operations and that their dynamics reveal nonlinear aspects of cortical processing.

Frontoparietal beta event characteristics are associated with early life stress and psychiatric symptoms in adults

Recent work has found that the presence of transient, oscillatory burst-like events, particularly within the beta band (15-29 Hz), is more closely tied to disease state and behavior across species than traditional electroencephalography (EEG) power metrics. This study sought to examine whether features of beta events over frontoparietal electrodes were associated with early life stress (ELS) and the related clinical presentation. Eighteen adults with documented ELS (n = 18; ELS + ) and eighteen adults without documented ELS (n = 18; ELS-) completed eyes-closed resting state EEG as part of their participation in a larger childhood stress study. The rate, power, duration, and frequency span of transient oscillatory events were calculated within the beta band at five frontoparietal electrodes. ELS variables were positively associated with beta event rate at Fp2 and beta event duration at Pz, in that greater ELS was associated with higher resting rates and longer durations. These beta event characteristics were used to successfully distinguish between ELS + and ELS- groups. In an independent clinical dataset (n = 25), beta event power at Pz was positively correlated with ELS. Beta events deserve ongoing investigation as a potential disease marker of ELS and subsequent psychiatric treatment outcomes.

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.

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.

Human Sensorimotor Beta Event Characteristics and Aperiodic Signal Are Highly Heritable

Individuals' phenotypes, including the brain's structure and function, are largely determined by genes and their interplay. The resting brain generates salient rhythmic patterns that can be characterized noninvasively using functional neuroimaging such as magnetoencephalography (MEG). One of these rhythms, the somatomotor (rolandic) beta rhythm, shows intermittent high amplitude "events" that predict behavior across tasks and species. Beta rhythm is altered in neurological disease. The aperiodic (1/f) signal present in electrophysiological recordings is also modulated by some neurological conditions and aging. Both sensorimotor beta and aperiodic signal could thus serve as biomarkers of sensorimotor function. Knowledge about the extent to which these brain functional measures are heritable could shed light on the mechanisms underlying their generation. We investigated the heritability and variability of human spontaneous sensorimotor beta rhythm events and aperiodic activity in 210 healthy male and female adult siblings' spontaneous MEG activity. The most heritable trait was the aperiodic 1/f signal, with a heritability of 0.87 in the right hemisphere. Time-resolved beta event amplitude parameters were also highly heritable, whereas the heritabilities for overall beta power, peak frequency, and measures of event duration remained nonsignificant. Human sensorimotor neural activity can thus be dissected into different components with variable heritability. We postulate that these differences partially reflect different underlying signal-generating mechanisms. The 1/f signal and beta event amplitude measures may depend more on fixed, anatomical parameters, whereas beta event duration and its modulation reflect dynamic characteristics, guiding their use as potential disease biomarkers.

HNN-core: A Python software for cellular and circuit-level interpretation of human MEG/EEG

HNN-core is a library for circuit and cellular level interpretation of non-invasive human magneto-/electro-encephalography (MEG/EEG) data. It is based on the Human Neocortical Neurosolver (HNN) software (Neymotin et al., 2020), a modeling tool designed to simulate multiscale neural mechanisms generating current dipoles in a localized patch of neocortex. HNN's foundation is a biophysically detailed neural network representing a canonical neocortical column containing populations of pyramidal and inhibitory neurons together with layer-specific exogenous synaptic drive (Figure 1 left). In addition to simulating network-level interactions, HNN produces the intracellular currents in the long apical dendrites of pyramidal cells across the cortical layers known to be responsible for macroscopic current dipole generation.

Bursting with Potential: How Sensorimotor Beta Bursts Develop from Infancy to Adulthood

Beta activity is thought to play a critical role in sensorimotor processes. However, little is known about how activity in this frequency band develops. Here, we investigated the developmental trajectory of sensorimotor beta activity from infancy to adulthood. We recorded EEG from 9-month-old, 12-month-old, and adult humans (male and female) while they observed and executed grasping movements. We analyzed "beta burst" activity using a novel method that combines time-frequency decomposition and principal component analysis. We then examined the changes in burst rate and waveform motifs along the selected principal components. Our results reveal systematic changes in beta activity during action execution across development. We found a decrease in beta burst rate during movement execution in all age groups, with the greatest decrease observed in adults. Additionally, we identified three principal components that defined waveform motifs that systematically changed throughout the trial. We found that bursts with waveform shapes closer to the median waveform were not rate-modulated, whereas those with waveform shapes further from the median were differentially rate-modulated. Interestingly, the decrease in the rate of certain burst motifs occurred earlier during movement and was more lateralized in adults than in infants, suggesting that the rate modulation of specific types of beta bursts becomes increasingly refined with age.SIGNIFICANCE STATEMENT We demonstrate that, like in adults, sensorimotor beta activity in infants during reaching and grasping movements occurs in bursts, not oscillations like thought traditionally. Furthermore, different beta waveform shapes were differentially modulated with age, including more lateralization in adults. Aberrant beta activity characterizes various developmental disorders and motor difficulties linked to early brain injury, so looking at burst waveform shape could provide more sensitivity for early identification and treatment of affected individuals before any behavioral symptoms emerge. More generally, comparison of beta burst activity in typical versus atypical motor development will also be instrumental in teasing apart the mechanistic functional roles of different types of beta bursts.

Pre-treatment frontal beta events are associated with executive dysfunction improvement after repetitive transcranial magnetic stimulation for depression: A preliminary report

Repetitive transcranial magnetic stimulation (rTMS) is an established clinical treatment for major depressive disorder (MDD) that has also been found to improve aspects of executive functioning. The objective of this study was to examine whether oscillatory burst-like events within the beta band (15-29 Hz) prior to treatment could predict subsequent change in self-reported executive dysfunction (EDF) across a clinical course of rTMS for MDD. Twenty-eight adults (64% female) with MDD completed the self-report Frontal Systems Behavior Scale (FrSBe) and provided eyes-closed resting-state electroencephalography (EEG) before and after a clinical course of rTMS therapy for primary MDD. The rate, power, duration, and frequency span of transient EEG measured oscillatory beta events were calculated. Events within delta/theta and alpha bands were examined to assess for beta specificity. After controlling for improvement in primary depressive symptoms, a lower rate of beta events at F3, Fz, F4, and Cz prior to rTMS treatment was associated with a larger improvement in EDF after rTMS treatment. In addition, a decrease in beta event rate at Fz pre-to-post treatment was associated with a larger improvement in EDF after treatment. Results were largely specific to the beta band. In this study, the rate of frontrocentral beta events prior to treatment significantly predicted the likelihood of subsequent improvement in EDF symptoms following a clinical course of rTMS for MDD. These preliminary findings suggest the potential utility of EEG measured beta events and rTMS for targeting EDF across an array of neuropsychiatric disorders.

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.

Diverse beta burst waveform motifs characterize movement-related cortical dynamics

Classical analyses of induced, frequency-specific neural activity typically average band-limited power over trials. More recently, it has become widely appreciated that in individual trials, beta band activity occurs as transient bursts rather than amplitude-modulated oscillations. Most studies of beta bursts treat them as unitary, and having a stereotyped waveform. However, we show there is a wide diversity of burst shapes. Using a biophysical model of burst generation, we demonstrate that waveform variability is predicted by variability in the synaptic drives that generate beta bursts. We then use a novel, adaptive burst detection algorithm to identify bursts from human MEG sensor data recorded during a joystick-based reaching task, and apply principal component analysis to burst waveforms to define a set of dimensions, or motifs, that best explain waveform variance. Finally, we show that bursts with a particular range of waveform motifs, ones not fully accounted for by the biophysical model, differentially contribute to movement-related beta dynamics. Sensorimotor beta bursts are therefore not homogeneous events and likely reflect distinct computational processes.

Fronto-central resting-state 15-29 Hz transient beta events change with therapeutic transcranial magnetic stimulation for posttraumatic stress disorder and major depressive disorder

Repetitive transcranial magnetic stimulation (rTMS) is an established treatment for major depressive disorder (MDD) and shows promise for posttraumatic stress disorder (PTSD), yet effectiveness varies. Electroencephalography (EEG) can identify rTMS-associated brain changes. EEG oscillations are often examined using averaging approaches that mask finer time-scale dynamics. Recent advances show some brain oscillations emerge as transient increases in power, a phenomenon termed "Spectral Events," and that event characteristics correspond with cognitive functions. We applied Spectral Event analyses to identify potential EEG biomarkers of effective rTMS treatment. Resting 8-electrode EEG was collected from 23 patients with MDD and PTSD before and after 5 Hz rTMS targeting the left dorsolateral prefrontal cortex. Using an open-source toolbox ( https://github.com/jonescompneurolab/SpectralEvents ), we quantified event features and tested for treatment associated changes. Spectral Events in delta/theta (1-6 Hz), alpha (7-14 Hz), and beta (15-29 Hz) bands occurred in all patients. rTMS-induced improvement in comorbid MDD PTSD were associated with pre- to post-treatment changes in fronto-central electrode beta event features, including frontal beta event frequency spans and durations, and central beta event maxima power. Furthermore, frontal pre-treatment beta event duration correlated negatively with MDD symptom improvement. Beta events may provide new biomarkers of clinical response and advance the understanding of rTMS.

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.

Fronto-central resting-state 15-29Hz transient beta events change with therapeutic transcranial magnetic stimulation for posttraumatic stress disorder and major depressive disorder

Repetitive transcranial magnetic stimulation (rTMS) is an established treatment for major depressive disorder (MDD) and shows promise for posttraumatic stress disorder (PTSD), yet effectiveness varies. Electroencephalography (EEG) can identify rTMS-associated brain changes. EEG oscillations are often examined using averaging approaches that mask finer time-scale dynamics. Recent advances show some brain oscillations emerge as transient increases in power, a phenomenon termed "Spectral Events," and that event characteristics correspond with cognitive functions. We applied Spectral Event analyses to identify potential EEG biomarkers of effective rTMS treatment. Resting 8-electrode EEG was collected from 23 patients with MDD and PTSD before and after 5Hz rTMS targeting the left dorsolateral prefrontal cortex. Using an open-source toolbox ( https://github.com/jonescompneurolab/SpectralEvents ), we quantified event features and tested for treatment associated changes. Spectral Events in delta/theta (1-6 Hz), alpha (7-14 Hz), and beta (15-29 Hz) bands occurred in all patients. rTMS-induced improvement in comorbid MDD PTSD were associated with pre-to post-treatment changes in fronto-central electrode beta event features, including frontal beta event frequency spans and durations, and central beta event maxima power. Furthermore, frontal pre-treatment beta event duration correlated negatively with MDD symptom improvement. Beta events may provide new biomarkers of clinical response and advance the understanding of rTMS.

β-Bursts over Frontal Cortex Track the Surprise of Unexpected Events in Auditory, Visual, and Tactile Modalities

One of the fundamental ways in which the brain regulates and monitors behavior is by making predictions about the sensory environment and adjusting behavior when those expectations are violated. As such, surprise is one of the fundamental computations performed by the human brain. In recent years, it has been well established that one key aspect by which behavior is adjusted during surprise is inhibitory control of the motor system. Moreover, because surprise automatically triggers inhibitory control without much proactive influence, it can provide unique insights into largely reactive control processes. Recent years have seen tremendous interest in burst-like β frequency events in the human (and nonhuman) local field potential-especially over (p)FC-as a potential signature of inhibitory control. To date, β-bursts have only been studied in paradigms involving a substantial amount of proactive control (such as the stop-signal task). Here, we used two cross-modal oddball tasks to investigate whether surprise processing is accompanied by increases in scalp-recorded β-bursts. Indeed, we found that unexpected events in all tested sensory domains (haptic, auditory, visual) were followed by low-latency increases in β-bursting over frontal cortex. Across experiments, β-burst rates were positively correlated with estimates of surprise derived from Shannon's information theory, a type of surprise that represents the degree to which a given stimulus violates prior expectations. As such, the current work clearly implicates frontal β-bursts as a signature of surprise processing. We discuss these findings in the context of common frameworks of inhibitory and cognitive control after unexpected events.

Spectral slowing in chronic stroke reflects abnormalities in both periodic and aperiodic neural dynamics

Decades of electrophysiological work have demonstrated the presence of "spectral slowing" in stroke patients - a prominent shift in the power spectrum towards lower frequencies, most evident in the vicinity of the lesion itself. Despite the reliability of this slowing as a marker of dysfunctional tissue across patient groups as well as animal models, it has yet to be explained in terms of the pathophysiological processes of stroke. To do so requires clear understanding of the neural dynamics that these differences represent, acknowledging the often overlooked fact that spectral power reflects more than just the amplitude of neural oscillations. To accomplish this, we used a combination of frequency domain and time domain measures to disambiguate and quantify periodic (oscillatory) and aperiodic (non-oscillatory) neural dynamics in resting state magnetoencephalography (MEG) recordings from chronic stroke patients. We found that abnormally elevated low frequency power in these patients was best explained by a steepening of the aperiodic component of the power spectrum, rather than an enhancement of low frequency oscillations, as is often assumed. However, genuine oscillatory activity at higher frequencies was also found to be abnormal, with patients showing alpha slowing and diminished oscillatory activity in the beta band. These aperiodic and periodic abnormalities were found to covary, and could be detected even in the un-lesioned hemisphere, however they were most prominent in perilesional tissue, where their magnitude was predictive of cognitive impairment. This work redefines spectral slowing as a pattern of changes involving both aperiodic and periodic neural dynamics and narrows the gap in understanding between non-invasive markers of dysfunctional tissue and disease processes responsible for altered neural dynamics.

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.