Zanieczyszczenie powietrza a pneumokoki. PM2,5 wiązało się z większym ryzykiem ciężkich zakażeń

Analiza 59 017 przypadków inwazyjnej choroby pneumokokowej w Republice Południowej Afryki wykazała dodatni związek między ekspozycją na pył zawieszony PM_2,5 a ryzykiem zachorowania. Wielkość i czas wystąpienia tego związku nie były jednak jednakowe dla wszystkich pacjentów i wszystkich pneumokoków. Istotne znaczenie miały wiek, postać kliniczna zakażenia, serotyp oraz linia genomowa Streptococcus pneumoniae. Badanie opublikowane w Nature Microbiology pokazuje, że ocena wpływu zanieczyszczenia powietrza na zakażenia może wymagać uwzględnienia także struktury populacji bakteryjnej.

Jak przeprowadzono analizę 59 017 przypadków IChP?

Badanie przeprowadzili Sophie Belman, Cebile Lekhuleni, Jackie Kleynhans i współpracownicy z zespołów związanych m.in. z Barcelona Supercomputing Center (BSC), Wellcome Sanger Institute, National Institute for Communicable Diseases of the National Health Laboratory Service w Johannesburgu, University of the Witwatersrand oraz innymi ośrodkami. Współautorami seniorskimi były Anne von Gottberg i Rachel Lowe.

Badacze analizowali inwazyjną chorobę pneumokokową (IChP; ang. invasive pneumococcal disease, IPD), czyli zakażenia, w których Streptococcus pneumoniae przedostaje się do miejsc fizjologicznie jałowych. W analizowanym zbiorze obejmowało to przede wszystkim bakteriemię, zapalenie opon mózgowo-rdzeniowych oraz inne postacie zakażenia inwazyjnego, w tym ropniaka i zapalenie osierdzia.

Podstawę epidemiologiczną stanowił południowoafrykański system nadzoru GERMS-SA. Do zbioru włączono 59 017 rozpoznanych przypadków IChP z 531 szpitali, zebranych w latach 2005-2023. Dane pochodziły ze wszystkich dziewięciu prowincji Republiki Południowej Afryki. W analizach dotyczących struktury genomowej pneumokoków wykorzystano 4 350 izolatów poddanych sekwencjonowaniu genomowemu z okresu do 2019 roku.

Wśród wszystkich przypadków 62% dotyczyło osób w wieku 15-64 lat, 23% dzieci poniżej 6. roku życia, a 4,8% osób w wieku powyżej 64 lat. Mediana wieku wynosiła 32 lata. Badanie obejmowało zatem bardzo zróżnicowaną wiekowo populację i nie było ograniczone ani do pediatrii, ani do dorosłych.

Autorzy łączyli dane epidemiologiczne z informacjami o warunkach meteorologicznych i jakości powietrza. Analizowano m.in. temperaturę, wilgotność względną i bezwzględną, opady, prędkość wiatru oraz stężenia PM_2,5, PM₁₀ i SO₂. Szczególną uwagę poświęcono drobnym cząstkom pyłu zawieszonego PM_2,5.

Modele uwzględniały przestrzenne i czasowe zróżnicowanie zachorowań oraz potencjalne czynniki zakłócające, w tym gęstość zaludnienia, okres wdrażania szczepień przeciw pneumokokom, temperaturę i wilgotność. Zastosowano modele pozwalające analizować nie tylko wielkość związku z ekspozycją, lecz także jego opóźnienie – od tygodnia ekspozycji do kolejnych tygodni.

PM2,5 wiązał się ze wzrostem ryzyka inwazyjnej choroby pneumokokowej

W modelach dotyczących jakości powietrza punkt odniesienia dla PM_2,5 ustalono na 40 µg/m³, odpowiadający przyjętemu w analizie południowoafrykańskiemu progowi jakości powietrza. Przy średnim tygodniowym stężeniu PM_2,5 wynoszącym 50 µg/m³ skumulowane w ciągu ośmiu tygodni ryzyko IChP było o około 3,4% większe – RR 1,03; 95% CI 1,02-1,05.

Znacznie większą różnicę obserwowano przy bardzo wysokiej ekspozycji. Dla średniego tygodniowego stężenia PM_2,5 wynoszącego 100 µg/m³ skumulowane ośmiotygodniowe RR wyniosło 1,18 (95% CI 1,12-1,23), czyli odpowiadało około 18-procentowemu wzrostowi względnego ryzyka względem przyjętego punktu odniesienia.

Związek przy wysokiej ekspozycji był widoczny już w tygodniu ekspozycji – dla 100 µg/m³ RR w tygodniu 0 wyniosło 1,02 (95% CI 1,01-1,04) – a podwyższone ryzyko utrzymywało się w modelach przez około sześć tygodni.

Nie tylko PM_2,5 wykazywał dodatnią zależność z IChP. Przy średnim tygodniowym stężeniu PM₁₀ równym 100 µg/m³ skumulowane ośmiotygodniowe RR wyniosło około 1,06 (95% CI 1,04-1,07). Dla SO₂ przy stężeniu 150 µg/m³ autorzy obserwowali wzrost ryzyka pojawiający się z większym opóźnieniem – od około trzeciego tygodnia po ekspozycji.

Temperatura i wilgotność również modyfikowały ryzyko

Zależności środowiskowe nie ograniczały się do zanieczyszczenia powietrza. W ośmiotygodniowym okresie obserwacji umiarkowana wilgotność względna, około 33-49%, oraz niskie temperatury minimalne, około 4-10°C, wiązały się odpowiednio z około 5- i 4-procentowym wzrostem ryzyka IChP. Jednocześnie wysokie temperatury maksymalne, w zakresie około 27-38°C, były związane z krótkoterminowym wzrostem ryzyka sięgającym około 10%.

Wyniki wskazują więc na złożoną, nieliniową zależność między środowiskiem a zachorowaniami. Nie można jej sprowadzić do prostego stwierdzenia, że wyłącznie zimno albo wyłącznie zanieczyszczenie zwiększa częstość zakażeń.

Wiek i postać kliniczna wpływały na wielkość związku z PM2,5

Jednym z istotnych wyników było zróżnicowanie efektu PM_2,5 w zależności od wieku. Największe skumulowane ryzyko obserwowano u osób w wieku co najmniej 65 lat. Przy średnim tygodniowym stężeniu 50 µg/m³ PM_2,5 RR w tej grupie wyniosło 1,16 (95% CI 1,10-1,22).

Inny wzorzec czasowy stwierdzono u osób w wieku 15-64 lat. W tej grupie wpływ wysokiego stężenia PM_2,5 był bardziej natychmiastowy i wzrost ryzyka pojawiał się już w tym samym tygodniu. Autorzy zwracają uwagę, że interpretacja tej obserwacji wymaga uwzględnienia specyfiki badanej populacji, w tym obciążenia zakażeniem HIV w Republice Południowej Afryki.

U dzieci w wieku do 5 lat związek również był widoczny. Przy 50 µg/m³ PM_2,5 skumulowane ośmiotygodniowe RR wyniosło 1,03 (95% CI 1,01-1,05). Wynik wskazuje, że zanieczyszczenie powietrza nie jest potencjalnie istotnym czynnikiem wyłącznie dla starszych osób.

Zapalenie opon mózgowo-rdzeniowych i inne postacie zakażenia

Autorzy analizowali również różne manifestacje IChP. Przy tygodniowym PM_2,5 równym 50 µg/m³ RR dla zapalenia opon mózgowo-rdzeniowych wyniosło 1,05 (95% CI 1,03-1,06), podczas gdy dla bakteriemii 1,01 (95% CI 1,00-1,03).

Jeszcze większą wartość uzyskano dla grupy innych zakażeń inwazyjnych, do których należały m.in. ropniak i zapalenie osierdzia – RR 1,12 (95% CI 1,08-1,16). Dla przypadków IChP z opornością na penicylinę przy tym samym poziomie PM_2,5 RR wyniosło 1,09 (95% CI 1,03-1,15).

Dane te nie oznaczają, że pył zawieszony wywołuje określoną postać zakażenia. Pokazują natomiast, że obserwowana zależność ekspozycja-choroba nie była jednakowa dla wszystkich analizowanych fenotypów klinicznych.

Serotyp pneumokoka zmieniał czas odpowiedzi na zanieczyszczenie powietrza

Streptococcus pneumoniae nie jest jednorodnym biologicznie patogenem. Znanych jest ponad 100 serotypów pneumokoka, a różnice genomowe pozwalają dodatkowo wyróżniać liczne linie określane jako Global Pneumococcal Sequence Clusters (GPSC).

W badaniu serotyp okazał się jednym z najważniejszych czynników modyfikujących zależność między PM_2,5 a IChP. Dla serotypów 4, 8 i 23F obserwowano wzrost ryzyka pojawiający się już bezpośrednio po wysokiej ekspozycji na PM_2,5. W przypadku innych serotypów najwyższe ryzyko pojawiało się dopiero po upływie kilku tygodni.

Analiza wskazała ponadto wyższe punktowe oszacowania skumulowanego ryzyka przy 50 µg/m³ dla choroby wywołanej przez serotypy 14, 19A i 8 w porównaniu z modelem obejmującym wszystkie przypadki. Interpretacja tego elementu wymaga jednak ostrożności. Dla części porównań przedziały ufności obejmowały wartość 1 – dla serotypu 19A RR ratio wyniosło 1,04 (95% CI 0,99-1,09), a dla serotypu 8 1,05 (95% CI 0,99-1,11). Wyniki te wskazują zatem na zróżnicowanie kierunku i wielkości efektu, ale nie wszystkie różnice między serotypami mają taką samą precyzję statystyczną.

GPSC21 i serotyp 19F

Szczególnie interesujące wyniki uzyskano dla linii genomowej GPSC21, która w analizowanym materiale była niemal całkowicie związana z serotypem 19F. Gdy GPSC21 występowała w danym miejscu i czasie z dużą częstością, wysokie stężenie PM_2,5 wiązało się ze wzrostem ryzyka IChP już w tym samym tygodniu. Przy małej częstości GPSC21 wzrost ryzyka pojawiał się dopiero około pięciu tygodni po ekspozycji.

Autorzy podkreślają, że mechanizm tego zjawiska pozostaje niewyjaśniony. Wynik sugeruje jednak, że przy analizie związku zanieczyszczenia powietrza z zakażeniami znaczenie może mieć nie tylko obecność konkretnego gatunku bakterii, ale również to, jakie linie genomowe danego patogenu aktualnie krążą w populacji.

Jak zanieczyszczenie powietrza może wpływać na pneumokoki i gospodarza?

Badanie miało charakter epidemiologiczny i nie było eksperymentem mechanistycznym, dlatego nie pozwala określić bezpośredniej przyczyny obserwowanych zależności. Autorzy omawiają jednak kilka biologicznie wiarygodnych mechanizmów opisanych wcześniej w literaturze.

Po stronie gospodarza zanieczyszczenie pyłowe może zwiększać podatność nabłonka dróg oddechowych na kolonizację, zaburzać klirens śluzowo-rzęskowy i uszkadzać barierę nabłonkową. Po stronie bakterii wcześniejsze badania wskazywały m.in. na możliwość zmian adhezji do komórek nabłonka, tworzenia biofilmu, obciążenia bakteryjnego oraz zdolności do transformacji genetycznej.

Autorzy rozważają także możliwość istnienia mechanizmów charakterystycznych dla poszczególnych serotypów. Jako przykład wskazują serotyp 4, w którym określone cechy wirulencji mogłyby potencjalnie oddziaływać z uszkodzeniem nabłonka związanym z PM_2,5. Jest to jednak hipoteza biologiczna wymagająca bezpośredniego potwierdzenia eksperymentalnego, a nie mechanizm wykazany w omawianym badaniu.

W szerszej perspektywie wyniki pokazują, że relacja między środowiskiem a zakażeniem może być wynikiem interakcji co najmniej trzech komponentów: podatności gospodarza, poziomu ekspozycji środowiskowej oraz właściwości biologicznych aktualnie krążącej populacji patogenu.

Ograniczenia badania są istotne dla interpretacji wyników

Mimo dużej liczebności i 19-letniego okresu obserwacji publikacja ma kilka ważnych ograniczeń. Przede wszystkim jest to badanie obserwacyjne oparte na modelowaniu przestrzenno-czasowym, dlatego wykazane związki nie dowodzą, że zanieczyszczenie powietrza bezpośrednio powodowało konkretne przypadki IChP.

Co wyniki mogą oznaczać dla nadzoru epidemiologicznego?

Najważniejszym elementem pracy jest pokazanie, że wpływ czynników środowiskowych na zakażenia pneumokokowe może być modyfikowany przez strukturę populacji samego patogenu. Dwa miasta o podobnym poziomie PM_2,5 nie muszą zatem wykazywać identycznego wzorca czasowego IChP, jeżeli różni je rozpowszechnienie serotypów i linii genomowych S. pneumoniae.

Autorzy proponują, aby w przyszłości łączyć dane o jakości powietrza i meteorologii z genomicznym nadzorem nad pneumokokami. Taki system mógłby potencjalnie ułatwiać prognozowanie okresów zwiększonego obciążenia IChP, planowanie zapotrzebowania na świadczenia oraz dokładniejsze modelowanie strategii szczepień.

Nie oznacza to jednak, że na podstawie stężenia pyłu można obecnie przewidywać zakażenie u pojedynczej osoby ani że wyniki powinny zmieniać indywidualne leczenie pacjentów. Potencjalne zastosowanie dotyczy przede wszystkim epidemiologii, zdrowia publicznego i planowania systemowego.

Badanie wzmacnia także argumenty za ograniczaniem ekspozycji populacji na PM_2,5. W analizie kontrfaktycznej dotyczącej prowincji Gauteng autorzy oszacowali, że przy utrzymaniu stężeń na niższych poziomach część przypadków IChP mogłaby teoretycznie nie wystąpić. Takie wyliczenia są jednak wynikiem modelowania, a nie bezpośrednio zaobserwowanym efektem interwencji środowiskowej.

W praktyce klinicznej szczególnie istotne jest rozróżnienie między epidemiologiczną zależnością a indywidualnym ryzykiem. Publikacja nie ustanawia nowego testu prognostycznego ani nie wskazuje nowej metody leczenia. Poszerza natomiast wiedzę o tym, jak środowisko, gospodarz i biologiczne zróżnicowanie pneumokoków mogą współtworzyć czasowe i przestrzenne wzorce ciężkich zakażeń.