$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In Bezug auf die Rekrutierungsergebnisse wurden die Teilnehmer in erster Linie durch das Versenden von Rekrutierungsschreiben und Nachverfolgungsanrufe auf der Grundlage der dargelegten Vorschriften des Atlanta VA Healthcare System rekrutiert. Das Studienteam rekrutierte insgesamt 50 Teilnehmer, was die Wirksamkeit der verwendeten Methoden zur Erreichung des Rekrutierungsziels belegt (siehe Abbildung 2). Die Verwendung der neuen diagnostischen Kriterien für klinische Fibromyalgie ermöglichte es dem Studienteam, Personen, die die Fibromyalgie-Kriterien32 nicht erfüllten, ordnungsgemäß auszusortieren. Die Befragung möglicher Teilnehmer zu einer Diagnose von Fibromyalgie ist kein so robustes Maß wie das zusätzliche Screening und hätte zu falsch geplanten Baseline-Besuchen oder Studienteilnahmen führen können. Achtundvierzig Teilnehmer wurden nach dem Zufallsprinzip in aktive und Scheingruppen eingeteilt; Zwei Teilnehmer wurden ausgeschlossen, weil sie die Zulassungskriterien auf der Grundlage der Einschlusstests für die Studie nicht erfüllten.

Abbildung 2: Flussdiagramm für die Rekrutierung. Ein Bericht und ein Flussdiagramm der Studienrekrutierung, Randomisierung und Zuweisung von Interventionen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Drei Endpunkte wurden bei der Berechnung der Stichprobengröße verwendet: DVPRS (klinischer Schmerz), 30 s Stuhlstandtest (Funktion) und DMN-SMN-Konnektivität (rs-fcMRI). Alle Leistungsberechnungen basierten auf vorläufigen Daten. Klinische Schmerzveränderungen mittels DVPRS wurden als primäres Ergebnis von Interesse gewählt. Der 30 s Stuhlstandtest (30sCST) wurde als repräsentatives funktionelles Ergebnis gewählt, da er die geringste Veränderung zwischen den Gruppen aufwies. Die DMN-SMN-Konnektivität wurde als sekundärer Endpunkt von Interesse als Neuroimaging-Biomarker für klinische Schmerzen und das Ansprechen auf die Behandlung ausgewählt. Die Stichprobengrößenanalysen wurden unter der Annahme einer Signifikanz von 1 % und 80 % Trennschärfe durchgeführt (2 Stichproben, 1-seitig).
Die Samen für diese Analyse wurden auf der Grundlage vorläufiger Daten sowie der Literatur zu Fibromyalgie, Schmerzen und CES 16,24,26,63,64,65,66,67,68,69 ausgewählt. Basierend auf den vorläufigen Daten beträgt der Mittelwert (SD) der Veränderung der Konnektivität des linken primären sensomotorischen Kortex (L-S1M1) nach der Behandlung zur linken posterioren cingulären Kortex (L-PCC) 0,041 (0,079) für die Behandlungsgruppe und -0,026 (0,049) für die Standardbehandlungsgruppe; Die beobachtete Differenzeffektgröße zwischen den Gruppen beträgt 1,03 (siehe Tabelle 2: 16,21,26,62,63,64,70,71,72,73). Die Studie würde 20 Probanden in jeder Gruppe benötigen, um eine Schärfe von 80 % zu erreichen, um den Unterschied zwischen der CES-Gruppe und der Standardbehandlungsgruppe in ihrer Veränderung der Konnektivität von L-S1M1 zu L-PCC nach der Behandlung auf dem Signifikanzniveau von 0,01 unter Verwendung eines zweiseitigen t-Tests zu erkennen, unter der Annahme, dass die Differenzeffektgröße zwischen den Gruppen 1,03 beträgt, wie in den Pilotdaten beobachtet. Obwohl in der Pilotstudie eine Abnutzungsrate von 17 % bei 12 Probanden in unserer vorherigen Studie zur aurikulären Neuromodulation beobachtet wurde (alle schlossen ihre Nachbeobachtungs-MRT ab, aber zwei gingen nach 8 Wochen und 12 Wochen für die Nachbeobachtung verloren), um eine konservative Schätzung für die Berechnung der Stichprobengröße beizubehalten, ging diese Studie von einer Abnutzungsrate von 20 % aus. Bei einer erwarteten Fluktuation von 20 % beim Nachbehandlungsbesuch musste die Studie 20/0,8 = 25 Probanden pro Gruppe rekrutieren.
| DMN-Saatgut (x,y,z) | SMN-Seeds (x,y,z) | SN-Seeds (x,y,z) |
| Medialer präfrontaler Kortex62,67 | Rechts Putamen 64,71 | Rechter dorsolateraler präfrontaler Kortex62 |
| Rechts PCC70 | Links M170 | Linke vordere Insula62 |
| Links PCC16,67 | Rechts M170 | Rechte vordere Insula62 |
| Precuneus71 | Rechts S1-Hand16,72 | Hintere Insulalinks 64 |
| Links S1-Hand16,72 | Insula rechts posterior63 |
| Thalamus21 | Dorsaler anteriorer cingulärer Kortex72,73 |
| | Interner temporoparietaler Übergangrechts 62 |
Tabelle 2: Saatgut für die Analyse. DMN-, SMN- und SN-Seeds, die für die Analyse auf der Grundlage von A-priori-Hypothesen ausgewählt wurden. Jedes Saatgut wird mit Verweisen auf frühere Literatur präsentiert, die seine Tests bei Schmerzsyndromen unterstützt.
Basierend auf früheren Forschungen zu CES bei zivilen Probanden mit Fibromyalgie sollten insgesamt 50 Probanden (n = 25 Schein und n = 25 wahr) eine Aussagekraft von 80 % erreichen, um einen Unterschied in den Schmerzwerten zwischen den beiden Gruppen zu erkennen17 (siehe Tabelle 3). Die Berechnungen der Stichprobengröße für diese Studie wurden unter Verwendung von sealedenvelope.com durchgeführt und basierten auf vorläufigen Daten.
| Steuerung | Eingriff | N pro Gruppe |
| Mittlere Veränderung | SD | Mittlere Veränderung | SD |
| Klinischer Schmerz (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Funktion (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcMRI (S1M1 auf PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tabelle 3: Berechnung des Stichprobenumfangs. Berechnungen im Zusammenhang mit der Stichprobengröße der Studie.
Importierte fMRIPrep-Funktionsdaten wurden mittels räumlicher Faltung mit einem Gaußschen Kern von 8 mm voller Breite halbes Maximum (FWHM) geglättet (siehe Abbildung 3) zeigt die funktionale Ausgabe von fMRIPrep, normiert in den MNI152NLin2009cAsym-Vorlagenraum (links) und das geglättete Funktionsbild aus der CONN Toolbox (rechts). Dies führt zu einer Erhöhung des Signal-Rausch-Verhältnisses, was wiederum die Detektion von blutsauerstoffabhängigen (BOLD) Signalen verbessert.

Abbildung 3: Einzelnes Motiv beim Vergleich eines ungeglätteten Funktionsbildes im MNI-Raum (links) mit seinem geglätteten Gegenstück bei 8 mm FWHM. Bitte klicken Sie hier, um eine größere Version dieser Abbildung zu sehen.
Die Daten wurden dann unter Verwendung einer Standard-Rauschunterdrückungspipeline74 entrauscht, einschließlich der Regression potenzieller Störeffekte, die durch Zeitreihen der weißen Substanz (5 CompCor-Rauschkomponenten), CSF-Zeitreihen (5 CompCor-Rauschkomponenten), Bewegungsparameter und deren Ableitungen erster Ordnung (12 Faktoren)75, Ausreißerscans (unter 295 Faktoren)48 und lineare Trends (2 Faktoren) innerhalb jedes funktionalen Laufs gekennzeichnet sind. gefolgt von einer Bandpass-Frequenzfilterung der BOLD-Zeitreihe76 zwischen 0,008 Hz und 0,09 Hz. CompCor49,77 Rauschkomponenten innerhalb der weißen Substanz und des Liquors wurden durch Berechnung des durchschnittlichen BOLD-Signals sowie der größten Hauptkomponenten orthogonal zum BOLD-Durchschnitt, der Bewegungsparameter und der Ausreißer-Scans innerhalb der erodierten Segmentierungsmasken jedes Probanden geschätzt (siehe Abbildung 4). Aus der Anzahl der Rauschterme, die in diese Rauschunterdrückungsstrategie einbezogen wurden, wurden die effektiven Freiheitsgrade des BOLD-Signals nach der Rauschunterdrückung auf 33 bis 240,6 (Durchschnitt 173,4) für alle Probanden geschätzt. Die Rauschunterdrückung führte zu einer Verringerung des physiologischen und anderen Fremdrauschens aus den Daten, das zu Störeffekten hätte führen können.

Abbildung 4: Qualitätsprüfungen. Qualitätsprüfungsdiagramme aus der CONN Toolbox, die die Auswirkungen der Rauschunterdrückung auf die funktionale Konnektivität (FC), das mittlere globale Signal und die maximale Bewegung anzeigen. Hinweis: (A) Das oberste Diagramm zeigt Daten für ein einzelnes Fach und eine einzelne Sitzung an, (B,C) Die Diagramme B und C sind die Ergebnisse der Rauschunterdrückung auf Gruppenebene. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Analysen auf Gruppenebene wurden unter Verwendung eines General Linear Model (GLM)78 durchgeführt. Für jedes einzelne Voxel wurde ein separater GLM geschätzt, wobei die Konnektivitätsmaße der ersten Ebene an diesem Voxel als abhängige Variablen (eine unabhängige Stichprobe pro Subjekt und eine Messung pro Aufgabe oder Versuchsbedingung, falls zutreffend) und Gruppen oder andere Identifikatoren auf Probandenebene als unabhängige Variablen verwendet wurden. Hypothesen auf Voxel-Ebene wurden unter Verwendung multivariater parametrischer Statistiken mit zufälligen Effekten über Probanden hinweg und Schätzung der Kovarianzschätzung der Stichprobe über mehrere Messungen hinweg bewertet. Die Inferenzen wurden auf der Ebene einzelner Cluster (Gruppen von zusammenhängenden Voxeln) durchgeführt. Die Inferenzen auf Clusterebene basierten auf parametrischen Statistiken der Gaußschen Zufallsfeldtheorie62,79. Die Ergebnisse wurden mit einer Kombination aus einem clusterbildenden p < 0,005 Voxel-Schwellenwert und einem familienbezogenen fehlerkorrigierten p < 0,0016 Clustergrößen-Schwellenwert80 ermittelt. Das Befolgen dieser Schritte führt zu funktionalen Konnektivitätswerten auf Gruppenebene, bei denen die CES- und Scheinbedingungen basierend auf den relevanten Regionen (ROI) verglichen werden. Diese Ergebnisse können auf vielfältige Weise über die Ergebnis-Explorer-GUI in der CONN Toolbox visualisiert werden. Eine Visualisierung der Volumenanzeige eines Beispiels für Ergebnisse auf Gruppenebene, wobei der rote Bereich Regionen mit größerer positiver Konnektivität mit den ROIs und der blaue Bereich Regionen mit größerer negativer Konnektivität mit dem ROI anzeigt, siehe Abbildung 5.

Abbildung 5: Volumendarstellung des postcingulären Saatguts. Das Bild zeigt eine größere positive Konnektivität mit dem Gyrus cinguli anterior (rot) und eine größere negative Konnektivität mit dem Precuneus (blau) für die True-Bedingung relativ zum Sham, pFWEc < 0,05. Bei dieser Zahl handelt es sich um partielle Daten auf Gruppenebene (n = 34). Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Ergänzende Datei 1: Bestellanweisungen für CES-Geräte Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 2: CES-Geräteprotokoll Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 3: fMRIPrep Boilerplate Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 4: CES CONN-Anweisungen Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 5: CES R-Code-Plots Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 6: R-Code CES eddy-qc Anova Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung 1: CONN-Arbeitsmappe. (A) Kovariaten der zweiten Ebene. (B) Versuchsbedingungen. (c) Entrauschung. (D) Ergebnisse. (E) Strukturdaten. Bitte klicken Sie hier, um diese Abbildung herunterzuladen.