Viene presentato un protocollo per valutare l’evoluzione temporale dell’ingresso neurale verso stimoli ripetitivi esterni. Le registrazioni di stato fisso della stessa condizione sperimentale vengono acquisite e mediate nel dominio temporale. Le dinamiche dello stato fisso vengono analizzate tracciando l’ampiezza della risposta in funzione del tempo.
L’intrappolamento neurale si riferisce alla sincronizzazione dell’attività neurale con la periodicità degli stimoli sensoriali. Questa sincronizzazione definisce la generazione di risposte evocate allo stato costante (cioè oscillazioni nell’elettroencefalogramma che si bloccano in fase agli stimoli di guida). L’interpretazione classica dell’ampiezza delle risposte allo stato costante presuppone una risposta neurale stereotipata invariante nel tempo più fluttuazioni di fondo casuali, in modo tale che la media su ripetute presentazioni dello stimolo recuperi il risposta stereotipata. Questo approccio ignora la dinamica dello stato stazionario, come nel caso dell’adattamento suscitato da esposizioni prolungate allo stimolo. Per analizzare la dinamica delle risposte allo stato stazionario, si può presumere che l’evoluzione temporale dell’ampiezza di risposta sia la stessa in diversi cicli di stimolazione separati da pause sufficientemente lunghe. Sulla base di questo presupposto, viene presentato un metodo per caratterizzare l’evoluzione temporale delle risposte in stato stazionario. Un numero sufficientemente elevato di registrazioni viene acquisito in risposta alla stessa condizione sperimentale. Le esecuzioni sperimentali (registrazioni) sono mediate a colonne (ad esempio, le esecuzioni sono mediate, ma l’epoca all’interno delle registrazioni non viene mediata con i segmenti precedenti). La media a colonne consente l’analisi delle risposte allo stato fisso nelle registrazioni con rapporti segnale-rumore notevolmente elevati. Pertanto, il segnale mediato fornisce una rappresentazione accurata dell’evoluzione temporale della risposta allo stato costante, che può essere analizzata sia nei domini di tempo che in quelli di frequenza. In questo studio, viene fornita una descrizione dettagliata del metodo, utilizzando potenziali evocati visivamente allo stato costante come esempio di risposta. Vantaggi e avvertenze vengono valutati sulla base di un confronto con metodi a prova singola progettati per analizzare l’ingresso neurale.
Quando viene registrata dal cuoio capelluto, l’attività elettrica del cervello è osservata come cambiamenti continui e regolari nelle tensioni nel tempo. Questa attività elettrica è chiamata elettroencefalogramma (EEG) ed è stata descritta per la prima volta da Hans Berger alla fine degli anni venti del secolo scorso1. Successivi studi seminali hanno descritto l’EEG come una serie temporale composta, in cui si possono osservare diversi modelli ritmici o ripetitivi2,3,4. Al giorno d’oggi, l’EEG è tipicamente diviso in cinque bande di frequenza ben consolidate, delta, theta, alfa, beta e gamma, che sono associate al diverso processo sensoriale e cognitivo.
Per anni, lo studio delle oscillazioni cerebrali utilizzando EEG è stato limitato a due analisi dello spettro nell’attività in corso o ai cambiamenti nell’attività oscillatoria suscitati da eventi sensoriali non periodici. Negli ultimi decenni sono state implementate diverse metodologie per modulare le oscillazioni EEG in corso ed esplorare gli effetti di tali modulazioni sui processi percettivi e cognitivi, compresa la presentazione della stimolazione sensoriale ritmica per indurre l’intrappolamento neurale. Il termine rientro neurale si riferisce alla sincronizzazione dell’attività neurale con le proprietà periodiche degli stimoli sensoriali. Questo processo porta alla generazione di potenziali evocati a stato costante (cioè oscillazioni EEG bloccate alle proprietà periodiche degli stimoli di guida). I potenziali evocati da stato costante sono più comunemente provocati dalla stimolazione visiva, uditiva e vibrotilativa, usando sia stimoli transitori presentati a un ritmo costante o stimolazione continua modulata in ampiezza alla frequenza di interesse. Mentre i potenziali somatosensoriali con rimorsi costanti (SSSEP) sono registrati in risposta alla stimolazione tattile ripetitiva5,6, i potenziali sottoposti a residuo visivamente evocati (SSVEP) sono generalmente suscitati dal periodico presentazione di sfarfallii di luminanza, immagini e volti7,8. Le risposte uditiche a stato stabile (ASSR) sono solitamente generate da treni di stimoli acustici transitori o dalla presentazione continua di toni modulati in ampiezza9,10.
L’estrazione di potenziali evocati dallo stato stazionario EEG si basa essenzialmente sulla media di ePochs EEG successivamente acquisita, bloccate nel tempo allo stimolo11. A causa della periodicità delle risposte, possono essere analizzate sia nei domini di tempo che in quelli di frequenza. Dopo la trasformazione frequenza-dominio, la risposta sensoriale viene osservata come picchi di ampiezza alla velocità di presentazione o alla frequenza di modulazione degli stimoli esterni e le armoniche corrispondenti. Queste procedure (media tra domini temporali e successiva trasformazione frequenza-dominio) sono state essenziali per lo sviluppo di un test dell’udito basato sul rilevamento di metodi ASSR con scopi clinici12,13,14 ,15,16.
Inoltre, la media classica del dominio temporale delle epoche EEG è stata estremamente utile per analizzare i processi fisiologici come la generazione e l’estinzione di SSVEP17,18. Presentando treni consecutivi di luci sfarfallio e media di epoche successive all’interno di una registrazione, Wacker et al.19 ha osservato che l’indice di blocco di fase del SSVEP è aumentato rapidamente durante i primi 400 ms di stimolazione ed è rimasto alto in seguito . Hanno anche riferito che è stato stabilito un robusto intrappolamento visivo tra 700-1 100 ms dopo l’insorgenza di stimoli. Un certo grado di intrappolamento è rimasto efficace dopo l’offset del treno di stimolazione, che è durato circa tre periodi della risposta oscillatoria17,19. Tali comportamenti sono stati interpretati come l’effetto coinvolgente/disinnesto delle oscillazioni osservate, che è una conseguenza dell’elaborazione delle informazioni non lineari nel sistema visivo umano17. In alternativa, è noto che in determinate condizioni sperimentali, la stimolazione dello sfarfallio può suscitare risposte all’inizio e off-risposta alla fine dei treni di stimolazione invece di intrappolamento neurale18.
L’ipotesi principale di epoche EEG acquisite consecutivamente nella media è che il segnale EEG rappresenta una combinazione lineare della risposta sensoriale e del rumore di fondo20. Inoltre, si presume che l’ampiezza, la frequenza e la fase della risposta oscillatoria siano stazionarie, mentre il rumore di fondo è considerato un’attività casuale. Tuttavia, nei casi in cui questa ipotesi non è soddisfatta, l’ampiezza di risposta calcolata dopo diverse epoche non corrisponde necessariamente all’ampiezza istantanea del potenziale evocato.
È stato recentemente riportato che l’ASSR generato nel tronco encefalico dei ratti si adatta alla continua presentazione di toni modulati con ampiezza (cioè, l’ampiezza di risposta esponenzialmente nel tempo)21,22. L’adattamento è stato interpretato come un meccanismo neurale che riflette la perdita di novità di uno stimolo sensoriale monotonamente ripetitivo, aumentando la sensibilità alle fluttuazioni rilevanti nell’ambiente acustico23,24. Nel percorso uditivo, l’adattamento può migliorare la comprensione del linguaggio in ambienti rumorosi. Inoltre, questo processo può essere una parte dei meccanismi esistenti per monitorare il feedback uditivo della propria voce per controllare la produzione vocale.
Analizzando l’evoluzione temporale dell’ASSR a 40 Hz negli esseri umani, Van Eeckhoutte et al.25 ha osservato una significativa ma piccola diminuzione dell’ampiezza di risposta nel tempo (circa -0.0002 V/s basata sull’analisi di gruppo, quando assumendo una diminuzione lineare nel tempo). Di conseguenza, questi autori hanno concluso che l’ASSR a 40 Hz nell’uomo non si adatta alla stimolazione. Nell’uomo, sono stati osservati comportamenti non stazionari durante l’analisi della stabilità dell’SSVEP26. Questi autori hanno osservato che l’ampiezza della frequenza fondamentale e la seconda armonica dell’SSVEP erano stazionarie solo nel 30% e nel 66,7% dei soggetti testati, rispettivamente. Le fasi di entrambi i componenti di frequenza SSVEP, anche se relativamente stabili nel tempo, hanno mostrato piccole derive26.
Pertanto, anche se la media classica del dominio temporale delle epoche successivamente acquisite consente di esplorare le proprietà stazionarie dell’interioramento neurale, questa metodologia deve essere rivista quando la dinamica a lungo termine dell’intrappolamento è al centro del ricerca, o quando la media delle dinamiche a breve termine è corrotta dall’insorgenza di dinamiche a lungo termine. Per caratterizzare i comportamenti non stazionari delle risposte allo stato stazionario, la risposta evocata calcolata in un determinato intervallo di tempo non deve essere compromessa da quelle calcolate nei segmenti EEG precedenti. In altre parole, il potenziale evocato dovrebbe essere estratto dal rumore di fondo senza che le epoche siano mediate nel dominio temporale con i segmenti EEG precedenti.
In questo studio, viene presentato un metodo per valutare la dinamica dell’inserimento neurale. Le risposte dello stato fisso vengono registrate ripetutamente in risposta alla stessa stimolazione, in cui le registrazioni consecutive vengono interfogliate da un intervallo di riposo di tre volte la lunghezza della corsa sperimentale. Considerando che se l’evoluzione temporale della risposta fisiologica è la stessa in diverse esecuzioni sperimentali indipendenti (registrazioni indipendenti), le registrazioni sono in media a colonne. In altre parole, vengono mediate le epoche corrispondenti alla stessa posizione nelle diverse registrazioni, senza una media delle epoche all’interno di una registrazione. Di conseguenza, l’ampiezza di risposta calcolata a qualsiasi intervallo di stimolazione corrisponderà all’ampiezza istantanea del potenziale evocato. Le risposte sensoriali possono essere analizzate nel periodo temporale o trasformate nel dominio di frequenza, a seconda dello scopo dell’esperimento. In ogni caso, le ampiezza possono essere tracciate in funzione del tempo per analizzare l’evoluzione del tempo della risposta allo stato costante. Generazione ed estinzione dei potenziali evocati allo stato costante possono essere valutati limitando l’analisi alla prima e all’ultima epoche delle registrazioni.
La dinamica dell’ingresso neurale può essere analizzata utilizzando altri approcci, come il filtraggio a banda stretta di misurazioni a singola prova intorno alla frequenza di interesse e il calcolo dell’inviluppo del segnale di alimentazione utilizzando il filtraggio a passa basso25 e il Trasformazione di Hilbert27. Rispetto a queste metodologie, la media a colonne delle epoche consente di calcolare i parametri a stato costante in base ai segnali con il più alto rapporto segnale-rumore (SNR). Recentemente, il filtraggio Kalman è emerso come una tecnica promettente per la stima di ampiezza ASSR 40-Hz28,29,30. L’implementazione del filtraggio Kalman può migliorare il rilevamento di risposte allo stato costante più vicino alla soglia elettrofisiologica e ridurre il tempo del test dell’udito29. Inoltre, non è necessario assumere risposte fisse quando si utilizza un approccio di filtraggio Kalman per stimare l’ampiezza dell’ASSR30. Tuttavia, solo uno studio ha analizzato l’evoluzione temporale degli ASSR utilizzando Kalman filtering25. La conclusione dello studio è che l’ampiezza ASSR a 40 Hz è stabile nell’intervallo di stimolazione. Pertanto, il filtraggio Kalman deve essere testato in condizioni in cui l’ASRR non è stazionario.
Anche se richiede molto tempo, il metodo di media per colonna è privo di modelli e non richiede valori di inizializzazione e/o definizioni a priori del comportamento del rumore. Inoltre, poiché non comporta tempi di convergenza, la media a colonne può fornire una rappresentazione più affidabile dell’insorgenza dell’insorgenza dell’insorgenza del rientro neurale. Pertanto, i risultati ottenuti con il metodo di media a colonne può essere considerato come la verità di base per analizzare la dinamica dell’intrappolamento neurale utilizzando il filtraggio Kalman.
Questa descrizione del protocollo si basa su un esempio di SSVEP. Tuttavia, è importante notare che il metodo qui presentato è indipendente dalla modalità, in modo che possa essere utilizzato anche per analizzare l’evoluzione temporale di SSSEP e ASSR.
Questo lavoro descrive una procedura sperimentale per analizzare la dinamica delle risposte cerebrali oscillatorie. Tale metodologia consiste nell’acquisire un numero sufficiente di esecuzioni sperimentali indipendenti della stessa condizione sperimentale e di epoche di media del settore temporale corrispondenti allo stesso intervallo temporale nelle diverse registrazioni (colonne-wise in media in Figura 1B). L’ampiezza calcolata nei dati medi rappresenta l’ampiezza istantanea della risposta oscillatoria. Tracciare queste ampiezza in funzione del tempo (o del numero di colonne nel set di dati) consente di analizzare l’evoluzione temporale della risposta oscillatoria bloccata dalla stimolazione. Questa metodologia è una modifica di quella proposta da Ritter et al.23 per analizzare l’adattamento dei potenziali corticali transitori evocati. Il metodo è stato utilizzato per analizzare la dinamica dei potenziali uditivi evocati sia in casi umani24 che in modelli animali20,21.
Da un punto di vista metodologico, la combinazione di parametri utilizzati per suscitare la risposta allo stato costante e quelli implementati per estrarre la risposta neurale dal rumore di fondo è fondamentale per analizzare l’evoluzione del tempo dei potenziali evocati allo stato costante 22. La lunghezza dello stimolo utilizzata nell’esperimento qui presentato (40 s) è stata selezionata sulla base dei risultati ottenuti in uno studio pilota. Questa lunghezza di stimolo è stata sufficiente per analizzare l’adattamento dell’ASSR generato nel tronco encefalico del ratto21,22. Inoltre, la lunghezza dello stimolo dovrebbe superare il tempo in cui viene raggiunto il potere anotoso della banda istantanea degli SSVIP (Figura 1 in Labecki et al.26). Tuttavia, il potere anodo anotico della banda istantanea degli SSVIP può essere raggiunto oltre gli anni ’60 in alcuni casi (Figura 2 in Labecki et al.26). Pertanto, si raccomanda di eseguire uno studio pilota di piccoli campioni per definire la lunghezza dello stimolo della stimolazione. In caso contrario, si raccomanda una lunghezza di stimolo più lunga di 90 s per ottenere una rappresentazione completa dell’evoluzione temporale della risposta. L’utilizzo di pause adeguatamente lunghe tra registrazioni consecutive implica considerare le esecuzioni sperimentali consecutive come statisticamente indipendenti (cioè diverse misure indipendenti della stessa variabile). Per quanto ne sappiamo, non sono stati effettuati esperimenti per analizzare la pausa ottimale tra le esecuzioni (pausa minima necessaria per effettuare esecuzioni indipendenti l’una dall’altra). Il criterio dell’uso delle pause di almeno 3 volte superiore alla lunghezza dello stimolo è sufficientemente conservativo da garantire che la risposta allo stato costante registrata in una determinata esecuzione non sia influenzata dalla stimolazione precedente.
Recentemente, stimoli alternati (condizioni sperimentali) è stato proposto come una scelta per ridurre la pausa tra le esecuzioni sperimentali, evitando un effetto di adattamento supplementare25. Allo stesso modo, il numero di esecuzioni sperimentali (30) implementate in questo protocollo sperimentale è conservativo, poiché la RNL asintotica e il pSNR sono in genere raggiunti dopo una media di 20 esecuzioni sperimentali, approssimativamente. Quando gli stimoli rientrano nella regione medio-alta della gamma dinamica della risposta (alti livelli di sensazione), è probabile che siano necessari un numero inferiore di corse per analizzare le dinamiche della risposta evocata. Tuttavia, nei casi in cui vengono testate diverse condizioni sperimentali, avere lo stesso numero di esecuzioni sperimentali è fondamentale per fare confronti tra le condizioni (cioè diversi livelli di sensazione).
Oltre alla media a colonne delle epoche, la dinamica dei potenziali evocati oscillatori è stata analizzata filtrando le misurazioni a prova singola in una banda di frequenza stretta intorno alla frequenza di interesse e calcolando l’involucro della potenza segnale utilizzando il filtro passa-basso26. Allo stesso modo, è stata attuata un’analisi di prova singola per caratterizzare il periodo di transizione che precede la regione stabile di SSVEP48e i cambiamenti di ampiezza e fase dell’SSVEP durante la regione stabile della risposta49. Mentre le singole analisi di prova consentono la discriminazione delle fluttuazioni relativamente rapide dell’ampiezza di risposta, i progetti sperimentali per analizzare la risposta media in blocchi separati con un determinato intervallo interblocco rappresentano solo variazioni a lungo termine ampiezza del potenziale evocato50,51. La media a colonne delle epoche si distingue tra queste due opzioni. La conversione del segnale medio al dominio di frequenza utilizzando il FFT implica analizzare la dinamica della risposta con una risoluzione uguale alla lunghezza dell’epoca. Nell’esempio qui presentato, l’SSVEP è stato segnalato ogni 4 s. Sebbene 4 s di risoluzione sia sufficiente a descrivere dinamiche che si verificano a intervalli di tempo superiori a decine di secondi, come quella dell’SSVEP26, le epoche parzialmente sovrapposte nelle registrazioni originali permette di descrivere l’evoluzione temporale del risposta allo stato costante in modo più raffinato25.
Le dinamiche delle risposte allo stato stazionario ottenute dopo la media a colonne delle epoche rappresentano principalmente l’evoluzione dell’attività oscillatoria sincronizzata tra i segmenti EEG mediati (quelli che sopravvivono alla media). Pertanto, un problema importante per quanto riguarda la fattibilità della metodologia è la possibile attenuazione delle ampiezza di risposta a causa delle variazioni nella fase delle oscillazioni neurali da una corsa sperimentale indipendente a un’altra (cioè tra le registrazioni). Questo argomento deve essere affrontato sperimentalmente. Tuttavia, le prove indicano che la fase delle risposte oscillatorie cerebrali è meno variabile del previsto. Infatti, diversi studi hanno segnalato una regolarità nella fase prevista dell’ASSR umano 80 Hz47,48,49. Quando le latenze sono stimate in base alla fase dell’attività oscillatoria, è stato osservato l’effetto prevedibile dell’intensità e la frequenza portante degli stimoli acustici sulla latenza delle risposte uditive (cioè la diminuzione della latenza l’intensità e l’aumento della frequenza dei portini)52,53,54. Inoltre, sono stati osservati cambiamenti maturazionali tipici dell’ampiezza e dell’asimmetria da sinistra a destra nei livelli uditivi, quando le latenze sono stimate dalla fase dell’ASSR47,55,56 , 57 del sistema , 58. Nel descrivere l’evoluzione temporale di SSVEP mediante analisi a prova singola, Labecki et al.26 ha osservato che, sebbene la variabilità inter-prova delle ampiezza di risposta all’interno dello stesso argomento fosse considerevolmente elevata, la variabilità del fase è stata significativamente meno pronunciata.
Sulla base delle loro osservazioni, Labecki et al.26 suggerirono che si dovrebbe calcolare una media di un minimo di 50 prove per ottenere una stima affidabile dell’involucro di potenza media della risposta. Questi risultati indicano che, anche quando l’ampiezza della risposta viene calcolata in singole prove, è necessaria una media (di buste in quel caso) per riportare risultati affidabili. Inoltre, la variabilità inter-prova nell’ampiezza di SSVEP riportata da Labecki et al.26 suggerisce che il calcolo di questo parametro in singole prove può essere fortemente influenzato dal rumore di fondo. Considerando l’evoluzione del rapporto segnale-rumore presentato nella Figura 2, il calcolo della risposta nel segnale medio invece di singole prove riduce significativamente il numero di segmenti EEG necessari per ottenere affidabili Misure. Inoltre, la bassa variabilità in fase ottenuta da Labecki et al.26 sostiene l’idea che la media a colonne delle epoche qui presentata sia una valida procedura per calcolare le dinamiche dei potenziali evocati oscillatori.
La media dei dati a diversi livelli porta a un’interpretazione diversa dei risultati. Per quanto riguarda i potenziali evocati oscillatori, calcolare l’ampiezza di risposta dopo la media del dominio temporale di esecuzioni indipendenti implica l’analisi solo delle oscillazioni bloccate nel tempo (cioè quelle che sopravvivono alla media). Questa procedura può filtrare le informazioni rilevanti relative alle dinamiche della risposta nelle singole prove. Tuttavia, consente un rapporto segnale-rumore sufficientemente elevato delle misurazioni. Questo aspetto potrebbe essere significativo quando le risposte sono vicine alla soglia elettrofisiologica, una condizione in cui il rilevamento dell’ingresso può essere compromesso a causa del basso rapporto segnale-rumore della misurazione.
The authors have nothing to disclose.
Gli autori riconoscono con gratitudine Lucàa Epeda, Grace A. Whitaker e Nicolas Nieto per il loro contributo alla produzione video. Questo lavoro è stato sostenuto in parte dai programmi CONICYT BASAL FB0008, MEC 80170124 e dalla borsa di dottorato 21171741, nonché dal National Institute on Deafness and Other Communication Disorders dei National Institutes of Health con il numero di premio P50DC015446. Il contenuto è di esclusiva responsabilità degli autori e non rappresenta necessariamente le opinioni ufficiali dei National Institutes of Health.
Active electrodes | Biosemi | P32-1020-32ACMS (ABC) | for channels 1-32 |
Active electrodes | Biosemi | P32-1020-32A (ABC) | for channels 33-64 |
Active electrodes | Biosemi | 8 x TP FLAT | external electrodes |
Active-Two adquisition system | Biosemi | version 7.0 | EEG adquisition system |
alcohol | Salcobrand | Code: 3309011 | for cleaning the scalp |
Electrode cap 64 channels | Biosemi | CAP MS xx yy | cap |
Electrode cap 64 channels | Biosemi | CAP ML xx yy | cap |
gel | Biosemi | SIGNA BOX12 | conductive gel |
Laptop | Asus | Core i7 1TB DD + 128GB SSD 8GB RAM | computer for stimulation |
Laptop | Asus | Core i7 1TB DD + 128GB SSD 8GB RAM | computer for recording |
LED screen | in-house production | – | The screen consists of four light-emitting diodes (LEDs) situated on the center of a 50×50 cm black screen, as vertexes of a square of 5×5 cm |
sterile gauze | Salcobrand | Code: 8730277 | for cleaning the scalp |