Signal to noise constraints on maximum length sequence auditory brain stem responses

Unfold: an integrated toolbox for overlap correction, non-linear modeling, and regression-based EEG analysis

Electrophysiological research with event-related brain potentials (ERPs) is increasingly moving from simple, strictly orthogonal stimulation paradigms towards more complex, quasi-experimental designs and naturalistic situations that involve fast, multisensory stimulation and complex motor behavior. As a result, electrophysiological responses from subsequent events often overlap with each other. In addition, the recorded neural activity is typically modulated by numerous covariates, which influence the measured responses in a linear or non-linear fashion. Examples of paradigms where systematic temporal overlap variations and low-level confounds between conditions cannot be avoided include combined electroencephalogram (EEG)/eye-tracking experiments during natural vision, fast multisensory stimulation experiments, and mobile brain/body imaging studies. However, even "traditional," highly controlled ERP datasets often contain a hidden mix of overlapping activity (e.g., from stimulus onsets, involuntary microsaccades, or button presses) and it is helpful or even necessary to disentangle these components for a correct interpretation of the results. In this paper, we introduce unfold, a powerful, yet easy-to-use MATLAB toolbox for regression-based EEG analyses that combines existing concepts of massive univariate modeling ("regression-ERPs"), linear deconvolution modeling, and non-linear modeling with the generalized additive model into one coherent and flexible analysis framework. The toolbox is modular, compatible with EEGLAB and can handle even large datasets efficiently. It also includes advanced options for regularization and the use of temporal basis functions (e.g., Fourier sets). We illustrate the advantages of this approach for simulated data as well as data from a standard face recognition experiment. In addition to traditional and non-conventional EEG/ERP designs, unfold can also be applied to other overlapping physiological signals, such as pupillary or electrodermal responses. It is available as open-source software at http://www.unfoldtoolbox.org.

Noise Attenuation Estimation for Maximum Length Sequences in Deconvolution Process of Auditory Evoked Potentials

The use of maximum length sequence (m-sequence) has been found beneficial for recovering both linear and nonlinear components at rapid stimulation. Since m-sequence is fully characterized by a primitive polynomial of different orders, the selection of polynomial order can be problematic in practice. Usually, the m-sequence is repetitively delivered in a looped fashion. Ensemble averaging is carried out as the first step and followed by the cross-correlation analysis to deconvolve linear/nonlinear responses. According to the classical noise reduction property based on additive noise model, theoretical equations have been derived in measuring noise attenuation ratios (NARs) after the averaging and correlation processes in the present study. A computer simulation experiment was conducted to test the derived equations, and a nonlinear deconvolution experiment was also conducted using order 7 and 9 m-sequences to address this issue with real data. Both theoretical and experimental results show that the NAR is essentially independent of the m-sequence order and is decided by the total length of valid data, as well as stimulation rate. The present study offers a guideline for m-sequence selections, which can be used to estimate required recording time and signal-to-noise ratio in designing m-sequence experiments.

Baseline correction of overlapping event-related responses using a linear deconvolution technique

The time it takes an event-related response (ERR) to subside is often longer than the interval between successive events so that the response to a new event interferes with a baseline formed by responses to preceding events. Thus, without proper baseline correction, the interpretation of an event-triggered average (ETA) of recorded data can be problematic. As the spectral compositions of ERR and baseline typically overlap, filtering the ETA is not always an adequate solution. The approach introduced here exploits that ETA and ERR are linearly related. Unless the series of events is exactly periodic, the ERR can be derived from the ETA by linear deconvolution. The performance of the method is illustrated with simulated examples as well as data from an auditory evoked field (AEF) study. It is also outlined how to handle experiments with two or more different events. Intriguing applications beyond the scope of baseline correction arise from the fact that the ERR is invariably estimated for a time window longer than the mean interval between successive events. The method may help, for example, to better understand the relationship between transient and steady-state responses or to delineate the component structure of a specific ERR.

A simulation study assessing the efficiency of deriving evoked responses using high stimulus rate

This study investigates the effects of stimulus rate on the recording efficiency of evoked-response scenario using a computer generated data simulating auditory evoked potentials (AEPs). We examine AEPs derived from three paradigms: 1) Conventional stimuli for ensemble averaging; 2) continuous loop averaging deconvolution (CLAD); and 3) maximum length sequence (MLS) stimuli, in which deconvolution techniques are incorporated to retrieve the overlapped responses. The performance is evaluated by correlation coefficients between the ideal and the derived responses. The results show that there is no distinct improvement of SNR for high rate stimulation paradigms, suggesting that the techniques of deconvolution in high rate will not provide more efficient recording approach. Further, responses from the CLAD paradigm tend to deteriorate compared to the other two methods due to the lower jittering arrangement of the sequence.

The effects of electrode montage on the amplitude of wave V in the auditory brainstem response to maximum length sequence stimuli

The use of maximum length sequence (MLS) stimuli to elicit an auditory brainstem response (ABR) has been limited, in part, by the observation that these stimuli reduce ABR wave amplitudes. This study recorded ABR waveforms from 14 normally hearing adults using MLS click stimuli (maximum stimulus rate = 250 clicks per second) at stimulus levels of 70, 60, 50, 40, 30 and 20 dB nHL, with a vertical and then an ipsilateral electrode montage. The vertical electrode montage produced significantly larger (p < 0.05) wave V amplitudes, with no change in wave V latencies (p > 0.05), at all stimulus levels. This result suggests a vertical electrode montage could be used to counter some of the loss in wave V amplitude observed when using MLS stimuli.

Signal to noise ratio analysis of maximum length sequence deconvolution of overlapping evoked potentials

In this study a general formula for the signal to noise ratio (SNR) of the maximum length sequence (MLS) deconvolution averaging is developed using the frequency domain framework of the generalized continuous loop averaging deconvolution procedure [Ozdamar and Bohórquez, J. Acoust. Soc. Am. 119, 429-438 (2006)]. This formulation takes advantage of the well known equivalency of energies in the time and frequency domains (Parseval's theorem) to show that in MLS deconvolution, SNR increases with the square root of half of the number of stimuli in the sweep. This increase is less than that of conventional averaging which is the square root of the number of sweeps averaged. Unlike arbitrary stimulus sequences that can attenuate or amplify phase unlocked noise depending on the frequency characteristics, the MLS deconvolution attenuates noise in all frequencies consistently. Furthermore, MLS and its zero-padded variations present optimal attenuation of noise at all frequencies yet they present a highly jittered stimulus sequence. In real recordings of evoked potentials, the time advantage gained by noise attenuation could be lost by the signal amplitude attenuation due to neural adaptation at high stimulus rates.

Signal-to-noise ratio and frequency analysis of continuous loop averaging deconvolution (CLAD) of overlapping evoked potentials

In this study, a frequency domain formulation of continuous loop averaging deconvolution (CLAD) of overlapping evoked potentials is developed and applied for the extraction of transient responses from recordings obtained at high stimulation rates. This formulation allows for a faster execution of CLAD by using fast Fourier transform algorithms. The frequency characteristics of the deconvolution filter depends exclusively on the stimulus sequence and determines whether the noncoherent noise is amplified or attenuated in different frequencies. A formula for calculating the signal-to-noise ratio (SNR) achieved by the deconvolution process is developed. The newly developed theory and the methodology is applied to the extraction of the auditory brainstem and middle latency responses using various sequences. The effects of the sequence used and the number of sweeps averaged in ongoing acquisition on SNR are examined by using single sweep recordings. The results verify the deconvolution theory and the methodology and show its limitations. Depending on the frequency characteristics of the sequence, the deconvolution process can amplify or attenuate the EEG noise. Proper selection of the stimulus sequence can increase the SNR enhancement obtained with conventional averaging.

The combined effects of forward masking by noise and high click rate on monaural and binaural human auditory nerve and brainstem potentials

OBJECTIVE

To study effects of forward masking and rapid stimulation on human monaurally- and binaurally-evoked brainstem potentials and suggest their relation to synaptic fatigue and recovery and to neuronal action potential refractoriness.

METHODS

Auditory brainstem evoked potentials (ABEPs) were recorded from 12 normally- and symmetrically hearing adults, in response to each click (50 dB nHL, condensation and rarefaction) in a train of nine, with an inter-click interval of 11 ms, that followed a white noise burst of 100 ms duration (50 dB nHL). Sequences of white noise and click train were repeated at a rate of 2.89 s(-1). The interval between noise and first click in the train was 2, 11, 22, 44, 66 or 88 ms in different runs. ABEPs were averaged (8000 repetitions) using a dwell time of 25 micros/address/channel. The binaural interaction components (BICs) of ABEPs were derived and the single, centrally located equivalent dipoles of ABEP waves I and V and of the BIC major wave were estimated.

RESULTS

The latencies of dipoles I and V of ABEP, their inter-dipole interval and the dipole magnitude of component V were significantly affected by the interval between noise and clicks and by the serial position of the click in the train. The latency and dipole magnitude of the major BIC component were significantly affected by the interval between noise and clicks. Interval from noise and the click's serial position in the train interacted to affect dipole V latency, dipole V magnitude, BIC latencies and the V-I inter-dipole latency difference. Most of the effects were fully apparent by the first few clicks in the train, and the trend (increase or decrease) was affected by the interval between noise and clicks.

CONCLUSIONS

The changes in latency and magnitude of ABEP and BIC components with advancing position in the click train and the interactions of click position in the train with the intervals from noise indicate an interaction of fatigue and recovery, compatible with synaptic depletion and replenishing, respectively. With the 2 ms interval between noise and the first click in the train, neuronal action potential refractoriness may also be involved.

Deconvolution of evoked responses obtained at high stimulus rates

Continuous loop averaging deconvolution (CLAD) is a new general mathematical theory and method developed to deconvolve overlapping auditory evoked responses obtained at high stimulation rates. Using CLAD, arbitrary stimulus sequences are generated and averaged responses deconvolved. Until now, only a few special stimulus series such as maximum length sequences (MLS) and Legendre sequences (LGS) were capable of performing this task. A CLAD computer algorithm is developed and implemented in an evoked potential averaging system. Computer simulations are used to verify the theory and methodology. Auditory brainstem responses (ABR) and middle latency responses (MLR) are acquired from subjects with normal hearing at high stimulation rates to validate and show the feasibility of the CLAD technique.

The cumulative effect of high click rate on monaural and binaural processing in the human auditory brainstem

OBJECTIVE

The objective of the present study was to compare the effects of high stimulus rate and click position in the train on monaurally and binaurally evoked activities in the human auditory brainstem and suggest their possible physiological mechanism.

METHODS

Auditory brainstem evoked potentials (ABEPs) were recorded from 15 normally and symmetrically hearing adults from 3 channels, in response to 50dB nHL, alternating polarity clicks, presented at a rate of 21/s as well as separately to each click in a train of 10 with an interstimulus interval of 11ms. Click trains were presented at a rate of 5.13/s. The binaural interaction components (BICs) of ABEPs were derived by subtracting the response to binaural clicks from the algebraic sum of monaural responses. Single, centrally located equivalent dipoles were estimated as concise measures of the surface-summated activity of ABEPs and BICs generators.

RESULTS

A significant effect of click position in the train on equivalent dipole latency of ABEP component V and on equivalent dipole magnitude of III were found. Latency was prolonged and amplitude was increased the later the click's position in the train. A significant effect of click position in the train on equivalent dipole latencies of all components of BICs was found. Latencies were prolonged if the click's position occurred later in the train, with most of the latency shift achieved by the third click in the train for the first major BIC and by the seventh click for other BIC components. No significant effects on equivalent dipole magnitudes of BICs were found. No significant effect of click position in the train on orientation of any of the equivalent dipoles of ABEP or BIC was found.

CONCLUSIONS

The progressive prolongation of latency of ABEP and BIC components with advancing position in the train may be attributed to cumulatively decreased synaptic efficacy at high stimulus rates, resulting in prolonged synaptic delays along the auditory pathway. The paradoxic enhancement of ABEP dipole III magnitude with advancing click position in the train may reflect higher sensitivity of inhibitory brainstem neurons to increased stimulus rate, resulting in disinhibition. The absence of significant effects on BIC dipole magnitudes may reflect the amplifying effect of divergence in the ascending auditory pathway, as has been observed for the monaurally evoked ABEP components from the upper pons.

The use of high stimulus rate auditory brainstem responses in the estimation of hearing threshold

This normative study investigates the efficiency of using the maximum length sequence (MLS) technique applied to auditory brainstem evoked response (ABR) testing to estimate hearing thresholds. Using a commercially available system, ABRs were recorded in sixteen subjects at two conventional rates--9.1 and 33.3 clicks/s--and six MLS rates between 88.8 and 1000 clicks/s. Each subject was tested at five stimulus levels from 60 down to 10 dBnHL. The wave JV amplitude input-output (I/O) functions, relative signal to noise ratio (SNR) and speed of test were calculated for all conditions. The JV amplitude and detectability decrease as the stimulus rate increases and level decreases. The latency of JV increases as the stimulus rate increases and the intensity decreases. While the slope of the amplitude I/O function was maximal at 200 clicks/s, at 300 clicks/s it was comparable with that obtained at conventional rates. At higher rates, the slope of the I/O function decreases. When compared with the conventional recording rate of 33.3 clicks/s there is a small improvement in SNR for MLS rates between 200 and 600 clicks/s at levels above 30 dBnHL. The calculated speed improvement at 300 clicks/s is a factor between 1.4 to 1.6 at a screening level of 30-40 dBnHL. It is felt therefore that there may be a small advantage to using MLS in screening and that the optimal rate for this lies at around 200 to 300 clicks/s. However even at these rates, as a consequence of the adaptation of the response with both rate and level, the improvement in SNR or speed of test would be modest when estimating threshold.

Nonlinear functional modeling of scalp recorded auditory evoked responses to maximum length sequences

The purpose of this study was to model the adult human's scalp recorded evoked response to auditory pulses separated by varying inter pulse intervals (IPIs). The responses modeled probably reflect auditory nerve and brainstem generators. The subjects were 10 young adult humans with normal hearing. They were presented pseudo random sequences of pulses (maximum length sequences, MLSs) in order to characterize their system response. For the stimuli and the responses modeled accounting for temporal nonlinearities (interactions among the pulses) improved model performance only marginally. Nonlinear contributions to the models decreased with increasing interval between the input pulses. Increasing the memory of the model beyond 20 ms did not increase modeled performance dramatically. Model performance varied as a function of minimum IPI (MIPI) of the MLSs. At the shortest MIPI overall model performance deteriorated (due, in part, to a decrease in SNR), but nonlinear effects became relatively more important. At the longest MIPI performance also deteriorated, possibly due to the increasing influence of longer latency, more variable evoked potential components. Modeled performance generalized to responses recorded in the same recording session to the same and different MLSs. This study confirms the similarity between MLS linear kernels and conventionally averaged evoked responses--both are adapted responses reflecting the IPIs of the evoking stimuli.

Responses of chopper units in the ventral cochlear nucleus of the anaesthetised guinea pig to clicks-in-noise and click trains

Auditory brainstem responses (ABRs) have been measured with clicks, clicks masked by noise, click trains and pseudorandom maximum length sequences (MLS) of clicks. To investigate the neuronal populations contributing to the ABR under these stimulation conditions, we measured the extracellular responses of ventral cochlear nucleus (VCN) units in the urethane-anaesthetised guinea pig. We studied 23 chopper, 7 primary-like and 7 onset units. This report focuses on the responses from chopper units. The probability of discharge for chopper units increased with increasing click level reaching nearly 100% in many units, over a range of about 20-30 dB. Following each response to a click there was a 5-10 ms suppression of the spontaneous or noise evoked activity. As the level of the noise was increased over a range of 20-30 dB, the response to the clicks gradually decreased leading to a complete abolition of the click response at high noise levels. In a few units, low level noise produced a facilitation of the response to single clicks. In response to constant level equally spaced click trains, discharge probability increased with increasing minimum pulse interval (MPI), approaching 100% for MPIs of 4-8 ms in some units. The recovery afforded by the gaps in the MLS train often resulted in higher discharge probability for MLS than click trains with the same MPI, while response probabilities for MLS and click trains were similar when compared at equivalent average click rates. At short MPIs (0.5 and 1.0 ms), peri stimulus time histograms in response to click trains resembled those to best frequency (BF) tones and noisebursts, with chopping peaks unrelated to unit BF. VCN units show highly synchronised and reliable responses to click trains, MLS trains and clicks masked by noise. The decrease in discharge rate and increase in latency of chopper units with decreasing click level, increasing click rate and increasing masker level parallel the peak amplitude and latency changes observed in the auditory brainstem response.

Auditory evoked brainstem and middle latency responses in Macaca mulatta and humans

Early (ABRs) and middle (MLRs) surface-recorded auditory evoked potentials were compared in eight adult monkeys (Macaca mulatta) and eight adult humans. Responses whose probable generators were the cochlear nucleus and lateral lemniscus were of shorter latency and larger amplitude in monkeys. Relative to humans, ABR response latencies in monkeys were less affected by stimulus intensity, stimulus rate, and masker level. In contrast, monkey amplitudes were relatively more affected by those same stimulus parameters. The most prominent MLR wave was longer in latency and greater in amplitude in humans than the homologous wave in monkeys. The reduction in amplitude of that wave with increasing rate was greater for humans than monkeys. Temporal interactions (the effect of prior stimuli on the response to current stimulation) were investigated from a non-linear systems identification framework using maximum length sequences (MLSs). Both monkey and human auditory systems were second and probably third-order systems at the levels assessed. As the separations between the stimulus pulses decreased, evidence for temporal interactions became more prominent, reached a maximum, and then decreased with further decreases in stimulus pulse separation. At the highest stimulus rates presented, variations in temporal spacing among stimuli had less of an effect on monkey than human evoked responses.