Auditory cortex modelled as a dynamical network of oscillators: understanding event-related fields and their adaptation.

Adaptation, the reduction of neuronal responses by repetitive stimulation, is a ubiquitous feature of auditory cortex (AC). It is not clear what causes adaptation, but short-term synaptic depression (STSD) is a potential candidate for the underlying mechanism. In such a case, adaptation can be directly linked with the way AC produces context-sensitive responses such as mismatch negativity and stimulus-specific adaptation observed on the single-unit level. We examined this hypothesis via a computational model based on AC anatomy, which includes serially connected core, belt, and parabelt areas. The model replicates the event-related field (ERF) of the magnetoencephalogram as well as ERF adaptation. The model dynamics are described by excitatory and inhibitory state variables of cell populations, with the excitatory connections modulated by STSD. We analysed the system dynamics by linearising the firing rates and solving the STSD equation using time-scale separation. This allows for characterisation of AC dynamics as a superposition of damped harmonic oscillators, so-called normal modes. We show that repetition suppression of the N1m is due to a mixture of causes, with stimulus repetition modifying both the amplitudes and the frequencies of the normal modes. In this view, adaptation results from a complete reorganisation of AC dynamics rather than a reduction of activity in discrete sources. Further, both the network structure and the balance between excitation and inhibition contribute significantly to the rate with which AC recovers from adaptation. This lifetime of adaptation is longer in the belt and parabelt than in the core area, despite the time constants of STSD being spatially homogeneous. Finally, we critically evaluate the use of a single exponential function to describe recovery from adaptation.

Explaining event-related fields by a mechanistic model encapsulating the anatomical structure of auditory cortex.

Event-related fields of the magnetoencephalogram are triggered by sensory stimuli and appear as a series of waves extending hundreds of milliseconds after stimulus onset. They reflect the processing of the stimulus in cortex and have a highly subject-specific morphology. However, we still have an incomplete picture of how event-related fields are generated, what the various waves signify, and why they are so subject-specific. Here, we focus on this problem through the lens of a computational model which describes auditory cortex in terms of interconnected cortical columns as part of hierarchically placed fields of the core, belt, and parabelt areas. We develop an analytical approach arriving at solutions to the system dynamics in terms of normal modes: damped harmonic oscillators emerging out of the coupled excitation and inhibition in the system. Each normal mode is a global feature which depends on the anatomical structure of the entire auditory cortex. Further, normal modes are fundamental dynamical building blocks, in that the activity of each cortical column represents a combination of all normal modes. This approach allows us to replicate a typical auditory event-related response as a weighted sum of the single-column activities. Our work offers an alternative to the view that the event-related field arises out of spatially discrete, local generators. Rather, there is only a single generator process distributed over the entire network of the auditory cortex. We present predictions for testing to what degree subject-specificity is due to cross-subject variations in dynamical parameters rather than in the cortical surface morphology.

Simultaneous spatio-temporal matching pursuit decomposition of evoked brain responses in MEG.

We present a novel approach to the spatio-temporal decomposition of evoked brain responses in magnetoencephalography (MEG) aiming at a sparse representation of the underlying brain activity in terms of spatio-temporal atoms. Our approach is characterized by three attributes which constitute significant improvements with respect to existing approaches: (1) the spatial and temporal decomposition is addressed simultaneously rather than sequentially, with the benefit that source loci and corresponding waveforms can be unequivocally allocated to each other, and, hence, allow a plausible physiological interpretation of the parametrized data; (2) it is free from severe a priori assumptions about the solution space; (3) it comprises an optimization technique for the use of very large spatial and temporal subdirectories to greatly reduce the otherwise enormous computational cost by making use of the Cauchy-Schwarz inequality. We demonstrate the efficiency of the approach with simulations and real MEG data obtained from a subject exposed to a simple auditory stimulus.

Averaging auditory evoked magnetoencephalographic and electroencephalographic responses: a critical discussion.

In the analysis of data from magnetoencephalography (MEG) and electroencephalography (EEG), it is common practice to arithmetically average event-related magnetic fields (ERFs) or event-related electric potentials (ERPs) across single trials and subsequently across subjects to obtain the so-called grand mean. Comparisons of grand means, e.g. between conditions, are then often performed by subtraction. These operations, and their statistical evaluation with parametric tests such as ANOVA, tacitly rely on the assumption that the data follow the additive model, have a normal distribution, and have a homogeneous variance. This may be true for single trials, but these conditions are rarely met when ERFs/ERPs are compared between subjects, meaning that the additive model is seldom the correct model for computing grand mean waveforms. Here, we summarize some of our recent work and present new evidence, from auditory-evoked MEG and EEG results, that the non-normal distributions and the heteroscedasticity observed instead result because ERFs/ERPs follow a mixed model with additive and multiplicative components. For peak amplitudes, such as the auditory M100 and N100, the multiplicative component dominates. These findings emphasize that the common practice of simply subtracting arithmetic means of auditory-evoked ERFs or ERPs is problematic without prior adequate transformation of the data. Application of the area sinus hyperbolicus (asinh) transform to data following the mixed model transforms them into the requested additive model with its normal distribution and homogeneous variance. We therefore advise checking the data for compliance with the additive model and using the asinh transform if required.

Optimising a computational model of human auditory cortex with an evolutionary algorithm.

We demonstrate how the structure of auditory cortex can be investigated by combining computational modelling with advanced optimisation methods. We optimise a well-established auditory cortex model by means of an evolutionary algorithm. The model describes auditory cortex in terms of multiple core, belt, and parabelt fields. The optimisation process finds the optimum connections between individual fields of auditory cortex so that the model is able to reproduce experimental magnetoencephalographic (MEG) data. In the current study, this data comprised the auditory event-related fields (ERFs) recorded from a human subject in an MEG experiment where the stimulus-onset interval between consecutive tones was varied. The quality of the match between synthesised and experimental waveforms was 98%. The results suggest that neural activity caused by feedback connections plays a particularly important role in shaping ERF morphology. Further, ERFs reflect activity of the entire auditory cortex, and response adaptation due to stimulus repetition emerges from a complete reorganisation of AC dynamics rather than a reduction of activity in discrete sources. Our findings constitute the first stage in establishing a new non-invasive method for uncovering the organisation of the human auditory cortex.

Application of modern tests for stationarity to single-trial MEG data: transferring powerful statistical tools from econometrics to neuroscience.

Stationarity is a crucial yet rarely questioned assumption in the analysis of time series of magneto- (MEG) or electroencephalography (EEG). One key drawback of the commonly used tests for stationarity of encephalographic time series is the fact that conclusions on stationarity are only indirectly inferred either from the Gaussianity (e.g. the Shapiro-Wilk test or Kolmogorov-Smirnov test) or the randomness of the time series and the absence of trend using very simple time-series models (e.g. the sign and trend tests by Bendat and Piersol). We present a novel approach to the analysis of the stationarity of MEG and EEG time series by applying modern statistical methods which were specifically developed in econometrics to verify the hypothesis that a time series is stationary. We report our findings of the application of three different tests of stationarity--the Kwiatkowski-Phillips-Schmidt-Schin (KPSS) test for trend or mean stationarity, the Phillips-Perron (PP) test for the presence of a unit root and the White test for homoscedasticity--on an illustrative set of MEG data. For five stimulation sessions, we found already for short epochs of duration of 250 and 500 ms that, although the majority of the studied epochs of single MEG trials were usually mean-stationary (KPSS test and PP test), they were classified as nonstationary due to their heteroscedasticity (White test). We also observed that the presence of external auditory stimulation did not significantly affect the findings regarding the stationarity of the data. We conclude that the combination of these tests allows a refined analysis of the stationarity of MEG and EEG time series.

Reducing the Number of MEG/EEG Trials Needed for the Estimation of Brain Evoked Responses: A Bootstrap Approach.

OBJECTIVE: A common problem in magnetoencephalographic (MEG) and electroencephalographic (EEG) experimental paradigms relying on the estimation of brain evoked responses is the lengthy time of the experiment, which stems from the need to acquire a large number of repeated recordings. Using a bootstrap approach, we aim at reliably reducing the number of these repeated trials. METHODS: To this end, we assessed five variants of non-parametric bootstrapping based on the classical signal-plus-noise model constituting the foundation of signal averaging in MEG/EEG. We explain which of these approaches should and which should not be used for the aforementioned purpose, and why. RESULTS: We present results for two advocated bootstrap variants applied to auditory MEG data. The ensuing trial-averaged magnetic fields served as input to the estimation of cortical source generators, with spatio-temporal matching pursuit as an example of an inverse solution technique. We propose, for a wide range of trial numbers, a general framework to evaluate the statistical properties of the parameter estimates for source locations and related time courses. CONCLUSION: The proposed bootstrap framework offers a systematic approach to reduce the number of trials required to estimate the evoked response. The general validity of our findings is neither bound to any particular type of MEG/EEG data nor to any specific source localization method. SIGNIFICANCE: Practical implications of this work relate to the optimization of acquisition time of MEG/EEG experiments, thus reducing stress for the subjects (especially for patients) and minimizing related artifacts.

Why do humans have unique auditory event-related fields? Evidence from computational modeling and MEG experiments.

Auditory event-related fields (ERFs) measured with magnetoencephalography (MEG) are useful for studying the neuronal underpinnings of auditory cognition in human cortex. They have a highly subject-specific morphology, albeit certain characteristic deflections (e.g., P1m, N1m, and P2m) can be identified in most subjects. Here, we explore the reason for this subject-specificity through a combination of MEG measurements and computational modeling of auditory cortex. We test whether ERF subject-specificity can predominantly be explained in terms of each subject having an individual cortical gross anatomy, which modulates the MEG signal, or whether individual cortical dynamics is also at play. To our knowledge, this is the first time that tools to address this question are being presented. The effects of anatomical and dynamical variation on the MEG signal is simulated in a model describing the core-belt-parabelt structure of the auditory cortex, and with the dynamics based on the leaky-integrator neuron model. The experimental and simulated ERFs are characterized in terms of the N1m amplitude, latency, and width. Also, we examine the waveform grand-averaged across subjects, and the standard deviation of this grand average. The results show that the intersubject variability of the ERF arises out of both the anatomy and the dynamics of auditory cortex being specific to each subject. Moreover, our results suggest that the latency variation of the N1m is largely related to subject-specific dynamics. The findings are discussed in terms of how learning, plasticity, and sound detection are reflected in the auditory ERFs. The notion of the grand-averaged ERF is critically evaluated.

Effects of the task of categorizing FM direction on auditory evoked magnetic fields in the human auditory cortex.

We examined effects of the task of categorizing linear frequency-modulated (FM) sweeps into rising and falling on auditory evoked magnetic fields (AEFs) from the human auditory cortex, recorded by means of whole-head magnetoencephalography. AEFs in this task condition were compared with those in a passive condition where subjects had been asked to just passively listen to the same stimulus material. We found that the M100-peak latency was significantly shorter for the task condition than for the passive condition in the left but not in the right hemisphere. Furthermore, the M100-peak latency was significantly shorter in the right than in the left hemisphere for the passive and the task conditions. In contrast, the M100-peak amplitude did not differ significantly between conditions, nor between hemispheres. We also analyzed the activation strength derived from the integral of the absolute magnetic field over constant time windows between stimulus onset and 260 ms. We isolated an early, narrow time range between about 60 ms and 80 ms that showed larger values in the task condition, most prominently in the right hemisphere. These results add to other imaging and lesion studies which suggest a specific role of the right auditory cortex in identifying FM sweep direction and thus in categorizing FM sweeps into rising and falling.

A multilevel and cross-modal approach towards neuronal mechanisms of auditory streaming.

We report first results of a multilevel, cross-modal study on the neuronal mechanisms underlying auditory sequential streaming, with the focus on the impact of visual sequences on perceptually ambiguous tone sequences which can either be perceived as two separate streams or one alternating stream. We combined two psychophysical experiments performed on humans and monkeys with two human brain imaging experiments which allow to obtain complementary information on brain activation with high spatial (fMRI) and high temporal (MEG) resolution. The same acoustic paradigm based on the pairing of tone sequences with visual stimuli was used in all human studies and, in an adapted version, in the psychophysical study on monkeys. Our multilevel approach provides experimental evidence that the pairing of auditory and visual stimuli can reliably introduce a bias towards either an integrated or a segregated perception of ambiguous sequences. Thus, comparable to an explicit instruction, this approach can be used to control the subject's perceptual organization of an ambiguous sound sequence without the need for the subject to directly report it. This finding is of particular importance for animal studies because it allows to compare electrophysiological responses of auditory cortex neurons to the same acoustic stimulus sequence eliciting either a segregated or integrated percept.

Improving Accuracy and Temporal Resolution of Learning Curve Estimation for within- and across-Session Analysis.

Estimation of learning curves is ubiquitously based on proportions of correct responses within moving trial windows. Thereby, it is tacitly assumed that learning performance is constant within the moving windows, which, however, is often not the case. In the present study we demonstrate that violations of this assumption lead to systematic errors in the analysis of learning curves, and we explored the dependency of these errors on window size, different statistical models, and learning phase. To reduce these errors in the analysis of single-subject data as well as on the population level, we propose adequate statistical methods for the estimation of learning curves and the construction of confidence intervals, trial by trial. Applied to data from an avoidance learning experiment with rodents, these methods revealed performance changes occurring at multiple time scales within and across training sessions which were otherwise obscured in the conventional analysis. Our work shows that the proper assessment of the behavioral dynamics of learning at high temporal resolution can shed new light on specific learning processes, and, thus, allows to refine existing learning concepts. It further disambiguates the interpretation of neurophysiological signal changes recorded during training in relation to learning.

Persistent neural activity in auditory cortex is related to auditory working memory in humans and nonhuman primates.

Working memory is the cognitive capacity of short-term storage of information for goal-directed behaviors. Where and how this capacity is implemented in the brain are unresolved questions. We show that auditory cortex stores information by persistent changes of neural activity. We separated activity related to working memory from activity related to other mental processes by having humans and monkeys perform different tasks with varying working memory demands on the same sound sequences. Working memory was reflected in the spiking activity of individual neurons in auditory cortex and in the activity of neuronal populations, that is, in local field potentials and magnetic fields. Our results provide direct support for the idea that temporary storage of information recruits the same brain areas that also process the information. Because similar activity was observed in the two species, the cellular bases of some auditory working memory processes in humans can be studied in monkeys.

Variance stabilization for computing and comparing grand mean waveforms in MEG and EEG.

Grand means of time-varying signals (waveforms) across subjects in magnetoencephalography (MEG) and electroencephalography (EEG) are commonly computed as arithmetic averages and compared between conditions, for example, by subtraction. However, the prerequisite for these operations, homogeneity of the variance of the waveforms in time, and for most common parametric statistical tests also between conditions, is rarely met. We suggest that the heteroscedasticity observed instead results because waveforms may differ by factors and additive terms and follow a mixed model. We propose to apply the asinh-transformation to stabilize the variance in such cases. We demonstrate the homogeneous variance and the normal distributions of data achieved by this transformation using simulated waveforms, and we apply it to real MEG data and show its benefits. The asinh-transformation is thus an essential and useful processing step prior to computing and comparing grand mean waveforms in MEG and EEG.

Stimulation-history effects on the M100 revealed by its differential dependence on the stimulus onset interval.

The amplitudes of the most prominent component of auditory evoked magnetic fields and electrical potentials, the M100 and N100, recorded from the human scalp depend on the duration of the stimulus onset interval (SOI). Here, we show, using magnetoencephalography, that the SOI dependence of the M100 amplitude strongly depends upon whether stimuli with different SOIs are presented in a conventional block design or in a random manner. This differential dependence reveals that the M100 is affected not only by the stimulus evoking it and by its predecessor, but by a longer-term history of stimulation. We provide a parsimonious model that accounts for our findings with both designs in a quantitative manner. It assumes a transient, temporally asymmetric reduction in the excitability of a fraction of potentially excitable neurons. A rather stereotyped response function may therefore underlie the stimulation-history effects in the human auditory cortex.

The M100 component of evoked magnetic fields differs by scaling factors: implications for signal averaging.

MEG and EEG studies of event-related responses often involve comparisons of grand averages, requiring homogeneity of the variances. Here, we examine the possibility, implied by the nature of neural sources and the measuring principles involved, that the M100 component of auditory-evoked magnetic fields of different subjects, hemispheres, to different stimuli, and at different sensors differs by scaling factors. Such a multiplicative model predicts a linear increase in the standard deviation with the mean, and thus would have important implications for averaging and comparing such data. Our analyses, at the sensor and the source level, clearly show that the multiplicative model applies. We therefore propose geometric, rather than arithmetic, averaging of the M100 component across subjects and suggest a novel and superior normalization procedure. Our results question the justification of the common practice of subtracting arithmetic grand averages.

Single-trial evoked brain responses modeled by multivariate matching pursuit.

We present a new approach to the analysis of brain evoked electromagnetic potentials and fields. Multivariate version of the matching pursuit algorithm (MMP) performs an iterative, exhaustive search for waveforms, which optimally fit to signal structures, persistent in all the responses (trials) with the same time of occurrence, frequency, phase, and time width, but varying amplitude. The search is performed in a highly redundant time--frequency dictionary of Gabor functions, i.e., sines modulated by Gaussians. We present the feasibility of such a single-trial MMP analysis of the auditory M100 response, using an illustrative dataset acquired in a magnetoencephalographic (MEG) measurement with auditory stimulation with sinusoidal 1-kHz tones. We find that the morphology of the M100 estimate obtained from simple averaging of single trials can be very well explained by the average reconstruction with a few Gabor functions that parametrize those single trials. The M100 peak amplitude of single-trial reconstructions is observed to decrease with repetitions, which indicates habituation to the stimulus. This finding suggests that certain waveforms fitted by MMP could possibly be related to physiologically distinct components of evoked magnetic fields, which would allow tracing their dynamics on a single-trial level.