$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In termen van wervingsresultaten werden deelnemers voornamelijk geworven via het mailen van wervingsbrieven en vervolgtelefoontjes op basis van de geschetste voorschriften van het Atlanta VA Healthcare System. Het onderzoeksteam rekruteerde in totaal 50 deelnemers, wat de effectiviteit bewijst van de methoden die zijn gebruikt om het wervingsdoel te bereiken (zie figuur 2). Het gebruik van de nieuwe klinische diagnostische criteria voor fibromyalgie stelde het onderzoeksteam in staat om personen die niet aan de fibromyalgiecriteria voldeden goed te screenen32. Het vragen van mogelijke deelnemers naar een diagnose van fibromyalgie is niet zo'n robuuste maatregel als de aanvullende screening en had kunnen leiden tot onjuist geplande basislijnbezoeken of deelname aan het onderzoek. Achtenveertig deelnemers werden gerandomiseerd in actieve en schijngroepen; Twee deelnemers werden uitgesloten omdat ze niet voldeden aan de geschiktheidscriteria op basis van de inclusietest voor het onderzoek.

Figuur 2: Stroomschema voor werving. Een rapport en stroomdiagram van studierekrutering, randomisatie en toewijzing van interventie. Klik hier om een grotere versie van deze figuur te bekijken.
Drie uitkomsten werden gebruikt bij de berekening van de steekproefomvang: DVPRS (klinische pijn), 30 s stoelstandtest (functie) en DMN-SMN-connectiviteit (rs-fcMRI). Alle vermogensberekeningen waren gebaseerd op voorlopige gegevens. Klinische pijnveranderingen met behulp van DVPRS werden gekozen als de primaire uitkomst die van belang is. De 30 s stoelstandtest (30sCST) werd gekozen als de representatieve functionele uitkomst omdat deze de kleinste verandering tussen de groepen vertoonde. DMN-SMN-connectiviteit werd gekozen als de secundaire uitkomst die van belang is als biomarker voor neuroimaging voor klinische pijn en respons op de behandeling. Er werden analyses van de steekproefomvang uitgevoerd, uitgaande van een significantie van 1% en 80% vermogen (2-steekproef, 1-zijdig).
De zaden voor deze analyse werden gekozen op basis van voorlopige gegevens, evenals de literatuur over fibromyalgie, pijn en CES 16,24,26,63,64,65,66,67,68,69. Op basis van de voorlopige gegevens is het gemiddelde (SD) van de verandering na de behandeling in de connectiviteit van de linker primaire sensomotorische cortex (L-S1M1) naar de linker posterieure cingulate cortex (L-PCC) 0,041 (0,079) voor de behandelingsgroep en -0,026 (0,049) voor de standaardbehandelingsgroep; de waargenomen verschilgrootte van het effect tussen de groepen is 1,03 (zie Tabel 2 16,21,26,62,63,64,70,71,72,73). De studie zou 20 proefpersonen in elke groep nodig hebben om 80% vermogen te bereiken om het verschil tussen de CES-groep en de standaardbehandelingsgroep te detecteren in hun verandering na de behandeling in connectiviteit voor L-S1M1 naar L-PCC op het significantieniveau van 0,01 met behulp van een tweezijdige t-toets, ervan uitgaande dat de verschileffectgrootte tussen de groep 1,03 is, zoals waargenomen in de pilotgegevens. Hoewel de pilotstudie een uitvalpercentage van 17% waarnam voor 12 proefpersonen in onze eerdere studie van auriculaire neuromodulatie (ze voltooiden allemaal hun follow-up MRI, maar twee gingen verloren voor follow-up na 8 weken en 12 weken), om een conservatieve schatting te behouden voor berekeningen van de steekproefomvang, ging deze studie uit van een uitvalpercentage van 20%. Met een verwacht verloop van 20% bij het bezoek na de behandeling, moest de studie 20/0,8 = 25 proefpersonen per groep rekruteren.
| DMN zaden (x,y,z) | SMN zaden (x,y,z) | SN zaden (x,y,z) |
| Mediale prefrontale cortex 62,67 | Rechts Putamen 64,71 | Rechter dorsolaterale prefrontale cortex62 |
| Rechts PCC70 | Links M170 | Linker voorste insula62 |
| Links PCC16,67 | Rechts M170 | Rechter voorste insula62 |
| Precuneus71 | Rechts S1-Hand16,72 | Linker posterieure insula64 |
| Links S1-wijzer16,72 | Rechter posterieure insula63 |
| Thalamus21 | Dorsale anterieure cingulate cortex72,73 |
| | Rechter temporopariëtale overgang62 |
Tabel 2: Zaden voor analyse. DMN-, SMN- en SN-zaden gekozen voor analyse op basis van a priori hypothesen. Elk zaadje wordt gepresenteerd met verwijzingen naar eerdere literatuur die het testen bij pijnsyndromen ondersteunt.
Op basis van eerder onderzoek naar CES bij civiele proefpersonen met fibromyalgie, zouden in totaal 50 proefpersonen (n = 25 schijnvertoning en n = 25 waar) 80% vermogen moeten bereiken om een verschil in pijnscores tussen de twee groepen te detecteren17 (zie tabel 3). De berekeningen van de steekproefomvang voor deze studie werden uitgevoerd met behulp van sealedenvelope.com en waren gebaseerd op voorlopige gegevens.
| Beheersen | Interventie | N per groep |
| Gemiddelde verandering | SD | Gemiddelde verandering | SD |
| Klinische pijn (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Functie (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcMRI (S1M1 naar PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tabel 3: Berekening van de steekproefomvang. Berekeningen met betrekking tot de omvang van de studiesteekproef.
Geïmporteerde fMRIPrep-functionele gegevens werden gladgestreken met behulp van ruimtelijke convolutie met een Gaussiaanse kern van 8 mm over de volledige breedte, half maximum (FWHM) (zie figuur 3), toont de functionele uitvoer van fMRIPrep genormaliseerd in MNI152NLin2009cAsym-sjabloonruimte (links) en het afgevlakte functionele beeld van CONN Toolbox (rechts). Dit resulteert in een verhoging van de signaal-ruisverhouding, wat op zijn beurt de detectie van bloed-zuurstofniveau-afhankelijke (BOLD) signalen verbetert.

Figuur 3: Enkel onderwerp waarin een niet-gladgestreken functioneel beeld in de MNI-ruimte (links) wordt vergeleken met zijn afgevlakte tegenhanger op 8 mm FWHM. Klik hier om een grotere versie van deze figuur te bekijken.
De gegevens werden vervolgens ontnoise met behulp van een standaard ruisonderdrukkingspijplijn74, inclusief de regressie van mogelijke verstorende effecten die worden gekenmerkt door tijdreeksen van witte stof (5 CompCor-ruiscomponenten), CSF-tijdreeksen (5 CompCor-ruiscomponenten), bewegingsparameters en hun eerste-orde-derivaten (12 factoren)75, uitbijterscans (minder dan 295 factoren)48 en lineaire trends (2 factoren) binnen elke functionele run, gevolgd door bandpass-frequentiefiltering van de BOLD-tijdreeks76 tussen 0,008 Hz en 0,09 Hz. CompCor49,77 ruiscomponenten binnen witte stof en CSF werden geschat door het gemiddelde BOLD-signaal te berekenen, evenals de grootste hoofdcomponenten loodrecht op het BOLD-gemiddelde, bewegingsparameters en uitbijterscans binnen de geërodeerde segmentatiemaskers van elke proefpersoon (zie figuur 4). Op basis van het aantal ruistermen dat in deze ruisonderdrukkingsstrategie was opgenomen, werd de effectieve vrijheidsgraad van het BOLD-signaal na ruisonderdrukking geschat op 33 tot 240,6 (gemiddeld 173,4) voor alle proefpersonen. De ruisonderdrukking resulteerde in een vermindering van fysiologische en andere externe ruis uit de gegevens die tot verstorende effecten had kunnen leiden.

Figuur 4: Kwaliteitscontroles. Grafieken voor kwaliteitscontroles van CONN Toolbox die de effecten van ruisonderdrukking op functionele connectiviteit (FC), gemiddeld wereldwijd signaal en maximale beweging weergeven. Opmerking (A) de bovenste grafiek geeft gegevens weer voor een enkel onderwerp en sessie, (B,C) grafieken B en C zijn de resultaten van ruisonderdrukking op groepsniveau. Klik hier om een grotere versie van deze figuur te bekijken.
Analyses op groepsniveau werden uitgevoerd met behulp van een Algemeen Lineair Model (GLM)78. Voor elke individuele voxel werd een afzonderlijke GLM geschat, met connectiviteitsmetingen op het eerste niveau bij deze voxel als afhankelijke variabelen (één onafhankelijke steekproef per proefpersoon en één meting per taak of experimentele omstandigheid, indien van toepassing), en groepen of andere identificatiemiddelen op onderwerpniveau als onafhankelijke variabelen. Hypothesen op Voxel-niveau werden geëvalueerd met behulp van multivariate parametrische statistieken met willekeurige effecten bij proefpersonen en schatting van de steekproefcovariantie bij meerdere metingen. Gevolgtrekkingen werden uitgevoerd op het niveau van individuele clusters (groepen van aaneengesloten voxels). Gevolgtrekkingen op clusterniveau waren gebaseerd op parametrische statistieken uit de Gaussiaanse toevalsveldentheorie 62,79. De resultaten werden gedrempeld met behulp van een combinatie van een clustervormende p < een voxelniveaudrempel van 0,005 en een gezinsgewijze fout gecorrigeerde p < 0,0016 clustergroottedrempel80. Het volgen van deze stappen resulteert in functionele connectiviteitswaarden op groepsniveau die de CES- en schijnomstandigheden vergelijken op basis van interessegebieden (ROI). Deze resultaten kunnen op verschillende manieren worden gevisualiseerd via de resultatenverkenner GUI binnen CONN Toolbox. Zie afbeelding 5 voor een visualisatie van volumeweergave van voorbeeldresultaten op groepsniveau, waarbij het rode gebied regio's aangeeft met een grotere positieve connectiviteit met de ROI's en het blauwe gebied regio's zijn met een grotere negatieve verbinding met de ROI.

Figuur 5: Volumeweergave van postcingulaat zaad. Het beeld vertoont een grotere positieve connectiviteit met de gyrus cingularis anterior en (rood) en een grotere negatieve connectiviteit met de precuneus (blauw) voor de ware toestand ten opzichte van de schijnvertoning, pFWEc < 0,05. Dit cijfer vertegenwoordigt gedeeltelijke gegevens op groepsniveau (n = 34). Klik hier om een grotere versie van deze figuur te bekijken.
Aanvullend bestand 1: Instructies voor het bestellen van CES-apparaten Klik hier om dit bestand te downloaden.
Aanvullend bestand 2: CES-apparaatlogboek Klik hier om dit bestand te downloaden.
Aanvullend bestand 3: fMRIPrep Boilerplate Klik hier om dit bestand te downloaden.
Aanvullend bestand 4: CES CONN instructies Klik hier om dit bestand te downloaden.
Aanvullend bestand 5: CES R-codeplots Klik hier om dit bestand te downloaden.
Aanvullend bestand 6: R-code CES eddy-qc Anova Klik hier om dit bestand te downloaden.
Aanvullende figuur 1: CONN-werkboek. (A) Covariabelen op het tweede niveau. (B) Experimentele omstandigheden. (C) Ruisonderdrukking. d) Resultaten. (E) Structurele gegevens. Klik hier om deze afbeelding te downloaden.