Selective processing of auditory evoked responses with iterative-randomized stimulation and averaging: A strategy for evaluating the time-invariant assumption

Multi-response deconvolution of auditory evoked potentials in a reduced representation space

The estimation of auditory evoked potentials requires deconvolution when the duration of the responses to be recovered exceeds the inter-stimulus interval. Based on least squares deconvolution, in this article we extend the procedure to the case of a multi-response convolutional model, that is, a model in which different categories of stimulus are expected to evoke different responses. The computational cost of the multi-response deconvolution significantly increases with the number of responses to be deconvolved, which restricts its applicability in practical situations. In order to alleviate this restriction, we propose to perform the multi-response deconvolution in a reduced representation space associated with a latency-dependent filtering of auditory responses, which provides a significant dimensionality reduction. We demonstrate the practical viability of the multi-response deconvolution with auditory responses evoked by clicks presented at different levels and categorized according to their stimulation level. The multi-response deconvolution applied in a reduced representation space provides the least squares estimation of the responses with a reasonable computational load. matlab/Octave code implementing the proposed procedure is included as supplementary material.

The hunt for hidden hearing loss in humans: From preclinical studies to effective interventions

Many individuals experience hearing problems that are hidden under a normal audiogram. This not only impacts on individual sufferers, but also on clinicians who can offer little in the way of support. Animal studies using invasive methodologies have developed solid evidence for a range of pathologies underlying this hidden hearing loss (HHL), including cochlear synaptopathy, auditory nerve demyelination, elevated central gain, and neural mal-adaptation. Despite progress in pre-clinical models, evidence supporting the existence of HHL in humans remains inconclusive, and clinicians lack any non-invasive biomarkers sensitive to HHL, as well as a standardized protocol to manage hearing problems in the absence of elevated hearing thresholds. Here, we review animal models of HHL as well as the ongoing research for tools with which to diagnose and manage hearing difficulties associated with HHL. We also discuss new research opportunities facilitated by recent methodological tools that may overcome a series of barriers that have hampered meaningful progress in diagnosing and treating of HHL.

Subspace-constrained deconvolution of auditory evoked potentials

Auditory evoked potentials can be estimated by synchronous averaging when the responses to the individual stimuli are not overlapped. However, when the response duration exceeds the inter-stimulus interval, a deconvolution procedure is necessary to obtain the transient response. The iterative randomized stimulation and averaging and the equivalent randomized stimulation with least squares deconvolution have been proven to be flexible and efficient methods for deconvolving the evoked potentials, with minimum restrictions in the design of stimulation sequences. Recently, a latency-dependent filtering and down-sampling (LDFDS) methodology was proposed for optimal filtering and dimensionality reduction, which is particularly useful when the evoked potentials involve the complete auditory pathway response (i.e., from the cochlea to the auditory cortex). In this case, the number of samples required to accurately represent the evoked potentials can be reduced from several thousand (with conventional sampling) to around 120. In this article, we propose to perform the deconvolution in the reduced representation space defined by LDFDS and present the mathematical foundation of the subspace-constrained deconvolution. Under the assumption that the evoked response is appropriately represented in the reduced representation space, the proposed deconvolution provides an optimal least squares estimation of the evoked response. Additionally, the dimensionality reduction provides a substantial reduction of the computational cost associated with the deconvolution. matlab/Octave code implementing the proposed procedures is included as supplementary material.

Optimizing Parameters for Using the Parallel Auditory Brainstem Response to Quickly Estimate Hearing Thresholds

OBJECTIVES

Timely assessments are critical to providing early intervention and better hearing and spoken language outcomes for children with hearing loss. To facilitate faster diagnostic hearing assessments in infants, the authors developed the parallel auditory brainstem response (pABR), which presents randomly timed trains of tone pips at five frequencies to each ear simultaneously. The pABR yields high-quality waveforms that are similar to the standard, single-frequency serial ABR but in a fraction of the recording time. While well-documented for standard ABRs, it is yet unknown how presentation rate and level interact to affect responses collected in parallel. Furthermore, the stimuli are yet to be calibrated to perceptual thresholds. Therefore, this study aimed to determine the optimal range of parameters for the pABR and to establish the normative stimulus level correction values for the ABR stimuli.

DESIGN

Two experiments were completed, each with a group of 20 adults (18-35 years old) with normal-hearing thresholds (≤20 dB HL) from 250 to 8000 Hz. First, pABR electroencephalographic (EEG) responses were recorded for six stimulation rates and two intensities. The changes in component wave V amplitude and latency were analyzed, as well as the time required for all responses to reach a criterion signal-to-noise ratio of 0 dB. Second, behavioral thresholds were measured for pure tones and for the pABR stimuli at each rate to determine the correction factors that relate the stimulus level in dB peSPL to perceptual thresholds in dB nHL.

RESULTS

The pABR showed some adaptation with increased stimulation rate. A wide range of rates yielded robust responses in under 15 minutes, but 40 Hz was the optimal singular presentation rate. Extending the analysis window to include later components of the response offered further time-saving advantages for the temporally broader responses to low-frequency tone pips. The perceptual thresholds to pABR stimuli changed subtly with rate, giving a relatively similar set of correction factors to convert the level of the pABR stimuli from dB peSPL to dB nHL.

CONCLUSIONS

The optimal stimulation rate for the pABR is 40 Hz but using multiple rates may prove useful. Perceptual thresholds that subtly change across rate allow for a testing paradigm that easily transitions between rates, which may be useful for quickly estimating thresholds for different configurations of hearing loss. These optimized parameters facilitate expediency and effectiveness of the pABR to estimate hearing thresholds in a clinical setting.

Latency-dependent filtering and compact representation of the complete auditory pathway response

Auditory evoked potentials (AEPs) include the auditory brainstem response (ABR), middle latency response (MLR), and cortical auditory evoked potentials (CAEPs), each one covering a specific latency range and frequency band. For this reason, ABR, MLR, and CAEP are usually recorded separately using different protocols. This article proposes a procedure providing a latency-dependent filtering and down-sampling of the AEP responses. This way, each AEP component is appropriately filtered, according to its latency, and the complete auditory pathway response is conveniently represented (with the minimum number of samples, i.e., without unnecessary redundancies). The compact representation of the complete response facilitates a comprehensive analysis of the evoked potentials (keeping the natural continuity related to the neural activity transmission along the auditory pathway), which provides a new perspective in the design and analysis of AEP experiments. Additionally, the proposed compact representation reduces the storage or transmission requirements when large databases are manipulated for clinical or research purposes. The analysis of the AEP responses shows that a compact representation with 40 samples/decade (around 120 samples) is enough for accurately representing the response of the complete auditory pathway and provides appropriate latency-dependent filtering. MatLab/Octave code implementing the proposed procedure is included in the supplementary materials.

Matrix-based formulation of the iterative randomized stimulation and averaging method for recording evoked potentials

The iterative randomized stimulation and averaging (IRSA) method was proposed for recording evoked potentials when the individual responses are overlapped. The main inconvenience of IRSA is its computational cost, associated with a large number of iterations required for recovering the evoked potentials and the computation required for each iteration [involving the whole electroencephalogram (EEG)]. This article proposes a matrix-based formulation of IRSA, which is mathematically equivalent and saves computational load (because each iteration involves just a segment with the length of the response, instead of the whole EEG). Additionally, it presents an analysis of convergence that demonstrates that IRSA converges to the least-squares (LS) deconvolution. Based on the convergence analysis, some optimizations for the IRSA algorithm are proposed. Experimental results (configured for obtaining the full-range auditory evoked potentials) show the mathematical equivalence of the different IRSA implementations and the LS-deconvolution and compare the respective computational costs of these implementations under different conditions. The proposed optimizations allow the practical use of IRSA for many clinical and research applications and provide a reduction of the computational cost, very important with respect to the conventional IRSA, and moderate with respect to the LS-deconvolution. matlab/Octave implementations of the different methods are provided as supplementary material.

Comparison of maximum length sequence and randomized stimulation and averaging methods on the bottlenose dolphin auditory brainstem response

The purpose of the present study was to compare auditory brainstem responses (ABRs) using two approaches that allow the use of high stimulation rates, but with different temporal variability in the interstimulus interval: maximum length sequences (MLS) and iterative randomized stimulation and averaging (I-RSA). ABRs were obtained to click stimuli in six bottlenose dolphins (Tursiops truncatus). In experiment 1, click level was held constant and click rate varied from 25 to 1250 Hz. For MLS, interstimulus intervals varied by a factor of 6 at each rate, while for I-RSA the interstimulus intervals varied by ± 0.5 ms regardless of rate. In experiment 2, stimulus rates ranged from 100 to 1000 Hz and click level varied from 105 to 135 dB re: 1 μPa. For experiment 1, MLS and I-RSA showed similar decreases in ABR peak amplitudes and increases in ABR peak latencies and interwave intervals with increasing rate. For experiment 2, there was an increase in peak latency and a decrease in peak amplitude with decreasing click level; however, the effects of click level were reduced at higher rates. The results indicate that the greater jitter for MLS compared to I-RSA does not substantially affect the dolphin ABR.

A simple method for EEG guided transcranial electrical stimulation without models

OBJECTIVE

There is longstanding interest in using EEG measurements to inform transcranial Electrical Stimulation (tES) but adoption is lacking because users need a simple and adaptable recipe. The conventional approach is to use anatomical head-models for both source localization (the EEG inverse problem) and current flow modeling (the tES forward model), but this approach is computationally demanding, requires an anatomical MRI, and strict assumptions about the target brain regions. We evaluate techniques whereby tES dose is derived from EEG without the need for an anatomical head model, target assumptions, difficult case-by-case conjecture, or many stimulation electrodes.

APPROACH

We developed a simple two-step approach to EEG-guided tES that based on the topography of the EEG: (1) selects locations to be used for stimulation; (2) determines current applied to each electrode. Each step is performed based solely on the EEG with no need for head models or source localization. Cortical dipoles represent idealized brain targets. EEG-guided tES strategies are verified using a finite element method simulation of the EEG generated by a dipole, oriented either tangential or radial to the scalp surface, and then simulating the tES-generated electric field produced by each model-free technique. These model-free approaches are compared to a 'gold standard' numerically optimized dose of tES that assumes perfect understanding of the dipole location and head anatomy. We vary the number of electrodes from a few to over three hundred, with focality or intensity as optimization criterion.

MAIN RESULTS

Model-free approaches evaluated include (1) voltage-to-voltage, (2) voltage-to-current; (3) Laplacian; and two Ad-Hoc techniques (4) dipole sink-to-sink; and (5) sink to concentric. Our results demonstrate that simple ad hoc approaches can achieve reasonable targeting for the case of a cortical dipole, remarkably with only 2-8 electrodes and no need for a model of the head.

SIGNIFICANCE

Our approach is verified directly only for a theoretically localized source, but may be potentially applied to an arbitrary EEG topography. For its simplicity and linearity, our recipe for model-free EEG guided tES lends itself to broad adoption and can be applied to static (tDCS), time-variant (e.g., tACS, tRNS, tPCS), or closed-loop tES.