Composition of oropharyngeal microbiota in patients with COVID-19 of different pneumonia severity

Full Text

Список сокращений

ГКС – глюкокортикостероиды

ИВЛ – искусственная вентиляция легких

ИПП – ингибиторы протонной помпы

КПУ – индекс количественного показателя кариозных, пломбированных и удаленных зубов

КТ – компьютерная томография

НПВС – нестероидные противовоспалительные средства

ХОБЛ – хроническая обструктивная болезнь легких

ASV (amplicon sequence variants) – последовательность ампликонов

NMDS (Non-Metric Multidimensional Scaling) – неметрическое многомерное шкалирование

PHP – индекс эффективности гигиены

PMA – папиллярно-маргинально-альвеолярный индекс

Введение

Течение COVID-19 зачастую осложняют сопутствующие бактериальные инфекции, которые могут определять клинический исход. Часть вторичных инфекций вызывают микроорганизмы, попадающие в легкие из верхних дыхательных путей. Не меньшую роль играют бактерии ротовой полости: с током слюны они попадают в глотку, откуда при вдыхании могут попасть в нижние дыхательные пути. Известно, что микробиота полости рта является второй по величине уровня α-разнообразия после микробиоты желудочно-кишечного тракта, включая более 700 видов бактерий, более 100 видов грибов и простейших [1]. Имеются сведения о влиянии состояния твердых тканей зубов и слизистой оболочки полости рта на частоту развития аспирационных пневмоний [2].

Особый интерес в свете анализа связи SARS-CoV-2 с бактериями верхних дыхательных путей и ротовой полости представляет ротоглотка, так как именно она чаще всего является входными воротами возбудителей респираторных инфекций, в том числе SARS-CoV-2 [3]. Предполагается, что состав микробиоты полости рта и ротоглотки может оказывать влияние на развитие и течение COVID-19, либо выступая в качестве барьера, тем или иным способом подавляя репликацию вируса, либо являясь источником вторичной бактериальной инфекции, что в свою очередь может влиять на тяжесть течения болезни и клинический исход.

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

Материалы и методы

Клиническая группа. В исследование включены 156 пациентов с лабораторно подтвержденной новой коронавирусной инфекцией SARS-CoV-2 и клиническими проявлениями COVID-19, которые проходили лечение на базе многопрофильного клинического медицинского центра «Кусково» ФГБОУ ВО «МГМСУ им. А.И. Евдокимова» в период с апреля по июнь 2021 г. Возраст пациентов составил от 20 до 82 лет, из них 68 пациентов мужского пола и 88 – женского пола. Госпитализированы 77 пациентов с поражением легких менее 25% по данным компьютерной томографии (КТ) – КТ-1, 66 – с поражением 35–50% (КТ-2) и 13 – с поражением 50–75% (КТ-3), при этом 53 человека нуждались в принудительной кислородотерапии (рис. 1).

 

Рис. 1. Распределение числа пациентов с диагностированным COVID-19 по возрасту, полу, степени поражения легких (KT) и потребности в дополнительной подаче кислорода.

 

Критерии включения в изучаемую группу:

• письменное согласие участника на участие в исследовании;
• подтвержденный тестом с помощью метода полимеразной цепной реакции диагноз COVID-19;
• возраст старше 18 лет.

Критерии исключения:

• тяжелое состояние с угрозой для жизни или противопоказания, препятствующие забору биоматериала;
• онкологические заболевания вне стадии ремиссии.

Также на основе опросников собраны данные о сопутствующих заболеваниях. Врачом-стоматологом проведены осмотры полости рта пациентов с заполнением индекса количественного показателя кариозных, пломбированных и удаленных зубов (КПУ) и папиллярно-маргинально-альвеолярного индекса (PMA) для оценки состояния твердых тканей зубов и слизистой оболочки полости рта. Также проводилась оценка гигиены пациентов с помощью индекса эффективности гигиены (PHP). Согласно данным опросникам 7 пациентов имели диагноз «хроническая обструктивная болезнь легких» (ХОБЛ), 30 – сахарный диабет, 23 – воспалительные заболевания кишечника и 98 – сердечно-сосудистые заболевания.

По стоматологическому статусу 63 пациента имели высокую резистентность эмали зубов к кариесу и здоровый пародонт, 93 имели множественные кариозные поражения зубов и заболевания пародонта. При измерении индекса PMA у 130 человек наблюдался гингивит легкой степени тяжести и у 26 – гингивит средней степени тяжести. Индексом КПУ измеряли интенсивность кариеса: 65 человек с высокой, 80 – со средней и 11 – с низкой интенсивностью кариеса. При оценке гигиены полости рта пациентов использовался индекс PHP, по результатам которого 54 пациента имели хороший индекс гигиены, 83 – средний и 19 – неудовлетворительный.

Кроме того, 78 пациентов принимали антибиотики в течение месяца до госпитализации.

Забор биоматериала осуществляли при помощи ватного тампона с пораженных участков ротоглотки – миндалин, не касаясь слизистой языка и полости рта. При необходимости биоматериал также забирали с дужек мягкого неба, небного язычка и задней стенки глотки. После забора образцов 2 зонда с ватным тампоном помещали в промаркированную криопробирку объемом 1,8 мл и подвергали немедленной заморозке. Хранение биоматериала осуществляли при температуре -70°C.

Метагеномное секвенирование. Выделение тотальной ДНК проводили с помощью набора MagMAX DNA Multi-Sample Ultra 2.0 Kit (Roche) и автоматической выделительной станции KingFisher Flex. Подготовку библиотек и их секвенирование осуществляли в соответствии с протоколом 16S Metagenomic Sequencing Library preparation для MiSeq, Illumina.

Приготовление геномных библиотек проводили набором Tersus PCR kit (Evrogen, Россия), для амплификации вариабельного региона V3–V4 гена 16s рРНК использовали праймеры 341F и 801R, для баркодирования образцов использовали уникальные комбинации индексирующих праймеров, являющихся аналогами праймеров из набора Nextera XT Index kit v2.

Процедуру секвенирования осуществляли с помощью инструмента HiSeq 2500 с использованием набора расходных материалов HiSeq Rapid SBS Kit v2 (500) (Illumina) согласно рекомендациям производителя.

Биоинформатическая обработка. Фильтрацию низкокачественных последовательностей проводили с помощью метода filterAndTrim из пакета DADA2. Для получения вариантов последовательностей ампликонов использовали пакет DADA2 [4].

Полученные последовательности ампликонов таксономически аннотировали методом idTaxa из пакета DECIPHER (версия 11.7) с использованием базы последовательностей рибосомальных РНК Silva (версия 138.1, nr99).

Для идентификации и удаления контаминантных последовательностей использовали пакет Decontam, использующий в качестве критерия обратную зависимость концентрации ДНК образца и представленности последовательностей [5]. Полная схема биоинформатической обработки представлена на рис. 2.

 

Рис. 2. Схема биоинформатической обработки последовательностей ампликонов гена 16s рРНК образцов пациентов.

 

Статистический анализ представленностей последовательностей ампликонов выполняли на языке программирования R (версия 4.0.2).

Для расчета индекса β-разнообразия значения представленности вариантов последовательностей ампликонов (ASV, Amplicon Sequence Variant) нормализовали при помощи функции transform пакета microbiome (версия 1.10.0). Матрицу попарных расстояний рассчитывали по метрике Брея–Кертиса с помощью функции dist пакета phyloseq (версия 1.32.0).

Значение меры α-разнообразия рассчитывали как индекс Шеннона при помощи функции plot_richness пакета phyloseq (версия 1.32.0) [6].

Для определения достоверности попарного различия сообществ по составу микробиоты и вклада метаданных пациентов в данное различие использовали метод непараметрического дисперсионного анализа (функция adonis пакета vegan, версия 2.5-7) с оценкой значимости при 999 перестановках. Анализ дифференциальной представленности для групп сравнения (КТ-1 и КТ≥2) выполнен с использованием инструмента Deseq 2 (версия 1.28.1) [7].

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

Для получения устойчивых, воспроизводимых результатов анализ проводили в 25 итерациях, в каждой из которых случайным образом выбирали 80% образцов. Полученные данные отфильтрованы в соответствии с выбранными нами порогами воспроизводимости: «строгим» (ASV, достоверные различия в представленности которых между группами воспроизводятся минимум 20 раз из 25) и «умеренным» (ASV, достоверные различия в представленности которых между группами воспроизводятся минимум 15 раз из 25). При проверке статистических гипотез минимальный уровень α составил 0,05 (по методу Бенджамини–Хохберга), а размер эффекта log2foldchange – 1.5.

Результаты

Организация сбора клинического материала и метаданных. Исследуемая выборка сформирована из пациентов, поступивших на лечение в многопрофильный клинический медицинский центр «Кусково» с лабораторно подтвержденной новой коронавирусной инфекцией SARS-CoV-2. При соответствии пациента критериям включения/исключения от него забраны мазки из ротоглотки на метагеномное тестирование. Деперсонализированные метаданные внесены в информационную систему, разработанную под задачи настоящего исследования с использованием международных практик и основанную на стандарте медицинских данных HL7 FHIR. Применение электронной системы обеспечило сбор информации о пациентах в едином формате, что является необходимым условием проведения дальнейшего биоинформатического анализа данных.

Метагеномное секвенирование. Собранные у поступивших в стационар пациентов мазки ротоглотки подвергнуты метагеномному секвенированию вариабельного участка V3–V4 гена 16s рРНК. С использованием генетического анализатора Illumina HiSeq 2500 для всей совокупности образцов получены прочтения длиной около 550 нуклеотидов. Среднее число прочтений на образец составило 45465.68.

Таксономический состав микробиоты ротоглотки пациентов с COVID-19. Исследование общего таксономического состава микробиоты ротоглотки пациентов с различной степенью поражения легких по данным КТ выявило 768 ASV (amplicon sequence variants – вариантов последовательностей ампликонов), принадлежащих 97 бактериальным родам (рис. 3). Наиболее распространенными родами являлись Streptococcus, Prevotella, Veillonella, Rothia, Actinomyces и Leptotrichia, что соответствует литературным данным о таксономическом составе микробиоты ротоглотки человека [8].

 

Рис. 3. Тепловая карта представленностей родов прокариот в образцах ротоглотки 156 пациентов с различной степенью поражения легких (обозначено желтым, оранжевым и фиолетовым цветом в верхней панели; КТ). Цвета тепловой карты отражают значения представленностей соответствующих бактериальных родов после преобразования CLR (centered log ratio transform).

 

Оценка β-разнообразия ASV сообществ образцов пациентов с разной формой пневмонии, развившейся на фоне COVID-19, проводилась при помощи неметрического многомерного шкалирования (Non-Metric Multidimensional Scaling – NMDS) по метрике Брея–Кертиса (Bray–Curtis); рис. 4, а. Расстояния между точками на графике соответствуют степени сходства таксономического состава образцов. При оценке β-разнообразия бактериальных сообществ исследуемых групп существенные различия не выявлены.

При оценке значений α-разнообразия микробных сообществ на основании сравнения значений индекса разнообразия Шеннона (рис. 4, b) не выявлено значимых различий между группами пациентов с различной степенью поражения легких (KT-1 и КТ≥2).

 

Рис. 4: а – график сходства таксономического состава (b-разнообразия) метагеномов пациентов с разной степенью поражения легких по KT: KT-1 (n=77), KT-2 (n=66), KT-3 (n=13) на основании NMDS по видовому разнообразию (метрика различия: мера Брея–Кертиса) методом главных компонент; b – таксономическое разнообразие (a-разнообразие) микробиоты пациентов с разной степенью поражения легких по KT: KT-1 (n=77), KT-2 (n=66), KT-3 (n=13), рассчитанное как индекс Шеннона.

 

Анализ вклада метаданных в вариацию дисперсии состава микробиома

Дисперсия состава микробиоты, объясняемая метаданными, рассчитана методом перестановочного многомерного дисперсионного анализа с использованием матриц расстояний, реализованного в функции adonis. При выполнении анализа собранные метаданные объединены в группы, такие как стоматологические показатели: стоматологический статус, индекс PMA, индекс PHP, индекс КПУз; состояние пациентов: степень тяжести пневмонии по КТ, потребность в дополнительном O₂, прием антибиотиков в течение месяца до госпитализации; физиологические показатели: пол, возраст; хронические заболевания сердечно-сосудистой системы, ХОБЛ, сахарный диабет, воспалительные заболевания кишечника (рис. 5). По результатам анализа данные показатели объяснили 12% дисперсии таксономического состава микробиоты, при этом наибольший вклад внесли стоматологические показатели и показатели состояния пациентов (4,5 и 3,1% вариации дисперсии соответственно).

 

Рис. 5. Вклад групп метаданных в вариацию дисперсии таксономического состава микробиоты образцов пациентов. Физиологические показатели: пол, возраст. Хронические заболевания: хронические болезни сердечно-сосудистой системы, ХОБЛ, сахарный диабет, воспалительные заболевания кишечника. Состояние пациентов: степень тяжести пневмонии по данным КТ, потребность в дополнительном O₂, прием антибиотиков в течение месяца до госпитализации. Стоматологические показатели: стоматологический статус, индекс PMA, индекс PHP, индекс КПУз.

 

Различия в таксономическом составе микробиоты ротоглотки пациентов с разной степенью поражения легких на фоне COVID-19

В результате анализа дифференциальной представленности последовательностей ампликонов гена 16s рРНК образцов пациентов с различной степенью поражения легких выявлены последовательности, предположительно ассоциированные с более тяжелой (КТ≥2) или более легкой (КТ-1) формой пневмонии, развившейся вне стационарного лечения (рис. 6).

 

Рис. 6. Анализ дифференциальной представленности 16s ампликонов бактерий, ассоциированных с тяжелой (КТ>2, n=79) и легкой (КТ=1, n=77) степенями поражения легких пациентов. Анализ проводился с использованием метода DeSeq2: a – ASV, воспроизводящиеся в ≥20 итерациях из 25; b – ASV, воспроизводящиеся в ≥15 итерациях из 25. ASV бактерий, которые ассоциированы с высокой степенью поражения легких у пациентов, показаны столбиками, направленными вправо; с легкой степенью поражения легких у пациентов – влево. Справа и слева от столбцов отображен уровень статистической значимости (a), рассчитанный в результате анализа для каждого дифференциально представленного ASV.

 

Результаты оценивались на двух порогах воспроизводимости: «строгом» (ASV, воспроизводящиеся в ≥20 итерациях из 25); см. рис. 6, а; и «умеренном» (ASV, воспроизводящиеся в ≥15 итерациях из 25); см. рис. 6, b.

После фильтрации результатов по «строгому» порогу воспроизводимости выявлено 18 вариантов последовательностей ампликонов (ASV), предположительно ассоциированных с меньшей степенью поражения легких. Они соответствовали таксонам класса Negativicutes (Veillonella), типа Bacteroidetes (Capnocytophaga, Porphyromonas, Alloprevotella и Prevotella), а также классам α- (Lautropia, Neisseria) и γ-протеобактерий (Haemophilus). Также выявлены ASV, относящиеся к роду Actinomyces и к филуму Saccharibacteria (ранее известному как TM7). При фильтрации результатов по «умеренному» порогу воспроизводимости дополнительно выявлено 15 ASV, соответствующих таксонам классов Negativicutes (Anaeroglobus, Selenomonas), γ-протеобактерий (Aggregatibacter) и ε-протеобактерий (Campylobacter), типа Actinobacteria (Atopobium) и порядку Fusobacteria (Leptotrichia).

Последовательности 5 ампликонов, предположительно ассоциированные с более сильным поражением легких, полученные при фильтрации по «строгому» порогу воспроизводимости, относились к грамположительным родам типов Firmicutes (Limosilactobacillus, Anoxybacillus, Mogibacterium) и Actinobacteria (Atopobium, Alloscardovia). При фильтрации результатов по «умеренному» порогу воспроизводимости дополнительно выявлены 5 ASV, соответствующих таксонам типов Firmicutes (Streptococcus, Bulleidia, Howardella) и Actinobacteria (Rothia), а также роду, относящемуся к типу Bacteroidetes (Alloprevotella).

На основании данных результатов можно заключить, что ASV, ассоциированные с меньшей степенью поражения легких, в основном относятся к классу грамотрицательных фирмикут (Negativicutes), к различным классам протеобактерий, а также к порядку Fusobacteria. В свою очередь ASV, ассоциированные с большей степенью поражения легких, относятся преимущественно к грамположительным фирмикутам классов Bacilli и Clostridia. Следует отметить, что разные ASV, отнесенные к родам Atopobium и Alloprevotella, находятся в числе ассоциированных как с легкой, так и с тяжелой формами пневмонии у пациентов с COVID-19.

Клиническое наблюдение за отобранной группой пациентов на фоне терапевтического лечения в условиях стационара

В зависимости от характера течения заболевания в период госпитализации включенные в исследование 156 пациентов разделены на 2 группы. Первую группу (n=70) составили пациенты со стабильным течением заболевания и положительной динамикой в течение стационарного лечения. Вторая группа (n=86) представлена больными с эпизодами отрицательной динамики течения заболевания в период госпитализации (прогрессирование пневмонии по данным КТ, развитие осложнений заболевания, перевод пациентов в отделение реанимации и интенсивной терапии, летальный исход, клинические состояния, требующие усиления фармакотерапии и применения интервенционных методов лечения, включая дополнительную инсуфляцию кислорода и искусственную вентиляцию легких – ИВЛ).

Сравнительный анализ между группами показал, что в терапии пациентов 2-й группы достоверно более часто использовались антибактериальные препараты, ингибиторы протонной помпы (ИПП), глюкокортикостероиды (ГКС), а также ИВЛ и дополнительная инсуфляция кислородом (табл. 1).

 

Таблица 1. Применяемые методы терапии в выделенных группах пациентов Table 1. Methods of therapy in selected groups of patients
Параметр, абс. (%)	1-я группа (без отрицательной динамики, n=70; 44,9%)	2-я группа (с отрицательной динамикой, n=86; 55,1%)	p (критерий Фишера)
Антибактериальные препараты (любые)	11 (15,7)	86 (100)	0,000000
Цефтриаксон	8 (11,4)	58 (67,4)	0,000000
Пенициллин	0	22 (25,6)	0,000000
Фторхинолоны	0	11 (12,8)	0,000000
Макролиды	3 (4,3)	8 (9,3)	0,347137
Антикоагулянты (любые)	64 (91,4)	82 (95,3)	0,345555
Эноксапарин	9 (12,9)	21 (24,4)	0,101270
Гепарин	54 (77,1)	62 (72,1)	0,543720
Аторвастатин	1 (1,4)	3 (3,5)	0,628167
ИПП	43 (61,4)	76 (88,4)	0,000118
ГКС	62 (88,6)	84 (97,7)	0,043502
Блокаторы Са2+-каналов	4 (5,7)	5 (5,8)	0,999999
β-Адреноблокаторы	10 (14,2)	10 (11,6)	0,638431
НПВС	28 (40,0)	37 (43,0)	0,745697
Противовирусные препараты	32 (45,7)	34 (39,5)	0,515036
Муколитики	20 (28,6)	27 (31,4)	0,728881
Ингибиторы ангиотензинпревращаю- щего фермента	15 (21,4)	15 (17,4)	0,546824
Антагонисты рецепторов альдостерона	6 (8,6)	5 (5,8)	0,543482
Гипогликемические средства	3 (4,3)	6 (6,9)	0,732106
ИВЛ	0 (0)	7 (8,1)	0,000000
Потребность в дополнительной инсуфляции кислорода	16 (22,9)	69 (80,2)	0,000000

 

Антибиотикотерапия проводилась в 1-й группе у 11/70 (15,7%) пациентов, муколитики назначены 20/70 (28,6%) пациентам, потребность в дополнительном О₂ выявлена у 16/70 (22,9%) пациентов. Во 2-й группе антибиотики (в том числе комбинированная схема приема) назначены всем 86/86 (100%) пациентам, муколитики – 27/86 (31,4%) пациентам, в дополнительном О₂ нуждались 69/86 (80,2%) пациентов, потребность в подключении к аппарату ИВЛ возникла у 7/86 (8,1%) пациентов. Данные показатели достоверно позволяют считать, что на фоне выраженного ухудшения состояния возникает потребность в купировании отрицательной динамики респираторной системы пациентов.

Назначение антикоагулянтов в 1-й группе отмечено у 64/70 (91,4%) пациентов, а во 2-й – у 82/86 (95,3%) – данные препараты назначаются практически всем пациентам в качестве превентивных мер повышенного тромбообразования и не являются достоверным критерием оценки изменения тяжести состояния. Точно так же с дополнительным назначением ГКС (дексаметазон), где в 1-й группе его назначили 62/70 (88,6%) пациентам и во 2-й группе – 84/86 (97,7%) пациентам – это служит предупредительной мерой развития отечных явлений легких, улучшает дыхательную динамику, тем самым предупреждая развитие осложнений.

Назначение нестероидных противовоспалительных средств (НПВС) связано с потребностью в снижении температуры тела, общим аналгизирующим эффектом на фоне возможного появления мышечных и головных болей. Так, в 1-й группе он назначен 28/70 (40%) пациентам, а во 2-й группе – 37/86 (43%) пациентам. Назначение ИПП можно связать с превенцией риска развития НПВC-гастропатии. В 1-й группе ИПП назначались 43/70 (61,4%) пациентам, а во 2-й – 76/86 (88,4%).

Сопоставление данных клинических наблюдений с первичным распределением пациентов на основании КТ легких (табл. 2) продемонстрировало, что в 1-й группе (КТ-1, n=77) незначительно преобладали пациенты без отрицательной динамики в ходе терапии (46/77; 59,7%), тогда как среди пациентов с более тяжелой степенью поражения легких на момент поступления (КТ-2–3, n=79) ожидаемо доминировали пациенты с эпизодами отрицательной динамики в период нахождения в стационаре (56/79; 70,9%).

 

Таблица 2. Сопоставление тяжести течения COVID-19 с формой пневмонии при поступлении пациентов в стационар Table 2. Comparison of the severity of the course of COVID-19 with the form of pneumonia upon admission of patients to the hospital
	КТ-1 (n=77)	КТ-2, КТ-3 (n=79)	Итого
Без отрицательной динамики (n=70)	46	24	70
С отрицательной динамикой (n=86)	31	55	86
Итого	77	79	156

 

Обсуждение

В данном исследовании мы рассмотрели возможную связь между микробиотой ротоглотки и развитием пневмонии на фоне COVID-19 и выявили ряд таксонов, предположительно ассоциированных с высокой и низкой степенями поражения легких.

В исследование намеренно включены мазки ротоглотки от пациентов с COVID-19, отобранные при их поступлении в стационар: такой подход позволил исключить потенциальное присоединение инфекций, связанных с оказанием медицинской помощи, позволяя рассматривать развившуюся у пациентов пневмонию как внебольничную.

Согласно нашим данным, с высокой степенью поражения легких ассоциирован ряд последовательностей, относящихся к грамположительным фирмикутам классов Bacilli и Clostridia, в числе которых последовательности рода Mogibacterium. Для данного микроорганизма, ранее выделявшегося из образцов пациентов с респираторными инфекциями, в недавнем исследовании показана провоспалительная активность in vitro [9].

С меньшей степенью поражения легких ассоциирован ряд последовательностей, относящихся к грамотрицательным фирмикутам (Negativicutes), протеобактериям и к порядку Fusobacteria (в частности, к роду Leptotrichia).

Представленность как протеобактерий, так и Leptotrichiа и Haemophilus у пациентов с меньшей степенью поражения легких можно предположительно связать с сиаловыми кислотами, находящимися на поверхностях клеток данных бактерий. Сиаловые кислоты поверхностей эпителиальных клеток человека являются дополнительной мишенью для связывания вирусных частиц SARS-CoV-2 помимо ACE2 рецептора [10]. Протеобактерии и Fusobacteria, к которым относится род Leptotrichia, способны к синтезу сиаловых кислот de novo, в то время как бактерии рода Haemophilus, которые также ассоциированы с меньшей степенью поражения легких, катаболизируют сиаловые кислоты с поверхности клеток человека [11]. Ранее высказывалось предположение, что бактерии верхних дыхательных путей, имеющие на своей поверхности сиаловые кислоты, могут выполнять барьерную функцию, являясь конкурентной мишенью для связывания SARS-CoV-2, что приводит к снижению вирусной нагрузки на нижние дыхательные пути [12]. Это может быть одним из объяснений обнаруженной нами ассоциации бактерий родов Leptotrichia и Haemophilus с более легким течением болезни.

О сниженной представленности рода Leptotrichia в носоглотке больных COVID-19 уже неоднократно сообщалось ранее [13, 14], при этом для данной бактерии показана положительная корреляция с метиловым эфиром хлорогеновой кислоты [14]. Исследователи предполагали, что данное соединение, предположительно производимое комменсальными микроорганизмами рода Leptotrichia, может играть роль в защите организма от проявлений COVID-19 [14]. В одном из предшествующих исследований для данного метаболита показана противовоспалительная активность на животной модели воспаления легких, осуществляющаяся за счет регуляции сигнального пути COX-2/NLRP3/NF-κB [15].

Несколько ASV, ассоциированных с меньшей степенью поражения легких, относятся к недавно выделенному филуму эпибионтов Saccharibacteria (также известному как TM7), а также к актинобактериям, в том числе Actinomyces и Schaalia, предположительно являющихся базибионтами (хозяевами) данных бактерий. Роль бактерий филума TM7 в организме человека до конца не определена. Ранее их ассоциировали с дисбиотическими микробиомами, выделяемыми у пациентов с периодонтитом [16]. Однако в недавнем исследовании in vivo показано, что TM7 способны регулировать патогенность актинобактерий-хозяев, тем самым снижая вызванный их активностью воспалительный ответ [17]. Это согласуется с данными другого исследования in vitro, в котором ассоциированные с хозяином TM7x способны снижать в макрофагах экспрессию фактора некроза опухоли α, вызванную активностью бактерии-хозяина [18]. На основе этого наблюдения высказывалось предположение, что TM7 способны регулировать иммунный ответ в организме человека [18, 19]. На настоящий момент понимание роли TM7 в микробиоте дыхательных путей и ротовой полости затруднено, в том числе ввиду ограниченности информации о спектре предположительных хозяев данных эпибионтов.

Ретроспективный анализ клинических данных показал, что среди пациентов, имевших легкую форму пневмонии при поступлении в стационар (группа КТ-1, n=76), в большинстве случаев (46/76; 60,5%) за время пребывания в стационаре не наблюдалось отрицательной динамики, тогда как в группе пациентов с более тяжелыми формами пневмонии (группа КТ-2–3, n=80) таких случаев значительно меньше (24/80; 30%). Однако обращает на себя внимание тот факт, что в обеих группах встречались пациенты с отрицательной динамикой во время пребывания в стационаре, что еще раз подчеркивает сложность прогнозирования клинического исхода COVID-19. В последующих исследованиях предстоит выяснить, насколько наблюдаемые отличия в картине течения болезни обусловлены составом микробиоты ротоглотки.

Заключение

На основании как полученных нами результатов, так и литературных данных можно предположить два механизма, посредством которых микробиота ротоглотки может оказывать влияние на течение вирусных респираторных инфекций, в частности вызванной SARS-CoV-2: образование «барьера», позволяющего задерживать вирусные частицы и таким образом снижать вирусную нагрузку на нижние дыхательные пути, а также модуляция иммунного ответа за счет регуляции экспрессии как про-, так и противовоспалительных цитокинов.

Раскрытие интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.

Disclosure of interest. The authors declare that they have no competing interests.

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

Authors’ contribution. The authors declare the compliance of their authorship according to the international ICMJE criteria. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work.

Источник финансирования. Эта работа была поддержана грантом РНФ №21-15-00431 и ГЗ №ЕГИСУ НИОКТР 122040100011-0.

About the authors

Elizaveta V. Starikova

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-6582-210X

науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Julia S. Galeeva

Research Institute of Systemic Biology and Medicine

Author for correspondence.
Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-6304-4607

мл. науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Dmitry N. Andreev

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-4007-7112

канд. мед. наук, доц., доц. каф. пропедевтики внутренних болезней и гастроэнтерологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Philipp S. Sokolov

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0003-2813-6498

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19, ст. лаборант каф. ЮНЕСКО «Здоровый образ жизни – залог успешного развития» ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Dmitry E. Fedorov

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-8468-7011

мл. науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Aleksander I. Manolov

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0003-3912-429X

науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Alexander V. Pavlenko

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-9549-0289

науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Ksenia M. Klimina

Federal Scientific and Clinical Center for Physical and Chemical Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-5563-644X

канд. биол. наук, ст. науч. сотр. ФГБУ ФНКЦ ФХМ

Russian Federation, Moscow

Vladimir A. Veselovsky

Federal Scientific and Clinical Center for Physical and Chemical Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-4336-9452

мл. науч. сотр. ФГБУ ФНКЦ ФХМ

Russian Federation, Moscow

Andrew V. Zaborovsky

Yevdokimov Moscow State University of Medicine and Dentistry, Moscow

Email: olgagaleeva546@gmail.com

D. Sci. (Med.)

Russian Federation, Moscow

Vladimir V. Evdokimov

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-9281-579X

д-р мед. наук, проф., нач. управления науки ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Nikolai G. Andreev

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-5136-0140

канд. мед. наук, доц., доц. каф.пропедевтики внутренних болезней и гастроэнтерологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Mikhail K. Devkota

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-3736-4196

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19, преподаватель каф. фармакологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Aleksei K. Fomenko

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-1794-7263

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19, преподаватель каф. фармакологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Vadim A. Khar'kovskii

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-8659-3502

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19 ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Pavel O. Asadulin

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-5236-1770

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19 ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Sergey A. Kucher

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-7981-1786

врач-специалист отд-ния микробиологического анализа Клинического центра COVID-19 ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Aleksandra S. Cheremushkina

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-1089-4322

студентка лечебного фак-та ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Oleg O. Yanushevich

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-4293-8465

акад. РАН, д-р мед. наук, проф., ректор, зав. каф. пародонтологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Igor V. Maev

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0001-6114-564X

акад. РАН, д-р мед. наук, проф., зав. каф. пропедевтики внутренних болезней и гастроэнтерологии лечебного фак-та ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Natella I. Krikheli

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0002-8035-0638

д-р мед. наук, проф., зав. каф. клинической стоматологии ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Oleg V. Levchenko

Yevdokimov Moscow State University of Medicine and Dentistry

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0003-0857-9398

д-р мед. наук, проф. каф. нейрохирургии и нейрореанимации ФГБОУ ВО «МГМСУ им. А.И. Евдокимова»

Russian Federation, Moscow

Elena N. Ilina

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0003-0130-5079

чл.-кор. РАН, проф. РАН, д-р биол. наук, гл. науч. сотр. ФБУН НИИ СБМ

Russian Federation, Moscow

Vadim M. Govorun

Research Institute of Systemic Biology and Medicine

Email: olgagaleeva546@gmail.com
ORCID iD: 0000-0003-0837-8764

акад. РАН, д-р биол. наук, проф., врио дир. ФБУН НИИ СБМ

Russian Federation, Moscow

References

Supplementary files