Niniejsze badanie retrospektywne zostało przeprowadzone zgodnie z wymogami instytucjonalnymi oraz Deklaracją Helsińską i zostało zatwierdzone przez Komisję Etyczną Szpitala Nanjing Tongren powiązanego z Wydziałem Medycznym Uniwersytetu Południowo-Wschodniego (numer zatwierdzenia 2024-03-039-K001). Ze względu na retrospektywny charakter badania oraz wykorzystanie zanonimizowanych danych, zrezygnowano z wymogu uzyskania indywidualnej świadomej zgody.
Dane kliniczne
W ramach grupy badawczej uwzględniono łącznie 172 dzieci z MPP leczonych między kwietniem 2018 r. a kwietniem 2021 r., a także 40 zdrowych dzieci jako grupę kontrolną. Grupa badawcza składała się ze 103 chłopców i 69 dziewcząt, średnia wieku wynosiła 7,15 ± 2,97 lat; grupa kontrolna składała się z 24 chłopców i 16 dziewcząt, średnia wieku wynosiła 7,02 ± 2,31 lat. Wiek i płeć wyjściowa nie różniły się istotnie między grupami (P > 0,05). Zgodnie z opisanymi kryteriami10, MPP klasyfikowano jako ciężkie, gdy wystąpiły co najmniej dwa z następujących objawów: wyraźna duszność, częstość akcji serca >120 uderzeń/min, gorączka trwająca >10 dni, zajęcie wielu płatów płuc lub rozległe nacieki płucne. Pozostałe przypadki zakwalifikowano jako MPP o przebiegu łagodnym do umiarkowanego.
Kryteria włączenia
Wszystkie dzieci z MPP spełniały kliniczne kryteria diagnostyczne dla MPP11. Rozpoznanie potwierdzono mianem przeciwciał przeciw MP >1:160, pozytywnym wynikiem testu na DNA MP lub czterokrotnym lub większym wzrostem poziomu swoistych dla MP IgG w surowicy między fazą ostrą a fazą rekonwalescencji. Objawy kliniczne obejmowały kaszel, duszność, tachykardię, sinicę lub dyspnoe, a obrazowanie klatki piersiowej wykazało jednostronne lub obustronne zmiany w płucach. Do badania zakwalifikowano uczestników w wieku od 8 miesięcy do 10 lat.
Dzieci wykluczono z badania, jeśli stwierdzono u nich wirusowe, bakteryjne lub grzybicze zapalenie płuc; ciężką niewydolność głównych narządów, w tym serca, mózgu, wątroby lub nerek; inną chorobę płuc; nowotwór złośliwy; wrodzoną dysfunkcję odporności lub niewystarczające informacje kliniczne do wymaganych analiz.
Dzieci w grupie kontrolnej były w wieku od 8 miesięcy do 10 lat, były zdrowe i miały prawidłowe wyniki badania fizykalnego. Grupa kontrolna została wykorzystana wyłącznie do porównań biomarkerów i nie brała udziału w analizach prognostycznych ani regresyjnych.
Gromadzenie danych
Przeanalizowano szpitalny system informacyjny pod kątem wieku, płci, czasu trwania gorączki, rodzaju i lokalizacji zmian, stanu funkcji odpornościowej, poziomu CRP, powikłań pozapłucnych, informacji o leczeniu oraz wyników obserwacji. Dysfunkcję odpornościową określono na podstawie współczesnej dokumentacji medycznej przed oceną wyników, w oparciu o wartości immunoglobulin, które znajdowały się poniżej przedziałów referencyjnych laboratorium szpitalnego dla danej grupy wiekowej. Dane zostały początkowo wyekstrahowane przez Bolinga Hana, a następnie zweryfikowane z oryginalną dokumentacją medyczną przez Rui Wang. Wszelkie rozbieżności rozstrzygano w drodze dyskusji między obwojem badaczy. Przed analizą zweryfikowano kompletność danych. Pełne definicje zmiennych, reguły kodowania oraz zanonimizowany zestaw danych kohorty znajdują się w Pliku uzupełniającym 1.
Pomiar wskaźników funkcji odpornościowej i markerów stanu zapalnego
W ciągu 24 godzin od zapisu każdego uczestnika pobrano 4 mL krwi żylnej na czczo. Próbki wirowano przy 1,500 × g przez 10 min w temperaturze 4°C, a otrzymane surowice przechowywano zgodnie z instytucjonalnym protokołem laboratoryjnym do czasu analizy. Stężenia IgA, IgG i IgM w surowicy zmierzono metodą immunoturbidymetrii, a poziom CRP oznaczono za pomocą enzymatycznego testu immunoabsorpcyjnego (ELISA). Próbki surowicy przechowywano w temperaturze –80 °C nie dłużej niż przez 3 miesiące przed analizą, unikając wielokrotnych cykli zamrażania i rozmrażania.
Protokół leczenia
Leczenie było dostosowane indywidualnie do stopnia ciężkości stanu klinicznego oraz obowiązującej praktyki instytucjonalnej. Standardowe postępowanie obejmowało terapię przeciwdrobnoustrojową oraz opiekę wspomagającą, w tym tlenoterapię, inhalacje nebulizowane, podawanie płynów i utrzymanie równowagi wodno-elektrolitowej, leczenie przeciwgorączkowe oraz inne środki objawowe, zgodnie ze wskazaniami klinicznymi. Dzieci z ciężką postacią MPP leczono azytromycyną, a metyloprednizolon stosowano jako pomocniczą terapię przeciwzapalną we wskazaniach klinicznych. Chociaż wszyscy pacjenci byli leczeni zgodnie z tymi samymi zasadami instytucjonalnymi, konkretne środki wspomagające były dobierane indywidualnie, a nie ściśle standaryzowane. Ze względu na retrospektywny charakter badania, dokładny czas rozpoczęcia leczenia oraz szczegóły schematów terapeutycznych nie były jednolicie dostępne dla wszystkich uczestników.
Obserwacja prognostyczna
Dzieci z ciężkim przebiegiem MPP były monitorowane co miesiąc podczas wizyt ambulatoryjnych oraz poprzez kontakt telefoniczny do kwietnia 2022 roku. Rokowanie uznano za dobre, gdy temperatura ciała wróciła do normy, objawy oddechowe uległy poprawie, a zmiany płucne zostały wchłonięte lub ustąpiły. Rokowanie złe zdefiniowano jako wynik złożony, obejmujący utrzymującą się gorączkę, nasilenie kaszlu lub duszności, brak poprawy w obrazowaniu radiologicznym lub pogorszenie obrazu radiologicznego, przeniesienie do innego szpitala lub zgon12. Składniki te mogły się nakładać i nie były analizowane jako wzajemnie wykluczające się wyniki. Wszystkie 114 dzieci z ciężkim MPP ukończyło okres obserwacji; nie odnotowano utraty kontaktu z żadnym z uczestników.
Analiza statystyczna
Analizy przeprowadzono przy użyciu programu IBM SPSS Statistics wersja 22.0. Zmienne ciągłe przedstawiono jako średnia ± odchylenie standardowe. Do porównań dwóch grup zastosowano test t dla prób niezależnych, a do porównań więcej niż dwóch grup — jednoczynnikową analizę wariancji. Zmienne kategoryczne przedstawiono jako n (%) i porównano za pomocą testu chi-kwadrat. Krzywe ROC wykorzystano do indywidualnej oceny IgA, IgG, IgM oraz CRP, raportując wartości AUC, 95% przedziały ufności, wartości odcięcia, czułość i swoistość. Czterech biomarkerów nie łączono w wieloczynnikowy model predykcyjny. Zmienne związane z rokowaniem w analizach jednowymiarowych wprowadzono do eksploracyjnej krokowej binarnej regresji logistycznej. Wyniki regresji przedstawiono jako współczynniki beta, błędy standardowe, statystyki Walda, ilorazy szans (ORs), 95% przedziały ufności oraz wartości P. Ze względu na ograniczoną liczbę zdarzeń o złym rokowaniu, regresję zinterpretowano jako eksploracyjną analizę asocjacji, a nie jako zwalidowany model predykcyjny. Za statystycznie istotne uznano dwustronne P < 0,05. Normalność rozkładu zmiennych ciągłych oceniono testem Shapiro-Wilka, a jednorodność wariancji oceniono testem Levene'a. W zmiennych uwzględnionych w analizach nie wystąpiły braki danych; w związku z tym nie przeprowadzono imputacji i zastosowano analizy przypadków pełnych.