W odniesieniu do wyników rekrutacji, uczestnicy byli pozyskiwani głównie poprzez wysyłkę listów rekrutacyjnych oraz następujące po nich rozmowy telefoniczne, zgodnie z określonymi regulacjami Atlanta VA Healthcare System. Zespół badawczy zrekrutował łącznie 50 uczestników, co potwierdziło skuteczność zastosowanych metod w osiągnięciu celu rekrutacyjnego (patrz Rysunek 2). Zastosowanie nowych klinicznych kryteriów diagnostycznych fibromialgii pozwoliło zespołowi badawczemu na prawidłowe odfiltrowanie osób, które nie spełniały kryteriów fibromialgii32. Pytanie potencjalnych uczestników o rozpoznanie fibromialgii nie jest tak wiarygodną miarą jak dodatkowy screening i mogłoby doprowadzić do niewłaściwego zaplanowania wizyt wyjściowych lub błędnego włączenia do badania. Czterdziestu ośmiu uczestników zostało przydzielonych losowo do grupy aktywnej i grupy pozorowanej; dwóch uczestników zostało wykluczonych, ponieważ nie spełniali kryteriów kwalifikacji na podstawie testów włączenia do badania.

Rysunek 2: Schemat rekrutacji. Raport i schemat blokowy rekrutacji do badania, randomizacji oraz przydziału do interwencji. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Do obliczeń wielkości próby wykorzystano trzy punkty końcowe: DVPRS (ból kliniczny), test wstawania z krzesła przez 30 s (funkcja) oraz łączność DMN-SMN (rs-fcMRI). Wszystkie obliczenia mocy oparto na danych wstępnych. Jako główny punkt końcowy wybrano zmiany bólu klinicznego mierzone skalą DVPRS. Test wstawania z krzesła przez 30 s (30sCST) wybrano jako reprezentatywny punkt końcowy funkcjonalny, ponieważ wykazał on najmniejszą różnicę między grupami. Łączność DMN-SMN wybrano jako wtórny punkt końcowy, służący jako biomarker neuroobrazowy dla bólu klinicznego i odpowiedzi na leczenie. Analizy wielkości próby przeprowadzono przy założeniu poziomu istotności 1% i mocy 80% (test 2-próbkowy, jednostronny).
Nasiona dla tej analizy wybrano w oparciu o dane wstępne oraz literaturę dotyczącą fibromialgii, bólu i CES16,24,26,63,64,65,66,67,68,69. Na podstawie danych wstępnych średnia (SD) zmiany po leczeniu w łączności między lewą pierwotną korą czuciowo-ruchową (L-S1M1) a lewą tylną korą obręczy (L-PCC) wyniosła 0,041 (0,079) dla grupy leczonej i -0,026 (0,049) dla grupy standardowego leczenia; obserwowana wielkość efektu różnicy między grupami wyniosła 1,03 (patrz Tabela 216,21,26,62,63,64,70,71,72,73). W badaniu potrzeba byłoby 20 osób w każdej grupie, aby osiągnąć moc 80% w celu wykrycia różnicy między grupą CES a grupą standardowego leczenia w zmianie łączności między L-S1M1 a L-PCC po leczeniu, przy poziomie istotności 0,01 przy użyciu dwustronnego testu t, zakładając, że wielkość efektu różnicy między grupami wynosi 1,03, zgodnie z obserwacjami z danych pilotażowych. Chociaż w badaniu pilotażowym odnotowano 17% wskaźnik rezygnacji wśród 12 osób w naszym poprzednim badaniu nad neuromodulacją małżowiny usznej (wszyscy ukończyli kontrolne badanie MRI, ale dwie osoby odpadły z obserwacji w 8. i 12. tygodniu), w celu zachowania konserwatywnego szacunku dla obliczeń wielkości próby, w niniejszym badaniu przyjęto 20% wskaźnik rezygnacji. Przy oczekiwanym 20% wskaźniku rezygnacji podczas wizyty po leczeniu, badanie wymagało zrekrutowania 20/0,8 = 25 osób na grupę.
| Ziarna DMN (x,y,z) | Ziarna SMN (x,y,z) | Ziarna SN (x,y,z) |
| Przyśrodkowa kora przedczołowa62,67 | Prawe jądro skorupowe64,71 | Prawa grzbietowo-boczna kora przedczołowa62 |
| Prawy PCC70 | Lewy M170 | Lewa przednia część wyspy62 |
| Lewy PCC16,67 | Prawy M170 | Prawa przednia część wyspy62 |
| Kuneus71 | Prawy S1-dłoń16,72 | Lewa tylna część wyspy64 |
| Lewy S1-dłoń16,72 | Prawa tylna część wyspy63 |
| Wzgórze21 | Grzbietowa przednia kora obręczy72,73 |
| | Prawy złącze skroniowo-ciemieniowe62 |
Tabela 2: Regiony seed do analizy. Regiony seed DMN, SMN i SN wybrane do analizy na podstawie hipotez a priori. Przy każdym regionie seed podano odniesienia do wcześniejszej literatury uzasadniającej jego badanie w zespołach bólowych.
W oparciu o wcześniejsze badania nad CES u osób cywilnych z fibromialgią, łączna liczba 50 uczestników (n = 25 w grupie pozorowanej i n = 25 w grupie rzeczywistej) powinna zapewnić 80% mocy statystycznej niezbędnej do wykrycia różnicy w ocenach bólu między dwiema grupami17 (patrz Tabela 3). Obliczenia wielkości próby dla niniejszego badania zostały przeprowadzone przy użyciu serwisu sealedenvelope.com i oparto je na danych wstępnych.
| Kontrola | Interwencja | liczba N w grupie |
| Średnia zmiana | OD | Średnia zmiana | Odchylenie standardowe |
| Ból kliniczny (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Funkcja (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcMRI (S1M1 do PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tabela 3: Obliczanie wielkości próby. Obliczenia związane z wielkością próby badania.
Zaimportowane funkcjonalne dane fMRIPrep zostały wygładzone przy użyciu splotu przestrzennego z jądrem Gaussa o pełnej szerokości połówki maksimum (FWHM) wynoszącej 8 mm (patrz Rysunek 3), który przedstawia funkcjonalny wynik z fMRIPrep znormalizowany do przestrzeni szablonu MNI152NLin2009cAsym (po lewej) oraz wygładzony obraz funkcjonalny z CONN Toolbox (po prawej). Przynosi to wzrost stosunku sygnału do szumu, co z kolei poprawia detekcję sygnałów zależnych od poziomu utlenowania krwi (BOLD).

Rycina 3: Pojedynczy badany; porównanie niewygładzonego obrazu funkcjonalnego w przestrzeni MNI (lewo) z jego wygładzonym odpowiednikiem przy 8 mm FWHM. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Następnie dane odszumiono za pomocą standardowego potoku odszumiania74, obejmującego regresję potencjalnych efektów zakłócających określonych przez szeregi czasowe istoty białej (5 komponentów szumu CompCor), szeregi czasowe CSF (5 komponentów szumu CompCor), parametry ruchu i ich pochodne pierwszego rzędu (12 czynników)75, skany odstające (poniżej 295 czynników)48 oraz trendy liniowe (2 czynniki) w każdym przebiegu funkcjonalnym, a następnie filtrowanie pasmowe szeregów czasowych BOLD76 w zakresie od 0.008 Hz do 0.09 Hz. Komponenty szumu CompCor49,77 w obrębie istoty białej i CSF oszacowano poprzez obliczenie średniego sygnału BOLD, a także największych głównych komponentów ortogonalnych do średniej BOLD, parametrów ruchu i skanów odstających w obrębie erodowanych masek segmentacji każdego badanego (patrz Rysunek 4). Biorąc pod uwagę liczbę składników szumu uwzględnionych w tej strategii odszumiania, oszacowano, że efektywne stopnie swobody sygnału BOLD po odszumianiu wynosiły od 33 do 240.6 (średnio 173.4) we wszystkich badanych. Odszumianie doprowadziło do redukcji szumu fizjologicznego i innych szumów zewnętrznych w danych, które mogłyby spowodować efekty zakłócające.

Rysunek 4: Kontrola jakości. Wykresy kontroli jakości z CONN Toolbox przedstawiające wpływ odszumiania na łączność funkcjonalną (FC), średni sygnał globalny oraz maksymalny ruch. Uwaga (A) najwyższy wykres przedstawia dane dla pojedynczego badanego i sesji, (B,C) wykresy B i C są wynikami odszumiania na poziomie grupowym. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Analizy na poziomie grupowym przeprowadzono przy użyciu ogólnego modelu liniowego (GLM)78. Dla każdego pojedynczego woksela oszacowano oddzielny model GLM, przyjmując miary łączności na pierwszym poziomie dla tego woksela jako zmienne zależne (jedna niezależna próba na badanego i jeden pomiar na zadanie lub warunek eksperymentalny, jeśli dotyczy) oraz grupy lub inne identyfikatory na poziomie badanych jako zmienne niezależne. Hipotezy na poziomie wokseli oceniano za pomocą wielowymiarowej statystyki parametrycznej z efektami losowymi między badanymi oraz szacowaniem kowariancji próby dla wielu pomiarów. Wnioskowanie przeprowadzono na poziomie poszczególnych klastrów (grup sąsiednich wokseli). Wnioskowanie na poziomie klastrów oparto na statystyce parametrycznej z teorii gaussowskich pól losowych62,79. Wyniki poddano progowaniu, stosując kombinację progu na poziomie wokseli p < 0.005 dla tworzenia klastrów oraz progu wielkości klastra p < 0.0016 z korekcją błędu rodzinnego (familywise error)80. Postępowanie zgodnie z tymi krokami prowadzi do uzyskania wartości łączności funkcjonalnej na poziomie grupowym, porównujących warunki CES i sham w oparciu o obszary zainteresowania (ROI). Wyniki te można wizualizować na wiele sposobów za pomocą interfejsu GUI eksploratora wyników w CONN Toolbox. Aby zapoznać się z wizualizacją wolumetryczną przykładowych wyników na poziomie grupowym, gdzie obszar czerwony wskazuje regiony o większej dodatniej łączności z ROI, a obszar niebieski regiony o większej ujemnej łączności z ROI, należy przejść do Rysunku 5.

Rycina 5: Wyświetlanie objętościowe ziarna z zakrętu obręczy tylnego. Obraz przedstawia silniejszą dodatnią łączność z zakrętem obręczy przedniego (kolor czerwony) oraz silniejszą ujemną łączność z przedkliniem (kolor niebieski) w warunku True w porównaniu do warunku Sham, pFWEc < 0.05. Rycina ta przedstawia częściowe dane na poziomie grupy (n = 34). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Plik uzupełniający 1: Instrukcje zamawiania urządzenia CES Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Dziennik urządzenia CES Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 3: Szablon fMRIPrep Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 4: Instrukcje CES CONN Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 5: Kod R do wykresów CES Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 6: Kod R CES eddy-qc Anova Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 1: Arkusz roboczy CONN.(A) Współzmienne drugiego poziomu. (B) Warunki eksperymentalne. (C) Odszumianie. (D) Wyniki. (E) Dane strukturalne. Kliknij tutaj, aby pobrać tę rycinię.