$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Considerações éticas e população estudada
Este estudo retrospectivo foi aprovado pelo comitê de ética institucional. O consentimento informado foi renunciado, pois todas as informações de saúde protegidas foram anonimizadas. A população do estudo consistiu em pacientes retirados do banco de dados de imagens do First People's Hospital de Nantong, que realizaram ressonância magnética da coluna lombar devido a dores lombares entre janeiro de 2022 e dezembro de 2023. As características clínicas dos pacientes também foram coletadas (Tabela 1).
Critérios de inclusão e exclusão
O padrão diagnóstico atual para DLBP segue o método de discografia da International Association for the Study of Pain de 1995, que provoca dor por meio do aumento da pressão, mas não é amplamente aceito devido à sua altainvasividade 8. Portanto, este estudo adotou um teste restritivo de provocação de dor discal intraoperatório como método diagnóstico para DLBP, consistente com os princípios da discografia e minimizando danos e complicações discaisadicionais 31. O procedimento envolveu pacientes em posição prona sob anestesia espinhal, com o disco lombar alvo localizado usando guiagem de raio-X em C-arm. Uma agulha de 18 G (<22 G) foi inserida por meio de uma abordagem posterolateral no núcleo pulposo central, evitando as raízes nervosas e o saco dural. Após confirmação da colocação da agulha, foi injetado soro fisiológico ou contraste não iônico (por exemplo, iohexol) a uma taxa não superior a 0,5 mL/min, com pressão abaixo de 50 psi e volume total de injeção inferior a 3 mL, para simular o aumento da pressão intradiscal e provocar sintomas familiares de dor lombar, com uma escala visual analógica (VAS) ≥ 7.
Critérios de inclusão do grupo DLBP: Os pacientes aqui incluídos passaram por exame de ressonância magnética, apresentam dor lombar recorrente por mais de 3 meses, com tratamento conservador fracassado, com ou sem dormência no membro inferior ou dor irradiada, e são positivos para o teste de provocação de dor discal intraoperatória.
Critérios de inclusão em grupos sem DLBP: Pacientes aqui incluídos passaram por ressonância magnética; não têm histórico de dor lombar em 3 meses ou estão pessoas saudáveis passando por exame físico; não apresentam anomalias na ressonância magnética; e possuem uma avaliação padronizada com o Índice de Deficiência (ODI) de Oswestry (ODI) <10 e VAS ≤2.
Critérios de Exclusão: Exclua pacientes que tenham outras causas de dor lombar, como hérnia de disco significativa que comprime nervos, fraturas, infecções na coluna, espondilolistese, tumores, osteoporose ou doenças metabólicas ósseas; histórico de cirurgia antes do exame; imagens pouco claras ou de baixa qualidade; incapacidade de identificar a região de interesse (ROI).
No total, foram incluídos 243 pacientes, compreendendo 81 pacientes com DLBP e 162 controles.
Parâmetros da ressonância magnética
Todos os pacientes incluídos neste estudo passaram por ressonâncias magnéticas de 3,0 T, utilizando sequências que incluíam imagem ponderada sagital e axial T1 (T1WI) e imagem ponderada T2 (T2WI). Três máquinas diferentes de ressonância magnética foram usadas no estudo: Siemens Verio, Siemens Prisma e Philips Ingenia CX. Os parâmetros de varredura foram definidos da seguinte forma: para T2WI sagital, TR variou de 2000 a 4597 ms, TE de 90 a 120 ms, espessura de corte de 4,0 a 4,8 mm, com 15 fatias incluídas, largura de banda de 250 Hz a 340 Hz, tamanho da matriz de 384 × 384 ou 512 × 512, percentual de campo de visão de fase de 100%, e um campo de visão de leitura de 300 mm.
Medição de características e estatísticas de imagem do HIZ
Este estudo utilizou o software 3D Slicer (versão 5.6.1, https://download.slicer.org/?version=5.6.1) para analisar manualmente imagens ponderadas em T2 de ressonância magnética lombar de 243 pacientes e controles. O fluxo de trabalho de processamento de dados era o seguinte: as imagens eram importadas usando a função Adicionar Dados DICOM no 3D Slicer. As medições foram realizadas usando a função Markups por dois estudantes de pós-graduação em pesquisa da coluna sob a orientação de radiologistas e cirurgiões de coluna com mais de 10 anos de experiência clínica. As discrepâncias nas medições foram resolvidas por meio de consulta com os dois médicos. Todos os valores das medições foram medidos por dois estatísticos para garantir a precisão. Em imagens sagitais, os dois diâmetros máximos mutuamente perpendiculares da zona de alta intensidade (HIZ) foram medidos, definidos como Hizh (direção quase vertical) e Hizw (direção quase horizontal). A função de segmentação era usada para delimitar a região HIZ, e a área do HIZ mais proeminente no plano sagital era calculada, denotada como Hizarea. Em imagens axiais, o comprimento máximo do HIZ foi medido, chamado Hizl. Além disso, três variáveis binárias foram definidas: Hiz (presença de HIZ), Outra (presença de HIZ multisegmental) e Posição (se o HIZ cruza a linha média posterior) (Tabela 2).
Radiômica apresenta extração e padronização
Para integrar a operação específica do Volume Escalar de Reamostragem no contexto do processamento de imagens ponderadas em T2 de ressonância magnética lombar de 243 pacientes e controles usando o 3D Slicer, foram seguidos os seguintes passos detalhados: O processo começou com a importação de imagens ponderadas em T2 para a ressonância magnética lombar para o 3D Slicer. Para garantir consistência e reduzir o viés de heterogeneidade, todas as imagens foram reamostradas para um tamanho voxel de 0,6 × 0,6 × 0,6 mm usando o módulo Resample Scalar Volume . Os passos específicos foram: navegar até a seção Módulos e selecionar Reamostrar Volume Escalar. Em Conjunto de Parâmetros, certifique-se de que o Volume Escalar de Reamostragem esteja selecionado. Na seção Parâmetros de Reamostragem , defina o espaçamento para 0,6, 0,6, 0,6 para definir as dimensões voxel alvo. Selecione um método de Interpolação apropriado entre opções como linear, vizinho mais próximo, bspline, hamming, cosseno, welch, lanczos ou blackman, tendo linear como padrão. Para Volume de Saída, selecione ou crie um novo volume para armazenar os dados reamostrados. Após verificar as configurações, clique em Aplicar para executar o processo de reamostragem.
Após a reamostragem, dois estudantes de pós-graduação em pesquisa da coluna, orientados por radiologistas e cirurgiões com mais de 10 anos de experiência, realizaram a delimitação semiautomática da região de interesse (ROI). As discrepâncias foram resolvidas por meio de consulta. Os passos específicos foram: selecionar a camada (reamostrar o novo volume), criar uma nova camada de segmentação e usar a função Desenhar para delimitação camada por camada. Use a função Preenchimento entre fatias para o recheio entre camadas. Suavize o ROI formado usando o método de suavização mediana com tamanho de grão de 3,0 mm, 5 x 5 pixels. Usando a biblioteca PyRadiomics, 107 características radiômicas foram extraídas das imagens reamostradas, incluindo características de forma, características estatísticas de primeira ordem, matriz de coocorrência de nível cinza (GLCM), matriz de comprimento de execução em nível cinza (GLRLM), matriz de zona de tamanho em nível cinza (GLSZM), matriz vizinha de diferença de tons de cinza (NGTDM) e características da matriz de dependência de nível cinza (GLDM) (Tabela Suplementar 1). Por fim, os dados extraídos de características foram padronizados usando o método Z-score para transformar seu alcance natural em um intervalo padronizado.
Agrupamento baseado em características
Os grupos foram distribuídos da seguinte forma:
d0 (Grupo base): Características clínicas incluídas, n(d0) = 5.
d1 (Grupo de ajuste fino-base): Características clínicas incluídas e características de imagem do HIZ, n(d1) = 12.
d2 (Grupo modelo): Características clínicas e radiômicas incluídas, n(d2) = 112.
D (Grupo de ajuste fino modelo): Características clínicas incluídas, características de imagem HIZ e características de radiômica, n(D) = 119.
Leitura e pré-processamento de dados
Os dados dos respectivos grupos eram lidos usando a função read_excel no software R (versão 4.3.1, https://www.r-project.org/, plataforma: x86_64-w64-mingw32/x64 [64 bits]) com codificação UTF-8 (padrão do sistema), que suporta a maioria das linguagens escritas globalmente. A função select era usada para separar a variável alvo das variáveis de características.
Seleção de recursos
Para garantir a reprodutibilidade da divisão dos dados, os quatro conjuntos de dados (D, d2, d1, d0) seguiram o mesmo fluxo de processamento de processamento, com uma semente aleatória fixa de 80. Uma coluna de rótulos foi adicionada ao quadro de dados, convertida para um tipo fator, e rótulos aleatórios (1 e 0) foram gerados usando uma distribuição binomial para dividir os dados em conjuntos de treinamento e teste em uma proporção de 8:2 (80% de probabilidade para o conjunto de treinamento). Devido ao número limitado de recursos, os grupos d0 e d1 não exigiam seleção de características, enquanto os grupos d2 e D passavam por seleção de características:
Primeiro, o teste Mann-Whitney U foi aplicado ao conjunto de treinamento para selecionar características com valores p <0,05. Em seguida, foi realizada regressão Lasso com validação cruzada 10 vezes e regularização L1 para determinar o coeficiente de penalidade ótimo e selecionar características de uma única iteração. Uma semente aleatória de 1 a 100 foi definida, e os passos acima (divisão de dados, teste Mann-Whitney U e regressão Lasso) foram repetidos em um ciclo. Recursos que aparecem mais de 50 vezes em 100 iterações foram selecionados para modelagem subsequente para garantir confiabilidade na construção do modelo e na previsão diagnóstica. As características de modelagem para cada grupo estão listadas na Tabela 3.
Busca em grade e ajuste de modelos
Para otimizar os modelos, a busca em grade foi aplicada para explorar e avaliar sistematicamente diferentes combinações de hiperparâmetros. Por meio de ajustes iterativos e avaliação de desempenho, este estudo identificou o conjunto de hiperparâmetros que forneceram os melhores valores de AUC ROC para os conjuntos de trem e teste, otimizando assim o desempenho do modelo.
Desenvolvimento e avaliação do modelo
Todo o desenvolvimento de modelos e análise de dados foram realizados em R (versão 4.3.1). Este estudo aplicou vários algoritmos de aprendizado de máquina, incluindo Random Forest (RF), Support Vector Machine (SVM), Decision Tree (TREE), K-Nearest Neighbors (KNN) e Logistic Regression (LOG), para desenvolver 20 modelos nos quatro grupos (d0, d1, d2 e D) para prever o diagnóstico de DLBP. O desempenho do modelo foi avaliado usando as seguintes métricas: ROC AUC, PR AUC, precisão, sensibilidade, especificidade, valor preditivo positivo (PPV), valor preditivo negativo (NPV) e escore F1.
Os valores ROC AUC e PR AUC foram gerados diretamente a partir das curvas ROC e PR, enquanto as métricas restantes (precisão, sensibilidade, especificidade, PPV, NPV e pontuação F1) foram calculadas e estatisticamente analisadas usando funções correspondentes em R. Entre os 20 modelos desenvolvidos, este estudo selecionou 8 modelos representativos para avaliação com base nas métricas de desempenho acima.
As fórmulas para cada métrica são as seguintes:




Q

TP: Verdadeiro Positivo; TN: Verdadeiro Negativo; FP: Falso Positivo; FN: Falso Negativo.
Análise de interpretabilidade SHAP
O modelo de Floresta Aleatória no grupo D demonstrou o melhor desempenho na avaliação do modelo. Para interpretar melhor as previsões do modelo, este estudo reconstruiu o modelo em Python (versão 3.7.9) e realizou análise de interpretabilidade SHAP. Para manter a consistência, os hiperparâmetros de treino do modelo Random Forest em Python eram idênticos aos usados em R, e o fluxo de trabalho de processamento de dados também era consistente. Especificamente, usamos o interpretador de modelo em árvore para analisar a contribuição de cada característica para as previsões diagnósticas do modelo nos conjuntos de treinamento e teste. Para melhor visualização, geramos gráficos de distribuição de valores SHAP, gráficos de importância de características e gráficos de força SHAP para previsões individuais.
Análise estatística
Todas as análises de dados foram realizadas usando R (versão 4.3.1, https://www.r-project.org/, plataforma: x86_64-w64-mingw32/x64 (64 bits)) e Python (versão 3.7.9, https://www.python.org/downloads/release/python-379/). Variáveis contínuas foram descritas como média ± desvio padrão, enquanto variáveis categóricas foram expressas como frequência e porcentagem. A comparação entre DLBP e não-DLBP foi analisada usando o teste U, com um valor P menor que 0,05 considerado estatisticamente significativo.