$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In termini di risultati del reclutamento, i partecipanti sono stati reclutati principalmente tramite l'invio di lettere di reclutamento e telefonate di follow-up in base alle normative delineate dell'Atlanta VA Healthcare System. Il team di studio ha reclutato un totale di 50 partecipanti, dimostrando l'efficacia dei metodi utilizzati per raggiungere l'obiettivo di reclutamento (vedi Figura 2). L'uso dei nuovi criteri diagnostici clinici per la fibromialgia ha permesso al team di studio di escludere correttamente gli individui che non soddisfacevano i criteri di fibromialgia32. Chiedere ai possibili partecipanti una diagnosi di fibromialgia non è una misura robusta come lo screening aggiuntivo e avrebbe potuto portare a visite di base programmate in modo improprio o alla partecipazione allo studio. Quarantotto partecipanti sono stati randomizzati in gruppi attivi e fittizi; Due partecipanti sono stati esclusi perché non soddisfacevano i criteri di ammissibilità basati sul test di inclusione per lo studio.

Figura 2: Diagramma di flusso del reclutamento. Un rapporto e un diagramma di flusso del reclutamento dello studio, della randomizzazione e dell'assegnazione dell'intervento. Clicca qui per visualizzare una versione più grande di questa figura.
Nei calcoli della dimensione del campione sono stati utilizzati tre risultati: DVPRS (dolore clinico), test di 30 s sulla sedia (funzione) e connettività DMN-SMN (RS-FCMRI). Tutti i calcoli della potenza si sono basati su dati preliminari. Le variazioni cliniche del dolore utilizzando DVPRS sono state scelte come esito primario di interesse. Il test in piedi della sedia 30 s (30sCST) è stato scelto come risultato funzionale rappresentativo poiché ha mostrato il più piccolo cambiamento tra i gruppi. La connettività DMN-SMN è stata scelta come risultato secondario di interesse come biomarcatore di neuroimaging per il dolore clinico e la risposta al trattamento. Le analisi delle dimensioni del campione sono state condotte assumendo una significatività dell'1% e dell'80% di potenza (2 campioni, 1 lato).
I semi per questa analisi sono stati scelti sulla base di dati preliminari, nonché della letteratura su fibromialgia, dolore e CES 16,24,26,63,64,65,66,67,68,69. Sulla base dei dati preliminari, la media (DS) della variazione post-trattamento della connettività tra corteccia sensomotoria primaria sinistra (L-S1M1) e corteccia cingolata posteriore sinistra (L-PCC) è 0,041 (0,079) per il gruppo di trattamento e -0,026 (0,049) per il gruppo di trattamento standard; la dimensione dell'effetto della differenza tra i gruppi osservata è 1,03 (vedi Tabella 2 16,21,26,62,63,64,70,71,72,73). Lo studio avrebbe bisogno di 20 soggetti in ciascun gruppo per raggiungere l'80% di potenza per rilevare la differenza tra il gruppo CES e il gruppo di trattamento standard nel loro cambiamento post-trattamento nella connettività per L-S1M1 a L-PCC al livello di significatività di 0,01 utilizzando un t-test a due lati, supponendo che la dimensione dell'effetto della differenza tra i gruppi sia 1,03 come osservato nei dati pilota. Sebbene lo studio pilota abbia osservato un tasso di abbandono del 17% per 12 soggetti nel nostro precedente studio sulla neuromodulazione auricolare (tutti hanno completato la risonanza magnetica di follow-up, ma due sono stati persi al follow-up a 8 settimane e 12 settimane), al fine di mantenere una stima conservativa per i calcoli della dimensione del campione, questo studio ha ipotizzato un tasso di abbandono del 20%. Con un attrito previsto del 20% alla visita post-trattamento, lo studio ha dovuto reclutare 20/0,8 = 25 soggetti per gruppo.
| Semi DMN (x,y,z) | Semi SMN (x,y,z) | Semi SN (x,y,z) |
| Corteccia Prefrontale Mediale62,67 | Putamen destro 64,71 | Corteccia prefrontale dorsolaterale destra62 |
| Destra PCC70 | Sinistra M170 | Insula anteriore sinistra62 |
| Sinistra PCC16,67 | Destra M170 | Insula anteriore destra62 |
| Precuneo71 | Destra S1-Mano16,72 | Insula posteriore sinistra64 |
| S1 Sinistro16,72 | Insula posteriore destra63 |
| Talamo21 | Corteccia cingolata anteriore dorsale72,73 |
| | Giunzione temporoparietale destra62 |
Tabella 2: Semi per l'analisi. Semi DMN, SMN e SN scelti per l'analisi basata su ipotesi a priori. Ogni seme è presentato con riferimenti alla letteratura precedente che supporta la sua sperimentazione nelle sindromi dolorose.
Sulla base di precedenti ricerche sulla CES in soggetti civili con fibromialgia, un totale di 50 soggetti (n = 25 sham e n = 25 vero) dovrebbero raggiungere l'80% di potenza per rilevare una differenza nei punteggi del dolore tra i due gruppi17 (vedi Tabella 3). I calcoli della dimensione del campione per questo studio sono stati condotti utilizzando sealedenvelope.com e si sono basati su dati preliminari.
| Controllo | Intervento | N per gruppo |
| Cambiamento medio | SD | Cambiamento medio | SD |
| Dolore clinico (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Funzione (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcMRI (da S1M1 a PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tabella 3: Calcolo della dimensione del campione. Calcoli relativi alla dimensione del campione di studio.
I dati funzionali di fMRIPrep importati sono stati levigati utilizzando la convoluzione spaziale con un kernel gaussiano di 8 mm di larghezza totale metà massima (FWHM) (vedere la Figura 3) mostra l'output funzionale di fMRIPrep normalizzato nello spazio del modello MNI152NLin2009cAsym (a sinistra) e l'immagine funzionale levigata da CONN Toolbox (a destra). Ciò si traduce in un aumento del rapporto segnale/rumore, che a sua volta migliora il rilevamento dei segnali dipendenti dal livello di ossigeno nel sangue (BOLD).

Figura 3: Singolo soggetto che confronta un'immagine funzionale non levigata nello spazio MNI (a sinistra) con la sua controparte levigata a 8 mm FWHM. Clicca qui per visualizzare una versione più grande di questa figura.
I dati sono stati quindi denoised utilizzando una pipeline di denoising standard74, inclusa la regressione dei potenziali effetti confondenti caratterizzati da serie temporali della sostanza bianca (5 componenti di rumore CompCor), serie temporali CSF (5 componenti di rumore CompCor), parametri di movimento e le loro derivate del primo ordine (12 fattori)75, scansioni anomale (inferiori a 295 fattori)48 e tendenze lineari (2 fattori) all'interno di ciascuna esecuzione funzionale, seguito dal filtraggio della frequenza passa-banda della serie temporale BOLD76 tra 0,008 Hz e 0,09 Hz. Le componenti di rumore CompCor49,77 all'interno della sostanza bianca e del CSF sono state stimate calcolando il segnale BOLD medio e le componenti principali più grandi ortogonali alla media BOLD, i parametri di movimento e le scansioni anomale all'interno delle maschere di segmentazione erose di ciascun soggetto (vedi Figura 4). Dal numero di termini di rumore inclusi in questa strategia di denoising, i gradi di libertà effettivi del segnale BOLD dopo la denoising sono stati stimati variare da 33 a 240,6 (media 173,4) in tutti i soggetti. La riduzione del rumore ha comportato una riduzione del rumore fisiologico e di altri rumori estranei dai dati che avrebbero potuto provocare effetti confondenti.

Figura 4: Controlli di qualità. Grafici di controllo qualità da CONN Toolbox che mostrano gli effetti del denoising sulla connettività funzionale (FC), sul segnale globale medio e sul movimento massimo. Nota (A) il grafico più in alto mostra i dati per un singolo argomento e sessione, (B,C) i grafici B e C sono i risultati della riduzione del rumore a livello di gruppo. Clicca qui per visualizzare una versione più grande di questa figura.
Le analisi a livello di gruppo sono state eseguite utilizzando un modello lineare generale (GLM)78. Per ogni singolo voxel, è stato stimato un GLM separato, con misure di connettività di primo livello a questo voxel come variabili dipendenti (un campione indipendente per soggetto e una misurazione per compito o condizione sperimentale, se applicabile) e gruppi o altri identificatori a livello di soggetto come variabili indipendenti. Le ipotesi a livello di voxel sono state valutate utilizzando statistiche parametriche multivariate con effetti casuali tra i soggetti e stima della covarianza del campione su più misurazioni. Le inferenze sono state eseguite a livello di singoli cluster (gruppi di voxel contigui). Le inferenze a livello di cluster si basavano sulle statistiche parametriche della teoria gaussiana dei campi casuali62,79. I risultati sono stati valutati utilizzando una combinazione di una soglia p < 0,005 a livello di voxel e una soglia p corretta per errore familiare < 0,0016 soglia80. Seguendo questi passaggi, si ottengono valori di connettività funzionale a livello di gruppo che confrontano il CES e le condizioni fittizie in base alle regioni di interesse (ROI). Questi risultati possono essere visualizzati in una moltitudine di modi attraverso l'interfaccia grafica di esplorazione dei risultati all'interno di CONN Toolbox. Per visualizzare una visualizzazione del volume di un esempio di risultati a livello di gruppo, vedere la Figura 5 con l'area rossa che indica le regioni con una maggiore connettività positiva con le ROI e le regioni blu con una maggiore connettività negativa con la ROI.

Figura 5: Visualizzazione del volume del seme post cingolato. L'immagine mostra una maggiore connettività positiva con il giro cingolato anteriore e (rosso) e una maggiore connettività negativa con il precuneo (blu) per la condizione True relativa allo Sham, pFWEc < 0,05. Questa cifra rappresenta dati parziali a livello di gruppo (n = 34). Clicca qui per visualizzare una versione più grande di questa figura.
File supplementare 1: Istruzioni per l'ordine del dispositivo CES Fare clic qui per scaricare questo file.
File supplementare 2: Registro del dispositivo CES Fare clic qui per scaricare questo file.
File supplementare 3: fMRIPrep Boilerplate Fare clic qui per scaricare questo file.
File supplementare 4: Istruzioni CES CONN Fare clic qui per scaricare questo file.
File supplementare 5: Grafici del codice R CES Fare clic qui per scaricare questo file.
File supplementare 6: Codice R CES eddy-qc Anova Fare clic qui per scaricare questo file.
Figura supplementare 1: Cartella di lavoro CONN. (A) Covariate di secondo livello. (B) Condizioni sperimentali. (C) Denoising. (d) Risultati. (E) Dati strutturali. Clicca qui per scaricare questa figura.