Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster.

BACKGROUND: . The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern. RESULTS: . The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes). CONCLUSIONS: . We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

Insight into Population Structure and Drug Resistance of Pediatric Tuberculosis Strains from China and Russia Gained through Whole-Genome Sequencing.

This study aimed to determine phenotypic and genotypic drug resistance patterns of Mycobacterium tuberculosis strains from children with tuberculosis (TB) in China and Russia, two high-burden countries for multi/extensively-drug resistant (MDR/XDR) TB. Whole-genome sequencing data of M. tuberculosis isolates from China (n = 137) and Russia (n = 60) were analyzed for phylogenetic markers and drug-resistance mutations, followed by comparison with phenotypic susceptibility data. The Beijing genotype was detected in 126 Chinese and 50 Russian isolates. The Euro-American lineage was detected in 10 Russian and 11 Chinese isolates. In the Russian collection, the Beijing genotype and Beijing B0/W148-cluster were dominated by MDR strains (68% and 94%, respectively). Ninety percent of B0/W148 strains were phenotypically pre-XDR. In the Chinese collection, neither of the Beijing sublineages was associated with MDR/pre-XDR status. MDR was mostly caused by low fitness cost mutations (rpoB S450L, katG S315T, rpsL K43R). Chinese rifampicin-resistant strains demonstrated a higher diversity of resistance mutations than Russian isolates (p = 0.003). The rifampicin and isoniazid resistance compensatory mutations were detected in some MDR strains, but they were not widespread. The molecular mechanisms of M. tuberculosis adaptation to anti-TB treatment are not unique to the pediatric strains, but they reflect the general situation with TB in Russia and China.

Impact of pathobiological diversity of Mycobacterium tuberculosis on clinical features and lethal outcome of tuberculosis.

BACKGROUND: Mycobacterium tuberculosis population in Russia is dominated by the notorious Beijing genotype whose major variants are characterized by contrasting resistance and virulence properties. Here we studied how these strain features could impact the progression of pulmonary tuberculosis (TB) concerning clinical manifestation and lethal outcome. RESULTS: The study sample included 548 M. tuberculosis isolates from 548 patients with newly diagnosed pulmonary TB in Omsk, West Siberia, Russia. Strains were subjected to drug susceptibility testing and genotyping to detect lineages, sublineages, and subtypes (within Beijing genotype). The Beijing genotype was detected in 370 (67.5%) of the studied strains. The strongest association with multidrug resistance (MDR) was found for epidemic cluster Beijing B0/W148 (modern sublineage) and two recently discovered MDR clusters 1071-32 and 14717-15 of the ancient Beijing sublineage. The group of patients infected with hypervirulent and highly lethal (in a mouse model) Beijing 14717-15 showed the highest rate of lethal outcome (58.3%) compared to Beijing B0/W148 (31.4%; P = 0.06), Beijing Central Asian/Russian (29.7%, P = 0.037), and non-Beijing (15.2%, P = 0.001). The 14717-15 cluster mostly included isolates from patients with infiltrative but not with fibrous-cavernous and disseminated TB. In contrast, a group infected with low virulent 1071-32-cluster had the highest rate of fibrous-cavernous TB, possibly reflecting the capacity of these strains for prolonged survival and chronicity of the TB process. CONCLUSIONS: The group of patients infected with hypervirulent and highly lethal in murine model 14717-15 cluster had the highest proportion of the lethal outcome (58.3%) compared to the groups infected with Beijing B0/W148 (31.4%) and non-Beijing (15.2%) isolates. This study carried out in the TB high-burden area highlights that not only drug resistance but also strain virulence should be considered in the implementation of personalized TB treatment.

First insight into the whole-genome sequence variations in Mycobacterium bovis BCG-1 (Russia) vaccine seed lots and their progeny clinical isolates from children with BCG-induced adverse events.

BACKGROUND: The only licensed live Bacille Calmette-Guérin (BCG) vaccine used to prevent severe childhood tuberculosis comprises genetically divergent strains with variable protective efficacy and rates of BCG-induced adverse events. The whole-genome sequencing (WGS) allowed evaluating the genome stability of BCG strains and the impact of spontaneous heterogeneity in seed and commercial lots on the efficacy of BCG-vaccines in different countries. Our study aimed to assess sequence variations and their putative effects on genes and protein functions in the BCG-1 (Russia) seed lots compared to their progeny isolates available from immunocompetent children with BCG-induced disease (mainly, osteitis). RESULTS: Based on the WGS data, we analyzed the links between seed lots 361, 367, and 368 used for vaccine manufacture in Russia in different periods, and their nine progeny isolates recovered from immunocompetent children with BCG-induced disease. The complete catalog of variants in genes relative to the reference genome (GenBank: CP013741) included 4 synonymous and 8 nonsynonymous single nucleotide polymorphisms, and 3 frameshift deletions. Seed lot 361 shared variants with 2 of 6 descendant isolates that had higher proportions of such polymorphisms in several genes, including ppsC, eccD5, and eccA5 involved in metabolism and cell wall processes and reportedly associated with virulence in mycobacteria. One isolate preserved variants of its parent seed lot 361 without gain of further changes in the sequence profile within 14 years. CONCLUSIONS: The background genomic information allowed us for the first time to follow the BCG diversity starting from the freeze-dried seed lots to descendant clinical isolates. Sequence variations in several genes of seed lot 361 did not alter the genomic stability and viability of the vaccine and appeared accumulated in isolates during the survival in the human organism. The impact of the observed variations in the context of association with the development of BCG-induced disease should be evaluated in parallel with the immune status and host genetics. Comparative genomic studies of BCG seed lots and their descendant clinical isolates represent a beneficial approach to better understand the molecular bases of efficacy and adverse events during the long-term survival of BCG in the host organism.

Mycobacterium tuberculosis RD-Rio Strain in Kazakhstan.

Mycobacterium tuberculosis RD-Rio strains are still rare in the former Soviet Union countries and Asia. We describe a strain in Kazakhstan that belongs to the RD-Rio secondary branch, which is endemic to northwest Russia and eastern Europe. Although RD-Rio strains are frequently multidrug resistant, this heterogeneous branch included only drug-susceptible isolates.

Simple Assay for Detection of the Central Asia Outbreak Clade of the Mycobacterium tuberculosis Beijing Genotype.

The Central Asia outbreak (CAO) clade is a branch of the Mycobacterium tuberculosis Beijing genotype that is associated with multidrug resistance, increased transmissibility, and epidemic spread in parts of the former Soviet Union. Furthermore, migration flows bring these strains far beyond their areas of origin. We aimed to find a specific molecular marker of the Beijing CAO clade and develop a simple and affordable method for its detection. Based on the bioinformatics analysis of the large M. tuberculosis whole-genome sequencing (WGS) data set (n = 1,398), we identified an IS6110 insertion in the Rv1359-Rv1360 intergenic region as a specific molecular marker of the CAO clade. We further designed and optimized a multiplex PCR method to detect this insertion. The method was validated in silico with the recently published WGS data set from Central Asia (n = 277) and experimentally with M. tuberculosis isolates from European and Asian parts of Russia, the former Soviet Union, and East Asia (n = 319). The developed molecular assay may be recommended for rapid screening of retrospective collections and for prospective surveillance when comprehensive but expensive WGS is not available or practical. The assay may be especially useful in high multidrug-resistant tuberculosis (MDR-TB) burden countries of the former Soviet Union and in countries with respective immigrant communities.

Whole-Genome Analysis of Mycobacterium tuberculosis from Patients with Tuberculous Spondylitis, Russia.

Whole-genome analysis of Mycobacterium tuberculosis isolates collected in Russia (N = 71) from patients with tuberculous spondylitis supports a detailed characterization of pathogen strain distributions and drug resistance phenotype, plus distinguished occurrence and association of known resistance mutations. We identify known and novel genome determinants related to bacterial virulence, pathogenicity, and drug resistance.

Emerging resistant clones of Mycobacterium tuberculosis in a spatiotemporal context.

Objectives: We assessed the genetic structure of the Mycobacterium tuberculosis population in Estonia with a special focus on major epidemic/endemic clones and drug resistance determinants. We investigated the hypothesis of the decisive impact of massive human influx on the locally circulating genotypes. Estonia received a mass immigration from Russia during 1945-90 followed by enhanced interaction with the EU since 1991. Methods: The study sample included M. tuberculosis isolates from patients newly diagnosed with TB in 2014 in North Estonia (including the capital Tallinn). The isolates were subjected to first- and second-line drug susceptibility testing, detection of mutations in rpoB, katG, inhA, rrs, embB and gyrA and lineage/clone-specific genotyping. Results: Of the M. tuberculosis isolates, 39.8% were assigned to the Beijing genotype; 56.8% of them were MDR. In contrast, all three major non-Beijing genotypes (LAM, Haarlem and Ural) were mainly drug susceptible. MDR was more prevalent among Beijing B0/W148-cluster isolates (81.8%) compared with other Beijing isolates (20.0%; P = 0.0007). The pre-XDR phenotype was found in eight isolates, of which six belonged to Beijing B0/W148. All rifampicin-resistant and ofloxacin-resistant and 97% of isoniazid-resistant isolates harboured resistance mutations in rpoB, gyrA and katG. The rpoB S531L, katG S315T and embB M306V mutations were the most prevalent. Conclusions: The major pool of the Beijing isolates was brought to Estonia before 1990. However, an active circulation of the most hazardous MDR-associated Beijing B0/W148-cluster started only in the last 20 years and its significantly increased circulation presents the major threat to TB control in Estonia. The overwhelming prevalence of the rpoB531 and katG315 mutations in the MDR-associated Beijing isolates requires attention.

Emerging peak on the phylogeographic landscape of Mycobacterium tuberculosis in West Asia: Definitely smoke, likely fire.

To date, a major attention was justly given to the global lineages of Mycobacterium tuberculosis. Here, we demonstrated an importance of the minor ones, on an example of intriguing and underestimated NEW-1 family that belongs to Euro-American lineage (lineage 4). Analysis of the global WGS/NGS datasets (5715 strains) identified 2235 strains of Lineage 4 and 66 strains of sublineage L4.5. This latter is marked with RD122 genomic deletion and includes NEW-1 family. Phylogenomic analysis confirmed a separate position of the NEW-1 family that we tentatively designate L4.5.1/Iran. We propose an evolution/migration scenario starting with origin of L4.5 1000-1300 ya in China, subsequent origin of the pre-NEW-1 intermediate genotype in Tibet, further migration to Xinjiang/Uyghur, and finally to Iran since 800 ya (origin of NEW-1), possibly, via expansion of the Mongol Yuan empire. Analysis of longitudinal phylogeographic datasets revealed a sharp increase in prevalence of NEW-1 strains in Iran and its eastwards neighbors in the last 20years; most alarmingly, it is accompanied with significant association with multidrug resistance (MDR). Ongoing migration, especially, Afghan refugees flows to developed countries emphasize a risk of the wider spread of the epidemic MDR subtype within NEW-1 family that we coin as emerging resistant clone of M. tuberculosis in West Asia.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography.

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

Tuberculous spondylitis in Russia and prominent role of multidrug-resistant clone Mycobacterium tuberculosis Beijing B0/W148.

Extrapulmonary and, in particular, spinal tuberculosis (TB) constitutes a minor but significant part of the total TB incidence. In spite of this, almost no studies on the genetic diversity and drug resistance of Mycobacterium tuberculosis isolates from spinal TB patients have been published to date. Here, we report results of the first Russian and globally largest molecular study of M. tuberculosis isolates recovered from patients with tuberculous spondylitis (TBS). The majority of 107 isolates were assigned to the Beijing genotype (n = 80); the other main families were T (n = 11), Ural (n = 7), and LAM (n = 4). Multidrug resistance (MDR) was more frequently found among Beijing (90.5%) and, intriguingly, Ural (71.4%) isolates than other genotypes (5%; P < 0.001). The extremely drug-resistant (XDR) phenotype was exclusively found in the Beijing isolates (n = 7). A notable prevalence of the rpoB531 and katG315 mutations in Beijing strains that were similarly high in both TBS (this study) and published pulmonary TB (PTB) samples from Russia shows that TBS and PTB Beijing strains follow the same paradigm of acquisition of rifampin (RIF) and isoniazid (INH) resistance. The 24-locus mycobacterial interspersed repetitive unit-variable-number tandem-repeat (MIRU-VNTR) subtyping of 80 Beijing isolates further discriminated them into 24 types (Hunter Gaston index [HGI] = 0.83); types 100-32 and 94-32 represented the largest groups. A genotype of Russian successful clone B0/W148 was identified in 30 of 80 Beijing isolates. In conclusion, this study highlighted a crucial impact of the Beijing genotype and the especially prominent role of its MDR-associated successful clone B0/W148 cluster in the development of spinal MDR-TB in Russian patients.

Trends in molecular epidemiology of drug-resistant tuberculosis in Republic of Karelia, Russian Federation.

BACKGROUND: Russian Republic of Karelia is located at the Russian-Finnish border. It contains most of the historical Karelia land inhabited with autochthonous Karels and more recently migrated Russians. Although tuberculosis (TB) incidence in Karelia is decreasing, it remains high (45.8/100 000 in 2014) with the rate of multi-drug resistance (MDR) among newly diagnosed TB patients reaching 46.5 %. The study aimed to genetically characterize Mycobacterium tuberculosis isolates obtained at different time points from TB patients from Karelia to gain insight into the phylogeographic specificity of the circulating genotypes and to assess trends in evolution of drug resistant subpopulations. METHODS: The sample included 150 M. tuberculosis isolates: 78 isolated in 2013-2014 ("new" collection) and 72 isolated in 2006 ("old" collection). Drug susceptibility testing was done by the method of absolute concentrations. Spoligotyping was used to test genotype-specific markers of a Latin-American-Mediterranean (LAM) family and its sublineages as well as a Beijing B0/W148-cluster. RESULTS: The largest spoligotypes were SIT1 (Beijing family, n = 42) and SIT40 (T family, n = 5). Beijing family was the largest (n = 43) followed by T (n = 11), Ural (n = 10) and LAM (n = 8). Successful Russian clone, Beijing В0/W148, was identified in 15 (34.9 %) of 43 Beijing isolates; all В0/W148 isolates were drug-resistant. Seven of 8 LAM isolates belonged to the RD115/LAM-RUS branch, 1 - to the LAM RD174/RD-Rio sublineage. MDR was found in Beijing (32/43), Ural (3/10), and LAM (3/8). In contrast, all T isolates were pansusceptible. Comparison of drug resistant subgroups of the new and old collections showed an increasing prevalence of the B0/W148 clonal cluster, from 18.0 % (mainly polyresistant) in 2006 to 32.6 % in 2014 (mainly MDR and pre-XDR). The West-east increasing gradient is observed for the Ural genotype that may be defined a 'Russian' strain. In contrast, the spoligotype SIT40 of the T family appears to be a historical Karelian strain. CONCLUSIONS: Circulation of the MDR M. tuberculosis isolates of the Beijing genotype and its B0/W148 cluster continues to critically influence the current situation with the MDR-TB control in northwestern Russia including the Republic of Karelia. Revealed phylogeographic patterns of some genotypes reflect a complex demographic history of Karelia within the course of the 20(th) century.

Mycobacterium tuberculosis Latin American-Mediterranean family and its sublineages in the light of robust evolutionary markers.

Mycobacterium tuberculosis has a clonal population structure, and the Latin American-Mediterranean (LAM) family is one of the largest and most widespread within this species, showing evidence for remarkable pathobiology and a confusing phylogeny. Here, we applied robust phylogenetic markers to study the evolution of the LAM family and its major sublineages circulating in Russia and neighboring countries. A total of 250 M. tuberculosis isolates were confirmed to belong to the LAM family based on the analysis of the LAM-specific single-nucleotide polymorphisms (SNPs) in the Rv3062 and Rv0129c genes. At this stage, the family status was rectified for 121 isolates misleadingly assigned by CRISPR spoligotyping to non-LAM families (T1- or T5-RUS1). Consequently, the reestimated LAM prevalence rate increased 2-fold in Russia and Kazakhstan and 4-fold in Belarus. The majority (91.8 to 98.7%) of the LAM isolates from all three countries belonged to the LAM-RUS sublineage. In contrast, the Ibero-American LAM RD-Rio sublineage was identified in only 7 Russian isolates. Taken together, our findings and further analyses suggest a monophyletic origin of LAM-RUS: at a historically distant time, in Russia, in a small founding bacterial/human population. Its dissemination pattern and high prevalence rate in Northern Eurasia may indicate a long-term coexistence of the LAM-RUS sublineage and local human populations hypothetically leading to coadaptation and reduced pathogenicity of the relatively more ancient clones, such as spoligotype international type 254 (SIT254), compared to the more recent SIT252 and SIT266 clones. In contrast, rare LAM RD-Rio isolates were likely brought to Russia through occasional human contact. The spread of RD-Rio strains is not as global as commonly claimed and is determined largely by human migration flows (rather than by pathobiological properties of these strains). Consequently, a host population factor appears to play a major role in shaping the in situ dissemination pattern of the imported strains in an autochthonous population.

Genome-wide Mycobacterium tuberculosis variation (GMTV) database: a new tool for integrating sequence variations and epidemiology.

BACKGROUND: Tuberculosis (TB) poses a worldwide threat due to advancing multidrug-resistant strains and deadly co-infections with Human immunodeficiency virus. Today large amounts of Mycobacterium tuberculosis whole genome sequencing data are being assessed broadly and yet there exists no comprehensive online resource that connects M. tuberculosis genome variants with geographic origin, with drug resistance or with clinical outcome. DESCRIPTION: Here we describe a broadly inclusive unifying Genome-wide Mycobacterium tuberculosis Variation (GMTV) database, (http://mtb.dobzhanskycenter.org) that catalogues genome variations of M. tuberculosis strains collected across Russia. GMTV contains a broad spectrum of data derived from different sources and related to M. tuberculosis molecular biology, epidemiology, TB clinical outcome, year and place of isolation, drug resistance profiles and displays the variants across the genome using a dedicated genome browser. GMTV database, which includes 1084 genomes and over 69,000 SNP or Indel variants, can be queried about M. tuberculosis genome variation and putative associations with drug resistance, geographical origin, and clinical stages and outcomes. CONCLUSIONS: Implementation of GMTV tracks the pattern of changes of M. tuberculosis strains in different geographical areas, facilitates disease gene discoveries associated with drug resistance or different clinical sequelae, and automates comparative genomic analyses among M. tuberculosis strains.

Real-time PCR assay for rapid detection of epidemiologically and clinically significant Mycobacterium tuberculosis Beijing genotype isolates.

Mycobacterium tuberculosis Beijing genotype strains are rapidly disseminating, frequently hypervirulent, and multidrug resistant. Here, we describe a method for their rapid detection by real-time PCR that targets the specific IS6110 insertion in the dnaA-dnaN genome region. The method was evaluated with a geographically and genetically diverse collection representing areas in East Asia and the former Soviet Union in which the Beijing genotype is endemic and epidemic (i.e., major foci of its global propagation) and with clinical specimens.

Genetic Diversity and Primary Drug Resistance of Mycobacterium tuberculosis Beijing Genotype Strains in Northwestern Russia.

The Beijing genotype is the main family of Mycobacterium tuberculosis in Russia. We analyzed its diversity and drug resistance in provinces across Northwestern Russia to identify the epidemiologically relevant Beijing strains. The study collection included 497 isolates from newly-diagnosed tuberculosis (TB) patients. Bacterial isolates were subjected to drug-susceptibility testing and genotyping. The Beijing genotype was detected in 57.5% (286/497); 50% of the Beijing strains were multidrug-resistant (MDR). Central Asian/Russian and B0/W148 groups included 176 and 77 isolates, respectively. MDR was more frequent among B0/W148 strains compared to Central Asian/Russian strains (85.7% vs. 40.3%, p < 0.0001). Typing of 24 minisatellite loci of Beijing strains revealed 82 profiles; 230 isolates were in 23 clusters. The largest Central Asian/Russian types were 94-32 (n = 75), 1065-32 (n = 17), and 95-32 (n = 12). B0/W148 types were 100-32 (n = 59) and 4737-32 (n = 5). MDR was more frequent in types 1065-32 (88.2%), 100-32 (83.1%), and 4737-32 (100%). In contrast, type 9391-32 (n = 9) included only drug-susceptible strains. To conclude, M. tuberculosis Beijing genotype is dominant in Northwestern Russia, and an active transmission of overwhelmingly MDR B0/W148 types explains the reported increase of MDR-TB. The presence of MDR-associated minor variants (type 1071-32/ancient Beijing and Central Asia Outbreak strain) in some of the studied provinces also requires attention.

A multifaceted interplay between virulence, drug resistance, and the phylogeographic landscape of Mycobacterium tuberculosis.

Latin-American Mediterranean (LAM) family is one of the most significant and global genotypes of Mycobacterium tuberculosis. Here, we used the murine model to study the virulence and lethality of the genetically and epidemiologically distinct LAM strains. The pathobiological characteristics of the four LAM strains (three drug resistant and one drug susceptible) and the susceptible reference strain H37Rv were studied in the C57BL/6 mouse model. The whole-genome sequencing was performed using the HiSeq Illumina platform, followed by bioinformatics and phylogenetic analysis. The susceptible strain H37Rv showed the highest virulence. Drug-susceptible LAM strain (spoligotype SIT264) was more virulent than three multidrug-resistant (MDR) strains (SIT252, SIT254, and SIT266). All three MDR isolates were low lethal, while the susceptible isolate and H37Rv were moderately/highly lethal. Putting the genomic, phenotypic, and virulence features of the LAM strains/spoligotypes in the context of their dynamic phylogeography over 20 years reveals three types of relationships between virulence, resistance, and transmission. First, the most virulent and more lethal drug-susceptible SIT264 increased its circulation in parts of Russia. Second, moderately virulent and pre-XDR SIT266 was prevalent in Belarus and continues to be visible in North-West Russia. Third, the low virulent and MDR strain SIT252 previously considered as emerging has disappeared from the population. These findings suggest that strain virulence impacts the transmission, irrespective of drug resistance properties. The increasing circulation of susceptible but more virulent and lethal strains implies that personalized TB treatment should consider not only resistance but also the virulence of the infecting M. tuberculosis strains. IMPORTANCE The study is multidisciplinary and investigates the epidemically/clinically important and global lineage of Mycobacterium tuberculosis, named Latin-American-Mediterranean (LAM), yet insufficiently studied with regard to its pathobiology. We studied different LAM strains (epidemic vs endemic and resistant vs susceptible) in the murine model and using whole-genome analysis. We also collected long-term, 20-year data on their prevalence in Eurasia. The findings are both expected and unexpected. (i) We observe that a drug-susceptible but highly virulent strain increased its prevalence. (ii) By contrast, the multidrug-resistant (MDR) but low-virulent, low-lethal strain (that we considered as emerging 15 years ago) has almost disappeared. (iii) Finally, an intermediate case is the MDR strain with moderate virulence that continues to circulate. We conclude that (i) the former and latter strains are the most hazardous and require close epidemiological monitoring, and (ii) personalized TB treatment should consider not only drug resistance but also the virulence of the infecting strains and development of anti-virulence drugs is warranted.

Russian "successful" clone B0/W148 of Mycobacterium tuberculosis Beijing genotype: a multiplex PCR assay for rapid detection and global screening.

We describe a multiplex PCR assay to detect the Mycobacterium tuberculosis Beijing genotype variant B0/W148, which is considered a "successful" clone of M. tuberculosis, widespread in Russia. Specificity and sensitivity of the assay were 100% based on the analysis of a collection of 516 M. tuberculosis isolates of different genotypes and origins. This assay may be used for accurate and simple detection and surveillance of this clinically and epidemiologically important variant of M. tuberculosis.

Population structure of Mycobacterium tuberculosis from referral clinics in Western Siberia, Russia: Before and during the Covid-19 pandemic.

The dramatic change in global health imposed by the Covid-19 pandemic has also impacted TB control. The TB incidence decreased dramatically not because of the improved situation but due to undertesting, reduced resources, and ultimately, substantially reduced detection rate. We hypothesized that multiple and partly counteracting factors could influence changes in the local Mycobacterium tuberculosis population. To test this hypothesis, we analyzed M. tuberculosis isolates collected in Western Siberia, Russia, before and during the Covid-19 pandemic. A total of 269 M. tuberculosis isolates from patients admitted at referral clinics were studied. The pre-pandemic and pandemic collections included 179 and 90 isolates, respectively. Based on genotyping, both pre-pandemic and pandemic samples are heavily dominated by the Beijing genotype isolates (95% and 88%) that were mostly MDR (80 and 68%). The high proportion of MDR isolates is due to the specific features of the studied collections biased towards patients with severe TB admitted at the National referral center in Novosibirsk. While no dramatic change was observed in the M. tuberculosis population structure in the survey area in Western Siberia during the Covid-19 pandemic in 2020-2021 compared to the pre-pandemic collection, still we note a certain decrease of the Beijing genotype and an increase in the proportion and diversity of the non-Beijing isolates. However, the transmissible and MDR Beijing B0/W148 did not increase its prevalence rate during the pandemic. More generally, the high prevalence rate of the Beijing genotype and its strong association with MDR both before and during the pandemic are alarming features of this region in Western Siberia, Russia.