Популяция исследования и характеристики образцов
В итоговую когорту анализа вошли 12 субъектов, из которых 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-часовое окно. Критерии контроля качества (К) на уровне образцов свидетельствовали о высокой достоверности результатов. Дифференциальный подсчет цитологических показателей был получен на основе подсчета ≥40 клеток во всех 112/12 образцах (10%). Данные проточной цитометрии соответствовали заранее установленному критерию ≥1 × 105 событий CD45⁺ во всех 112/12 образцах (10%), при этом медиана жизнеспособности живых клеток составила 89,4% (IQR 85,6–92,8%). Все мультиплексные планшеты для анализа цитокинов соответствовали критериям КК после одного повторного запуска. Для всех 13 аналитов внутрисерийный коэффициент вариации (CV) составил 4,7% (2,1–9,3%, все <10%), межпланшетный CV пулированных образцов КК составил 1,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: Демографические, клинические характеристики и параметры обработки БАЛЖ когорты исследования. Сравниваются исходные демографические характеристики, сопутствующие заболевания, показания к бронхоскопии, возбудители инфекции, параметры процедуры и показатели лабораторных анализов периферической крови у 64 участников группы с инфекцией и 48 участников контрольной группы без инфекции. Непрерывные переменные представлены в виде медианы [межквартильный размах] и сравнивались с помощью U -критерия Манна–Уитни; категориальные переменные представлены как n (%) и сравнивались с помощью критерия хи-квадрат или точного критерия Фишера, где это применимо. Сокращения: БАЛ = бронхоальвеолярный лаваж; БАЛЖ = бронхоальвеолярный лаважный жидкость; ИМТ = индекс массы тела; IQR = межквартильный размах; WBC = количество лейкоцитов.
Визуализация рекрутирования нейтрофилов в БАЛЖ
Как показано на Рисунке 2, у инфицированных субъектов наблюдался отчетливый сдвиг соотношения популяций от преобладания макрофагов к преобладанию нейтрофилов. Препараты, полученные с помощью цитоцентрифугирования, были окрашены по методу Райта — Гимзы; окрашивание выявило явное различие между двумя группами: в контрольной группе преобладали альвеолярные макрофаги, присутствовали некоторые лимфоциты и несколько нейтрофилов, в то время как в инфицированной группе доминировали нейтрофилы с частым наличием внутриклеточных бактерий, дегенерированными ядрами и токсической грануляцией. Для подтверждения различий использовался дифференциальный подсчет не менее 40 клеток, в результате чего доля нейтрофилов увеличилась с медианного значения 3.2% (IQR: 1.8–5.6%) в контрольной группе до 68.4% (IQR: 49.2–82.5%) в инфицированной группе (p < 0.01), тогда как доля альвеолярных макрофагов снизилась с 8.2% (IQR: 82.4–92.5%) до 2.6% (IQR: 12.3–38.5%) (p < 0.01). Процент лимфоцитов был схожим в обеих группах (7.4% [4.8–1.3%] против 6.5% [3.4–1.2%], p = 0.42).
Анализ гейтирования и снижение размерности подтвердили и дополнили данные о клеточных фенотипах, полученные с помощью микроскопии. После исключения клеточного дебриса, дуплетов и мертвых клеток доля CD45⁺CD6b⁺ гранулоцитов составила в медиане 65,8% (IQR: 46,2–79,4%) от всех CD45⁺ лейкоцитов в группе инфекции по сравнению с 4,1% (IQR: 2,3–6,8%) в контрольной группе (p < 0.01). В популяции гранулоцитов зрелые нейтрофилы CD15⁺CD16⁺CD1b⁺ составили 92,4% (IQR: 87,5–95,8%) всех гранулоцитов в группе инфекции против 8,6% (IQR: 82,3–93,5%) в контрольной группе, что указывает на то, что пул привлеченных нейтрофилов состоял преимущественно из зрелых, полностью дифференцированных клеток, а не из ранних предшественников. Проекции UMAP и t-SNE на основе матрицы поверхностных белков с arcsinh-преобразованием показали почти идентичные топологии, в которых случаи инфекции занимали выраженный «остров» зрелых нейтрофилов CD15⁺CD16⁺CD1b⁺, который был слабо заселен в контрольных образцах, тогда как «остров» HLA-DR⁺CD14⁺ альвеолярных макрофагов демонстрировал противоположную картину: плотное заселение в контроле и редкое в группе инфекции. Мета-кластеризация с помощью FlowSOM выявила 13 субпопуляций, среди которых наибольшее межгрупповое расширение продемонстрировали субпопуляция зрелых нейтрофилов (38,4% [2,6–56,2%] против 2,8% [1,5–4,6%], p < 0.001) и кластер активированных CD16-dim CD1b-high нейтрофилов (2,6% [1,4–32,5%] против 0,4% [0,1–1,2%], p < 0.001).
Что касается абсолютного числа нейтрофилов, в группе инфекции наблюдалось примерно 50-кратное увеличение количества нейтрофилов в БАЛЖ по сравнению с контрольной группой. Общее количество ядерных клеток на µL БАЛЖ увеличилось со 156 cells/µL (IQR 8–232) в контрольной группе до 384 cells/µL (IQR 180–712) в группе инфекции (p < 0.01), а число нейтрофилов возросло с 5 cells/µL (IQR 2–13) в контрольной группе до 262 cells/µL (IQR 89–578) в группе инфекции (p < 0.01). Попарные различия в процентном содержании и абсолютном количестве нейтрофилов между группами оставались высокозначимыми после поправки Бенджамини–Хохберга (BH) (q < 0.01).

Рисунок 2Визуализация привлечения нейтрофилов в БАЛЖ. Репрезентативные поля цитоспинов, окрашенных по методу Райта — Гимзы, полученные от контрольного участника и инфицированного участника, иллюстрируют сдвиг клеточного профиля БАЛЖ от преобладания макрофагов к преобладанию нейтрофилов. Четырехэтапная стратегия гейтирования при проточной цитометрии (синглеты FSC → живые CD45⁺ лейкоциты → CD6b⁺ гранулоциты после исключения CD14⁺HLA-DR⁺ → CD15⁺CD16⁺CD1b⁺ зрелые нейтрофилы), UMAP-проекция данных проточной цитометрии с arcsinh-преобразованием и наложением 13 метакластеров FlowSOM, расцвеченных по биологической аннотации, а также диаграммы размаха процента нейтрофилов в БАЛЖ и абсолютной концентрации нейтрофилов (клеток/µL, логарифмическая шкала) представлены рядом для прямого сравнения клеточных фенотипов, выявленных в группах инфекции и контроля. Сокращения: BALF = бронхоальвеолярный лаважный fluid (БАЛЖ); 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,01 для всех 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–1,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; 12 [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.01 для пяти ведущих дискриминаторов), при этом точечные оценки были несколько выше, а не ниже медиан всей когорты (например, CXCL8/IL-8 2 860 [980–7 10] pg/mL и IL-1β 92 [34–210] pg/mL), что согласуется с частичным подавлением цитокинов дыхательных путей из-за применения антибиотиков перед БАЛ. Стратификация группы инфекции по подтипу пневмонии показала одинаковый порядок ранжирования хемокинов и провоспалительных цитокинов для внебольничной пневмонии (n = 35), нозокомиальной пневмонии (n = 19) и вентилятор-ассоциированной пневмонии (n = 10), при этом в подгруппе вентилятор-ассоциированной пневмонии наблюдались самые высокие медианные уровни CXCL8, CXCL1, IL-6 и G-CSF; факт проведения искусственной вентиляции легких не изменил направление какого-либо из межгрупповых сравнений.

Рисунок 3Визуализация профилей воспалительных цитокинов в БАЛЖ. Тепловая карта иерархической кластеризации log₁₀-преобразованных и z-стандартизированных концентраций 13 измеренных аналитов для всех 112 участников; цветные аннотационные полосы над тепловой картой указывают принадлежность к группе (контрольная группа, инфекция), подтип пневмонии (внебольничная, госпитальная или вентилятор-ассоциированная пневмония) и категорию возбудителя (грамотрицательная, грамположительная, атипичная, вирусная или смешанная инфекция). Диаграммы размаха для наиболее значимых дифференциально экспрессируемых цитокинов (CXCL8/IL-8, IL-1β, IL-6, CXCL1/GRO-α, G-CSF и TNF-α; шкала log₁₀, пг/мл) для сравнения групп инфекции и контроля. *** p < 0,01 после коррекции ложноположительных результатов по методу Бенджамини–Хохберга. Сокращения: BALF = бронхоальвеолярный лаваж; BH = Бенджамини–Хохберг; CAP = внебольничная пневмония; HAP = госпитальная пневмония; VAP = вентилятор-ассоциированная пневмония. Нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Совместная визуализация рекрутирования нейтрофилов и профилей цитокинов
Сопоставление клеточных и растворимых компонентов выявило сильную биологическую связь между привлечением нейтрофилов и системой хемокинов, как показано на Рисунке 4. Двумерное представление было получено с помощью анализа главных компонентов (PCA) стандартизированной матрицы клеток и цитокинов, где PC1 составил 51,1% общей дисперсии, а PC2 — еще 12,7%, что в сумме дает 63,8%. PC1 отделил группу с инфекцией от контрольной группы и характеризовался высокими положительными нагрузками на нейтрофилы БАЛЖ, CXCL8, CXCL1, IL-6, IL-1β и G-CSF, а также высокими отрицательными нагрузками на макрофаги. PC2 разделил траекторию с высоким содержанием хемокинов (высокие положительные нагрузки на 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.01). Второй уровень значимых корреляций был обнаружен между количеством нейтрофилов и провоспалительным каскадом (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 признаками 9 соответствовали совместному критерию |ρ| ≥ 0.30 и скорректированному по BH p < 0.05 и были сохранены. Граф сети корреляций продемонстрировал четкую трехблочную модульную архитектуру, которая повторяла кластеры, полученные в результате иерархической кластеризации: самые плотные и сильные связи наблюдались между модулем хемокинов и количеством и процентом нейтрофилов в БАЛ; связи средней интенсивности связывали провоспалительный модуль как с нейтрофилами, так и с модулем хемокинов, а слабые связи из регуляторного модуля соединяли блоки хемокинов и провоспалительных факторов. Узел процента макрофагов располагался напротив блока нейтрофилов и хемокинов и был связан с ним веером ребер с отрицательным весом. В совокупности PCA, корреляционная матрица и сетевой граф выявили согласованную закономерность, при которой инфильтрация нейтрофилов в БАЛ и ответ хемокинов сильно связаны на уровне каждого отдельного участника, а провоспалительные и регуляторные пути служат четко определенными, но взаимосвязанными модулями вокруг этой оси нейтрофилов и хемокинов.

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