Protokół diagnostyczny oparty na obrazowaniu metodą rezonansu magnetycznego oraz przewidywanie ryzyka wystąpienia lipomatozy znadtwardówkowej w przebiegu zwyrodnieniowej spondylolistezy lędźwiowej

Podsumowanie

Niniejszy protokół opisuje ustandaryzowaną metodę opartą na obrazowaniu metodą rezonansu magnetycznego w celu diagnozowania lipomatozy znadtwardówkowej oraz przewidywania ryzyka jej wystąpienia u pacjentów z degeneracyjnym spondylolistezą lędźwiową.

Streszczenie

Niniejszy protokół przedstawia ustandaryzowane podejście diagnostyczne oparte na obrazowaniu metodą rezonansu magnetycznego (MRI) w ocenie lipomatozy nadtwardówkowej u pacjentów z degeneracyjnym spondylolistezą lędźwiową (DLS) oraz opisuje opracowanie klinicznego modelu przewidywania ryzyka. Kolejni pacjenci z DLS zostali poddani ustandaryzowanej ocenie MRI, obejmującej sekwencje T1 w celu wykrycia lipomatozy nadtwardówkowej, ocenę stopnia nacieku tłuszczowego mięśnia wielodzielnego według skali Goutalliera oraz ocenę radiologiczną. Kohortę badawczą podzielono na grupę modelowania (n = 248) i grupę walidacji (n = 106). Niezależne predyktory zidentyfikowano za pomocą analizy jednowymiarowej, a następnie wielowymiarowej regresji logistycznej binarnej. W celu zapewnienia przejrzystości raportowania, wpływ wieku oraz wskaźnika masy ciała (BMI) wyrażono odpowiednio na każde 10 lat i na każde 5 kg/m² wzrostu. Równanie predykcyjne do użytku przyłóżkowego przedstawiono z wykorzystaniem algebraicznie równoważnych współczynników jednostek surowych, aby wiek (lata) i BMI (kg/m²) można było wprowadzać bezpośrednio. Zidentyfikowano pięć niezależnych predyktorów: wiek (iloraz szans [OR] = 1,510 na każde 10 lat wzrostu), płeć żeńską (OR = 2,354), BMI (OR = 1,874 na każde 5 kg/m² wzrostu), zajęcie segmentu L5 (OR = 3,766) oraz stopień 3-4 w skali Goutalliera (OR = 3,184). Model osiągnął wartości powierzchni pod krzywą charakterystyki ROC wynoszące 0,834 w kohorcie modelowania i 0,815 w kohorcie walidacyjnej. Analizy kalibracji i krzywych decyzyjnych potwierdziły przydatność modelu jako eksploracyjnego narzędzia wspomagania decyzji klinicznych; jednak przed szerokim wdrożeniem klinicznym wymagana jest walidacja zewnętrzna. Niniejszy protokół dostarcza powtarzalnych ram oceny lipomatozy nadtwardówkowej w oparciu o MRI wraz z zintegrowanym narzędziem predykcyjnym, które łączy parametry kliniczne i radiologiczne w celu ułatwienia ustandaryzowanej oceny i zindywidualizowanego przewidywania ryzyka u pacjentów z DLS.

Wprowadzenie

Zwyrodnieniowa spondylolisteza lędźwiowa (DLS) dotyczy 4%–8% populacji ogólnej, przy czym jej rozpowszechnienie gwałtownie wzrasta po piątej dekadzie życia^1,2. Definiowana jako przednie przemieszczenie kręgu względem kręgu pod nim przy zachowaniu nienaruszonego tylnego łuku neuronalnego, DLS często prowadzi do zwężenia kanału kręgowego i chromania neurogennego, co znacznie upośledza mobilność pacjenta, siłę mięśniową oraz jakość życia³. Wśród powiązanych cech obrazowych objaw lipomatozy nadtwardówkowej — charakteryzujący się pasmowym, hiperintensywnym sygnałem tłuszczowym w obrazowaniu rezonansem magnetycznym (MRI) w płaszczyźnie strzałkowej T1 na poziomie przesunięcia — jest istotnym wskaźnikiem radiologicznym, który może ułatwić rozpoznanie choroby i ocenę jej stopnia zaawansowania. Lipomatozę nadtwardówkową tradycyjnie wiązano z podawaniem glikokortykosteroidów, ekspozycją na sterydy wewnątrzoponowe lub zaburzeniami endokrynologicznymi, takimi jak zespół Cushinga^4,5,6. Jednakże coraz częściej jest ona dokumentowana u osób z otyłością lub zespołem metabolicznym bez ekspozycji na sterydy^7,8,9,10, co sugeruje, że czynniki metaboliczne i biomechaniczne odgrywają istotną rolę etiologiczną. Niedawny przegląd systematyczny podsumował różnorodne etiologie i wyniki kliniczne w zgłoszonych przypadkach¹¹, a udokumentowano również poprawę objawów po redukcji masy ciała¹². Opisano wcześniejsze metody stopniowania w oparciu o MRI i ilościowej oceny lipomatozy nadtwardówkowej kręgosłupa, w tym podejścia oparte na stosunku tłuszczu nadtwardówkowego do worka oponowego w przekroju osiowym oraz systemy stopniowania lokalno-regionalnego^13,14. Niemniej jednak nie ustalono jeszcze zstandardyzowanego protokołu MRI przeznaczonego specjalnie do systematycznej oceny lipomatozy nadtwardówkowej u pacjentów z DLS. W związku z tym interpretacja obrazów pozostaje zależna od operatora, brakuje spójnych kryteriów diagnostycznych dostosowanych do tej populacji pacjentów, determinanty kliniczne odróżniające osoby dotknięte tą dolegliwością od osób zdrowych nie zostały w pełni scharakteryzowane, a klinicyści nie dysponują praktycznym narzędziem do szacowania ryzyka podczas wstępnej oceny.

Aby wypełnić te luki, niniejszy protokół ustanawia zestandaryzowane, klinicznie stosowalne ramy oceny lipomatozy nadtwardówkowej w oparciu o MRI u pacjentów z DLS. Protokół integruje parametry kliniczne, zdefiniowane kryteria diagnostyczne MRI ze specyficznymi parametrami akwizycji obrazu, ocenę mięśni przykręgosłupowych z wykorzystaniem dostosowanej do kręgosłupa skali Goutalliera, osiową ocenę potwierdzającą ucisk worka oponowego oraz model prognostyczny do stratyfikacji ryzyka przy łóżku pacjenta. System gradingowy Manjila jest uznawany za istotny punkt odniesienia w zakresie lokalizacji regionalnej, ponieważ ocenia zakres kraniokaudalny w obrazowaniu strzałkowym oraz nasilenie w przekroju osiowym przy użyciu oceny opartej na siatce 3 × 3, co może pomóc w planowaniu chirurgicznym i ułatwić rozróżnienie między grzbietowym, brzusznym i obwodowym rozkładem tłuszczu w przestrzeni nadtwardówkowej^14,15,16. Niniejszy protokół nie ma na celu zastąpienia istniejących systemów gradingowych lipomatozy nadtwardówkowej kręgosłupa; zamiast tego dostarcza on ukierunkowany schemat postępowania dla poziomów przesuniętych związanych z DLS, dokumentując jednocześnie cechy osiowe i lokalno-regionalne, gdy są one istotne klinicznie. Podejście to ma na celu promowanie bardziej powtarzalnej oceny lipomatozy nadtwardówkowej, wspierając zestandaryzowaną ocenę obrazową oraz indywidualną szacunkową ocenę ryzyka klinicznego u pacjentów z DLS.

Protokół

Niniejszy protokół MRI został opracowany dla nadprzewodzącego skanera o natężeniu pola 1,5 Tesla, wyposażonego w dedykowaną 16-kanałową cewkę powierzchniową do obrazowania kręgosłupa. Protokół można dostosować do innych systemów 1,5 lub 3,0 T o równoważnej specyfikacji sprzętowej. Wszystkie procedury zostały przeprowadzone zgodnie z Deklaracją Helsińską i zostały zatwierdzone przez Instytucjonalną Komisję Rewizyjną (Institutional Review Board) Tianjin Union Medical Center (numer zatwierdzenia: 2024-IRB-089).

1. Wybór pacjentów i ocena kliniczna

1. Przed włączeniem pacjentów do badania należy uzyskać zgodę Instytucjonalnej Rady Przeglądowej (Institutional Review Board) oraz Komisji Etycznej. Należy upewnić się, że wszystkie procedury z udziałem ludzi są zgodne z Deklaracją Helsińską. Numer zgody (2024-IRB-089) musi zostać odnotowany we wszystkich dokumentach badania.
   ​UWAGA: Wymóg uzyskania indywidualnej świadomej zgody może zostać zniesiony w przypadku, gdy badanie wykorzystuje w pełni zanonimizowane, istniejące dane kliniczne oraz zestawy danych obrazowych, pod warunkiem że zniesienie to jest zgodne z polityką instytucjonalną oraz przepisami krajowymi regulującymi wtórne wykorzystanie danych klinicznych. Należy uzyskać wyraźną zgodę komisji etycznej na zniesienie tego wymogu. Przed przystąpieniem do rekrutacji pacjentów należy uzyskać zgodę komisji etycznej.
2. Należy przeprowadzić przesiew kolejnych pacjentów zgłaszających się do oddziału chirurgii kręgosłupa z potwierdzoną diagnozą DLS. Należy przeanalizować zdjęcia rentgenowskie w projekcji bocznej w pozycji stojącej lub w pozycjach zgięcia i wyprostu, aby potwierdzić przednie przesunięcie kręgu o ≥3 mm lub ≥10% szerokości trzonu kręgu.
   1. Zweryfikuj zmiany zwyrodnieniowe na odpowiadającym poziomie, w tym degenerację krążka międzykręgowego oraz przerost stawów facetowych¹⁷Potwierdź brak defektu pars interarticularis w obrazach rentgenowskich w projekcji bocznej oraz w obrazowaniu MRI w płaszczyźnie strzałkowej.
   2. Należy zakwalifikować pacjentów w wieku 18–85 lat, u których wykonano standardowe badanie MRI odcinka lędźwiowego kręgosłupa (sekwencje T1-zależne w płaszczyźnie strzałkowej, T2-zależne w płaszczyźnie strzałkowej oraz T2-zależne w płaszczyźnie osiowej) w ciągu 2 tygodni od oceny klinicznej. Należy wykluczyć pacjentów, u których w ciągu ostatnich 6 miesięcy przeprowadzono operację kręgosłupa lub wykonano nadtwardówkowe wstrzyknięcie sterydów.
   3. Zmierz przednie przesunięcie kręgu na zdjęciach rentgenowskich w projekcji bocznej w pozycji stojącej lub w podczas zgięcia i wyprostu, stosując metodę linii tylnej krawędzi trzonu kręgu. Narysuj pionową linię referencyjną wzdłuż tylnej kory trzonu kręgu położonego niżej, a następnie zmierz odległość prostopadłą od tej linii do tylno-dolnego rogu przesuniętego kręgu położonego wyżej. Zapisz wartość przesunięcia w milimetrach oraz, jeśli jest to wymagane, jako procent szerokości przednio-tylnej trzonu kręgu położonego niżej.
3. Wykluczyć pacjentów z niedegeneracyjnym spondylolisthesis (istmicznym, urazowym, patologicznym, pooperacyjnym lub rozwojowym lub z imitacjami anatomicznymi, takimi jak hipoplazja trzonu kręgu L5 i pseudolisthesis L5–S1), stosującymi glikokortykosteroidy w celu >łącznie 3 miesiące w ciągu poprzedniego roku lub potwierdzony endogenny hiperkortyzolizm (zespół Cushinga, gruczolak nadnerczy lub przysadkowy guz wydzielający hormon adrenokortykotropowy).
   1. Należy wykluczyć pacjentów, którzy w ciągu ostatnich 6 miesięcy otrzymywali sterydy anaboliczne lub terapię zastępczą testosteronem. Należy wykluczyć badania MRI o niewystarczającej jakości obrazu, w tym z artefaktami ruchowymi stopnia II lub wyższego bądź z niepełnym zakresem obrazowania odcinka lędźwiowo-krzyżowego.
   2. Należy wykluczyć pacjentów z ciężką niewydolnością krążeniowo-oddechową lub nerkową (klasa IV–V w skali stanu fizycznego American Society of Anesthesiologists), aktywną infekcją lub przeciwwskazaniami do badania MRI.
      ​UWAGA: Przed wejściem do skanera należy przeprowadzić pełny screening bezpieczeństwa MRI. Należy wykluczyć pacjentów z wszczepionymi urządzeniami niekompatybilnymi z MRI.
   3. Oceń artefakty ruchowe przy użyciu czterostopniowej skali jakości obrazu: stopień 0 = brak widocznego artefaktu; stopień I = łagodny artefakt bez upośledzenia diagnostycznego; stopień II = umiarkowany artefakt przesłaniający granice anatomiczne lub punkty orientacyjne pomiaru; oraz stopień III = ciężki artefakt sprawiający, że sekwencja jest nienadająca się do diagnostyki. Wyklucz badania z artefaktami stopnia II lub III w jakiejkolwiek sekwencji wymaganej do klasyfikacji wyniku lub powtórz dotkniętą sekwencję przed włączeniem badania do badania.
4. Wyodrębnić dane kliniczne z bazy danych elektronicznej dokumentacji medycznej (EMR) oraz systemu archiwizacji i przesyłania obrazów (PACS), korzystając ze zstandaryzowanego formularza zbierania danych REDCap wypełnionego przez dwóch przeszkolonych asystentów badawczych.
   1. Zapisz wiek (lata), płeć (mężczyzna/kobieta), wskaźnik masy ciała (BMI; kg/m²)²), oraz przesuniętego segmentu (L3, L4 lub L5).
   2. Należy odnotować czas trwania objawów, występowanie chromania neurogennego oraz wcześniej zastosowane leczenie zachowawcze.
   3. Dwóch przeszkolonych asystentów badawczych powinno niezależnie wyodrębnić zmienne kliniczne z EMR oraz PACS do powielonych formularzy REDCap. Po ekstrakcji danych należy porównać oba formularze i rozwiązać rozbieżności, odwołując się do rekordów źródłowych. Nierozstrzygnięte rozbieżności należy przekazać do rozstrzygnięcia starszemu badaczowi.
   4. Zidentyfikuj brakujące lub niepełne zmienne podczas sprawdzania zakresu i kompletności danych w systemie REDCap. Zweryfikuj brakujące wartości w oparciu o EMR, PACS oraz oryginalne raporty obrazowe. Wyklucz pacjenta z opracowywania modelu tylko wtedy, gdy nie można odzyskać wymaganego predyktora, zmiennej wynikowej lub niezbędnej sekwencji MRI. Nie stosuj imputacji statystycznej dla wymaganych zmiennych modelu.
5. Kwalifikujących się pacjentów należy losowo przydzielić do kohorty modelowania (n = 248) oraz kohorty walidacyjnej (n = 106), stosując stosunek alokacji 7:3 w programie SPSS wersja 27.0 z ziarnem losowym 42.
   1. Pacjentów w kohorcie modelowej należy podzielić na grupę z lipomatozą znadtwardówkową oraz grupę bez lipomatozy, zgodnie z konsensualną diagnozą dwóch starszych radiologów specjalizujących się w obrazowaniu kręgosłupa z >10 lat doświadczenia, zgodnie z protokołem interpretacji MRI opisanym w kroku 3.
   2. W programie SPSS w wersji 27.0 ustaw ziarno losowości na 42, a następnie wybierz Transform (Przekształć) > Oblicz zmienną i wygeneruj jednostajną zmienną losową za pomocą funkcji RV.UNIFORM(0,1). Posortuj przypadki w kolejności rosnącej według tej zmiennej, przypisz pierwsze 70% kwalifikujących się przypadków do kohorty modelowej, a pozostałe 30% do kohorty walidacyjnej, oraz zachowaj zmienną alokacji bez zmian dla wszystkich kolejnych analiz.
   3. Należy zastosować prostą randomizację zamiast randomizacji warstwowej. Nie należy wymuszać równowagi zmiennych demograficznych, radiologicznych ani obrazów MRI w momencie przydzielania do grup.
   4. Po losowym przydziale należy porównać wyjściowe zmienne demograficzne, kliniczne, radiologiczne oraz zmienne z obrazowania MRI pomiędzy kohortą modelującą a kohortą walidacyjną, wykorzystując test t Studenta t-test lub test chi-kwadrat, w zależności od potrzeb. Przypisanie uznaje się za dopuszczalne, gdy nie stwierdzono klinicznie istotnych dysproporcji, a ocenę porównywalności kohort należy przedstawić w Tabela 1.
      UWAGA: Przydział do kohort należy przeprowadzić po zweryfikowaniu wszystkich kryteriów kwalifikowalności.

Zmienna	Kohorta modelowania (n = 248)	Kohorta walidacyjna (n = 106)	t/χ2	Wartość P
Wiek (lata)	67.87 ± 11.39	68.16 ± 11.93	0.216	0.829
Płeć, n (%)			0.002	0.966
Mężczyzna	96 (38.71)	42 (39.62)
Kobieta	152 (61.29)	64 (60.38)
Wskaźnik masy ciała (kg/m²)	24.88 ± 3.62	24.91 ± 3.48	0.077	0.939
Ssegment przesunięty, n (%)			0.564	0.754
L3	31 (12.50)	14 (13.21)
L4	149 (60.08)	67 (63.21)
L5	68 (27.42)	25 (23.58)
Stopień wg Pfirrmanna, n (%)			0.104	0.991
Stopień II	18 (7.26)	8 (7.55)
Stopień III	74 (29.84)	33 (31.13)
Stopień IV	93 (37.50)	38 (35.85)
Stopień V	63 (25.40)	27 (25.47)
Stopień wg Goutalliera, n (%)			3.022	0.388
Stopień 1	78 (31.45)	27 (25.47)
Stopień 2	86 (34.68)	33 (31.13)
Stopień 3	52 (20.97)	28 (26.42)
Stopień 4	32 (12.90)	18 (16.98)
Kąt stawu facetowego dystalnego (°)	55.26 ± 6.42	55.18 ± 6.39	0.111	0.912
Wysięk w stawie facetowym dystalnym (mm)	0.94 ± 0.21	0.93 ± 0.20	0.444	0.657
Wysokość krążka dystalnego (mm)	9.16 ± 1.92	9.20 ± 1.88	0.193	0.847
Kąt stawu facetowego proksymalnego (°)	49.95 ± 6.03	50.03 ± 5.99	0.115	0.908
Wysięk w stawie facetowym proksymalnym (mm)	0.84 ± 0.26	0.83 ± 0.25	0.237	0.813
Wysokość krążka proksymalnego (mm)	8.05 ± 1.52	8.08 ± 1.50	0.156	0.876
Incydencja miednicy (PI, °)	52.18 ± 7.94	52.35 ± 7.85	0.188	0.851
Nachylenie miednicy (PT, °)	24.35 ± 7.17	24.28 ± 7.23	0.080	0.936
Kąt nachylenia kości krzyżowej (SS, °)	44.91 ± 6.83	45.01 ± 6.79	0.129	0.898
Kifoza piersiowa (TK, °)	21.18 ± 7.67	21.25 ± 7.62	0.075	0.940

Tabela 1: Charakterystyka wyjściowa kohort modelowania i walidacji. Zmienne ciągłe przedstawiono jako średnia ± odchylenie standardowe (SD), a zmienne kategoryczne jako liczba (%). Porównano demograficzne, kliniczne, radiologiczne i obrazowanie rezonansem magnetycznym cechy wyjściowe pomiędzy kohortą modelowania (n = 248) a kohortą walidacji (n = 106) w celu oceny porównywalności kohort. Skróty: BMI, wskaźnik masy ciała; PI, incydencja miedniczna; PT, nachylenie miednicy; SS, kąt nachylenia kości krzyżowej; TK, kifoza piersiowa.

2. Protokół akwizycji obrazów MRI

1. Przeprowadź badanie przesiewowe pod kątem bezpieczeństwa MRI przy użyciu standaryzowanego kwestionariusza. Zidentyfikuj przeciwwskazania, w tym ciążę, wszczepione urządzenia elektroniczne, ferromagnetyczne ciała obce, klaustrofobię oraz historię reakcji alergicznych na środki kontrastowe na bazie gadolinu.
   1. Poinstruuj pacjenta, aby usunął wszystkie metalowe przedmioty, w tym biżuterię, zegarki, kolczyki i aparaty słuchowe. Ułóż pacjenta w pozycji leżącej na plecach na stole do MRI.
   2. Umieść miękką poduszkę pod kolanami, aby uzyskać zgięcie wynoszące około 20°–30° i zmniejszyć fizjologiczną lordozę lędźwiową. Ustaw cewkę kręgosłupa centralnie na poziomie kręgu L3, aby zapewnić obrazowanie od T12 do S3.
   3. Zabezpiecz cewkę tak, aby miała bezpośredni kontakt ze skórą, eliminując widoczne przerwy powietrzne. Umieść podkładki piankowe po obu stronach, aby zminimalizować ruchy pacjenta, i zapewnij zgodne z MRI urządzenie do wezwania pomocy w sytuacjach awaryjnych.
   4. Nie podawaj dożylnych środków kontrastowych na bazie gadolinu. Procedura diagnostyczna obejmuje sekwencje bezkontrastowe: sagitalne T1-zależne, sagitalne T2-zależne oraz aksjalne T2-zależne.
2. Nabyj sekwencję lokalizacyjną (scout) w trzech płaszczyznach, stosując sekwencję gradient-echo z TR = 8 ms, TE = 4 ms, grubością warstwy = 8 mm i polem widzenia (FOV) = 400 mm × 400 mm.
   1. Zweryfikuj pełny zakres obrazowania kręgosłupa lędźwiowego od T12/L1 do kości krzyżowej. Skoryguj ułożenie pacjenta, jeśli przejście lędźwiowo-krzyżowe nie jest w pełni widoczne.
3. Nabyj sagitalne obrazy T1-zależne turbo spin-echo, stosując TR = 500 ms, TE = 12 ms, pole widzenia = 280 × 280 mm, macierz = 512 × 256, grubość warstwy = 4 mm, odstęp między warstwami = 0,4 mm oraz NEX = 2.
   1. Wyznacz warstwy równolegle do wyrostków kolczystych. Nabyj minimum 11 warstw sagitalnych, aby zapewnić pełny zakres obrazowania od lewej do prawej strony.
      ​UWAGA: Obrazowanie T1-zależne jest podstawową sekwencją do identyfikacji lipomatozy nadtwardówkowej, ponieważ tkanka tłuszczowa wykazuje naturalnie wysoką intensywność sygnału w stosunku do struktur otaczających. Wybrane parametry akwizycji maksymalizują kontrast T1, utrzymując jednocześnie odpowiedni stosunek sygnału do szumu i rozdzielczość przestrzenną. Typowy czas akwizycji wynosi około 3–4 min.
   2. Podczas adaptacji protokołu do równoważnej platformy 1,5 lub 3,0 T, zachowaj wymagany zakres anatomiczny, typ sekwencji, orientację warstw, grubość warstw, odstępy między warstwami oraz płaszczyzny diagnostyczne. Zoptymalizuj pole widzenia, macierz, liczbę wzbudzeń, długość pociągu ech i współczynnik obrazowania równoległego w razie potrzeby, pod warunkiem, że tłuszcz nadtwardówkowy, brzegi worka oponowego i anatomia poziomu przesunięcia pozostaną wyraźnie widoczne do interpretacji diagnostycznej.
4. Nabyj sagitalne obrazy T2-zależne turbo spin-echo, stosując TR = 3 500-4 000 ms, TE = 100-120 ms, pole widzenia = 280 mm × 280 mm, macierz = 512 × 256, grubość warstwy = 4 mm, odstęp między warstwami = 0,4 mm oraz NEX = 2.
   1. Dostosuj rozmieszczenie warstw do akwizycji sagitalnej T1-zależnej, aby umożliwić bezpośrednie porównanie obrazów.
      ​UWAGA: Obrazowanie T2-zależne dostarcza uzupełniających informacji anatomicznych poprzez poprawę wizualizacji płynu mózgowo-rdzeniowego (PMR), worka oponowego i morfologii krążków międzykręgowych. Sekwencja ta jest również wykorzystywana do oceny stopnia degeneracji krążków międzykręgowych według klasyfikacji Pfirrmanna. Typowy czas akwizycji wynosi około 4–5 min.
5. Nabyj aksjalne obrazy T2-zależne fast spin-echo na każdym poziomie krążka międzykręgowego odcinka lędźwiowego od L1/L2 do L5/S1, stosując TR = 4 000 ms, TE = 112 ms, pole widzenia = 180 mm × 180 mm, macierz = 320 × 256, grubość warstwy = 4 mm, odstęp między warstwami = 0,4 mm oraz NEX = 3.
   1. Wyznacz warstwy aksjalne równolegle do każdej przestrzeni międzykręgowej. Nabyj 6-8 warstw w obrębie segmentu przesuniętego, aby zapewnić pełny zakres anatomiczny.
      ​UWAGA: Aksjalne obrazowanie T2-zależne pozwala na ocenę rozkładu tłuszczu nadtwardówkowego i morfologii worka oponowego. Połączenie sekwencji sagitalnych T1-zależnych, sagitalnych T2-zależnych i aksjalnych T2-zależnych dostarcza komplementarnych informacji anatomicznych do kompleksowej oceny. Całkowity czas akwizycji dla wszystkich poziomów lędźwiowych wynosi około 15–20 min. Całe badanie, wliczając ułożenie pacjenta i obrazowanie lokalizacyjne, trwa około 30–35 min.
6. Oceń jakość obrazu przed zakończeniem badania. Potwierdź brak artefaktów ruchowych stopnia II lub wyższego, pełną wizualizację przejścia lędźwiowo-krzyżowego, odpowiedni stosunek sygnału do szumu oraz brak istotnych artefaktów zawijania (wrap-around).
   1. Powtórz każdą sekwencję, która nie spełnia kryteriów jakości obrazu. Eksportuj wszystkie obrazy w formacie DICOM (Digital Imaging and Communications in Medicine) do systemu PACS z zachowaniem metadanych.
   2. Zdefiniuj akceptowalną jakość obrazu jako pełny zakres obrazowania od T12/L1 do kości krzyżowej, ostrość obrazowania tylnych krawędzi korowych ciał kręgów i brzegów worka oponowego, brak artefaktu zawijania wnikającego do kanału kręgowego oraz wystarczający stosunek sygnału do szumu, aby odróżnić tłuszcz nadtwardówkowy od PMR i mięśni przykręgosłupowych. Powtórz każdą sekwencję, która nie spełnia tych kryteriów.
   3. Użyj tego samego czterostopniowego systemu oceny artefaktów ruchowych opisanego w kroku 1.3.3. Powtórz sekwencje z artefaktami ruchowymi stopnia II lub III przed zwolnieniem pacjenta lub odrzuć badanie, jeśli ponowne obrazowanie nie jest możliwe.
      UWAGA: Zweryfikuj jakość obrazu przed wypuszczeniem pacjenta ze skanera.

3. Interpretacja obrazu MRI w przypadku lipomatozy nadtwardówkowej

UWAGA: Lipomatoza nadtwardówkowa jest definiowana jako patologiczny rozrost nieotoczonej torebką dojrzałej tkanki tłuszczowej w przestrzeni nadtwardówkowej kręgosłupa. W przypadku DLS stan ten przyczynia się do zwężenia kanału kręgowego i może wymagać połączonej dekompresji podczas operacji stabilizacji. Dokładna diagnoza wymaga systematycznej oceny w wielu sekwencjach i płaszczyznach obrazowania.

1. Przeprowadź interpretację obrazów na dedykowanej stacji roboczej PACS wyposażonej w diagnostyczny monitor medyczny (minimum 3 megapiksele; DICOM Part 14 Grayscale Standard Display Function). Zachowaj kontrolowane oświetlenie otoczenia (≤50 lux).
   1. Wczytaj kompletne badanie MRI. Wyświetl obok siebie obrazy T1-zależne w płaszczyźnie strzałkowej oraz obrazy T2-zależne w płaszczyźnie strzałkowej w górnym oknie podglądu oraz obrazy T2-zależne w płaszczyźnie osiowej w dolnym oknie podglądu, z włączoną zsynchronizowaną referencją anatomiczną.
   2. Użyj instytucyjnego przeglądarki PACS z narzędziami pomiarowymi skalibrowanymi zgodnie ze standardem DICOM i włączoną zsynchronizowaną referencją strzałkową-osiową. Stosuj spójne ustawienia okna (window) i poziomu (level) podczas porównywania obrazów obok siebie w każdej sesji odczytu. Zapisz producenta systemu PACS oraz wersję oprogramowania w Tabeli materiałów.
2. Zidentyfikuj poziom przesunięcia kręgu na obrazie T1-zależnym w płaszczyźnie środkowej strzałkowej, obserwując przednie przemieszczenie kręgu. Narysuj linię referencyjną wzdłuż tylnej kory przyległego poniżej ciała kręgu.
   1. Zmierz prostopadłą odległość od linii referencyjnej do tylno-dolnego rogu przesuniętego kręgu. Potwierdź, że zmierzone przemieszczenie wynosi ≥3 mm lub ≥10% szerokości ciała kręgu.
      ​UWAGA: Dokładna lokalizacja poziomu przesunięcia jest kluczowa, ponieważ lipomatozę nadtwardówkową ocenia się konkretnie na poziomie przesunięcia kręgu.
   2. Wykonaj pomiary przesunięcia kręgu za pomocą wbudowanego, skalibrowanego elektronicznego suwmiarki na stacji roboczej PACS. Przed pomiarem zweryfikuj kalibrację pikseli w metadanych DICOM i zapisz wszystkie odległości z dokładnością do 0,1 mm.
3. Zbadaj tylną przestrzeń nadtwardówkową na poziomie przesunięcia na obrazie T1-zależnym w płaszczyźnie środkowej strzałkowej. Zidentyfikuj pasmowaty lub półksiężycowaty sygnał hiperintensywny znajdujący się pomiędzy workiem oponowym z przodu a więzadłem żółtym lub wyrostkiem kolczystym z tyłu.
   1. Porównaj bezpośrednio intensywność sygnału zmiany z tkanką tłuszczową podskórną na tym samym obrazie. Jeśli dostępne są sekwencje z tłumieniem tłuszczu, potwierdź tłumienie sygnału, aby zweryfikować tłuszczowy charakter zmiany.
      ​UWAGA: Odróżnij lipomatozę nadtwardówkową od innych hiperintensywnych w T1 zmian nadtwardówkowych, w tym od podostrego krwiaka nadtwardówkowego i zbiorników płynu białkowego. Potwierdzaj lipomatozę nadtwardówkową tylko wtedy, gdy zmiana wykazuje charakterystykę sygnału identyczną z tłuszczem podskórnym i dopasowuje się do przestrzeni nadtwardówkowej.
   2. Stosuj zdefiniowane w protokole kryteria strzałkowe jako podstawową metodę przesiewową poziomu przesunięcia i weryfikuj wyniki dodatnie lub niejednoznaczne na obrazach T2-zależnych w płaszczyźnie osiowej. Podczas wdrożenia prospektywnego dokumentuj rozkład tłuszczu nadtwardówkowego, deformację worka oponowego oraz zgodność z uznanymi metodami oceny opartymi na MRI, w tym podejściami opartymi na stosunku tłuszczu nadtwardówkowego do worka oponowego w płaszczyźnie osiowej oraz lokoregionalnym systemie stopniowania Manjila. W kohorcie retrospektywnej wykorzystywanej do opracowania modelu nie stosowano tych uznanych systemów stopniowania jako niezależnych standardów referencyjnych, a pierwotny wynik binarny nie był retrospektywnie reklasyfikowany.
4. Zmierz zasięg kraniokaudalny złogu tłuszczu nadtwardówkowego za pomocą elektronicznej suwmiarki na stacji roboczej PACS. Umieść suwmiarkę na górnej i dolnej krawędzi ciągłego złogu tłuszczu nadtwardówkowego i zapisz maksymalny pomiar.
   1. Potwierdź, że zasięg kraniokaudalny wynosi ≥5 mm w tej samej płaszczyźnie strzałkowej. Wykorzystuj widoczność na co najmniej dwóch sąsiednich przekrojach strzałkowych wyłącznie do potwierdzenia ciągłości w płaszczyźnie mediolateralnej i redukcji artefaktów częściowej objętości; nie używaj sąsiednich przekrojów strzałkowych do szacowania zasięgu kraniokaudalnego.
   2. Zweryfikuj zmianę na obrazach T2-zależnych w płaszczyźnie strzałkowej, aby ocenić charakterystykę sygnału. Przeanalizuj odpowiadające im obrazy T2-zależne w płaszczyźnie osiowej, aby ocenić obwodowy rozkład tłuszczu nadtwardówkowego oraz redukcję przednio-tylną worka oponowego.
      ​UWAGA: Przyjmij zasięg kraniokaudalny ≥5 mm jako operacyjny próg przesiewowy dla poziomu przesunięcia w tym protokole. Zapisz średnicę przednio-tylną oraz pole przekroju poprzecznego worka oponowego na obrazach osiowych jako ciągłe pomiary potwierdzające. Nie stosuj stałych wartości odcięcia w płaszczyźnie osiowej jako niezależnych standardów referencyjnych w niniejszej kohorcie retrospektywnej.
   3. Jeśli sygnał tłuszczowy jest widoczny tylko na jednym przekroju strzałkowym, zaklasyfikuj znalezisko jako niejednoznaczne, chyba że obrazy osiowe wykazują odpowiadającą akumulację tłuszczu nadtwardówkowego, a zasięg kraniokaudalny w tej samej płaszczyźnie nadal wynosi ≥5 mm. Kryterium to odróżnia główny pomiar kraniokaudalny od wspomagającej oceny ciągłości mediolateralnej.
   4. Zmierz średnicę przednio-tylną worka oponowego na obrazie T2-zależnym w płaszczyźnie osiowej wykazującym maksymalną kompresję na poziomie przesunięcia, mierząc odległość od przedniej do tylnej krawędzi worka oponowego. Zmierz pole przekroju poprzecznego, obrysowując wewnętrzną granicę oponytwardej za pomocą narzędzia region-of-interest w systemie PACS lub narzędzia do wyboru wielokąta w programie ImageJ, i zapisz powierzchnię w mm².
   5. Jakościowo porównaj zasięg strzałkowy i stopień nasilenia w płaszczyźnie osiowej z ramami stopniowania Manjila. W tym protokole skoncentrowanym na DLS stosuj system stopniowania Manjila jako lokoregionalny punkt odniesienia do dokumentacji, a nie jako główną definicję wyniku, ponieważ punktem końcowym badania jest obecność lipomatozy nadtwardówkowej na poziomie przesunięcia w modelu predykcyjnym.
5. Zastosuj cztery pierwotne kryteria poziomu przesunięcia użyte do klasyfikacji kohorty: (1) pasmowaty lub półksiężycowaty sygnał hiperintensywny w T1 w tylnej przestrzeni nadtwardówkowej na poziomie przesunięcia; (2) intensywność sygnału identyczna z tłuszczem podskórnym; (3) zasięg kraniokaudalny ≥5 mm w tej samej płaszczyźnie strzałkowej; oraz (4) powtarzalność na sąsiednich przekrojach strzałkowych. W przypadku prospektywnego wdrożenia klinicznego nie polegaj wyłącznie na kryteriach strzałkowych. Potwierdź odpowiadającą akumulację tłuszczu nadtwardówkowego i deformację konturu worka oponowego na obrazach osiowych oraz udokumentuj, czy rozkład tłuszczu jest grzbietowy, brzuszny czy obwodowy.
   1. Zapisz ilościowe wyniki w płaszczyźnie osiowej, w tym średnicę przednio-tylną i pole przekroju poprzecznego worka oponowego, wraz z jakościową zgodnością z metodami oceny opartymi na stosunku tłuszczu nadtwardówkowego do worka oponowego i lokoregionalnym systemem stopniowania Manjila.
      UWAGA: Pierwotne badanie nie zostało zaprojektowane jako badanie dokładności diagnostycznej z użyciem uznanego standardu referencyjnego MRI, a oryginalne obrazy DICOM nie były ponownie analizowane podczas rewizji manuskryptu. W związku z tym nie przeliczano czułości, swoistości ani zgodności między metodami. Te dodatkowe deskryptory osiowe mają na celu wsparcie prospektywnego wdrożenia klinicznego i nie zmieniają pierwotnego wyniku badania ani modelu predykcyjnego.
6. Zleć dwóm doświadczonym radiologom kręgosłupa (≥10 lat doświadczenia) niezależną interpretację wszystkich badań. Wypełnij ustandaryzowany, ustrukturyzowany formularz raportowania, dokumentujący obecność lub brak każdego kryterium diagnostycznego, poziom przesunięcia kręgu, maksymalny zasięg kraniokaudalny (mm) oraz diagnozę ogólną.
   1. Rozstrzygnij rozbieżności poprzez konsensus z trzecim doświadczonym radiologiem kręgosłupa (>15 lat doświadczenia), który nie zna wyników wstępnych ocen. Oblicz zgodność między obserwatorami za pomocą statystyki kappa Cohena.
   2. Użyj ustrukturyzowanego formularza raportowania, aby udokumentować poziom przesunięcia, rozkład tłuszczu nadtwardówkowego (grzbietowy, brzuszny lub obwodowy), maksymalny zasięg kraniokaudalny, pomiar strzałkowy w tej samej płaszczyźnie, widoczność na sąsiednich przekrojach, potwierdzenie osiowe, średnicę worka oponowego, pole przekroju poprzecznego worka oponowego, ostateczną diagnozę oraz pewność interpretatora. Przed formalną interpretacją obrazów skalibruj wszystkich obserwatorów, wykorzystując 20 reprezentatywnych przypadków szkoleniowych, obejmujących badania negatywne, graniczne i dodatnie.
      UWAGA: Dokończ wszystkie interpretacje konsensualne przed przystąpieniem do analizy statystycznej.

4. Ocena degeneracji krążka międzykręgowego

1. Przejrzyj obrazy MRI w płaszczyźnie strzałkowej z obrazowaniem T2-zależnym w miejscu przesuniętego segmentu. Przewiń wszystkie przekroje strzałkowe i wybierz obraz, który zapewnia najwyraźniejszą wizualizację jądra miażdżystego (nucleus pulposus).
   1. Określ stopień degeneracji krążka międzykręgowego, korzystając z systemu stopniowania Pfirrmanna¹⁸. Przypisz stopień I do krążków z jednorodnym, hiperintensywnym sygnałem identycznym z płynem mózgowo-rdzeniowym, prawidłową wysokością krążka oraz wyraźnym rozróżnieniem między jądrem a pierścieniem włóknistym.
   2. Przypisz stopień II do krążków z niejednorodnym sygnałem hiperintensywnym zawierającym poziome pasma niskiego sygnału, zachowaną wysokością krążka i wyraźnym rozróżnieniem między jądrem a pierścieniem włóknistym. Przypisz stopień III do krążków z sygnałem w odcieniu szarym, niewyraźną granicą między jądrem a pierścieniem włóknistym oraz prawidłową lub lekko zmniejszoną wysokością krążka.
   3. Przypisz stopień IV do krążków z niejednorodnym, hipointensywnym sygnałem w kolorze ciemnoszarym, całkowitym brakiem rozróżnienia między jądrem a pierścieniem włóknistym oraz umiarkowanie lub znacznie zmniejszoną wysokością krążka. Przypisz stopień V do krążków z niejednorodnym sygnałem hipointensywnym, zapadniętą przestrzenią krążka, całkowitym brakiem rozróżnienia między jądrem a pierścieniem włóknistym oraz wyraźnym zwężeniem przestrzeni krążka.
2. Zleć dwóm certyfikowanym radiologom z co najmniej 5-letnim doświadczeniem w obrazowaniu układu mięśniowo-szkieletowego niezależne przypisanie stopni Pfirrmanna. Zapisz przypisany stopień dla każdego ocenianego krążka.
   1. Rozstrzygnij rozbieżności w stopniowaniu poprzez arbitraż trzeciego certyfikowanego radiologa. Oblicz zgodność między obserwatorami, korzystając z ważonej statystyki kappa Cohena.
   2. Zastosuj metodę ślepej próby dla radiologów w odniesieniu do objawów klinicznych, historii leczenia, stanu lipomatozy znadtwardówkowej, zmiennych modelu predykcyjnego oraz ocen pozostałych badaczy podczas niezależnego stopniowania według Pfirrmanna.
   3. Przeprowadź stopniowanie według Pfirrmanna, korzystając z tych samych monitorów o jakości diagnostycznej i ustandaryzowanych warunków oglądania opisanych w kroku 3.1.
      UWAGA: Całkowite uzgodnienie stopniowania musi zostać zakończone przed przejściem do oceny nacieku tłuszczowego mięśni przykręgosłupowych.

5. Ocena nacieku tłuszczowego mięśni przykręgosłupowych

1. Przejrzyj osiowe obrazy T2-zależne fast spin-echo na poziomie środka przesuniętego krążka międzykręgowego. Zidentyfikuj obustronnie mięśnie wielodzielne jako głębokie mięśnie przykręgosłupowe położone bezpośrednio bocznie względem wyrostka kolczystego i łuku kręgu.
   1. Przejrzyj obrazy strzałkowe, aby określić stronę przeważającego przesunięcia kręgu. Wybierz mięsień wielodzielny po odpowiadającej stronie do oceny.
   2. Zastosuj podejście do stopniowania według Goutalliera dostosowane do kręgosłupa, wywodzące się z oryginalnej klasyfikacji tłuszczowej degeneracji mięśni i wcześniejszych zastosowań w odniesieniu do mięśnia wielodzielnego odcinka lędźwiowego^19,20. Przypisz Stopień 0 mięśniom wykazującym jednorodną niską intensywność sygnału, całkowity brak smug tłuszczowych i dobrze zachowaną masę mięśniową. Przypisz Stopień 1 mięśniom zawierającym małe, liniowe, wysokosygnałowe smug tłuszczowe zajmujące <5% pola przekroju poprzecznego.
   3. Przypisz Stopień 2 mięśniom wykazującym wyraźnie widoczną infiltrację tłuszczową obejmującą 5%–50% pola przekroju poprzecznego. Przypisz Stopień 3 mięśniom wykazującym infiltrację tłuszczową ≥50% z zachowaniem resztkowej tkanki mięśniowej.
   4. Przypisz Stopień 4 mięśniom, w których cały przekrój poprzeczny został zastąpiony przez tłuszcz, bez widocznej tkanki mięśniowej. Na potrzeby analizy statystycznej dokonaj dychotomizacji wyników na Stopnie 0-2 oraz Stopnie 3-4.
      ​UWAGA: Mięsień wielodzielny jest największym i najbardziej przyśrodkowym mięśniem przykręgosłupowym i pełni rolę głównego dynamicznego stabilizatora segmentu lędźwiowego. Silna infiltracja tłuszczowa (Stopnie 3-4) wskazuje na znaczne zastąpienie tłuszczowe i jest traktowana w modelowaniu jako kategoria klinicznie ciężka, ponieważ odzwierciedla próg, przy którym tłuszcz zajmuje co najmniej około połowy pola przekroju poprzecznego mięśnia.
   5. Określ, że oryginalny system stopniowania Goutalliera został dostosowany do oceny mięśnia wielodzielnego lędźwiowego w MRI osiowym. Dokonaj dychotomizacji Stopni 0–2 i Stopni 3–4, aby odróżnić brak lub umiarkowaną infiltrację tłuszczową od silnego zastąpienia tłuszczowego, zmniejszając tym samym złożoność modelu przy zachowaniu odpowiedniego stosunku zdarzeń do zmiennych w modelu predykcyjnym. Progi procentowe użyte w tym protokole traktuj jako operacyjne kryteria oparte na MRI, a nie jako dosłowne odwzorowanie oryginalnej klasyfikacji Goutalliera dla tomografii komputerowej barku.
   6. Określ stronę przeważającego przesunięcia kręgu, analizując obrazy osiowe i strzałkowe pod kątem asymetrycznej translacji kręgu, przemieszczenia rotacyjnego lub większego zwężenia recesowania bocznego. Jeśli żadna strona nie dominuje, oceń stronę wykazującą większą infiltrację tłuszczową mięśnia wielodzielnego. Jeśli obie strony są symetryczne, oceń prawy mięsień wielodzielny i udokumentuj tę decyzję.
2. Wyeksportuj osiowy obraz DICOM T2-zależny fast spin-echo z systemu PACS. Otwórz obraz w oprogramowaniu ImageJ (wersja 1.53t).
   1. Wybierz narzędzie zaznaczania wielokątnego i obrysuj granicę powięziową docelowego mięśnia wielodzielnego. Zidentyfikuj granicę powięziową jako cienką ciemną linię otaczającą mięsień.
   2. Wybierz Image > Adjust > Threshold. Ustaw dolny próg na 120 jednostek arbitralnych, a górny próg na maksymalną wartość piksela.
   3. Wybierz kolor czerwony jako kolor wyświetlania, aby zweryfikować identyfikację obszarów tłuszczowych. Kliknij Apply, aby wygenerować maskę binarną.
   4. Wartość progu 120 jednostek arbitralnych traktuj jako punkt wyjścia specyficzny dla skanera i sekwencji, a nie jako próg uniwersalny. Ponownie skalibruj próg, używając reprezentatywnych obrazów lokalnych, ilekroć platforma skanera, cewka radiowej częstotliwości, natężenie pola magnetycznego lub parametry akwizycji obrazu różnią się od tych stosowanych w protokole bazowym.
   5. Nie interpretuj wartości progu jako zwalidowanej histologicznie dla wszystkich systemów obrazowania. Używaj jej jako powtarzalnego pomocnika w przetwarzaniu obrazu i wizualnie weryfikuj każdą maskę binarną względem oryginalnego obrazu T2-zależnego w płaszczyźnie osiowej przed wykonaniem pomiarów ilościowych.
3. Otwórz Analyze > Measure w ImageJ. Przed wykonaniem pomiaru włącz opcje Area, Area Fraction, Limit to Threshold oraz Display Label.
   1. Zapisz wartość %Area jako ilościowy procent tłuszczu wewnątrzmięśniowego.
   2. Powtórz pomiar po upływie co najmniej 24 h przez tego samego operatora. Oblicz współczynnik korelacji wklasowej (ICC), aby ocenić niezawodność wewnątrzoperatorową.
      ​UWAGA: ICC >0,90 wskazuje na doskonałą niezawodność wewnątrzoperatorową.
   3. Podczas powtórnych pomiarów zapewnij zaślepienie operatora względem pomiaru początkowego poprzez przechowywanie oryginalnej wartości w zablokowanym polu REDCap i ponowne otwarcie obrazu DICOM jako nowej sesji pomiarowej po upływie co najmniej 24 h.
4. Zatrudnij dwóch certyfikowanych radiologów do niezależnego przypisywania stopni według Goutalliera i wykonywania ilościowych pomiarów w programie ImageJ. Przed przystąpieniem do ocen w badaniu zapewnij zestandaryzowane szkolenie z wykorzystaniem referencyjnego zestawu obrazów.
   1. Rozwiązuj rozbieżności poprzez konsensus z trzecim, starszym radiologiem muskuloszkieletowym. Oblicz kapę Cohena dla klasyfikacji Goutalliera oraz współczynnik korelacji wklasowej dla ilościowego procentu tłuszczu.
   2. Przed oceną w badaniu wykorzystaj zestaw treningowy składający się z 20 przypadków, w tym po pięć przypadków reprezentujących minimalną, lekką, umiarkowaną i silną infiltrację tłuszczową. Przeprowadź sesję kalibracyjną w celu przeanalizowania obrazów referencyjnych, rozmieszczenia obszarów zainteresowania, dostrojenia progu oraz procedur rozwiązywania rozbieżności przed rozpoczęciem niezależnego stopniowania obrazów.
      UWAGA: Jeśli kappa Cohena między obserwatorami spadnie poniżej 0,70 lub współczynnik korelacji wklasowej spadnie poniżej 0,85, przeprowadź sesję rekalibracyjną przed wznowieniem analizy.

6. Opracowanie i zastosowanie modelu predykcji ryzyka

1. Zgromadzić wszystkie zmienne badania w ustrukturyzowanym zbiorze danych przy użyciu systemu REDCap. Zdefiniować zmienną wynikową jako binarną: 1 = obecność lipomatozy zewnątrzoponowej (spełnienie wszystkich czterech kryteriów diagnostycznych z kroku 3.5) oraz 0 = brak lipomatozy zewnątrzoponowej.
   1. Zakoduj zmienne predyktora w następujący sposób: wiek (lata); płeć (0 = mężczyzna, 1 = kobieta); BMI (kg/m²).²); przesunięty segment (0 = L3/L4, 1 = L5) oraz stopień według Goutalliera (0 = stopnie 0–2; 1 = stopnie 3–4, reprezentujące ciężką infiltrację tłuszczową). W raporcie należy przedstawić efekt wieku na każde 10 lat wzrostu oraz efekt BMI na każde 5 kg/m²² wzrost. W równaniu predykcyjnym stosowanym przy łóżku pacjenta należy użyć współczynników jednostek surowych, które są algebraicznie równoważne raportowanym współczynnikom opartym na przyrostach.
   2. Przeprowadź czyszczenie danych poprzez identyfikację brakujących wartości, weryfikację wartości odstających oraz potwierdzenie, że wszystkie zmienne mieszczą się w oczekiwanych zakresach.
   3. Wyklucz pacjentów z brakującymi wymaganymi zmiennymi wynikowymi lub predyktorami po zweryfikowaniu dokumentacji źródłowej. Sprawdź zmienne ciągłe, korzystając z analizy zakresów i wykresów rozrzutu, popraw zweryfikowane błędy wprowadzania danych na podstawie dokumentacji źródłowej, zachowaj biologicznie prawdopodobne wartości odstające i nie stosuj imputacji dla wymaganych zmiennych modelu.
2. Otwórz program SPSS w wersji 27.0 i ustaw ziarno generatora liczb losowych na 42, aby zapewnić powtarzalność wyników. Przeprowadź przesiewową analizę jednowymiarową przy użyciu testu t Studenta Proszę podać tekst do przetłumaczenia.-test dla zmiennych ciągłych oraz test chi-kwadrat dla zmiennych kategorycznych.
   1. Wybierz zmienne z P < 0,05 dla analizy wielowymiarowej. Dopasowano wielowymiarowy model regresji logistycznej binarnej z zastosowaniem wstecznej eliminacji krokowej w oparciu o minimalizację kryterium informacyjnego Akaike (AIC).
   2. Oceń dopasowanie modelu za pomocą testu Hosmera-Lemeshow. Oblicz współczynnik inflacji wariancji (VIF) dla wszystkich zachowanych predyktorów i potwierdź, że wszystkie wartości VIF są <5.
   3. Oblicz wartość wykładniczą współczynników regresji, aby uzyskać skorygowane ilorazy szans (OR) z 95-procentowymi przedziałami ufności (CI). Wygeneruj wykres lasu (forest plot) przy użyciu pakietu forestplot w programie R, aby zwizualizować wielkości efektów (Rysunek 1). Należy zgłosić efekt wieku na każde 10 lat wzrostu oraz efekt BMI na każde 5 kg/m²² wzrost
   4. W celu opracowania modelu eksploracyjnego dla tego zbioru danych z jednego ośrodka należy zastosować jednoczynnikową analizę przesiewową, a następnie procedurę wstecznej eliminacji krokowej. Końcowy model należy interpretować jako narzędzie wspomagania decyzji klinicznych służące do generowania hipotez, ponieważ selekcja krokowa może prowadzić do niestabilnych oszacowań współczynników oraz zbyt optymistycznych wyników oceny wydajności modelu. Przed wdrożeniem klinicznym wydajność modelu musi zostać potwierdzona poprzez walidację zewnętrzną.
   5. Wyeksportuj oczyszczony zbiór danych SPSS jako plik wartości rozdzielonych przecinkami (.csv), usuwając etykiety zmiennych i zachowując kodowanie. Wygeneruj wykres leśny (forest plot) w programie R w wersji 4.3.1, korzystając z pakietu forestplot. Przed sporządzeniem wykresu potwierdź, że kodowanie zmiennych w programie R jest zgodne ze słownikiem kodowania SPSS.
3. Oceń dyskryminację modelu poprzez obliczenie powierzchni pod krzywą charakterystyki operacyjnej odbiornika (AUC) dla kohorty modelowej. Oszacuj błąd standardowy, stosując test DeLonga.
   1. Oceń kalibrację modelu za pomocą wykresów kalibracyjnych, które porównują prawdopodobieństwa przewidywane i obserwowane w grupach przedziałów przewidywanego ryzyka zarówno w kohorcie modelowania, jak i walidacji (Rysunek 2). Krzywą kalibracyjną kohorty walidacyjnej należy interpretować z ostrożnością, ponieważ w skrajnych obszarach przewidywanego ryzyka dostępnych jest mniej obserwacji.
   2. Przeprowadź analizę krzywej decyzyjnej (decision curve analysis), wykorzystując pakiet rmda (wersja 1.6) w programie R w wersji 4.3.1. Porównaj model predykcyjny ze strategiami „lecz wszystkich” (treat-all) i „nie lecz nikogo” (treat-none) w zakresie prawdopodobieństw progowych (Rycina 3).
   3. Zweryfikuj ostateczny model z wykorzystaniem wydzielonej wewnętrznej kohorty walidacyjnej (n = 106). Przedstaw w formie opisowej AUC walidacji, wykres kalibracji oraz krzywą decyzyjną. Zgodność między kohortą modelowania a kohortą walidacyjną należy interpretować jako dowód spójności wewnętrznej, a nie jako dowód zewnętrznej trafności.
   4. Przeprowadź analizę charakterystyki operacyjnej odbiornika (ROC), wykreślenie kalibracji oraz analizę krzywej decyzyjnej w programie R w wersji 4.3.1. Zarchiwizuj skrypty R, wersje pakietów, zestaw danych wejściowych, słownik kodowania oraz otrzymane wykresy wraz z dokumentacją analizy statystycznej.
   5. Nie należy stosować stałych progów, takich jak różnica AUC ≤0,05 lub nachylenie kalibracji w zakresie 0,8–1,2, jako formalnych kryteriów walidacji. Wartości te należy raportować opisowo jako wskaźniki wewnętrznej stabilności modelu oraz stwierdzić, że przed szerokim zastosowaniem klinicznym wymagana jest wieloośrodkowa walidacja zewnętrzna.
      UWAGA: Przed zastosowaniem modelu do zewnętrznych zbiorów danych klinicznych należy przeprowadzić pełną walidację wewnętrzną.
4. Zebierz pięć zmiennych predyktora dla nowego pacjenta z DLS: wiek (lata), płeć, BMI (kg/m²²), przesunięty segment (L3/L4 lub L5) oraz stopień Goutalliera mięśnia wielodzielnego na poziomie przesunięcia.
   1. Oblicz przewidywane prawdopodobieństwo, korzystając z równania regresji logistycznej wyrażonego w surowych jednostkach klinicznych:
      
      Tutaj, z = -5,521 + (0,0412 × wiek) + (0,856 × płeć) + (0,1256 × BMI) + (1,326 × przesunięty segment) + (1,158 × stopień według Goutalliera). Wprowadzić wiek w latach; zakodować płeć jako 0 = mężczyzna i 1 = kobieta; wprowadzić BMI w kg/m²²; zakodowano przesunięty segment jako 0 = L3/L4 oraz 1 = L5; a stopień wg Goutalliera zakodowano jako 0 = stopnie 0–2 oraz 1 = stopnie 3–4. Współczynniki dla jednostek surowych w przypadku wieku i BMI są algebraicznie równoważne współczynnikom wyrażonym na 10 lat i na 5 kg/m²² odpowiednio zwiększają.
   2. Zinterpretuj przewidywane prawdopodobieństwo opisowo w następujący sposób: <wartość 20% wskazuje na niskie szacowane prawdopodobieństwo, 20%–50% wskazuje na pośrednie szacowane prawdopodobieństwo, a >Wartość 50% wskazuje na wysokie szacowane prawdopodobieństwo wystąpienia lipomatozy nadtwardówkowej. Kategorii tych należy używać do komunikacji ryzyka i ustalania priorytetów w procesie diagnostycznym, a nie do wyboru metody leczenia.
      UWAGA: Modelu predykcyjnego należy używać wyłącznie jako narzędzia wspomagania decyzji klinicznych. Model nie może zastąpić oceny klinicznej ani interpretacji radiologicznej. Nie należy stosować modelu w populacjach spoza kohorty pochodzenia, w tym u pacjentów z spondylolistezą istmową, pooperacyjnymi stanami kręgosłupa lub steroidową lipomatozą zewnątrzoponową, bez przeprowadzenia dodatkowej walidacji.
   3. Wartości wieku i BMI wprowadza się bezpośrednio do równania predykcyjnego. Algebraiczne przeskalowanie współczynników wieku i BMI zmienia jedynie prezentację równania i nie wpływa na przewidywane prawdopodobieństwa, ilorazy szans, krzywe charakterystyki operacyjnej odbiorcy (ROC), wykresy kalibracyjne ani wyniki analizy krzywej decyzyjnej.
   4. Ilorazy szans należy raportować konsekwentnie w następujący sposób: wiek na każde 10 lat wzrostu; BMI na każde 5 kg/m²² wzrost; płeć, żeńska versus męska; przesunięty segment, L5 versus L3/L4; oraz stopień wg Goutalliera, stopnie 3–4 versus stopnie 0–2.
   5. Należy wykorzystać kategorie niskiego, średniego i wysokiego prawdopodobieństwa jako pragmatyczne grupy komunikacji klinicznej, wyłonione na podstawie przewidnianego przez model rozkładu prawdopodobieństwa oraz analizy krzywych decyzyjnych. Kategorii tych nie należy interpretować jako zwalidowanych progów terapeutycznych, ponieważ do ustalenia optymalnych progów decyzyjnych wymagana jest prospektywna walidacja zewnętrzna.

Rycina 1. Wykres leśny (forest plot) niezależnych predyktorów lipomatozy zewnątrzoponowej u pacjentów z degeneracyjnym spondylolistezą lędźwiową. Wykres leśny przedstawiający skorygowane ilorazy szans (OR) oraz 95% przedziały ufności (CI) dla pięciu niezależnych predyktorów zidentyfikowanych za pomocą wieloczynnikowej binarnej regresji logistycznej. Czerwone kwadraty oznaczają skorygowane OR, poziome niebieskie linie reprezentują odpowiadające im 95% CI, a pionowa przerywana linia wskazuje wartość zerową (OR = 1). Efekt wieku podano w przeliczeniu na wzrost o 10 lat, a efekt wskaźnika masy ciała (BMI) podano w przeliczeniu na wzrost o 5 kg/m². Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rysunek 2. Krzywe kalibracyjne modelu predykcyjnego regresji logistycznej. (A) Krzywa kalibracyjna dla kohorty modelowania (n = 248). (B) Krzywa kalibracyjna dla kohorty walidacyjnej (n = 106). Oś x reprezentuje przewidywane prawdopodobieństwo lipomatozy zewnątrzoponowej, a oś y reprezentuje zaobserwowaną częstość. Ciągła czarna linia diagonalna oznacza idealną kalibrację, kolorowe linie z symbolami reprezentują pozorny wynik modelu, a obszary zacieniowane oznaczają 95% przedziały ufności (C.I.). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 3. Krzywa charakterystyki operacyjnej odbiornika i analiza krzywej decyzyjnej modelu predykcyjnego. (A) Krzywe charakterystyki operacyjnej odbiornika (ROC) dla kohorty modelowej (pole pod krzywą [AUC] = 0.834) i kohorty walidacyjnej (AUC = 0.815). Szara linia przekątna reprezentuje linię odniesienia dla klasyfikatora niedyskryminującego. (B) Analiza krzywej decyzyjnej pokazująca korzyść netto modelu predykcyjnego w zależności od prawdopodobieństw progowych. Czerwona linia ciągła reprezentuje kohortę modelową, niebieska linia przerywana reprezentuje kohortę walidacyjną, czarna linia kropka-kreska reprezentuje strategię leczenia wszystkich, a szara linia kropkowana reprezentuje strategię nieleczenia nikogo. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Wyniki

Protokół zastosowano u 354 kolejnych pacjentów z DLS. Losowy przydział pozwolił na utworzenie kohorty modelowania (n = 248) oraz kohorty walidacyjnej (n = 106). Porównanie demograficznych, klinicznych, radiograficznych oraz charakterystyk MRI w punkcie wyjściowym nie wykazało statystycznie istotnych różnic między kohortami, co wskazuje na porównywalny rozkład pacjentów w kontekście późniejszego opracowania modelu i walidacji wewnętrznej (Tabela 1).

W grupie modelowej u 87 pacjentów stwierdzono lipomatozę nadtwardówkową, a u 161 pacjentów nie stwierdzono lipomatozy nadtwardówkowej, zgodnie ze zstandaryzowanym protokołem MRI. Analiza jednowymiarowa wykazała różnice między grupami pod względem wieku (P = 0,014), BMI (P < 0,001), segmentu przesuniętego (P = 0,010) oraz stopnia według Goutalliera (P = 0,024). Płeć żeńska znajdowała się na nominalnym progu istotności (P = 0,050). Nie zaobserwowano statystycznie istotnych różnic w przypadku stopnia według Pfirrmanna ani pozostałych parametrów kręgosłupowo-miednicznych i stawów międzywyrostkowych (Tabela 2).

Zmienna	Grupa z lipomatozą (n = 87)	Grupa bez lipomatozy (n = 161)	t/χ2	P wartość
Wiek (lata)	70.28 ± 10.84	66.57 ± 11.50	2.47	0.014
Płeć, n (%)			3.844	0.05
Mężczyzna	26 (29.89)	70 (43.48)
Kobieta	61 (70.11)	91 (56.52)
Wskaźnik masy ciała (kg/m²)	26.84 ± 3.68	23.82 ± 3.12	6.815	<0.001
Segment przesunięty, n (%)			9.304	0.01
L3	6 (6.90)	25 (15.53)
L4	48 (55.17)	101 (62.73)
L5	33 (37.93)	35 (21.74)
Stopień według Pfirrmanna, n (%)			0.99	0.804
Stopień II	5 (5.75)	13 (8.07)
Stopień III	24 (27.59)	50 (31.06)
Stopień IV	34 (39.08)	59 (36.65)
Stopień V	24 (27.59)	39 (24.22)
Stopień według Goutalliera, n (%)			9.446	0.024
Stopień 1	22 (25.29)	56 (34.78)
Stopień 2	26 (29.89)	60 (37.27)
Stopień 3	21 (24.14)	31 (19.25)
Stopień 4	18 (20.69)	14 (8.70)
Kąt stawu międzywyrostkowego dystalnego (°)	55.38 ± 6.52	55.20 ± 6.39	0.208	0.835
Wysięk ze stawu międzywyrostkowego dystalnego (mm)	0.96 ± 0.22	0.93 ± 0.20	1.091	0.277
Wysokość krążka dystalnego (mm)	9.08 ± 1.96	9.20 ± 1.90	0.469	0.64
Kąt stawu międzywyrostkowego proksymalnego (°)	49.82 ± 6.17	50.02 ± 5.97	0.252	0.802
Wysięk ze stawu międzywyrostkowego proksymalnego (mm)	0.85 ± 0.27	0.83 ± 0.25	0.584	0.56
Wysokość krążka proksymalnego (mm)	8.02 ± 1.56	8.07 ± 1.50	0.247	0.805
Inklunacja miednicy (PI, °)	52.46 ± 8.12	52.03 ± 7.86	0.404	0.687
Nachylenie miednicy (PT, °)	24.58 ± 7.33	24.22 ± 7.10	0.375	0.708
Kąt nachylenia kości krzyżowej (SS, °)	44.78 ± 6.91	44.98 ± 6.80	0.224	0.823
Kifoza piersiowa (TK, °)	21.08 ± 7.79	21.24 ± 7.62	0.158	0.875

Tabela 2: Analiza jednowymiarowa czynników powiązanych z lipomatozą znadtwardówkową w kohorcie modelowej. Zmienne ciągłe przedstawiono jako średnia ± odchylenie standardowe (SD), a zmienne kategoryczne jako liczba (%). Porównania przeprowadzono pomiędzy grupą z lipomatozą (n = 87) a grupą bez lipomatozy (n = 161) w celu zidentyfikowania zmiennych kandydujących do wieloczynnikowej analizy regresji logistycznej. Skróty: BMI, wskaźnik masy ciała; PI, inklinacja miedniczna; PT, tilt miedniczny; SS, kąt nachylenia kości krzyżowej; TK, kifoza piersiowa.

Wieloczynnikowa binarna regresja logistyczna pozwoliła zidentyfikować pięć niezależnych predyktorów lipomatozy znadtwardówkowej: wiek (iloraz szans [OR] = 1,510 na każde 10 lat wzrostu), płeć żeńska (OR = 2,354), BMI (OR = 1,874 na każde 5 kg/m² wzrostu), przesunięcie segmentu L5 (OR = 3,766) oraz stopnie Goutalliera 3–4 (OR = 3,184). Odpowiednie współczynniki regresji, błędy standardowe, statystyki Walda, wartości OR oraz 95% przedziały ufności (CI) podsumowano w Tabeli 3. Względne wielkości efektów oraz przedziały CI dla wybranych predyktorów przedstawiono na wykresie leśnym (Rysunek 1).

Zmienna	β	SE	χ-kwadrat Walda2	Wartość p	LUB	95% CI (przedział ufności)
Wiek	0.412	0.156	6.975	0.008	1.51	1.112-2.050
Płeć (żeńska)	0.856	0.368	5.408	0.02	2.354	1.144-4.842
Wskaźnik masy ciała	0.628	0.184	11.648	0.001	1.874	1.307-2.688
Przesunięcie segmentu (L5)	1.326	0.342	15.034	<0.001	3.766	1.926-7.362
stopień Goutalliera (stopnie 3-4)	1.158	0.316	13.421	<0.001	3.184	1.714-5.915
Stała	-5.521	1.248	19.574	<0.001	0.004	—

Tabela 3: Wieloczynnikowa analiza regresji logistycznej binaryjnej w celu wyłonienia niezależnych predyktorów lipomatozy nadtwardówkowej. Wyniki przedstawiono jako współczynnik regresji (β), błąd standardowy (SE), statystykę Wald χ², iloraz szans (OR) oraz odpowiadający mu 95% przedział ufności (CI) dla niezależnych predyktorów utrzymanych w końcowym wieloczynnikowym modelu regresji logistycznej binaryjnej. Podany współczynnik dla wieku odpowiada wzrostowi o 10 lat, a podany współczynnik dla wskaźnika masy ciała (BMI) odpowiada wzrostowi o 5 kg/m². Do bezpośredniego wprowadzenia surowych wartości klinicznych do równania predykcyjnego zastosowano współczynniki równoważne wynoszące 0,0412 na rok wieku oraz 0,1256 na kg/m² BMI. Skróty: β, współczynnik regresji; SE, błąd standardowy; OR, iloraz szans; CI, przedział ufności.

Wykresy kalibracji wykazały przybliżoną zgodność między prawdopodobieństwami przewidywanymi a obserwowanymi w kohortach modelowania i walidacji, przy czym w kohorcie walidacyjnej odnotowano większą niepewność przy wyższych przewidywanych prawdopodobieństwach (Rysunek 2A,B).

Dyskryminację modelu oceniono za pomocą analizy charakterystyki operacyjnej odbiornika (ROC). Model predykcyjny osiągnął wartości AUC wynoszące 0,834 w kohorcie modelowania i 0,815 w kohorcie walidacyjnej (Rycina 3A). Analiza krzywej decyzyjnej sugerowała potencjalną korzyść netto w wybranych zakresach prawdopodobieństwa progowego w porównaniu ze strategiami leczenia wszystkich (treat-all) i leczenia żadnego (treat-none) (Rycina 3B). Ponieważ kohorta walidacyjna została wydzielona wewnętrznie z jednej instytucji, wyniki te wskazują na spójność wewnętrzną i nie powinny być interpretowane jako dowód na szeroką generalizowalność zewnętrzną.

Zbiory danych wygenerowane i analizowane w ramach niniejszego badania są dostępne na uzasadnioną prośbę skierowaną do autora korespondencyjnego (R.T.). Ze względu na prywatność pacjentów oraz instytucjonalne zasady zarządzania danymi, kliniczne i obrazowe dane na poziomie indywidualnym nie mogą zostać udostępnione publicznie. Zanonimizowane dane podsumowujące, współczynniki regresji, równanie predykcyjne, przykładowy arkusz kodowania oraz kod analizy statystycznej mogą zostać udostępnione wykwalifikowanym badaczom po uzyskaniu zgody Rady Przeglądu Instytucjonalnego i Komisji Etyki Tianjin Union Medical Center oraz podpisaniu umowy o wykorzystaniu danych.

Dyskusja

Niniejszy protokół integruje standaryzowane pozyskiwanie obrazów MRI, ustrukturyzowaną interpretację obrazu oraz ilościowe modelowanie ryzyka w spójne narzędzie wspomagania decyzji klinicznych służące do identyfikacji lipomatozy nadtwardówkowej u pacjentów z DLS. Proces obejmuje cztery następujące po sobie komponenty: systematyczną charakterystykę pacjenta, obejmującą zmienne demograficzne i antropometryczne; standaryzowane pozyskiwanie obrazów MRI w płaszczyźnie strzałkowej w obrazowaniu T1, zoptymalizowane pod kątem wizualizacji tłuszczu nadtwardówkowego; ustrukturyzowaną jakościową ocenę odkładania się tłuszczu nadtwardówkowego oraz nacieku tłuszczowego mięśni przykręgosłupowych przy użyciu zdefiniowanych kryteriów diagnostycznych; oraz zastosowanie wielowymiarowego równania predykcyjnego w celu wygenerowania indywidualnych szacunków prawdopodobieństwa. Krytyczne kroki zapewniające powtarzalność protokołu obejmują wykorzystanie obrazowania T1 w płaszczyźnie strzałkowej jako głównej sekwencji diagnostycznej, próg zasięgu w kierunku czaszowo-ogonowym ≥5 mm w tej samej płaszczyźnie w celu zdefiniowania powtarzalnego złogu tłuszczu nadtwardówkowego na poziomie przesunięcia, potwierdzenie w płaszczyźnie osiowej ucisku worka oponowego, systematyczną ocenę nacieku tłuszczowego mięśnia wielodzielnego według skali Goutalliera dostosowanej do kręgosłupa oraz integrację tych parametrów w ilościowej ocenie ryzyka.

Protokół można dostosować do różnych warunków klinicznych i platform obrazowania. W systemach 3.0 T możliwe może być zastosowanie cieńszych warstw; jednak czas echa należy starannie dostosować, aby zminimalizować artefakty prześwitu T1 (T1 shine-through), które mogą być bardziej wyraźne przy wyższych natężeniach pola. W systemach 1.5 T standardowe parametry akwizycji zapewniają odpowiednią jakość diagnostyczną; w przypadku nadmiernego szumu obrazu należy zwiększyć liczbę wzbudzeń z 2 do 3 lub zastosować współczynniki przyspieszenia obrazowania równoległego nieprzekraczające 2. W przypadku nietypowych wzorców rozmieszczenia tkanki tłuszczowej, w tym rozmieszczenia obwodowego zamiast grzbietowego lub złogów plamistych zamiast pasmowych, należy zastosować kryterium ≥5 mm do największego ciągłego lub półciągłego złogu w tej samej płaszczyźnie strzałkowej i potwierdzić znalezisko na obrazach osiowych ważonych T2. Jeśli nadal utrzymuje się niepewność co do różnicy między patologiczną lipomatozą nadtwardówkową a fizjologiczną tkanką tłuszczową w przestrzeni nadtwardówkowej, należy porównać złóg tłuszczowy z średnicą przednio-tylną i powierzchnią przekroju poprzecznego worka oponowego na obrazach osiowych. W niniejszym protokole ilościowe pomiary osiowe służą wsparciu potwierdzenia i dokumentacji, a nie jako niezależny standard referencyjny do redefinicji pierwotnego wyniku modelu. Aby zminimalizować zmienność między obserwatorami w skali Goutalliera, należy zestandaryzować rozmieszczenie obszaru zainteresowania (ROI) na warstwie osiowej w połowie wysokości przestrzeni wysuniętego krążka międzykręgowego. Rozbieżności między dwoma niezależnymi obserwatorami należy rozstrzygać poprzez konsensus z trzecim, starszym stażem obserwatorem.

Należy wziąć pod uwagę kilka ograniczeń metodologicznych. Po pierwsze, diagnoza opiera się na jakościowej ocenie wizualnej charakterystyk sygnału MRI oraz zasięgu kraniokaudalnego, a nie na ilościowym pomiarze objętościowym. Chociaż podejście to jest praktyczne i szeroko stosowane, nie obejmuje ono trójwymiarowej objętości tłuszczu nadtwardówkowego, która może korelować bardziej bezpośrednio z nasileniem objawów i pilnością operacyjną. Ustanowione podejścia do graduacji w MRI, w tym metody oparte na stosunku tłuszczu nadtwardówkowego do worka oponowego oraz system graduacji lokoregionalnej Manjila, dostarczają uzupełniających informacji na temat ucisku kanału w przekroju osiowym, zasięgu strzałkowego oraz grzbietowego, brzusznego lub okrężnego rozmieszczenia tłuszczu^13,14. Cechy te mogą być szczególnie istotne na poziomie L5-S1, gdzie geometria kanału kręgowego i anatomia recesji bocznej różnią się od tych na innych poziomach lędźwiowych. Niniejszy protokół prospektywnie dokumentuje ucisk osiowy i rozmieszczenie tłuszczu, ale nie wykorzystuje tych systemów graduacji jako niezależnych standardów referencyjnych, ponieważ oryginalne badanie nie zostało zaprojektowane jako badanie dokładności diagnostycznej względem standardu referencyjnego. W konsekwencji nie przeliczano czułości, swoistości ani zgodności między metodami, aby uniknąć post hoc reklasyfikacji wyników oryginalnego badania. Zaawansowane ilościowe techniki MRI, w tym separacja tłuszczu i wody w oparciu o metodę Dixon oraz segmentacja wolumetryczna, mogą zapewnić dokładniejszą kwantyfikację tłuszczu, ale wymagają specjalistycznego oprogramowania i dodatkowego czasu akwizycji. Po drugie, system graduacji Goutalliera pozostaje półilościowy i jest podatny na zmienność międzyobserwacyjną. Choć pięciostopniowy schemat zapewnia praktyczną klasyfikację kliniczną, nie oferuje on ciągłego pomiaru nacieku tłuszczowego, który był powiązany z ekspresją cytokin zapalnych w mięśniach i tkance tłuszczowej nadtwardówkowej²¹ oraz ze zmienionymi pasywnymi właściwościami mechanicznymi mięśnia wielodzielnego^22,23. W badaniu tym zastosowano adaptację koncepcji Goutalliera do kręgosłupa, a stopnie 0–2 pogrupowano oddzielnie od stopni 3–4, aby odróżnić łagodne od ciężkiego zastąpienia tłuszczowego, zachowując przy tym parsymonię modelu. Po trzecie, model predykcyjny obejmuje tylko zmienne kliniczne i obrazowe i nie zawiera markerów biochemicznych, takich jak profile lipidowe surowicy, markery zapalne czy adipokiny, które mogłyby poprawić wydajność predykcyjną w przypadkach o podłożu metabolicznym. Po czwarte, zastosowanie przesiewowej analizy jednowymiarowej, po której następuje krokowa regresja logistyczna, może wprowadzić błąd selekcji zmiennych i niestabilność współczynników; dlatego model należy uznać za eksploracyjny i jedynie wewnętrznie zwalidowany. Wreszcie, przed szerokim wdrożeniem klinicznym konieczna pozostaje zewnętrzna walidacja w różnych warunkach klinicznych, na różnych platformach obrazowania i w różnych populacjach pacjentów. Ponieważ walidacja została przeprowadzona na kohorcie z jednej instytucji, przyszłe wieloośrodkowe badania prospektywne powinny potwierdzić wydajność modelu i ustalić klinicznie użyteczne progi prawdopodobieństwa.

Niniejszy protokół oferuje kilka zalet w porównaniu z ustrukturyzowaną praktyką kliniczną. Tradycyjna praktyka często opiera się na przypadkowym rozpoznaniu tłuszczowatości nadtwardówkowej przez radiologa bez zastosowania zestandaryzowanych kryteriów diagnostycznych, co prowadzi do niespójności w wykrywaniu tej zmiany. W przeciwieństwie do tego, niniejszy protokół dostarcza wyraźnych kryteriów diagnostycznych i ustrukturyzowanego schematu interpretacji, który może być stosowany konsekwentnie przez różnych lekarzy i w różnych placówkach. W porównaniu z bardziej złożonymi strategiami oceny multimodalnej, obejmującymi zaawansowane obrazowanie lub inwazyjne techniki diagnostyczne, obecne podejście jest praktyczne z punktu widzenia operacyjnego, ponieważ wykorzystuje standardowe sekwencje MRI odcinka lędźwiowego oraz rutynowo dostępne dane kliniczne. Jego głównym wkładem jest integracja dostosowanej do kręgosłupa oceny Goutalliera w zakresie naciekania tłuszczowego mięśnia multifidus z ramami przewidywania ryzyka tłuszczowatości nadtwardówkowej związanej z DLS. W ocenie degeneracji krążków międzykręgowych system stopniowania Pfirrmanna pozostaje praktyczną i powszechnie akceptowaną klasyfikacją morfologiczną w MRI¹⁸. Elastografia MR pojawiła się jako uzupełniająca technika ilościowa, ponieważ sztywność krążka wzrasta wraz z postępującą degeneracją i może dostarczać obiektywnych informacji biomechanicznych wykraczających poza porządkowe stopnie Pfirrmanna²⁴. Niniejszy protokół zachowuje stopniowanie według Pfirrmanna, ponieważ jest ono łatwo stosowalne w rutynowych badaniach MRI odcinka lędźwiowego, podczas gdy przyszłe wersje mogą włączyć elastografię MR w miarę zwiększania dostępności tej technologii.

Klinicznie szacunki prawdopodobieństwa mogą pomóc w uświadamianiu o ryzyku, analizie obrazowej, dyskusjach dotyczących planowania chirurgicznego oraz w poradnictwie pacjentów. Wcześniejsze prace wykazały porównywalne wyniki dwuletnie w przypadku lipomatozy nadtwardówkowej i stenozy degeneracyjnej po dekompresji^25,26. Jednakże niniejszy protokół nie ocenia wyników leczenia ani nie porównuje strategii dekompresji i nie powinien być wykorzystywany do nakazywania dekompresji metodą otwartą, małoinwazyjną lub endoskopową. Zamiast tego model ma za zadanie ostrzec klinicystów przed możliwością wystąpienia klinicznie istotnego gromadzenia się tłuszczu w przestrzeni nadtwardówkowej oraz zachęcać do uważnej analizy kompresji osiowej, grzbietowego lub brzusznego rozkładu tłuszczu oraz ucisku worka oponowego. Hipoplazja trzonu kręgu L5 oraz pseudolisthesis L5-S1 są potencjalnymi imitacjami rozwojowymi lub anatomicznymi, ponieważ stosunkowo mały wymiar przednio-tylny trzonu kręgu L5 oraz tylna tkanka tłuszczowo-włóknista mogą symulować pozorne przesunięcie kręgów. Aby zredukować to źródło błędnej klasyfikacji, protokół wyklucza niedegeneracyjną i rozwojową spondylolistezę, potwierdza DLS za pomocą zdjęć rentgenowskich w pozycji stojącej oraz braku defektów łuku kręgu i wymaga korelacji poziomu przesunięcia pomiędzy obrazami MRI w płaszczyźnie strzałkowej i osiowej. W poradnictwie pacjentów indywidualne szacunki prawdopodobieństwa mogą ułatwić dyskusję nad wynikami obrazowania i przewidywanymi kwestiami śródoperacyjnymi. Po operacji protokół może pomóc w zidentyfikowaniu pacjentów wymagających ściślejszej obserwacji, choć zastosowanie to wymaga dalszej walidacji. Poza DLS, opracowany schemat może być adaptowalny do innych degeneracyjnych schorzeń kręgosłupa, w których gromadzenie się tłuszczu w przestrzeni nadtwardówkowej przyczynia się do objawów, w tym do stenozy kanału kręgowego odcinka lędźwiowego oraz spondylozy lędźwiowej bez listhesis.

Podsumowując, niniejszy protokół zapewnia ustrukturyzowane i praktyczne klinicznie podejście do oceny lipomatozy nadtwardówkowej u pacjentów z DLS. Dzięki integracji znormalizowanej interpretacji MRI, systematycznej oceny mięśni przykręgosłowych, prospektywnego potwierdzenia ucisku worka twardówki w przekroju osiowym oraz wieloczynnikowego modelowania ryzyka, protokół ten umożliwia powtarzalną analizę obrazową, komunikację ryzyka, dyskusje nad planowaniem chirurgicznym oraz poradnictwo dla pacjentów przy użyciu standardowego klinicznego MRI bez konieczności stosowania specjalistycznego sprzętu. Pierwotna definicja wyniku diagnostycznego oraz zgłoszona wydajność modelu zostały zachowane podczas rewizji. Model należy interpretować jako wewnętrznie zwalidowane narzędzie wspomagania decyzji klinicznych, a nie jako diagnostyczny standard odniesienia lub regułę decyzji chirurgicznej. Przyszłe badania powinny obejmować porównanie z uznanymi standardami odniesienia MRI, wieloośrodkową walidację zewnętrzną, zautomatyzowaną analizę obrazu oraz zaawansowane ilościowe techniki MRI, takie jak obrazowanie Dixon i elastografia MR.

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
Skaner do obrazowania metodą rezonansu magnetycznego 1.5 Tesla	Sprzęt	Siemens Healthineers	Magnetom Aera
16-kanałowa cewka powierzchniowa do kręgosłupa	Sprzęt	Siemens Healthineers	Konfiguracja standardowa
Monitor diagnostyczny zgodny ze standardem DICOM Part 14 (3 megapiksele)	Sprzęt	Barco NV	Nio Color 3MP, MDNC-3421CN; kalibracja DICOM za pomocą Barco QAWeb Enterprise
Pakiet forestplot do środowiska R	Oprogramowanie	R Foundation for Statistical Computing	Wersja 3.1.3
ImageJ	Oprogramowanie	National Institutes of Health	Wersja 1.53t
Stacja robocza PACS	Sprzęt	Dell Inc.	Precision 3660 Tower Workstation; system operacyjny Windows 10, 64-bitowy
System archiwizacji i transmisji obrazów (PACS)	Oprogramowanie	INFINITT Healthcare Co., Ltd.	INFINITT PACS 7.0
R	Oprogramowanie	R Foundation for Statistical Computing	Wersja 4.3.1
REDCap	Oprogramowanie	Vanderbilt University	Wersja 13.1
Pakiet rmda do środowiska R	Oprogramowanie	R Foundation for Statistical Computing	Wersja 1.6
SPSS Statistics	Oprogramowanie	IBM	Wersja 27.0