Популяция исследования и характеристики образцов
В окончательную когорту анализа вошли 112 субъектов, из которых 64 человека находились в группе инфекции и 48 человек — в контрольной группе без инфекции. Исходные демографические показатели были схожи в обеих группах, в то время как лабораторные параметры, связанные с инфекцией, были выше в группе инфекции, как показано в Таблице 1. Инфекции были вызваны преимущественно грамотрицательными бактериями, основными бактериальными патогенами которых являлись Klebsiella pneumoniae, Pseudomonas aeruginosa и Acinetobacter baumannii. Другими патогенами были Streptococcus pneumoniae, Staphylococcus aureus, атипичные агенты, вирусы и смешанные инфекции.
Показатели восстановления БАЛЖ соответствовали заранее установленному порогу ≥40% для всех сохраненных образцов, при этом показатели в группе инфекции были немного ниже, чем в контрольной группе (медиана 48,3% против 53,3%, p = 0,03), что согласуется с усилением воспаления дыхательных путей и застойными явлениями в слизистой оболочке инфицированных дыхательных путей. Поскольку восстановленный объем также различался между группами (58 mL против 64 mL, p = 0,04), все сообщаемые концентрации цитокинов могут быть подвержены незначительной межгрупповой разнице в коэффициенте разведения; предварительно заданный анализ чувствительности, ограниченный образцами с показателем восстановления ≥50% (n = 42 в группе инфекции, n = 34 в контрольной группе), воспроизвел направление и величину каждого межгруппового сравнения (максимальное изменение величины эффекта < 8%), что указывает на маловероятность влияния ошибки разведения на наблюдаемые различия. Время между сбором и обработкой составило 65 min для группы инфекции и 62 min для контрольной группы, оставаясь в пределах 2 ч окна. Критерии контроля качества (КК) на уровне образцов указали на высокую достоверность результатов. Цитологический дифференциальный подсчет был получен на основе подсчета ≥400 клеток во всех 112/112 образцах (100%). Результаты проточной цитометрии соответствовали заранее установленному критерию ≥1 × 105 событий CD45⁺ во всех 112/112 образцах (100%) при медианной жизнеспособности живых клеток 89,4% (IQR 85,6–92,8%). Все мультиплексные планшеты для анализа цитокинов соответствовали критериям КК после одного повторного запуска. Для всех 13 аналитов внутрисерийный коэффициент вариации (CV) составил 4,7% (2,1–9,3%, все <10%), межпланшетный CV объединенных образцов КК составил 11,2% (6,8–18,4%, все <20%), а медианное число бусин составило 78 событий/аналит/лунка (53–142, все ≥ 50).
| Переменная | Инфекция (n = 64) | Контрольная группа (n = 48) | p-значение |
| Возраст, лет, медиана [МКР] | 62 [54–71] | 60 [52–68] | 0.41 |
| Мужчины, n (%) | 38 (59.4) | 26 (54.2) | 0.59 |
| ИМТ, кг/м²2, медиана [МКР] | 23.8 [21.4–26.2] | 23.5 [21.6–25.7] | 0.66 |
| Анамнез курения, n (%) | | | 0.52 |
| Никогда | 28 (43.8) | 25 (52.1) | |
| Предыдущий | 22 (34.4) | 16 (33.3) | |
| Текущий | 14 (21.9) | 7 (14.6) | |
| Артериальная гипертензия, n (%) | 27 (42.2) | 18 (37.5) | 0.62 |
| Сахарный диабет, n (%) | 17 (26.6) | 10 (20.8) | 0.49 |
| Ишемическая болезнь сердца, n (%) | 11 (17.2) | 6 (12.5) | 0.51 |
| Хроническая болезнь почек, n (%) | 5 (7.8) | 3 (6.3) | 0.76 |
| Показание к бронхоскопии, n (%) | | | — |
| Внебольничная пневмония | 35 (54.7) | — | |
| Госпитальная пневмония | 19 (29.7) | — | |
| Вентилятор-ассоциированная пневмония | 10 (15.6) | — | |
| Неопределенный легочный узел | — | 23 (47.9) | |
| Стойкая рентгенологическая патология | — | 15 (31.3) | |
| Диагностика хронического кашля | — | 10 (20.8) | |
| Возбудитель, n (%) | | — | |
| Грамотрицательные бактерии | 28 (43.8) | | |
| Грамположительные бактерии | 16 (25.0) | | |
| Атипичные патогены | 7 (10.9) | | |
| Респираторные вирусы | 8 (12.5) | | |
| Смешанная инфекция | 5 (7.8) | | |
| ИВЛ при БАЛ, n (%) | 14 (21.9) | 0 (0) | <0.001 |
| Прием антибиотиков в течение 72 ч до проведения БАЛ, n (%) | 41 (64.1) | 0 (0) | <0.001 |
| Введенный объем, мл, медиана [МКР] | 120 [120–140] | 120 [120–140] | 0.84 |
| Восстановленный объем, мл, медиана [МКР] | 58 [49–68] | 64 [55–73] | 0.04 |
| Процент восстановления БАЛЖ, %, медиана [МКР] | 48.3 [42.5–55.6] | 53.3 [46.2–60.5] | 0.03 |
| Время обработки, мин, медиана [МКР] | 65 [50–85] | 62 [48–80] | 0.55 |
| Лейкоциты крови, × 10⁹/л, медиана [МР] | 11.8 [8.6–14.5] | 6.4 [5.3–7.8] | <0.001 |
| Нейтрофилы крови, %, медиана [МКИ] | 78.5 [70.2–86.4] | 60.3 [54.6–66.8] | <0.001 |
| С-реактивный белок, мг/л, медиана [МР] | 78.5 [42.6–128.4] | 4.2 [1.8–8.6] | <0.001 |
| Прокальцитонин, нг/мл, медиана [МКИ] | 0.86 [0.32–2.45] | 0.05 [0.03–0.08] | <0.001 |
Таблица 1: Демографические, клинические и характеристики обработки BALF когорты исследования. Сравниваются исходные демографические характеристики, сопутствующие заболевания, показания к бронхоскопии, возбудители инфекций, параметры процедур и лабораторные показатели периферической крови между 64 участниками группы с инфекцией и 48 участниками контрольной группы без инфекции. Непрерывные переменные представлены в виде медианы [межквартильный размах] и сравнивались с использованием U-критерия Манна-Уитни; категориальные переменные представлены как n (%) и сравнивались с использованием критерия хи-квадрат или точного критерия Фишера, где это применимо. Сокращения: BAL = бронхоальвеолярный лаваж; BALF = бронхоальвеолярная лаважная жидкость; BMI = индекс массы тела; IQR = межквартильный размах; WBC = количество лейкоцитов.
Визуализация рекрутирования нейтрофилов в БАЛЖ
У инфицированных субъектов наблюдался явный сдвиг соотношения популяций от преобладания макрофагов к преобладанию нейтрофилов, как показано на Рисунке 2. Препараты, полученные с помощью цитоцентрифуги, окрашивали по методу Райта — Гимзы; окрашивание выявило четкое различие между двумя группами: в контрольной группе преобладали альвеолярные макрофаги, присутствовали некоторые лимфоциты и несколько нейтрофилов, тогда как в инфицированной группе доминировали нейтрофилы с частыми внутриклеточными бактериями, дегенерацией ядер и токсической грануляцией. Для подтверждения различий использовали дифференциальный подсчет не менее 400 клеток, в результате чего доля нейтрофилов увеличилась с медианного значения 3.2% (IQR: 1.8–5.6%) в контрольной группе до 68.4% (IQR: 49.2–82.5%) в инфицированной группе (p < 0.001), тогда как доля альвеолярных макрофагов снизилась с 88.2% (IQR: 82.4–92.5%) до 22.6% (IQR: 12.3–38.5%) (p < 0.001). Процентное содержание лимфоцитов было схожим в обеих группах (7.4% [4.8–11.3%] против 6.5% [3.4–11.2%], p = 0.42).
Анализ гейтирования и снижение размерности подтвердили и расширили данные о клеточных фенотипах, наблюдаемых при микроскопии. После исключения детрита, дуплетов и мертвых клеток доля CD45⁺CD66b⁺ гранулоцитов составила в медиане 65.8% (IQR: 46.2–79.4%) от всех CD45⁺ лейкоцитов в группе инфекции по сравнению с 4.1% (IQR: 2.3–6.8%) в контрольной группе (p < 0.001). В популяции гранулоцитов CD15⁺CD16⁺CD11b⁺ зрелые нейтрофилы составили 92.4% (IQR: 87.5–95.8%) всех гранулоцитов в группе инфекции против 88.6% (IQR: 82.3–93.5%) в контрольной группе, что указывает на то, что пул привлеченных нейтрофилов состоял преимущественно из зрелых, полностью дифференцированных клеток, а не из ранних предшественников. Проекции UMAP и t-SNE, основанные на матрице поверхностных белков с arcsinh-преобразованием, показали почти идентичную топологию, при которой образцы с инфекцией занимали выраженный «остров» CD15⁺CD16⁺CD11b⁺ зрелых нейтрофилов, который был слабо заселен в контрольных образцах, в то время как «остров» HLA-DR⁺CD14⁺ альвеолярных макрофагов демонстрировал противоположную картину: высокая плотность в контроле и низкая в группе инфекции. Мета-кластеризация с помощью FlowSOM выявила 13 субпопуляций, среди которых наибольшее межгрупповое расширение продемонстрировали субпопуляция зрелых нейтрофилов (38.4% [22.6–56.2%] против 2.8% [1.5–4.6%], p < 0.001) и кластер активированных CD16-dim CD11b-high нейтрофилов (22.6% [11.4–32.5%] против 0.4% [0.1–1.2%], p < 0.001).
Что касается абсолютного числа нейтрофилов, в группе инфекции наблюдалось примерно 50-кратное увеличение количества нейтрофилов в БАЛЖ по сравнению с контрольной группой. Общее количество ядерных клеток на µL БАЛЖ увеличилось со 156 cells/µL (IQR 88–232) в контрольной группе до 384 cells/µL (IQR 180–712) в группе инфекции (p < 0.001), а число нейтрофилов выросло с 5 cells/µL (IQR 2–13) в контрольной группе до 262 cells/µL (IQR 89–578) в группе инфекции (p < 0.001). Попарные различия в проценте и абсолютном количестве нейтрофилов между группами оставались высокозначимыми после поправки Бенджамини — Хохберга (BH) (q < 0.001).

Рисунок 2: Визуализация привлечения нейтрофилов в БАЛЖ. Представительные поля цитоспина, окрашенные по методу Райта — Гимзы, от участника контрольной группы и инфицированного участника иллюстрируют переход от клеточного профиля БАЛЖ с преобладанием макрофагов к профилю с преобладанием нейтрофилов. Рядом представлены четырехэтапная стратегия гейтирования при проточной цитометрии (синглеты FSC → живые CD45⁺ лейкоциты → CD66b⁺ гранулоциты после исключения CD14⁺HLA-DR⁺ → CD15⁺CD16⁺CD11b⁺ зрелые нейтрофилы), UMAP-вложение данных проточной цитометрии с arcsinh-преобразованием и наложенными 13 метакластерами FlowSOM, раскрашенными по биологической аннотации, а также диаграммы размаха процентов нейтрофилов в БАЛЖ и абсолютной концентрации нейтрофилов (клетки/µL, логарифмическая шкала), что позволяет напрямую сравнить клеточные фенотипы, выявленные в группах инфекции и контроля. Сокращения: БАЛЖ = бронхоальвеолярный лаважный жидкость; FSC = прямое рассеяние; UMAP = равномерная аппроксимация и проекция многообразия. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Визуализация профилей воспалительных цитокинов в БАЛЖ
Характер скоординированных ответов по нескольким цитокиновым осям позволил отличить группу инфекции от контрольной группы; при этом цитокины были разделены на три устойчивых функциональных модуля, как показано на Рисунке 3. Кластеризация нормализованных концентраций цитокинов, преобразованных по логарифмической шкале log₁₀, с использованием метода Ward.D2 привела к созданию четкой дендрограммы образцов, разделенной на две группы, где образцы с инфекцией и контрольные образцы располагались на противоположных концах, при этом было отмечено всего семь ошибочно классифицированных участников (шесть контрольных образцов сгруппировались с образцами инфекции, в основном случаи хронического кашля с легким нейтрофильным воспалением дыхательных путей, и один образец инфекции сгруппировался с контролем, что соответствовало ранней атипичной пневмонии). Дендрограмма аналитов определила три кластера цитокинов. Хемокиновый кластер, включающий CXCL8, CXCL1, CXCL10, CCL2 и CCL3, продемонстрировал наибольшее кратное изменение между группой инфекции и контрольной группой, а также самую высокую внутримодульную корреляцию. Провоспалительный кластер, содержащий IL-1β, IL-6, IL-17A и TNF-α, показал схожие направленные повышения, хотя и несколько менее выраженные. Регуляторный кластер, включающий IL-10, IL-1RA, IFN-γ и G-CSF, также был повышен в образцах с инфекцией, что указывает на то, что ответ не был чисто провоспалительным, а параллельно задействовал контррегуляторные и противовирусные механизмы.
Медианы (МКИ) всех исследованных аналитов были выше в группе с инфекцией по сравнению с контрольной группой, а скорректированное значение p было меньше 0,001 для всех 13 аналитов после коррекции ложноположительных результатов по методу Бенджамини — Хохберга (BH-FDR). Среди аналитов CXCL8/IL-8 имел самую высокую концентрацию и наиболее выраженное кратное изменение (с 68 [32–142] pg/mL в контрольной группе до 2 485 [842–6 520] pg/mL в группе с инфекцией, что соответствует примерно 36-кратному увеличению). Аналогичная закономерность наблюдалась для CXCL1/GRO-α с примерно 21-кратным увеличением (78 [36–148] против 1 672 [524–3 890] pg/mL). Для IL-1β наблюдалось примерно 13-кратное увеличение (5,8 [2,6–11,4] против 78 [29–184] pg/mL). IL-6, IL-1RA и G-CSF продемонстрировали увеличение в 13–16 раз (16,4 [7,8–30,5] против 248 [82–612] pg/mL; 112 [48–258] против 1 748 [658–3 964] pg/mL; 24 [10–52] против 384 [125–912] pg/mL соответственно). Другие, менее выраженные увеличения включали TNF-α (~8-кратное), CXCL10/IP-10 (~9-кратное), CCL2/MCP-1 (~8-кратное), CCL3/MIP-1α (~10-кратное), IL-10 (~10-кратное), IFN-γ (~6-кратное) и IL-17A (~4-кратное).
Анализ ROC-кривых внутри когорты позволил определить IL-8 (AUC 0,94), IL-1β (AUC 0,92), IL-6 (AUC 0,91), CXCL1 (AUC 0,91) и G-CSF (AUC 0,89) как аналиты с наибольшей дискриминирующей способностью одиночных маркеров при разграничении группы инфекции и контрольной группы; данные значения AUC отражают дискриминацию внутри анализируемой когорты в отсутствие отдельного внешнего валидационного набора. Анализ чувствительности с использованием множественной импутации методом последовательных уравнений (m = 20) дал оценки в пределах 5% от медиан для полных случаев по каждому аналиту и не изменил ранговый порядок величин эффекта.
Для оценки влияния предшествующей антибиотикотерапии и подтипа пневмонии на наблюдаемые различия в уровнях цитокинов были проведены заранее определенные анализы чувствительности. При повторном сравнении основных групп после исключения 41 участника группы инфекции (64,1%), получавших антибиотики в течение 72 ч до проведения БАЛ, все 13 аналитов оставались значимо более высокими в оставшихся 23 случаях инфекции, чем в контрольной группе, после коррекции BH-FDR (скорректированное p < 0,001 для пяти основных дискриминаторов); при этом точечные оценки были несколько выше, а не ниже медиан всей когорты (например, CXCL8/IL-8 2 860 [980–7 110] pg/mL и IL-1β 92 [34–210] pg/mL), что согласуется с частичным подавлением цитокинов дыхательных путей вследствие приема антибиотиков перед БАЛ. Стратификация группы инфекции по подтипу пневмонии показала одинаковый порядок ранжирования хемокинов и провоспалительных цитокинов для внебольничной пневмонии (n = 35), госпитальной пневмонии (n = 19) и вентилятор-ассоциированной пневмонии (n = 10), причем в подгруппе с вентилятор-ассоциированной пневмонией наблюдались самые высокие медианные уровни CXCL8, CXCL1, IL-6 и G-CSF; статус механической вентиляции легких не изменил направление ни одного межгруппового сравнения.

Рисунок 3: Визуализация профилей воспалительных цитокинов в БАЛЖ. Тепловая карта иерархической кластеризации лог₁₀-преобразованных и z-стандартизированных концентраций 13 измеренных аналитов у всех 112 участников; цветные полосы аннотаций над тепловой картой указывают на принадлежность к группе (контроль, инфекция), подтип пневмонии (внебольничная, госпитальная или вентилятор-ассоциированная пневмония) и категорию возбудителя (грамотрицательная, грамположительная, атипичная, вирусная или смешанная инфекция). Диаграммы размаха для наиболее значимых дифференциально экспрессируемых цитокинов (CXCL8/IL-8, IL-1β, IL-6, CXCL1/GRO-α, G-CSF и TNF-α; лог₁₀ шкала, pg/mL) позволяют сравнить группы инфекции и контроля. *** p < 0.001 после коррекции по методу Бенджамини–Хохберга (BH) для контроля ложных открытий. Сокращения: БАЛЖ = бронхоальвеолярный лаважный жидкость; BH = Бенджамини–Хохберг; CAP = внебольничная пневмония; HAP = госпитальная пневмония; VAP = вентилятор-ассоциированная пневмония. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Совместная визуализация рекрутирования нейтрофилов и профилей цитокинов
Сопоставление клеточных и растворимых фракций выявило сильную биологическую связь между привлечением нейтрофилов и системой хемокинов, как показано на рисунке 4. Двумерное представление было получено с помощью метода главных компонентов (PCA) стандартизированной матрицы клеток и цитокинов, где ГК1 (PC1) составила 51,1% общей дисперсии, а ГК2 (PC2) — еще 12,7%, что в совокупности дает 63,8%. ГК1 разделила группу инфекции и контрольную группу; группа инфекции характеризовалась высокими положительными значениями для нейтрофилов в БАЛЖ, CXCL8, CXCL1, IL-6, IL-1β и G-CSF, и высокими отрицательными значениями для макрофагов. ГК2 разграничила траекторию с высоким содержанием хемокинов (высокие положительные значения для CXCL8 и CXCL1) и траекторию с преобладанием регуляторных факторов (положительные значения для IL-10, IL-1RA и IFN-γ), таким образом, что внутри группы инфекции участники имели тенденцию разделяться по степени тяжести заболевания: случаи вентилятор-ассоциированной пневмонии с применением искусственной вентиляции легких попадали в квадрант с высоким уровнем хемокинов.
Попарный корреляционный анализ рангов Спирмена между 17 клеточными и растворимыми признаками выявил сильную положительную корреляцию между количеством нейтрофилов в БАЛЖ и общим хемокиновым каскадом (CXCL8 ρ = 0.78; CXCL1 ρ = 0.74; CCL3 ρ = 0.69; CCL2 ρ = 0.61; CXCL10 ρ = 0.54; все скорректированные p < 0.001). Второй уровень значимых корреляций был обнаружен между количеством нейтрофилов и провоспалительным каскадом (G-CSF ρ = 0.65; IL-1β ρ = 0.62; IL-6 ρ = 0.58; TNF-α ρ = 0.51), за которым следовали умеренные положительные корреляции с регуляторным каскадом (IL-10 ρ = 0.45; IL-1RA ρ = 0.48). Процент макрофагов демонстрировал зеркальную отрицательную корреляцию аналогичной силы как с хемокиновым, так и с провоспалительным каскадами. Внутрикаскадные корреляции были существенно сильнее, чем межкаскадные (медиана |ρ| внутри каскада = 0.68 против медианы |ρ| между каскадами = 0.33), что подтверждает концепцию общих вышестоящих регуляторных путей, контролирующих хемокиновый, провоспалительный и регуляторный каскады.
Сеть корреляций отобразила эту связь на графе, в котором нейтрофилы БАЛ являлись наиболее центральным узлом, причем все 13 ребер, ведущих к цитокинам, сохранили значимость после коррекции Бенджамина-Хохберга (BH) для множественного тестирования. Из 136 возможных попарных связей между 17 признаками 99 соответствовали совместному критерию |ρ| ≥ 0.30 и скорректированному значению p BH < 0.05 и были сохранены. Граф сети корреляций продемонстрировал четкую трехблочную модульную архитектуру, которая повторяла кластеры, полученные в результате иерархической кластеризации: самые плотные и сильные связи наблюдались между модулем хемокинов и количеством и процентом нейтрофилов в БАЛ; связи средней интенсивности соединяли провоспалительный модуль как с нейтрофилами, так и с модулем хемокинов; а слабые связи регуляторного модуля служили мостиком между блоками хемокинов и провоспалительных факторов. Узел с процентом макрофагов располагался напротив блока «нейтрофилы — хемокины» и был связан с ним веером ребер с отрицательным весом. В совокупности PCA, корреляционная матрица и сетевой граф выявили согласованную закономерность, при которой инфильтрация нейтрофилов в БАЛ и ответ хемокинов тесно взаимосвязаны на уровне отдельных участников, а провоспалительные и регуляторные пути функционируют как четко определенные, но взаимосвязанные модули вокруг этой оси «нейтрофилы — хемокины».

Рисунок 4: Совместная визуализация рекрутирования нейтрофилов и профилей воспалительных цитокинов. Биплот анализа главных компонентов стандартизированной матрицы клеточных и цитокиновых данных с цветовой разметкой участников по группам, тепловая карта корреляции Спирмена и сетевой граф всех 17 клеточных и растворимых признаков (размер узлов соответствует степени связности, ребра масштабированы по |ρ| и отфильтрованы по BH-скорректированному p < 0.05 и |ρ| ≥ 0.30), выделяющий хаб нейтрофилы–хемокины. Сокращения: BH = Бенджамини–Хохберг; PC = главная компонента. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
ДОСТУПНОСТЬ ДАННЫХ:
Обезличенный набор данных на уровне участников и код анализа на языке R, использованные в данном исследовании, представлены в качестве Дополнительная таблица S1 и Дополнительный файл 1Другие данные, подтверждающие результаты данного исследования, могут быть предоставлены автором, ответственным за переписку, по обоснованному запросу.
Дополнительная таблица S1: Набор данных BALF. Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 1: ZIP-архив, содержащий коды для анализа. Пожалуйста, нажмите здесь, чтобы скачать этот файл.