Введение
Бактерия Y. pestis является возбудителем особо опасного природно-очагового заболевания — чумы, способного к эпидемическому и пандемическому распространению. Глобальный характер проблемы чумы определяется существованием большого числа природных очагов чумы на различных континентах и высокой вирулентностью штаммов Y. pestis. В природе возбудитель чумы циркулирует на широком круге носителей и переносчиков, который включает 363 вида млекопитающих и 280 видов и подвидов блох [1, 2]. Считается, что географические и климатические факторы, а также природные катаклизмы, оказывают значительное влияние на эпизоотическую и эпидемическую активность природных очагов чумы [3, 4].
Штаммы Y. pestis вследствие их широкого распространения и адаптации к условиям разных природных биоценозов претерпевают различные изменения генома в ответ на воздействие окружающей среды. Эта изменчивость лежит в основе существующего генетического разнообразия штаммов Y. pestis в природных очагах чумы. Ранее проведенные молекулярно-генетические исследования показали, что в очагах чумы России и других стран СНГ распространены высоковирулентные и эпидемически значимые штаммы Y. pestis основного подвида, относящиеся к средневековому и античному биоварам [5]. В эндемичных очагах Кавказа, Горного Алтая и Центральной Азии встречаются штаммы Y. pestis неосновных подвидов, которые избирательно вирулентны для мелких мышевидных грызунов и зайцеобразных.
По генетической номенклатуре штаммы Y. pestis делят на ряд филогенетических линий и ветвей [6—8]. По современной классификации вид Y. pestis включает семь подвидов: кавказский (0.PE2), ангольский (0.PE3), центральноазиатский (0.PE4), улегейский (0.PE5), тибетский (0.PE7), цинхайский (0.PE10), основной. Последний включает 4 биовара — античный, представленный 5 филогенетическими линиями 0.ANT-4.ANT; средневековый — 2.MED; восточный — 1.ORI; промежуточный — 1.INT [5]. В XX веке широкое распространение в очагах Восточной Европы и Центральной Азии получили штаммы наиболее молодой филогенетической ветви 2.MED1 средневекового биовара, которые встречаются в 33 из 45 очагов чумы стран СНГ, в том числе в 7 из 11 природных очагов России. Штаммы 2.MED1 занимают 93,3% очаговых территорий СНГ с общей площадью 1959965 км2, в том числе 6 очагов сусликового (221347 км2), 25 очагов песчаночьего типа (1728676 км2), 2 очага сурочьего (9942 км2) типов. Эти штаммы распространены в очагах Прикаспия, Кавказа и Центральной Азии [9, 10]. В Центрально-Кавказском высокогорном очаге наряду с ветвью 2.MED1 средневекового биовара циркулирует и другая, наиболее ранняя из всех ветвей средневекового биовара 2.MED0. Также в первой половине XX века в очагах Северного и Северо-Западного Прикаспия циркулировали штаммы ветви 2.MED4, которые во второй половине этого века уже не выделялись. Штаммы средневекового биовара еще двух филогенетических ветвей 2.MED2 и 2.MED3 встречаются в Китае [8].
Распространение штаммов филогенетической ветви 2.MED1 в XX веке в Восточной Европе и Центральной Азии привело к многочисленным вспышкам чумы с высокими показателями летальности среди населения [9]. Ранее проведенные исследования показали, что штаммы 2.MED1 из очагов Восточной Европы и Центральной Азии в свою очередь делятся на две ветви — каспийскую (CS) (очаги Кавказа, Прикаспия) и центрально-азиатскую (CA) — очаги Центральной Азии (преимущественно Прибалхашский, Муюнкумский, Таукумский, Кызылкумский пустынные очаги) [9].
Штаммы 2.MED1 генетически однородны и их дифференциация является не простой, но актуальной для мониторинга очагов чумы России и других стран СНГ задачей. Для исследования штаммов Y. pestis средневекового биовара, имеющих важное эпидемиологическое значение, и для повышения эффективности эпидемиологического мониторинга этих территорий требуется наличие высокоразрешающих методов молекулярного типирования [5]. Как показал опыт молекулярно-генетических исследований, одним из наиболее высокоразрешающих методов типирования и определения филогеографической принадлежности штаммов Y. pestis является метод VNTR (от англ. variable number tandem repeats) и его мультилокусный вариант MLVA (от англ. multiple locus VNTR analysis) [12, 13]. Метод основан на анализе повторяющихся тандемных последовательностей, количество которых может изменяться в одной и той же области генома при сохранении структуры окружающих последовательностей. MLVA удовлетворяет большинству критериев, необходимых для эффективного типирования таких генетически мономорфных возбудителей как Y. pestis, Bacillus anthracis, Francisella tularensis, включая простоту использования, высокую дискриминационную способность и хорошую воспроизводимость. Метод MLVA25 успешно используется для молекулярно-эпидемиологических исследований Y. pestis; а также для клональной кластеризации определенных филогенетических групп; решения задач, связанных с микроэволюцией и таксономией возбудителя чумы [14—16]. Однако эффективность различных VNTR локусов и их значимость в генотипировании штаммов средневекового биовара остается не оцененной. Важной задачей является оптимизация метода MLVA25 для быстрой и высокоразрешающей дифференциации эпидемически значимых штаммов Y. pestis средневекового биовара, широко распространенных в природных очагах РФ и сопредельных стран на территории Восточной Европы и Центральной Азии.
Цель работы: оценка эффективности различных VNTR локусов и оптимизация метода MLVA25 для дифференциации штаммов Y. pestis средневекового биовара из природных очагов РФ и других стран СНГ, расположенных в Восточной Европе и Центральной Азии.
Материал и методы
В работе использованы полногеномные последовательности 220 штаммов Y. pestis средневекового биовара основного подвида из 31 очага чумы России и других стран СНГ, полученных в 1917—2003 гг. от носителей, переносчиков, больных людей, а также полногеномные последовательности 31 штамма других филогенетических линий Y. pestis из очагов чумы мира (табл. 1). Мультилокусное типирование MLVA25 проводили по 25 VNTR локусам с исключением из анализа трудно идентифицируемого локуса ms09 [13]. Праймеры для локусов ms46 и ms62 использовали согласно описанным [17]. Поиск участков VNTR осуществляли в полногеномных последовательностях штаммов Y. pestis с применением авторского скрипта, принцип действия которого основан на поиске последовательностей, фланкирующих VNTR локусы праймеров, с последующем вырезанием найденных участков и сохранением их в multi-FASTA файле. Подсчет количества тандемных повторов проводили в программе TRF [18], полученные значения округляли в меньшую сторону до целого числа. Степень вариабельности локусов оценивали посредством расчета индекса аллельного полиморфизма (h) с использованием формулы разнообразия Нея (Nei’s Diversity Index). Также был рассчитан соответствующий 95% доверительный интервал (CI от англ. Confidence Interval) [19].
Таблица 1. Штаммы Y. pestis, использованные в MLVA25 анализе
| Природный очаг (№) | Штаммы | Источник, год выделения | Филогенетическая ветвь |
| 01 — Центрально-Кавказский высокогорный | С627 | Горный сурок; 1986 г. | 2.MED0 |
| 02 — Терско-Сунженский низкогорный | 1116Д | Малый суслик; 1970 г. | 2.MED1 |
| 03 — Дагестанский равнинно-предгорный | С754, С791 | Малый суслик, блохи; 1999—2003 гг. | 2.MED1 |
| 06 — Зангезуро-Карабахский горный a/o | 146, 149 | Человек; 1931 г. | 2.MED4, 2.MED1 |
| 07 — Приараксинский низкогорный | 818 | Блохи; 1968 г. | 2.MED1 |
| 09 — Кобыстанский равнинно-предгорный | 261, 806, 807, 808, 812 ,813, 814, 815 | Краснохвостая песчанка, домовая мышь, серый хомячок, общественная полевка, блохи; 1953—1955 гг. | 2.MED1 |
| 10 — Мильско-Карабахский равнинно-предгорный | 1240 | Краснохвостая песчанка; 1967 г. | 2.MED1 |
| 11 — Джейранчёльский равнинно-предгорный | 7937, KM872, С36, С48 | Краснохвостая песчанка, блохи; 1956—1984 г. | 2.MED1 |
| 13 — Йоркский равнинно-предгорный | С293 | Краснохвостая песчанка; 1971 г. | 2.MED1 |
| 14 — Прикаспийский Северо-Западный степной | 9, 128, KM570, 1408, M1355, 372, 27 | Человек, блохи, серая крыса, малый суслик; 1923—1986 гг. | 2.MED1, 2.MED4 |
| 15 — Волго-Уральский степной | 107, 122, 165, 438, 3192, M577-M579 | Человек, малый суслик, обыкновенная полевка; 1929—1979 г. | 2.MED1 |
| 16 — Волго-Уральский песчаный | 4, 31, 768, 5, 3, 8, 15, 111, 174, 232, 260, 521, 546, 1077, KM639, KM642, 411, 484, 633, M591, M688, M696, M697, M1217, M1218, M1245, M1274, M1443-M1445, M1495, M1497, M1499-M1501, M1503-M1505, M1722, M1773, 106, 626, 685 | Человек, верблюд, полуденная песчанка, обыкновенная полевка, гребенщиковая песчанка, серая крыса, песчаный суслик, степная петрушка, домовая мышь, мохноногие тушканчики, блохи; 1917—2002 гг. | 2.MED4, 2.MED1 |
| 17 — Урало-Уильский степной | А1793, M1489, М1446, М960, M961, M963, M965, M1448, M1450, M1491, M1493 | Большая песчанка, блохи, малый суслик; 1971—1992 г. | 2.MED1 |
| 18 — Урало-Эмбенский пустынный | M1452, М1453, M1459, M1460, M1461, M1463, M1464, 48, 52, 7702, 6645, 4033 | Человек, блохи, большая песчанка; 1966—1990 гг. | 2.MED1 |
| 19 — Предустюртский пустынный | 224, M1467, 299 | Человек, верблюд, большая песчанка; 1967—1990 гг. | 2.MED1 |
| 20 — Устюртский пустынный | 85, 452, M549, М545 | Большая песчанка, верблюд; 1927—1975 гг. | 2.MED1 |
| 21 — Северо-Приаральский пустынный | 580, 578, 244, 247, 927-929, 621, 617, 605, 586, 246, 241, 600, 592, 609, 602, 242 | Человек, верблюд, блохи, полуденная песчанка, большая песчанка; 1945—1967 гг. | 2.MED1 |
| 23 — Мангышлакский пустынный | 381, 173, А1771, А1773 | Человек, сайга, большая песчанка; 1967—1978 гг. | 2.MED1 |
| 24 — Приаральско-Каракумский пустынный | 615, А1763, 4635, КМ938, 930, 550, 932, А1764, 1252, 6661, 1065, А1764 | Человек, нимфы клещей, блохи, блохи без определения, желтый суслик, малый суслик, полуденная песчанка, большая песчанка; 1945—1974 гг. | 2.MED1 |
| 25 — Каракумский пустынный | КМ816, М1283, М1285-М1287, М1295, М1299, М1301, М1302, М1304, 153, 650, 693, 868, M525, M557 | Человек, блохи, краснохвостая песчанка, большая песчанка; 1949—1985 гг. | 2.MED1 |
| 26 — Копетдагский пустынный | М521, M519 | Большая песчанка; 1967—1968 гг. | 2.MED1 |
| 27 — Кызылкумский пустынный | 19, 20, KM803 | Человек, блохи; 1924—1982 гг. | 2.MED1 |
| 28 — Мойынкумский пустынный | А1753, А1756, 1001, М1120, KM806, KM671, M974-M977, 1099, 135, М1197-М1199, М1201, М1507- М1524, М1614, М1615 | Блохи, блохи без определения, краснохвостая песчанка, большая песчанка; 1964—1993 гг. | 2.MED1 |
| 29 — Таукумский пустынный | 139, 141 | Большая песчанка; 1965—1966 гг. | 2.MED1 |
| 30 — Прибалхашский пустынный | 40, 505, KM805, KM936, KM939, A1920, A1760, A1778, 468, 1246, A197, A144 | Человек, блохи без определения, большая песчанка; 1938—1989 г. | 2.MED1 |
| 31 — Сарыджаский высокогорный a/o | 1564, 915, 333 | блохи, серый сурок; 1948—1961 гг. | 0.ANT5, 0.ANT3 |
| 32 — Верхненарынский высокогорный a/o | А-1694 | Блохи; 1928 | 0.ANT5 |
| 33 — Аксайский высокогорный а/о | А-1734 | Серый сурок; 1972 г. | 2.MED1 |
| 38 — Забайкальский-степной | 17 | Тарбаган № 9; 1923 г. | 2.MED1 |
| 40 — Таласский высокогорный | A1809, SCPM-O-B-7074, SCPM-O-B-7019 | Блохи; 1980 г. | 2.MED1, 0.PE4t |
| 43 — Прикаспийский песчаный | 258, М1220 | Человек, блохи; 1930—1999 гг. | 2.MED1 |
| Курдистан (Иран) | KIM10 | 1968 г. | 2.MED1 |
| Китай | 2501, B42003004, F1991016, Z176003, 91001 | Песчанка без определения, гималайский сурок, полевка брандта, желтогрудая крыса, серый сурок; 1976—2005 гг. | 2.MED1, 0.ANT2, 1.ORI2, 1.IN2, 0.PE4m (microtus) |
| Мадагаскар | IP275, MG05-1020 | Человек; 1995—2005 гг. | 1.ORI3 |
| Африка | Antiqua | Человек; 1965 г. | 1.ANT |
| Уганда | UG05-0454 | Человек; 2004 г. | 1.ANT |
| США, Angola | CO92, Dodson, Shasta, CA88-4125, A1122, Angola | Человек; 1939—1991 гг. | 1.ORI1, 0.PE3 |
| Монголия | I2422 | Блохи; 1974 г. | 0.PE5 |
| Киргизия | A1815 | Блохи; 1980 г. | 0.PE4 |
| Таджикистан | A1249 | Полевка без определения; 1970 г. | 0.PE4 |
| Армения | 3544 | Обыкновенная полевка; 1979 г. | 0.PE2 |
| NA* | Pestoides F, SCPM-O-B-6212, SCPM-O-B-6304, I2998, Pestoides A, EV76-CN | Монгольская пищуха; 1972—1982 гг. | 0.PE2, 0.PE5, 0.PE4, 1.ORI3 |
Примечание: NA* (Not Available) — не известно.
Полученные последовательности VNTR локусов вносили в программу Bionumerics 7.6.3. Построение филогенетического дерева осуществляли методами невзвешенного кластирования с арифметическим усреднением (UPGMA) и MST (Minimum spanning tree, минимальное остовное дерево). Оценку дискриминационной возможности системы типирования (Discrimination Power — DP) производили по индексу Симпсона [20].
Результаты и обсуждение
Оценка эффективности различных VNTR локусов в MLVA25 типировании штаммов Y. pestis средневекового биовара
В работе использованы полногеномные последовательности 220 штаммов Y. pestis средневекового биовара основного подвида из 31 очага чумы России и других стран СНГ, полученных в 1923—2003 гг., а также полногеномные последовательности 31 штамма других филогенетических линий (см. табл. 1). Поиск фрагментов VNTR осуществляли в полногеномных последовательностях штаммов Y. pestis с использованием авторского скрипта. Подсчет количества тандемных повторов по VNTR локусам проводили с помощью программы TRF (табл. 2). Современные методы секвенирования нового поколения генерируют относительно короткие последовательности [21], которые могут не охватывать весь VNTR локус. Это относится к VNTR локусам, которые могут иметь очень большое число копий [22]. По этой причине из анализа был исключен локус ms09. Для каждого VNTR локуса была определена степень вариабельности посредством расчета индекса аллельного полиморфизма (h) с использованием формулы разнообразия Нея (Nei’s Diversity Index) [23]. Для 24 использованных VNTR локусов количество аллелей, наблюдаемых по каждому локусу, варьировало от 1 до 18 (см. табл. 2).
Таблица 2. Анализ аллельного полиморфизма 24 локусов VNTR у штаммов Y. pestis, использованных в этой работе
| Локус | Штаммы средневекового биовара | Вся выборка штаммов | ||||||
| n* | h** | CI*** | Аллель*1 | n | h | CI | Аллель*1 | |
| ms01 | 2 | 0,239 | 0,173—0,305 | 7,6 | 5 | 0,337 | 0,271—0,403 | 7,6,8,9,10 |
| ms04 | 1 | — | — | 5 | 6 | 0,174 | 0,11—0,238 | 5,7,4,3,8,6 |
| ms05 | 3 | 0,106 | 0,050—0,162 | 9,10,8 | 4 | 0,217 | 0,15—0,284 | 9,10,8,3 |
| ms06 | 2 | 0,082 | 0,032—0,132 | 3,4 | 6 | 0,21 | 0,142—0,278 | 3,4,8,6,7,9 |
| ms07 | 4 | 0,203 | 0,132—0,274 | 7,9,6,8 | 5 | 0,22 | 0,152—0,288 | 7,9,8,6,2 |
| ms15 | 3 | 0,074 | 0,025—0,123 | 9,8,7 | 3 | 0,114 | 0,06—0,168 | 9,8,7 |
| ms20 | 1 | — | — | 6 | 2 | 0,178 | 0,118—0,238 | 6,7 |
| ms21 | 3 | 0,064 | 0,018—0,11 | 9,5,8 | 4 | 0,101 | 0,049—0,153 | 9,8,5,10 |
| ms35 | 2 | 0,056 | 0,013—0,099 | 10,9 | 3 | 0,12 | 0,066—0,174 | 10,9,11 |
| ms38 | 1 | — | — | 6 | 5 | 0,209 | 0,143—0,275 | 6,8,7,3,9 |
| ms40 | 1 | — | — | 8 | 5 | 0,101 | 0,049—0,153 | 8,7,9,2,11 |
| ms41 | 1 | — | — | 4 | 3 | 0,212 | 0,149—0,275 | 4,6,5 |
| ms44 | 1 | — | — | 7 | 1 | — | — | 7 |
| ms45 | 1 | — | — | 5 | 2 | 0,136 | 0,078—0,194 | 5,6 |
| ms46 | 16 | 0,896 | 0,878—0,914 | 8,14,9,13,15,7,12,5,11, 10,17,18,16,20,28,27 | 18 | 0,905 | 0,891—0,919 | 8,14,5,9,7,15,13,12,11,10,6,17,18,16,4,20,27,28 |
| ms51 | 1 | — | — | 4 | 1 | — | — | 4 |
| ms54 | 1 | — | — | 2 | 1 | — | — | 2 |
| ms56 | 4 | 0,523 | 0,503—0,543 | 8,9,7,6 | 5 | 0,583 | 0,554—0,612 | 8,9,7,6,5 |
| ms62 | 13 | 0,803 | 0,758—0,848 | 4,14,7,12,10,11,8,6,9,5, 13,15,17 | 14 | 0,837 | 0,802—0,873 | 4,7,14,12,10,8,11,9,6,5,13,15,17,3 |
| ms69 | 1 | — | — | 5 | 1 | — | — | 5 |
| ms70 | 7 | 0,702 | 0,66—0,744 | 6,7,9,4,8,5,10 | 9 | 0,746 | 0,709—0,783 | 6,7,9,4,8,5,3,10,12 |
| ms71 | 2 | 0,082 | 0,03—0,134 | 4,5,3 | 4 | 0,109 | 0,055—0,163 | 4,5,3,7 |
| ms73 | 2 | 0,013 | 0,002—0,024 | 4,5 | 3 | 0,194 | 0,132—0,256 | 4,5,2 |
| ms74 | 4 | 0,14 | 0,114—0,166 | 7,8,6,3 | 6 | 0,36 | 0,287—0,433 | 7,5,8,6,4,3 |
Примечание: n* — число аллелей; h** — индекс аллельного полиморфизма; CI*** — доверительный интервал; *1 — количество повторов по конкретному аллелю, представлено от большей частоты к меньшей.
Как показал проведенный MLVA анализ, наименее вариабельными по количеству аллелей (2-3) оказались локусы ms15, ms20, ms35, ms41, ms45, ms73, в то время как наиболее вариабельными по количеству аллелей (6-18) были ms04, ms06, ms46, ms62, ms70, ms74. Индекс аллельного полиморфизма, рассчитанный на основе частот аллелей, у исследованных штаммов варьировал в широких пределах от 0,101 до 0,905, за исключением тех локусов, где обнаружен только 1 аллель (см. табл. 2). По данным литературы по индексу аллельного полиморфизма вариабельность VNTR локусов делят на три типа: высоко вариабельные (h>0,25), умеренно вариабельные (0,11<h<0,25) и низко вариабельные (0,01<h<0,11) [24].
По итогам проведенного статистического анализа дискриминационной силы каждого локуса для штаммов средневекового биовара выявлены локусы с наибольшим генетическим полиморфизмом: ms01, ms07, ms46, ms56, ms62, ms70 и ms74 (рис. 1). Данные локусы включают от 2 до 16 аллелей с индексом аллельного полиморфизма 0,14-0,896. Наиболее вариабельные локусы с наибольшей эффективностью разделяют близкородственные штаммы по месту или времени происхождения, поэтому их включение в оптимизированную схему MLVA типирования средневекового биовара было целесообразным.
Рис. 1. График дискриминационной способности различных VNTR локусов для 220 штаммов Y.pestis средневекового биовара, выделенных в 1917—2003 гг. в 31 очаге чумы России и других стран СНГ.
DP — индекс дискриминационной способности.
Скорость мутаций значительно варьирует в зависимости от таксона, локуса и даже аллеля [12]. Участки последовательностей повторов во многих умеренно вариабельных и низко вариабельных VNTR в связи с низкой скоростью мутирования этих локусов долго остаются неизменными среди групп штаммов, выделенных в один период времени на одной территории. Поэтому даже незначительные изменения в количестве этих повторов могут повлиять на распределение штаммов на дендрограмме. Так, умеренно и низко вариабельные локусы ms05, ms06, ms15, ms21, ms71 по результатам проведенного нами статистического анализа оказали наибольшее влияние на разделение штаммов некоторых популяций средневекового биовара, таких как 2.MED4, а также 2.MED1 из Каракумского пустынного, Мангышлакского пустынного и Урало-Уильского степного очагов чумы. Из этого следует, что использование высоко вариабельных локусов — ms01, ms07, ms46, ms56, ms62, ms70 и ms74 и таких локусов с низкой скоростью накопления мутаций как ms05, ms06, ms15, ms21, ms71 позволяет добиться высокого разрешения в анализе и разделении штаммов средневекового биовара на генотипы. Эти VNTR локусы были использованы нами при разработке оптимизированной схемы MLVA12 типирования для штаммов Y. pestis средневекового биовара. При использовании в MLVA анализе других низко вариабельных локусов — ms35 и ms73 никакой зависимости топологии дерева штаммов средневекового биовара от их присутствия выявлено не было. Использование еще 10 других локусов ms04, ms20, ms38, ms40, ms41, ms44, ms45, ms51, ms54, ms69 из системы MLVA25 типирования в оптимизируемой схеме типирования средневекового биовара было также не целесообразным, поскольку каждый из них у штаммов средневекового биовара был представлен только одной аллелью.
Анализ эффективности использования меньшего набора VNTR локусов MLVA12 для дифференциации штаммов Y. pestis средневекового биовара
Для дифференциации штаммов Y. pestis средневекового биовара были отобраны локусы с наибольшим генетическим полиморфизмом (ms01, ms07, ms46, ms56, ms62, ms70 и ms74), а также некоторые умеренно и низко вариабельные локусы, оказывающие наибольшее влияние на распределение штаммов на дендрограмме (ms05, ms06, ms15, ms21, ms71).
Для оценки эффективности оптимизированного набора VNTR локусов — MLVA12 — и проверки дискриминирующей способности метода для штаммов средневекового биовара построена дендрограмма методом UPGMA с использованием программы Bionumerics 7.6.3 (рис. 2, см. https://mediasphera.ru/upload/medialibrary/files/Mol_genetika_2024_04_022_add.zip). Распределение штаммов на дендрограмме на основе 12 VNTR локусов и на дендрограмме, построенной на основе 24 локусов (данные не приводятся), было идентично (DP =0,98).
Штаммы 2.MED1 из очагов Восточной Европы и Центральной Азии делятся на две ветви — каспийскую (CS) и центральноазиатскую (CA). Список генотипов CA и CS представлен в табл. 3.
Таблица 3. Список MLVA12 генотипов 220 штаммов Y. pestis средневекового биовара из природных очагов Восточной Европы и Центральной Азии, использованных в работе
| Генотип | Штамм (Природный очаг, год выделения) |
| CA1 | 40 (Прибалхашский, 1961 г.) |
| CA2 | 1001 (Мойынкумский, 1965 г.) |
| CA3 | А-1763 (Приаральско-Каракумский, 1974 г.) |
| CA4 | 241, 242, 244, 246, 247 — (Северо-Приаральский, 1967 г.) |
| CA5 | КМ671 (Мойынкумский, 1966 г.); 1099 (Мойынкумский, 1964 г.); 650 (Каракумский, 1968 г.) |
| CA6 | 135 (Мойынкумский, 1966 г.) |
| CA7 | 261 (Кобыстанский, 1955 г.); 1240 (Мильско-Карабахский, 1967 г.) |
| CA8 | 505 (Прибалхашский, 1939 г.); 17 (Забайкальский, 1923 г.) |
| CA9 | 468 (Прибалхашский, 1938 г.) |
| CA10 | КМ816 (Каракумский, 1985 г.) |
| CA11 | М1283, М1285, М1286, М1287, М1295, М1301, М1302, М1304, М1299 — (Каракумский, 1984 г.) |
| CA12 | КМ806 (Мойынкумский, 1984 г.) |
| CA13 | М1519 (Мойынкумский, 1993 г.) |
| CA14 | М1220 (Прикаспийский песчаный, 1999 г.) |
| CA15 | М974, М977, М1198, М1507, М1508, М1510, М1515, М1520, М1521, М1522, М1197, М1509, М1511, М1516, М1518 — (Мойынкумский, 1993 г.); |
| CA16 | М975 (Мойынкумский, 1990 г.) |
| CA17 | М1523, М1524 — (Мойынкумский, 1993 г.) |
| CA18 | М1199, М1201, М1517 — (Мойынкумский, 1993 г.) |
| CA19 | М1513 (Мойынкумский, 1993 г.) |
| CA20 | М1514, М1614, М1512, М1615 — (Мойынкумский, 1993 г.) |
| CA21 | М976 (Мойынкумский, 1990 г.) |
| CA22 | КМ805 (Прибалхашский, 1963 г.); 141 (Таукумский, 1965 г.); М1489 (Урало-Уильский, 1992 г.); КМ803 (Кызылкумский, 1982 г.) |
| CA23 | А1760, А1778 — (Прибалхашский, 1972 г.) |
| CA24 | А1753 (Мойынкумский, 1973 г.) |
| CA25 | 1246 (Прибалхашский, 1970 г.) |
| CA26 | А1734 (Аксайский, 1972 г.) |
| CA27 | А-144, А-197 — (Прибалхашский, 1965 г.) |
| CA28 | А1920 (Прибалхашский, 1988 г.) |
| CA29 | А1809 (Таласский, 1980 г.) |
| CA30 | А-1764 (Приаральско-Каракумский, 1973 г.); КМ936, КМ939 — (Прибалхашский, 1989 г.) |
| CA31 | 2501 (Китай, 2005) |
| CA32 | 139 (Таукумский, 1966 г.) |
| CA33 | 19 (Каракумский, 1924 г.) |
| CA34 | 20 (Каракумский, 1924 г.) |
| CS1 | 106 (Волго-Уральский песчаный, 1928 г.) |
| CS2 | 149 (Зангезуро-Карабахский, 1931 г.) |
| CS3 | 521 (Волго-Уральский песчаный , 1944 г.) |
| CS4 | 633 (Волго-Уральский песчаный, 1945 г.) |
| CS5 | 52 (Урало-Эмбенский, 1966 г.) |
| CS6 | 626 (Волго-Уральский песчаный, 1945 г.); 9 (Прикаспийский Северо-Западный, 1923 г.) |
| CS7 | 372 (Прикаспийский Северо-Западный, 1937 г.), 165 (Волго-Уральский степный, 1932 г.) |
| CS8 | 411 (Волго-Уральский песчаный, 1938 г.) |
| CS9 | 812, 813 — (Кобыстанский, 1953 г.) |
| CS10 | 107 (Волго-Уральский степной, 1929 г.); 546 (Волго-Уральский песчаный, 1945 г.) |
| CS11 | 229 (Предустюртский, 1967 г.) |
| CS12 | 1408, КМ570 — (Прикаспийский Северо-Западный, 1972 г.) |
| CS13 | 818 (Приараксинский, 1968 г.) |
| CS14 | 815 (Кобыстанский, 1953 г.) |
| CS15 | С293 (Йоркский, 1971 г.) |
| CS16 | 1116Д (Терско-Сунженский, 1970 г.) |
| CS17 | 85 (Устюртский, 1927 г.) |
| CS18 | 685 (Волго-Уральский песчаный, 1947 г.) |
| CS19 | 5 (Волго-Уральский песчаный, 1917 г.) |
| CS20 | 128 (Прикаспийский Северо-Западный, 1929 г.) |
| CS21 | 3 (Волго-Уральский песчаный, 1917 г.) |
| CS22 | 8 (Волго-Уральский песчаный, 1922 г.) |
| CS23 | 15 (Волго-Уральский песчаный, 1923 г.) |
| CS24 | С36 (Джейранчёльский, 1956 г.) |
| CS25 | С48 (Джейранчёльский 1958 г.) |
| CS26 | М1453 (Урало-Эмбенский, 1990 г.) |
| CS27 | 768 (Волго-Уральский песчаный, 1968 г.); М696 (Волго-Уральски песчаный, 1981 г.) |
| CS28 | М963 (Урало-Уильский, 1971 г.); М1493 (Урало-Уильский, 1992 г.) |
| CS29 | 617, 621 — (Северо-Приаральский, 1945); 224 (Предустюртский, 1967 г.) |
| CS30 | М960, М1443, М1444, М1445, М1495, М1497, М1499, М1500, М1501, М1503, М1504 — (Урало-Уильский, 1989—1992 гг.) |
| CS31 | 48 (Урало-Эмбенский, 1966 г.) |
| CS32 | 111, 232, 1077 — (Волго-Уральский песчаный, 1962—1971 г.) |
| CS33 | КМ639 (Волго-Уральский песчаный, 1980 г.) |
| CS34 | 806, 807, 808 — (Кобыстанский, 1949—1953 гг.); 693 (Каракумский, 1949 г.) |
| CS35 | 814 (Кобыстанский, 1953 г.) |
| CS36 | М1452 (Урало-Эмбенский, 1990 г.) |
| CS37 | 7702 (Урало-Эмбенский, 1966 г.) |
| CS38 | 3192 (Волго-Уральский степной, 1959 г.) |
| CS39 | М688 (Волго-Уральский песчаный, 1980 г.) |
| CS40 | М697 (Волго-Уральский песчаный, 1981 г.) |
| CS41 | 258 (Прикаспийский песчаный, 1930 г.) |
| CS42 | 4635 (Урало-Уильский, 1959 г.) |
| CS43 | М965, М1446 — (Урало-Уильский, 1989 г.) |
| CS44 | М100 (Волго-Уральский степной, 1937 г.) |
| CS45 | М577, М578, М579, М591 — (Волго-Уральский степной, 1979 г.) |
| CS46 | КМ642 (Урало-Уильский, 1980 г.) |
| CS47 | М961 (Волго-Уральский песчаный, 1989 г.) |
| CS48 | М1217, М1467, М1218, М1245 — (Волго-Уральский песчаный, 1990—1999 гг.) |
| CS49 | М1773 (Волго-Уральский песчаный, 2002) |
| CS50 | М1448 (Урало-Уильский, 1990) |
| CS51 | А-1793 (Урало-Уильский, 1978 г.) |
| CS52 | М1450, М1491 — (Каракумский, 1989—1991 гг.) |
| CS53 | 868 (Каракумский, 1953 г.) |
| CS54 | 930 (Приаральско-Каракумский, 1955 г.) |
| CS55 | 1065 (Приаральско-Каракумский, 1966 г.) |
| CS56 | 1252 (Приаральско-Каракумский, 1966 г.); М519 (Копетдагский, 1967 г.) |
| CS57 | 6661 (Приаральско-Каракумский, 1966 г.) |
| CS58 | 550, КМ938 (Северо-Приаральский, 1955—1966 гг.); А1756 (Мойынкумский, 1973 г.) |
| CS59 | 928, 929 — (Северо-Приаральский, 1955 г.) |
| CS60 | 927 (Приаральско-Каракумский, 1955 г.); М556 (Волго-Уральский песчаный, 1945 г.); 153 (Каракумский, 1964 г.) |
| CS61 | 932 (Приаральско-Каракумский, 1965 г.) |
| CS62 | М1355 (Прикаспийский Северо-Западный, 1986 г.) |
| CS63 | С754, С791 — (Дагестанский, 1999—2003 гг.); КМ872 (Джейранчёльский, 1984 г.) |
| CS64 | М1274 (Волго-Уральский песчаный, 2001 г.) |
| CS65 | М1722 (Волго-Уральский песчаный, 1977 г.) |
| CS66 | 580, 615, 605, 600, 592, 602 — (Северо-Приаральский, 1945 г.) |
| CS67 | 578 (Северо-Приаральский, 1945 г.) |
| CS68 | 586 (Северо-Приаральский, 1945 г.) |
| CS69 | 609 (Северо-Приаральский, 1945 г.) |
| CS70 | 6645, 4033 — (Урало-Эмбенский, 1968 г.) |
| CS71 | М485, М487 (Мангышлакский, 1974 г.) |
| CS72 | 438 (Волго-Уральский степной, 1968 г.) |
| CS73 | 484 (Волго-Уральский песчаный, 1968 г.) |
| CS74 | М545, М549 — (Устюртский, 1975 г.) |
| CS75 | 7937 (Джейранчёльский, 1966 г.) |
| CS76 | 452 (Устюртский, 1962 г.) |
| CS77 | 173 (Мангышлакский, 1978 г.) |
| CS78 | М1459, М1460, М1463, М1464 — (Урало-Эмбенский, 1990 г.) |
| CS79 | М1461 (Урало-Эмбенский, 1990 г.) |
| CS80 | 122 (Волго-Уральский степной, 1929 г.) |
| CS81 | 260 (Волго-Уральский песчаный, 1924 г.) |
| CS82 | 174 (Волго-Уральский песчаный, 1932 г.) |
| CS83 | KIM10 (Иран, 1968 г.) |
| CS84 | М521 (Копетдагский, 1968 г.) |
| CS85 | M525 (Каракумский, 1970 г.) |
| CS86 | M557 (Каракумский, 1968 г.); 381 (Мангышлакский, 1967 г.) |
| 2.MED0 | M919 (Центрально-Кавказский, 1986 г.) |
| 2.MED4 | 146 (Зангезуро-Карабахский, 1931 г.), 27 (Прикаспийский Северо-Западный, 1924 г.), 4 (Волго-Уральский песчаный, 1917 г.), 31 (Волго-Уральский песчаный, 1924 г.) |
Анализ дендрограммы показал, что каждый очаг чумы имел свои уникальные типы MLVA. Большинство изолятов Y. pestis средневекового биовара (88%, 193/220) можно разделить на 4 основных кластера, №№1—4. Первый из них (№1) включает 58 штаммов. В этот кластер вошли штаммы из высокогорных, низкогорных и равнинно-предгорных природных очагов чумы Кавказа (CA7, CS2, CS9, CS13, CS14, CS15, CS16, CS24, CS25), Прибалхашского (CA1, CA8, CA9), Муюнкумского (CA2, CA5, CA6), Северо-Приаральского (CA4), Каракумского (CA11, CA12) пустынных очагов, а также Волго-Уральского степного (CS10) и Волго-Уральского песчаного (CS3, CS4, CS6, CS8, CS18, CS19, CS21, CS22, CS23) очагов чумы, выделенных в основном в первой половине XX века. Штаммы (27 штаммов) центральноазиатской ветви из смежных очагов чумы западного Казахстана выделены в 2 отдельные группы (CA1-CA9 и CA10-CA11). Данный кластер образован за счет отличий в аллелях по локусам ms01, ms07, ms46, ms56, ms74 и разделен на 37 MLVA типов. Кластер №2 (CA12-CA23) состоит из 36 штаммов центральноазиатской ветви, в основном из Мойынкумского пустынного очага чумы (CA12-CA13, CA15-CA21), выделенных в период 1984—1993 гг. Этот кластер выделен благодаря локусам: ms01, ms56, ms70 и ms74. Кластер №3 содержал наибольшее количество штаммов (79), разделяющихся на 42 MLVA генотипа, из которых 38 относятся к каспийской, а 4 к центральноазиатской ветвям эволюции средневекового биовара. Данный кластер в основном состоит из групп штаммов каспийской ветви, выделенных на территории Волго-Уральского междуречья, Астраханской и Волгоградской областей (CS28, CS30, CS41- CS44, CS46, CS50- CS51, CS62) Российской Федерации, а также Атырауской (CS27, CS31, CS32, CS33, CS36-CS40, CS45-CS49, CS64-CS65) и других западных областей (CA24- CA27, CS29, CS34, CS52-CS61) Республики Казахстан, из 14, 15, 16, 17, 18, 20 21, 24, 23, 43 очагов чумы (см. табл. 1). Также в третьем кластере отдельную группу составили пять штаммов центральноазиатской ветви (CA24-CA27), выделенные (1965—1973 гг.) в Прибалхашском, Муюнкумскм и Аксайском очагах чумы. Эти штаммы включены в №3 кластер за счет аллелей в локусах ms46, ms56 и ms62. Кластер №4 состоит из 20 штаммов (CS66-CS77) каспийской ветви, выделенных в период 1945-1978 гг. из 11, 15, 16, 18, 20, 21, 23, 24 очагов чумы (см. табл. 1). Этот кластер разделен на 12 MLVA типов за счет отличий в аллелях по локусам: ms46, ms62, ms70.
К линии средневекового биовара относятся также штаммы (12%), не вошедшие в эти 4 кластера. К ним принадлежат 2 штамма (генотип CA33, CA34) из Кызылкумского природного очага чумы, выделенные в 1924 г. во время вспышки чумы в с. Ак-Камыш. Их MLVA генотип отличается от других изолятов средневекового биовара по локусам ms07, ms21, ms70 и ms74. Также отдельно расположены штамм филогенетической линии 2.MED0 из Центрально-Кавказского высокогорного очага чумы, 4 штамма филогенетической линии 2.MED4 и 20 штаммов 2.MED1 (CA28-CA32, CS78-CS86), выделенные в период 1924—1970 гг. в Волго-Уральском степном и Волго-Уральском песчаном (CS80-CS82), Урало-Эмбенском (CS78-CS79), Каракумском (CS85-CS86) и Копетдагском (CS84) природных очагах чумы. Данные штаммы имели отличные повторы от остальных штаммов средневекового биовара в локусах ms05, ms46, ms62, ms70, ms74. Штаммы филогенетической линии 2.MED4 также отличались по локусам ms05, ms15, ms56.
Анализ кластеризации штаммов на дендрограмме свидетельствует о приуроченности MLVA генотипов к определенным территориям природных очагов чумы. Расположение на дендрограмме штаммов, входящих в кластерные комплексы и формирующих индивидуальные MLVA генотипы, доказывает их длительное сохранение в отдельных экологических нишах.
В результате статистического анализа выбранных нами 12 VNTR локусов у штаммов средневекового биовара обнаружены 126 индивидуальных генотипов, что соответствует DP=0,98. Это свидетельствует о высоком разрешении метода MLVA12 и о том, что комбинация этих 12 локусов может быть использована в качестве метода типирования Y. pestis с общепринятой вероятностью ошибки 5%. Данный показатель коррелирует с кластерным анализом, при котором исследуемые штаммы разделились на 126 MLVA генотипов. Использование предложенного набора VNTR локусов для штаммов средневекового биовара, циркулирующих на территории Восточной Европы и Центральной Азии, позволяет получить идентичный результат типирования, как при использовании полного набора локусов в MLVA25 анализе.
Локусы, у которых индекс полиморфизма равен нулю, не влияют на генотипирование штаммов, т.к. индивидуальные генотипы формируются за счет уникальных комбинаций чисел количества повторов во всех локусах. Из этого следует, что индекс дискриминационной возможности Симпсона (0,98) остается неизменным при исключении из выборки всех не вариабельных локусов (h=0). По данным индекса частот аллелей, посчитанного нами для штаммов средневекового биовара, были выявлены 10 локусов с нулевым диапазоном частот тандемных повторов — их аллели представлены 1 вариантом повтора (см. табл. 2). Такими локусами оказались: ms04, ms20, ms38, ms40, ms41, ms44, ms45, ms51, ms54, ms69. Независимо от их наличия или отсутствия в анализе количество MLVA генотипов средневекового биовара сохранялось неизменным и равнялось 126. Таким образом, использование этих локусов для внутрибиоварной дифференциации штаммов средневекового биовара, циркулирующих на территории Восточной Европы и Центральной Азии, методом MLVA нецелесообразно. У этих локусов единый вариант VNTR повтора присутствует в группе средневекового биовара, который также изредка встречается и во взятых нами в исследование штаммов других филогенетических линий. Четыре из этих 10 локусов — ms44, ms51, ms54 и ms69 — оказались мономорфными для всех групп штаммов, использованных в работе. Однако другими исследователями в локусах ms44, ms51, ms54 и ms69 были выявлены и другие варианты аллелей, отличавшиеся от наших, что связано с использованием в анализе большего числа штаммов Y.pestis других филогенетических линий [15]. Процент встречаемости единого аллеля средневекового биовара у оставшихся 6 из 10 локусов (ms04, ms20, ms38, ms40, ms41, ms45) у штаммов других биоваров составляет 22,5, 19,3, 3,1, 67,7, 0 и 9,3% соответственно, из чего следует, что применение комплекса этих 6 локусов обеспечивает дифференциацию средневекового биовара от других филогенетических групп Y. pestis. Анализ данных литературы подтверждает сделанное нами заключение о том, что у штаммов античного, восточного биоваров и других филогенетических линий, варианты аллелей этих 6 VNTR локусов отличаются от единого варианта аллеля штаммов средневекового биовара [15, 25—28]. У штаммов античного, восточного и промежуточного биоваров, а также кавказского подвида (0.PE2) и биовара микротус центральноазиатского подвида (0.PE4m), количество повторов варьировало в локусах: ms04 — 2-10; ms20 — 2-11; ms38 — 5-11; ms40 — 7-11; ms41 — 5-7; ms45 — 2-10. При комплексном использовании данных локусов, вероятность перекрытия по аллельному признаку минимальна, так как встречаемость доминирующего аллеля средневекового биовара в локусах, приведенных выше, единична. Из этого следует, что комбинация локусов ms04, ms20, ms38, ms40, ms41, ms45 дает возможность эффективно отделять штаммы Y. pestis средневекового биовара от античного и восточного биоваров, а также от штаммов других филогенетических линий (рис. 3).
Рис. 3. Дифференциация штаммов Y. pestis средневекового биовара на основе 6 локусов VNTR.
а — дендрограмма 251 штамма Y. pestis разных филогенетических линий с использованием 24 локусов VNTR; б — дендрограмма 251 штамма Y. pestis с использованием 6 маркерных локусов VNTR: ms04, ms20, ms38, ms40, ms41, ms45. Minimum Spanning tree (MST), программа Bionumerics 7.6.3.
Для визуализации результата использования 6 локусов VNTR (ms04, ms20, ms38, ms40, ms41, ms45) в MLVA анализе для быстрой дифференциации штаммов Y. pestis средневекового биовара, построена дендрограмма методом MST с использованием программы Bionumerics 7.6.3 (см. рис. 3).
На дендрограмме, построенной методом MST, MLVA типы представлены в виде кругов. Размер каждого круга определяется количеством изолятов с данным VNTR генотипом. На дендрограмме (см. рис. 3 а) штаммы сформировали отдельные группы в зависимости от принадлежности к филогенетической ветви или биовару. Большая часть исследуемых изолятов состоит из штаммов Y. pestis средневекового биовара основного подвида. При использовании в MLVA анализе 6 локусов VNTR (ms04, ms20, ms38, ms40, ms41, ms45) (см. рис. 3 б) штаммы средневекового биовара объединяются в отдельный кластер. При этом другие штаммы формируют другие кластеры, которые не воспроизводят порядок эволюции, но кластеризуют штаммы по филогеографической принадлежности и могут использоваться для дифференцирования штаммов средневекового биовара основного подвида от других штаммов Y. pestis.
Заключение
Таким образом, по результатам проведенного MLVA анализа при исследовании 220 штаммов Y. pestis средневекового биовара, полученных на территории Восточной Европы и Центральной Азии в период 1917—2003 гг., выявлены наиболее эффективные локусы с высокой дискриминационной силой, оказывающие наибольшее влияние на дифференциацию штаммов по генотипам. Для каждого локуса подсчитан индекс аллельного полиморфизма, в качестве интервальной оценки статистического параметра рассчитан доверительный интервал.
По результатам кластерного анализа установлено, что для MLVA типирования штаммов Y. pestis средневекового биовара основного подвида нецелесообразно использовать 12 локусов (ms04, ms20, ms35, ms38, ms40, ms41, ms44, ms45, ms51, ms54, ms69 и ms73) из-за их низкой дискриминационной силы.
По данным MLVA анализа выявлены сходные MLVA генотипы у штаммов, выделенных на смежных территориях. Полученные генотипы коррелируют с местом и временем выделения исследуемых штаммов. Эта особенность позволяет использовать данный метод для кластеризации штаммов, выделенных как внутри одного очага, так и для разделения штаммов, полученных на разных очаговых территориях. Метод MLVA дает возможность проследить микроэволюционную изменчивость Y. pestis в отдельных географических регионах, что может быть полезно для углубленного понимания популяционной структуры Y. pestis в разных очагах чумы, а также для повышения эффективности молекулярно-эпидемиологического мониторинга очаговых по чуме территорий.
Авторы заявляют об отсутствии конфликта интересов.