Search

The amplitude envelope of speech carries crucial low-frequency acoustic information that assists linguistic decoding at multiple time scales. Neurophysiological signals are known to track the amplitude envelope of adult-directed speech (ADS), particularly in the theta-band. Acoustic analysis of infant-directed speech (IDS) has revealed significantly greater modulation energy than ADS in an amplitude-modulation (AM) band centered on ∼2 Hz. Accordingly, cortical tracking of IDS by delta-band neural signals may be key to language acquisition. Speech also contains acoustic information within its higher-frequency bands (beta, gamma). Adult EEG and MEG studies reveal an oscillatory hierarchy, whereby low-frequency (delta, theta) neural phase dynamics temporally organize the amplitude of high-frequency signals (phase amplitude coupling, PAC). Whilst consensus is growing around the role of PAC in the matured adult brain, its role in the development of speech processing is unexplored.Here, we examined the presence and maturation of low-frequency (Graphical abstractHighlightsLongitudinal EEG study in which 4, 7- & 11-month infants listened to nursery rhymesWe demonstrate cortical speech tracking via delta & theta neural signals (mTRF)Periodogram (PSD) analysis revealed stimulus related delta & theta PSD peaksDelta and theta driven phase amplitude coupling (PAC) was found at all agesGamma frequency amplitudes displayed stronger PAC to low frequency phases than beta

Interactions between different brain regions can be revealed by dependencies between their neuronal oscillations. We examined the sensitivity of different oscillatory connectivity measures in revealing interhemispheric interactions between primary motor cortices (M1s) during unilateral finger movements. Based on frequency, amplitude, and phase of the oscillations, a number of metrics have been developed to measure connectivity between brain regions, and each metric has its own strengths, weaknesses, and pitfalls. Taking advantage of the well-known movement-related modulations of oscillatory amplitude in M1s, this study compared and contrasted a number of leading connectivity metrics during distinct phases of oscillatory power changes. Between M1s during unilateral movements, we found that phase-based metrics were effective at revealing connectivity during the beta (15–35 Hz) rebound period linked to movement termination, but not during the early period of beta desynchronization occurring during the movement itself. Amplitude correlation metrics revealed robust connectivity during both periods. Techniques for estimating the direction of connectivity had limited success. Granger Causality was not well suited to studying these connections because it was strongly confounded by differences in signal-to-noise ratio linked to modulation of beta amplitude occurring during the task. Phase slope index was suggestive but not conclusive of a unidirectional influence between motor cortices during the beta rebound. Our findings suggest that a combination of amplitude and phase-based metrics is likely required to fully characterize connectivity during task protocols that involve modulation of oscillatory power, and that amplitude-based metrics appear to be more sensitive despite the lack of directional information.

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

The neurophysiological response to visual stimulation in both humans and animals is characterized by an increase in high frequency amplitude peaking in the gamma range (40–100Hz) and a suppression of low frequency amplitude peaking in the alpha range (10–16Hz). Due to the large number of studies linking amplitude and peak frequency to perception and neurological disorders, there is great interest in understanding the basis of inter-subject variability in gamma and alpha responses. To address this, we measured gamma and alpha amplitude and peak frequency of response to visual stimulation in 42 healthy humans. Using FMRI to delineate active cortical tissue in the same subjects, we correlated these neurophysiological metrics with two structural metrics: distance from active cortex to electrode, and dipole cancellation over active cortex. We find that distance strongly predicted inter-subject gamma amplitude, but had little effect on alpha amplitude, while cancellation had little effect on gamma or alpha amplitude. Neither alpha peak frequency nor gamma peak frequency correlated with our structural metrics. These results suggest that inter-subject variability in gamma amplitude may reflect gross morphology rather than neurophysiological variability, and should be interpreted with caution, while peak frequency may serve as a more sensitive metric of differences in neuronal activity across subjects.

Electroencephalographic (EEG) and magnetoencephalographic (MEG) signals can often be exposed to strong power line interference at 50 or 60 Hz. A widely used method to remove line noise is the notch filter, but it comes with the risk of potentially severe signal distortions. Among other approaches, the Discrete Fourier Transform (DFT) filter and CleanLine have been developed as alternatives, but they may fail to remove power line noise of highly fluctuating amplitude. Here we introduce spectrum interpolation as a new method to remove line noise in the EEG and MEG signal. This approach had been developed for electromyographic (EMG) signals, and combines the advantages of a notch filter, while synthetic test signals indicate that it introduces less distortion in the time domain. The effectiveness of this method is compared to CleanLine, the notch (Butterworth) and DFT filter. In order to quantify the performance of these three methods, we used synthetic test signals and simulated power line noise with fluctuating amplitude and abrupt on- and offsets that were added to an MEG dataset free of line noise. In addition, all methods were applied to EEG data with massive power line noise due to acquisition in unshielded settings. We show that spectrum interpolation outperforms the DFT filter and CleanLine, when power line noise is nonstationary. At the same time, spectrum interpolation performs equally well as the notch filter in removing line noise artifacts, but shows less distortions in the time domain in many common situations.

Power (amplitude) and frequency are two important characteristics of EEG alpha oscillations (8–12 Hz). There is an extensive literature showing that alpha power can be modulated in a goal-oriented manner to either enhance or suppress sensory information processing. Only a few studies to date have examined the task-dependent modulation of alpha frequency. Instead, alpha frequency is often viewed as a trait variable, and used to characterize individual differences in cognitive functioning. We performed two experiments to examine the task-dependent modulation of alpha frequency and its functional significance. In the first experiment, high-density EEG was recorded from 21 participants performing a Sternberg working memory task. The results showed that: (1) during memory encoding, alpha frequency decreased with increasing memory load, whereas during memory retention and retrieval, alpha frequency increased with increasing memory load, (2) higher alpha frequency prior to the onset of probe was associated with longer reaction time, and (3) higher alpha frequency prior to the onset of cue or probe was associated with weaker early cue-evoked or probe-evoked neural responses. In the second experiment, simultaneous EEG-fMRI was recorded from 59 participants during resting state. An EEG-informed fMRI analysis revealed that the spontaneous fluctuations of alpha frequency, but not alpha power, were inversely associated with BOLD activity in the visual cortex. Taken together, these findings suggest that alpha frequency is task-dependent, may serve as an indicator of cortical excitability, and along with alpha power, provides more comprehensive indexing of sensory gating.

Typically, time-frequency analysis (TFA) of electrophysiological data is aimed at isolating narrowband signals (oscillatory activity) from broadband non-oscillatory (1/f) activity, so that changes in oscillatory activity resulting from experimental manipulations can be assessed. A widely used method to do this is to convert the data to the decibel (dB) scale through baseline division and log transformation. This procedure assumes that, for each frequency, sources of power (i.e., oscillations and 1/f activity) scale by the same factor relative to the baseline (multiplicative model). This assumption may be incorrect when signal and noise are independent contributors to the power spectrum (additive model). Using resting-state EEG data from 80 participants, we found that the level of 1/f activity and alpha power are not positively correlated within participants, in line with the additive but not the multiplicative model. Then, to assess the effects of dB conversion on data that violate the multiplicativity assumption, we simulated a mixed design study with one between-subject (noise level, i.e., level of 1/f activity) and one within-subject (signal amplitude, i.e., amplitude of oscillatory activity added onto the background 1/f activity) factor. The effect size of the noise level × signal amplitude interaction was examined as a function of noise difference between groups, following dB conversion. Findings revealed that dB conversion led to the over- or under-estimation of the true interaction effect when groups differing in 1/f levels were compared, and it also led to the emergence of illusory interactions when none were present. This is because signal amplitude was systematically underestimated in the noisier compared to the less noisy group. Hence, we recommend testing whether the level of 1/f activity differs across groups or conditions and using multiple baseline correction strategies to validate results if it does. Such a situation may be particularly common in aging, developmental, or clinical studies.

The predominant activity of slow wave sleep is cortical slow oscillations (SOs), thalamic spindles and hippocampal sharp wave ripples. While the precise temporal nesting of these rhythms was shown to be essential for memory consolidation, the coordination mechanism is poorly understood. Here we develop a minimal hippocampo-cortico-thalamic network that can explain the mechanism underlying the SO-spindle-ripple coupling indicating of the succession of regional neuronal interactions. Further we verify the model predictions experimentally in naturally sleeping rodents showing our simple model provides a quantitative match to several experimental observations including the nesting of ripples in the spindle troughs and larger duration but lower amplitude of the ripples co-occurring with spindles or SOs compared to the isolated ripples. The model also predicts that the coupling of ripples to SOs and spindles monotonically enhances by increasing the strength of hippocampo-cortical connections while it is stronger at intermediate values of the cortico-hippocampal projections.

Highlights • Comprehensive investigation of local PAC in human resting state MEG. • Pipeline to distinguish genuine non-harmonic PAC from harmonic PAC. • Observation of wide-spread harmonic local PAC, but no evidence for non-harmonic local PAC. • Spurious PAC due to physiological artifacts even after application of state-of-the-art artifact rejection methods., Phase-amplitude coupling (PAC) has been hypothesized to coordinate cross-frequency interactions of neuronal activity in the brain. However, little is known about the distribution of PAC across the human brain and the frequencies involved. Furthermore, it remains unclear to what extent PAC may reflect spurious cross-frequency coupling induced by physiological artifacts or rhythmic non-sinusoidal signals with higher harmonics. Here, we combined MEG, source-reconstruction and different measures of cross-frequency coupling to systematically characterize local PAC across the resting human brain. We show that cross-frequency measures of phase-amplitude, phase-phase, and amplitude-amplitude coupling are all sensitive to signals with higher harmonics. In conjunction, these measures allow to distinguish harmonic and non-harmonic PAC. Based on these insights, we found no evidence for non-harmonic local PAC in resting-state MEG. Instead, we found cortically and spectrally wide-spread PAC driven by harmonic signals. Furthermore, we show how physiological artifacts and spectral leakage cause spurious PAC across wide frequency ranges. Our results clarify how different measures of cross-frequency interactions can be combined to characterize PAC, and cast doubt on the presence of prominent non-harmonic phase-amplitude coupling in human resting-state MEG.

In EEG data acquired in the presence of fMRI, gradient-related spike artifacts contaminate the signal following the common preprocessing step of average artifact subtraction. Spike artifacts compromise EEG data quality since they overlap with the EEG signal in frequency, thereby confounding frequency-based inferences on activity. As well, spike artifacts can inflate or deflate correlations among time series, thereby confounding inferences on functional connectivity. We present Schrödinger filtering, which uses the Schrödinger equation to decompose the spike-containing input. The basis functions of the decomposition are localized and pulse-shaped, and selectively capture the various input peaks, with the spike components clustered at the beginning of the spectrum. Schrödinger filtering automatically subtracts the spike components from the data. On real and simulated data, we show that Schrödinger filtering (1) simultaneously accomplishes high spike removal and high signal preservation without affecting evoked activity, and (2) reduces spurious pairwise correlations in spontaneous activity. In these regards, Schrödinger filtering was significantly better than three other despiking techniques: median filtering, amplitude thresholding, and wavelet denoising. These results encourage the use of Schrödinger filtering in future EEG-fMRI pipelines, as well as in other spike-related applications (e.g., fMRI motion artifact removal or action potential extraction).

Motor imagery modulates specific neural oscillations like actual movement does. Representatively, suppression of the alpha power (e.g., event-related desynchronization [ERD]) is the typical pattern of motor imagery in the motor cortex. However, in addition to this amplitude-based feature, the coupling across frequencies includes important information about the brain functions and the existence of such complex information has been reported in various invasive studies. Yet, the interaction across multiple frequencies during motor imagery processing is still unclear and has not been widely studied, particularly concerning the non-invasive signals. In this study, we provide empirical evidence of the comodulation between the phase of alpha rhythm and the amplitude of high gamma rhythm during the motor imagery process. We used electroencephalography (EEG) in our investigation during the imagination of left- or right-hand movement recorded from 52 healthy subjects, and quantified the ERD of alpha and phase-amplitude coupling (PAC) which is a relative change of modulation index to the base line period (before the cue). As a result, we found that the coupling between the phase of alpha (8–12 Hz) and the amplitude of high gamma (70–120 Hz) and this PAC decreases during motor imagery and then rebounds to the baseline like alpha ERD (r = 0.29 to 0.42). This correlation between PAC and ERD was particularly stronger in the ipsilateral area. In addition, trials that demonstrated higher alpha power during the ready period (before the cue) showed a larger ERD during motor imagery and similarly, trials with higher modulation index during the ready period yielded a greater decrease in PAC during imagery. In the classification analysis, we found that the effective phase frequency that showed better decoding accuracy in left and right-hand imagery, varied across subjects. Motivated by result, we proposed a weighted cross-frequency coupling (WCFC) method that extracts the maximal discriminative feature by combining band power and CFC. In the evaluation, WCFC with only two electrodes yielded a performance comparable to the conventional algorithm with 64 electrodes in classifying left and right-hand motor imagery. These results indicate that the phase-amplitude frequency plays an important role in motor imagery, and that optimizing this frequency ranges is crucial for extracting information features to decode the motor imagery types.

Functional connectivity (FC), thought to provide a window into neural communication, has become a core focus in the study of brain function and cognition. However, there is no consensus on how to conceptualize large-scale FC in electrophysiology. Phase coupling (PhC), defined as coupling between the phases of two signals, reflects the synchronization of rhythmic oscillation cycles. Conversely, amplitude coupling (AmpC), defined as coupling between the envelopes of two signals, reflects correlation of activation amplitude. Despite quantifying different electrophysiological properties, the relationship between PhC and AmpC remains largely unknown. We assessed spatial and temporal correspondence between PhC and AmpC over 5 canonical frequency bands during a cue-based motor task using electrocorticography (ECoG) in 18 patients (8 females) undergoing presurgical monitoring. Significant correspondence between the spatial pattern of PhC and AmpC was detected during stimulus processing across all subjects and frequency bands (R ​≈ ​0.50 for theta, decreasing with increasing frequency). The cross-measure spatial correlation vanished almost entirely when accounting for the portion of FC equally present during pre- and post-stimulus intervals, suggesting that the spatial correlations reflect intrinsic FC independent of stimulus processing. Stimulus-related processing modulated both PhC and AmpC, however in a spatially independent manner. Examining the linear temporal correlation, we found no evidence for linear dependence between PhC and AmpC. Supporting the robustness of our findings, results extended to a verb generation task in a second ECoG dataset of 6 subjects. We conclude that PhC and AmpC reflect intrinsic FC similarly across space, but exhibit divergent stimulus-related FC changes over space and time.

Over the recent years, significant advances in Spin-Echo (SE) Echo-Planar (EP) Diffusion MRI (dMRI) have enabled improved fiber tracking conspicuity in the human brain. At the same time, pushing the spatial resolution and using higher b-values inherently expose the acquired images to further eddy-current-induced distortion and blurring. Recently developed data-driven correction techniques, capable of significantly mitigating these defects, are included in the reconstruction pipelines developed for the Human Connectome Project (HCP) driven by the NIH BRAIN initiative. In this case, however, corrections are derived from the original diffusion-weighted (DW) magnitude images affected by distortion and blurring. Considering the complexity of k-space deviations in the presence of time varying high spatial order eddy currents, distortion and blurring may not be fully reversed when relying on magnitude DW images only. An alternative approach, consisting of iteratively reconstructing DW images based on the actual magnetic field spatiotemporal evolution measured with a magnetic field monitoring camera, has been successfully implemented at 3T in single band dMRI (Wilm et al., 2017, 2015). In this study, we aim to demonstrate the efficacy of this eddy current correction method in the challenging context of HCP-style multiband (MB ​= ​2) dMRI protocol. The magnetic field evolution was measured during the EP-dMRI readout echo train with a field monitoring camera equipped with 16 19F NMR probes. The time variation of 0th, 1st and 2nd order spherical field harmonics were used to reconstruct DW images. Individual DW images reconstructed with and without field correction were compared. The impact of eddy current correction was evaluated by comparing the corresponding direction-averaged DW images and fractional anisotropy (FA) maps. 19F field monitoring data confirmed the existence of significant field deviations induced by the diffusion-encoding gradients, with variations depending on diffusion gradient amplitude and direction. In DW images reconstructed with the field correction, residual aliasing artifacts were reduced or eliminated, and when high b-values were applied, better gray/white matter delineation and sharper gyri contours were observed, indicating reduced signal blurring. The improvement in image quality further contributed to sharper contours and better gray/white matter delineation in mean DW images and FA maps. In conclusion, we demonstrate that up-to-2nd-order-eddy-current-induced field perturbation in multiband, in-plane accelerated HCP-style dMRI acquisition at 7T can be corrected by integrating the measured field evolution in image reconstruction.

Neural oscillations dominate electrophysiological measures of macroscopic brain activity and fluctuations in these rhythms offer an insightful window on cortical excitation, inhibition, and connectivity. However, in recent years the ‘classical’ picture of smoothly varying oscillations has been challenged by the idea that many ‘oscillations’ may actually be formed from the recurrence of punctate high-amplitude bursts in activity, whose spectral composition intersects the traditionally defined frequency ranges (e.g. alpha/beta band). This finding offers a new interpretation of measurable brain activity, however neither the methodological means to detect bursts, nor their link to other findings (e.g. connectivity) have been settled. Here, we use a new approach to detect bursts in magnetoencephalography (MEG) data. We show that a time-delay embedded Hidden Markov Model (HMM) can be used to delineate single-region bursts which are in agreement with existing techniques. However, unlike existing techniques, the HMM looks for specific spectral patterns in timecourse data. We characterise the distribution of burst duration, frequency of occurrence and amplitude across the cortex in resting state MEG data. During a motor task we show how the movement related beta decrease and post movement beta rebound are driven by changes in burst occurrence. Finally, we show that the beta band functional connectome can be derived using a simple measure of burst overlap, and that coincident bursts in separate regions correspond to a period of heightened coherence. In summary, this paper offers a new methodology for burst identification and connectivity analysis which will be important for future investigations of neural oscillations.

Non-invasive brain stimulation to target specific network activity patterns, e.g. transcranial alternating current stimulation (tACS), has become an essential tool to understand the causal role of neuronal oscillations in cognition and behavior. However, conventional sinusoidal tACS limits the ability to record neuronal activity during stimulation and lacks spatial focality. One particularly promising new tACS stimulation paradigm uses amplitude-modulated (AM) high-frequency waveforms (AM-tACS) with a slow signal envelope that may overcome the limitations. Moreover. AM-tACS using high-frequency carrier signals is more tolerable than conventional tACS, e.g. in terms of skin irritation and occurrence of phosphenes, when applied at the same current intensities (e.g. 1–2 mA). Yet, the fundamental mechanism of neuronal target-engagement by AM-tACS waveforms has remained unknown. We used a computational model of cortex to investigate how AM-tACS modulates endogenous oscillations and compared the target engagement mechanism to the case of conventional (unmodulated) low-frequency tACS. Analysis of stimulation amplitude and frequency indicated that cortical oscillations were phase-locked to the envelope of the AM stimulation signal, which thus exhibits the same target engagement mechanism as conventional (unmodulated) low frequency tACS. However, in the computational model substantially higher current intensities were needed for AM-tACS than for low-frequency (unmodulated) tACS waveforms to achieve pronounced phase synchronization. Our analysis of the carrier frequency suggests that there might be a trade-off between the use of high-frequency carriers and the stimulation amplitude required for successful entrainment. Together, our computational simulations support the use of slow-envelope high frequency carrier AM waveforms as a tool for noninvasive modulation of brain oscillations. More empirical data will be needed to identify the optimal stimulation parameters and to evaluate tolerability and safety of both, AM- and conventional tACS.

When combined with source modeling, magneto‐ (MEG) and electroencephalography (EEG) can be used to study long-range interactions among cortical processes non-invasively. Estimation of such inter-areal connectivity is nevertheless hindered by instantaneous field spread and volume conduction, which artificially introduce linear correlations and impair source separability in cortical current estimates. To overcome the inflating effects of linear source mixing inherent to standard interaction measures, alternative phase‐ and amplitude-correlation based connectivity measures, such as imaginary coherence and orthogonalized amplitude correlation have been proposed. Being by definition insensitive to zero-lag correlations, these techniques have become increasingly popular in the identification of correlations that cannot be attributed to field spread or volume conduction. We show here, however, that while these measures are immune to the direct effects of linear mixing, they may still reveal large numbers of spurious false positive connections through field spread in the vicinity of true interactions. This fundamental problem affects both region-of-interest-based analyses and all-to-all connectome mappings. Most importantly, beyond defining and illustrating the problem of spurious, or “ghost” interactions, we provide a rigorous quantification of this effect through extensive simulations. Additionally, we further show that signal mixing also significantly limits the separability of neuronal phase and amplitude correlations. We conclude that spurious correlations must be carefully considered in connectivity analyses in MEG/EEG source space even when using measures that are immune to zero-lag correlations.Highlights✓Reliable estimation of neuronal coupling with MEG and EEG is challenged by signal mixing✓A number of coupling techniques attempt to overcome this limitation by excluding zero-lag interactions✓Contrary to what is commonly admitted, our simulations illustrate that such interaction metrics will still yield false positives✓Spurious, or “ghost”, interactions are generally detected between sources in the vicinity of true phase-lagged interacting sources✓Signal mixing also severely affects the mutual separability of phase and amplitude correlations

Natural environments elicit both phase-locked and non-phase-locked neural responses to the stimulus in the brain. The interpretation of the BOLD signal to date has been based on an association of the non-phase-locked power of high-frequency local field potentials (LFPs), or the related spiking activity in single neurons or groups of neurons. Previous studies have not examined the prediction of the BOLD signal by phase-locked responses. We examined the relationship between the BOLD response and LFPs in the same nine human subjects from multiple corresponding points in the auditory cortex, using amplitude modulated pure tone stimuli of a duration to allow an analysis of phase locking of the sustained time period without contamination from the onset response. The results demonstrate that both phase locking at the modulation frequency and its harmonics, and the oscillatory power in gamma/high-gamma bands are required to predict the BOLD response. Biophysical models of BOLD signal generation in auditory cortex therefore require revision and the incorporation of both phase locking to rhythmic sensory stimuli and power changes in the ensemble neural activity.

A method is presented for correcting the effects of stimulated and indirect echoes on quantitative T2 mapping data acquired with multiple spin echo techniques, such as turbo spin echo. In contrast to similar correction techniques proposed in the literature, the method does not require a priori knowledge of the radio frequency (RF) pulse profiles. In a first step, for the T2 mapping protocol under investigation, signal decay curves S(TE) are simulated for a range of different RF pulse profiles. The actual signal decay S(TE) is then measured on a phantom with known T2, so the approximate RF pulse profiles can be derived via comparison with the simulated decay curves. In a second step, with the RF pulses obtained from step one, signal decay curves S(TE) are simulated for different T2 values and fitted mono-exponentially, thus allowing to deduce the relationship between true T2 and the apparent T2 (T2app) values. Results show that this relationship is approximately linear, allowing for a direct correction of T2app maps. If the amplitude of the transmitted RF field (B1) does not exceed the nominal value by more than 10%, it is shown that a B1-independent correction of T2app maps yields sufficiently accurate results for T2. A B1-dependent version is also presented. The method is tested in vitro on a phantom with different T2 values and in vivo on healthy subjects.

Natural scenes share a consistent distribution of energy across spatial frequencies (SF) known as the 1/fα amplitude spectrum (α≈0.8–1.5, mean 1.2). This distribution is scale-invariant, which is a fractal characteristic of natural scenes with statistically similar structure at different spatial scales. While the sensitivity of the visual system to the 1/f properties of natural scenes has been studied extensively using psychophysics, relatively little is known about the tuning of cortical responses to these properties. Here, we use fMRI and retinotopic mapping techniques to measure and analyze BOLD responses in early visual cortex (V1, V2, and V3) to synthetic noise images that vary in their 1/fα amplitude spectra (α=0.25 to 2.25, step size: 0.50) and contrast levels (10% and 30%) (Experiment 1). To compare the dependence of the BOLD response between the photometric (intensity based) and geometric (fractal) properties of our stimuli, in Experiment 2 we compared grayscale noise images to their binary (thresholded) counterparts, which contain only black and white regions. In both experiments, early visual cortex responded maximally to stimuli generated to have an input 1/f slope corresponding to natural 1/fα amplitude spectra, and lower BOLD responses were found for steeper or shallower 1/f slopes (peak modulation: 0.59% for 1.25 vs. 0.31% for 2.25). To control for changing receptive field sizes, responses were also analyzed across multiple eccentricity bands in cortical surface space. For most eccentricity bands, BOLD responses were maximal for natural 1/fα amplitude spectra, but importantly there was no difference in the BOLD response to grayscale stimuli and their corresponding thresholded counterparts. Since the thresholding of an image changes its measured 1/f slope (α) but not its fractal characteristics, this suggests that neuronal responses in early visual cortex are not strictly driven by spectral slope values (photometric properties) but rather their embedded geometric, fractal-like scaling properties.

Somatosensory stimulation intensity and behavioral detection are positively related, and both correlate with neural responses. However, it is still controversial as to what extent stimulus intensity and early somatosensory evoked potentials (SEP) predict detection and how these parameters interact with pre-stimulus brain oscillatory states, which also influence sensory processing. Here we investigated how early SEP components encode stimulation intensity, how pre-stimulus alpha- and beta-band amplitudes interact with SEPs, and which neural markers predict stimulus detection. To this end, we randomly presented electrical finger nerve stimulation with various intensities distributed along the individual psychometric response function (including catch trials) while recording the EEG. Participants reported stimulus presence on a trial-by-trial basis (one-alternative-forced-choice). For the lowest (imperceptible) intensities, participants showed zero (behavioral) sensitivity despite measurable early cortical processing reflected by the P50 component. The P50 amplitude scaled with increasing stimulation intensities but was not predictive of stimulus detection. Instead, detection was associated with the later negative N150 component, as well as with pre-stimulus lowered somatosensory alpha- and increased frontal beta-band amplitudes. Our results give evidence for a serial representation of stimulus intensity and detection, as reflected by the P50 and N150 amplitude, respectively. Furthermore, stimulus detection seems to depend on the current brain state, rendering upcoming stimulation being reportable or not.

Spatio-temporal patterns in electroencephalography (EEG) can be described by microstate analysis, a discrete approximation of the continuous electric field patterns produced by the cerebral cortex. Resting-state EEG microstates are largely determined by alpha frequencies (8-12 Hz) and we recently demonstrated that microstates occur periodically with twice the alpha frequency. To understand the origin of microstate periodicity, we analyzed the analytic amplitude and the analytic phase of resting-state alpha oscillations independently. In continuous EEG data we found rotating phase patterns organized around a small number of phase singularities which varied in number and location. The spatial rotation of phase patterns occurred with the underlying alpha frequency. Phase rotors coincided with periodic microstate motifs involving the four canonical microstate maps. The analytic amplitude showed no oscillatory behaviour and was almost static across time intervals of 1-2 alpha cycles, resulting in the global pattern of a standing wave. In n=23 healthy adults, time-lagged mutual information analysis of microstate sequences derived from amplitude and phase signals of awake eyes-closed EEG records showed that only the phase component contributed to the periodicity of microstate sequences. Phase sequences showed mutual information peaks at multiples of 50 ms and the group average had a main peak at 100 ms (10 Hz), whereas amplitude sequences had a slow and monotonous information decay. This result was confirmed by an independent approach combining temporal principal component analysis (tPCA) and autocorrelation analysis. We reproduced our observations in a generic model of EEG oscillations composed of coupled non-linear oscillators (Stuart-Landau model). Phase-amplitude dynamics similar to experimental EEG occurred when the oscillators underwent a supercritical Hopf bifurcation, a common feature of many computational models of the alpha rhythm. These findings explain our previous description of periodic microstate recurrence and its relation to the time scale of alpha oscillations. Moreover, our results corroborate the predictions of computational models and connect experimentally observed EEG patterns to properties of critical oscillator networks.

Item does not contain fulltext Modulations of beta oscillatory power serve a predictive role, in preparation of future actions. It is well known that beta amplitude decreases prior to voluntary movements and expected tactile stimuli. Paradoxically, recent studies have reported a beta amplitude increase prior to expected visual and auditory stimuli. Moreover, it has been suggested that, in isochronic stimulus series, the rising beta slope is adjusted to the duration of the interstimulus interval. We investigated the characteristics of such timing related pre-stimulus beta power increases using visual stimulus sequences that were presented at three regular rates (0.61, 0.74 and 0.95Hz). EEG was recorded from twenty participants while they attended the sequences by performing a clock reading task. Time-frequency analyses showed a consistent pattern of beta modulation: the post-stimulus beta power decrease was followed by a steep increase. Contrary to recent views, we found that the peaks of beta power, for the three presentation rates, were reached at a similar latency post-stimulus, instead of a fixed interval preceding the next stimulus. This demonstrates that, at interstimulus intervals between 1-2s, beta synchronization slopes are not modulated by timing mechanisms related to prediction of upcoming stimuli. We reconcile the discrepant results by proposing that when shorter interval durations are used, as in previous studies, beta resynchronization is interrupted by the presentation of a new stimulus, making it seem as if beta power peaks prior to upcoming stimuli. We emphasize caution with respect to the notion that the timing of beta synchronization is an expression of predictive timing.

Functional imaging investigations into the brain's resting state interactions have yielded a wealth of insight into the intrinsic and dynamic neural architecture supporting cognition and behavior. Electrophysiological studies however have highlighted the fact that synchrony across large-scale cortical systems is composed of spontaneous interactions occurring at timescales beyond the traditional resolution of fMRI, a feature that limits the capacity of fMRI to draw inference on the true directional relationship between network nodes. To approach the question of directionality in resting state signals, we recorded resting state functional MRI (rsfMRI) and electrocorticography (ECoG) from four human subjects undergoing invasive epilepsy monitoring. Using a seed-point based approach, we employed phase-amplitude coupling (PAC) and biPhase Locking Values (bPLV), two measures of cross-frequency coupling (CFC) to explore both outgoing and incoming connections between the seed and all non-seed, site electrodes. We observed robust PAC between a wide range of low-frequency phase and high frequency amplitude estimates. However, significant bPLV, a CFC measure of phase-phase synchrony, was only observed at specific narrow low and high frequency bandwidths. Furthermore, the spatial patterns of outgoing PAC connectivity were most closely associated with the rsfMRI connectivity maps. Our results support the hypothesis that PAC is relatively ubiquitous phenomenon serving as a mechanism for coordinating high-frequency amplitudes across distant neuronal assemblies even in absence of overt task structure. Additionally, we demonstrate that the spatial distribution of a seed-point rsfMRI sensorimotor network is strikingly similar to specific patterns of directional PAC. Specifically, the high frequency activities of distal patches of cortex owning membership in a rsfMRI sensorimotor network were most likely to be entrained to the phase of a low frequency rhythm engendered from the neural populations at the seed-point, suggestive of greater directional coupling from the seed out to the site electrodes.

Alignment of neural oscillations with temporally regular input allows listeners to generate temporal expectations. However, it remains unclear how behavior is governed in the context of temporal variability: What role do temporal expectations play, and how do they interact with the strength of neural oscillatory activity? Here, human participants detected near-threshold targets in temporally variable acoustic sequences. Temporal expectation strength was estimated using an oscillator model and pre-target neural amplitudes in auditory cortex were extracted from magnetoencephalography signals. Temporal expectations modulated target-detection performance, however, only when neural delta-band amplitudes were large. Thus, slow neural oscillations act to gate influences of temporal expectation on perception. Furthermore, slow amplitude fluctuations governed linear and quadratic influences of auditory alpha-band activity on performance. By fusing a model of temporal expectation with neural oscillatory dynamics, the current findings show that human perception in temporally variable contexts relies on complex interactions between multiple neural frequency bands.

© 2015. In resting-state functional connectivity magnetic resonance imaging (fcMRI) studies, measures of functional connectivity are often calculated after the removal of a global mean signal component. While the application of the global signal regression approach has been shown to reduce the influence of physiological artifacts and enhance the detection of functional networks, there is considerable controversy regarding its use as the method can lead to significant bias in the resultant connectivity measures. In addition, evidence from recent studies suggests that the global signal is linked to neural activity and may carry clinically relevant information. For instance, in a prior study we found that the amplitude of the global signal was negatively correlated with EEG measures of vigilance across subjects and experimental runs. Furthermore, caffeine-related decreases in global signal amplitude were associated with increases in EEG vigilance. In this study, we extend the prior work by examining measures of global signal amplitude and EEG vigilance under eyes-closed (EC) and eyes-open (EO) resting-state conditions. We show that changes (EO minus EC) in the global signal amplitude are negatively correlated with the associated changes in EEG vigilance. The slope of this EO-EC relation is comparable with the slope of the previously reported relation between caffeine-related changes in the global signal amplitude and EEG vigilance. Our findings provide further support for a basic relationship between global signal amplitude and EEG vigilance.

Coupling of neuronal oscillations may reflect and facilitate the communication between neuronal populations. Two primary neuronal coupling modes have been described: phase-coupling and amplitude-coupling. Theoretically, both coupling modes are independent, but so far, their neuronal relationship remains unclear. Here, we combined MEG, source-reconstruction and simulations to systematically compare cortical amplitude-coupling and phase-coupling patterns in the human brain. Importantly, we took into account a critical bias of amplitude-coupling measures due to phase-coupling. We found differences between both coupling modes across a broad frequency range and most of the cortex. Furthermore, by combining empirical measurements and simulations we ruled out that these results were caused by methodological biases, but instead reflected genuine neuronal amplitude coupling. Our results show that cortical phase- and amplitude-coupling patterns are non-redundant, which may reflect at least partly distinct neuronal mechanisms. Furthermore, our findings highlight and clarify the compound nature of amplitude coupling measures., Highlights • Systematic comparison of cortical phase- and amplitude-coupling patterns. • Demonstration of genuine amplitude coupling independent of phase coupling bias. • Amplitude- and phase coupling patterns differ across many cortical regions and frequencies.

Mu, beta, and gamma rhythms increase and decrease in amplitude during movement. This event-related synchronization (ERS) and desynchronization (ERD) can be readily recorded non-invasively using magneto- and electro-encephalography (M/EEG). In addition, event-related potentials and fields (i.e., evoked responses) can be elucidated during movement. There is some evidence that the frequency, amplitude and latency of the movement-related ERS/ERD changes with ageing, however the evidence surrounding this topic comes mainly from studies in sample sizes on the order of tens of participants. The objective of this study was to examine a large open-access MEG dataset for age-related changes in movement-related ERS/ERD and evoked responses. MEG data acquired at the Cambridge Centre for Ageing and Neuroscience during cued button pressing was used from 567 participants between the ages of 18 and 88 years. The characteristics movement-related ERD/ERS and evoked responses were calculated for each individual participant. Based on linear regression analysis, significant relationships were found between participant age and some response characteristics, although the predictive value of these relationships was low. Specifically, we conclude that peak beta rebound frequency and amplitude decreased with age, peak beta suppression amplitude increased with age, movement-related gamma burst amplitude decreased with age, and peak motor-evoked response amplitude increased with age. Given our current understanding of the underlying mechanisms of these responses, our findings suggest the existence of age-related changes in the neurophysiology of thalamocortical loops and local circuitry in the primary somatosensory and motor cortices.

Increasing evidence suggests that neuronal communication is a defining property of functionally specialized brain networks and that it is implemented through synchronization between population activities of distinct brain areas. The detection of long-range coupling in electroencephalography (EEG) and magnetoencephalography (MEG) data using conventional metrics (such as coherence or phase-locking value) is by definition contaminated by spatial leakage. Methods such as imaginary coherence, phase-lag index or orthogonalized amplitude correlations tackle spatial leakage by ignoring zero-phase interactions. Although useful, these metrics will by construction lead to false negatives in cases where true zero-phase coupling exists in the data and will underestimate interactions with phase lags in the vicinity of zero. Yet, empirically observed neuronal synchrony in invasive recordings indicates that it is not uncommon to find zero or close-to-zero phase lag between the activity profiles of coupled neuronal assemblies. Here, we introduce a novel method that allows us to mitigate the undesired spatial leakage effects and detect zero and near zero phase interactions. To this end, we propose a projection operation that operates on sensor-space cross-spectrum and suppresses the spatial leakage contribution but retains the true zero-phase interaction component. We then solve the network estimation task as a source estimation problem defined in the product space of interacting source topographies. We show how this framework provides reliable interaction detection for all phase-lag values and we thus refer to the method as Phase Shift Invariant Imaging of Coherent Sources (PSIICOS). Realistic simulations demonstrate that PSIICOS has better detector characteristics than existing interaction metrics. Finally, we illustrate the performance of PSIICOS by applying it to real MEG dataset recorded during a standard mental rotation task. Taken together, using analytical derivations, data simulations and real brain data, this study presents a novel source-space MEG/EEG connectivity method that overcomes previous limitations and for the first time allows for the estimation of true zero-phase coupling via non-invasive electrophysiological recordings.

Traditionally, EEG is understood as originating from the synchronous activation of neuronal populations that generate rhythmic oscillations in specific frequency bands. Recently, new neuronal dynamics regimes have been identified (e.g. neuronal avalanches) characterized by irregular or arrhythmic activity. In addition, it is starting to be acknowledged that broadband properties of EEG spectrum (following a 1 / f law) are tightly linked to brain function. Nevertheless, there is still no theoretical framework accommodating the coexistence of these two EEG phenomenologies: rhythmic/narrowband and arrhythmic/broadband. To address this problem, we present a new framework for EEG analysis based on the relation between the Gaussianity and the envelope of a given signal. EEG Gaussianity is a relevant assessment because if EEG emerges from the superposition of uncorrelated sources, it should exhibit properties of a Gaussian process, otherwise, as in the case of neural synchronization, deviations from Gaussianity should be observed. We use analytical results demonstrating that the coefficient of variation of the envelope (CVE) of Gaussian noise (or any of its filtered sub-bands) is the constant 4 π − 1 ≈ 0.523 , thus enabling CVE to be a useful metric to assess EEG Gaussianity. Furthermore, a new and highly informative analysis space (envelope characterization space) is generated by combining the CVE and the envelope average amplitude. We use this space to analyze rat EEG recordings during sleep-wake cycles. Our results show that delta, theta and sigma bands approach Gaussianity at the lowest EEG amplitudes while exhibiting significant deviations at high EEG amplitudes. Deviations to low-CVE appeared prominently during REM sleep, associated with theta rhythm, a regime consistent with the dynamics shown by the synchronization of weakly coupled oscillators. On the other hand, deviations to high-CVE, appearing mostly during NREM sleep associated with EEG phasic activity and high-amplitude Gaussian waves, can be interpreted as the arrhythmic superposition of transient neural synchronization events. These two different manifestations of neural synchrony (low-CVE/high-CVE) explain the well-known spectral differences between REM and NREM sleep, while also illuminating the origin of the EEG 1 / f spectrum.

Diffusion-weighted (DW) functional magnetic resonance imaging (fMRI) signal changes have been noted as a promising marker of neural activity. Although there is no agreement on the signal origin, the blood oxygen level dependent (BOLD) effect has figured as one of the most likely sources. In order to investigate possible BOLD and non-BOLD contributions to the signal, DW fMRI was performed on normal volunteers using a sequence with two echo-planar acquisitions after pulsed-gradient spin-echo. Along with the changes to the signal amplitude (ΔS/S) measured at both echo-times, this sequence allowed changes to the transverse relaxation rate (ΔR 2 ) to be estimated for multiple b-values during hypercapnia (HC) and visual stimulation (VS). ΔS/S and ΔR 2 observed during HC were relatively insensitive to increasing b-value. On the other hand, ΔS/S demonstrated a clear dependence on b-value at both echo-times for VS. In addition, ΔR 2 during the latter half of VS was significantly more negative at b = 1400 s/mm 2 than for the time-courses at lower b-value, but ΔR 2 during the post-stimulus undershoot was independent of b-value. The results have been discussed in terms of two models: the standard intravascular-extravascular model for fMRI and a three-compartment model (one intra- and two extravascular compartments). Within these interpretations the results suggest that the majority of the response is linked to changes in transverse relaxation, but possible contributions from other sources may not be ruled out.

Behavioral inhibition often is studied by comparing the electroencephalographic responses to stop and to go signals. Most studies simply assess amplitude differences of the N200 and P300 event-related potentials, which seem to best correspond to increased activity in the theta and delta frequency bands, respectively. However, neither have reliable indicators for successful behavioral inhibition been identified nor have the causal dependencies of stop-related neurocognitive processes been addressed yet. By studying functional and effective connectivity underlying stopping behavior, this study opens new directions for the investigation of behavioral inhibition. Group independent component analysis was used to infer functionally coherent networks from electroencephalographic data, which were recorded from healthy human participants during processing of a stop signal task. Then, the temporal dynamics of causal dependencies between independent components were identified by means of Bayesian network estimations. The mean clustering coefficient and the characteristic path length measure indicated time windows between 130 and 180 ms and between 420 and 500 ms to express significantly different connectivity profiles between conditions. Three components showed significant correlations between 120 and 260 ms with stop signal reaction times and the number of failed stops. Two of these components acted as sources of causal flow, one capturing P300/delta characteristics while the other was characterized by alpha power depletion putatively representing the evaluation or processing of stimulus features. Although results suggest that the P300 and associated delta activity seem to be statistically dependent on earlier processes associated with behavioral inhibition, the time window critical for inhibition coincides with early changes in causal patterns and largely precedes peak amplitude differences between go and stop trials. Altogether, utilizing the analysis of stopping-related connectivity, previously undetected patterns emerged that warrant further investigation.

Spontaneous (or resting-state) brain activity has attracted a growing body of neuroimaging research over the last/ndecades.Whole-brain networkmodels have proved helpful to investigate the source of slow(b0.1 Hz) correlated/nhemodynamic fluctuations revealed in fMRI during rest. However, the mechanisms mediating resting-state/nlong-distance correlations and the relationship with the faster neural activity remain unclear. Novel insights/ncoming from MEG studies have shown that the amplitude envelopes of alpha- and beta-frequency oscillations/n(8–30 Hz) display similar correlation patterns as the fMRI signals./nIn thiswork, we combine experimental and theoreticalwork to investigate the mechanisms of spontaneousMEG/nfunctional connectivity. Using a simple model of coupled oscillators adapted to incorporate realisticwhole-brain/nconnectivity and conduction delays, we explore how slow and structured amplitude envelopes of band-pass/nfiltered signals – fairly reproducing MEG data collected from 10 healthy subjects at rest – are generated spontaneously/nin the space-time structure of the brain network./nOur simulation results show that the large-scale neuroanatomical connectivity provides an optimal network/nstructure to support a regimewith metastable synchronization. In this regime, different subsystems may temporarily/nsynchronize at reduced collective frequencies (falling in the 8–30 Hz range due to the delays) while the/nglobal system never fully synchronizes. This mechanism modulates the frequency of the oscillators on a slow/ntime-scale (b0.1 Hz) leading to structured amplitude fluctuations of band-pass filtered signals. Taken overall,/nour results reveal that the structured amplitude envelope fluctuations observed in resting-state MEG data may/noriginate from spontaneous synchronization mechanisms naturally occurring in the space-time structure of/nthe brain. The research reported herein was supported by the ERC Advanced/nGrant: DYSTRUCTURE (n. 295129), by the FET Flagship Human Brain/nProject, by the Spanish Research Project SAF2010-16085, by the/nCONSOLIDER-INGENIO 2010 CSD2007-00012, by the BrainNRG through/nthe James S. McDonnell Foundation, by the FP7-ICT BrainScales, by the/nRCUK Digital Economy – Centre for Doctoral Training in Healthcare/nInnovation, by theMINDLab Investment Capital for University Research/nFund and by the TrygFonden Charitable Foundation.

Investigating the neural generators of auditory steady-state responses (ASSRs), i.e., auditory evoked brain responses, with a wide range of screening and diagnostic applications, has been the focus of various studies for many years. Most of these studies employed a priori assumptions regarding the number and location of neural generators. The aim of this study is to reconstruct ASSR sources with minimal assumptions in order to gain in-depth insight into the number and location of brain regions that are activated in response to low- as well as high-frequency acoustically amplitude modulated signals. In order to reconstruct ASSR sources, we applied independent component analysis with subsequent equivalent dipole modeling to single-subject EEG data (young adults, 20–30 years of age). These data were based on white noise stimuli, amplitude modulated at 4, 20, 40, or 80 Hz. The independent components that exhibited a significant ASSR were clustered among all participants by means of a probabilistic clustering method based on a Gaussian mixture model. Results suggest that a widely distributed network of sources, located in cortical as well as subcortical regions, is active in response to 4, 20, 40, and 80 Hz amplitude modulated noises. Some of these sources are located beyond the central auditory pathway. Comparison of brain sources in response to different modulation frequencies suggested that the identified brain sources in the brainstem, the left and the right auditory cortex show a higher responsiveness to 40 Hz than to the other modulation frequencies.

Visually induced motion sickness (VIMS) can occur via prolonged exposure to visual stimulation that generates the illusion of self-motion (vection). Not everyone is susceptible to VIMS and the neural mechanism underlying susceptibility is unclear. This study explored the differences of electroencephalographic (EEG) signatures between VIMS-susceptible and VIMS-resistant groups. Thirty-two-channel EEG data were recorded from 12 VIMS-susceptible and 15 VIMS-resistant university students while they were watching two patterns of moving dots: (1) a coherent rotation pattern (vection-inducing and potentially VIMS-provoking pattern), and (2) a random movement pattern (non-VIMS-provoking control). The VIMS-susceptible group exhibited a significantly larger increase in the parietal N2 response when exposed to the coherent rotating pattern than when exposed to control patterns. In members of the VIMS-resistant group, before vection onset, global connectivity from all other EEG electrodes to the right-temporal-parietal and to the right-central areas increased, whereas after vection onset the global connectivity to the right-frontal area reduced. Such changes were not observed in the susceptible group. Further, the increases in N2 amplitude and the identified phase synchronization index were significantly correlated with individual motion sickness susceptibility. Results suggest that VIMS susceptibility is associated with systematic impairment of dynamic cortical coordination as captured by the phase synchronization of cortical activities. Analyses of dynamic EEG signatures could be a means to unlock the neural mechanism of VIMS.

Analyzing single trial brain activity remains a challenging problem in the neurosciences. We gain purchase on this problem by focusing on globally synchronous fields in within-trial evoked brain activity, rather than on localized peaks in the trial-averaged evoked response (ER). We analyzed data from three measurement modalities, each with different spatial resolutions: magnetoencephalogram (MEG), electroencephalogram (EEG) and electrocorticogram (ECoG). We first characterized the ER in terms of summation of phase and amplitude components over trials. Both contributed to the ER, as expected, but the ER topography was dominated by the phase component. This means the observed topography of cross-trial phase will not necessarily reflect the phase topography within trials. To assess the organization of within-trial phase, traveling wave (TW) components were quantified by computing the phase gradient. TWs were intermittent but ubiquitous in the within-trial evoked brain activity. At most task-relevant times and frequencies, the within-trial phase topography was described better by a TW than by the trial-average of phase. The trial-average of the TW components also reproduced the topography of the ER; we suggest that the ER topography arises, in large part, as an average over TW behaviors. These findings were consistent across the three measurement modalities. We conclude that, while phase is critical to understanding the topography of event-related activity, the preliminary step of collating cortical signals across trials can obscure the TW components in brain activity and lead to an underestimation of the coherent motion of cortical fields.

EEG recordings made during concurrent fMRI are confounded by the pulse artefact (PA), which although smaller than the gradient artefact is often more problematic because of its variability over multiple cardiac cycles. A better understanding of the PA is needed in order to generate improved methods for reducing its effect in EEG–fMRI experiments. Here we performed a study aimed at identifying the relative contributions of three putative sources of the PA (cardiac-pulse-driven head rotation, the Hall effect due to pulsatile blood flow and pulse-driven expansion of the scalp) to its amplitude and variability. EEG recordings were made from 6 subjects lying in a 3 T scanner. Accelerometers were fixed on the forehead and temple to monitor head motion. A bite-bar and vacuum cushion were used to restrain the head, thus greatly attenuating the contribution of cardiac-driven head rotation to the PA, while an insulating layer placed between the head and the EEG electrodes was used to eliminate the Hall voltage contribution. Using the root mean square (RMS) amplitude of the PA averaged over leads and time as a measure of the PA amplitude, we found that head restraint and insulating layer reduced the PA by 61% and 42%, respectively, when compared with the PA induced with the subject relaxed, indicating that cardiac-pulse-driven head rotation is the dominant source of the PA. With both the insulating layer and head restraint in place, the PA was reduced in RMS amplitude by 78% compared with the relaxed condition, the remaining PA contribution resulting from scalp expansion or residual head motion. The variance of the PA across cardiac cycles was more strongly reduced by the insulating layer than the head restraint, indicating that the flow-induced Hall voltage makes a larger contribution than pulse-driven head rotation to the variability of the PA., Highlights ► We experimentally isolate the putative sources of the pulse artefact (PA) in humans. ► Head rotation is the dominant source its removal produces a 61% reduction of the PA. ► Removing Hall effect produces a considerable, 42%, reduction in PA. ► Hall effect produces largest contribution to PA variance across cardiac cycles.

There is an increasing demand for source analysis of neonatal EEG, but currently there is inadequate knowledge about i) the spatial patterning of neonatal scalp EEG and hence ii) the number of electrodes needed to capture neonatal EEG in full spatial detail. This study addresses these issues by using a very high density (2.5 mm interelectrode spacing) linear electrode array to assess the spatial power spectrum, by using a high density (64 electrodes) EEG cap to assess the spatial extent of the common oscillatory bouts in the neonatal EEG and by using a neonatal size spherical head model to assess the effects of source depth and skull conductivities on the spatial frequency spectrum. The linear array recordings show that the spatial power spectrum decays rapidly until about 0.5–0.8 cycles per centimeter. The dense array EEG recordings show that the amplitude of oscillatory events decays within 4–6 cm to the level of global background activity, and that the higher frequencies (12–20 Hz) show the most rapid spatial decline in amplitude. Simulation with spherical head model showed that realistic variation in skull conductivity and source depths can both introduce orders of magnitude difference in the spatial frequency of the scalp EEG. Calculation of spatial Nyquist frequencies from the spatial power spectra suggests that an interelectrode distance of about 6–10 mm would suffice to capture the full spatial texture of the raw EEG signal at the neonatal scalp without spatial aliasing or under-sampling. The spatial decay of oscillatory events suggests that a full representation of their spatial characteristics requires an interelectrode distance of 10–20 mm. The findings show that the conventional way of recording neonatal EEG with about 10 electrodes ignores most spatial EEG content, that increasing the electrode density is necessary to improve neonatal EEG source localization and information extraction, and that prospective source models will need to carefully consider the neonatally relevant ranges of tissue conductivities and source depths when source localizing cortical activity in neonates.

Phase/amplitude coupling (PAC) is emerging as an important electrophysiological measure of local and long-distance neuronal communication. Current techniques for calculating PAC provide a numerical index that represents an average value across an arbitrarily long time period. This requires researchers to rely on block design experiments and temporal concatenation at the cost of the sub-second temporal resolution afforded by electrophysiological recordings. Here we present a method for calculating event-related phase/amplitude coupling (ERPAC) designed to capture the temporal evolution of task-related changes in PAC across events or between distant brain regions that is applicable to human or animal electromagnetic recording.

Corticomuscular interactions are studied mostly with EEG/EMG coherence, which, however, does not allow quantification of amplitude dynamics of sensorimotor oscillations. Here, we investigated the amplitude dynamics of sensorimotor EEG beta oscillations during an isometric task and their relation to corticomuscular coherence (CMC). We used amplitude envelopes of beta oscillations, derived from multichannel EEG and EMG recordings, as a measure of local cortical and spinal-cord synchronization. In general, we showed that the amplitude of cortical beta oscillations can influence CMC in two ways. First, we showed that the signal-to-noise ratio of pre-stimulus beta oscillations affects CMC. Second, we demonstrated that the attenuation of beta oscillations upon imperative stimulus correlated with the CMC strength. Attenuation of cortical beta oscillations was previously hypothesized to reflect increased motor cortex excitability. Consequently, this correlation might indicate that high cortical excitability, produced by imperative stimulus, facilitates the recruitment of neuronal networks responsible for establishing reliable corticospinal control manifested in larger CMC. Critically, we demonstrated that the amplitude envelopes of beta oscillations in EEG and EMG are positively correlated on time scales ranging from 50 to 1000 ms. Such correlations indicate that the amplitude of cortical beta oscillations might relate to the rhythmic spiking output of both corticospinal neurons and their spinal targets. Compared to CMC, however, amplitude-envelope correlations were detected in fewer cases, which might relate to a higher susceptibility of these correlations to signal-to-noise ratio. We conclude that EEG beta oscillations, originating from the sensorimotor cortex, can transmit not only their phase but also amplitude dynamics through the spinal motoneurons down to peripheral effectors.

With increasing age, cognitive control processes steadily decline. Prior research suggests that healthy older adults have a generally intact performance monitoring system, but show specific deficits in error awareness, i.e., the ability to detect committed errors. We examined the neural processing of errors across the adult lifespan (69 participants; age range 20–72 years) by analysing the error (-related) negativity (Ne/ERN) and the error positivity (Pe) using an adapted version of the Go/Nogo task. At a stable overall error rate, higher age was associated with a greater proportion of undetected errors. While the Ne/ERN was associated with the processing of errors in general, the Pe amplitude was modulated by detected errors only. Furthermore, the Pe amplitude for detected errors was significantly smaller in older adults, in contrast to the Ne/ERN amplitude which did not show age-related changes. Structural path models suggested that through those age-related changes in Pe amplitude, an indirect effect on the performance was observed. Our results confirm and extend previous extreme-group based findings about specific deficits in error detection associated with higher age using age as a continuous predictor. Age-related reductions in Pe amplitude, associated with more undetected errors, are independent of early error processing, as evidenced by the preserved Ne/ERN.

Encoding higher spatial resolution in simultaneous multi-slice (SMS) EPI is highly dependent on gradient performance, high density receiver coil arrays and pulse sequence optimization. We simulate gradient amplitude and slew rate determination of EPI imaging performance in terms of minimum TE, echo spacing (ES) and spatial resolution. We discuss the effects of image zooming in pulse sequences that have been used for sub-millimeter resolutions and the trade-offs in using partial Fourier and parallel imaging to reduce TE, PSF and ES. Using optimizations for SMS EPI pulse sequences with available gradient and receiver hardware, experimental results in ultra-high resolution (UHR) (0.45-0.5mm isotropic) SMS-EPI fMRI and mapping ocular dominance columns (ODC) in human brain at 0.5 mm isotropic resolution are demonstrated. We discuss promising future directions of UHR fMRI.

A fundamental feature of the temporal organization of neural activity is phase-amplitude coupling between brain rhythms at different frequencies, where the amplitude of a higher frequency varies according to the phase of a lower frequency. Here, we show that this rule extends to brain-organ interactions. We measured both the infra-slow (~0.05 Hz) rhythm intrinsically generated by the stomach – the gastric basal rhythm – using electrogastrography, and spontaneous brain dynamics with magnetoencephalography during resting-state with eyes open. We found significant phase-amplitude coupling between the infra-slow gastric phase and the amplitude of the cortical alpha rhythm (10–11 Hz), with gastric phase accounting for 8% of the variance of alpha rhythm amplitude fluctuations. Gastric-alpha coupling was localized to the right anterior insula, and bilaterally to occipito-parietal regions. Transfer entropy, a measure of directionality of information transfer, indicates that gastric-alpha coupling is due to an ascending influence from the stomach to both the right anterior insula and occipito-parietal regions. Our results show that phase-amplitude coupling so far only observed within the brain extends to brain-viscera interactions. They further reveal that the temporal structure of spontaneous brain activity depends not only on neuron and network properties endogenous to the brain, but also on the slow electrical rhythm generated by the stomach., Highlights • Spontaneous alpha rhythm fluctuations locked to gastric phase. • Gastric phase explains 8% of alpha variance. • Gastric-alpha coupling mostly to due to ascending influences from stomach to brain. • The stomach as an external oscillator constraining spontaneous brain activity.

Movements result from a complex interplay of multiple brain regions. These regions are assembled into distinct functional networks depending on the specific properties of the action. However, the nature and details of the dynamics of this complex assembly process are unknown. In this study, we sought to identify key markers of the neural processes underlying the preparation and execution of motor actions that always occur irrespective of differences in movement initiation, hence the specific neural processes and functional networks involved. To this end, EEG activity was continuously recorded from 18 right-handed healthy participants while they performed a simple motor task consisting of button presses with the left or right index finger. The movement was performed either in response to a visual cue or at a self-chosen, i.e., non-cued point in time. Despite these substantial differences in movement initiation, dynamic properties of the EEG signals common to both conditions could be identified using time–frequency and phase locking analysis of the EEG data. In both conditions, a significant phase locking effect was observed that started prior to the movement onset in the δ – θ frequency band (2–7 Hz), and that was strongest at the electrodes nearest to the contralateral motor region (M1). This phase locking effect did not have a counterpart in the corresponding power spectra (i.e., amplitudes), or in the event-related potentials. Our finding suggests that phase locking in the δ – θ frequency band is a ubiquitous movement-related signal independent of how the actual movement has been initiated. We therefore suggest that phase-locked neural oscillations in the motor cortex are a prerequisite for the preparation and execution of motor actions.