$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Ethische overwegingen en studiepopulatie
Deze retrospectieve studie werd goedgekeurd door de ethische commissie van de instelling. Geïnformeerde toestemming werd opgegeven omdat alle beschermde gezondheidsinformatie werd geanonimiseerd. De onderzoekspopulatie bestond uit patiënten die waren teruggevonden uit de beeldvormingsdatabase van het First People's Hospital of Nantong, die tussen januari 2022 en december 2023 een lumbale wervelkolom-MRI ondergingen vanwege lage rugpijn. Klinische kenmerken van de patiënten werden ook verzameld (Tabel 1).
Inclusie- en uitsluitingscriteria
De huidige diagnostische standaard voor DLBP volgt de discografiemethode van de International Association for the Study of Pain uit 1995, die pijn veroorzaakt door verhoogde druk maar niet breed wordt geaccepteerd vanwege de hoge invasieve aard8. Daarom nam deze studie een restrictieve intraoperatieve test voor schijfpijn aan als diagnostische methode voor DLBP, in overeenstemming met discografieprincipes en met minimale verdere schijfschade en complicaties31. De procedure betrof patiënten in buikligging onder spinale anesthesie, waarbij de doellendenwervelschijf werd gelokaliseerd met behulp van C-arm röntgenbegeleiding. Een 18 G (<22 G) naald werd via een posterolaterale benadering in de centrale nucleus pulposus ingebracht, waarbij zenuwwortels en de durale zak werden vermeden. Na bevestiging van naaldplaatsing werd zoutoplossing of niet-ionisch contrast (bijv. iohexol) geïnjecteerd met een snelheid van niet meer dan 0,5 mL/min, met een druk onder de 50 psi en een totaal injectievolume van minder dan 3 mL, om verhoogde intradiscale druk te simuleren en bekende symptomen van lage rugpijn op te roepen, met een visuele analoge schaal (VAS) score ≥ 7.
DLBP-groepsinklusiecriteria: De hier opgenomen patiënten ondergingen een MRI-onderzoek, hebben terugkerende lage rugpijn van meer dan 3 maanden, met mislukte conservatieve behandeling, met of zonder gevoelloosheid in de onderledemaat of uitstralende pijn, en zijn positief voor de intraoperatieve test voor schijfpijnproccatie.
Inclusiecriteria voor de niet-DLBP-groep: Patiënten die hier zijn opgenomen, ondergingen een MRI; geen geschiedenis van lage rugpijn binnen 3 maanden of gezonde personen die een lichamelijk onderzoek ondergaan; geen MRI-afwijkingen; en een gestandaardiseerde beoordeling hebben met Oswestry Disability Index (ODI) <10 en VAS ≤2.
Uitsluitingscriteria: Sluit patiënten uit met andere oorzaken van lage rugpijn, zoals significante hernia die zenuwen samendrukt, breuken, wervelkolominfecties, spondylolisthesis, tumoren, osteoporose of metabole botziekten; Voorgeschiedenis van de operatie voorafgaand aan het onderzoek; onduidelijke of slechte kwaliteit afbeeldingen; onvermogen om het interessegebied (ROI) te identificeren.
Uiteindelijk werden 243 patiënten opgenomen, bestaande uit 81 DLBP-patiënten en 162 controles.
MRI-parameters
Alle patiënten die in deze studie waren opgenomen, ondergingen 3,0 T MRI-scans, met sequenties die sagittale en axiale T1-gewogen beeldvorming (T1WI) en T2-gewogen beeldvorming (T2WI) omvatten. In de studie werden drie verschillende MRI-machines gebruikt: Siemens Verio, Siemens Prisma en Philips Ingenia CX. De scanparameters werden als volgt ingesteld: voor sagittale T2WI varieerde TR van 2000 tot 4597 ms, TE van 90 tot 120 ms, slicedikte van 4,0 tot 4,8 mm, met 15 slices inbegrepen, bandbreedte van 250 Hz tot 340 Hz, matrixgrootte van 384 × 384 of 512 × 512, faseveldpercentage van 100%, en een uitleesveld van 300 mm.
HIZ-beeldvormingsfuncties meten en statistieken
Deze studie gebruikte 3D Slicer-software (versie 5.6.1, https://download.slicer.org/?version=5.6.1) om handmatig lumbale MRI T2-gewogen beelden van 243 patiënten en controles te analyseren. De gegevensverwerkingsworkflow was als volgt: Afbeeldingen werden geïmporteerd met de functie Add DICOM Data in 3D Slicer. De metingen werden uitgevoerd met behulp van de Markups-functie door twee promovendi in wervelkolomonderzoek onder begeleiding van radiologen en wervelkolomchirurgen met meer dan 10 jaar klinische ervaring. Verschillen in metingen werden opgelost door overleg met de twee artsen. Alle meetwaarden werden door twee statistici gemiddeld om nauwkeurigheid te waarborgen. Op sagittale beelden werden de twee onderling loodrechte maximale diameters van de hoogintensiteitszone (HIZ) gemeten, gedefinieerd als Hizh (bijna-verticale richting) en Hizw (bijna-horizontale richting). De Segmentatiefunctie werd gebruikt om het HIZ-gebied af te bakenen, en het gebied van de meest prominente HIZ op het sagittale vlak werd berekend, aangeduid als Hizarea. Op axiale beelden werd de maximale lengte van de HIZ gemeten, genaamd Hizl. Daarnaast werden drie binaire variabelen gedefinieerd: Hiz (aanwezigheid van HIZ), Overig (aanwezigheid van multisegmentale HIZ) en Positie (of HIZ de posterieure middellijn kruist) (Tabel 2).
Radiomics kenmerken extractie en standaardisatie
Om de specifieke werking van Resampling Scalar Volume te integreren in de context van het verwerken van lumbale MRI T2-gewogen beelden van 243 patiënten en controles met behulp van de 3D Slicer, werden de volgende gedetailleerde stappen gevolgd: Het proces begon met het importeren van lumbale MRI T2-gewogen beelden in de 3D Slicer. Om consistentie te waarborgen en heterogeniteitsbias te verminderen, werden alle beelden opnieuw gesampled tot een voxelgrootte van 0,6 × 0,6 × 0,6 mm met behulp van de Resample Scalar Volume-module . De specifieke stappen waren: Ga naar de sectie Modules en selecteer Resample Scalar Volume. Onder Parameter Set zorgt u ervoor dat Resample Scalar Volume is geselecteerd. In de sectie Resampling Parameters stelt u de afstand in op 0,6, 0,6, 0,6 om de doel-voxeldimensies te definiëren. Kies een geschikte interpolatiemethode uit opties zoals lineair, dichtstbijzijnde buur, bspline, hamming, cosinus, welch, lanczos of blackman, met lineair als standaard. Voor Output Volume selecteer of maak je een nieuw volume om de geresamplede data op te slaan. Na het verifiëren van de instellingen klik je op Toepassen om het herbemonsteringsproces uit te voeren.
Na het opnieuw bemonsteren voerden twee promovendi wervelkolomonderzoek, begeleid door radiologen en wervelkolomchirurgen met meer dan 10 jaar ervaring, een semi-automatische afbakening uit van het interessegebied (ROI). Onenigheid werd opgelost door overleg. De specifieke stappen waren: Selecteer de laag (opnieuw gesampled nieuw volume), maak een nieuwe segmentatielaag aan en gebruik de Draw-functie voor laag-voor-laag afbakening. Gebruik de functie Vullen tussen slices voor het vullen van lagen tussen lagen. Glad de gevormde ROI met de mediaan gladmakingsmethode met een kernelgrootte van 3,0 mm, 5 x 5 pixels. Met behulp van de PyRadiomics-bibliotheek werden 107 radiomische kenmerken geëxtraheerd uit de herbemonsterde beelden, waaronder vormkenmerken, eerste-orde statistische kenmerken, grijs-niveau co-occurent-matrix (GLCM), grijs-niveau run-length matrix (GLRLM), grijs-niveau grootte zone-matrix (GLSZM), aangrenzende grijs-toon verschil matrix (NGTDM) en grijs-niveau afhankelijkheidsmatrix (GLDM) kenmerken (Aanvullende Tabel 1). Ten slotte werden de geëxtraheerde featuregegevens gestandaardiseerd met behulp van de Z-score methode om hun natuurlijke bereik om te zetten in een gestandaardiseerd bereik.
Groeperen op basis van kenmerken
De groepen werden als volgt toegewezen:
d0 (basisgroep): Opgenomen klinische kenmerken, n(d0) = 5.
d1 (Basis finetunegroep): Inclusief klinische kenmerken en HIZ-beeldvormende kenmerken, n(d1) = 12.
d2 (modelgroep): Inclusief klinische kenmerken en radiomische kenmerken, n(d2) = 112.
D (Model finetune groep): Inclusief klinische kenmerken, HIZ-beeldvormingsfuncties en radiomische kenmerken, n(D) = 119.
Gegevenslezen en voorverwerking
De gegevens voor de respectievelijke groepen werden gelezen met de read_excel functie in R-software (versie 4.3.1, https://www.r-project.org/, platform: x86_64-w64-mingw32/x64 [64-bit]) met UTF-8-codering (systeemstandaard), die de meeste geschreven talen wereldwijd ondersteunt. De select-functie werd gebruikt om de doelvariabele te scheiden van de featurevariabelen.
Selectie van functies
Om reproduceerbaarheid van datasplitsing te waarborgen, volgden de vier datasets (D, d2, d1, d0) dezelfde verwerkingsworkflow, met een vaste willekeurige seed van 80. Een labelkolom werd toegevoegd aan het dataframe, omgezet naar een factortype, en willekeurige labels (1 en 0) werden gegenereerd met behulp van een binomiale verdeling om de data te splitsen in trainings- en testsets in een verhouding van 8:2 (80% kans voor de trainingsset). Vanwege het beperkte aantal functies hoefden d0- en d1-groepen geen feature-selectie te ondergaan, terwijl d2- en D-groepen feature-selectie ondergingen:
Eerst werd de Mann-Whitney U-test toegepast op de trainingsset om kenmerken met p-waarden <0,05 te selecteren. Vervolgens werd lasso-regressie uitgevoerd met 10-voudige kruisvalidatie en L1-regularisatie om de optimale strafcoëfficiënt te bepalen en kenmerken uit één enkele iteratie te selecteren. Er werd een willekeurige seed van 1 tot 100 gezet, en de bovenstaande stappen (datasplitsing, Mann-Whitney U-test en Lasso-regressie) werden herhaald in een lus. Kenmerken die meer dan 50 keer in 100 iteraties voorkwamen, werden geselecteerd voor latere modellering om betrouwbaarheid in modelconstructie en diagnostische voorspelling te waarborgen. Modelleerkenmerken voor elke groep worden vermeld in Tabel 3.
Rasterzoektocht en modelafstemming
Om de modellen te optimaliseren, werd grid search toegepast om systematisch verschillende hyperparametercombinaties te verkennen en te evalueren. Door iteratief afstellen en prestatiebeoordeling identificeerde deze studie de set hyperparameters die de beste ROC AUC-waarden leverden voor de trein- en testsets, waardoor de prestaties van het model werden geoptimaliseerd.
Modelontwikkeling en beoordeling
Alle modelontwikkeling en data-analyse werden uitgevoerd in R (versie 4.3.1). Deze studie paste verschillende machine learning-algoritmen toe, waaronder Random Forest (RF), Support Vector Machine (SVM), Decision Tree (TREE), K-Nearest Neighbors (KNN) en Logistic Regression (LOG), om 20 modellen te ontwikkelen over de vier groepen (d0, d1, d2 en D) voor het voorspellen van DLBP-diagnose. De modelprestaties werden geëvalueerd met de volgende maatstaven: ROC AUC, PR AUC, nauwkeurigheid, gevoeligheid, specificiteit, positieve voorspellende waarde (PPV), negatieve voorspellende waarde (NPV) en F1-score.
De ROC AUC- en PR-AUC-waarden werden direct gegenereerd uit de ROC- en PR-curves, terwijl de overige meetwaarden (nauwkeurigheid, gevoeligheid, specificiteit, PPV, NPV en F1-score) werden berekend en statistisch geanalyseerd met behulp van overeenkomstige functies in R. Van de 20 ontwikkelde modellen selecteerde deze studie 8 representatieve modellen voor beoordeling op basis van bovenstaande prestatie-indicatoren.
De formules voor elke metriek zijn als volgt:




Q

TP: Echt positief; TN: Echt negatief; FP: Vals-positief; FN: Vals-negatief.
SHAP-interpretabiliteitsanalyse
Het Random Forest-model in de D-groep toonde de beste prestaties in modelbeoordeling. Om de voorspellingen van het model verder te interpreteren, heeft deze studie het model opnieuw opgebouwd in Python (versie 3.7.9) en SHAP-interpretabiliteitsanalyse uitgevoerd. Om consistentie te behouden, waren de train-hyperparameters van het Random Forest-model in Python identiek aan die in R, en was de dataverwerkingsworkflow ook consistent. Specifiek gebruikten we de boommodelinterpreter om de bijdrage van elk kenmerk aan de diagnostische voorspellingen van het model in de trein- en testsets te analyseren. Voor betere visualisatie hebben we SHAP-waardeverdelingsplots, feature-importance-plots en SHAP-krachtplots voor individuele voorspellingen gegenereerd.
Statistische analyse
Alle data-analyses werden uitgevoerd met R (versie 4.3.1, https://www.r-project.org/, platform: x86_64-w64-mingw32/x64 (64-bit)) en Python (versie 3.7.9, https://www.python.org/downloads/release/python-379/). Continue variabelen werden beschreven als gemiddelde ± standaarddeviatie, terwijl categorische variabelen werden uitgedrukt als frequentie en percentage. De vergelijking tussen DLBP en Non-DLBP werd geanalyseerd met behulp van de U-test, waarbij een P-waarde van minder dan 0,05 als statistisch significant werd beschouwd.