Auditory cortex iEEG responses to pure tone stimuli in normal-hearing and deafened cochlear implant rats.

A, Schematic of 60-channel cortical surface electrode arrays, covering 3.25 x 3.25 mm on the surface of auditory cortex to record tone-evoked responses in normal hearing (NH) trained rats.

B, Some animals were bilaterally deafened and fitted with a unilateral cochlear implant (CI).

C, Behavioral training of subset of animals on auditory go/no-go frequency recognition task. Animals were first trained when normal-hearing (NH, N=4, d’: 2.9±0.1) and then re-trained to respond to cochlear implant stimulation (CI, N=3, d’: 0.9±0.2).

D, Examples of trial-averaged single-site iEEG responses (top, ERPs; middle, ERP spectrograms; bottom, HG). Clear transients were evoked by tones in normal-hearing animals (NH, left column) and in cochlear implant rats (CI, right column).

E, iEEG response magnitude was similar between normal-hearing (NH, ERP amplitude: 1.9±0.1 µV; HG amplitude: 0.6±0.1 a.u.) and cochlear-implant rats (CI, ERP amplitude: 2.0±0.2 µV, p=0.43 compared to normal-hearing ERPs, Student’s paired two-tailed t-test; HG amplitude: 0.6±0.1 a.u., p=0.33 compared to normal-hearing HG).

Tone-evoked and cochlear implant-evoked iEEG measurements are spatially organized.

A, F, Trial-averaged tone-evoked event-related potentials (ERPs, panel A) and high gamma (HG, panel F) were spatially restricted to regions along the iEEG recording array; these active regions shifted as a function of stimulus.

B, G, Spatial correlations of evoked response areas decreased monotonically as a function of increasing stimulus separation. Evoked iEEG responses as assessed by ERPs and HG were nonrandom for normal-hearing and implanted rats (p<10-8 compared to all shuffled responses), but were locally tonotopic in normal-hearing rats (p<10-4) but not implanted rats (ERPs: p=0.5, HG: p=0.7).

C, H, Determining the preferred stimulus of spatially-evoked activity revealed smooth gradients shifting from high-to-low to high tone frequencies or cochlear implant channels (ERPs, panel C; HG, panel H).

D, I, Local tonotopic gradients for each recording site plotted on a unit circle (blue) as a function of magnitude (strength of gradient) and angle (direction of gradient). The vectors across each recording site were averaged to produce an overall tonotopic vector (black) whose magnitude represents a metric of overall topography (ERPs, panel D; HG, panel I). The mean vector was plotted against the mean vector computed from n=1,000 shuffled maps (ERPs, panel D: normal-hearing z-score: 6.6, p<10-10; implanted z-score: 7.9, p<10-14; HG, panel I: normal-hearing z-score: 2.1, p=0.02; implanted z-score: 1.5, p=0.07).

E, J, Z-scored magnitudes of mean tonotopic vectors across animals (ERPs, panel E: normal-hearing mean z-score: 4.7±1.0, implanted mean z-score: 2.6±1.0, p=0.22, Student’s paired two-tailed t-test; HG, panel J: normal-hearing mean z-score: 3.0±0.9, implanted mean z-score: 2.1±0.9, p=0.96).

Trial-by-trial cochleotopic encoding is more variable in cochlear implant vs normal-hearing rats.

A, C, Trial-by-trial iEEG measurements of both tone-evoked (NH) and cochlear implant-evoked (CI) neural activity revealed peaks across time (left) and space (right) (ERP, A; HG, C)

B, D, Variability of iEEG measurements across trials (root mean square, rms) was consistently higher for cochlear implant-evoked compared to tone-evoked activity, for ERPs in panel B (top, ERP temporal NH mean rms: 1.6±0.3 over all animals, CI mean rms: 1.9±0.5, Student’s two-tailed paired t-test for animals monitored before and after deafening, p=0.12; bottom, ERP spatial NH mean rms: 0.32±0.01, CI mean rms: 0.33±0.02, p=0.05) and HG in panel D (top, HG temporal NH mean rms: 2.4±0.3 over all animals, CI mean rms: 3.3±0.9, p=0.18; bottom, HG spatial NH mean rms: 1.2±0.2, CI mean rms: 1.3±0.3, p=0.32).

Trial-by-trial iEEG measurements encode stimulus identity.

A, Trial-by-trial ERPs evoked by 1.4 kHz (left), 8 kHz (middle), and 32 kHz tones (right) are plotted for three recording sites (blue) reveal discernable evoked transients and spatially restricted patterns of evoked magnitudes (stimulus onsets depicted by dotted line).

B, Trial-by-trial iEEG measurements concatenated such that columns represent recording sites (spatial) and post-stimulus sampling (temporal), rows by stimulus trials. Data then reorganized by PCA and reduced to the 15 components according to magnitude of explained variance.

C, Classification matrices are plotted means across bootstrapped repeated (N=1,000) versions of LDA classifiers trained using 13 randomly selected trials for each stimulus, and classification predictions of single and remaining trials reveal significant prediction of stimulus identity (dashed line: chance-level error distance; solid lines: actual mean error distances).

D, Stimulus prediction probabilities plotted across animals (black: mean, grey: s.e.m.) and as function of either octaves (normal-hearing) or channels (cochlear implant) from actual stimulus, reflecting a gradient in which adjacent stimuli share more encoded features than other stimuli.

E, Decoder performance (mean error distance) across individual animals was somewhat worse for cochlear implant ERPs compared to tone-evoked ERPs (ERP mean error rate relative to chance for normal-hearing: 0.74±0.05, cochlear implant: 0.78±0.06; for animals assessed both first when normal-hearing and then after cochlear implantation p=0.02, Student’s two-tailed paired t-test), but not for HG (HG mean error rate relative to chance for normal-hearing: 0.90±0.03, cochlear implant: 0.91±0.04; for animals assessed both first when normal-hearing and then after cochlear implantation p=0.11).

Trial-by-trial decoding from either spatial or temporal aspects of iEEG signals.

A, Raw trial-by-trial ERPs for a single animal are plotted for 3 channels.

B-C, Trial-by-trial iEEG measurements are reduced into either spatial only (B) or temporal only (C) by extracting either the magnitude of evoked activity or averaging across all recording sites.

D, Re-training PCA-LDA decoders on spatial-only (middle) or temporal-only (right) does not abolish the encoding of stimulus identity.

E, Decoder errors (mean error distance re. chance) across animals are plotted for decoders trained on spatial+temporal, spatial-only, and temporal-only iEEG measurements. Mean error distance for spatial+temporal ERPs, normal-hearing: 0.74±0.04, implant: 0.86±0.06; for spatial+temporal HG, normal-hearing: 0.90±0.03, implant: 0.97±0.04. Mean error distance for spatial-only ERPs, normal-hearing: 0.77±0.04, implant: 0.77±0.08; for HG, normal-hearing: 0.89±0.04, implant: 0.87±0.05. Mean error distance for temporal-only ERPs, normal-hearing: 0.82±0.05, implant: 0.81±0.07; for HG, normal-hearing: 0.87±0.04, implant: 0.93±0.05.

TCA-based decoding of single-trial iEEG measurements.

A, Schematic of TCA with 3-dimensional tensors with orthogonal dimensions of spatial, temporal and trial factors.

B, TCA-reduced iEEG measurements of stimulus-evoked ERPs in normal-hearing rats are plotted across 15 components of spatial (top row), temporal (middle row) and trial factors (bottom row).

C, Classification matrices of mean decoder predictions across bootstrapped repeated (N=1,000) versions LDA classifiers trained using 13 randomly selected trials for each stimulus indicated that iEEG measurements encode stimulus identity from single stimulus presentations (dashed line: chance-level error distance; solid lines: actual mean error distances).

D, Stimulus prediction probabilities are plotted across animals (black: mean, grey: std. error) and as a function of either octaves or channels in away from actual stimulus, normal-hearing and implanted rats. Correct prediction probabilities are high across all conditions and prediction errors coalescing of predictions toward actual stimuli (orange: mean and std. error of decoder performance on shuffled trials; grey: mean error distance, dashed line: chance level error distance).

E, Decoder errors (mean error distance) across individual animals were slightly smaller for evoked iEEG measurements in normal-hearing compared to implanted rats (ERP mean error rate re. chance normal-hearing: 0.77±0.04, implanted: 0.80±0.06; HG mean error rate re. chance normal-hearing: 0.90±0.03, implanted: 0.88±0.05). Animals with both normal-hearing and implanted iEEG measurements only showed significant differences for ERP (top, Student’s paired t-test: p=0.01) but not HG (bottom, p=0.09).

TCA data reductions reveal latent spatial factors are topographically organized.

A, TCA-based analysis of evoked neural activity, divided into 15 unique components where each component is separated into orthogonal dimensions of spatial factors (top row), temporal factors (middle row) and trial factors (bottom row).

B, Spatial maps for tone-evoked or electrode-evoked responses recreated from TCA reduced data.

C, Reducing re-organized spatial factors revealed a TCA map (left) with topographical gradients in same location and direction as the tonotopic map reduced from raw measurements (right).

D, Spatial correlations of evoked response areas decreased monotonically as a function of increasing stimulus separation. TCA-reduced models of evoked iEEG responses as assessed by ERPs and HG were nonrandom for normal-hearing and implanted rats (p<10-8 compared to all shuffled responses). TCA-reduced models were locally tonotopic only for HG in normal-hearing rats (normal-hearing ERPs: p=0.18, normal-hearing HG: p<10-4, cochlear implant ERPs: p=0.43, cochlear implant HG: p=0.40).

E, Z-scored magnitudes of mean tonotopic vectors across animals were similar (ERPs, normal-hearing mean z-score: 6.2±3.6, cochlear implant mean z-score: 3.9±2.1, p=0.28, Student’s paired two-tailed t-test; HG, normal-hearing mean z-score: 2.9±1.1, implanted mean z-score: 2.5±0.7, p=0.88).

Lack of information transfer between acoustic and electrical stimulation representations in the same animals.

A, Evoked iEEG measurements from a normal-hearing rat reduced using TCA.

B, Evoked iEEG measurements in cochlear implant rat modeled using TCA, constraining the model with the spatial and temporal factors from the TCA model optimized on normal-hearing data, leaving only the trials as the optimizable variables.

C, Linear discriminant analysis classifiers are trained on the trial factors extracted from the TCA-reduced models of evoked data in normal-hearing conditions and then used to predict stimulus identity from the trial factors from the TCA-reduced models of implant-evoked measurements.

D, Classification matrices are plotted means of decoder predictions across bootstrapped repeated (N=1,000) versions of linear-discriminant analysis (LDA) classifiers reveal little information transfer (IT for animal 1; for ERPs: 1.1%, for HG: 0.3%).

E, Classification matrices of normal-hearing trained and implant tested decoders reveal little-to- no information transfer in three other animals with iEEG measurements of both tone-evoked and implant-evoked responses (IT for animal 2, ERPs: 5.9%, HG: 0.2%; for animal 3, ERPs: 1.0%, HG: 0.1%; for animal 4, ERPs: 0.9%, HG: 0.1%).