Charakterystyka pacjentów
Podstawowa charakterystyka 153 pacjentów z CRC, u których przeprowadzono resekcję radykalną, została szczegółowo przedstawiona w Tabeli 1. Mediana wieku w kohorcie wynosiła 66 lat (IQR: 17), a mediana BMI wyniosła 23,1 kg/m2 (IQR: 4,7). Palacze stanowili 39,9% (n = 61) pacjentów, natomiast 31,4% (n = 48) zgłosiło spożycie alkoholu. Rozkład stopnia zaawansowania nowotworu przedstawiał się następująco: stopień I, 22,2% (n = 34); stopień II, 47,1% (n = 72); stopień III, 22,9% (n = 35) oraz stopień IV, 7,8% (n = 12). Mediana czasu hospitalizacji wynosiła 8 dni (IQR: 3).
Stosunek czynników zapalnych w okresie okołooperacyjnym
W Tabeli 2 przedstawiono mediany i rozstępy międzykwartylowe (IQR) dla stosunków wartości pooperacyjnych do przedoperacyjnych zmierzonych czynników zapalnych. Analiza wykazała, że poziom IL–6 wzrósł najsilniej po operacji (mediana stosunku = 7,23, IQR = 28,31), a następnie IL–10 (mediana stosunku = 2,13, IQR = 4,18). W przeciwieństwie do nich, mediany stosunków dla IFN–γ (mediana stosunku = 0,89, IQR = 1,99), IL–12p70 (mediana stosunku = 0,88, IQR = 0,91), IL–17 (mediana stosunku = 0,82, IQR = 0,97) oraz TNF–α (mediana stosunku = 0,87, IQR = 1,69) były mniejsze niż 1,0, co wskazuje na spadki pooperacyjne. Warto zauważyć, że poziomy IL–1β, IL–2 i IL–8 pozostały stosunkowo stabilne, przy zgłoszonych medianach stosunków wynoszących 1,00 (IQR odpowiednio 0,79, 0,95 i 1,96).
Jedno- i wieloczynnikowe modele analizy regresji Coxa
Jednoczynnikowa analiza regresji proporcjonalnego hazardu Coxa wykazała, że kilka parametrów klinicznych jest istotnie powiązanych z długoterminowym przeżyciem. Konkretnie, wiek pacjenta (HR = 1,05, 95% CI = 1,00–1,10, p = 0,031), stopień zaawansowania guza (HR = 3,61, 95% CI = 2,18–6,00, p < 0,001), poziom IFN–α (HR = 1,00, 95% CI = 1,00–1,00, p = 0,023) oraz stosunek IL–1β (HR = 1,00, 95% CI = 1,00–1,01, p = 0,012) wykazało statystycznie istotne powiązania (p < 0,05). Ponadto poziom TNF–α (HR = 1,00, 95% CI = 1,00–1,00, p = 0,076) spełnił przyjęty wcześniej próg istotności (p < 0,1).
Zmienne z wartościami p < 0,1 w analizie jednowymiarowej (wiek, czas pobytu w szpitalu, stopień zaawansowania guza, IFN‑α_ratio, IL‑1β_ratio oraz TNF‑α_ratio; łącznie 21 zdarzeń) zostały następnie włączone do wielowymiarowego modelu regresji. Wskaźnik IL‑1β modelowano jako ciągły składnik liniowy bez transformacji w obu analizach. Po skorygowaniu o potencjalne czynniki zakłócające wyniki analizy wielowymiarowej potwierdziły, że zarówno wskaźnik IL–1β (skorygowany HR = 1,05, 95% CI = 1,00–1,10, p = 0,037), jak i stopień zaawansowania guza (skorygowany HR = 3,83, 95% CI = 2,27–6,46, p < 0,001) zachowały niezależną wartość prognostyczną (próg istotności p < 0,05). Każdy wzrost wskaźnika IL–1β o 1 jednostkę wiązał się z 5% wzrostem ryzyka śmiertelności (HR = 1,05, 95% CI = 1,00–1,10, p = 0,037). Jednowymiarowa wartość HR wynosząca 1,00 (95% CI 1,00–1,01) wynika z zaokrąglenia do dwóch miejsc po przecinku (rzeczywiste β ≈ 0,004), a zmiana na 1,05 w modelu wielowymiarowym jest przypisywana korekcie czynników zakłócających, a nie niespójności w specyfikacji modelu. Co istotne, choć IFN–α (p = 0,023) i wiek pacjenta (p = 0,031) wykazały istotność statystyczną w analizie jednowymiarowej, nie utrzymały one istotności w modelu wielowymiarowym. Ponadto żaden z pozostałych markerów zapalnych nie wykazał statystycznie istotnej wartości prognostycznej w żadnej z analiz. Szczegółowe wyniki statystyczne przedstawiono w Tabeli 3.
Wpływ stosunku Il-1β na przeżycie specyficzne dla nowotworu
Analiza Kaplana-Meiera u pacjentów podzielonych według stosunku IL–1β (z podziałem dychotomicznym na poziomie 90. percentyla wynoszącego 5,35, jak pokazano na Rysunku 2) wykazała trwałą przewagę przeżycia w grupie o niskim stosunku (<5,35; n = 136) w porównaniu z grupą o wysokim stosunku (≥5,35; n = 17) (log–rank p = 0,026), przy krzywych, które sukcesywnie się rozchodziły. Liczby te spadły odpowiednio do 121 i 12 w 40. miesiącu obserwacji (co odpowiada 16 zgonom w grupie o niskim stosunku i 5 zgonom w grupie o wysokim stosunku). Podczas obserwacji odnotowano łącznie 21 zgonów z powodu nowotworu; zgony niezwiązane z nowotworem zostały wykluczone zgodnie z kryteriami badania i nie były wliczane do punktu końcowego. Aby dodatkowo uwzględnić potencjalne zakłócenia związane ze stadium choroby, w badaniu przeprowadzono dodatkową post-hoc analizę przeżycia specyficznego dla nowotworu ograniczoną do pacjentów w stadiach I–III (n = 141, 16 zgonów). W tej podgrupie skorygowane HR dla stosunku IL–1β wynosiło 1,04 (95% CI 0,99–1,10, p = 0,11), co wskazuje na spójną, ale statystycznie nieistotną tendencję, odzwierciedlającą prawdopodobnie zmniejszoną moc statystyczną ze względu na mniejszą liczbę zdarzeń.
DOSTĘPNOŚĆ DANYCH:
Dane wspierające wyniki niniejszego badania są dostępne na żądanie autora korespondencyjnego. Zanonimizowane dane znajdują się w Tabeli uzupełniającej 1, a kod źródłowy jest dostępny w Pliku uzupełniającym 1.

Rysunek 1: Schemat blokowy selekcji pacjentów i analizy danych. Początkowo przesiewowo przebadano łącznie 429 pacjentów z rakiem jelita grubego (CRC), którzy zostali poddani resekcji radykalnej. Po wykluczeniu 276 pacjentów (193 z niepełnym monitorowaniem, w tym 11 zgonów niezwiązanych z rakiem, oraz 83 z brakiem danych dotyczących IL‑1β przed lub po operacji), do ostatecznej analizy włączono 153 pacjentów. Rysunek został przygotowany w programie Microsoft PowerPoint (Microsoft 365). Copyright 2025 Autorzy. Jest to oryginalny rysunek stworzony przez autorów. Wszelkie prawa zastrzeżone. Skrót: IL-1β = interleukina‑1 beta. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 2: Analiza Kaplana–Meiera przeżycia specyficznego dla nowotworu u pacjentów z CRC (n = 153) w podziale według stosunku IL‑1β po operacji do wartości przedoperacyjnej. Stosunek ten został podzielony dychotomicznie przy użyciu 90. percentyla wynoszącego 5,35: grupa o wysokim stosunku (≥5,35, n = 17) kontra grupa o niskim stosunku (<5,35, n = 136). Test log-rank wykazał p = 0,026, co wskazuje na istotną przewagę przeżycia w grupie o niskim stosunku. Zacieniowane obszary reprezentują 95% przedziały ufności (CI). Skróty: IL-1β = interleukina-1 beta; CRC = rak jelita grubego. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Czynniki | n = 153 |
| Wiek, rok | 66 (17) |
| BMI, kg/m²2 | 23.1 (4.7) |
| Palenie tytoniu | 61 (39.9%) |
| Picie | 48 (31.4%) |
| Stopień zaawansowania nowotworu | |
| I | 34 (22.2%) |
| II | 72 (47.1%) |
| III | 35 (22.9%) |
| IV | 12 (7.8%) |
| Hospitalizacje, liczba dni | 8 (3) |
Tabela 1: Informacje wyjściowe całej kohorty. Zmienne wyrażono jako medianę (IQR = Q3–Q1) lub n (%). Skróty: BMI = wskaźnik masy ciała; IQR = rozstęp międzykwartylny.
| Zmienna | Mediana | IQR (rozstęp międzykwartylowy) |
| stosunek IFN–α | 1.02 | 1.81 |
| stosunek_IFN–γ | 0.89 | 1.99 |
| stosunek IL–10 | 2.13 | 4.18 |
| stosunek_IL–12p70 | 0.88 | 0.91 |
| stosunek IL–17 | 0.82 | 0.97 |
| stosunek IL–1β | 1 | 0.79 |
| stosunek_IL–2 | 1 | 0.95 |
| stosunek_IL–4 | 0.85 | 2.23 |
| stosunek IL–5 | 1.06 | 1.17 |
| stosunek IL–6 | 7.23 | 28.31 |
| stosunek IL–8 | 1 | 1.96 |
| stosunek TNF-α | 0.87 | 1.69 |
Tabela 2: Stosunek czynników zapalnych w okresie okołooperacyjnym. Uwaga: IQR oznacza różnicę między Q3 a Q1. Skróty: IQR = rozstęp międzykwartylny; IFN–β = interferon-beta; IFN–γ = interferon-gamma; IL–1β = interleukina-1 beta; TNF–α = czynnik martwicy nowotworów alfa; IL–2_ratio = stosunek poziomu IL–2 po operacji do poziomu przed operacją.
| Zmienna | Analiza jednowymiarowa | Analiza wielowymiarowa |
| HR (95% CI) | wartość p | HR (95% CI) | wartość p |
| Wiek | 1.05 (1.00–1.10) | 0.031 | 1.04 (0.99–1.09) | 0.09 |
| BMI | 0.95 (0.83–1.09) | 0.471 | | |
| Palenie tytoniu | 0.92 (0.38–2.23) | 0.858 | | |
| Alkohol | 1.12 (0.45–2.78) | 0.805 | | |
| Stopień zaawansowania nowotworu | 3.61 (2.18–6.00) | < 0.001 | 3.83 (2.27–6.46) | < 0.001 |
| Pobyt w szpitalu | 1.03 (1.00–1.06) | 0.079 | 1.03 (0.99–1.07) | 0.151 |
| IFN–α_ratio | 1.00 (1.00–1.00) | 0.023 | 0.98 (0.96–1.01) | 0.254 |
| IFN–γ_ratio | 1.00 (0.98–1.02) | 0.981 | | |
| IL–10_ratio | 0.99 (0.96–1.02) | 0.524 | | |
| IL–12p70_ratio | 0.97 (0.83–1.13) | 0.683 | | |
| IL–17_ratio | 0.92 (0.77–1.11) | 0.385 | | |
| IL–1β_ratio | 1.00 (1.00–1.01) | 0.012 | 1.05 (1.00–1.10) | 0.037 |
| IL–2_ratio | 0.91 (0.72–1.16) | 0.451 | | |
| IL–4_ratio | 1.00 (0.99–1.01) | 0.955 | | |
| IL–5_ratio | 0.98 (0.92–1.04) | 0.524 | | |
| IL–6_ratio | 1.00 (1.00–1.00) | 0.614 | | |
| IL–8_ratio | 0.99 (0.94–1.04) | 0.593 | | |
| TNF–α_ratio | 1.00 (1.00–1.00) | 0.076 | 1.00 (0.99–1.00) | 0.323 |
Tabela 3: Analiza regresji Coxa dla wartości przeżycia specyficznego dla nowotworu. Skróty: HR = hazard ratio; 95% CI, 95% przedział ufności. Wartości pogrubione są istotne statystycznie (p < 0,1) lub marginalnie istotne (p < 0,05). Model wielowymiarowy obejmował sześć współzmiennych: wiek, czas pobytu w szpitalu, stadium guza, IFN‑α_ratio, IL‑1β_ratio oraz TNF‑α_ratio. Wskaźnik IL‑1β wprowadzono jako ciągły człon liniowy bez transformacji logarytmicznej.
Tabela uzupełniająca 1: Surowe dane. Zanonimizowany zbiór danych wykorzystany w badaniu.Prosimy kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 1: Kod źródłowy. Kod źródłowy wykorzystany w badaniu. Kliknij tutaj, aby pobrać ten plik.