Thalamocortical Impulse Propagation and Information Transfer in EEG and MEG

First activity and interactions in thalamus and cortex using raw single-trial EEG and MEG elicited by somatosensory stimulation

Introduction

One of the primary motivations for studying the human brain is to comprehend how external sensory input is processed and ultimately perceived by the brain. A good understanding of these processes can promote the identification of biomarkers for the diagnosis of various neurological disorders; it can also provide ways of evaluating therapeutic techniques. In this work, we seek the minimal requirements for identifying key stages of activity in the brain elicited by median nerve stimulation.

Methods

We have used a priori knowledge and applied a simple, linear, spatial filter on the electroencephalography and magnetoencephalography signals to identify the early responses in the thalamus and cortex evoked by short electrical stimulation of the median nerve at the wrist. The spatial filter is defined first from the average EEG and MEG signals and then refined using consistency selection rules across ST. The refined spatial filter is then applied to extract the timecourses of each ST in each targeted generator. These ST timecourses are studied through clustering to quantify the ST variability. The nature of ST connectivity between thalamic and cortical generators is then studied within each identified cluster using linear and non-linear algorithms with time delays to extract linked and directional activities. A novel combination of linear and non-linear methods provides in addition discrimination of influences as excitatory or inhibitory.

Results

Our method identifies two key aspects of the evoked response. Firstly, the early onset of activity in the thalamus and the somatosensory cortex, known as the P14 and P20 in EEG and the second M20 for MEG. Secondly, good estimates are obtained for the early timecourse of activity from these two areas. The results confirm the existence of variability in ST brain activations and reveal distinct and novel patterns of connectivity in different clusters.

Discussion

It has been demonstrated that we can extract new insights into stimulus processing without the use of computationally costly source reconstruction techniques which require assumptions and detailed modeling of the brain. Our methodology, thanks to its simplicity and minimal computational requirements, has the potential for real-time applications such as in neurofeedback systems and brain-computer interfaces.

Time-Frequency Analysis of Somatosensory Evoked High-Frequency (600 Hz) Oscillations as an Early Indicator of Arousal Recovery after Hypoxic-Ischemic Brain Injury

Cardiac arrest (CA) remains the leading cause of coma, and early arousal recovery indicators are needed to allocate critical care resources properly. High-frequency oscillations (HFOs) of somatosensory evoked potentials (SSEPs) have been shown to indicate responsive wakefulness days following CA. Nonetheless, their potential in the acute recovery phase, where the injury is reversible, has not been tested. We hypothesize that time-frequency (TF) analysis of HFOs can determine arousal recovery in the acute recovery phase. To test our hypothesis, eleven adult male Wistar rats were subjected to asphyxial CA (five with 3-min mild and six with 7-min moderate to severe CA) and SSEPs were recorded for 60 min post-resuscitation. Arousal level was quantified by the neurological deficit scale (NDS) at 4 h. Our results demonstrated that continuous wavelet transform (CWT) of SSEPs localizes HFOs in the TF domain under baseline conditions. The energy dispersed immediately after injury and gradually recovered. We proposed a novel TF-domain measure of HFO: the total power in the normal time-frequency space (NTFS) of HFO. We found that the NTFS power significantly separated the favorable and unfavorable outcome groups. We conclude that the NTFS power of HFOs provides earlier and objective determination of arousal recovery after CA.

Directionality of the injected current targeting the P20/N20 source determines the efficacy of 140 Hz transcranial alternating current stimulation (tACS)-induced aftereffects in the somatosensory cortex

Interindividual anatomical differences in the human cortex can lead to suboptimal current directions and may result in response variability of transcranial electrical stimulation methods. These differences in brain anatomy require individualized electrode stimulation montages to induce an optimal current density in the targeted area of each individual subject. We aimed to explore the possible modulatory effects of 140 Hz transcranial alternating current stimulation (tACS) on the somatosensory cortex using personalized multi-electrode stimulation montages. In two randomized experiments using either tactile finger or median nerve stimulation, we measured by evoked potentials the plasticity aftereffects and oscillatory power changes after 140 Hz tACS at 1.0 mA as compared to sham stimulation (n = 17, male = 9). We found a decrease in the power of oscillatory mu-rhythms during and immediately after tactile discrimination tasks, indicating an engagement of the somatosensory system during stimulus encoding. On a group level both the oscillatory power and the evoked potential amplitudes were not modulated by tACS neither after tactile finger stimulation nor after median nerve stimulation as compared to sham stimulation. On an individual level we could however demonstrate that lower angular difference (i.e., differences between the injected current vector in the target region and the source orientation vector) is associated with significantly higher changes in both P20/N20 and N30/P30 source activities. Our findings suggest that the higher the directionality of the injected current correlates to the dipole orientation the greater the tACS-induced aftereffects are.

Reconstructing subcortical and cortical somatosensory activity via the RAMUS inverse source analysis technique using median nerve SEP data

This study concerns reconstructing brain activity at various depths based on non-invasive EEG (electroencephalography) scalp measurements. We aimed at demonstrating the potential of the RAMUS (randomized multiresolution scanning) technique in localizing weakly distinguishable far-field sources in combination with coinciding cortical activity. As we have shown earlier theoretically and through simulations, RAMUS is a novel mathematical method that by employing the multigrid concept, allows marginalizing noise and depth bias effects and thus enables the recovery of both cortical and subcortical brain activity. To show this capability with experimental data, we examined the 14-30 ms post-stimulus somatosensory evoked potential (SEP) responses of human median nerve stimulation in three healthy adult subjects. We aim at reconstructing the different response components by evaluating a RAMUS-based estimate for the primary current density in the nervous tissue. We present source reconstructions obtained with RAMUS and compare them with the literature knowledge of the SEP components and the outcome of the unit-noise gain beamformer (UGNB) and standardized low-resolution brain electromagnetic tomography (sLORETA). We also analyzed the effect of the iterative alternating sequential technique, the optimization technique of RAMUS, compared to the classical minimum norm estimation (MNE) technique. Matching with our previous numerical studies, the current results suggest that RAMUS could have the potential to enhance the detection of simultaneous deep and cortical components and the distinction between the evoked sulcal and gyral activity.

Individually optimized multi-channel tDCS for targeting somatosensory cortex

OBJECTIVE

Transcranial direct current stimulation (tDCS) is a non-invasive neuro-modulation technique that delivers current through the scalp by a pair of patch electrodes (2-Patch). This study proposes a new multi-channel tDCS (mc-tDCS) optimization method, the distributed constrained maximum intensity (D-CMI) approach. For targeting the P20/N20 somatosensory source at Brodmann area 3b, an integrated combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis is used with individualized skull conductivity calibrated realistic head modeling.

METHODS

Simulated electric fields (EF) for our new D-CMI method and the already known maximum intensity (MI), alternating direction method of multipliers (ADMM) and 2-Patch methods were produced and compared for the individualized P20/N20 somatosensory target for 10 subjects.

RESULTS

D-CMI and MI showed highest intensities parallel to the P20/N20 target compared to ADMM and 2-Patch, with ADMM achieving highest focality. D-CMI showed a slight reduction in intensity compared to MI while reducing side effects and skin level sensations by current distribution over multiple stimulation electrodes.

CONCLUSION

Individualized D-CMI montages are preferred for our follow up somatosensory experiment to provide a good balance between high current intensities at the target and reduced side effects and skin sensations.

SIGNIFICANCE

An integrated combined MEG and EEG source analysis with D-CMI montages for mc-tDCS stimulation potentially can improve control, reproducibility and reduce sensitivity differences between sham and real stimulations.

A novel method for calibrating head models to account for variability in conductivity and its evaluation in a sphere model

The accuracy in electroencephalography (EEG) and combined EEG and magnetoencephalography (MEG) source reconstructions as well as in optimized transcranial electric stimulation (TES) depends on the conductive properties assigned to the head model, and most importantly on individual skull conductivity. In this study, we present an automatic pipeline to calibrate head models with respect to skull conductivity based on the reconstruction of the P20/N20 response using somatosensory evoked potentials and fields. In order to validate in a well-controlled setup without interplay with numerical errors, we evaluate the accuracy of this algorithm in a 4-layer spherical head model using realistic noise levels as well as dipole sources at different eccentricities with strengths and orientations related to somatosensory experiments. Our results show that the reference skull conductivity can be reliably reconstructed for sources resembling the generator of the P20/N20 response. In case of erroneous assumptions on scalp conductivity, the resulting skull conductivity parameter counterbalances this effect, so that EEG source reconstructions using the fitted skull conductivity parameter result in lower errors than when using the standard value. We propose an automatized procedure to calibrate head models which only relies on non-invasive modalities that are available in a standard MEG laboratory, measures under in vivo conditions and in the low frequency range of interest. Calibrated head modeling can improve EEG and combined EEG/MEG source analysis as well as optimized TES.

Parametrizing the Conditionally Gaussian Prior Model for Source Localization with Reference to the P20/N20 Component of Median Nerve SEP/SEF

In this article, we focused on developing the conditionally Gaussian hierarchical Bayesian model (CG-HBM), which forms a superclass of several inversion methods for source localization of brain activity using somatosensory evoked potential (SEP) and field (SEF) measurements. The goal of this proof-of-concept study was to improve the applicability of the CG-HBM as a superclass by proposing a robust approach for the parametrization of focal source scenarios. We aimed at a parametrization that is invariant with respect to altering the noise level and the source space size. The posterior difference between the gamma and inverse gamma hyperprior was minimized by optimizing the shape parameter, while a suitable range for the scale parameter can be obtained via the prior-over-measurement signal-to-noise ratio, which we introduce as a new concept in this study. In the source localization experiments, the primary generator of the P20/N20 component was detected in the Brodmann area 3b using the CG-HBM approach and a parameter range derived from the existing knowledge of the Tikhonov-regularized minimum norm estimate, i.e., the classical Gaussian prior model. Moreover, it seems that the detection of deep thalamic activity simultaneously with the P20/N20 component with the gamma hyperprior can be enhanced while using a close-to-optimal shape parameter value.

Electrodiagnostic applications of somatosensory evoked high-frequency EEG oscillations: Technical considerations

INTRODUCTION

High frequency oscillations (HFOs) embedded within the somatosensory evoked potential (SEP) are not routinely recorded/measured as part of standard clinical SEPs. However, HFOs could provide important additional diagnostic/prognostic information in various patient groups in whom SEPs are tested routinely. One area is the management of patients with hypoxic ischaemic encephalopathy (HIE) in the intensive care unit (ICU). However, the sensitivity of standard clinical SEP recording techniques for detecting HFOs is unknown.

METHODS

SEPs were recorded using routine clinical methods in 17 healthy subjects (median nerve stimulation; 0.5 ms pulse width; 5 Hz; maximum 4000 stimuli) in an unshielded laboratory. Bipolar EEG recordings were acquired (gain 50 k; bandpass 3Hz-2 kHz; sampling rate 5 kHz; non-inverting electrode 2 cm anterior to C3/C4; inverting electrode 2 cm posterior to C3/C4). Data analysis was performed in MATLAB.

RESULTS

SEP-HFOs were detected in 65% of controls using standard clinical recording techniques. In 3 controls without significant HFOs, experiments were repeated using a linear electrode array with higher spatial sampling frequency. SEP-HFOs were observed in all 3 subjects.

CONCLUSIONS

Currently standard clinical methods of recording SEPs are not sufficiently sensitive to permit the inclusion of SEP-HFOs in routine clinical diagnostic/prognostic assessments. Whilst an increase in the number/density of EEG electrodes should improve the sensitivity for detecting SEP-HFOs, this requires confirmation. By improving and standardising clinical SEP recording protocols to permit the acquisition/analysis of SEP-HFOs, it should be possible to gain important insights into the pathophysiology of neurological disorders and refine the management of conditions such as HIE.

Cortical somatosensory evoked high-frequency (600Hz) oscillations predict absence of severe hypoxic encephalopathy after resuscitation

OBJECTIVE

Following cardiac arrest (CA), hypoxic encephalopathy (HE) frequently occurs and hence reliable neuroprognostication is crucial to decide on the extent of intensive care. Several investigations predict severe HE leading to persistent unresponsive wakefulness or death, with high specificity. Only few studies attempted to predict absence of severe HE. Cortical somatosensory evoked high-frequency (600Hz) oscillation (HFO) bursts indicate the presence of highly synchronized spiking activity in the primary somatosensory cortex. Since global neuronal damage characterizes severe HE preserved cortical HFOs may early exclude severe HE.

METHODS

We determined amplitudes of early and late HFO bursts in 302 comatose CA patients after median nerve somatosensory evoked potential (SSEPs) and clinical outcome upon intensive care unit discharge using the cerebral performance category (CPC) scale.

RESULTS

We detected significant early HFO bursts in 146 patients and late HFO bursts in 95 patients. Only one of 27 unresponsive wakefulness patients had a late HFO burst amplitude above 70nV and all seventeen patients who died despite higher amplitudes died from non-neurological causes.

CONCLUSIONS

High-frequency SSEP components can reliably be studied in comatose CA patients using standard equipment.

SIGNIFICANCE

Late HFO burst amplitudes above 70nV largely exclude severe HE incompatible with regaining consciousness.