$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
En ce qui concerne les résultats du recrutement, les participants ont été principalement recrutés par l’envoi de lettres de recrutement et d’appels téléphoniques de suivi basés sur les réglementations définies du système de santé de l’Atlanta VA. L’équipe de recherche a recruté un total de 50 participants, ce qui prouve l’efficacité des méthodes utilisées pour atteindre l’objectif de recrutement (voir la figure 2). L’utilisation des nouveaux critères de diagnostic clinique de la fibromyalgie a permis à l’équipe de l’étude d’éliminer correctement les personnes qui ne répondaient pas aux critères32 de la fibromyalgie. Interroger les participants potentiels sur un diagnostic de fibromyalgie n’est pas une mesure aussi robuste que le dépistage supplémentaire et aurait pu conduire à des visites de base ou à une participation à l’étude mal planifiées. Quarante-huit participants ont été répartis au hasard en groupes actifs et simulés ; Deux participants ont été exclus parce qu’ils ne répondaient pas aux critères d’admissibilité basés sur le test d’inclusion de l’étude.

Figure 2 : Organigramme de recrutement. Un rapport et un organigramme du recrutement de l’étude, de la randomisation et de l’allocation de l’intervention. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Trois critères de jugement ont été utilisés dans le calcul de la taille de l’échantillon : DVPRS (douleur clinique), test sur chaise debout de 30 s (fonction) et connectivité DMN-SMN (rs-fcMRI). Tous les calculs de puissance ont été basés sur des données préliminaires. Les changements cliniques de la douleur à l’aide de la DVPRS ont été choisis comme critère de jugement principal. Le test du support de chaise 30 s (30sCST) a été choisi comme résultat fonctionnel représentatif car il présentait le plus petit changement entre les groupes. La connectivité DMN-SMN a été choisie comme critère de jugement secondaire d’intérêt en tant que biomarqueur de neuroimagerie pour la douleur clinique et la réponse au traitement. Des analyses de la taille de l’échantillon ont été effectuées en supposant une signification de 1 % et 80 % de puissance (2 échantillons, 1 face).
Les graines de cette analyse ont été choisies sur la base de données préliminaires, ainsi que de la littérature sur la fibromyalgie, la douleur et CES 16,24,26,63,64,65,66,67,68,69. D’après les données préliminaires, la moyenne (ET) de la variation post-traitement de la connectivité du cortex sensorimoteur primaire gauche (L-S1M1) au cortex cingulaire postérieur gauche (L-PCC) est de 0,041 (0,079) pour le groupe de traitement et de -0,026 (0,049) pour le groupe de traitement standard ; l’ampleur de l’effet de la différence observée entre les groupes est de 1,03 (voir tableau 2: 16, 21, 26, 62, 63, 64, 70, 71, 72, 73). L’étude aurait besoin de 20 sujets dans chaque groupe pour atteindre une puissance de 80 % afin de détecter la différence entre le groupe CES et le groupe de traitement standard dans leur changement post-traitement de la connectivité pour L-S1M1 à L-PCC au niveau de signification de 0,01 à l’aide d’un test t bilatéral, en supposant que la taille de l’effet de la différence entre les groupes est de 1,03, comme observé dans les données pilotes. Bien que l’étude pilote ait observé un taux d’attrition de 17 % chez 12 sujets dans notre étude antérieure sur la neuromodulation auriculaire (tous ont terminé leur IRM de suivi, mais deux ont été perdus de vue à 8 semaines et 12 semaines), afin de maintenir une estimation prudente pour le calcul de la taille de l’échantillon, cette étude a supposé un taux d’attrition de 20 %. Avec un taux d’attrition attendu de 20 % lors de la visite post-traitement, l’étude devait recruter 20/0,8 = 25 sujets par groupe.
| Graines de DMN (x,y,z) | Graines SMN (x,y,z) | Graines SN (x,y,z) |
| Cortex préfrontal médial62,67 | Putamen droit 64,71 | Cortex préfrontal dorsolatéral droit62 |
| PCC70 droit | Gauche M170 | Insula antérieure gauche62 |
| Gauche PCC16,67 | Droite M170 | insula antérieure droite62 |
| Precuneus71 | Main droite S116,72 | Insula postérieure gauche64 |
| Main gauche S116,72 | insula postérieure droite63 |
| Thalamus21 | Cortex cingulaire antérieur dorsal72,73 |
| | Jonction temporo-pariétale droite62 |
Tableau 2 : Graines pour l’analyse. Les graines DMN, SMN et SN sont choisies pour l’analyse sur la base d’hypothèses a priori. Chaque graine est présentée avec des références à la littérature antérieure soutenant ses tests dans les syndromes douloureux.
Sur la base de recherches antérieures sur le CES chez des sujets civils atteints de fibromyalgie, un total de 50 sujets (n = 25 simulé et n = 25 vrai) devrait atteindre une puissance de 80% pour détecter une différence dans les scores de douleur entre les deux groupes17 (voir tableau 3). Les calculs de la taille de l’échantillon pour cette étude ont été effectués à l’aide de sealedenvelope.com et ont été fondés sur des données préliminaires.
| Contrôle | Intervention | N par groupe |
| Changement moyen | SD | Changement moyen | SD |
| Douleur clinique (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Fonction (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcIRM (S1M1 à PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tableau 3 : Calcul de la taille de l’échantillon. Calculs liés à la taille de l’échantillon de l’étude.
Les données fonctionnelles fMRIPrep importées ont été lissées à l’aide d’une convolution spatiale avec un noyau gaussien de 8 mm pleine largeur demi-maximum (FWHM) (voir la figure 3) montre la sortie fonctionnelle de fMRIPrep normalisée dans l’espace de modèle MNI152NLin2009cAsym (à gauche) et l’image fonctionnelle lissée de CONN Toolbox (à droite). Cela se traduit par une augmentation du rapport signal/bruit, ce qui améliore la détection des signaux dépendants du niveau d’oxygène dans le sang (BOLD).

Figure 3 : Sujet unique comparant une image fonctionnelle non lissée dans l’espace MNI (à gauche) à son homologue lissée à 8 mm FWHM. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Les données ont ensuite été débruitées à l’aide d’un pipeline de débruitage standard74, y compris la régression des effets confusionnels potentiels caractérisés par des séries chronologiques de matière blanche (5 composantes de bruit CompCor), des séries chronologiques LCR (5 composantes de bruit CompCor), des paramètres de mouvement et leurs dérivées de premier ordre (12 facteurs)75, des balayages aberrants (inférieurs à 295 facteurs)48 et des tendances linéaires (2 facteurs) au sein de chaque exécution fonctionnelle. suivi du filtrage de fréquence passe-bande de la série temporelle BOLD76 entre 0,008 Hz et 0,09 Hz. Les composantes de bruit CompCor49,77 dans la substance blanche et le LCR ont été estimées en calculant le signal BOLD moyen ainsi que les composantes principales les plus grandes orthogonales à la moyenne BOLD, les paramètres de mouvement et les balayages aberrants dans les masques de segmentation érodés de chaque sujet (voir Figure 4). D’après le nombre de termes de bruit inclus dans cette stratégie de débruitage, les degrés de liberté effectifs du signal BOLD après le débruitage ont été estimés entre 33 et 240,6 (moyenne de 173,4) pour tous les sujets. Le débruitage a entraîné une réduction des bruits physiologiques et autres bruits étrangers par rapport aux données, ce qui aurait pu entraîner des effets confondants.

Figure 4 : Contrôles de qualité. Graphiques de contrôles de qualité de CONN Toolbox affichant les effets du débruitage sur la connectivité fonctionnelle (FC), le signal global moyen et le mouvement maximal. Le graphique le plus haut de la note affiche les données d’un seul sujet et d’une seule session, les graphiques B et C sont les résultats du débruitage au niveau du groupe. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Des analyses au niveau du groupe ont été effectuées à l’aide d’un modèle linéaire général (GLM)78. Pour chaque voxel individuel, un GLM distinct a été estimé, avec des mesures de connectivité de premier niveau à ce voxel en tant que variables dépendantes (un échantillon indépendant par sujet et une mesure par tâche ou condition expérimentale, le cas échéant), et des groupes ou d’autres identificateurs au niveau du sujet en tant que variables indépendantes. Les hypothèses au niveau du voxel ont été évaluées à l’aide de statistiques paramétriques multivariées avec des effets aléatoires entre les sujets et une estimation de la covariance de l’échantillon entre plusieurs mesures. Les inférences ont été effectuées au niveau des grappes individuelles (groupes de voxels contigus). Les inférences au niveau du cluster étaient basées sur des statistiques paramétriques de la théorie gaussienne des champs aléatoires62,79. Les résultats ont été seuillés à l’aide d’une combinaison d’un seuil de taille de voxel p formant un agrégat << 0,005 et d’un seuil de taille de cluster de80 corrigé par erreur familiale corrigée par famille. En suivant ces étapes, on obtient des valeurs de connectivité fonctionnelle au niveau du groupe comparant le CES et les conditions fictives en fonction des régions d’intérêt (ROI). Ces résultats peuvent être visualisés de multiples façons via l’interface graphique de l’explorateur de résultats dans CONN Toolbox. Pour voir une visualisation de l’affichage du volume d’un exemple de résultats au niveau du groupe, la zone rouge indiquant les régions avec la plus grande connectivité positive avec les zones d’intérêt et le bleu sont les régions avec une plus grande connectivité négative avec le retour sur investissement, consultez la figure 5.

Figure 5 : Affichage du volume des graines post-cingulaires. L’image montre une plus grande connectivité positive avec le gyrus cingulaire antérieur et (en rouge) et une plus grande connectivité négative avec le précuneus (en bleu) pour la condition True par rapport au Sham, pFWEc < 0,05. Ce chiffre représente des données partielles au niveau du groupe (n = 34). Veuillez cliquer ici pour voir une version agrandie de cette figure.
Fichier supplémentaire 1 : Instructions de commande de l’appareil CES Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 2 : Journal de l’appareil CES Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 3 : fMRIPrep Boilerplate Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 4 : Instructions CES CONN Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 5 : Tracés du code R du CES Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 6 : Code R CES eddy-qc Anova Veuillez cliquer ici pour télécharger ce fichier.
Figure supplémentaire 1 : Classeur CONN. (A) Covariables de deuxième niveau. (B) Conditions expérimentales. (C) Débruitage. (D) Résultats. E) Données structurelles. Veuillez cliquer ici pour télécharger cette figure.