CC BY
4
ПОЛИМОРФИЗМЫ RS17713054 И RS1800629 ГЕНОВ LZTFL1 И TNFАССОЦИИРОВАНЫ С ТЯЖЕСТЬЮ ТЕЧЕНИЯ COVID-19
А. А. Траспов1, М. М. Минашкин1, С. В. Поярков1, А. Г. Комаров2, И. А. Штинова2, Г. И. Спешилов2, И. А. Карбышев2, Н. В. Позднякова1 М. А. Годков3
1 «Система-БиоТех», Москва, Россия
2 Диагностический центр (Центр лабораторных исследований) Департамента здравоохранения города Москвы, Москва, Россия
3 Научно-исследовательский институт скорой помощи имени Н. В. Склифосовского
Генетические и негенетические факторы ответственны за высокую межиндивидуальную вариабельность ответа на SARS-CoV-2. Хотя многочисленные генетические полиморфизмы были идентифицированы как факторы риска тяжелого течения COVID-19, они остаются недостаточно изученными в российской популяции. Целью данного исследования было выявить генетические детерминанты, ассоциированные с тяжелым течением COVID-19 на выборке пациентов из Российской Федерации. Проведена оценка связи генетических полиморфизмов rs17713054 в гене LZTFL1 и rs1800629 в гене TNF (фактор некроза опухоли) с тяжестью COVID-19. Были исследованы образцы ДНК 713 пациентов (324 мужчины и 389 женщины) возрастом 18-95 лет с COVID-19, протекавшей с разной степенью тяжести. Идентифицированы rs1800629 TNF (OR = 1,5; p = 0,02) и rs17713054 LZTFL1 (OR = 1,60; p = 0,0043) как факторы риска тяжелого течения. Полиморфизмы TNF rs1800629 и LZTFL1 rs17713054 могут быть потенциальным предиктивным и предиктивным биомаркерами. Показана сильная ассоциация rs17713054 G > A с тяжелым течением. Полученные данные могут в дальнейшем стать основой для разработки тест-систем предсказания рисков тяжести течения вирусных заболеваний респираторного тракта.
Ключевые слова: COVID-19, SARS-CoV-2, SNP, LZTFL1, TNF, тяжесть течения, генетические детерминанты
Финансирование: исследование профинансировано главным акционером компании ООО «Система-БиоТех», АФК «Система». Исследование профинансировано в рамках двойного слепого клинического исследования Департаментом здравоохранения г. Москвы
Благодарности: авторы благодарят заместителя руководителя Департамента здравоохранения г. Москвы Ю. О. Антипову за организацию исследования; за предоставление образцов биологических материалов пациентов ГБУЗ «Диагностический центр (Центр лабораторных исследований) г. Москвы».
Вклад авторов: Н. В. Позднякова, С. В. Поярков — концепция и дизайн исследования; М. М. Минашкин — выполнение молекулярно-генетических работ, лабораторная постановка; С. В. Поярков, А. А. Траспов — обзор литературы, написание статьи; А. А. Траспов — статистическая обработка и редактирование статьи; А. Г. Комаров — руководство, планирование эксперимента, анализ данных; И. А. Штинова, Г. И. Спешилов и И. А. Карбышев — предоставление клинического материала и метаданных; М. А. Годков — контроль за сбором образцов, предоставление обезличенных данных.
Соблюдение этических стандартов: одобрение этического комитета не требовалось в связи с работой с обезличенными данными пациентов и двойным слепым форматом исследования.
[>3 Для корреспонденции: Наталья Вячеславовна Позднякова
1-я Курьяновская ул., д. 34, кор.11, Москва, Россия; [email protected]
Статья получена: 18.11.2022 Статья принята к печати: 14.12.2022 Опубликована онлайн: 28.12.2022
DOI: 10.24075/vrgmu.2022.065
THE RS17713054 AND RS1800629 POLYMORPHISMS OF GENES LZTFL1 AND TNF ARE ASSOCIATED WITH COVID-19 SEVERITY
Traspov AA1, Minashkin MM1, Poyarkov SV1, Komarov AG2, Shtinova IA2, Speshilov GI2, Karbyshev IA2, Pozdnyakova NV1H, Godkov MA3
1 Sistema BioTech LLC, Moscow, Russia
2 Diagnostic Center (Center of Laboratory Testing) of the Moscow Department of Health, Moscow, Russia
3 Sklifosovsky Research Institute for Emergency Medicine, Moscow, Russia
Both genetic and non-genetic factors are responsible for high interindividual variability in response to SARS-CoV-2. Despite the fact that multiple genetic polymorphisms have been identified as risk factors of severe COVID-19, such polymorphisms are still insufficiently studied in the Russian population. The study was aimed to identify genetic determinants associated with severe COVID-19 in the sample of patients from the Russian Federation. The correlation of the rs17713054 polymorphism in gene LZTFL1 and rs1800629 polymorphism in gene TNF (tumor necrosis factor) with the COVID-19 severity was assessed. DNA samples obtained from 713 patients (324 males and 389 females) aged 18-95 with COVID-19 of varying severity were analyzed. The rs1800629 polymorphism of gene TNF (OR = 1.5; p = 0.02) and rs17713054 polymorphism of gene LZTFL1 (OR = 1.60; p = 0.0043) were identified as risk factors of severe disease. The TNF polymorphism rs1800629 and LZTFL1 polymorphism rs17713054 could be considered as potential predictive biomarkers. The rs17713054 G > A polymorphism was strongly associated with severe disease. In the future the findings may provide the basis for the development of test-systems for prediction of the risk of severe viral respiratory diseases.
Keywords: COVID-19, SARS-CoV-2, SNP, LZTFL1, TNF, severity, genetic determinants
Funding: the study was supported by the major Sistema BioTech LLC shareholder, the Sistema Public Joint Stock Financial Corporation. The study was funded by the Moscow Department of Health as part of the double blind clinical trial.
Acknowledgements: the authors express their gratitude to Antipova YuO, Deputy Head of the Moscow Department of Health, for organizing the blind clinical trial of COVID-19 severity and the patients of the Diagnostic Center (Center of Laboratory Testing) of the Moscow Department of Health for the provided biomaterial samples.
Author contribution: Pozdnyakova NV, Poyarkov SV — study concept and design; Minashkin MM — molecular genetic research, laboratory tests; Poyarkov SV, Traspov AA — literature review, manuscript writing; Traspov AA — statistical processing of the results, manuscript editing; Komarov AG — project management, planning the experiment, data analysis; Shtinova IA, Speshilov GI, Karbyshev IA — providing clinical data and metadata; Godkov MA — control over clinical sample collection at the Sklifosovsky Research Institute for Emergency Medicine, providing anonymized data.
Compliance with ethical standards: the study did not need to be approved by the Ethics Committee because of dealing with anonymized patient data and the double blind study format.
gg Correspondence should be addressed: Natalia V Pozdnyakova
1 Kurjanovskaja, 34, korp. 11, Moscow, Russia; [email protected]
Received: 18.11.2022 Accepted: 14.12.2022 Published online: 28.12.2022
DOI: 10.24075/brsmu.2022.065
За последние десятилетия было показано, что вариации генома человека способствуют возникновению гетерогенности в ответ на инфекционные заболевания.
Пандемия тяжелого острого респираторного синдрома, вызванного коронавирусом SARS-CoV-2, привела к гибели миллионов людей во всем мире.
Для заболевания, вызванного новой коронавирусной инфекцией (COVID-19), характерна высокая вариабельность клинических проявлений (Всемирная организация здравоохранения, 2021 г.). Большинство пациентов, а именно 81%, имеют бессимптомное или легкое течение, тогда как 14% страдают тяжелым и 5% — критическим течением заболевания [1]. Наиболее распространенные симптомы включают лихорадку, сухой кашель и утомляемость; наблюдаются также агевзия, аносмия и желудочно-кишечные симптомы. Для тяжелой формы COVID-19 характерна дыхательная недостаточность, требующая искусственной вентиляции легких или интенсивной подачи кислорода.
Идет постоянный поиск факторов риска, которые могут быть связаны с болезнью. Изначально предполагалось, что люди пожилого возраста, мужского пола, имевшие ранее сердечно-сосудистые заболевания, заболевают в тяжелой форме. Но оказалось, что генетический бэкграунд вносит не меньший вклад. В результате многочисленных исследований было показано, что полиморфизмы в генах факторов организма-хозяина, необходимых для реализации жизненного цикла вируса, например ACE2, TMPRSS2, AR, а также в генах, задействованных во врожденном иммунном ответе, например TNF, TLR7, могут быть связаны с тяжестью COVID-19 [2].
Первое полногеномное ассоциативное исследование (GWAS) тяжести течения COVID-19, в котором сравнивались 1980 пациентов с тяжелым течением из Италии и Испании с контрольной популяцией с неизвестным статусом инфекции SARS-CoV-2, выявил два значимых для всего генома локуса, картируемых в области 3p21.31, включающей шесть генов (SLC6A26, LZTFL1, CXCR6, CCR1, CCR3, CCR9) и в области 9q34.2, которая содержит локус группы крови ABO [3].
Сигнал в локусе 3p21.31 остается наиболее устойчивым и сильным во множестве исследований, что связано как с восприимчивостью к инфекции, так и с тяжестью заболевания. Аллель C варианта rs10490770 связан с самым сильным увеличением риска тяжелого протекания пневмонии COVID-19 в данном локусе [4].
В полногеномно-ассоциативных исследованиях (Genome-Wide Association Studies, GWAS) было выявлено, что область 3p21.31 связана с двукратным увеличением риска появления дыхательной недостаточности [5].
Роль факторов хозяина, определяющих как чувствительность так и тяжесть течения инфекции, показана для многих патогенов, включая микобактерии туберкулеза, HIV, Candida albicans и многих других [6].
Известна ассоциация риска тяжелого течения инфекции с генами, вовлеченными в иммунный и воспалительный ответы. Так, известна роль минорных аллелей в TLR7 [7] и в генах системы интерферона [8].
Помимо аллелей риска, существуют варианты, которые обладают защитным эффектом. Так, генотип ТТ SNP rs5443 в гене GNB3 ассоциирован с защитой от летального исхода при COVID-19 [9]. Полиморфизмы в регуляторных регионах генов провоспалительных цитокинов, которые могут влиять на уровень экспрессии мРНК этих генов, например C-572G rs1800796 в промоторном регионе гена
IL6, ассоциирован с защитным эффектом от развития тяжелого течения и смерти от COVID-19 в популяции народов Азии [10].
ПАЦИЕНТЫ И МЕТОДЫ
Клинико-демографические характеристики обследованных пациентов.
В исследование были включены пациенты отделения реанимации и интенсивной терапии (ОРИТ) НИИ скорой помощи им. Н. В. Склифосовского, ГКБ № 1 ООО «Медси», ГКБ № 40, а также сотрудники АФК «Система» и лаборатории «Система-БиоТех». Было проведено ретроспективное исследование (исход заболевания был известен), в которое включили 713 пациентов (n = 713). Критерии включения: коронавирусная инфекция различной степени тяжести. Критерии исключения: наличие различных типов коморбидности, способных резко повлиять на общее состояние пациента (онкология, ИБС, дефекты иммунитета).
Для выявления ассоциации исследуемых вариантов с тяжелым и крайне тяжелым течением пациенты были разделены на две группы: группу контроля, куда были включены две категории (пациенты с легким и средне-тяжелым течением), и опыта (пациенты с тяжелым течением) в соответствии со шкалой визуальной оценки легких по результатам компьютерной томографии (или МСКТ), разработанной во время пандемии новой коронавирусной инфекции специалистами из Центра диагностики и телемедицины США, после изучения КТ-исследований 13 003 человек, которые составили основную выборку. Данная классификация получила широкое распространение и в российской практике [24]. Критерием степени тяжести служили данные компьютерной томографии (КТ) легких: для контрольной группы — КТ1 и 2, для тяжелой — КТ3 и 4. Основные признаки исследуемой когорты представлены в табл. 1.
Выделение геномной ДНК
В качестве биологического материала для выделения геномной ДНК использовали венозную кровь, стабилизированную ЭДТА.
ДНК выделяли с помощью набора для выделения ДНК из цельной крови DiaGene («Диа-М»; Россия). Чистоту выделенной НК проверяли на спектрофотометре NanoDrop OneC (Themo FS; США). Соотношение A260/280 колебалось от 1,8 до 1,91, а соотношение A260/230 — от 1,62 до 2,28. Концентрацию ДНК измеряли с помощью набора dsDNA BR на флуориметре Qubit Flex (Thermo FS; США). Ее значения менялись в пределах 15-300 нг/мкл. Концентрацию всех образцов ДНК доводили до 2 нг/мкл.
Генотипирование ДНК
Генотипирование ДНК в области исследуемых полиморфных маркеров проводили методом ПЦР в реальном времени с флуоресцентной детекцией с использованием наборов для генотипического анализа TaqMan SNP (Thermo FS; США). ПЦР проводили на амплификаторе CFX96 (BioRad; США).
Генотипировали ряд маркеров, который в дальнейшем подтверждали секвенированием по методу Сэнгера. Праймеры были разработаны с использованием PrimerBlast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/ (по состоянию на 24 сентября 2020 г.)). В реакции
Таблица 1. Основные параметры возраста и пола, взятые для исследования. Классифицировано на две группы по степени тяжести (группа «контроль» и группа «тяжелые»), по полу и возрасту
n = 713
Группа «контроль» Группа «тяжелые»
Всего 373 (52,31%) 340 (47,69%)
МУЖЧИНЫ 154 (41%) 170 (50%)
ЖЕНЩИНЫ 219 (59%) 170 (50%)
Возраст 21-94 18-95
циклического секвенирования BigDye Terminator v3.1 использовали набор для циклического секвенирования BigDye (Thermo FS; США); секвенирование проводили на Genetic Analyzer 3500 (Thermo FS; США).
Статистическая обработка полученных данных
Для расчета статистической значимости различий частот встречаемости использовали точный двусторонний критерий Фишера и критерий х2 (хи-квадрат). Расчеты проводили как для минорного аллеля (АА против Аа+аа; доминантная модель), так и для гомозиготных генотипов по минорному аллелю (АА+Аа против аа; рецессивная модель). Для групп, в которых были обнаружены значимые различия, рассчитывали отношение шансов (OR), 95%-й доверительный интервал (CI) и критерий статистической значимости. Расчеты проводились с использованием ресурса SNPStats (https://www.snpstats.net/ (по состоянию на 20 сентября 2021 г.)), которое предназначено для выявления связи между однонуклеотидными полиморфизмами и риском заболевания. Для оценки влияния множественных сравнений для значений OR статистически значимых полиморфизмов была проведена поправка Холма-Бонферрони. Были проанализированы выборки, включая все полиморфизмы, на достаточность размера выборки. Анализ проводили в пакете Statsoft Statistica 12. Чтобы проанализировать требуемый размер выборки при заданной силе 80, мы использовали независимый f-тест с 2 выборками для сравнения двух средних значений между двумя популяциями (нулевая гипотеза р1 = р2). Учитывали особенности пакета SNPstats и используемые в нем величины: средние с критериями Акаике (AIC) и Байесовским информационным критерием (Bayesian information criterion, BIC).
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ
В табл. 1 описана исследуемая выборка с целью оценки ее разнообразия.
Для расчета статистической значимости разницы в частотах встречаемости использовали точный двухсторонний критерий Фишера и х2-критерий. Статистические расчеты проводили в соответствии с доминантной моделью (с идентификацией аллеля риска) и рецессивной моделью (с идентификацией генотипа риска). Кроме того, для групп, в которых были обнаружены значимые различия, рассчитывали отношение шансов (OR), с 95%-м доверительным интервалом и критерием статистической значимости р.
Для всех полиморфизмов был сделан (в зависимости от разделения пациентов на группы) расчет статистической значимости с использованием двух моделей: доминантной и рецессивной. Модель базировалась на оценке риска перехода легкой формы в более тяжелую; частоту сравнивали у пациентов контрольной группы, с одной стороны, и у пациентов с тяжелой формой, с другой (табл. 2).
Поскольку в эксперименте проверяли большое количество гипотез на одном и том же наборе данных, из-за эффекта множественных сравнений была внесена поправка Холма-Бонферрони, чтобы избежать ошибок I типа.
Для двух маркеров наличие минорной аллели и гомозиготного генотипа было связано с повышенным риском тяжелого течения заболевания. Это rs1800629 в гене TNF (OR = 1,5; р = 0,02), и rs17713054 в гене LZTFL1 (OR = 1,60; р = 0,0043), однако rs17713054 показал патогенные свойства и в рецессивной модели (OR = 4,56; р = 0,0025).
Таблица 2. БИР, используемые в модели оценки риска тяжести течения в двух группах (контрольная и опытная)
Ген/SNP Доминантная модель (AA против Aa+aa) Рецессивная модель (AA+Aa против aa)
OR (95% CI) ДИ P-value OR (95% CI) ДИ P-value
TN/F/rs1800629 1,5 1,06-2,12 0,02 1,1 0,35-3,44 0,87
IFH1/rs1990760 1,2 0,88-1,63 0,26 1,4 0,98-2,01 0,07
IFITM3/rs12252 0,87 0,55-1,38 0,56 1,66 0,46-5,92 0,43
CCR2/rs1799864 1,27 0,91-1,78 0,16 1,33 0,57-3,11 0,51
STAT3/rs744166 0,94 0,70-1,27 0,71 0,53 0,32-0,88 0,01
STAT6/rs324011 0,93 0,68-1,26 0,62 1,17 0,75-1,81 0,48
TLR2/rs1898830 0,72 0,53-0,98 0,04 0,73 0,48-1,11 0,13
C3AR1/rs7842 1,13 0,84-1,52 0,41 1,37 0,77-2,42 0,28
TMPRSS2/rs75603675 1,23 0,90-1,68 0,19 0,85 0,57-1,27 0,43
LZTFL1/rs17713054 1,6 1,16-2,21 0,0043 4,56 1,51-13,76 0,0025
Примечание: значения ДИ присутствуют только для отношений шансов.
Таблица 3. Расчет частот встречаемости 10 изученных маркеров по доминантной и рецессивной моделям (Allele frequency), частот генотипов HOM1/HET/ HOM2 (Genotype distribution), а также достоверность отклонения по Харди-Вайнбергу (HWE p-value). Все исследование проводили на экспериментальной (n = 340), контрольной (n = 373) и объединенной группах (n =713)
Группа «контроль» в сравнении с группой «тяжелые»
Частота аллеля Распределение генотипов Отклонение по Харди-Вайнбергу, P-value
Ген/Аллельные варианты все контроль опыт все контроль опыт все контроль опыт
TNF (G/A) 0,87/0,13 0,89/0,11 0,85/0,15 0,02/0,22/0,76 0,02/0,18/0,8 0,02/0,26/0,73 0,86 0,42 0,83
IFIH1 (T/C) 0,57/0,43 0,59/0,41 0,54/0,46 0,21/0,44/0,35 0,18/0,45/0,37 0,24/0,43/0,33 0,0047 0,13 0,016
IFITM3 (A/G) 0,93/0,07 0,93/ 0,07 0,94/0,06 0,88/0,11/0,01 0,87/0,12/0,01 0,89/0,09/0,02 0,00053 0,081 0,0012
CCR2 (G/A) 0,85/0,15 0,87/0,13 0,84/0,16 0,03/0,23/0,74 0,03/0,21/0,76 0,04/0,25/0,71 0,049 0,12 0,22
STAT3 (A/G) 0,66/0,34 0,64/0,36 0,68/0,32 0,42/0,47/0,11 0,42/0,45/0,13 0,43/0,49/0,08 0,24 0,65 0,019
STAT6 (C/T) 0,62/ 0,38 0,62/0,38 0,62/0,38 0,37/0,5/0,13 0,36/0,51/0,12 0,38/0,48/0,14 0,13 0,078 0,82
TLR2 (A/G) 0,59/0,41 0,57/0,43 0,62/0,38 0,34/0,52/0,15 0,3/0,53/0,17 0,37/0,5/0,13 0,062 0,11 0,25
C3AR1 (T/C) 0,72/0,28 0,73/ 0,27 0,71/0,29 0,07/0,42/0,51 0,06/0,41/0,53 0,08/0,42/0,5 0,46 0,43 0,9
TMPRSS2 (C/A) 0,59/0,41 0,6/0,4 0,58/ 0,42 0,16/0,5/0,34 0,17/ 0,46/0,36 0,15/0,53/0,32 0,54 0,52 0,094
LZTFL1 (C/T) 0,84/0,16 0,87/0,13 0,8/0,2 0,7/0,27/0,03 0,75/0,24/0,01 0,65/0,3/0,05 0,68 0,37 0,39
В маркере TLR2 защитными свойствами обладал гомозиготный генотип по минорному аллелю в доминантной модели.
Все выявленные маркеры являются факторами риска как для восприимчивости к заболеванию, так и для перехода заболевания в более тяжелую форму, но статистическая значимость была достигнута только для вариантов, проиллюстрированных в таблицах выше.
Кроме того, был также проведен анализ исследуемой выборки на отклонение от распределения равновесия Харди-Вайнберга (п = 713), для 10 исследованных маркеров (табл. 3).
На основе суммирования коэффициентов отклонения от равновесия Харди-Вайнберга и их уровней значимости P было выявлено отклонение от равновесия для IFIH1 p = 0,0047), ССЯ2 (p = 0,049) и значительное — для IFITM3 (p = 0,00053).
Таблица 4. Расчет уровней статистической значимости в различиях между
Таким образом, варианты rs1800629 в гене TNF и rs17713054 в гене LZTFL1 могут быть рассмотрены как потенциальные кандидаты для дальнейшего анализа на расширенных выборках (табл. 4).
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
В нашем исследовании мы демонстрируем подтверждающие доказательства участия локуса 3p21.31 человека в патофизиологии заболевания COVID-19 с использованием независимой когорты пациентов и контрольной группы населения Российской Федерации.
Ранее нами были проанализированы распределение 10 SNP и их ассоциация с тяжестью течения [11]. В настоящей работе, проведенной на большей выборке и с расширенным набором полиморфизмов, исключен полиморфизм в TLR7 и добавлен полиморфизм LZTFL1 и
тпами опыта и контроля в 10 исследованных полиморфизмах
Ген/SNP Аллель Опыт Контроль P-value Отношения шансов ДИ
TNF/rs1800629 G 665 581 0,03801 1,398585 1,009831-1,941277
A 81 99
/F/H1/rs1990760 T 442 370 0,06881 1,217995 0,9818877-1,5111842
C 304 310
IFITM3/rs12252 A 695 636 0,8319 0,9427985 0,6060591-1,4617779
G 51 44
CCR2/rs1799864 G 647 571 0,1536 1,247354 0,9192725-1,6940554
A 99 109
S7AT3/rs744166 A 478 460 0,1627 0,8531148 0,6806029-1,0687811
G 268 220
STAT6/rs1898830 C 464 423 1 0,999678 0,8020877-1,2457747
T 282 257
TLR2/rs1898830 A 423 424 0,03094 0,7908369 0,6357947-0,9831671
G 323 256
C3AR1/rs7842 T 547 481 0,2877 1,137109 0,8959515-1,4432459
C 199 199
TMPRSS2/ rs75603675 C 445 397 0,628 1,053809 0,8482685-1,3092215
A 301 283
LZTF1/rs17713054 C 649 546 0,0007003 1,641442 1,223191-2,208574
T 97 134
обнаружено, что всего два полиморфизма в генах LZTFL1 и TNF показывают сильную ассоциацию с тяжелым течением COVID-19.
Среди проанализированных полиморфизмов только два убедительно ассоциациированы с тяжестью течения. Интересно отметить, что после расширения выборки до 713 пациентов результаты для части полиморфизмов из раннего исследования [11] оказались статистически незначимыми. Тем не менее сильный сигнал был обнаружен в rs17713054 (ген LZTFL1, локус 3p21.3). Вторым по значимости был rs1800629 в гене TNF. Известно, что оба полиморфизма ассоциированы с тяжелым течением COVID-19 и носители минорных аллелей имеют повышенный риск развития тяжелого течения COVID-19 [12]. Полиморфизм rs17713054 располагается в энхансере, который регулирует экспрессию генов в этом локусе, включая LZTFL1 и SLC6A20, а также гены хемокинов.
Этот полиморфизм приводит к появлению нового сайта связывания для транскрипционного фактора C/EBP beta, что приводит к повышению экспрессии LZTFL1 и соседних генов в этом локусе [13]. Ген LZTFL1 участвует в цилиарной функции клеток эпителия легких, что важно для очистки дыхательных путей от вирусов [14].
LZTFL1 широко экспрессируется в эпителиальных клетках легких, включая реснитчатые эпителиальные клетки, которые были идентифицированы как одна из основных клеточных мишеней для инфекции SARS-CoV-2. Кроме того, гомозиготная потеря LZTFL1 вызывает классическую цилиопатию: синдром Барде-Бидля [15]; известно, что респираторные вирусы могут воздействовать на мукоцилиарный клиренс. LZTFL1 кодирует цитозольный белок с мотивом лейциновой застежки, который связывается c E-кадгерином (эпителиальным маркером) и участвует в переносе многочисленных сигнальных молекул. Известно также, что активация LZTFL1 в контексте злокачественных новообразований ингибирует путь эпителиально-мезенхимального перехода (EMT), который, как известно, является частью как механизма заживления ран, так и иммунного ответа [16]. Исследование посмертных биоптатов легких пациентов, умерших от осложнений COVID-19, демонстрирует широко распространенную эпителиальную дисфункцию с признаками EMT [13]. Сигнал в локусе 3p21.31 остается наиболее устойчивым и сильным по данным нескольких исследований, что связано как с восприимчивостью к инфекции, так и с тяжестью заболевания [17].
Показано, что аллель риска rs17713054 A в энхансере гена LZTFL1, в значительной степени ответственен за двукратное увеличение риска дыхательной недостаточности от COVID-19, связанного с 3p21.31 [13]. LZTFL1 широко экспрессируется в эпителиальных клетках легких, включая реснитчатые эпителиальные клетки, которые были идентифицированы как одна из основных клеточных мишеней для инфекции SARS-CoV-2 [18].
Анализ других полиморфизмов гена LZTFL1, таких как rs11385942, показал повышенный риск госпитализации р < 0,01; OR = 5,73; 95% CI: 1,2-26,5 по аллельному тесту) в колумбийской [19] и в латышской [20] популяциях; интронный вариант rs35280891 (p = 6,88 х 10-7; OR = 19,846, 95% CI: 5,728-68,761,) ассоциирован с тяжелым течением в сербской популяции [21].
Интересно отметить, что минорный аллель А rs17713054 является частью расширенного гаплотипа, унаследованного от неандертальцев. На сегодняшний
день этот гаплотип представляет собой основной генетический фактор риска, связанный с тяжелыми симптомами после заражения SARS-CoV-2. Различия в частоте гаплотипа риска у населения Южной и Восточной Азии привели к предположению, что некое селективное давление, возможно, связанное с холерой, привело к его распространению среди жителей Южной Азии [22].
В нашем исследовании мы подтвердили сильную ассоциацию rs17713054 с тяжелым течением (OR = 4,56; р = 0,0025) в рецессивной модели.
Было также показано, что во всех вариантах выборок как минорный аллель варианта rs1800629 в гене TNF — Т, так и минорный генотип ТТ являются фактором риска развития тяжелой формы при значениях OR = 1,5; р = 0,02.
Известно, что вариант rs1800629 гена TNF связан с необходимостью респираторной поддержки и ее большей продолжительностью у пациентов с COVID-19 [23]. Кроме патогенных вариантов нами был также обнаружен rs1898830 в гене TLR2 с протективным действием (OR = 0,72).
Проведя исследование ассоциации 10 SNP, мы обнаружили, что два полиморфизма обладают достаточной силой для оценки рисков тяжелого течения.
Проведено исследование ассоциации LZTFL1 rs17713054 в РФ на сбалансированной выборке пациентов с тяжелым и легким течением, показавшим убедительную ассоциацию с тяжестью течения. Это подтверждает гипотезу о функциональном значении данного полиморфизма.
Простой, быстрый и недорогой тест, направленный на прогнозирование рисков тяжелого течения COVID-19 на основе индивидуальных полиморфизмов ДНК, был бы полезен как для стратификации пациентов с высоким и низким риском осложнений, а также, что более важно, для оценки тяжести заболевания в популяции здоровых людей в случаях инфицирования и последующего заболевания как перспективный прогностический тест.
ВЫВОДЫ
В генетическом исследовании выборки из 713 пациентов с подтвержденным диагнозом COVID-19 были найдены ключевые SNP в генах TNF и LZTFL1, связанные с тяжелым течением заболевания. Охарактеризован новый генетический маркер в гене LZTFL1, прогностическая ценность которого составила более 91%. Отклонения частот ассоциированных с рисками тяжелого течения COVID-19 генетических полиморфизмов от равновесия Харди-Вайнберга имеют важное эпидемиологическое значение и требуют дальнейшего изучения. Варианты rs1800629 в гене TNF и rs17713054 в гене LZTFL1 могут быть рассмотрены как потенциальные кандидаты для дальнейшего анализа на расширенных выборках. На основании проведенного клинического исследования можно сделать ряд выводов для дальнейшей клинической практики. Общий поток пациентов, поступающих в клиники с диагнозом острой респираторной коронавирусной инфекции, гетерогенен и на ранних этапах может быть разделен на группы пациентов с потенциально легким течением и потенциально тяжелым течением заболевания по факту наличия или отсутствия в генотипе пациента выбранных полиморфизмов TNF rs1800629 и LZTFL1 rs17713054. Полиморфизмы в генах LZTFL1 и TNF могут быть использованы не только как прогностические биомаркеры, но и как предиктивные и представлять собой научное обоснование новых подходов в генотип-ориентированной терапии больных с тяжелым поражением
легких, связанным с СОУЮ-19. Настоящее исследование предполагает перспективную возможность организации
Литература
1. Hu B, et al., Characteristics of SARS-CoV-2 and COVID-19. Nat Rev Microbiol. 2021; 19 (3): 141-54.
2. Gupta K, et al., Systematic review and meta-analysis of human genetic variants contributing to COVID-19 susceptibility and severity. Gene. 2022; 844: 146790.
3. Severe Covid GG, et al., Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N Engl J Med. 2020; 383 (16): 1522-34.
4. Colona VL, Vasiliou V, Watt J, Novelli G, Reichardt JKV. Update on human genetic susceptibility to COVID-19: susceptibility to virus and response. Hum Genomics. 2021; 15 (1): 57.
5. Nakanishi T, Pigazzini S, Degenhardt F, Cordioli M, Butler-Laporte G, Maya-Miles D, et al. Age-dependent impact of the major common genetic risk factor for COVID-19 on severity and mortality. J Clin Invest. 2021; 131 (23): e152386.
6. Kwok AJ, Mentzer A, Knight JC. Host genetics and infectious disease: new tools, insights and translational opportunities. Nat Rev Genet. 2021; 22 (3): 137-53.
7. El-Hefnawy SM, et al, COVID-19 susceptibility, severity, clinical outcome and Toll-like receptor (7) mRNA expression driven by TLR7 gene polymorphism (rs3853839) in middle-aged individuals without previous comorbidities. Gene Rep. 2022; 27: 101612.
8. Casanova JL, Abel L. From rare disorders of immunity to common determinants of infection: Following the mechanistic thread. Cell. 2022; 185 (17): 3086-103.
9. Mohlendick B, et al. The GNB3 c.825C>T (rs5443) polymorphism and protection against fatal outcome of corona virus disease 2019 (COVID-19). Front Genet. 2022; 13: 960731.
10. Chen T, et al. A Low-Producing Haplotype of Interleukin-6 Disrupting CTCF Binding Is Protective against Severe COVID-19. mBio. 2021; 12 (5): e0137221.
11. Minashkin MM, et al. The Role of Genetic Factors in the Development of Acute Respiratory Viral Infection COVID-19: Predicting Severe Course and Outcomes. Biomedicines. 2022; 10 (3).
12. Fricke-Galindo I, Buendía-Roldán I, Ruiz A, Palacios Y, Pérez-Rubio G, de Jesus Hernández-Zenteno R, et al. TNFRSF1B and TNF Variants Are Associated With Differences in Levels of Soluble Tumor Necrosis Factor Receptors in Patients With Severe
скрининговых исследований с очевидной экономической доступностью и эффективностью.
COVID-19. J Infect Dis. 2022; 226 (5): 778-87.
13. Downes DJ, et al. Identification of LZTFL1 as a candidate effector gene at a COVID-19 risk locus. Nat Genet. 2021; 53 (11): 1606-15.
14. Fink-Baldauf IM, et al. CRISPRi links COVID-19 GWAS loci to LZTFL1 and RAVER1. EBioMedicine. 2022; 75: 103806.
15. Marion V, et al., Exome sequencing identifies mutations in LZTFL1, a BBSome and smoothened trafficking regulator, in a family with Bardet--Biedl syndrome with situs inversus and insertional Polydactyly. J Med Genet. 2012; 49 (5): 317-21.
16. Stewart CA, Gay CM, Ramkumar K, Cargill KR, Cardnell RJ, Nilsson MB, et al. Lung Cancer Models Reveal Severe Acute ResPiratory Syndrome Coronavirus 2-Induced EPithelial-to-Mesenchymal Transition Contributes to Coronavirus Disease 2019 Pathophysiology. J Thorac Oncol. 2021; 16 (11): 1821-39.
17. Redin C, CW. Thorball, and J. Fellay, Host genomics of SARS-CoV-2 infection. Eur J Hum Genet. 2022; 30 (8): 908-14.
18. Wei Q, Chen ZH, Wang L, Zhang T, Duan L, Behrens C, et al. LZTFL1 suppresses lung tumorigenesis by maintaining differentiation of lung epithelial cells. Oncogene. 2016; 35 (20): 2655-63.
19. Angulo-Aguado M, et al., Association Between the LZTFL1 rs11385942 Polymorphism and COVID-19 Severity in Colombian Population. Front Med (Lausanne). 2022; 9: 910098.
20. Rescenko R, et al. Replication of LZTFL1 Gene Region as a Susceptibility Locus for COVID-19 in Latvian Population. Virol Sin. 2021; 36 (5): 1241-4.
21. Zecevic M, et al. Genome-Wide Association Study of COVID-19 Outcomes Reveals Novel Host Genetic Risk Loci in the Serbian Population. Front Genet. 2022; 13: 911010.
22. Niemi MEK, Daly MJ, Ganna A. The human genetic epidemiology of COVID-19. Nat Rev Genet. 2022; 23 (9): 533-46.
23. Fishchuk L, et al. Modifying effects of TNF-alpha, IL-6 and VDR genes on the development risk and the course of COVID-19. Pilot study. Drug Metab Pers Ther. 2021; 37 (2): 133-9.
24. Kudryavtsev YS, Beregov MM, Berdalin AB, Lelyuk VG. Comparison of the Main Staging Systems for Assessing the Severity of Lung Injury in Patients with COVID-19 and Evaluation of Their Predictive Value. Journal of radiology and nuclear medicine. 2021; 102 (5): 296-303.
References
1. Hu B, et al., Characteristics of SARS-CoV-2 and COVID-19. Nat Rev Microbiol. 2021; 19 (3): 141-54.
2. Gupta K, et al., Systematic review and meta-analysis of human genetic variants contributing to COVID-19 susceptibility and severity. Gene. 2022; 844: 146790.
3. Severe Covid GG, et al., Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N Engl J Med. 2020; 383 (16): 1522-34.
4. Colona VL, Vasiliou V, Watt J, Novelli G, Reichardt JKV. Update on human genetic susceptibility to COVID-19: susceptibility to virus and response. Hum Genomics. 2021; 15 (1): 57.
5. Nakanishi T, Pigazzini S, Degenhardt F, Cordioli M, Butler-Laporte G, Maya-Miles D, et al. Age-dependent impact of the major common genetic risk factor for COVID-19 on severity and mortality. J Clin Invest. 2021; 131 (23): e152386.
6. Kwok AJ, Mentzer A, Knight JC. Host genetics and infectious disease: new tools, insights and translational opportunities. Nat Rev Genet. 2021; 22 (3): 137-53.
7. El-Hefnawy SM, et al, COVID-19 susceptibility, severity, clinical outcome and Toll-like receptor (7) mRNA expression driven by TLR7 gene polymorphism (rs3853839) in middle-aged individuals
without previous comorbidities. Gene Rep. 2022; 27: 101612.
8. Casanova JL, Abel L. From rare disorders of immunity to common determinants of infection: Following the mechanistic thread. Cell. 2022; 185 (17): 3086-103.
9. Mohlendick B, et al. The GNB3 c.825C>T (rs5443) polymorphism and protection against fatal outcome of corona virus disease 2019 (COVID-19). Front Genet. 2022; 13: 960731.
10. Chen T, et al. A Low-Producing Haplotype of Interleukin-6 Disrupting CTCF Binding Is Protective against Severe COVID-19. mBio. 2021; 12 (5): e0137221.
11. Minashkin MM, et al. The Role of Genetic Factors in the Development of Acute Respiratory Viral Infection COVID-19: Predicting Severe Course and Outcomes. Biomedicines. 2022; 10 (3).
12. Fricke-Galindo I, Buendía-Roldán I, Ruiz A, Palacios Y, Pérez-Rubio G, de Jesus Hernández-Zenteno R, et al. TNFRSF1B and TNF Variants Are Associated With Differences in Levels of Soluble Tumor Necrosis Factor Receptors in Patients With Severe COVID-19. J Infect Dis. 2022; 226 (5): 778-87.
13. Downes DJ, et al. Identification of LZTFL1 as a candidate effector gene at a COVID-19 risk locus. Nat Genet. 2021; 53 (11): 1606-15.
14. Fink-Baldauf IM, et al. CRISPRi links COVID-19 GWAS loci to LZTFL1 and RAVER1. EBioMedicine. 2022; 75: 103806.
15. Marion V, et al., Exome sequencing identifies mutations in LZTFL1, a BBSome and smoothened trafficking regulator, in a family with Bardet--Biedl syndrome with situs inversus and insertional Polydactyly. J Med Genet. 2012; 49 (5): 317-21.
16. Stewart CA, Gay CM, Ramkumar K, Cargill KR, Cardnell RJ, Nilsson MB, et al. Lung Cancer Models Reveal Severe Acute ResPiratory Syndrome Coronavirus 2-Induced EPithelial-to-Mesenchymal Transition Contributes to Coronavirus Disease 2019 Pathophysiology. J Thorac Oncol. 2021; 16 (11): 1821-39.
17. Redin C, CW. Thorball, and J. Fellay, Host genomics of SARS-CoV-2 infection. Eur J Hum Genet. 2022; 30 (8): 908-14.
18. Wei Q, Chen ZH, Wang L, Zhang T, Duan L, Behrens C, et al. LZTFL1 suppresses lung tumorigenesis by maintaining differentiation of lung epithelial cells. Oncogene. 2016; 35 (20): 2655-63.
19. Angulo-Aguado M, et al., Association Between the LZTFL1
rs11385942 Polymorphism and COVID-19 Severity in Colombian Population. Front Med (Lausanne). 2022; 9: 910098.
20. Rescenko R, et al. Replication of LZTFL1 Gene Region as a Susceptibility Locus for COVID-19 in Latvian Population. Virol Sin. 2021; 36 (5): 1241-4.
21. Zecevic M, et al. Genome-Wide Association Study of COVID-19 Outcomes Reveals Novel Host Genetic Risk Loci in the Serbian Population. Front Genet. 2022; 13: 911010.
22. Niemi MEK, Daly MJ, Ganna A. The human genetic epidemiology of COVID-19. Nat Rev Genet. 2022; 23 (9): 533-46.
23. Fishchuk L, et al. Modifying effects of TNF-alpha, IL-6 and VDR genes on the development risk and the course of COVID-19. Pilot study. Drug Metab Pers Ther. 2021; 37 (2): 133-9.
24. Kudryavtsev YS, Beregov MM, Berdalin AB, Lelyuk VG. Comparison of the Main Staging Systems for Assessing the Severity of Lung Injury in Patients with COVID-19 and Evaluation of Their Predictive Value. Journal of radiology and nuclear medicine. 2021; 102 (5): 296-303.
