Molecular detection and phylogeography of the hypervirulent subtype of Mycobacterium tuberculosis Beijing 14717-15
- Authors: Mokrousov I.1, Vyazovaya A.1, Badleeva M.1, Gerasimova A.1, Belopolskaya O.1, Masharsky A.1, Kostyukova I.1, Zhdanova S.1, Mudarisova R.1, Avadenii I.1, Solovieva N.1, Skiba Y.1, Zhuravlev V.1, Pasechnik O.1, Ogarkov O.1
-
Affiliations:
- Санкт-Петербургский научно-исследовательский институт эпидемиологии и микробиологии им. Пастера
- Section: Статьи
- URL: https://iimmun.ru/iimm/article/view/2082
- DOI: https://doi.org/10.15789/2220-7619-MDA-2082
Cite item
Full Text
Abstract
Aim. To develop a method for the detection of Mycobacterium tuberculosis strains of the highly lethal, hypervirulent, and multidrug-resistant Beijing 14717-15-cluster and its application for screening retrospective collections of M. tuberculosis from regions of Russia. The collection included M. tuberculosis DNA samples collected within population studies. Spoligotyping and genotyping of 24 variable number of tandem repeats loci were performed according to standard protocols. Phylogenetic analysis was based on the whole genome sequencing data. Polymorphism at the genomic position 2423040 A>G, specific for the Beijing 14717-15 cluster, was detected using PCR-RFLP and HhaI restriction enzyme. The bioinformatic and phylogenetic analysis of whole genome sequencing data of 205 strains of the early ancient sublineage of the M. tuberculosis Beijing genotype showed that strains with the VNTR 14717-15 profile (according to MIRU-VNTRplus.org) and those close to it were grouped into one monophyletic cluster of related strains, which we defined as Beijing 14717-15-cluster. Among the SNPs specific for the 14717-15 cluster, we chose a functionally neutral SNP at the genomic position 2423040 A>G (Rv2161c Val33Ala) and developed a method for its detection in the PCR-RFLP format. Next, we applied this method to screening DNA collections from the regions of the European and Asian parts of Russia and Asian countries. The results demonstrate the presence of strains of the Beijing 14717-15 cluster only in the Asian part of Russia, with a peak in Buryatia (18%). We have developed a method for the detection of the hypervirulent and highly lethal genetic cluster M. tuberculosis Beijing 14717-15 based on the detection of a specific mutation in the Rv2161c gene using PCR followed by HhaI restriction of the PCR product and agarose gel electrophoresis to discriminate between wild-type and mutant alleles. The advantages of the proposed method are the speed, simple and unambiguous interpretation of the results, the ability to detect contamination and mixed infection, suitability for efficient analysis of large collections of M. tuberculosis strains.
Full Text
Введение
Патобиология штамма (устойчивость к антибиотикам, вирулентность, выживаемость, имуногенность) является ключевым фактором формирования структуры популяции Mycobacterium tuberculosis, которая, в свою очередь, определяет эпидемическую ситуацию по туберкулезу. Клиническая значимость структуры популяции M. tuberculosis состоит в том, что отдельные генетические варианты имеют большее клиническое и/или эпидемиологическое значение чем другие. Вид M. tuberculosis включает несколько крупных филогенетических линий, и генотип Beijing (основной компонент Восточно-Азиатской линии) является одним из наиболее изученных [1-3]. Генотип Beijing разделяют на две ветви или сублинии - более гетерогенную предковую/древнюю и эволюционно более молодую современную [4-6]. Предыдущие исследования показали связь сублиний Beijing с различными патологическими и клиническими проявлениями, которые влияют на трансмиссивность конкретных штаммов в человеческих популяциях. В рамках генотипа Beijing его древние и современные штаммы различаются по вирулентности, что также частично зависит от происхождения штамма.
Современная сублиния распространена по всему миру и включает несколько хорошо известных эпидемических или эндемических генетических кластеров штаммов. Напротив, штаммы древней сублинии генотипа Beijing редко встречаются за пределами Восточной Азии. Недавно два генетических кластера этой сублинии были обнаружены в азиатской части России (типы 1071-32 и 14717-15, согласно номенклатуре ресурса MIRU-VNTRplus.org), и оба были ассоциированы с множественной лекарственной устойчивостью (МЛУ) [7]. Положение данных кластеров в рамках глобальной филогении генотипа Beijing показано на Рис. 1. Один из кластеров, обозначенный как 1071-32, относительно распространен в Омске, Сибирь (7%), и встречается спорадически в разных регионах России с середины 1990-х годов, также его единичные изоляты выявлены в Восточной и Южной Европе [8]. Второй кластер 14717-15 был обнаружен только в азиатской части России, где он показал нарастающий градиент от 2,6% в Омске, Западная Сибирь, до 16% в Бурятии, Дальний Восток [9].
В мышиной модели штаммы этих кластеров продемонстрировали противоположные свойства [10]. Штамм 396 (14717-15-кластер из Бурятии, Дальний Восток) приводил к очень высокой смертности, большей потере веса, более высокой бактериальной нагрузке и более тяжелой патологии легких. Cравнение с опубликованными данными по другим российским эпидемическим штаммах генотипа Beijing (B0/W148, CAO, Central Asian/Russian [11]), показало, что штамм 396 продемонстрировал наивысшую летальность таким образом являясь наиболее летальным российским штаммом M. tuberculosis. При этом ретроспективное исследование больных в Омской области показало, что группа пациентов, инфицированных кластером Beijing 14717-15, имела и самую высокую долю летального исхода по сравнению с другими группами [12].
Рис. 1. Схема эволюции генотипа Beijing Mycobacterium tuberculosis, включая основные российские субтипы.
Fig. 1. Scheme of the evolution of the Beijing Mycobacterium tuberculosis genotype, including the main Russian subtypes.
Сейчас становится очевидным, что при внедрении персонализированного лечения туберкулеза следует учитывать не только лекарственную устойчивость, но и вирулентность штамма M. tuberculosis, что требует разработки методов быстрой детекции значимых гипервирулентных геновариантов. Выявление штаммов кластера Beijing 14717-15 возможно с использованием полногеномного секвенирования или генотипирования по 24 локусам VNTR с последующим сравнением с базой данных www.MIRU-VNTRplus.org [7, 9].
Цель исследования - разработка простого метода для выявления штаммов M. tuberculosis, относящихся к высоколетальному, гипервирулентному и мультирезистентному кластеру Beijing 14717-15 и применение метода для скрининга ретроспективных коллекций ДНК M. tuberculosis из различных регионов России.
Материалы и методы
Изученная коллекция включала образцы ДНК M. tuberculosis, собранные в период с 2015 по 2022 годы в результате популяционных исследований. Выделение ДНК из культуры М. tuberculosis, выращенной на среде Левенштейна-Йенсена, проводили по van Embden et al. [13] с использованием лизоцима, протеиназы К, додецилсульфата натрия и цетилтриметиламмонийбромида с последующей экстракцией с использованием смеси фенол-хлороформ-изоамиловый спирт (25:24:1), центрифугированием, осаждением изопропанолом, промывкой 70% этанолом, и растворением осадка в 30-50 мкл буфера ТЕ х1.
Сполиготипирование проводили по стандартному протоколу с использованием мембраны с иммобилизованными 43 олигонуклеотидами, изготовленной в нашей лаборатории [14]. Генотип Beijing был идентифицирован экспериментально или in silico на основе определения делеции RD207 (позиции 3120521-3127920 в геноме H37Rv, NC_000962.3). Генотипирование 24 локусов MIRU-VNTR проводили по Supply et al. [15], полученные числовые профили сравнивали с ресурсом MIRU-VNTRplus.org.
Филогенетический анализ штаммов древней сублинии генотипа Beijing на основе данных полногеномного секвенирования, полученных в нашей лаборатории (депонированы в SRA NCBI под номерами проектов PRJNA822891, PRJNA489691) или извлеченных из архива коротких прочтений SRA NCBI, был описан ранее [9]. Пакет Geneious 9.0 (Biomatters Ltd) использовали для картирования ридов на геном референтного штамма M. tuberculosis H37Rv (NC_00962.3). Онлайн-ресурс SAM-TB (https://samtb.uni-medica.com/index) использовали для анализа геномных данных, определения SNP и филогенетического анализа. Пакет MEGA7 [16] использовали для визуализации дендрограмм.
Полиморфизм в геномной позиции 2423040 A>G (позиция в гене Rv2161c 98T>C, аминокислотная замена Val(s)33Ala GTG/GCG), специфический для кластера Beijing 14717-15 (Рис. 2) определяли с использованием метода ПЦР-ПДРФ и рестриктазы HhaI. В случае данной мутации возникает дополнительный сайт рестрикции HhaI (сайт рестрикции GCGC). Длина продукта ПЦР - 150 пн. Размеры рестриктов при аллеле дикого типа – 60 и 89 пн, при мутации - 21, 60 и 68 пн (Рис. 3). Для ПЦР использовали два праймера GTCCGGCAGCTCTCCACCG и TGCAGTTCGTCACCGACCTGACC. ПЦР проводили при следующем температурном режиме: 95°С, 5 минут, 35 циклов 95°С, 40 сек, 67°С, 20 сек, 72°С, 20 сек; 72°С, 5 мин. Продукт ПЦР обрабатывали рестриктазой HhaI (37 °С, 2 часа), продукты рестрикции разделяли в 1,4% агарозном геле (Рис. 3).
Результаты и обсуждение
Проведенный нами биоинформационный и филогенетический анализ данных полногеномного секвенирования 205 штаммов ранней древней сублинии генотипа Beijing M. tuberculosis показал, что штаммы из Омска и Бурятии с профилем VNTR 14717-15 и близкие ему также группировались в один кластер родственных штаммов и на полногеномном дереве, который мы определили как Beijing 14717-15-кластер (Рис. 2).
Следует отметить что ссй кластер был монофилетическим и сравнение его штаммов с соседними ветвями на дендрограмме позволило выявить 58 специфических полиморфизмов SNP. Среди SNP специфичных для кластера 14717-15 мы выбрали функционально нейтральный SNP (PAM1=9867) в геномной позиции 2423040 A>G (Rv2161c Val33Ala) и разработали метод его детекции в формате ПЦР-ПДРФ. Фактически эволюционная нейтральность данной замены снижает риск ее селекции в процессе эволюции и адаптации и, таким образом, делает ее наиболее пригодной в качестве диагностического молекулярного маркера. Данная мутация приводит к появлению дополнительного сайта рестриктазы HhaI и может быть выявлена с помощью анализа HhaI-ПДРФ амплифицированной ПЦР-области (Рис. 3). Метод был оптимизирован с использованием ДНК изолятов с имевшейся информацией по их полным геномам у которых в данной позиции был или аллель дикого типа или мутантный аллель. Длина продукта ПЦР составляет 150 п.н. Поскольку в данном участке гена присутствует дополнительный инвариантый сайт HhaI, размеры фрагментов рестрикции составляют 60 и 89 п.н. в случае аллеля дикого типа и 60, 21, 68 п.н. в случае мутации (Рис. 3), верхний фрагмент хорошо различим при агарозном гель-электрофорезе даже при небольшом времени электрофореза и малой длине геля, что удобно практически.
Рис. 2. Дендрограмма штаммов древней сублинии генотипа Beijing из России и Кореи, построенная на основе полногеномых SNP на основе алгоритма максимального правдоподобия с использованием ресурса SAM-TB. Ветвь, включающая штаммы Beijing 14717-15-кластера, отмечена звездочкой. Шкала показывает количество SNP. Цифры курсивом рядом с узлами показывают процент бутстрапа. KR – Корея (в номере штамма).
Fig. 2. Dendrogram of strains of the ancient subline of the Beijing genotype from Russia and Korea, built on the basis of whole genome SNPs based on the maximum likelihood algorithm using the SAM-TB resource. The branch containing strains of the Beijing 14717-15 cluster is marked with an asterisk. The scale shows the number of SNPs. The numbers in italics next to the nodes show the bootstrap percentage. KR - Korea (in the strain number).
Рис. 3. HhaI ПЦР-ПДРФ детекция замены в геномной позиции 2423040 A>G (Rv2161c Val(s)33Ala). Дорожки 1-5 – Beijing 14717-15-кластер. Дорожки 6-7 – другие генотипы. М – маркер молекулярных весов 100 bp ladder (GE Healthcare).
Fig. 3. HhaI PCR-RFLP detection of substitution at genomic position 2423040 A>G (Rv2161c Val(s)33Ala). Lanes 1-5 - Beijing 14717-15-cluster. Lanes 6-7 - other genotypes. M, 100 bp ladder molecular weight marker (GE Healthcare).
Анализ ПЦР-ПДРФ был дополнительно применен к изолятам генотипа Beijing, которые представляли разные сублинии Beijing и имели разные профили VNTR. Все изоляты были ранее генотипированы методом мультилокусного VNTR-анализа, позволяющим точно определить принадлежность штамма к тому или иному VNTR-кластеру в том числе к кластеру 14717-15. Эти валидационные коллекции включали изоляты из России, Центральной и Восточной Азии. Анализ ПЦР-ПДРФ корректно отнес все изоляты с известным профилем VNTR 14717-15 и родственными профилями к кластеру Beijing 14717-15. Метод имеет 100% чувствительность и 100% специфичность для обнаружения данного кластера.
Далее мы применили метод ПЦР-ПДРФ для скрининга всех имеющихся коллекций ДНК из регионов Европейской и Азиатской частей России и стран Азии. Эти коллекции были из Восточной Сибири и Дальнего Востока (регионы с различной долей штаммов Beijing и Beijing 14717-15 [17]), из Омской области Западной Сибири (регион с небольшим, но постоянным присутствием штаммов Beijing 14717-15 в локальной популяции [7]), из Китая, Вьетнама и Японии (страны с достаточно высокой долей штаммов древних сублиний генотипа Beijing [18-20]). Результаты по всем проверенным коллекциям приведены в Таблице 1 и демонстрируют четкий пик кластера 14717-15 в Бурятии.
Дополнительно мы рассмотрели географическое распространение основных российских кластеров современной сублинии генотипа Beijing (B0/W148 и Central Asian/Russian) и двух основных кластеров древней сублинии, включая Beijing 14717-15-кластер на основании наших и ранее опубликованных данных [21-25] (рис. 4). Это сравнение показывает доминирование кластеров современной сублинии генотипа Beijing и относительно более широкое распространение кластера Beijing 1071-32 (древняя сублиния), хотя и в небольшой доле, за исключением Омской области (рис. 4).
Рис 4. Ретроспективный анализ географического распространения штаммов Beijing 14717-15 и других субтипов Beijing в Российских регионах. Показан процент в локальной популяции. CAR – Central Asian Russian.
Figure 4. Retrospective analysis of the geographical distribution of Beijing 14717-15 strains and other Beijing subtypes in Russian regions. The percentage in the local population is shown. CAR - Central Asian Russian.
Таблица 1. Детекция штаммов Beijing 14717-15 в региональных коллекциях
Table 1. Detection of strains Beijing 14717-15 in regional collections
Регион РФ или страна | Всего штаммов* | Beijing** | Beijing 14717-15 кластер Количество штаммов и процент от всей коллекции |
Вологда | 82 | 51 | 0 |
Мурманск | 67 | 35 | 0 |
Карелия | 67 | 36 | 0 |
Калининград | 73 | 46 | 0 |
Омск | 208 | 138 | 4 (2%) |
Иркутская область | 198 | 143 | 4 (2%) |
Бурятия | 265 | 180 | 47 (18%) |
Забайкальский край | 62 | 41 | 5 (8%) |
Монголия | 143 | 107 | 0 |
Китай | 113 | 102 | 0 |
Казахстан | 97 | 71 | 0 |
Япония* | 87 | 87 | 0 |
Примечания.
*коллекция из Японии включала только штаммы генотипа Beijing.
**Генотип Beijing определяли на основании сполиготипирования или ВНТР-типирования. Штаммы кластера Beijing 14717-15 определяли на основании разработанного метода детекции специфического SNP. Большинство штаммов этого кластера имели профили 14717-15 по 24-локусной схеме MIRU-VNTRplus.org.
Выводы
Клинически значимые штаммы могут возникать среди менее заметных геновариантов, и примером может служить недавно обнаруженный гипервирулентный и высоколетальный кластер Beijing 14717-15. Нами разработан метод его детекции на основе выявления специфической мутации в гене Rv2161c с помощью ПЦР с последующей рестрикцией продукта ПЦР и электрофорезом в агарозном геле для дискриминации дикого и мутантного аллелей. Преимуществами предлагаемого способа являются быстрота, простая и однозначная интерпретация результатов, возможность выявления контаминации и микс-инфекции, пригодность для производительного анализа больших коллекций штаммов М. tuberculosis.
Следует отметить, что штаммы Beijing 14717-15 входят в крупную ветвь ранней древней сублинии генотипа Beijing с интактным локусом RD181 (Рис. 1), которая весьма гетерогенна и включает штаммы с разнообразными профилями по VNTR. В этом смысле, выявленный нами маркер является маркером только специфического для России кластера Beijing, но не маркером всей ветви RD181-intact.
Российские изоляты кластера 14717-15, включенные в данное исследование, были выделены от больных в азиатской части России. Анализ большего количества штаммов (геномов) из Восточной Азии позволить провести более надежную реконструкцию эволюционной истории этого клинически значимого геноварианта M. tuberculosis, в то время как дальнейшие омиксные исследования дадут более информированное понимание его патобиологии.
About the authors
I. Mokrousov
Санкт-Петербургский научно-исследовательский институт эпидемиологии и микробиологии им. Пастера
Email: igormokrousov@yahoo.com
A. Vyazovaya
Email: annavyazovaya@gmail.com
M. Badleeva
Email: mbadleeva@mail.ru
A. Gerasimova
Email: kantarelle@mail.ru
O. Belopolskaya
Email: o.belopolskaya@spbu.ru
A. Masharsky
Email: alexey.masharsky@spbu.ru
I. Kostyukova
Email: i.v.yarusova@mail.ru
S. Zhdanova
Email: svetnii73@gmail.com
R. Mudarisova
Email: mudarisova.regina@bk.ru
I. Avadenii
Email: avadanii2013@yandex.ru
N. Solovieva
Email: baclab@spbniif.ru
Yu. Skiba
Email: yuriy.skiba@gmail.com
V. Zhuravlev
Email: jouravlev-slava@mail.ru
O. Pasechnik
Email: opasechnik@mail.ru
O. Ogarkov
Author for correspondence.
Email: obogarkov@gmail.com
References
- Gagneux S. Ecology and evolution of Mycobacterium tuberculosis. Nat Rev Microbiol. 2018; 16(4):202-213.
- Luo T, Comas I, Luo D, Lu B, Wu J, Wei L, et al. Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese. Proc Natl Acad Sci USA. 2015;112:8136-41..
- Mokrousov I, Ly HM, Otten T, Lan NN, Vyshnevskyi B, Hoffner S, et al. Origin and primary dispersal of the Mycobacterium tuberculosis Beijing genotype: clues from human phylogeography. Genome Res. 2005;15:1357-1364.
- Mokrousov I, Narvskaya O, Otten T, Vyazovaya A, Limeschenko E, Steklova L, et al. Phylogenetic reconstruction within Mycobacterium tuberculosis Beijing genotype in northwestern Russia. Res. Microbiol. 2002;153: 629–637.
- Kremer K, Glynn JR, Lillebaek T, Niemann S, Kurepina NE, Kreiswirth BN, et al. Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers. J Clin Microbiol. 2004;42:4040-4049.
- Maeda S, Hang NT, Lien LT, Thuong PH, Hung NV, Hoang NP, et al. Mycobacterium tuberculosis strains spreading in Hanoi, Vietnam: Beijing sublineages, genotypes, drug susceptibility patterns, and host factors. Tuberculosis (Edinb). 2014;94:649-656.
- Mokrousov I, Vyazovaya A, Pasechnik O, Gerasimova A, Dymova M, Chernyaeva E, et al. Early ancient sublineages of Mycobacterium tuberculosis Beijing genotype: unexpected clues from phylogenomics of the pathogen and human history. Clin Microbiol Infect. 2019;25:1039.e1-1039.e6.
- Mokrousov I, Vyazovaya A, Sinkov V, Gerasimova A, Ioannidis P, Jiao W, Khromova P, Papaventsis D, Pasechnik O, Perdigão J, Rastogi N, Shen A, Skiba Y, Solovieva N, Suffys P, Tafaj S, Umpeleva T, Vakhrusheva D, Yarusova I, Zhdanova S, Zhuravlev V, Ogarkov O. Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster. Sci Rep. 2021;11(1):21392.
- Mokrousov I, Sinkov V, Vyazovaya A, Pasechnik O, Solovieva N, Khromova P, Zhuravlev V, Ogarkov O. Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia. Int J Antimicrob Agents. 2020;56(2):106036.
- Vinogradova T., Dogonadze M, Zabolotnykh N, Badleeva M, Yarusova I, Vyazovaya A, Gerasimova A, Zhdanova S, Vitovskaya M, Solovieva N, Pasechnik O, Ogarkov O, Mokrousov I. Extremely lethal and hypervirulent Mycobacterium tuberculosis strain cluster emerging in Far East, Russia. Emerg Microbes Infect. 2021;10(1):1691-1701.
- Bespyatykh JA, Vinogradova ТI, Manicheva OA, et al. In vivo virulence of Beijing genotype Mycobacterium tuberculosis. Russ J Infect Immun. 2019;9:173-182. (In Russian).
- Mokrousov I, Pasechnik O, Vyazovaya A, Yarusova I, Gerasimova A, Blokh A, Zhuravlev V. Impact of pathobiological diversity of Mycobacterium tuberculosis on clinical features and lethal outcome of tuberculosis. BMC Microbiol. 2022;22(1):50.
- van Embden J, Cave M, Crawford J, Dale JW, Eisenach KD, Gicquel B, et al. Strain identification on Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J. Clin. Microbiol. 1993;31:406-409.
- Kamerbeek J, L Schouls, A Kolk, M van Agterveld, D van Soolingen, S Kuijper, A Bunschoten, H Molhuizen, R Shaw, M Goyal, J van Embden. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology J Clin Microbiol. 1997;35(4):907-14.
- Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rüsch-Gerdes S, Willery E, et al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006;44:4498-4510.
- Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35(6):1547-1549.
- Жданова С.Н. Молекулярно-эпидемиологическое исследование распространения основных генотипов возбудителя туберкулеза в Северной Азии: дис. … д-ра мед. наук / С.Н. Жданова. - Иркутск, 2018. - 276 с.
- . Shamputa IC, Lee J, Allix-Béguec C, Cho EJ, Lee JI, Rajan V, et al. Genetic diversity of Mycobacterium tuberculosis isolates from a tertiary care tuberculosis hospital in South Korea. J Clin Microbiol. 2010;48:3873-3894.
- Wada T, Iwamoto T, Maeda S. Genetic diversity of the Mycobacterium tuberculosis Beijing family in East Asia revealed through refined population structure analysis. FEMS Microbiol Lett, 2009;291:35-43.
- Yin QQ, Liu HC, Jiao WW, Li QJ, Han R, Tian JL, et al. Evolutionary History and Ongoing Transmission of Phylogenetic Sublineages of Mycobacterium tuberculosis Beijing Genotype in China. Sci Rep. 2016;6:34353.
- Pasechnik O, Vyazovaya A, Vitriv S, Tatarintseva M, Blokh A, Stasenko V, Mokrousov I. Major genotype families and epidemic clones of Mycobacterium tuberculosis in Omsk region, Western Siberia, Russia, marked by a high burden of tuberculosis-HIV coinfection. Tuberculosis (Edinb). 2018;108:163-168.
- Zhdanova S, Mokrousov I, Orlova E, Sinkov V, Ogarkov O. Transborder molecular analysis of drug-resistant tuberculosis in Mongolia and Eastern Siberia, Russia. Transbound Emerg Dis. 2022 69(5):e1800-e1814.
- Mokrousov I, Vyazovaya A, Solovieva N, Sunchalina T, Markelov Y, Chernyaeva E, Melnikova N, Dogonadze M, Starkova D, Vasilieva N, Gerasimova A, Kononenko Y, Zhuravlev V, Narvskaya O. Trends in molecular epidemiology of drug-resistant tuberculosis in Republic of Karelia, Russian Federation. BMC Microbiol. 2015;15:279.
- Вязовая А.А., Лебедева И.А., Ушакова Н.Б., и др. Молекулярно-генетический анализ популяции Mycobacterium tuberculosis в Вологодской области — регионе с низкой заболеваемостью туберкулезом // Инфекция и иммунитет. - 2021. - Т. 11. - №3. - C. 497-505.
- Вязовая А.А., Гаврилова Н.Ю., Герасимова А.А., Бычкова А.О., Авадений И., Аникиева Е.В., Соловьева Н.С., Журавлев В.Ю., Мокроусов И.В., Нарвская О.В. Молекулярно-генетический мониторинг популяции Mycobacterium tuberculosis в Мурманской области // Молекулярная генетика, микробиология и вирусология. – 2022. - №40(2). – С.21–27.
Supplementary files
There are no supplementary files to display.
