W badaniu wykorzystano zanonimizowane dokumentacje kliniczne z Szpitala Ludowego nr 2 w Kunshan, a samo badanie zostało zatwierdzone przez Komisję Etyczną szpitala (numer zatwierdzenia ksehllsp2024-05). Badanie przeprowadzono zgodnie z zasadami Deklaracji Helsińskiej. Ze względu na to, że informacje o pacjentach zostały zanonimizowane przed analizą, komisja etyczna zwolniła badaczy z obowiązku uzyskania pisemnej świadomej zgody. Przepływ procesu selekcji pacjentów, ich grupowania oraz główna analiza zostały przedstawione w Rycina 1Sprzęt laboratoryjny, systemy baz danych oraz oprogramowanie statystyczne wykorzystane w niniejszym badaniu wymieniono w Tabela materiałów.

Rysunek 1Schemat blokowy przesiewu i grupowania pacjentów. Z systemu elektronicznej dokumentacji medycznej szpitala wstępnie zidentyfikowano łącznie 25 hospitalizowanych pacjentów z dokumentacją diagnostyczną dotyczącą podejrzenia lub potwierdzenia urazowego uszkodzenia mózgu (TBI) w latach 2017–2023. Po wykluczeniu pacjentów w wieku Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
Projekt badania i kohorta pacjentów
Przeanalizowano jednorodniową kohortę pacjentów hospitalizowanych z powodu TBI w okresie od stycznia 2017 do grudnia 2023 roku. Potencjalnie kwalifikujących się pacjentów zidentyfikowano w elektronicznym systemie dokumentacji medycznej szpitala. Przesiew oparto na zapisach diagnozy przy wypisie zawierających terminy diagnostyczne związane z TBI, a następnie zweryfikowano diagnozę za pomocą dokumentacji z hospitalizacji oraz wyników neuroobrazowania. Terminy diagnostyczne związane z TBI oparto na kodach Międzynarodowej Klasyfikacji Chorób, 10. rewizji (ICD-10): S06.0 (wstrząśnienie mózgu), S06.1 (urazowy obrzęk mózgu), S06.3 (ogniskowe uszkodzenie mózgu), S06.4 (krwiak nadtwardówkowy), S06.5 (urazowy krwiak podtwardówkowy), S06.6 (urazowy krwotok podpajęczynówkowy), S06.8 (uraz wewnątrzczaszkowy z przedłużonym śpiączką) oraz S06.9 (uraz wewnątrzczaszkowy, nieokreślony). Potwierdzoną pierwotną diagnozę TBI zdefiniowano jako obecność jednego lub więcej z tych kodów ICD-10 jako głównej diagnozy przy wypisie, popartej zgodnymi znaleziskami w tomografii komputerowej (CT) czaszki lub obrazowaniu metodą rezonansu magnetycznego (MRI). Łącznie pod kątem kwalifikacji oceniono 2 5 hospitalizowanych pacjentów z podejrzeniem lub potwierdzonym TBI.
W przypadku pacjentów z więcej niż jedną hospitalizacją w okresie badania, za hospitalizację indeksową uznano jedynie pierwsze kwalifikujące się przyjęcie. Krok ten zastosowano w celu uniknięcia powielonych rekordów na poziomie pacjenta. Do badania włączono pacjentów, którzy mieli ≥18 lat, potwierdzoną diagnozę pierwotnego TBI, dostępną liczbę neutrofili przy przyjęciu, wymagane do analizy bazowe zmienne kliniczne i laboratoryjne oraz potwierdzony status przeżycia po 14 dniach od przyjęcia. Zmiennymi bazowymi wymaganymi do włączenia do badania były: wiek, płeć, kategoria diagnostyczna, status operacyjny (kraniotomia lub dekompresyjna kraniektomia), czas protrombinowy (PT), czas aktywowanej częściowej tromboplastyny (APTT), fibrynogen (FIB), liczba płytek krwi (PLT), albuminy (ALB), cholesterol całkowity (TC), cholesterol lipoprotein o wysokiej gęstości (HDL), cholesterol lipoprotein o niskiej gęstości (LDL) oraz poziom glukozy we krwi na czczo (FBG).
Pacjentów wykluczono, jeśli mieli mniej niż 18 lat, byli hospitalizowani wielokrotnie, cierpieli na nietraumatyczny uraz mózgu lub mieli niejasną diagnozę, brakowało danych o liczbie neutrofili w momencie przyjęcia, brakowało kluczowych klinicznych lub laboratoryjnych zmiennych wyjściowych lub brakowało informacji o przeżywalności po 14 dniach. Ekstremalne wartości liczby neutrofili zidentyfikowano poprzez sprawdzenie zakresów, a następnie zweryfikowano z oryginalną dokumentacją laboratoryjną. Kryteria sprawdzania zakresu dla ekstremalnych liczb neutrofili określono jako wartości poniżej 0,5 × 109/L lub powyżej 40,0 × 109/L. Wartości wykluczano tylko wtedy, gdy potwierdzono, że wynikały one z błędów wprowadzania danych, błędów w raportach laboratoryjnych lub klinicznie nieprawdopodobnych błędów pomiarowych. Klinicznie nieprawdopodobne błędy pomiarowe zdefiniowano jako liczby neutrofili, które były albo (1) niezgodne z jednoczesną liczbą białych krwinek i rozmazem u pacjenta (np. liczba neutrofili przekraczająca całkowitą liczbę białych krwinek) lub (2) zapisane jako wartości nienumeryczne lub znaki zastępcze w surowym eksporcie z elektronicznej dokumentacji medycznej. Ostatecznie do końcowej kohorty analitycznej włączono 1361 pacjentów.
Pacjentów zaklasyfikowano zgodnie z główną kategorią diagnostyczną TBI odnotowaną podczas pierwszej hospitalizacji. Kategorie diagnostyczne obejmowały urazowy krwotok wewnątrzczaszkowy, uraz pnia mózgu oraz proste złamanie czaszki. W przypadku odnotowania więcej niż jednej diagnozy związanej z TBI, kategorię diagnostyczną przypisano zgodnie z główną diagnozą przy wypisie. Główną diagnozę przy wypisie zdefiniowano jako pierwszą wymienioną diagnozę w epikryzie, przypisaną przez prowadzącego neurochirurga na podstawie stanu klinicznie najistotniejszego, wymagającego najdłuższego czasu hospitalizacji lub największego wykorzystania zasobów opieki zdrowotnej podczas przyjęcia.
Włączeni pacjenci zostali następnie podzieleni na grupy zgodnie z tercylami liczby neutrofili w momencie przyjęcia. Niski tercyl obejmował wartości 1,650–7,057 × 109/L, średni tercyl 7,057–10,947 × 109/L, a wysoki tercyl 10,947–32,40 × 109/L. Punkty odcięcia tercyli określono poprzez posortowanie liczby neutrofili przy przyjęciu w kolejności rosnącej i podzielenie kohorty na trzy grupy o zbliżonej liczebności (odpowiednio n = 454, 453 i 454). W całym tekstem opracowania użyto terminu „tercyle”, ponieważ kohorta została podzielona na trzy grupy.
Gromadzenie danych dotyczących ekspozycji, wyników i danych klinicznych
Główną ekspozycją była bezwzględna liczba neutrofili zmierzona przy przyjęciu. Laboratoryjne wartości przy przyjęciu zdefiniowano jako pierwsze dostępne wyniki laboratoryjne uzyskane podczas wstępnej oceny na oddziale ratunkowym lub oceny przy przyjęciu. Próbki krwi pobierano zazwyczaj bezpośrednio po przyjęciu, a wyniki testów raportowano zazwyczaj w ciągu 2 h. Dopuszczalne okno czasowe dla zdefiniowania laboratoryjnych wartości przy przyjęciu wynosiło 6 h od przybycia na oddział ratunkowy lub przyjęcia do szpitala. Jeśli w tym początkowym okresie oceny dostępny był więcej niż jeden wynik, do analizy wykorzystano wynik najwcześniejszy. Liczbę neutrofili zapisano w jednostkach ×109/L i analizowano zarówno jako zmienną ciągłą, jak i w podziale na grupy tertylowe.
Głównym punktem końcowym była 14-dniowa śmiertelność z przyczyn ogólnych po indeksowej hospitalizacji. Zgon w ciągu 14 dni zakodowano jako 1, a przeżycie powyżej 14 dni zakodowano jako 0. Śmiertelność w czternastym dniu określono na podstawie dokumentacji zgonów pacjentów stacjonarnych, dokumentacji wypisowej oraz szpitalnego systemu monitorowania stanu zdrowia po wypisie. W przypadku rozbieżności między źródłami danych o śmiertelności (np. gdy pacjent został odnotowany jako żywy w momencie wypisu, ale jako zmarły w systemie monitorowania w ciągu 14 dni), rekord z systemu monitorowania uznano za standard odniesienia, ponieważ obejmował on śmiertelność po wypisie ze szpitala. Pacjenci, których status przeżycia w 14. dniu nie udało się potwierdzić, zostali wykluczeni z analizy końcowej.
Z elektronicznego systemu dokumentacji medycznej wyodrębniono zmienne demograficzne i kliniczne. Zmienne te obejmowały wiek, płeć, kategorię diagnostyczną, status operacyjny oraz status tracheostomii. Wynik w skali Glasgow (GCS), reakcja źrenic, klasyfikacja Marshalla w tomografii komputerowej (CT) oraz stopień urazów pozaczaszkowych nie zostały uwzględnione w analizie, ponieważ zmienne te nie były konsekwentnie rejestrowane w elektronicznym systemie dokumentacji medycznej w okresie badania (2017–2023), szczególnie w przypadku pacjentów z łagodnym TBI lub tych, którzy nie zostali przyjęci na oddział intensywnej terapii. W związku z tym zmienne te nie mogły zostać wiarygodnie wyodrębnione dla większości kohorty. Operację zdefiniowano jako kraniotomię w celu usunięcia krwiaka lub dekompresyjną kraniektomię wykonaną podczas indeksowej hospitalizacji. Tracheostomię zakodowano jako „tak” lub „nie” w zależności od tego, czy procedura została wykonana podczas tej samej hospitalizacji.
Zmienne laboratoryjne przy przyjęciu pobrano z pierwszej dostępnej oceny laboratoryjnej podczas wstępnej oceny ratunkowej lub przyjęcia. Zmienne koagulacyjne i hematologiczne obejmowały PT, APTT, międzynarodowy wskaźnik INR (INR), FIB, D-dimer, PLT, hemoglobinę (Hb) oraz liczbę neutrofili. Zmienne biochemiczne obejmowały ALB, TC, triglicerydy (TGs), HDL, LDL, apolipoproteinę A-I (Apo A-I), apolipoproteinę B (Apo B), cystatynę C (CYS), aminotransferazę alaninową (ALT), aminotransferazę asparaginianową (AST), kreatyninę (Cr), azot krwi z mocznikiem (BUN), kwas moczowy (UA) oraz FBG.
Wszystkie próbki krwi zostały przebadane w szpitalnym laboratorium klinicznym zgodnie ze standardowymi procedurami operacyjnymi. Oznaczenia morfologii krwi wykonano przy użyciu automatycznego analizatora hematologicznego, badania krzepliwości przy użyciu automatycznego analizatora krzepliwości, a badania biochemiczne przy użyciu automatycznego analizatora biochemicznego. Wszystkie testy zostały przeprowadzone przez przeszkolonych techników laboratoryjnych, przy rutynowej wewnętrznej kontroli jakości prowadzonej w ramach programu kontroli jakości laboratorium klinicznego.
Wstępne przetwarzanie danych i porównanie z poziomem bazowym
Przed analizą statystyczną wyodrębniony zestaw danych sprawdzono pod kątem powtórnych przyjęć, brakujących wartości, zdublowanych rekordów, niespójnych jednostek oraz wartości klinicznie nieplausible. Powtórne przyjęcia usunięto, zachowując jedynie pierwszą kwalifikującą się hospitalizację. Ujednolicono jednostki laboratoryjne w całym okresie badania. W szczególności liczbę neutrofili i liczbę płytek krwi ustandaryzowano do ×109/L, hemoglobinę do g/L, PT i APTT do sekund, INR pozostawiono jako miara bezwymiarowa, FIB do g/L, D-dimer do mg/L, ALB, Apo A-I i Apo B do g/L, TC, TGs, HDL, LDL oraz glukozę w surowicy na czczo (FBG) do mmol/L, CYS do mg/L, ALT i AST do U/L, kreatyninę (Cr) i kwas moczowy (UA) do µmol/L oraz BUN do mmol/L. Wartości nieplausible zweryfikowano z oryginalną dokumentacją medyczną lub raportami laboratoryjnymi i potwierdzono w oparciu o rekordy źródłowe. Wartości klinicznie nieplausible zdefiniowano jako wartości wykraczające poza biologicznie dopuszczalne zakresy dla odpowiadającego im parametru laboratoryjnego (np. PT <5 s lub >120 s, FIB <0,1 g/L lub >20 g/L, FBG <1,0 mmol/L lub >50 mmol/L lub ALT/AST >10 0 U/L). Dodatkowe wartości nieplausible obejmowały niespójności logiczne, takie jak liczba neutrofili przekraczająca całkowitą liczbę białych krwinek.
Główna analiza została przeprowadzona z wykorzystaniem kompletnym przypadków. Pacjenci zostali wykluczeni z końcowej kohorty analitycznej, jeśli brakowało danych dotyczących liczby neutrofili przy przyjęciu, statusu śmiertelności po 14 dniach lub współzmiennych wymaganych do głównej analizy skorygowanej. Współzmienne wymagane dla w pełni skorygowanego modelu (Model 4) to wiek, kategoria diagnostyczna, status operacyjny, APTT, FIB, PLT, LDL, HDL oraz INR. Po wyborze kohorty wszystkie zmienne uwzględnione w sekwencyjnych modelach regresji były dostępne dla 1 361 pacjentów włączonych do analizy końcowej.
Zmienne kategoryczne zakodowano w następujący sposób: płeć jako męska lub żeńska, operacja i tracheostomia jako tak lub nie, kategoria diagnostyczna jako urazowe krwawienie wewnątrzczaszkowe, urazowe uszkodzenie pnia mózgu lub proste złamanie czaszki oraz śmiertelność w ciągu 14 dni jako zgon lub przeżycie.
Zmienne ciągłe zostały zestawione w głównych tabelach opisowych jako średnia ± odchylenie standardowe. Zmienne kategoryczne przedstawiono jako liczebności i wartości procentowe. Porównania między grupami przeprowadzono za pomocą analizy wariancji lub nieparametrycznych testów rangowych, w zależności od zastosowania. Jednoczynnikowa analiza wariancji (ANOVA) została zastosowana dla zmiennych ciągłych spełniających założenia normalności rozkładu i jednorodności wariancji, ocenianych za pomocą testu Shapiro-Wilka (P > 0,05) oraz test Levene’a (P > odpowiednio. W przypadku zmiennych ciągłych, które nie spełniały tych założeń, zastosowano test Kruskala–Wallisa. Dla zmiennych kategorycznych, w zależności od oczekiwanej liczebności komórek, zastosowano test chi-kwadrat Pearsona lub dokładny test Fishera. Dokładny test Fishera stosowano w sytuacjach, gdy oczekiwana liczebność którejkolwiek komórki w tabeli kontyngencji była mniejsza niż 5. W tabeli danych wyjściowych, PWartość p została wykorzystana jako główna podstawa do porównań między grupami. Pwartość p* została zachowana jedynie jako uzupełniający punkt odniesienia dla odpowiadającego jej nieparametrycznego testu alternatywnego lub testu dla małych prób, jeżeli był on stosowany; istotność statystyczną w sekcji Wyniki oceniano zgodnie z Pwartość p oraz te P wartości nie zostały skorygowane pod kątem wielokrotnych porównań.
Aby ułatwić interpretację różnic wyjściowych, wygenerowano opisowe wykresy wielopanelowe dla wybranych zmiennych klinicznych i laboratoryjnych, w tym wieku, poziomu D-dimerów, glukozy w krwi na czczo, fibrynogenu, czasu protrombinowego, składu diagnostycznego, interwencji chirurgicznej, tracheostomii oraz śmiertelności w 14. dobie. Zmienne te wybrano do wizualizacji, ponieważ reprezentowały one istotne klinicznie cechy i wykazały znaczące różnice pomiędzy tercylami neutrofili w analizach wyjściowych. Wykresy te posłużyły do przedstawienia profili klinicznych i laboratoryjnych pacjentów w grupach niskiego, średniego i wysokiego tercyla neutrofili.
Analiza statystyczna i wizualizacja
W pierwszej kolejności zastosowano jednowymiarową regresję logistyczną w celu oceny powiązań między zmiennymi kandydującymi a śmiertelnością w 14. dobie. Przedstawiono ilorazy szans (OR), 95% przedziały ufności (CI) oraz wartości P. Zmienne rozważano pod kątem korekty wielowymiarowej, jeśli były istotne klinicznie, wykazały związek ze śmiertelnością w 14. dobie w analizie jednowymiarowej lub zmieniały OR dla liczby neutrofili o ≥10% po dodaniu do modelu lub usunięciu z niego. Istotność kliniczną określono a priori na podstawie uznanych czynników prognostycznych dla urazowego uszkodzenia mózgu oraz znanych mechanizmów patofizjologicznych opisanych w literaturze, w tym wieku, parametrów krzepnięcia, glukozy oraz zmiennych związanych z lipidami. Próg identyfikacji zmiennych związanych ze śmiertelnością w 14. dobie w analizie jednowymiarowej ustalono na poziomie P < 0,10, aby uniknąć wykluczenia potencjalnie ważnych czynników zakłócających z modeli wielowymiarowych.
Związek między liczbą neutrofili w momencie przyjęcia a śmiertelnością w ciągu 14 dni zbadano następnie za pomocą wieloczynnikowej regresji logistycznej. Liczbę neutrofili analizowano jako zmienną ciągłą w odniesieniu do każdego wzrostu o 1 ×109/L. Opracowano cztery sekwencyjne modele. Model 1 nie był korygowany. Model 2 skorygowano o wiek, rozpoznanie i operację. Model 3 skorygowano o wiek, rozpoznanie, operację, APTT, FIB, PLT oraz LDL. Model 4 skorygowano o wiek, rozpoznanie, operację, APTT, FIB, PLT, LDL, HDL oraz INR. Do oceny potencjalnej współliniowości zmiennych towarzyszących przed dopasowaniem modeli wieloczynnikowych zastosowano czynniki inflacji wariancji (VIF). Jako próg identyfikacji potencjalnie problematycznej współliniowości przyjęto wartość VIF równą 5.
Aby zbadać, czy zależność ta była nieliniowa, zastosowano uogólnione modele addytywne w celu oceny relacji ekspozycja-odpowiedź pomiędzy liczbą neutrofili a śmiertelnością w ciągu 14 dni. Jako funkcję wygładzającą wykorzystano sześcienne splajny regresyjne z wymiarem bazy (k) równym 4. Parametr wygładzania oszacowano przy użyciu ograniczonej maksymalnej wiarygodności (REML). Wygenerowano skorygowaną krzywą wygładzania, aby zwizualizować zależność między liczbą neutrofili w momencie przyjęcia a przewidywanym prawdopodobieństwem śmiertelności w ciągu 14 dni, po skorygowaniu o współzmienne uwzględnione w pełnym modelu dostosowania. Gdy krzywa sugerowała wzorzec nieliniowy, do oszacowania punktu przegięcia zastosowano dwuczęściową regresję liniową. Punkt przegięcia (K) zidentyfikowano za pomocą algorytmu rekurencyjnego, który oceniał potencjalne punkty przełamania w zakresie od 5. do 95. percentyla liczby neutrofili. Punkt przełamania związany z najlepszym dopasowaniem modelu, określone na podstawie kryterium informacyjnego Akaike, wybrano jako ostateczny punkt przegięcia. Jednowątkowy model liniowy oraz model dwuczęściowy porównano za pomocą testu ilorazu wiarygodności. Efekty progowe raportowano poniżej i powyżej oszacowanego punktu przegięcia.
Efekty progowe zwizualizowano za pomocą wykresu leśnego (forest plot), przedstawiającego szacowane wartości OR poniżej i powyżej punktu przegięcia K specyficznego dla danego modelu. Wartości K w Modelach 1–4 wynosiły odpowiednio 18,54, 3,82, 3,86 oraz 3,81 ×109/L. W przypadku zmiennych kategorycznych z bardzo małą liczbą lub brakiem zdarzeń końcowych w danej podgrupie, odpowiadające im wartości OR uznano za potencjalnie niestabilne i interpretowano z ostrożnością. Szacunki OR klasyfikowano jako niestabilne, jeśli podgrupa zawierała mniej niż pięć zdarzeń końcowych lub jeśli szerokość 95% przedziału ufności przekraczała 10,0.
Rysunki zostały przygotowane w celu przedstawienia procesu przesiewowego pacjentów, wyjściowych profili klinicznych i laboratoryjnych w podziale na tercyle neutrofili, wskaźników śmiertelności w ciągu 14 dni w podziale na tercyle, szacunków regresji z modeli sekwencyjnych, efektów progowych oraz skorygowanych krzywych wygładzających. Analizy statystyczne przeprowadzono z użyciem programu R w wersji 4.5.2 oraz platformy online EmpowerStats. Rysunki zostały wygenerowane przy użyciu pakietów R: ggplot2, patchwork, pdftools, png oraz scales. Dwustronną wartość P <0,05 uznano za istotną statystycznie.