Введение
Группа спирохет Borrelia burgdorferi sensu lato включает более 20 видов [1], из которых в России получены изоляты B. garinii, B. afzelii, B. burgdorferi sensu stricto, B. valaisiana, B. lusitaniae, B. spielmanii, B. finlandensis и B. bavariensis [2, 3]. Известно предложение о выделении спирохет группы B. burgdorferi sensu lato в отдельный род Borreliella [4, LPSN database, NCBI taxonomy database]. Поскольку оно остается остро дискуссионным [1], в данной работе мы используем прежнее родовое название спирохет этой группы.
B. bavariensis описана как самостоятельный вид рода Borrelia, благодаря мультилокусному сиквенс-типированию (MLST) генов домашнего хозяйства боррелий [5, 6]. Этим методом, а также мультилокусным сиквенс-анализом (MLSA) показано, что B. garinii подгруппы NT29, относятся к B. bavariensis [2, 6—9]. Этот вид широко распространен в Евразии и циркулирует среди различных переносчиков и резервуарных хозяев [10]. Различают европейскую и азиатскую популяции (по терминологии авторов [8, 10]) B. bavariensis. Далее мы называем эти популяции геномными подгруппами, исходя из понятий, которые были предложены ранее для описания внутривидового генетического полиморфизма боррелий [2, 11, 12].
На территории Европы и Восточной Азии изоляты B. bavariensis выделены из биоптатов кожи, спинномозговой и синовиальной жидкостей больных ИКБ [6, 8]. В России боррелии этого вида, выделенные от пациентов, ранее не описаны.
Природные очаги иксодовых клещевых боррелиозов (ИКБ) распространены в России в основном в лесной зоне. По уровню заболеваемости ИКБ занимают первое место среди всех природно-очаговых инфекций [13]. Показатели ее интенсивности в Пермском крае, где были получены охарактеризованные в данной работе изоляты, ежегодно значительно превышают среднероссийский уровень [14, 15]. Особенности клинических проявлений ИКБ описаны в ряде публикаций. В различных регионах страны до 80% заболеваний сопровождаются мигрирующей эритемой [13, 14, 16], а на территории Пермского края, по многолетним данным, — около 54% [15].
Цель нашего исследования — уточнение таксономического статуса одного из возбудителей эритемной формы ИКБ в Восточной Европе.
Материал и методы
Исследованы 23 изолята от людей, больных эритемной формой ИКБ, из музея боррелий [3] лаборатории переносчиков инфекций ФГБУ «Национального исследовательского центра эпидемиологии и микробиологии имени почетного академика Н.Ф. Гамалеи» Минздрава России (Москва), которые ранее по результатам изучения межгенного спейсера rrfA-rrlB отнесены к геномной подгруппе NT29 B. garinii [17]. Все изоляты получены от пациентов Краевой инфекционной больницы г. Пермь, который расположен на востоке Восточной Европы (58°00’ с.ш.; 56°15’ в.д.), путем посева на среду BSK биоптатов кожи из периферической части мигрирующей эритемы (17 изолятов) или плазмы крови (6 изолятов) больных ИКБ. Культивирование боррелий осуществляли при температуре 33 °C в течение 1 месяца [18, 19]. Клиническая картина эритемной формы ИКБ в Пермском крае, в частности, пациентов, от которых получены изоляты, описана ранее [14, 18]. В соответствии с принятой нами системой маркировки музейных изолятов перед их журнальными номерами проставляются первые буквы латинского названия источника изоляции. В частности, «Hs» означает «Homo sapiens», т.е. изолят, полученный от человека.
Выделение бактериальной ДНК проведено с помощью коммерческого набора «Проба-НК» («ДНК-технология», Москва). MLSA выполнен по схеме, предложенной D.Richter и соавт. [7]. Исследовали 6 генов домашнего хозяйства боррелий rrs, fla, groEL, hbb, recA, ospA, а также межгенный спейсер rrfA-rrlB. Размер локусов каждого изолята составил от 150 до 501 п.н. Полимеразная цепная реакция (ПЦР) проведена в термоциклере Eppendorf Mastercycler nexus (Eppendorf, Germany) в полном соответствии с описанной рекомендацией [7]. Секвенирование полученных ампликонов произведено с помощью набора реактивов ABI PRISM BigDye Terminator v. 3.1 с последующим анализом продуктов реакции на автоматическом секвенаторе Applied Biosystems 3730 DNA Analyser в Межинститутском Центре коллективного пользования «Геном» (Институт молекулярной биологии им. В.А. Энгельгардта РАН, Москва).
Информация об исследованных нами изолятах и сведения о штаммах из международных баз данных представлены в табл. 1. Нуклеотидные последовательности локусов анализируемых генов сравнивали с нуклеотидными последовательностями штаммов боррелий разных видов из баз данных GenBank (INSDC), а также PubMLST Borrelia spp, и анализировали с использованием пакета программ BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Выравнивание нуклеотидных последовательностей проведено с применением алгоритма MUSCLE. Дистанции рассчитывали с использованием поправки Джукса и Кантора [20]. Все объединенные нуклеотидные последовательности (2026 п.н.) исследованных изолятов (за исключением нуклеотидных последовательностей гена ospA штамма SZ и гена hbb изолята Hs-139, которые не удалось получить) использованы для построения дендрограммы методом невзвешенного попарного среднего (UPGMA) при величине bootstrap 1000 повторов с использованием пакета программ MEGA-X (версия 10.2.4).
Таблица 1. Информация об изолятах боррелий, исследованных в данной работе, а также из баз данных GenBank и PubMLST
Номер изолята | Источник изоляции | Страна изоляции штамма | Ссылка на результат генотипирования | Год изоляции |
Hs-1; Hs-2; Hs-3; Hs-4; Hs-6; Hs-7; Hs-9; Hs-10; Hs-11; Hs-12. Hs-52; Hs-55; Hs-58. Hs-61; Hs-64; Hs-65. Hs-80. | Кожный биоптат человека | Россия | [17] | 1992 1997 1998 1999 |
Hs-104; Hs-115; Hs-119; Hs-125; Hs-128; Hs-139 | Плазма крови человека | Россия | [17] | 2006 |
PBi | Ликвор человека | Германия | [5] | 1986 |
Lubl25 | Кожный биоптат человека | Словения | [21] | 1995 |
NT29 | Клещ I. persulcatus | Япония | [5] | 1994 |
SZ | Клещ Dermacentor sp. | Китай | [30] | 2010 |
Prm7564-11 | Клещ I. persulcatus | Россия | [21] | 2011 |
Prm7569-11 | Клещ I. persulcatus | Россия | [21] | 2011 |
Hiratsuka | Кожный биоптат человека | Япония | [21] | 2008 |
20047T | Клещ I. ricinus | Франция | [7] | 1997 |
IPT114 | Клещ I. ricinus | Франция | [7] | 2011 |
VS461T | Клещ I. ricinus | Швейцария | [7] | 1992 |
PGau | Кожный биоптат человека | Германия | [7] | 1988 |
A14S | Кожный биоптат человека | Нидерланды | [7] | 1992 |
PotiB2T | Клещ I. ricinus | Португалия | [7] | 1997 |
VS116T | Клещ I. ricinus | Швейцария | [7] | 1997 |
B31T | Клещ I. scapularis | США | [7] | 1984 |
IP1 | Ликвор человека | Франция | [7] | 1992 |
В базу данных Genbank депонированы 11 нуклеотидных последовательностей локусов исследованных генов размером более 200 п.н. (номера доступа MT713138-MT713148). В базу данных European Nucleotide Archive депонированы 5 нуклеотидных последовательностей размером менее 200 п. н. (номера доступа LR828307-LR828308, LR861813, LR862289-LR862290). Эти фрагменты участков генома характеризуют описанные геноварианты евразийской геномной подгруппы B. bavariensis и выбраны по дендрограмме (рис. 1) и табл. 2, исходя из показателей наименьшего сходства между ними.
Дендрограмма результатов MLSA локусов 7 генов, исследованных и использованных для сравнения изолятов.
Построена методом UPGMA в MEGA-X (см. Материалы и методы). Bootstrap поддержка 1000 повторов. На дендрограмме в круглых скобках приведено обозначение штамма, в квадратных — его сиквенс-типа [ST] из базы данных Borrelia spp. MLST (https://pubmlst.org/organisms/borrelia-spp).
Римские цифры в квадратных скобках: [I] — номера доступа нуклеотидных последовательностей в Genbank и ENA — MT713138-MT713142, LR862289, LR861813 — локусов всех исследованных нами генов изолята Hs-10; [II] — то же для изолята Hs-128: MT713144-MT713148, LR862290, LR828307; [III] — то же в ENA для межгенного спейсера rrfA-rrlB изолята Hs-58: LR828308; [IV] — то же в Genbank для гена groEL изолята Hs-55: MT713143.
Таблица 2. Сходство (в %) результатов MLSA изученных изолятов с результатами аналогичного анализа боррелий различных видов, имеющихся в базах данных GenBank и PubMLST
Hs-1 | Hs-6 | Hs-7 | Hs-55 | Hs-58 | Hs-80 | Hs-128 | Hs-139 | PBi | NT29 | SZ | Prm 7564-11 | Prm 7569-11 | Hiratsuka | 20047T | |
Hs-1 | 100 | ||||||||||||||
Hs-6 | 98,9 | 100 | |||||||||||||
Hs-7 | 98,9 | 99,9 | 100 | ||||||||||||
Hs-55 | 99,9 | 98,9 | 98,9 | 100 | |||||||||||
Hs-58 | 99,9 | 98,8 | 98,9 | 99,9 | 100 | ||||||||||
Hs-80 | 100,0 | 98,9 | 98,9 | 99,9 | 99,9 | 100 | |||||||||
Hs-128 | 98,9 | 100 | 99,9 | 98,8 | 98,8 | 98,9 | 100 | ||||||||
Hs-139 | 98,9 | 99,7 | 99,7 | 98,8 | 98,8 | 98,9 | 99,7 | 100 | |||||||
PBi | 98,4 | 98,1 | 98,2 | 98,3 | 98,3 | 98,4 | 98,1 | 97,8 | 100 | ||||||
NT29 | 100 | 98,8 | 98,8 | 99,9 | 99,9 | 100 | 98,8 | 98,8 | 98,4 | 100 | |||||
SZ | 99,9 | 99,2 | 99,3 | 99,9 | 99,9 | 99,9 | 99,2 | 99,3 | 98,7 | 99,9 | 100 | ||||
Prm 7564-11 | 98,9 | 100 | 99,9 | 98,8 | 98,8 | 98,8 | 100 | 99,7 | 98,1 | 98,8 | 99,2 | 100 | |||
Prm 7569-11 | 98,6 | 99,0 | 99,1 | 98,5 | 98,5 | 98,6 | 99,0 | 98,6 | 98,2 | 98,5 | 98,9 | 99,1 | 100 | ||
Hiratsuka | 98,9 | 99,9 | 100 | 98,9 | 98,9 | 98,9 | 99,5 | 99,7 | 98,2 | 98,8 | 99,3 | 99,5 | 99,1 | 100 | |
20047T | 97,6 | 97,7 | 97,8 | 97,5 | 97,5 | 97,6 | 97,7 | 97,3 | 98,1 | 97,6 | 98,5 | 97,7 | 97,8 | 97,8 | 100 |
Результаты и обсуждение
По результатам MLSA нуклеотидные последовательности исследованных нами изолятов боррелий от больных эритемной формой ИКБ отличались не менее чем на 2,2—2,7% от последовательностей штамма B. garinii 20047T (Sequence Type, ST 82 (Borrelia spp. MLST, https://pubmlst.org/organisms/borrelia-spp)). С типовым штаммом B. bavariensis PBi ST84 [6] они оказались схожи на 98,3—98,4% и отличались не менее чем на 5,4% от нуклеотидных последовательностей других видов боррелий, представленных в базе данных GenBank (см. табл. 2). В связи с ограничением объема статьи в этой таблице представлена информация о 8 изолятах из 23 исследованных нами. На дендрограмме, построенной по результатам MLSA, все изоляты Hs образуют общую ветвь с B. bavariensis и четко отделяются от сестринской ветви с B. garinii 20047T (рисунок). Это позволяет отнести изученные изоляты к B. bavariensis. Таким образом, впервые установлено, что этот вид боррелий вызывает на территории России эритемную форму ИКБ.
Все изоляты Hs оказались практически идентичными (сходство нуклеотидных последовательностей от 98,5 до 100%) штаммам NT29 [10], SZ ST131, Prm7564-11 ST433, Prm7569-11 ST374 и Hiratsuka ST128 [21] азиатской (по терминологии авторов [8, 10]) геномной подгруппы B. bavariensis (см. табл. 2) и группировались с ними на дендрограмме (см. рисунок). Геномная подгруппа B. garinii NT29, которая по современной таксономии рода Borrelia относится к B. bavariensis [6], обнаруженная у различных переносчиков и резервуарных хозяев [17, 22, 23], по особенностям распространения ранее была названа евразийской [2]. Поэтому мы предлагаем сохранить за этой геномной подгруппой B. bavariensis более точное название — «евразийская». Изоляция боррелий этой подгруппы, а также обнаружение маркеров ее ДНК у клещей Ixodes persulcatus, I. ricinus и I. pavlovskyi, мелких млекопитающих в Европейской [2, 22, 24, 25] и Азиатской [2, 22, 25—27] частях России говорят о широком распространении и циркуляции евразийской геномной подгруппы B. bavariensis, как это ранее было показано для геномной подгруппы B. garinii NT29 [28]. Боррелии этой геномной подгруппы, как показывают представленные нами данные, вызывают эритемную форму ИКБ у людей, заразившихся на территории восточной части Восточной Европы.
По результатам прежних молекулярно-биологических исследований изолятов, основанных на изучении полиморфизма длины фрагментов рестрикции межгенного спейсера rrfA-rrlB, принято считать, что наибольшее эпидемическое значение в России имеют B. garinii и B. afzelii, которые циркулируют среди широкого круга переносчиков и резервуарных хозяев [12, 13, 23, 24]. Данные, полученные нами, позволяют отнести B. bavariensis евразийской геномной подгруппы к числу этиологических агентов, вызывающих ИКБ на территории России. Это не исключает аналогичное значение европейской геномной подгруппы B. bavariensis, которая была обнаружена у клеща I. ricinus на территории Северо-Кавказского федерального округа России [2]. Результаты нашего исследования, полученные методом MLST (см. табл. 2, рисунок), подтверждают генетические отличия боррелий европейской (ST84, ST85) и евразийской (ST131, ST433, ST374, ST128 и др.) геномных подгрупп B. bavariensis, циркулирующих на территории России [22, 24].
Среди пермских изолятов евразийской геномной подгруппы выявлены два геноварианта, нуклеотидные последовательности которых отличались друг от друга более чем на 1% (см. табл. 2), а сходство между сиквенсами каждого геноварианта составляло от 99,7 до 100%. Большинство клинических изолятов (15 из биоптата кожи и 4 из плазмы крови пациентов) одного из геновариантов образовывали на дендрограмме (см. рисунок) общий кластер вместе с японским NT29 и китайским SZ ST131 штаммами евразийской геномной подгруппы B. bavariensis. Четыре изолята (по два из биоптата кожи и из плазмы крови пациентов) представляли другой геновариант, группируясь вместе с двумя пермскими штаммами B. bavariensis Prm7564-11, Prm7569-11 и японскимштаммом Hiratsuka (см. рисунок). Изоляты Hs-6 и Hs-128 группировались со штаммом Prm7564-11, относящимся к ST433, а изолят Hs-7 со штаммом Hiratsuka ST128. По результатам MLST оба сиквенс-типа отличаются всего 1 аллелем гена recG (база данных PubMLST) и обнаружены в России у клещей I. persulcatus. При этом ST433 найден пока только в Пермском крае, а ST128 в разных регионах России [22, 26]. На территории России ST131, с которым группировались большинство наших изолятов, не обнаружен у иксодовых клещей [22, 26], однако известно, что в Японии боррелии этого ST вызывают ИКБ у людей, а также обнаружены у грызунов [8].
Изоляты обоих геновариантов имеют определенные аллельные изменения в структуре генов. Например, выявлено, что они отличались по гену ospA пятью несинонимичными нуклеотидными заменами. Наши данные подтверждают сведения других исследователей [2, 5, 8, 18, 22, 26, 27] о большей генетической гетерогенности евразийской геномной подгруппы B. bavariensis. Возможное влияние различных геновариантов B. bavariensis на клинические проявления заболевания нуждается в дальнейшем изучении.
Совместная циркуляция на одной и той же территории различных видов боррелий, их геномных подгрупп и геновариантов выявлена ранее для B. garinii и B. afzelii [11, 22, 23, 25, 27]. Микст-изолят евразийской геномной подгруппы B. bavariensis (ранее B. garinii NT29) + B. garinii 20047T был изолирован от больного ИКБ в Пермском крае [19, 29]. Это дает основание предполагать возможность микст-инфицирования переносчиков и резервуарных хозяев двумя и более видами боррелий при участии B. bavariensis.
Эритемную форму ИКБ, протекающую как в серопозитивной, так и в серонегативной форме, могут вызывать различные виды боррелий [19], а их замедленный иммуногенез затрудняет раннюю серолого-этиологическую идентификацию инфекции [14]. Этому могут способствовать молекулярно-биологические методы лабораторной диагностики, особенно при серонегативном заболевании.
Заключение
Приведенные нами данные в совокупности с фактами, которые были опубликованы ранее другими исследователями, позволяют сделать вывод, что наряду с B. garinii и B. afzelii в Евразии широко распространена и имеет важное значение как возбудитель ИКБ B. bavariensis. На территории России циркулируют евразийская и европейская геномные подгруппы этого возбудителя. Евразийская геномная подгруппа B. bavariensis может вызывать эритемную форму ИКБ.
Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов.
The authors declare no conflicts of interest.
*Представлено на 14-м Международном симпозиуме по иксодовым клещам и клещевым заболеваниям. Германия, март 2021 г.
Presented at the 14-th International Symposium on Ticks and Tick-borne Diseases. Germany, March, 2021.