Ogólne informacje o pacjentach
W badaniu wzięło udział łącznie 101 pacjentów. Jak pokazano w Tabeli 1, pacjenci w grupie z niekorzystnym wynikiem byli istotnie starsi (P < 0,05) i mieli dłuższy czas trwania choroby (P < 0,05). Nie zaobserwowano statystycznie istotnych różnic między obiema grupami pod względem płci, BMI, historii palenia tytoniu, chorób współistniejących, poziomu operowanego ani dostępu chirurgicznego (P > 0,05).
| Zmienne | Grupa pomyślna
(n = 76) | Grupa niekorzystna
(n = 25) | Proszę o podanie tekstu źródłowego do tłumaczenia.χ²/Z | P |
| Wiek (lata) | 45.2 ± 11.8 | 51.6 ± 12.9 | 2.322 | 0.022 |
| Samiec | 42 (55.3) | 16 (64.0) | 0.587 | 0.444 |
| BMI (kg/m²)2) | 24.4 ± 3.0 | 25.4 ± 3.6 | 1.531 | 0.129 |
| Czas trwania choroby (miesiące) | 9 (8, 13) | 15 (12, 17) | 4.174 | <0.001 |
| Wywiad tytoniowy | 18 (23.7) | 9 (36.0) | 1.457 | 0.227 |
| Historia cukrzycy | 8 (10.5) | 5 (20.0) | 0.779 | 0.337 |
| Historia nadciśnienia tętniczego | 12 (15.8) | 6 (24.0) | 0.396 | 0.529 |
| Segment operacyjny | | | | |
| L3/4 | 11 (14.5) | 1 (4.0) | | |
| L4/5 | 40 (52.6) | 16 (64.0) | | |
| L5/S1 | 25 (32.9) | 8 (32.0) | 2.18 | 0.336 |
| Podejście operacyjne | | | | |
| Dostęp transforaminalny | 59 (77.6) | 19 (76.0) | | |
| Metoda interlaminarna | 17 (22.4) | 6 (24.0) | 0.029 | 0.866 |
Tabela 1: Charakterystyka wyjściowa [średnia±SD, n (%), M(Q1, Q3)]. Charakterystyka wyjściowa z podziałem na grupy prognostyczne. Dane przedstawiono jako średnia ± SD, n (%) lub mediana (Q1, Q3). Grupa niepomyślna (Unfavorable) była istotnie starsza (P = 0.022) i charakteryzowała się dłuższym czasem trwania choroby (P < 0.01). Pozostałe zmienne nie wykazały istotnych różnic międzygrupowych. BMI, wskaźnik masy ciała.
Porównanie przedoperacyjnych parametrów morfologicznych MRI
Przedoperacyjne parametry morfologiczne MRI u 101 pacjentów przedstawiono w Tabeli 2. Odnotowano istotne różnice (P < 0,05) między dwiema grupami w zakresie proporcji zmian Modica oraz typów przepukliny; wskaźnik DHI oraz powierzchnia otworu międzykręgowego w grupie korzystnej były istotnie wyższe niż w grupie niekorzystnej (P < 0,05), natomiast współczynnik zajęcia kanału kręgowego oraz maksymalna powierzchnia przepukliny w grupie korzystnej były istotnie niższe niż w grupie niekorzystnej (P < 0,05); nie stwierdzono statystycznie istotnych różnic (P > 0,05) w stopniu wg klasyfikacji Pfirrmanna, lokalizacji przepukliny oraz stopniu migracji. Wyniki sugerują, że powyższe parametry mogą być potencjalnymi wskaźnikami prognostycznymi dla wyników pooperacyjnych zabiegu PELD.
| Grupa korzystna
(n = 76) | Grupa niekorzystna (n = 25) | Proszę podać tekst źródłowy do tłumaczenia.χ² | P |
| skala Pfirrmanna | | | | |
| I-II | 28(36.8) | 5(20.0) | | |
| III | 32(42.1) | 10(40.0) | | |
| IV-V | 16(21.1) | 10(40.0) | 4.277 | 0.118 |
| Zmiany Modicia | | | | |
| Proszę podać tekst źródłowy do tłumaczenia. | 48(63.2) | 8(32.0) | | |
| Typ I | 12(15.8) | 3(12.0) | | |
| Typ II | 14(18.4) | 13(52.0) | | |
| Typ III | 2(2.6) | 1(4.0) | 11.529 | 0.009 |
| Typ przepukliny | | | | |
| Wypuklenie | 8(10.5) | 1(4.0) | | |
| Wypukłość | 44(57.9) | 7(28.0) | | |
| Wytłaczanie | 20(26.3) | 12(48.0) | | |
| Sekwestracja | 4(5.3) | 5(20.0) | 11.605 | 0.009 |
| Lokalizacja przepukliny | | | | |
| Centralny | 12(15.8) | 2(8.0) | | |
| Paracentralny | 48(63.2) | 18(72.0) | | |
| Otworowy | 12(15.8) | 4(16.0) | | |
| Pozakanałowy | 4(5.3) | 1(4.0) | 1.11 | 0.775 |
| Stopień migracji | | | | |
| Brak migracji (stopień 0) | 25(32.9) | 6(24.0) | | |
| Migracja niskiego stopnia (stopień 1–2) | 32(42.1) | 10(40.0) | | |
| Migracja wysokiego stopnia (stopień 3–4) | 19(25.0) | 9(36.0) | 1.326 | 0.515 |
| Stopień zajęcia kanału kręgowego (%) | 35.6 ± 6.9 | 47.2 ± 8.4 | 6.892 | <0.001 |
| Wskaźnik zajęcia kanału kręgowego > 40% | 30(39.5) | 14(56.0) | | |
| DHI (%) | 28.5 ± 5.2 | 24.5 ± 4.8 | 3.392 | 0.001 |
| DHI > 25% | 58(76.3) | 11(44.0) | | |
| Pole otworu (cm²)²) | 1.1 ± 0.3 | 0.8 ± 0.3 | 4.029 | <0.001 |
| Obszar otworu wyjścia nerwu > 1 cm² | 49(64.5) | 13(52.0) | | |
| Maksymalna powierzchnia przepukliny (mm2)²) | 59.8 ± 11.9 | 86.0 ± 15.7 | 8.764 | <0.001 |
| Maksymalna powierzchnia przepukliny > 75 mm² | 21(27.6) | 13(52.0) | | |
Tabela 2: Porównanie przedoperacyjnych morfologicznych parametrów MRI [średnia±SD, n (%)]. Porównanie parametrów MRI przed operacją pomiędzy grupami. Dane przedstawiono jako średnia ± SD lub n (%). Istotne różnice zaobserwowano w zakresie zmian Modica, typu przepukliny, stopnia migracji, współczynnika zajęcia kanału, DHI, pola otworu międzykręgowego oraz pola przepukliny (wszystkie P < 0,05), natomiast nie stwierdzono istotnych różnic w stopniu wg skali Pfirrmanna ani lokalizacji przepukliny. DHI, wskaźnik wysokości krążka międzykręgowego (Disc Height Index).
Zmiany we wskaźnikach wyników klinicznych
Tabela 3 przedstawia wyniki kliniczne 101 pacjentów ocenianych przy użyciu zmodyfikowanych kryteriów MacNaba podczas wizyty kontrolnej po 1 roku. Wyniki wykazały doskonałe efekty u 51 pacjentów (50,5%) i dobre efekty u 25 pacjentów (24,8%), co łącznie daje 76 pacjentów (75,2%) w grupie z wynikiem korzystnym. Wyniki umiarkowane zaobserwowano u 16 pacjentów (15,8%), a niekorzystne rokowania u 9 pacjentów (8,9%), co stanowi 25 pacjentów (24,8%) w grupie z wynikiem niekorzystnym, co wskazuje na ogólnie korzystną skuteczność kliniczną.
| Stopień rezultatu | Liczba przypadków (n) | Procent (%) | Definicja wyniku klinicznego |
| Doskonały | 51 | 50.50% | Grupa z korzystnym rokowaniem |
| Dobrze | 25 | 24.80% | Grupa korzystna |
| Uczciwy | 16 | 15.80% | Grupa niekorzystna |
| Słabe | 9 | 8.90% | Grupa niekorzystna |
Tabela 3: Roczna ocena skuteczności klinicznej po operacji na podstawie zmodyfikowanych kryteriów MacNaba. Wyniki po roku od operacji według zmodyfikowanych kryteriów MacNaba. Wyniki doskonałe/dobre stanowiły grupę korzystną (75,2%), natomiast wyniki przeciętne/słabe stanowiły grupę niekorzystną (24,8%).
Tabela 4 przedstawia dynamiczne zmiany w pooperacyjnych wynikach klinicznych. Przed operacją nie stwierdzono statystycznie istotnych różnic między dwiema grupami w zakresie wskaźnika VAS dla bólu odcinka lędźwiowego kręgosłupa, VAS dla bólu nogi ani w wynikach wskaźnika ODI (P > 0,05). We wszystkich punktach czasowych po operacji wyniki w obu grupach wykazały znaczną poprawę w porównaniu z wartościami przedoperacyjnymi (wewnątrzgrupowo P < 0,01), co wskazuje na ogólną skuteczność PELD. Jednak porównania międzygrupowe wykazały, że grupa korzystna osiągnęła istotnie lepsze wyniki VAS (plecy i noga) oraz ODI we wszystkich punktach czasowych po operacji w porównaniu z grupą niekorzystną (P < 0,01). Dodatkowo poprawa wskaźnika ODI w roku po zabiegu była istotnie większa w grupie korzystnej niż w grupie niekorzystnej. Wyniki te dowodzą, że choć sama operacja prowadzi do znacznej poprawy, pacjenci podzieleni na grupy na podstawie przedoperacyjnych cech morfologicznych w MRI wykazują trwałe różnice w trajektoriach powrotu do zdrowia, co dodatkowo potwierdza wartość prognostyczną cech MRI przedoperacyjnego dla długoterminowych wyników.
| Czas | Grupa z korzystnym wynikiem
(n = 76) | Grupa z niekorzystnym wynikiem
(n = 25) | P (wewnątrzgrupowo) | P (międzygrupowo) |
| VAS-plecy | T0 | 6(6, 7) | 7(6, 7) | - | 0.212 |
| T1 | 3(3, 4) | 4(4, 5) | <0.01 | <0.01 |
| T2 | 2(1, 2) | 3(2, 4) | <0.01 | <0.01 |
| T3 | 1(1, 1) | 2(2, 3) | <0.01 | <0.001 |
| VAS-noga | T0 | 8(7, 9) | 8(7, 8) | - | 0.486 |
| T1 | 3(2, 3) | 5(4, 5) | <0.001 | <0.01 |
| T2 | 2(1, 2) | 3(3, 4) | <0.01 | <0.01 |
| T3 | 1(1, 1) | 3(2, 3) | <0.01 | <0.01 |
| ODI (%) | T0 | 68.4 ± 9.6 | 70.2 ± 10.4 | - | 0.408 |
| T1 | 30.1 ± 7.4 | 42.4 ± 8.9 | <0.01 | <0.01 |
| T2 | 18.7 ± 5.5 | 29.5 ± 6.4 | <0.01 | <0.01 |
| T3 | 10.1 ± 4.0 | 21.8 ± 5.5 | <0.001 | <0.01 |
| Zmiana wyniku ODI | T3-T0 | 58.2 ± 10.9 | 48.5 ± 1.5 | - | <0.01 |
Tabela 4: Zmiany wskaźnikach wyników klinicznych przed i po operacji [średnia ± SD, M(Q1, Q3)]. Podłużne zmiany w punktacji VAS i ODI. Dane przedstawiono jako medianę (Q1, Q3) lub średnią ± SD. W obu grupach odnotowano istotną poprawę po operacji (wewnątrzgrupowo P < 0.01), jednak grupa Favorable wykazywała konsekwentnie lepsze wyniki we wszystkich punktach czasowych (międzygrupowo P < 0.01). Porównania wewnątrz grup przeprowadzono za pomocą testu znaków Wilcoxona, natomiast porównania między grupami za pomocą testu U Manna-Whitneya.
Analiza korelacji między morfologicznymi parametrami MRI a wynikami klinicznymi Tabela 5 przedstawia analizę korelacji Pearsona, która wykazała, że współczynnik zajęcia kanału kręgowego oraz maksymalna powierzchnia przepukliny były istotnie dodatnio skorelowane z pooperacyjnymi wynikami VAS dla bólu pleców, VAS dla bólu nóg oraz wynikami ODI (r: 0,375–0,538, P < 0,01). Wskazuje to, że większy stopień zwężenia kanału kręgowego i większa objętość przepukliny wiążą się z bardziej nasilonymi objawami resztkowymi i gorszą regeneracją funkcjonalną. Z kolei DHI oraz powierzchnia otworu międzykręgowego wykazały istotne ujemne korelacje ze wszystkimi miernikami wyników klinicznych (r: od -0,637 do -0,240, P < 0,05). spośród nich korelacja między powierzchnią otworu międzykręgowego a wynikiem ODI była bliska istotności (P=0,053). Wyniki te sugerują, że zachowana wysokość dysku międzykręgowego i odpowiednia przestrzeń otworu międzykręgowego ułatwiają pooperacyjną dekompresję nerwów i łagodzenie objawów. Powyższe wyniki korelacji dodatkowo określają siłę powiązań między parametrami MRI a rokowaniami. DHI wykazało najsilniejszą korelację z VAS-Back (r = -0,637), co podkreśla jego szczególną wartość prognostyczną dla pooperacyjnej regeneracji funkcjonalnej i stanowi ilościową podstawę do przedoperacyjnej oceny obrazowej.
| VAS-Back | VAS-noga | ODI |
| Parametry MRI | r | P | Proszę podać tekst źródłowy do tłumaczenia. | P | r | P |
| Wskaźnik zajętości kanału kręgowego | 0.501 | <0.001 | 0.52 | <0.001 | 0.404 | <0.001 |
| DHI | -0.637 | <0.001 | -0.608 | <0.001 | -0.323 | 0.001 |
| Obszar otworu wyjścia nerwu | -0.24 | 0.016 | -0.298 | 0.002 | -0.193 | 0.053 |
| Maksymalna powierzchnia przepukliny | 0.475 | <0.001 | 0.538 | <0.001 | 0.375 | <0.001 |
Tabela 5: Analiza korelacji Pearsona pomiędzy morfologicznymi parametrami MRI a wskaźnikami wyników klinicznych w 1. roku po operacji.Korelacje Pearsona pomiędzy parametrami MRI a wynikami w 1. roku. Współczynnik zajęcia kanału oraz powierzchnia przepukliny korelowały dodatnio z wynikami VAS/ODI (r = 0.457–0.653), natomiast DHI i powierzchnia otworu korelowały ujemnie (r = –0.627 do –0.31), wszystkie P < 0.01.
Wieloczynnikowa analiza regresji logistycznej
Przyjmując zmodyfikowane kryteria MacNaba w 1. roku po operacji jako zmienną zależną (prognoza korzystna = 1, prognoza niekorzystna = 0), do wieloczynnikowej analizy regresji logistycznej włączono zmienne, dla których P < 0,1 w analizie jednoczynnikowej parametrów morfologicznych MRI (zmiany Modica, typ przepukliny, zakres migracji, wskaźnik zajęcia kanału kręgowego, DHI, powierzchnia otworu międzykręgowego oraz maksymalna powierzchnia przepukliny). Ze względu na współliniowość między typem przepukliny a maksymalną powierzchnią przepukliny oraz fakt, że maksymalna powierzchnia przepukliny jako zmienna ciągła dostarcza więcej informacji, do modelu włączono maksymalną powierzchnię przepukliny zamiast typu przepukliny. Jak pokazano w Tabeli 6, po korekcie wieloczynnikowej zmiany Modica typu II (OR = 0,213, 95% CI: 0,062–0,728, P = 0,014) zostały zidentyfikowane jako niezależny czynnik ryzyka niekorzystnej prognozy, podczas gdy DHI > 25% (OR = 5,154, 95% CI: 1,62–16,376, P = 0,05) był niezależnym czynnikiem ochronnym dla korzystnej prognozy. Tabela 6 wskazuje, że wartości VIF dla wszystkich współzmiennych są mniejsze niż 2, co świadczy o braku istotnego problemu współliniowości.
| | | | | | 95% CI dla OR | |
| B | SE | Wald | P | LUB | Niższy | Górna | VIF |
| Zmiany Modica typu I | -0.542 | 0.853 | 0.404 | 0.525 | 0.581 | 0.109 | 3.095 | 1.082 |
| Zmiany Modica typu II | -1.547 | 0.627 | 6.082 | 0.014 | 0.213 | 0.062 | 0.728 |
| Zmiany Modica typu III | -0.802 | 1.452 | 0.305 | 0.581 | 0.449 | 0.026 | 7.72 |
| Migracja niskostopniowa | -0.539 | 0.707 | 0.583 | 0.445 | 0.583 | 0.146 | 2.329 | 1.073 |
| Migracja wysokiego stopnia | -0.953 | 0.738 | 1.669 | 0.196 | 0.386 | 0.091 | 1.637 |
| Współczynnik zajęcia kanału kręgowego > 40% | -0.813 | 0.563 | 2.084 | 0.149 | 0.444 | 0.147 | 1.337 | 1.016 |
| DHI > 25% | 1.64 | 0.59 | 7.729 | 0.005 | 5.154 | 1.622 | 16.376 | 1.026 |
| Obszar otworu wyjścia nerwu > 1 cm² | 0.84 | 0.582 | 2.084 | 0.149 | 2.316 | 0.741 | 7.245 | 1.047 |
| Maksymalna powierzchnia przepukliny > 75 mm² | -1.153 | 0.614 | 3.526 | 0.06 | 0.316 | 0.095 | 1.052 | 1.154 |
| Stała | 1.596 | 0.782 | 4.173 | 0.041 | 4.936 | - | - | |
Tabela 6: Wieloczynnikowa analiza regresji logistycznej pooperacyjnych wyników klinicznych. Wieloczynnikowa regresja logistyczna dla predyktorów korzystnej prognozy. Zmiany Modica typu II stanowiły niezależny czynnik ryzyka (OR = 0.213, P = 0.014), podczas gdy DHI > 25% był czynnikiem ochronnym (OR = 5.154, P = 0.05). Wartość VIF < 2 wskazała na brak współliniowości (Hosmer-Lemeshow P = 0.933). Test Hosmera-Lemeshow χ2=3.017, P=0.93, Nagelkerke R2 = 0.353.
Analiza prognostyczna pacjentów z różnymi typami zmian Modica
Tabela 7 przedstawia wyniki prognostyczne pacjentów w podziale na różne typy zmian Modica. Wyniki VAS w zakresie 4–6 wskazują na umiarkowany ból, a wyniki ODI w zakresie 21%–40% wskazują na umiarkowaną niepełnosprawność. Pacjenci ze zmianami Modica typu II wykazali najniższy odsetek korzystnej prognozy (51,9%) oraz stosunkowo wysoki odsetek niepełnosprawności funkcjonalnej (3,3%).
| Zmiany Modica | n | Korzystna prognoza | Umiarkowany ból pleców | Umiarkowany ból nogi | Umiarkowana niepełnosprawność |
| Brak | 56 | 48(85.7) | 2(3.6) | 1(1.8) | 5(8.9) |
| I | 15 | 12(80.0) | 0 | 0 | 2(13.3) |
| II | 27 | 14(51.9) | 2(7.4) | 2(7.4) | 8(3.3) |
| III | 3 | 2(6.7) | 0 | 0 | 1(33.3) |
| χ²/F | - | 1.529 | 1.608 | 2.68 | 6.617 |
| P | - | 0.09 | 0.658 | 0.46 | 0.085 |
Tabela 7: Porównanie wyników klinicznych w zależności od różnych typów zmian Modica [n (%)].Wyniki kliniczne w zależności od typu zmian Modica. Dane przedstawiono jako n (%). Grupa Modic II wykazała najniższy odsetek korzystnych rokowań (51,9%) i najwyższy odsetek umiarkowanej niepełnosprawności (32,0%, P = 0,048).
DOSTĘPNOŚĆ DANYCH:
Wszystkie dane potwierdzające wyniki niniejszego badania zostały zawarte w artykule oraz w jego Pliku uzupełniającym 1.
Plik uzupełniający 1: Surowe dane. Kliknij tutaj, aby pobrać ten plik.