Исследовательская статья

Интегрированная визуализация рекрутирования нейтрофилов и профилирование воспалительных цитокинов в бронхоальвеолярной лаважной жидкости

0 просмотров

DOI:

10.3791/72246

15 сентября 2026 г.

В этой статье

Краткое содержание

В данном исследовании представлен комплексный рабочий процесс визуализации, объединяющий стандартизированный сбор и обработку жидкости для бронхоальвеолярного лаважа, панель проточной цитометрии с семью маркерами, мультиплексный анализ цитокинов с 13 аналитами и единый конвейер анализа.

Аннотация

Жидкость для бронхоальвеолярного лаважа (БАЛ) является наиболее доступным матриксом для забора проб при изучении воспаления нижних дыхательных путей, однако традиционные анализы опираются на проточную цитометрию с использованием одного маркера в сочетании с иммуноферментным анализом (ELISA) для каждого отдельного аналита; такой подход ограничивает плотность получаемой информации и скрывает взаимосвязь между клетками и факторами. В данном исследовании представлен интегрированный рабочий процесс визуализации, сочетающий стандартизированный сбор и предварительную обработку с панелью для проточной цитометрии нейтрофилов и моноцитов по семи маркерам, мультиплексным анализом 13 цитокинов и унифицированным конвейером анализа на основе равномерного аппроксимированного проекционного отображения (UMAP), метакластеризации FlowSOM, иерархической кластеризации, сетей корреляции Спирмена и анализа главных компонентов (PCA).

Данный рабочий процесс был применен к одноцентровой проспективной когорте из 112 взрослых пациентов, перенесших бронхоскопию, которая включала 64 пациентов с микробиологически подтвержденной инфекцией легких и 48 лиц из контрольной группы без инфекции. Данные цитологического анализа и проточной цитометрии выявили выраженный сдвиг альвеолярного профиля от макрофагального к нейтрофильному при инфекции: процент нейтрофилов увеличился с медианы 3,2% до 68,4%, а абсолютная концентрация нейтрофилов возросла примерно в 50 раз. Мультиплексный профильный анализ выявил скоординированный трехмодульный цитокиновый ответ, при котором хемокиновый, провоспалительный и регуляторный модули были активированы; при этом CXCL8/IL-8, IL-1β, IL-6, CXCL1/GRO-α и G-CSF продемонстрировали наибольшую дискриминационную способность в качестве одиночных маркеров в исследуемой когорте (площадь под кривой [AUC] от 0,89 до 0,94). Интегративная визуализация с использованием PCA, корреляционной матрицы и сетевого графа выявила тесную связь оси «нейтрофилы — хемокины» на уровне участников, при этом количество нейтрофилов в БАЛЖ оказалось центральным узлом сети.

Данный рабочий процесс использует только материалы, регулярно получаемые в ходе клинически обоснованного БАЛ, не требует дополнительного объема образцов и предоставляет информативно насыщенную основу для анализа воспаления дыхательных путей при респираторных инфекциях. Это может быть реализовано в лабораториях специализированной медицинской помощи, оснащенных системами многопараметрической проточной цитометрии и платформами мультиплексного микросферного анализа, с возможностью применения в лонгитюдных исследованиях, анализах единичных клеток и оценке терапевтической эффективности. Более широкое внедрение отдельных аналитов в клиническую практику в качестве биомаркеров потребует проспективной валидации в независимых многоцентровых когортах.

Введение

Инфекционные смерти в результате инфекций нижних дыхательных путей остаются значительной причиной глобальной смертности, несмотря на достижения в стратегиях вакцинации и методах лечения; тем не менее, это особенно актуально для детского и пожилого населения1. Нейтрофилы являются основными клетками, регистрируемыми на ранних стадиях врожденного иммунного ответа в пораженной дыхательной системе, что жизненно важно для уничтожения вторгающегося патогена; однако неконтролируемое нейтрофильное воспаление приводит к разрушению альвеолярно-капиллярного барьера, окислительному повреждению тканей и развитию острого респираторного дистресс-синдрома (ОРДС)2,3. Нейтрофилы в воздушном пространстве легких также демонстрируют кинетику привлечения, особенности удерживания и эффекторное поведение, которые отличаются от других тканей, и эти различия влияют на то, разрешится ли инфекция или перейдет в тяжелую пневмонию4.

Жидкость для бронхоальвеолярного лаважа (БАЛЖ) является важным источником информации об альвеолярном микроокружении и обеспечивает наиболее точное представление о клиническом иммунитете в данной области5. Хемокин CXCL8 играет ключевую роль в привлечении нейтрофилов в легкие через его рецепторы CXCR1 и CXCR2, наряду с другими хемокинами, включая CXCL1, CXCL2 и CXCL56. Цитокины IL-1β, IL-6 и TNF-α дополнительно усиливают этот процесс, в то время как IL-10 действует как противодействующий цитокин. Было показано, что их уровни в БАЛЖ коррелируют с тяжестью инфекций легких и развитием осложнений в легких вследствие бактериальной, вирусной и атипичной пневмонии7, а мультиплексные профили биомаркеров БАЛЖ внедряются для более ранней диагностики и дифференциации патогенов при инфекциях нижних дыхательных путей8.

Стандартное исследование БАЛЖ обычно не позволяет полностью охватить биологические механизмы, лежащие в основе теста. Количество клеток оценивают путем ручного дифференцирования клеток в образцах после цитоспинирования, в то время как растворимые медиаторы измеряют с помощью отдельных иммуноферментных анализов (ELISA) и представляют в виде столбчатых диаграмм. Хемокины, цитокины и каждый тип клеток оцениваются по отдельности, при этом между градиентом концентрации конкретного хемокина и популяциями клеток, на которые он влияет в образце, корреляция не устанавливается. Сочетание цитокиновых модулей, зависимостей «доза-ответ» хемокинов и нейтрофилов, а также подгрупп пациентов с их собственными воспалительными профилями часто упускается из виду в отсутствие многомерного анализа9.

Развитие вычислительных мощностей делает такой холистический подход все более реалистичным. Сочетание многопараметрической и спектральной проточной цитометрии с методами снижения размерности, такими как t-распределенное стохастическое вложение соседей (t-SNE) или равномерная аппроксимация и проекция многообразий (UMAP), позволяет переносить многомерные иммунные фенотипы в двухмерное пространство, при этом UMAP демонстрирует лучшую сохранность глобальной топологии и воспроизводимость применительно к цитометрическим данным10. Мультианалитные панели цитокинов с последующим иерархическим кластеризатором и построением корреляционных сетей помогают выявить сигнатуры, коррелирующие с тяжестью заболевания, которые невозможно определить с помощью анализа отдельных цитокинов11. Несмотря на прогресс в данной области, до сих пор не было разработано рабочего процесса, объединяющего эти методы для изучения рекрутирования нейтрофилов и цитокиновых сигнатур из одних и тех же образцов БАЛЖ при респираторных инфекциях. В данной работе описывается интегрированный рабочий процесс визуализации для исследования рекрутирования нейтрофилов с помощью цитологии и многопараметрической проточной цитометрии с последующим UMAP-анализом и анализом цитокиновых сигнатур методом мультиплексного иммуноанализа, иерархической кластеризацией, построением корреляционных матриц и анализом главных компонентов. Этот рабочий процесс позволяет создать интегрированную визуализацию, объединяющую клеточные и растворимые иммунные показатели дыхательных путей на одном графике с использованием образцов БАЛЖ, которые рутинно получают во время клинически обоснованной бронхоскопии при инфекциях нижних дыхательных путей.

Протокол

Популяция исследования и этические аспекты

Это было одноцентровое проспективное обсервационное исследование, проведенное в многопрофильной университетской клинике в период с января 2026 г. по апрель 2026 г. Протокол был утвержден под номером 2026HL-038; от каждого участника было получено письменное информированное согласие перед началом любых процедур исследования.

В исследование последовательно включали взрослых пациентов в возрасте от 18 до 75 лет, которым в рамках стандартного клинического обследования проводилась диагностическая гибкая бронхоскопия с бронхоальвеолярным лаважем (БАЛ); участников распределяли в группу с инфекцией или в контрольную группу без инфекции на основании окончательных клинических и микробиологических диагнозов. Обзор интегрированного рабочего процесса визуализации, от набора пациентов до получения итоговых визуальных результатов, представлен на рисунке 1.

figure-protocol-1
Рисунок 1: Обзор интегрированного рабочего процесса визуализации для профилирования рекрутирования нейтрофилов и ответов воспалительных цитокинов на основе БАЛЖ. Пять этапов кратко описывают (слева направо): перспективный набор групп инфекции и контроля с помощью бронхоскопии и БАЛ, предварительную обработку БАЛЖ с разделением на клеточный осадок и супернатантную фракцию, параллельные измерения клеточного состава (цитология по Райту–Гимзе и многопараметрическая проточная цитометрия) и растворимых компонентов (13-плексный иммуноанализ на основе микросфер), единый конвейер визуализации (UMAP с мета-кластеризацией FlowSOM, иерархическая тепловая карта, анализ главных компонентов и сеть корреляции Спирмена), а также интегрированный визуальный результат. Сокращения: BAL = бронхоальвеолярный лаваж; BALF = БАЛЖ (жидкость при БАЛ); UMAP = равномерное приближение и проекция многообразия. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Критерии включения в группу инфицированных были следующими: (1) диагноз внебольничной пневмонии (CAP), госпитальной пневмонии (HAP) или вентилятор-ассоциированной пневмонии (VAP) согласно рекомендациям Общества инфекционных болезней Америки (IDSA) и Американского торакального общества (ATS); (2) обнаружение одного или нескольких патогенов с помощью посева БАЛЖ с титром ≥104 колониеобразующих единиц (CFU)/mL, мультиплексной полимеразной цепной реакции (PCR) респираторных проб, метагеномного секвенирования нового поколения БАЛЖ или специфических тестов на антиген/антитела в парных образцах респираторного тракта и сыворотки крови.

Критерии для контрольной группы были следующими: (1) неинфекционные показания к бронхоскопии, включающие оценку неопределенного узла в легком, аномальные результаты рентгенологического исследования, не указывающие на инфекцию, или обследование по поводу хронического кашля; (2) отсутствие острой легочной инфекции, подтвержденное отрицательными результатами микробиологического тестирования и отсутствием каких-либо признаков воспаления или системной инфекции на момент проведения бронхоскопии; и (3) отсутствие применения антибиотиков в течение 30 дней до участия в исследовании. Контрольная группа была сформирована из пациентов, проходящих бронхоскопию по другим клиническим показаниям, а не из здоровых добровольцев, что отражает популяцию, у которой в данном исследовательском центре обычно проводят BAL по клиническим показаниям.

Критерии исключения включали: (1) злокачественные новообразования и/или проведение цитотоксической химиотерапии в течение предшествующих шести месяцев; (2) серопозитивность по ВИЧ; (3) трансплантация солидных органов или гемопоэтических стволовых клеток; (4) длительная иммуносупрессивная терапия с ежедневной эквивалентной дозой преднизолона > 0.3 mg/kg на протяжении более 3 недель; (5) беременность или период лактации; (6) геморрагический диатез в анамнезе или любые противопоказания к гибкой бронхоскопии; (7) имеющиеся тяжелые хронические заболевания дыхательных путей (хроническая обструктивная болезнь легких [ХОБЛ] стадии III/IV по классификации Global Initiative for Chronic Obstructive Lung Disease [GOLD], идиопатический легочный фиброз, тяжелая персистирующая астма или клинически значимый бронхоэктаз); (8) активный туберкулез; (9) участие в другом интервенционном исследовании в течение последних 30 дней. БАЛ проводили во время первой бронхоскопии, в течение 72 h с момента поступления при наличии инфекций, а у пациентов на искусственной вентиляции легких с инфекциями — в течение 48 h после интубации.

Расчет размера выборки был сосредоточен на основном клеточном показателе: количестве нейтрофилов в БАЛЖ в разных группах. С учетом ожидаемой величины эффекта 0,6 SD между инфицированными и неинфицированными дыхательными путями и с учетом 15% допустимого выбывания участников, для обеспечения двустороннего α = 0,05 при статистической мощности 80% потребовался минимальный размер выборки 45 субъектов в группе. Из 156 подходящих кандидатов 27 не соответствовали критериям отбора (18 не соответствовали критериям включения, главным образом из-за невозможности определить патогенный агент, ответственный за заболевание, в установленный срок, или из-за несогласия с самой процедурой; 9 соответствовали критериям исключения, прежде всего из-за активного злокачественного новообразования или наличия в анамнезе химиотерапии, хронического применения кортикостероидов, тяжелого предшествующего хронического заболевания легких или геморрагического диатеза). Остальные 129 участников прошли бронхоскопию. Впоследствии из исследования были исключены 17 субъектов (7 из-за недостаточного выхода БАЛЖ < 40%; 4 из-за отзыва согласия во время госпитализации; и 6 из-за переклассификации субъектов после окончательной микробиологической оценки — 4 кандидата в контрольную группу с диагностированной скрытой инфекцией по результатам посева; и 2 кандидата с инфекцией, у которых после анализа посева не было обнаружено доказательств инфекции). Итоговый пул участников исследования составил 112 человек: 64 с инфекцией и 48 без нее.

Сбор и предварительная обработка БАЛЖ

Бронхоскопию проводили опытные операторы по стандартному протоколу учреждения, в соответствии с опубликованными техническими рекомендациями по бронхоальвеолярному лаважу12. Лидокаин применяли местно; пациентам проводили умеренную внутривенную седацию, при этом подавали дополнительный кислород для поддержания сатурации кислорода (SpO₂) >92%. Гибкий бронхоскоп вводили и позиционировали в бронхопульмонарный сегмент с радиологическими изменениями. В случаях, когда радиологические изменения были диффузными или нелокализованными, выбирали среднюю долю правого легкого или лингвилярные сегменты левого легкого. После каждого введения жидкости выполняли осторожное ручное аспирационное отсасывание с отрицательным давлением < 100 mmHg. Вводили от трех до пяти порций стерильного физраствора при 37 °C объемом 20–40 mL каждая, чтобы общий объем составил 100–150 mL. Коэффициент восстановления оценивали как отношение объема полученного конденсата к объему введенной жидкости; для анализа учитывали только образцы с восстановлением ≥40%.

Образцы, собранные в результате второго и последующих промываний, объединяли для снижения уровня бронхиального загрязнения и получения образца, обогащенного альвеолярными клетками. Объединенный пул жидкости для бронхоальвеолярного лаважа (БАЛЖ) хранили на льду и обрабатывали в течение 2 ч после сбора, учитывая имеющиеся данные об изменениях клеточного состава и дифференциального подсчета после длительного хранения и нагрева образцов БАЛЖ человека13. После фильтрования через стерильную нейлоновую сетку с размером пор 70 микрон для удаления слизи (фильтрацию проводили путем механического пропускания без добавления муколитиков или антиагглютинационных агентов, таких как дитиотреитол, так как предварительные испытания показали, что этого достаточно для удаления слизистых конгломератов без заметной потери целостности поверхностных эпитопов), жидкость центрифугировали при 400 × g в течение 10 мин при 4 °C. Супернатант сливали и распределяли по криопробиркам с низким уровнем связывания белков для анализа цитокинов, а клетки суспендировали в фосфатно-солевом буферном растворе (PBS) с 0,5% бычьим сывороточным альбумином и 2 mM этилендиаминтетрауксусной кислотой (EDTA), часть которых оставляли для цитологического исследования и проточной цитометрии. Образцы цитокинов замораживали в жидком азоте, хранили при –80 °C и подвергали не более чем одному циклу замораживания-оттаивания перед использованием.

Цитология нейтрофилов и количественное определение методом проточной цитометрии

Для цитопатологического исследования 100 µL клеточной суспензии BALF наносили на предметное стекло с помощью цитоцентрифуги (см. Таблицу материалов) при 300 × g в течение 5 мин, высушивали на воздухе и окрашивали по методу Райта-Гимзы. Дифференциальный подсчет проводил специалист по гематологии, не имевший информации об анамнезе пациента; при этом на каждом стекле подсчитывали не менее 400 клеток и рассчитывали процентное содержание макрофагов, лимфоцитов, нейтрофилов и эозинофилов. Процент нейтрофилов в BALF считали основным показателем цитопатологического анализа, поскольку имеющиеся данные свидетельствуют об их полезности для дифференциации инфекционных и неинфекционных заболеваний легких у пациентов с легочными инфильтратами14.

Протокол количественного анализа методом проточной цитометрии был основан на недавней валидации протокола дифференцировки лейкоцитов жидкости бронхоальвеолярного лаважа (БАЛ) с использованием четырех цветов15 и включал дополнительные маркеры для идентификации нейтрофилов. В частности, примерно 5 × 105 клеток, выделенных из БАЛ, блокировали с помощью реагента для блокирования Fc-рецепторов человека (см. Таблицу материалов) в течение 10 мин при 4 °C, а затем окрашивали в течение 20 мин при 4 °C антителами к CD45 (клон HI30), CD66b (клон G10F5), CD15 (клон HI98), CD16 (клон 3G8), CD11b (клон ICRF44), CD14 (клон M5E2) и HLA-DR (клон L243), а также фиксирующим красителем для определения жизнеспособности (см. Таблицу материалов). Для лизиса эритроцитов использовали коммерческий раствор для лизиса эритроцитов (см. Таблицу материалов) в течение 10 мин при комнатной температуре. После лизиса эритроцитов клетки дважды промывали PBS с добавлением 2% фетальной телячьей сыворотки и ресуспендировали в буфере для проточной цитометрии, после чего анализировали с помощью многоцветного проточного цитометра (см. Таблицу материалов). Компенсацию проводили с использованием компенсационных бусин с одним красителем (см. Таблицу материалов); для каждой панели включали пробы с одним отсутствующим красителем (FMO) и неокрашенные пробы.

Для каждого образца было зарегистрировано не менее 1 × 105 событий CD45⁺. Сначала на графике FSC-A/FSC-H были выделены синглеты, затем живые (viability dye⁻) CD45⁺ лейкоциты, после чего было проведено исключение CD14⁺HLA-DR⁺ моноцитов-макрофагов для отбора CD66b⁺ гранулоцитов с высоким боковым рассеянием. Затем нейтрофилы в популяции гранулоцитов определялись как события CD15⁺CD16⁺CD11b⁺. Это первичное окно для CD66b⁺ гранулоцитов было проведено широко вокруг всех событий с высоким боковым рассеянием, а не ограничивалось границей CD15⁺CD16⁺CD11b⁺, чтобы CD16-dim, CD11b-high активированные нейтрофилы и другие менее зрелые формы гранулоцитов оставались в общем автоматизированном подсчете нейтрофилов; более строгое определение CD15⁺CD16⁺CD11b⁺ применялось только для характеристики подмножества зрелых нейтрофилов, представленного отдельно в последующем кластерном анализе. Количество нейтрофилов в БАЛЖ (cells/µL) оценивали на основе процента нейтрофилов и общего количества ядерных клеток, определенного с помощью автоматического гематологического анализатора (см. Таблицу материалов). Матрицы компенсации, стратегия гейтирования и индивидуальные отчеты о контроле качества образцов были заархивированы, а аналитики оставались слепыми к распределению участников по группам на протяжении всего процесса сбора и анализа данных.

Мультиплексный анализ профиля воспалительных цитокинов

Супернатанты БАЛВ размораживали один раз на влажном льду, после чего подвергали кратковременному перемешиванию в вортексе и осветлению путем центрифугирования при 10,000 × g в течение 5 min при 4 °C непосредственно перед проведением анализа. Измерение 13 аналитов трех функциональных категорий осуществляли с помощью мультиплексной иммуноаналитической платформы на основе микросфер с использованием специализированной панели магнитных микросфер (см. Таблицу материалов). Провоспалительная группа включала IL-1β, IL-6, IL-17A и TNF-α; группа хемокинов включала CXCL8/IL-8, CXCL1/GRO-α, CXCL10/IP-10, CCL2/MCP-1 и CCL3/MIP-1α; регуляторная группа включала IL-10, IFN-γ, IL-1RA и G-CSF. Выбор аналитов основывался на данных, которые чаще всего описывались как информативные в мультиплексных анализах цитокинов БАЛВ или сыворотки крови в исследованиях пневмонии и тяжелых респираторных инфекций7,11. Коэффициент разведения (1:2 в аналитическом буфере) определяли в ходе пилотного эксперимента на репрезентативной подгруппе как инфицированных, так и контрольных образцов, чтобы обеспечить попадание >90% измерений в диапазон стандартной кривой.

Стандартные кривые определяли с использованием семиточечных серийных разведений рекомбинантных стандартов производителя, приготовленных в дубликатах на каждой планшете, включая одну холостую лунку для вычитания фона. Аппроксимация кривой выполнялась методом пятипараметрической логистической (5PL) регрессии с помощью специализированного программного обеспечения прибора для сбора и анализа данных (см. Таблица материалов). Нижний предел количественного определения (НПКО) рассчитывали как минимальное количество стандарта, обеспечивающее степень восстановления при обратном расчете от 80% до 120% и внутрисерийный коэффициент вариации (КВ) ≤ 20%. Значения, оказавшиеся ниже НПКО, заменяли значением НПКО/√2. Два объединенных образца жидкости для бронхоальвеолярного лаважа (БАЛ) для контроля качества, полученных от репрезентативного инфицированного пациента и здорового субъекта соответственно, были разделены на аликвоты и добавлены на каждую планшетку для отслеживания долгосрочной воспроизводимости результатов. Стандарты обеспечения качества основывались на рекомендациях по передовым методам, выработанных в ходе 12-летней многоцентровой программы по проверке квалификации в проведении мультиплексного анализа на основе магнитных частиц.16включая требования к коэффициенту вариации (CV) внутри анализа ≤ 10% для дубликатных лунок, CV между планшетами ≤ 20% для контроля качества и минимальному количеству микросфер ≥ 50 на один аналит в лунке. Планшеты, не соответствующие вышеуказанным критериям, были переанализированы. Анализ проводили одним батчем в течение четырех недель, используя одну серию наборов, с рандомизированным чередованием инфицированных образцов и контрольных проб.

Конвейер анализа визуализации

Файлы стандарта проточной цитометрии (FCS) сначала анализировали с помощью коммерческого программного обеспечения для проточной цитометрии (см. Таблицу материалов) для компенсации и выделения гейтом живых синглетов CD45⁺, после чего экспортировали для объективного анализа в R v4.4. Для каждого образца проводили даунсэмплинг и конкатенацию 10 000 событий CD45⁺, а интенсивность сигналов трансформировали с использованием функции arcsinh с кофактором 150. Снижение размерности и проекции UMAP выполняли с помощью пакета uwot при 15 ближайших соседях и минимальном расстоянии 0,1, а также параллельное вложение t-SNE с использованием пакета Rtsne с перплексией 3010. Фенотипические кластеры определяли с помощью FlowSOM со 100 узлами самоорганизующейся карты, которые затем разделяли на 12–15 кластеров посредством мета-кластеризации, используя основанный на оптимизации подход, который был недавно валидирован для многомерных данных цитометрии17. Каждый кластер аннотировали на основе медианной экспрессии маркеров линий и активации, после чего сверяли с результатами ручного выделения гейтами.

Для анализа данных мультиплексного определения цитокинов исходные значения сначала подвергали log₁₀-преобразованию, а затем стандартизировали между аналитами с помощью z-оценки перед проведением дальнейшего анализа. Для аналитов и пациентов выполняли иерархическую кластеризацию с использованием евклидова расстояния и метода связывания Ward.D2; результаты были представлены в виде теплокарты с кластеризацией по строкам и столбцам. Такой подход применялся для выявления коэкспрессируемых модулей цитокинов и классификации подгрупп пациентов в мультиплексных воспалительных когортах9. Матрицы корреляции для всех клеточных и растворимых биомаркеров были построены с использованием попарных ранговых корреляций Спирмена. Теплокарты матриц корреляции были представлены вместе с сетью, в которой узлы соответствуют биомаркерам, а ребра указывают на корреляции, сохранившиеся после поправки на множественное тестирование методом Бенджамини — Хохберга (BH). Ширина ребер масштабировалась в соответствии с величиной коэффициента корреляции |ρ|, и ребра включались только при |ρ| ≥ 0.30 и p < 0.05 после поправки. Анализ главных компонентов (PCA) стандартизированной матрицы, содержащей данные по клеткам и цитокинам, выполняли с помощью функции prcomp() с масштабированием единичной дисперсии. Участники были спроецированы на две первые главные компоненты с наложением групп инфекции и контроля. Последовательный даунсэмплинг и стандартизация всех видов анализа обеспечили сопоставимость показателей клеточных и растворимых биомаркеров на всех визуализациях.

Статистический анализ

Непрерывные переменные при отсутствии нормального распределения по результатам теста Шапиро–Уилка описывали с помощью медианы (межквартильный размах, МКР), в остальных случаях — с помощью среднего значения (стандартное отклонение, СО). Категориальные переменные описывали с помощью частот и долей. Сравнение групп проводили с использованием критерия Манна–Уитни U тест, когда непрерывные переменные имели распределение, отличное от нормального, тест Уэлча т—тест применялся при нормальном распределении непрерывных переменных, а критерий хи-квадрат или точный тест Фишера — в соответствующих случаях для категориальных переменных. Связи между непрерывными переменными оценивали с помощью коэффициента ранговой корреляции Спирмена. Для учета частоты ложноположительных результатов (FDR) при мультиплексном анализе цитокинов и корреляционном анализе значения p корректировали методом Бенджамини–Хохберга с пороговым значением q < 0,05. Двусторонний критерий p < Значение 0,05 считалось статистически значимым во всех основных анализах. По возможности случаи с отсутствующими данными по отдельным переменным сохранялись и анализировались методом анализа полных случаев для соответствующей переменной, в то время как для проведения анализа чувствительности применялось многократное импотирование с использованием цепочек уравнений (m = 20). Для отдельных цитокинов в рамках всей когорты исследования с помощью пакета pROC были построены рабочие характеристики приемника (ROC-кривые) и рассчитана соответствующая площадь под кривой (AUC); отдельная когорта для внешней валидации не использовалась. В рамках предварительно определенных анализов чувствительности основные сравнения групп были повторены после исключения участников, принимавших антибиотики в течение 72 часов до бронхоскопии, а также после стратификации группы инфекции по статусу искусственной вентиляции легких и по подтипу пневмонии (внебольничная, госпитальная или вентилятор-ассоциированная). Анализы проводились в R v4.4 с использованием пакетов rstatix, FlowSOM, uwot, ComplexHeatmap, igraph, pROC и mice. Рисунки были созданы с помощью коммерческого программного обеспечения для построения научных графиков (см. Таблица материалов).

Результаты

Популяция исследования и характеристики образцов

В окончательную когорту анализа вошли 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).

figure-results-1
Рисунок 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; статус механической вентиляции легких не изменил направление ни одного межгруппового сравнения.

figure-results-2
Рисунок 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, корреляционная матрица и сетевой граф выявили согласованную закономерность, при которой инфильтрация нейтрофилов в БАЛ и ответ хемокинов тесно взаимосвязаны на уровне отдельных участников, а провоспалительные и регуляторные пути функционируют как четко определенные, но взаимосвязанные модули вокруг этой оси «нейтрофилы — хемокины».

figure-results-3
Рисунок 4: Совместная визуализация рекрутирования нейтрофилов и профилей воспалительных цитокинов. Биплот анализа главных компонентов стандартизированной матрицы клеточных и цитокиновых данных с цветовой разметкой участников по группам, тепловая карта корреляции Спирмена и сетевой граф всех 17 клеточных и растворимых признаков (размер узлов соответствует степени связности, ребра масштабированы по |ρ| и отфильтрованы по BH-скорректированному p < 0.05 и |ρ| ≥ 0.30), выделяющий хаб нейтрофилы–хемокины. Сокращения: BH = Бенджамини–Хохберг; PC = главная компонента. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

ДОСТУПНОСТЬ ДАННЫХ:

Обезличенный набор данных на уровне участников и код анализа на языке R, использованные в данном исследовании, представлены в качестве Дополнительная таблица S1 и Дополнительный файл 1Другие данные, подтверждающие результаты данного исследования, могут быть предоставлены автором, ответственным за переписку, по обоснованному запросу.

Дополнительная таблица S1: Набор данных BALF. Пожалуйста, нажмите здесь, чтобы скачать этот файл.

Дополнительный файл 1: ZIP-архив, содержащий коды для анализа. Пожалуйста, нажмите здесь, чтобы скачать этот файл.

Обсуждение

Надежность данного протокола зависит от ряда решений, принятых до и во время анализа. Установление нижней границы степени восстановления на уровне 40% и объединение второй и последующих аликвот были преднамеренными решениями, так как первая аликвота содержит более высокую концентрацию материалов дыхательных путей бронхов, в то время как степень восстановления систематически варьирует в зависимости от возраста, пола, истории курения и целевого бронха18. Блокирование Fc-рецепторов и титрование антител против отдельных клонов для обеспечения правильных границ гейтирования между зрелыми CD15⁺CD16⁺CD11b⁺ нейтрофилами и CD14⁺HLA-DR⁺ альвеолярными макрофагами (доминирующими клетками в образце БАЛЖ с наиболее высоким уровнем неспецифического окрашивания и фоновой автофлуоресценции) имеют жизненно важное значение для поддержания согласованности между двумя группами19,20. Ограничение по количеству циклов замораживания-размораживания и 2-часовой временной интервал для мультиплексного анализа также не являются случайными. Данные по БАЛЖ человека показали значительное снижение процента нейтрофилов в период от 6 до 24 ч даже при хранении при 4 °C, в то время как процент макрофагов и лимфоцитов оставался стабильным13. Параметры алгоритмов t-SNE и UMAP, а также параметры сетки FlowSOM из 100 узлов и 12–15 метакластеров были определены до анализа групп, чтобы предотвратить подстройку параметров post hoc, что является одной из известных причин невоспроизводимости протоколов цитометрии16. Сведения об операторах и идентификаторы планшетов/запусков могут быть сохранены для каждого клеточного события с целью применения методов коррекции пакетного эффекта в будущих крупномасштабных анализах21.

Существует несколько распространенных сценариев неудач, которые не требуют внесения изменений в протокол. Низкий выход (<40%) обычно обусловлен неправильным позиционированием клина или чрезмерным всасыванием; эту проблему следует устранить, скорректировав положение бронхоскопа и снизив интенсивность всасывания до <100 mmHg перед аспирацией следующей аликвоты. Избыточную аутофлуоресценцию альвеолярных макрофагов, являющуюся основным мешающим фактором при анализе клеток легких методом проточной цитометрии, можно минимизировать, включая в каждую серию контроли из неокрашенных клеток легких и четко определяя гейт макрофагов по высокому боковому рассеянию, а не полагаясь исключительно на красители жизнеспособности19. Мультиплексные панели часто подвержены влиянию низкого количества микросфер и аномально высоких значений. Выполнение предварительной кривой разведений для подмножества инфицированных и контрольных образцов перед выбором разведения 1:2 позволяет избежать насыщения сигнала CXCL8, CXCL1 и G-CSF, которые для образцов при пневмонии находятся в верхней части стандартной кривой вследствие ранее обнаруженного матричного эффекта в образцах БАЛЖ22. Протокол может быть дополнительно скорректирован в зависимости от особенностей когорты. Если образцы преимущественно содержат один тип патогена, сокращение основной панели до хемокинов и провоспалительных медиаторов позволит сохранить достаточную статистическую мощность. В случае смешанной популяции, содержащей бактерии, вирусы и атипичные патогены, будет полезно измерить также регуляторные медиаторы, так как недавно было доказано существование специфического для патогена паттерна в альвеолярных нейтрофилах и макрофагах23, а комбинация про- и регуляторных цитокинов в БАЛЖ, как было показано, отражает тяжесть течения внебольничной пневмонии24.

Вклад данного рабочего процесса заключается не в отдельном наблюдении за притоком нейтрофилов или повышением уровня хемокинов и провоспалительных цитокинов при респираторной инфекции, что подробно описано в литературе, а в сочетании панели проточной цитометрии для нейтрофилов и моноцитов из семи маркеров и мультиплексного иммуноанализа из 13 аналитов в рамках общего конвейера визуализации, применяемого к одному и тому же образцу БАЛЖ. Внедрение панели из семи маркеров, нацеленной на нейтрофилы и моноциты, вместе с мультиплексным иммуноанализом из 13 аналитов позволяет зафиксировать фенотипическую и цитокиновую сложность; без этого анализ потребовал бы проведения нескольких независимых тестов с гораздо большим объемом БАЛЖ, чем можно получить при стандартной бронхоскопии. Первый метод хорошо согласуется с предыдущими работами, показывающими, что многоцветная проточная цитометрия столь же надежна, как и ручная микроскопия, для идентификации типов клеток в образцах БАЛЖ15, в то время как критерии приемлемости для второго метода отражают передовой опыт признанной многоцентровой программы оценки компетенций по мультиплексному анализу на основе микросфер16. В частности, комбинированная тепловая карта, сетевой график и PCA, построенные с использованием идентичных стандартизированных входных данных, выявляют взаимосвязи между типами клеток и связанными с ними факторами, которые ранее можно было только косвенно предположить, суммируя независимые столбчатые диаграммы. Это важно, поскольку хемокины CXCL8, CCL2 и CCL3, обнаруживаемые в БАЛЖ, совместно работают для привлечения нейтрофилов и моноцитов из общего циркуляторного русла, и их совместная структура несет информацию, которую невозможно извлечь из независимых показателей отдельных аналитов9,10. БАЛЖ уже рутинно получают в клинических целях, и дополнительный объем не требуется, поэтому данный рабочий процесс может быть позиционирован как трансляционное дополнение к существующим услугам бронхоскопии, а не как отдельная диагностическая процедура.

Практическое время выполнения всего рабочего процесса, включающего транспортировку образцов, окрашивание, регистрацию данных методом проточной цитометрии, инкубацию и считывание мультиплексного анализа, а также последующий вычислительный конвейер, составляет порядка 24 h в условиях исследовательской лаборатории. Таким образом, данный рабочий процесс позиционируется как инструмент глубокого фенотипирования для механистической характеристики и мониторинга ответа на лечение, а не как экспресс-метод прикроватной сортировки для принятия критических по времени решений, таких как немедленное различие между вентилятор-ассоциированной пневмонией и стерильным острым респираторным дистресс-синдромом, при которых стандартом медицинской помощи остаются более быстрые микробиологические анализы и тесты на единичные биомаркеры.

Необходимо упомянуть о нескольких ограничениях. БАЛЖ собирается из дыхательных путей и отражает непосредственные поверхностные взаимодействия в воздушном пространстве, однако он не содержит клеток, расположенных в альвеолярной стенке и перегородках, где локализуется множество активированных макрофагов и резидентных лимфоцитов ткани. Следовательно, между иммунным статусом воздушного пространства и паренхимы нет однозначного соответствия, что согласуется с ранее опубликованными данными25. Мультиплексное определение на основе микросфер подвержено перекрестной реактивности антител, а также матриксным эффектам, которые особенно выражены в среде БАЛЖ, содержащей муцин и сурфактант, по сравнению с сывороткой крови; кроме того, процесс разведения эпителиальной выстилающей жидкости физиологическим раствором невозможно компенсировать с помощью каких-либо эндогенных маркеров, таких как мочевина, что согласуется с недавно опубликованными критическими обзорами26. Объем полученного промывного содержимого был в среднем на 6 mL ниже в группе инфекции, чем в контрольной группе, поэтому сообщаемые здесь абсолютные концентрации подвержены небольшой разнице в коэффициенте разведения между группами, и эндогенный маркер разведения не применялся. Предварительный прием антибиотиков 41 из 64 участников группы инфекции (64,1%) может привести к снижению концентрации цитокинов в дыхательных путях и затруднить выделение патогена. Важно отметить, что протокол представляет собой лишь «мгновенный снимок» в одной временной точке во время клинически показанной бронхоскопии, что исключает анализ кинетики привлечения нейтрофилов, максимальной секреции хемокинов и их разрешения у отдельных пациентов. Значения AUC для одного аналита были получены из той же когорты, которая использовалась для идентификации биомаркеров, при отсутствии независимой валидационной когорты, и должны интерпретироваться как показатели дискриминации внутри когорты, а не как внешне валидированная диагностическая эффективность; внедрение любого порогового значения для одного маркера в клиническую практику потребует проспективной валидации в независимой многоцентровой когорте. Заранее определенные анализы чувствительности, ограниченные образцами с более высоким выходом, участниками, не получавшими антибиотиков, а также стратифицированные по статусу искусственной вентиляции легких и по подтипу пневмонии (внебольничная, нозокомиальная или вентилятор-ассоциированная), не изменили направление межгрупповых сравнений.

В будущем эти ограничения можно будет устранить, добавив секвенирование РНК отдельных клеток или CITE-seq к тем же осадкам БАЛЖ, которые использовались для количественного определения белков и клеток; данные методы были успешно применены в крупных когортах пациентов с бактериальной пневмонией и педиатрической пневмонией, вызванной Mycoplasma pneumoniae, для выявления подмножеств альвеолярных нейтрофилов и макрофагов, которые невозможно различить с помощью традиционной проточной цитометрии27,28. Дополнение информации, полученной при анализе воздушных пространств, данными масс-цитометрии изображений или пространственной транскриптомики соответствующих биопсий или даже материалов аутопсии позволило бы контекстуализировать эти данные29. Связывание лонгитюдного отбора проб БАЛЖ с оценкой действия антагонистов CXCR1/CXCR2 и блокады рецепторов IL-6, которые продемонстрировали преимущество в выживаемости или необходимости органосистемной поддержки при тяжелой пневмонии и остром респираторном дистресс-синдроме (ОРДС), также позволило бы протестировать ось «нейтрофилы-хемокины», выдвинутую на передний план в данном рабочем процессе, в качестве средства принятия терапевтических решений30,31.

В заключение, в данной работе описывается комплексный рабочий процесс визуализации, состоящий из стандартизированных процедур сбора и подготовки БАЛЖ, панели проточной цитометрии с семью маркерами для количественного определения состава нейтрофилов и моноцитов, мультиплексного анализа цитокинов для измерения 13 цитокинов, а также общей платформы визуализации, включающей UMAP, иерархическую кластеризацию, корреляционную сеть и PCA-анализ. При применении этого рабочего процесса к проспективной когорте из 112 взрослых пациентов, которым была показана бронхоскопия по клиническим показаниям, анализ выявил резкий переход в составе иммунных клеток альвеол от макрофагов к нейтрофилам у пациентов с инфекцией, обнаружил трехмодульную структуру цитокинового ответа и выявил сильно связанную ось нейтрофил–хемокин, которая не прослеживается при использовании измерений по одному маркеру. Данная процедура опирается исключительно на материал, собранный в ходе клинической диагностики, не требует дополнительного объема образца и может быть реализована в лабораториях третичного звена, оснащенных оборудованием для многопараметрической проточной цитометрии и мультиплексного анализа на основе микросфер. Более широкое внедрение и любое использование отдельных аналитов в качестве диагностических биомаркеров потребуют проспективной валидации в независимых многоцентровых когортах. Таким образом, исследование представляет воспроизводимую и насыщенную данными методологию для характеристики воспаления дыхательных путей при респираторных инфекциях и создает основу для ее дальнейшего развития в лонгитюдных исследованиях, одноклеточном и пространственном анализе, а также для терапевтической оценки иммуномодуляции, направленной на ось нейтрофил–хемокин.

Раскрытие информации

Авторы заявляют об отсутствии конфликта интересов.

Благодарности

Мы искренне благодарим всех пациентов и добровольцев, принявших участие в данном исследовании, за их неоценимое сотрудничество. Мы выражаем признательность Комитету по этике за этическое одобрение (номер одобрения 2026HL-038) и нормативное руководство на протяжении всего процесса исследования.

Мы благодарим медицинский персонал отделения респираторной медицины, центра бронхоскопии и клинической лаборатории за помощь в сборе клинических образцов, обработке жидкости при бронхоальвеолярном лаваже и микробиологическом обнаружении. Мы также выражаем признательность за техническую поддержку центральной лаборатории больницы при проведении экспериментов по многомерной проточной цитометрии и мультиплексному профилированию цитокинов.

Данное исследование не получило никаких специальных грантов от финансирующих организаций из государственного, коммерческого или некоммерческого секторов.

Мы благодарим всех исследователей, техников и средний медицинский персонал, принимавших участие в этом проекте, за их вклад в скрининг пациентов, управление образцами, проведение экспериментов и систематизацию данных.

Материалы

Список материалов, использованных в этой статье
ИмяКомпанияКаталожный номерКомментарии
Антитело к CD11b человека, Brilliant Violet 421, клон ICRF44BioLegendКат. № 301323Маркер активации нейтрофилов; фиолетовый лазерный канал
URL: https://www.biolegend.com
Антитела к CD14 человека, APC/Cyanine7, клон M5E2BioLegendКат. № 301820Маркер исключения моноцитов-макрофагов; красный лазерный канал
Поскольку вы не предоставили исходный текст для перевода, пожалуйста, вставьте текст или предоставьте ссылку на контент JoVE, который необходимо перевести. Я готов приступить к работе, следуя всем вашим инструкциям по соблюдению академического стиля, научной точности и терминологического единства. https://www.biolegend.com
Антитело к CD15 человека, PE, клон HI98BioLegendКат. № 301906Маркер линии нейтрофилов; синий лазерный канал
Пожалуйста, предоставьте исходный текст для перевода. https://www.biolegend.com
Антитело против CD16 человека, Brilliant Violet 510, клон 3G8BioLegendКат. № 302048Созревание нейтрофилов и маркер NK-клеток; фиолетовый лазерный канал
Поскольку текст в поле URL отсутствует, пожалуйста, предоставьте исходный текст для перевода. https://www.biolegend.com
Антитело к CD45 человека, PerCP, клон HI30BioLegendКат. № 304026 / RRID: AB_893338Панлейкоцитарный маркер; синий лазерный канал
Пожалуйста, предоставьте исходный текст для перевода. https://antibodyregistry.org/AB_893338
Антитело к CD66b человека, FITC, клон G10F5BioLegendКат. № 305103Маркер гранулоцитов; синий лазерный канал
Пожалуйста, предоставьте исходный текст для перевода. https://www.biolegend.com
Антитела к HLA-DR человека, PerCP/Cyanine5.5, клон L243BioLegendКат. № 307629Маркер антигенпрезентирующих клеток; синий лазерный канал
Поскольку вы не предоставили исходный текст для перевода (поле URL пусто), пожалуйста, введите текст или предоставьте ссылку на контент, который необходимо перевести на русский язык. Я готов приступить к работе в соответствии с вашими инструкциями по соблюдению академического стиля и научной точности. https://www.biolegend.com
Автоматический гематологический анализатор (Sysmex XN-1000)Sysmex CorporationПрибор (модель XN-1000)Используется для определения общего количества нуклеированных клеток с целью расчета абсолютного числа нейтрофилов
Пожалуйста, предоставьте текст для перевода. https://www.sysmex.com
Бычий сывороточный альбуминSigma-Aldrich (MilliporeSigma)Кат. № A7906Используется в буфере для окрашивания при проточной цитометрии (0,5% масс./об. в PBS с 2 мМ ЭДТА)
Пожалуйста, предоставьте исходный текст для перевода. https://www.sigmaaldrich.com/US/ru/product/sigma/a7906
Пакет ComplexHeatmap (R/Bioconductor)Проект Bioconductorv2.20 (с открытым исходным кодом)Используется для построения тепловых карт иерархической кластеризации данных мультиплексного анализа цитокинов
Пожалуйста, предоставьте исходный текст для перевода. https://bioconductor.org/packages/ComplexHeatmap
Цитоцентрифуга (Cytospin 4)Thermo Fisher ScientificКат. № A78300003Используется для приготовления цитологических препаратов из клеточной суспензии БАЛЖ (бронхоальвеолярного лаважа)
Пожалуйста, предоставьте исходный текст для перевода. Набор для выделения РНК (RNA Isolation Kit)
Раствор для лизиса эритроцитов (BD FACS Lysing Solution, 10-кратный концентрат)BD BiosciencesКат. № 349202Используется для лизиса эритроцитов после поверхностного окрашивания
Поскольку вы не предоставили исходный текст для перевода, пожалуйста, вставьте текст, который необходимо перевести с английского на русский язык. Я выполню перевод в соответствии с вашими инструкциями, соблюдая академический стиль и научную терминологию. https://www.bdbiosciences.com
Этилендиаминтетрауксусная кислота (ЭДТА), 0,5 M, pH 8,0Thermo Fisher Scientific (Invitrogen)Кат. № AM9260GИспользуется в буфере для окрашивания при проточной цитометрии в концентрации 2 мМ
Пожалуйста, предоставьте текст для перевода. Набор для выделения РНК из тканей и клеток, с колонкой, 50 препарирований
Фиксируемый краситель для определения жизнеспособности (набор Zombie NIR Fixable Viability Kit)BioLegendКат. № 423105Исключение живых и мертвых клеток перед гейтированием; канал красного лазера
Пожалуйста, предоставьте исходный текст для перевода. https://www.biolegend.com
Гибкий видеобронхоскопOlympus Corporation (BF-1TQ290)Прибор (модель BF-1TQ290)Используется для БАЛ при умеренной седации
Пожалуйста, предоставьте текст для перевода. https://www.olympus-global.com
Проточный цитометр, многоцветный (CytoFLEX LX)Beckman Coulter Life SciencesПрибор (модель CytoFLEX LX, шестилазерная конфигурация)Используется для сбора данных при применении панели из семи маркеров для проточной цитометрии нейтрофилов/моноцитов
Поскольку в вашем запросе отсутствует исходный текст для перевода (поле URL пустое), пожалуйста, предоставьте текст или ссылку на материалы JoVE, которые необходимо перевести на русский язык. Я готов приступить к работе, соблюдая все указанные вами требования к академическому стилю, научной терминологии и точности перевода. https://www.beckman.com
Программное обеспечение для анализа данных проточной цитометрии (FlowJo)BD Biosciences (FlowJo LLC)v10.10Используется для компенсации и ручного гейтирования FCS-файлов
Пожалуйста, предоставьте исходный текст для перевода. https://www.flowjo.com
Пакет FlowSOM (R/Bioconductor)Проект Bioconductorv2.12 (с открытым исходным кодом)Используется для кластеризации и метакластеризации данных проточной цитометрии с помощью самоорганизующихся карт
Поскольку вы не предоставили текст для перевода после «URL:», пожалуйста, вставьте исходный английский текст, который необходимо перевести. Я выполню перевод в соответствии с вашими строгими инструкциями для JoVE. https://bioconductor.org/packages/FlowSOM
GraphPad Prism (программное обеспечение для статистического анализа и построения графиков)Программное обеспечение GraphPad (Dotmatics)v10Используется для создания итоговых рисунков
Пожалуйста, предоставьте исходный текст для перевода. https://www.graphpad.com
Реагент для блокирования Fc-рецепторов человека (Human TruStain FcX)BioLegendКат. № 422301 / RRID: AB_2818986Блокировка Fc-рецепторов перед окрашиванием антителами
Пожалуйста, предоставьте текст для перевода. Вы прислали только «URL:», но не указали саму ссылку или исходный текст. Как только вы предоставите материал, я переведу его в соответствии с вашими инструкциями. https://antibodyregistry.org/AB_2818986
пакет igraph (R)CRANv2.0 (с открытым исходным кодом)Используется для построения и визуализации графа корреляционной сети
Пожалуйста, предоставьте исходный текст для перевода. https://cran.r-project.org/package=igraph
Криопробирки с низкой сорбцией белков, 2,0 млThermo Fisher Scientific (Nunc)Кат. № 375418Используется для аликвот супернатанта БАЛЖ для анализа цитокинов
Пожалуйста, предоставьте исходный текст для перевода. https://www.thermofisher.com/order/catalog/product/375418
Платформа для мультиплексного иммуноанализа на основе магнитных частиц Luminex (MAGPIX)Luminex Corporation (DiaSorin)Прибор (модель MAGPIX)Используется для сбора данных мультиплексного анализа цитокинов
Поскольку вы не предоставили исходный текст для перевода, пожалуйста, отправьте текст или ссылку на материал, который необходимо перевести с английского на русский язык. Я готов приступить к работе, строго соблюдая ваши инструкции по научной терминологии и академическому стилю. https://www.luminexcorp.com
Май-Грüраствор красителя по методу Ройнальда-ГимзыSigma-Aldrich (MilliporeSigma)Кат. № GS500Используется для дифференциального подсчета клеток методом цитоспина
Поскольку в вашем запросе отсутствует исходный текст для перевода, пожалуйста, предоставьте текст, который необходимо перевести с английского на русский язык. Я готов приступить к работе, соблюдая все ваши инструкции по научной терминологии, академическому стилю и точности перевода для материалов JoVE. https://www.sigmaaldrich.com/US/ru/product/sigma/gs500
пакет mice (R)CRANv3.16 (с открытым исходным кодом)Используется для множественной импутации с помощью цепочек уравнений в анализе чувствительности
Пожалуйста, предоставьте текст для перевода. Вы указали только «URL:», но не привели саму ссылку или исходный текст. Как только вы отправите содержимое, я приступлю к профессиональному переводу на русский язык в соответствии с вашими инструкциями. https://cran.r-project.org/package=mice
Мультиплексная панель цитокинов/хемокинов на магнитных микросферах (13-плексная, индивидуальная; MILLIPLEX MAP Human Cytokine/Chemokine)Merck / MilliporeSigmaКат. № HCYTOMAG-60KСпециализированная панель для измерения 13 аналитов, описанных в разделе 2.4
Пожалуйста, предоставьте исходный текст для перевода. Вы отправили только «URL:», но не указали саму ссылку или текст. Как только вы предоставите контент, я переведу его в соответствии с вашими инструкциями для JoVE. https://www.merckmillipore.com
сетчатый нейлоновый фильтр для клеток, 70 µmCorning (Falcon)Кат. № 352350Используется для механической фильтрации БАЛЖ с целью удаления слизистых сгустков
Пожалуйста, предоставьте текст для перевода. https://www.corning.com
Фосфатно-солевой буферный раствор (PBS), pH 7,4Thermo Fisher Scientific (Gibco)Кат. № 10010023Базовый буфер для ресуспендирования и окрашивания клеток
URL: https://www.thermofisher.com/order/catalog/product/10010023
пакет pROC (R)CRANv1.18 (с открытым исходным кодом)Используется для построения ROC-кривых и расчета значений AUC
Пожалуйста, предоставьте исходный текст для перевода. https://cran.r-project.org/package=pROC
R (среда статистических вычислений)R Core Teamv4.4 (с открытым исходным кодом)Платформа для статистического анализа и визуализации
Пожалуйста, предоставьте исходный текст для перевода. https://www.r-project.org
пакет rstatix (R)CRANv0.7 (с открытым исходным кодом)Используется для групповых сравнений и корреляционной статистики
Пожалуйста, предоставьте текст для перевода. https://cran.r-project.org/package=rstatix
Пакет Rtsne (R)CRANv0.17 (с открытым исходным кодом)Используется для вложений t-SNE
Пожалуйста, предоставьте исходный текст для перевода. https://cran.r-project.org/package=Rtsne
Компенсационные бусины с одним красителем (компенсационные бусины UltraComp eBeads)Thermo Fisher Scientific (Invitrogen)Кат. № 01-2222-42Используется для одноцветных контролей компенсации
Пожалуйста, предоставьте исходный текст для перевода. Вы отправили только «URL:» без самого контента. Как только вы пришлете текст, я переведу его в соответствии с вашими инструкциями. Набор для выделения РНК из тканей и клеток (RNA Pure), 50 реакций Набор RNA Pure предназначен для высокоэффективного выделения общей РНК из различных типов тканей и культур клеток. Данный метод обеспечивает получение высококачественной РНК, пригодной для последующего анализа с помощью ОТ-ПЦР (RT-PCR), северного блоттинга или синтеза кДНК. Особенности: - Высокая чистота и выход РНК. - Совместимость с широким спектром образцов. - Оптимизированный протокол для минимизации деградации РНК. Состав набора: - Раствор для лизиса. - Раствор для осаждения. - Растворы для промывки. - RNase-free вода. Условия хранения: 4°C.
Раствор хлорида натрия для инъекций, 0,9% (стерильный физиологический раствор)Baxter InternationalКат. № 2B1324Используется для введения препарата при БАЛ (100–150 мл, 20–аликоты по 40 мл при 37 °C)
Пожалуйста, предоставьте текст для перевода. https://www.baxter.com
пакет uwot (R)CRANv0.2 (с открытым исходным кодом)Используется для снижения размерности методом UMAP
Пожалуйста, предоставьте исходный текст для перевода. https://cran.r-project.org/package=uwot
xPONENT (программное обеспечение для сбора данных мультиплексного анализа и аппроксимации кривых)Luminex Corporation (DiaSorin)v4.3Специализированное программное обеспечение прибора для аппроксимации стандартной кривой по модели 5PL
Пожалуйста, предоставьте текст для перевода. https://www.luminexcorp.com

Ссылки

  1. Yu X, et al. Estimating the global and regional burden of lower respiratory infections attributable to leading pathogens and the protective effectiveness of immunization programs. Int J Infect Dis. 2024;149:107268.
  2. Rizo-Téllez SA, Filep JG. Beyond host defense and tissue injury: the emerging role of neutrophils in tissue repair. Am J Physiol Cell Physiol. 2024;326(3):C661-83.
  3. Zhou X, Jin J, Lv T, Song Y. A narrative review: the role of NETs in acute respiratory distress syndrome/acute lung injury. Int J Mol Sci. 2024;25(3):1464.
  4. Haynes ME, Sullivan DP, Muller WA. Neutrophil infiltration and function in the pathogenesis of inflammatory airspace disease. Am J Pathol. 2024;194(5):628-36.
  5. Davidson KR, Ha DM, Schwarz MI, Chan ED. Bronchoalveolar lavage as a diagnostic procedure: a review of known cellular and molecular findings in various lung diseases. J Thorac Dis. 2020;12(9):4991-5019.
  6. Cambier S, Gouwy M, Proost P. The chemokines CXCL8 and CXCL12: molecular and functional properties, role in disease and efforts towards pharmacological intervention. Cell Mol Immunol. 2023;20(3):217-51.
  7. Lin J, et al. Local and systemic cytokine profiles in children with pneumonia-associated lung consolidation. Front Immunol. 2025;16:1546730.
  8. Huang C, et al. Heparin-binding protein in bronchoalveolar lavage fluid as a biomarker for discriminating severe bacterial and viral pneumonia in critically ill children. Mediators Inflamm. 2023;2023:6123911.
  9. Polley DJ, et al. Identification of novel clusters of co-expressing cytokines in a diagnostic cytokine multiplex test. Front Immunol. 2023;14:1223817.
  10. Roca CP, et al. A cross entropy test allows quantitative statistical comparison of t-SNE and UMAP representations. Cell Rep Methods. 2023;3(1):100389.
  11. Hawerkamp HC, et al. Characterisation of the pro-inflammatory cytokine signature in severe COVID-19. Front Immunol. 2023;14:1170012.
  12. Meyer KC, et al. An official American Thoracic Society clinical practice guideline: the clinical utility of bronchoalveolar lavage cellular analysis in interstitial lung disease. Am J Respir Crit Care Med. 2012;185(9):1004-14.
  13. Shionoya Y, et al. Impact of storage conditions on bronchoalveolar lavage fluid analysis: a human study. Diagnostics (Basel). 2025;15(11):1386.
  14. Li J, et al. Diagnostic test accuracy of cellular analysis of bronchoalveolar lavage fluid in distinguishing pulmonary infectious and non-infectious diseases in patients with pulmonary shadow. Front Med (Lausanne). 2024;11:1496088.
  15. Bratke K, Weise M, Stoll P, Virchow JC, Lommatzsch M. Flow cytometry as an alternative to microscopy for the differentiation of BAL fluid leukocytes. Chest. 2024;166(4):793-801.
  16. Rountree W, et al. Sources of variability in Luminex bead-based cytokine assays: evidence from twelve years of multi-site proficiency testing. J Immunol Methods. 2024;531:113699.
  17. Tao W, Sinha A, Raddassi K, Pandit A. Parameter optimization for stable clustering using FlowSOM: a case study from CyTOF. Front Immunol. 2024;15:1414400.
  18. Shikano K, et al. What are the factors affecting the recovery rate of bronchoalveolar lavage fluid? Clin Respir J. 2022;16(2):142-51.
  19. Baumann Z, et al. Optimized full-spectrum flow cytometry panel for deep immunophenotyping of murine lungs. Cell Rep Methods. 2024;4(11):100885.
  20. Bonilla DL, et al. The power of reagent titration in flow cytometry. Cells. 2024;13(20):1677.
  21. Ogishi M, et al. Multibatch cytometry data integration for optimal immunophenotyping. J Immunol. 2021;206(1):206-13.
  22. Sudo K, et al. Case study observational research: inflammatory cytokines in the bronchial epithelial lining fluid of COVID-19 patients with acute hypoxemic respiratory failure. Crit Care. 2024;28(1):134.
  23. Williams JE, et al. Epigenetic regulation of neutrophils in ARDS. Cells. 2025;14(15):1151.
  24. Zhang Y, et al. Interleukin-6 in blood and bronchoalveolar lavage fluid of hospitalized children with community-acquired pneumonia. Front Pediatr. 2022;10:922143.
  25. Bailey JI, et al. Profibrotic monocyte-derived alveolar macrophages are expanded in patients with persistent respiratory symptoms and radiographic abnormalities after COVID-19. Nat Immunol. 2024;25(11):2097-109.
  26. Haeger S, et al. The bronchoalveolar lavage dilution conundrum: an updated view on a long-standing problem. Am J Physiol Lung Cell Mol Physiol. 2024;327(5):L807-13.
  27. Xiao K, et al. A pan-immune panorama of bacterial pneumonia revealed by a large-scale single-cell transcriptome atlas. Signal Transduct Target Ther. 2025;10(1):5.
  28. Shen X, et al. Single-cell transcriptome atlas revealed bronchoalveolar immune features related to disease severity in pediatric Mycoplasma pneumoniae pneumonia. MedComm (2020). 2024;5(10):e748.
  29. Megas S, Wilbrey-Clark A, Maartens A, Teichmann SA, Meyer KB. Spatial transcriptomics of the respiratory system. Annu Rev Physiol. 2025;87(1):447-70.
  30. Landoni G, et al. A multicenter phase 2 randomized controlled study on the efficacy and safety of reparixin in the treatment of hospitalized patients with COVID-19 pneumonia. Infect Dis Ther. 2022;11(4):1559-74.
  31. Sellarès-Nadal J, et al. Efficacy of tocilizumab for hospitalized patients with COVID-19 pneumonia and high IL-6 levels: a randomized controlled trial. Infection. 2025;53(5):1851-61.

Перепечатки и разрешения

Теги

UMAPFlowSOM