Qualità dei dati e fattibilità dell'acquisizione
L'applicazione dei criteri predefiniti di qualità dei dati ha portato all'esclusione di canali e partecipanti con qualità del segnale insufficiente. I risultati positivi si caratterizzano per un'elevata conservazione dei canali in tutte le condizioni, registrazioni emodinamiche stabili e modelli di connettività interpretabili che permettono una stima affidabile delle caratteristiche della rete. I risultati scadenti o subottimali includono bassa qualità del segnale, artefatti dovuti al movimento eccessivi, accoppiamento inadeguato tra optodi e cuoio capelluto, scarsa conservazione dei canali o stime di rete inconsistenti che limitano l'interpretazione. Come illustrato in Figura 3 e in Figura 4, il dataset finale dopo la valutazione della qualità includeva 19 partecipanti con disturbo depressivo maggiore (MDD; età media: 31,21 anni, DS: 8,99; 10 femmine) e 19 controlli sani (HC; età media: 33,5 anni, DS: 13,51; 7 femmine). Sette partecipanti (3 MDD e 4 HC) sono stati esclusi a causa della qualità insufficiente dei dati, principalmente legata ad artefatti dovuti al movimento o a un accoppiamento inadeguato tra optodi e cuoio capelluto.
Le mappe di qualità per canale hanno ulteriormente dimostrato distribuzioni spaziali eterogenee della qualità del segnale tra i partecipanti, con regioni localizzate di qualità ridotta in alcuni individui. Questi risultati evidenziano l'importanza di applicare soglie predeterminate basate sull'indice di accoppiamento del cuoio capelluto e sul picco dello spettro di potenza prima di ulteriori analisi. Queste procedure di controllo qualità garantiscono che le misure della rete successive siano calcolate a partire da dati affidabili e rappresentino un componente critico del protocollo. Gli output rappresentativi mostrati qui dipendono direttamente dal completamento con successo della verifica della qualità dei dati, del preelaboramento standardizzato, dell'estrazione delle serie temporali basata su ROI e dell'applicazione della procedura di selezione predeterminata basata sulla somma dei ranghi.
Pre-elaborazione e costruzione della rete
Il disegno sperimentale e il flusso di lavoro di preelaborazione sono mostrati in Figura 1 e Figura 2. Dopo il passaggio di controllo qualità (Figura 2), i segnali sono stati correttamente convertiti in concentrazione di emoglobina e corretti per artefatti di movimento e rumore fisiologico sistemico. Il corretto completamento di questi passaggi di preelaborazione è indicato dalla conservazione di canali di alta qualità, da serie temporali regionali stabili e da stime di connettività interpretabili, mentre esiti scadenti includono perdita eccessiva del segnale, artefatti non risolti o misure di rete non affidabili. L'estrazione delle serie temporali regionali dai nodi della rete della modalità predefinita ha permesso la costruzione di matrici di connettività funzionale per ogni condizione. Questi risultati rappresentativi dipendono direttamente dall'esecuzione corretta del controllo qualità dei dati17, del preelaborazione standardizzata18, della definizione delle ROI e dell'applicazione della procedura di classificazione predeterminata. Le matrici risultanti di connettività funzionale della rete della modalità predefinita costituiscono la base per confronti successivi basati sulla distanza.
Esiti rappresentativi a livello di rete
Figura 5A mostra un esempio del profilo di connettività funzionale della rete della modalità predefinita (DMN) specifico per condizione, relativo a un partecipante rappresentativo affetto da disturbo depressivo maggiore (MDD), affiancato al modello di riferimento della rete DMN derivato da soggetti sani. Queste matrici permettono una visualizzazione diretta dei modelli di connettività tra regioni in diverse condizioni e aiutano a distinguere un output del protocollo di successo da stime di rete meno interpretabili o di scarsa qualità. L'ispezione visiva rivela variazioni dipendenti dalla condizione nell'organizzazione della rete rispetto al modello, illustrando come la frequenza della stimolazione possa modulare la struttura della rete DMN all'interno di una singola sessione.
Le misure di distanza (Figura 5B) quantificano queste differenze attraverso diverse metriche (connettività, autovalori, autovettori e controllabilità media). Per ciascuna metrica, le condizioni possono essere ordinate in base alla loro vicinanza al modello di riferimento, con valori di distanza più bassi che indicano una somiglianza maggiore. In particolare, l'ordinamento fornito dalle diverse metriche di distanza varia a seconda delle condizioni, riflettendo le informazioni complementari catturate da ciascuna metrica sulle rappresentazioni della rete.
Selezione della frequenza basata sul consenso
Per integrare le informazioni provenienti da quattro metriche, i ranghi sono stati aggregati all'interno di ogni partecipante (Figura 5C). Una chiara separazione dei punteggi di rango cumulativo tra le condizioni indica una differenziazione efficace delle risposte della rete specifiche per frequenza, mentre punteggi di rango molto simili tra le condizioni suggeriscono una sensibilità limitata del protocollo o una scarsa separazione tra gli effetti della stimolazione. La condizione con il punteggio di rango cumulativo più basso è stata interpretata come la frequenza di stimolazione rappresentativa candidata per quel partecipante. Tabella 1 riassume queste selezioni tra i partecipanti, insieme ai corrispondenti valori di guadagno di stimolazione e al grado di accordo tra le metriche quantificato tramite il coefficiente W di Kendall.
Nel gruppo MDD, 12 partecipanti su 19 hanno mostrato il rango cumulativo più basso per una frequenza di stimolazione attiva piuttosto che per la condizione placebo. La frequenza di stimolazione selezionata è variata tra i partecipanti, come riassunto nella Tabella 1. L'accordo tra le metriche, quantificato utilizzando il coefficiente W di Kendall, differiva anche tra i singoli individui, riflettendo una variabilità nella coerenza dei ranghi. Valori più elevati del coefficiente W di Kendall indicano una maggiore concordanza tra le metriche e una maggiore affidabilità nel fatto che diverse misure della rete sostengano lo stesso ordinamento delle condizioni, mentre valori più bassi indicano un accordo più debole e una maggiore eterogeneità tra le metriche. Pertanto, la selezione finale deve essere interpretata come un'integrazione guidata dai dati di caratteristiche di rete complementari, piuttosto che come una decisione basata su una singola metrica.
Questi risultati rappresentativi e basati sulla fattibilità illustrano i tipi di output a livello di rete generati dal protocollo e dimostrano che il flusso di lavoro può essere implementato per confronti intra-soggetto tra diverse condizioni di stimolazione. I risultati hanno lo scopo di illustrare l'acquisizione del segnale, il pre-elaborazione, l'estrazione delle caratteristiche e la mappatura dei parametri basata sui ranghi, piuttosto che costituire una prova di beneficio terapeutico, efficacia clinica o validazione delle prestazioni, come mostrato in Figure 1–5 e Tabella 1.

Figura 1. Progettazione sperimentale e configurazione della sonda fNIRS. (A) Progettazione sperimentale within-subject. Ogni partecipante ha completato tutte le condizioni all'interno di una singola sessione, inclusi lo stato di riposo basale (RS0), la stimolazione sham, quattro frequenze di stimolazione attiva (2, 10, 25 e 40 Hz) e lo stato di riposo post-stimolazione (RS1). L'ordine della stimolazione attiva è stato randomizzato tra i partecipanti. Ogni blocco è durato 6 min e separato da intervalli di circa 2 min senza stimolazione. I dati di spettroscopia funzionale nel vicino infrarosso (fNIRS) sono stati registrati continuamente (circa 54 min totali) e suddivisi per condizione ai fini dell'analisi. (B) Disposizione della sonda fNIRS che mostra 32 sorgenti (rosse) e 30 detector (blu), formanti 89 canali a lunga separazione e 8 canali a breve separazione (evidenziati). I canali a breve separazione sono stati utilizzati per rimuovere il rumore fisiologico superficiale. (C) Profilo di sensibilità corticale della disposizione della sonda rispetto alla rete della modalità predefinita (DMN), basato sull'atlante Automated Anatomical Labeling (AAL), che illustra la sensibilità della misurazione dipendente dalla profondità. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 2. Pipeline di preelaborazione ed estrazione delle caratteristiche per l'analisi di rete basata su fNIRS.
I dati grezzi sull'intensità della luce sono stati suddivisi in blocchi specifici per condizione (RS0, placebo, frequenze di stimolazione e RS1). La qualità dei dati è stata valutata utilizzando l'indice di accoppiamento cutaneo (SCI > 0,8) e il picco dello spettro di potenza (PSP > 0,1), e sono stati mantenuti solo i canali di alta qualità. I segnali sono stati convertiti in densità ottica, seguiti dalla correzione del movimento mediante riparazione della distribuzione della derivata temporale (TDDR) e filtraggio wavelet. Sono stati calcolati i cambiamenti di concentrazione dell'emoglobina, ed è stata applicata una regressione con separazione breve per ridurre il rumore fisiologico sistemico. Per l'analisi della rete, le serie temporali sono state estratte da regioni di interesse (ROI) all'interno della rete della modalità predefinita (DMN), inclusa la corteccia prefrontale mediale (mPFC), il precuneo (Prec), il lobulo parietale inferiore (IPL) e il lobo temporale mediale (MTL). È stata calcolata la connettività funzionale tra le ROI. Le caratteristiche derivate hanno incluso valori di connettività, autovalori, autovettori e controllabilità media. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3. Valutazione della qualità dei dati nei partecipanti con disturbo depressivo maggiore (MDD).
(A) Numero di canali che soddisfano i criteri di qualità durante la condizione basale (RS0) per ciascun partecipante con MDD. Le barre indicano tutti i canali (grigio), i canali a lunga separazione (blu) e i canali a breve separazione (arancione) che soddisfano i criteri di qualità (≥70% del tempo di registrazione con SCI > 0,8 e PSP > 0,1). Le linee tratteggiate indicano il numero massimo di canali a lunga separazione disponibili (89), di canali a breve separazione (8) e la soglia per una qualità accettabile dei dati a livello individuale (≥70% del totale dei canali conservati). I partecipanti al di sotto della soglia sono evidenziati (D13, D18, D19). (B) Qualità dei dati per ciascun canale attraverso i partecipanti. La mappa termica mostra la percentuale di tempo di registrazione che soddisfa i criteri di qualità per ciascun canale e partecipante durante RS0. Un’ombreggiatura più scura indica una qualità inferiore dei dati. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 4. Valutazione della qualità dei dati nei partecipanti controlli sani. (A) Numero di canali che soddisfano i criteri di qualità durante la condizione basale (RS0) per ciascun partecipante controllo sano. Le barre indicano tutti i canali (grigio), i canali a lunga separazione (blu) e i canali a breve separazione (arancione) che soddisfano i criteri di qualità (≥70% del tempo di registrazione con SCI > 0,8 e PSP > 0,1). Le linee tratteggiate indicano il numero massimo di canali disponibili e la soglia per una qualità accettabile dei dati a livello di partecipante (≥ 70% dei canali totali mantenuti). I partecipanti al di sotto della soglia sono evidenziati (HC13, HC14, HC15, HC18). (B) Qualità dei dati per canale attraverso i partecipanti. La mappa termica mostra la percentuale di tempo di registrazione che soddisfa i criteri di qualità per ciascun canale e partecipante durante RS0. Un’ombreggiatura più scura indica una qualità inferiore dei dati. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 5. Selezione rappresentativa basata sulla distanza della frequenza di stimolazione in un singolo partecipante.
(A) Matrici di connettività funzionale che mostrano le interazioni a coppie tra quattro regioni della rete della modalità predefinita per un partecipante rappresentativo con disturbo depressivo maggiore (D11) in diverse condizioni (RS0, sham e frequenze di stimolazione). Le regioni includono la corteccia prefrontale mediale (mPFC), il lobulo parietale inferiore (IPL), il precuneo (Prec) e il lobo temporale mediale (MTL). La matrice centrale rappresenta il modello di controllo sano derivato dai dati di livello di gruppo della condizione RS0. (B) Misure di distanza tra ciascuna condizione e il modello di controllo sano calcolate secondo quattro metriche: distanza di connettività (D_conn), distanza degli autovalori (D_eig), distanza degli autovettori (D_evec) e distanza della controllabilità media (D_AC). Le barre rappresentano i valori di distanza; valori più bassi indicano una maggiore somiglianza con il modello. Vengono mostrati i ranghi per ciascuna condizione (rango 1 = distanza minore). (C) Rango di consenso ottenuto sommando i ranghi per ciascuna condizione su tutte le metriche. La condizione con il rango totale più basso è stata selezionata come frequenza di stimolazione candidata per questo partecipante. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| ID partecipante | Connettività | Autovalore | Autovettore | Controllabilità media | Consenso | Kendall-W |
| D01 | 40HZ (G=1.00) | 40HZ (G=1.00) | 40HZ (G=1.00) | 40HZ (G=0.29) | 40HZ | 0.85 |
| D02 | Sham (G=0.18) | Sham (G=0.00) | Sham (G=0.11) | Sham (G=-0.09) | Sham | 0.89 |
| D03 | 10HZ (G=0.20) | Sham (G=0.05) | 2HZ (G=0.56) | Sham (G=0.63) | 40HZ | 0.34 |
| D07 | Sham (G=0.29) | 25HZ (G=0.06) | 2HZ (G=0.33) | 25HZ (G=0.45) | 25HZ | 0.36 |
| D10 | Sham (G=1.00) | Sham (G=1.00) | Sham (G=1.00) | Sham (G=1.00) | Sham | 0.59 |
| D11 | 10HZ (G=0.71) | 10HZ (G=0.66) | 2HZ (G=0.73) | 40HZ (G=1.00) | 10HZ | 0.34 |
| D15 | 40HZ (G=0.58) | 40HZ (G=0.53) | 40HZ (G=0.21) | 40HZ (G=-0.27) | 40HZ | 0.96 |
| D16 | 25HZ (G=1.00) | 10HZ (G=1.00) | 25HZ (G=0.90) | 25HZ (G=0.40) | 25HZ | 0.64 |
| D17 | 2HZ (G=0.51) | 2HZ (G=0.75) | 25HZ (G=0.73) | 10HZ (G=0.82) | 10HZ | 0.21 |
| D20 | 25HZ (G=-0.08) | 25HZ (G=0.04) | 25HZ (G=-0.27) | Sham (G=0.31) | 25HZ | 0.84 |
| D21 | Sham (G=0.37) | 2HZ (G=0.08) | Sham (G=0.27) | 40HZ (G=-0.28) | Sham | 0.41 |
| D22 | 2HZ (G=-0.06) | 10HZ (G=0.19) | 2HZ (G=0.75) | 2HZ (G=1.00) | 2HZ | 0.47 |
| D23 | Sham (G=0.53) | 25HZ (G=-0.14) | Sham (G=0.71) | 40HZ (G=0.23) | Sham | 0.31 |
| D24 | 25HZ (G=0.59) | 25HZ (G=0.18) | 25HZ (G=0.18) | Sham (G=0.34) | 25HZ | 0.5 |
| D25 | 25HZ (G=1.00) | Sham (G=-0.01) | 10HZ (G=1.00) | Sham (G=0.19) | 25HZ | 0.16 |
| D26 | 10HZ (G=0.52) | Sham (G=0.01) | 10HZ (G=1.00) | 10HZ (G=0.71) | 10HZ | 0.93 |
| D27 | 10HZ (G=-0.11) | 10HZ (G=0.20) | 2HZ (G=-0.20) | 40HZ (G=1.00) | 2HZ | 0.05 |
| D28 | 10HZ (G=0.51) | 10HZ (G=0.04) | 10HZ (G=1.00) | 10HZ (G=-0.02) | 10HZ | 0.81 |
| D29 | 10HZ (G=0.81) | 10HZ (G=0.54) | Sham (G=1.00) | 10HZ (G=1.00) | 10HZ | 0.85 |
| Guadagno | 0.50, 0.52 | 0.32, 0.18 | 0.58, 0.73 | 0.46, 0.40 | | 0.55, 0.50 |
| (Media, Mediana) |
Tabella 1: Selezione della frequenza specifica per partecipante in base alle metriche di distanza.
Per ogni partecipante, viene riportata la condizione di stimolazione associata alla distanza minima rispetto al modello del gruppo di controllo sano per quattro metriche basate su reti: distanza di connettività (D_conn), distanza dell'autovalore (D_eig), distanza dell'autovettore (D_evec) e distanza della controllabilità media (D_AC). I valori tra parentesi (G) indicano il guadagno di stimolazione. La colonna Consensus mostra la frequenza di stimolazione proposta in base all'aggregazione dei ranghi ottenuti dalle quattro metriche. Il coefficiente W di Kendall quantifica l'accordo tra le metriche all'interno di ciascun partecipante. L'ultima riga riassume i valori medi e mediani del guadagno di stimolazione ottenuti tra i partecipanti per ciascuna metrica, insieme ai corrispondenti valori di W di Kendall.