Search

Abstract 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.

ABSTRACT 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.

Abstract Ancient sublineage of the Mycobacterium tuberculosis Beijing genotype is endemic and prevalent in East Asia and rare in other world regions. While these strains are mainly drug susceptible, we recently identified a novel clonal group Beijing 1071-32 within this sublineage emerging in Siberia, Russia and present in other Russian regions. This cluster included only multi/extensive drug resistant (MDR/XDR) isolates. Based on the phylogenetic analysis of the available WGS data, we identified three synonymous SNPs in the genes Rv0144, Rv0373c, and Rv0334 that were specific for the Beijing 1071-32-cluster and developed a real-time PCR assay for their detection. Analysis of the 2375 genetically diverse M. tuberculosis isolates collected between 1996 and 2020 in different locations (European and Asian parts of Russia, former Soviet Union countries, Albania, Greece, China, Vietnam, Japan and Brazil), confirmed 100% specificity and sensitivity of this real-time PCR assay. Moreover, the epidemiological importance of this strain and the newly developed screening assay is further stressed by the fact that all identified Beijing 1071-32 isolates were found to exhibit MDR genotypic profiles with concomitant resistance to additional first-line drugs due to a characteristic signature of six mutations in rpoB450, rpoC485, katG315, katG335, rpsL43 and embB497. In conclusion, this study provides a set of three concordant SNPs for the detection and screening of Beijing 1071-32 isolates along with a validated real-time PCR assay easily deployable across multiple settings for the epidemiological tracking of this significant MDR cluster.

Mycobacterium tuberculosis strains of the early ancient sublineage of the Beijing genotype are mostly drug susceptible and mainly circulate in East Asia. We have recently discovered two clusters of this sublineage emerging in the Asian part of Russia (VNTR-defined 1071-32 and 14717-15 types) and, to our surprise, both were strongly MDR/XDR-associated. Here, we evaluated their pathogenic features. The clinical isolates and reference laboratory strain H37Rv were investigated in the C57BL/6 mouse model to assess their virulence and lethality properties. The BACTEC MGIT 960 system was used to study the in vitro growth characteristics. In the murine model, strains 396 (14717-15-cluster, from Buryatia, Far East) and 6691 (1071-32-cluster, from Omsk, Siberia) demonstrated contrasting properties. The 396-infected group had significantly higher mortality, more weight loss, higher bacterial burden, and more severe lung pathology. Furthermore, compared to the previously published data on other Russian epidemic Beijing strains (B0/W148, CAO, Central Asian Russian), strain 396 demonstrated the highest mortality. Under the in vitro growth experiment, cluster 14717-15 isolates had significantly shorter lag-phase. To conclude, low-virulent MDR strain 6691 belongs to the Beijing 1071-32-cluster widespread across FSU countries but at low prevalence. This corresponds to common expectation that multiple drug resistance mutations reduce fitness and virulence. In contrast, highly lethal and hypervirulent MDR strain 396 represents an intriguing Beijing 14717-15 cluster predominant only in Buryatia, Far East (16%), sporadically found beyond it, but not forming clusters of transmission. Further in-depth study of this most virulent Russian Beijing cluster is warranted.

Abstract 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.

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.

We performed synthesis of new nitrofuranyl amides and investigated their anti-TB activity and primary genetic response of mycobacteria through whole-genome sequencing (WGS) of spontaneous resistant mutants. The in vitro activity was assessed on reference strain Mycobacterium tuberculosis H37Rv. The most active compound 11 was used for in vitro selection of spontaneous resistant mutants. The same mutations in six genes were detected in bacterial cultures grown under increased concentrations of 11 (2×, 4×, 8× MIC). The mutant positions were presented as mixed wild type and mutant alleles while increasing the concentration of the compound led to the semi-proportional and significant increase in mutant alleles. The identified genes belong to different categories and pathways. Some of them were previously reported as mediating drug resistance or drug tolerance, and counteracting oxidative and nitrosative stress, in particular: Rv0224c, fbiC, iniA, and Rv1592c. Gene-set interaction analysis revealed a certain weak interaction for gene pairs Rv1592–Rv1639c and Rv1592–Rv0224c. To conclude, this study experimentally demonstrated a multifaceted primary genetic response of M. tuberculosis to the action of nitrofurans. All three 11-treated subcultures independently presented the same six SNPs, which suggests their non-random occurrence and likely causative relationship between compound action and possible resistance mechanism.

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.

Mycobacterium tuberculosis is a highly studied pathogen due to public health importance. Despite this, problems like early drug resistance, diagnostics and treatment success prediction are still not fully resolved. Here, we analyze the incidence of point mutations widely used for drug resistance detection in laboratory practice and conduct comparative analysis of whole-genome sequence (WGS) for clinical M. tuberculosis strains collected from patients with pulmonary tuberculosis (PTB) and extra-pulmonary tuberculosis (XPTB) localization. A total of 72 pulmonary and 73 extrapulmonary microbiologically characterized M. tuberculosis isolates were collected from patients from 2007 to 2014 in Russia. Genomic DNA was used for WGS and obtained data allowed identifying major mutations known to be associated with drug resistance to first-line and second-line antituberculous drugs. In some cases previously described mutations were not identified. Using genome-based phylogenetic analysis we identified M. tuberculosis substrains associated with distinctions in the occurrence in PTB vs. XPTB cases. Phylogenetic analyses did reveal M. tuberculosis genetic substrains associated with TB localization. XPTB was associated with Beijing sublineages Central Asia (Beijing CAO), Central Asia Clade A (Beijing A) and 4.8 groups, while PTB localization was associated with group LAM (4.3). Further, the XPTB strain in some cases showed elevated drug resistance patterns relative to PTB isolates. HIV was significantly associated with the development of XPTB in the Beijing B0/W148 group and among unclustered Beijing isolates.

Perchlozone ([PCZ] 4-thioureido-iminomethylpyridinium perchlorate) is a new thiosemicarbazone approved for the treatment of multidrug-resistant tuberculosis (MDR-TB) in Russia and some other countries. The ethA and hadABC mutations may confer PCZ resistance. At the same time, ethA mutations are known to mediate resistance to ethionamide (ETH) and prothionamide (PTH). We aimed to study the genetic variation underlying Mycobacterium tuberculosis resistance to PCZ through whole genome sequencing (WGS) of consecutive isolates recovered during long-term treatment. This prospective study included patients admitted in 2018–2019 to the regional tuberculosis dispensary, Kaliningrad, Russia, whose treatment regimen included PCZ. Multiple M. tuberculosis isolates were recovered during PCZ treatment, and the bacterial DNA was subjected to WGS followed by bioinformatics analysis. We identified mutations in the genes putatively associated with PCZ resistance, ethA, and hadA. The most frequent one was a frameshift ethA 106 GA > G (seven of nine patients) and most of the other mutations were also likely present before PCZ treatment. In one patient, a frameshift mutation ethA 702 CT > C emerged after six months of PCZ treatment. A frequent presence of cross-resistance mutations to PCZ and ETH/PTH should be taken into consideration when PCZ is included in the treatment regimen of MDR-TB patients.