$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Em termos de resultados de recrutamento, os participantes foram recrutados principalmente por meio de cartas de recrutamento e telefonemas de acompanhamento com base nos regulamentos descritos do Atlanta VA Healthcare System. A equipe do estudo recrutou um total de 50 participantes, comprovando a eficácia dos métodos usados para atingir a meta de recrutamento (ver Figura 2). O uso dos novos critérios diagnósticos clínicos de fibromialgia permitiu que a equipe do estudo rastreasse adequadamente os indivíduos que não atendiam aos critérios de fibromialgia32. Perguntar a possíveis participantes sobre um diagnóstico de fibromialgia não é uma medida tão robusta quanto a triagem adicional e poderia ter levado a visitas iniciais agendadas incorretamente ou à participação no estudo. Quarenta e oito participantes foram randomizados em grupos ativos e simulados; Dois participantes foram excluídos por não atenderem aos critérios de elegibilidade com base no teste de inclusão para o estudo.

Figura 2: Fluxograma de recrutamento. Um relatório e fluxograma de recrutamento, randomização e alocação de intervenção do estudo. Clique aqui para ver uma versão maior desta figura.
Três desfechos foram usados nos cálculos do tamanho da amostra: DVPRS (dor clínica), teste de levantar da cadeira de 30 s (função) e conectividade DMN-SMN (rs-fcMRI). Todos os cálculos de potência foram baseados em dados preliminares. As alterações clínicas da dor usando DVPRS foram escolhidas como o desfecho primário de interesse. O teste de levantar da cadeira de 30 s (30sCST) foi escolhido como o resultado funcional representativo, pois exibiu a menor mudança entre os grupos. A conectividade DMN-SMN foi escolhida como o desfecho secundário de interesse como um biomarcador de neuroimagem para dor clínica e resposta ao tratamento. As análises do tamanho da amostra foram conduzidas assumindo um poder de significância de 1% e 80% (2 amostras, 1 lado).
As sementes para esta análise foram escolhidas com base em dados preliminares, bem como na literatura sobre fibromialgia, dor e CES 16,24,26,63,64,65,66,67,68,69. Com base nos dados preliminares, a média (DP) da mudança pós-tratamento na conectividade do córtex sensório-motor primário esquerdo (L-S1M1) para o córtex cingulado posterior esquerdo (L-PCC) é de 0,041 (0,079) para o grupo de tratamento e -0,026 (0,049) para o grupo de tratamento padrão; o tamanho do efeito de diferença entre os grupos observado é de 1,03 (ver Tabela 2 16,21,26,62,63,64,70,71,72,73). O estudo precisaria de 20 indivíduos em cada grupo para atingir 80% de poder para detectar a diferença entre o grupo CES e o grupo de tratamento padrão em sua mudança pós-tratamento na conectividade de L-S1M1 para L-PCC no nível de significância de 0,01 usando um teste t bilateral, assumindo que o tamanho do efeito da diferença entre os grupos é de 1,03, conforme observado nos dados piloto. Embora o estudo piloto tenha observado uma taxa de atrito de 17% para 12 indivíduos em nosso estudo anterior de neuromodulação auricular (todos completaram sua ressonância magnética de acompanhamento, mas dois foram perdidos no acompanhamento em 8 semanas e 12 semanas), a fim de manter uma estimativa conservadora para cálculos de tamanho de amostra, este estudo assumiu uma taxa de atrito de 20%. Com uma expectativa de 20% de atrito na visita pós-tratamento, o estudo precisou recrutar 20/0,8 = 25 indivíduos por grupo.
| Sementes DMN (x,y,z) | Sementes de SMN (x,y,z) | Sementes SN (x, y, z) |
| Córtex Pré-Frontal Medial62,67 | Putâmen Direito64,71 | Córtex pré-frontal dorsolateral direito62 |
| PCC70 direito | M170 esquerdo | Ínsula anterior esquerda62 |
| PCC esquerdo16,67 | Direito M170 | Ínsula anterior direita62 |
| Precuneus71 | Mão DireitaS1 16,72 | Ínsula posterior esquerda64 |
| Mão EsquerdaS1-Mão 16,72 | Ínsula posterior direita63 |
| Tálamo21 | Córtex cingulado anterior dorsal72,73 |
| | Junção temporoparietal direita62 |
Tabela 2: Sementes para análise. Sementes de DMN, SMN e SN escolhidas para análise com base em hipóteses a priori. Cada semente é apresentada com referências à literatura anterior que apóia seus testes em síndromes dolorosas.
Com base em pesquisas anteriores de CES em indivíduos civis com fibromialgia, um total de 50 indivíduos (n = 25 simulados e n = 25 verdadeiros) devem atingir 80% de poder para detectar uma diferença nos escores de dor entre os dois grupos17 (ver Tabela 3). Os cálculos do tamanho da amostra para este estudo foram realizados usando sealedenvelope.com e foram baseados em dados preliminares.
| Controle | Intervenção | N por grupo |
| Mudança média | SD | Mudança média | SD |
| Dor Clínica (DVPRS) | 0.375 | 1.493 | -1.833 | 2.229 | 10 |
| Função (30sCST) | -0.250 | 1.500 | 3.000 | 4.980 | 14 |
| rs-fcMRI (S1M1 para PCC) | -0.026 | 0.049 | 0.041 | 0.079 | 20 |
Tabela 3: Cálculo do tamanho da amostra. Cálculos relacionados ao tamanho da amostra do estudo.
Os dados funcionais importados do fMRIPrep foram suavizados usando convolução espacial com um kernel gaussiano de 8 mm de largura total e meio máximo (FWHM) (veja a Figura 3) mostra a saída funcional do fMRIPrep normalizada no espaço do modelo MNI152NLin2009cAsym (esquerda) e a imagem funcional suavizada do CONN Toolbox (direita). Isso resulta em um aumento na relação sinal-ruído, o que, por sua vez, melhora a detecção de sinais dependentes do nível de oxigênio no sangue (BOLD).

Figura 3: Sujeito único comparando uma imagem funcional não suavizada no espaço MNI (esquerda) com sua contraparte suavizada em 8 mm FWHM. Clique aqui para ver uma versão maior desta figura.
Os dados foram então removidos usando um pipeline de redução de ruído padrão74, incluindo a regressão de potenciais efeitos de confusão caracterizados por séries temporais de substância branca (5 componentes de ruído CompCor), séries temporais CSF (5 componentes de ruído CompCor), parâmetros de movimento e seus derivados de primeira ordem (12 fatores) 75, varreduras discrepantes (abaixo de 295 fatores) 48 e tendências lineares (2 fatores) dentro de cada execução funcional, seguido pela filtragem de frequência passa-banda da série temporal BOLD76 entre 0,008 Hz e 0,09 Hz. Os componentes de ruído CompCor49,77 na substância branca e no LCR foram estimados calculando o sinal BOLD médio, bem como os maiores componentes principais ortogonais à média BOLD, parâmetros de movimento e varreduras discrepantes dentro das máscaras de segmentação erodidas de cada sujeito (ver Figura 4). A partir do número de termos de ruído incluídos nesta estratégia de redução de ruído, os graus efetivos de liberdade do sinal BOLD após a redução de ruído foram estimados em 33 a 240,6 (média de 173,4) em todos os indivíduos. A redução de ruído resultou em uma redução de ruídos fisiológicos e outros ruídos estranhos dos dados que poderiam ter resultado em efeitos de confusão.

Figura 4: Verificações de qualidade. Gráficos de verificações de qualidade do CONN Toolbox exibindo os efeitos da redução de ruído na conectividade funcional (FC), sinal global médio e movimento máximo. Nota (A) o gráfico superior exibe dados para um único assunto e sessão, (B,C) os gráficos B e C são os resultados da redução de ruído em nível de grupo. Clique aqui para ver uma versão maior desta figura.
As análises em nível de grupo foram realizadas usando um Modelo Linear Geral (GLM)78. Para cada voxel individual, um GLM separado foi estimado, com medidas de conectividade de primeiro nível neste voxel como variáveis dependentes (uma amostra independente por sujeito e uma medição por tarefa ou condição experimental, se aplicável) e grupos ou outros identificadores de nível de sujeito como variáveis independentes. As hipóteses em nível de voxel foram avaliadas usando estatísticas paramétricas multivariadas com efeitos aleatórios entre os sujeitos e estimativa de covariância amostral em várias medições. As inferências foram realizadas no nível de clusters individuais (grupos de voxels contíguos). As inferências em nível de cluster foram baseadas em estatísticas paramétricas da teoria Gaussiana de Campos Aleatórios 62,79. Os resultados foram limiares usando uma combinação de um limite de nível de voxel 80. Seguir essas etapas resulta em valores de conectividade funcional em nível de grupo comparando as condições CES e simuladas com base nas regiões de interesse (ROI). Esses resultados podem ser visualizados de várias maneiras por meio da GUI do explorador de resultados no CONN Toolbox. Para ver uma visualização de exibição de volume de um exemplo de resultados em nível de grupo com a área vermelha indicando regiões com maior conectividade positiva com os ROIs e o azul são regiões com maior conectividade negativa com o ROI, consulte a Figura 5.

Figura 5: Exibição do volume da semente pós-cíngula. A imagem exibe maior conectividade positiva com o giro cingulado anterior e (vermelho) e maior conectividade negativa com o precuneus (azul) para a condição True em relação ao Sham, pFWEc < 0,05. Esta figura representa dados parciais em nível de grupo (n = 34). Clique aqui para ver uma versão maior desta figura.
Arquivo suplementar 1: Instruções de pedido de dispositivo CES Clique aqui para baixar este arquivo.
Arquivo suplementar 2: Registro do dispositivo CES Clique aqui para baixar este arquivo.
Arquivo Suplementar 3: Boilerplate do fMRIPrep Clique aqui para baixar este arquivo.
Arquivo Suplementar 4: Instruções CES CONN Clique aqui para baixar este arquivo.
Arquivo Suplementar 5: Plotagens de código CES R Clique aqui para baixar este arquivo.
Arquivo Suplementar 6: Código R CES eddy-qc Anova Clique aqui para baixar este arquivo.
Figura suplementar 1: Pasta de trabalho CONN. (A) Covariáveis de segundo nível. (B) Condições experimentais. (C) Redução de ruído. (d) Resultados. (E) Dados estruturais. Clique aqui para baixar esta figura.