Insights into the origin, emergence, and current spread of a successful Russian clone of Mycobacterium tuberculosis

Increasing circulation of multi-drug resistant tuberculosis strains in Buryatia, high-burden and ethnically diverse region in the Russian Far East

Buryatia is a multidrug-resistant tuberculosis (MDR-TB) high-burden region in the Russian Far East with ethnically diverse population (30 % Mongoloid Buryats and 65 % Russians). Two hundred M. tuberculosis strains from newly-diagnosed patients were subjected to phenotypic testing and genotyping. The Beijing genotype was more prevalent among Russians than Buryats (68 % vs 53 %; P = 0.055). European non-Beijing genotypes (LAM, Ural, Haarlem) were double more prevalent in Buryats vs Russians (39.2 % vs 20.5 %; P = 0.01). Higher prevalence of Beijing among former prison inmates (79 % vs 61 % in other patients, P = 0.1) suggests its increased transmissibility. The Russian epidemic cluster B0/W148 was in 9.5 %, double smaller than elsewhere in Siberia. The hypervirulent Beijing 14717-15-cluster was endemic in Buryatia but paradoxically enough, it was more frequently isolated from Russians than Buryats (9.1 % vs 3.9 %; P = 0.2). Beijing subtypes B0/W148, CAO, and 14717-15 were associated with poly/multi-drug resistance (P = 0.01-0.0001). HIV coinfection was more frequent in Russians than in Buryats: 35/141 (24.8 %) vs 5/51 (9.8 %), P = 0.03. To conclude, M. tuberculosis population structure in Buryatia retained its singularities compared to other parts of Russia and remains strikingly different from the neighboring Mongolia. A circulation of strongly MDR-associated Beijing subtypes and drug-resistant non-Beijing strains highlights a risk of their broader dissemination.

Genomic Interactions Between Mycobacterium tuberculosis and Humans

Mycobacterium tuberculosis is considered by many to be the deadliest microbe, with the estimated annual cases numbering more than 10 million. The bacteria, including Mycobacterium africanum, are classified into nine major lineages and hundreds of sublineages, each with different geographical distributions and levels of virulence. The phylogeographic patterns can be a result of recent and early human migrations as well as coevolution between the bacteria and various human populations, which may explain why many studies on human genetic factors contributing to tuberculosis have not been replicable in different areas. Moreover, several studies have revealed the significance of interactions between human genetic variations and bacterial genotypes in determining the development of tuberculosis, suggesting coadaptation. The increased availability of whole-genome sequence data from both humans and bacteria has enabled a better understanding of these interactions, which can inform the development of vaccines and other control measures.

Genetic diversity, evolution and drug resistance of Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis causes a chronic infectious disease called tuberculosis. Phylogenetic lineage 2 (L2) of M. tuberculosis, also known as the East Asian lineage, is associated with high virulence, increased transmissibility, and the spread of multidrug-resistant strains. This review article examines the genomic characteristics of the M. tuberculosis genome and M. tuberculosis lineage 2, such as the unique insertion sequence and spoligotype patterns, as well as MIRU-VNTR typing, and SNP-based barcoding. The review describes the geographical distribution of lineage 2 and its history of origin. In addition, the article discusses recent studies on drug resistance and compensatory mechanisms of M. tuberculosis lineage 2 and its impact on the pathogen's transmissibility and virulence. This review article discusses the importance of establishing a unified classification for lineage 2 to ensure consistency in terminology and criteria across different studies and settings.

Characteristic SNPs defining the major multidrug-resistant Mycobacterium tuberculosis clusters identified by EuSeqMyTB to support routine surveillance, EU/EEA, 2017 to 2019

BackgroundThe EUSeqMyTB project, conducted in 2020, used whole genome sequencing (WGS) for surveillance of drug-resistant Mycobacterium tuberculosis in the European Union/European Economic Area (EU/EEA) and identified 56 internationally clustered multidrug-resistant (MDR) tuberculosis (TB) clones.AimWe aimed to define and establish a rapid and computationally simple screening method to identify probable members of the main cross-border MDR-TB clusters in WGS data to facilitate their identification and track their future spread.MethodsWe screened 34 of the larger cross-border clusters identified in the EuSeqMyTB pilot study (2017-19) for characteristic single nucleotide polymorphism (SNP) signatures that could identify and define members of each cluster. We also linked this analysis with published clusters identified in previous studies and identified more distant genetic relationships between some of the current clusters.ResultsA panel of 30 characteristic SNPs is presented that can be used as an initial (routine) screen for members of each cluster. For four of the clusters, no unique defining SNP could be identified; three of these are closely related (within approximately 20 SNPs) to one or more other clusters and likely represent a single established MDR-TB clade composed of multiple recent subclusters derived from the previously described ECDC0002 cluster.ConclusionThe identified SNP signatures can be integrated into routine pipelines and contribute to the more effective monitoring, rapid and widespread screening for TB. This SNP panel will also support accurate communication between laboratories about previously identified internationally transmitted MDR-TB genotypes.

Prevalence of Beijing Central Asian/Russian Cluster 94-32 among Multidrug-Resistant M. tuberculosis in Kazakhstan

The Beijing genotype is the most distributed M. tuberculosis family in Kazakhstan. In this study, we identified dominant Beijing clusters in Kazakhstan and assessed their drug susceptibility profiles and association with the most widely spread mutation Ser531Leu of the rpoB gene and the mutation Ser315Thr of the katG gene associated with resistance to rifampicin and isoniazid, respectively. M. tuberculosis isolates (n = 540) from new TB cases were included in the study. MIRU-VNTR genotyping was performed for 540 clinical isolates to determine M. tuberculosis families using 24 loci. RD analysis was additionally performed for the Beijing isolates. The identification of mutations in the drug-resistance genes of M. tuberculosis was performed with allele-specific real-time PCR and Sanger sequencing. The Beijing genotype was identified in 60% (324/540) of the clinical isolates. Central Asian/Russian cluster 94-32 was the most distributed cluster among the Beijing isolates (50.3%; 163/324). Three other dominant Beijing clusters were identified as 94-33 (3.4%; 11/324), 100-32 (3.1%; 10/324) and 99-32 (3.1%; 10/324). The Beijing genotype was associated with drug-resistant TB (p < 0.0001), including multidrug-resistant TB (p < 0.0001), in our study. An association of the mutation Ser531Leu of the rpoB gene with the Beijing genotype was found (p < 0.0001; OR = 16.0000; 95%CI: 4.9161-52.0740). Among the Beijing isolates, cluster 94-32 showed an association with MDR-TB (p = 0.021). This is why the evaluation of the Beijing genotype and its clusters is needed to control MDR-TB in Kazakhstan.

The Global Success of Mycobacterium tuberculosis Modern Beijing Family Is Driven by a Few Recently Emerged Strains

Strains of the Mycobacterium tuberculosis complex (MTBC) Beijing family aroused concern because they were often found in clusters and appeared to be exceptionally transmissible. However, it was later found that strains of the Beijing family were heterogeneous, and the transmission advantage was restricted to sublineage L2.3 or modern Beijing. In this study, we analyzed the previously published genome sequences of 7,896 L2.3 strains from 51 different countries. Using BEAST software to approximate the temporal emergence of L2.3, our calculations suggest that L2.3 initially emerged in northern East Asia during the early 15th century and subsequently diverged into six phylogenetic clades, identified as L2.3.1 through L2.3.6. Using terminal branch length and genomic clustering as proxies for transmissibility, we found that the six clades displayed distinct population dynamics, with the three recently emerged clades (L2.3.4 to L2.3.6) exhibiting significantly higher transmissibility than the older three clades (L2.3.1 to L2.3.3). Of the Beijing family strains isolated outside East Asia, 83.1% belonged to the clades L2.3.4 to L2.3.6, which were also associated with more cross-border transmission. This work reveals the heterogeneity in sublineage L2.3 and demonstrates that the global success of Beijing family strains is driven by the three recently emerged L2.3 clades. IMPORTANCE The recent population dynamics of the global tuberculosis epidemic are heavily shaped by Mycobacterium tuberculosis complex (MTBC) strains with enhanced transmissibility. The infamous Beijing family strain stands out because it has rapidly spread throughout the world. Identifying the strains responsible for the global expansion and tracing their evolution should help to understand the nature of high transmissibility and develop effective strategies to control transmission. In this study, we found that the L2.3 sublineage diversified into six phylogenetic clades (L2.3.1 to L2.3.6) with various transmission characteristics. Clades L2.3.4 to L2.3.6 exhibited significantly higher transmissibility than clades L2.3.1 to L2.3.3, which helps explain why more than 80% of Beijing family strains collected outside East Asia belong to these three clades. We conclude that the global success of L2.3 was not caused by the entire L2.3 sublineage but rather was due to the rapid expansion of L2.3.4 to L2.3.6. Tracking the transmission of L2.3.4 to L2.3.6 strains can help to formulate targeted TB prevention and control.

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.

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.

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.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

Molecular characteristics of Mycobacterium tuberculosis drug-resistant isolates from HIV- and HIV+ tuberculosis patients in Russia

Background

High burden of drug-resistant (DR) tuberculosis (TB) is a significant threat to national TB control programs all over the world and in the Russian Federation. Different Mycobacterium tuberculosis (MTB) genotypes are hypothesized to have specific characteristics affecting TB control programs. For example, Beijing strains are supposed to have higher mutation rates compared to strains of other genotypes and subsequently higher capability to develop drug-resistance.

Results

Clinical MTB isolates from HIV- and HIV+ patients from four regions of Russia were analyzed for genotypes and mutations conferring resistance to Isoniazid, Rifampicin, Ethambutol, aminoglycosides, and fluoroquinolones. Analysis of genotypes and polymorphism of genomic loci according to the HIV status of the patients – sources of MTB isolates were performed. Studied MTB isolates from HIV- TB patients belonged to 15 genotypes and from HIV + TB patients – to 6 genotypes. Beijing clinical isolates dominated in HIV- (64,7%) and HIV+ (74,4%) groups. Other isolates were of LAM (including LAM1 and LAM9), Ural, and 4 minor groups of genotypes (including 5 subclones T). The spectrum of genotypes in the HIV- group was broader than in the HIV+ group. PR of B0/W148 Beijing was significantly lower than of other Beijing genotypes in susceptible and MDR-XDR isolates. Rates of isolates belonging to non-Beijing genotypes were higher than Beijing in susceptible isolates from HIV- patients.

Conclusions

Beijing genotype isolates prevailed in clinical isolates of all drug susceptibility profiles both from HIV- and HIV+ patients, although B0/W148 Beijing genotype did not dominate in this study. Genome loci and mutations polymorphisms were more pronounced in clinical isolates from HIV- patients, than from HIV+.

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.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster

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.

Genetic Profile of Linezolid-Resistant M. tuberculosis Clinical Strains from Moscow

Background: Linezolid, bedaquiline, and newer fluoroquinolones are currently placed as priority Group A drugs for the treatment of drug-resistant tuberculosis. The number of reported linezolid-resistant clinical strains is still low, and the correlation of molecular determinants with phenotype is not perfect. Methods: We determined the linezolid MICs for clinical isolates from the Moscow region and identified mutations in rplC and rrl genes. Results: All 16 linezolid-resistant isolates had previously reported mutations in the rplC or rrl loci, and 13 of them bore a RplC C154R substitution. Detection of this substitution in a heteroresistant state was not successful, probably, due to the more stable DNA secondary structure of the mutated fragment, which precludes its amplification in mixes with the wild-type DNA. Strains with an rplC mutation had higher linezolid MIC compared to isolates with rrl mutations. Conclusions: Linezolid resistance mostly emerged during treatment with the latest regimen. Three primary cases with linezolid resistance question the possible transmission of totally drug-resistant tuberculosis in the Moscow region, which demands further investigation.

MDR and Pre-XDR Clinical Mycobacterium tuberculosis Beijing Strains: Assessment of Virulence and Host Cytokine Response in Mice Infectious Model

Mycobacterium tuberculosis Beijing genotype associated with drug resistance is a growing public health problem worldwide. The aim of this study was the assessment of virulence for C57BL/6 mice after infection by clinical M. tuberculosis strains 267/47 and 120/26, which belong to the modern sublineages B0/W148 and Central Asia outbreak of the Beijing genotype, respectively. The sublineages were identified by the analysis of the strains' whole-genomes. The strains 267/47 and 120/26 were characterized as agents of pre-extensively drug-resistant (pre-XDR) and multidrug-resistant (MDR) tuberculosis, respectively. Both clinical strains were slow-growing in 7H9 broth compared to the M. tuberculosis H37Rv strain. The survival rates of C57BL/6 mice infected by 267/47, 120/26, and H37Rv on the 150th day postinfection were 10%, 40%, and 70%, respectively. Mycobacterial load in the lungs, spleen, and liver was higher and histopathological changes were more expressed for mice infected by the 267/47 strain compared to those infected by the 120/26 and H37Rv strains. The cytokine response in the lungs of C57BL/6 mice after infection with the 267/47, 120/26, and H37Rv strains was different. Notably, proinflammatory cytokine genes Il-1α, Il-6, Il-7, and Il-17, as well as anti-inflammatory genes Il-6 and Il-13, were downregulated after an infection caused by the 267/47 strain compared to those after infection with the H37Rv strain.

Genetic composition and evolution of the prevalent Mycobacterium tuberculosis lineages 2 and 4 in the Chinese and Zhejiang Province populations

Background

There are seven human-adaptation lineages of Mycobacterium tuberculosis (Mtb). Tuberculosis (TB) dissemination is strongly influenced by human movements and host genetics. The detailed lineage distribution evolution of Mtb in Zhejiang Province is unknown. We aim to determine how different sub-lineages are transmitted and distributed within China and Zhejiang Province.

Methods

We analysed whole-genome sequencing data for a worldwide collection of 1154 isolates and a provincial collection of 1296 isolates, constructed the best-scoring maximum likelihood phylogenetic tree. Bayesian evolutionary analysis was used to calculate the latest common ancestor of lineages 2 and 4. The antigenic diversity of human T cell epitopes was evaluated by calculating the pairwise dN/dS ratios.

Results

Of the Zhejiang isolates, 964 (74.38%) belonged to lineage 2 and 332 (25.62%) belonged to lineage 4. The distributions of the sub-lineages varied across the geographic regions of Zhejiang Province. L2.2 is the most ancient sub-lineage in Zhejiang, first appearing approximately 6897 years ago (95% highest posterior density interval (HDI): 6513–7298). L4.4 is the most modern sub-lineage, first appearing approximately 2217 years ago (95% HDI: 1864–2581). The dN/dS ratios showed that the epitope and non-epitope regions of lineage 2 strains were significantly (P < 0.001) more conserved than those of lineage 4.

Conclusions

An increase in the frequency of lineage 4 may reflect its successful transmission over the last 20 years. The recent common ancestors of the sub-lineages and their transmission routes are relevant to the entry of humans into China and Zhejiang Province. Diversity in T cell epitopes may prevent Mycobacterium tuberculosis from being recognized by the immune system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00673-7.

Genomic Variations in Drug Resistant Mycobacterium tuberculosis Strains Collected from Patients with Different Localization of Infection

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.

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

Genotyping, Assessment of Virulence and Antibacterial Resistance of the Rostov Strain of Mycobacterium tuberculosis Attributed to the Central Asia Outbreak Clade

The Central Asia Outbreak (CAO) clade is a growing public health problem for Central Asian countries. Members of the clade belong to the narrow branch of the Mycobacterium tuberculosis Beijing genotype and are characterized by multidrug resistance and increased transmissibility. The Rostov strain of M. tuberculosis isolated in Russia and attributed to the CAO clade based on PCR-assay and whole genome sequencing and the laboratory strain H37Rv were selected to evaluate the virulence on C57Bl/6 mice models by intravenous injection. All mice infected with the Rostov strain succumbed to death within a 48-day period, while more than half of the mice infected by the H37Rv strain survived within a 90-day period. Mice weight analysis revealed irreversible and severe depletion of animals infected with the Rostov strain compared to H37Rv. The histological investigation of lung and liver tissues of mice on the 30th day after injection of mycobacterial bacilli showed that the pattern of pathological changes generated by two strains were different. Moreover, bacterial load in the liver and lungs was higher for the Rostov strain infection. In conclusion, our data demonstrate that the drug-resistant Rostov strain exhibits a highly virulent phenotype which can be partly explained by the CAO-specific mutations.

Increased transmissibility of Russian successful strain Beijing B0/W148 of Mycobacterium tuberculosis: Indirect clues from history and demographics

The local situation with tuberculosis (TB) is shaped by the complex interplay of multiple factors related to both human host and Mycobacterium tuberculosis. We hypothesized that TB epidemiology in the rural regions in the Soviet Union was impacted by construction of the Gulag camps and significant incoming migration. This molecular M. tuberculosis study was conducted in 2017 in the Komi Republic in northern Russia, a region with high rate (26%) of primary multidrug-resistant (MDR) TB. MDR was detected in 30.8% (40/130) isolates; eight were extensively drug resistant. The Beijing genotype was predominant (56.2%). The main Beijing subtypes B0/W148 and 94-32 differed in the MDR rate, 83.3% and 27.2%, respectively. The non-Beijing isolates represented five genotypes (LAM, Ural, Haarlem, X, T). The proportion of Beijing B0/W148 in the "camp" cities (originated from Gulag camps) was twice as large as in other districts of the Komi Republic. To conclude, сirculation of the MDR-associated Beijing B0/W148 cluster critically influences the current situation with MDR-TB in this Russian region. The increased prevalence of B0/W148 in the urban setting on the whole, and in the "camp cities", in particular, indirectly points to the increased transmission capacity of this successful Russian strain of M. tuberculosis.

System OMICs analysis of Mycobacterium tuberculosis Beijing B0/W148 cluster

Mycobacterium tuberculosis Beijing B0/W148 is one of the most widely distributed clusters in the Russian Federation and in some countries of the former Soviet Union. Recent studies have improved our understanding of the reasons for the "success" of the cluster but this area remains incompletely studied. Here, we focused on the system omics analysis of the RUS_B0 strain belonging to the Beijing B0/W148 cluster. Completed genome sequence of RUS_B0 (CP020093.1) and a collection of WGS for 394 cluster strains were used to describe the main genetic features of the population. In turn, proteome and transcriptome studies allowed to confirm the genomic data and to identify a number of finds that have not previously been described. Our results demonstrated that expression of the whiB6 which contains cluster-specific polymorphism (a151c) increased almost 40 times in RUS_B0. Additionally, the level of ethA transcripts in RUS_B0 was increased by more than 7 times compared to the H37Rv. Start sites for 10 genes were corrected based on the combination of proteomic and transcriptomic data. Additionally, based on the omics approach, we identified 5 new genes. In summary, our analysis allowed us to summarize the available results and also to obtain fundamentally new data.

Population structure of drug-resistant Mycobacterium tuberculosis in Central Asia

Background

Drug-resistant tuberculosis (TB) is a major public health concern threathing the success of TB control efforts, and this is particularily problematic in Central Asia. Here, we present the first analysis of the population structure of Mycobacterium tuberculosis complex isolates in the Central Asian republics Uzbekistan, Tajikistan, and Kyrgyzstan.

Methods

The study set consisted of 607 isolates with 235 from Uzbekistan, 206 from Tajikistan, and 166 from Kyrgyzstan. 24-loci MIRU-VNTR (Mycobacterial Interspersed Repetitive Units - Variable Number of Tandem Repeats) typing and spoligotyping were combined for genotyping. In addition, phenotypic drug suceptibility was performed.

Results

The population structure mainly comprises strains of the Beijing lineage (411/607). 349 of the 411 Beijing isolates formed clusters, compared to only 33 of the 196 isolates from other clades. Beijing 94–32 (n = 145) and 100–32 (n = 70) formed the largest clusters. Beijing isolates were more frequently multidrug-resistant, pre-extensively resistant (pre-XDR)- or XDR-TB than other genotypes.

Conclusions

Beijing clusters 94–32 and 100–32 are the dominant MTB genotypes in Central Asia. The relative size of 100–32 compared to previous studies in Kazakhstan and its unequal geographic distribution support the hypothesis of its more recent emergence in Central Asia. The data also demonstrate that clonal spread of resistant TB strains, particularly of the Beijing lineage, is a root of the so far uncontroled MDR-TB epidemic in Central Asia.

Cross-Border Movement of Highly Drug-Resistant Mycobacterium tuberculosis from Papua New Guinea to Australia through Torres Strait Protected Zone, 2010-2015

In this retrospective study, we used whole-genome sequencing (WGS) to delineate transmission dynamics, characterize drug-resistance markers, and identify risk factors of transmission among Papua New Guinea residents of the Torres Strait Protected Zone (TSPZ) who had tuberculosis diagnoses during 2010-2015. Of 117 isolates collected, we could acquire WGS data for 100; 79 were Beijing sublineage 2.2.1.1, which was associated with active transmission (odds ratio 6.190, 95% CI 2.221-18.077). Strains were distributed widely throughout the TSPZ. Clustering occurred more often within than between villages (p = 0.0013). Including 4 multidrug-resistant tuberculosis isolates from Australia citizens epidemiologically linked to the TSPZ into the transmission network analysis revealed 2 probable cross-border transmission events. All multidrug-resistant isolates (33/104) belonged to Beijing sublineage 2.2.1.1 and had high-level isoniazid and ethionamide co-resistance; 2 isolates were extensively drug resistant. Including WGS in regional surveillance could improve tuberculosis transmission tracking and control strategies within the TSPZ.

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.

First insights on the genetic diversity of MDR Mycobacterium tuberculosis in Lebanon

Background

Lebanon hosts a heterogeneous population coming from underdeveloped and developing countries, resulting in increasing incidences of tuberculosis over the past years. The genetic heterogeneity and lineages associated with tuberculosis, along with their resistance determinants have not been studied at the genomic level previously in the region.

Methods

Isolates were recovered from the American University of Beirut Medical Center (AUBMC). Antimicrobial susceptibility profiles were determined using the MGIT automated system for the first-line drugs at AUBMC, while second-line drug susceptibility was tested at Mayo Clinic Laboratories. Whole Genome Sequencing (WGS) was performed to classify mycobacterial lineages and highlight single nucleotide mutations causing resistance to both 1st line and 2nd line antimicrobials. wgSNP analysis provided insights on the phylogeny of the isolates along with spoligotyping and core genomic SNVs, IS6110 insertion sites, and variable number tandem repeats (VNTR).

Results

The analyzed isolates carry distinct resistance determinants to isoniazid, rifampicin, ethambutol, quinolones, and streptomycin. The isolates belonged to different lineages including the Euro/American lineage (Lineage 4) (53.8%), M. bovis (15.4%) and Delhi/Central Asia (Lineage 1) (15.4%), Beijing/East Asia (Lineage 2) (7.7%), and East Africa/Indian Ocean lineage (Lineage 3) (7.7%) showing great phylogenetic differences at the genomic level.

Conclusions

The population diversity in Lebanon holds an equally diverse and uncharacterized population of drug resistant mycobacteria. To achieve the WHO “END-TB” milestones of 2025 and 2035, Lebanon must decrease TB incidences by 95% in the next decade. This can only be done through WGS-based patient centered diagnosis with higher throughput and genomic resolution to improve treatment outcomes and to monitor transmission patterns.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3626-3) contains supplementary material, which is available to authorized users.

Line-probe assay and molecular typing reveal a potential drug resistant clone of Mycobacterium tuberculosis in Ethiopia

Background

Antimicrobial resistance is a global concern of increasing significance. Multidrug resistant tuberculosis (MDR-TB) is spreading worldwide. It is important to monitor trends of antimycobacterial resistance. This is particularly true for high TB burden countries such as Ethiopia where disproportionally less drug sensitivity data are reported from.

Methods

The prevalence of drug resistance was assessed with the line probe assay GenoType MTBDRplus in a set of 161 M. tuberculosis strains that were selected from four common lineages and sub-lineages previously identified in Ethiopia. Most of the tested M. tuberculosis isolates had been genotyped by established Spoligotyping and MIRU-VNTR typing methods.

Results

The proportion of MDR-TB among the isolates was 3.1%. Mono-resistance was 1.2% to rifampicin and 4.3% to isoniazid, and resistance to either of the two first line drugs was 8.7%. Strains of Lineage 4 had the highest resistance rate (13.6%) followed by Lineage 3 (4.9%). None of the isolates representing Lineages 1 and Lineage 7 were drug resistant. Multidrug resistance among pulmonary TB and TB lymphadenitis clinical isolates was 2.8 and 3.7%, respectively. Drug resistance of strains carrying the most prevalent spoligotype in Ethiopia - SIT149 - was further explored. Stratification by MIRU-VNTR identified one genotype with a high rate of drug resistance against Rifampicin and Isoniazid and circulation of a potential MDR-TB clone is proposed.

Conclusion

Although the strain selection was not fully randomized, the overall M. tuberculosis drug resistance rate in this strain set was 8.7% while the rate of MDR was 3.1%. In parallel, we identified a sub-lineage that showed a high rate of resistance to both rifampicin and isoniazid. These resistant strains may belong to a clone of M. tuberculosis that is circulating in the highlands of Ethiopia.

Electronic supplementary material

The online version of this article (10.1186/s40794-018-0075-3) contains supplementary material, which is available to authorized users.

Longitudinal genotyping surveillance of Mycobacterium tuberculosis in an area with high tuberculosis incidence shows high transmission rate of the modern Beijing subfamily in Japan

Tuberculosis (TB) is a severe and wide-spread infectious disease worldwide. The modern Beijing subfamily, one lineage of M. tuberculosis, reportedly has high pathogenicity and transmissibility. This study used a molecular epidemiological approach to investigate the transmissibility of the modern Beijing subfamily in the Airin area of Osaka City, Japan. During 2006-2016, we collected 596 M. tuberculosis clinical isolates in the Airin area, Osaka city, Japan. We analyzed the 24-locus variable number of tandem repeats typing optimized for the Beijing family of isolates, M. tuberculosis lineage, and patient epidemiological data. The proportion of the modern Beijing subfamily was significantly higher not only than previously obtained data for the Airin area: it was also higher than the nationwide in Japan. The rate of recent clusters, defined as a variable number of tandem repeats profile identified within two years, of the modern Beijing subfamily was significantly higher than that the rate of recent clusters of the ancient Beijing subfamily. Results suggest that TB control measures formulated with attention to the modern Beijing subfamily might be an important benchmark to understanding recent TB transmission in the area.

MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis

Background

Tuberculosis (TB) resulted in an estimated 1.7 million deaths in the year 2016. The disease is caused by the members of Mycobacterium tuberculosis complex, which includes Mycobacterium tuberculosis, Mycobacterium bovis and other closely related TB causing organisms. In order to understand the epidemiological dynamics of TB, national TB control programs often conduct standardized genotyping at 24 Mycobacterial-Interspersed-Repetitive-Units (MIRU)-Variable-Number-of-Tandem-Repeats (VNTR) loci. With the advent of next generation sequencing technology, whole-genome sequencing (WGS) has been widely used for studying TB transmission. However, an open-source software that can connect WGS and MIRU-VNTR typing is currently unavailable, which hinders interlaboratory communication. In this manuscript, we introduce the MIRU-profiler program which could be used for prediction of MIRU-VNTR profile from WGS of M. tuberculosis.

Implementation

The MIRU-profiler is implemented in shell scripting language and depends on EMBOSS software. The in-silico workflow of MIRU-profiler is similar to those described in the laboratory manuals for genotyping M. tuberculosis. Given an input genome sequence, the MIRU-profiler computes alleles at the standard 24-loci based on in-silico PCR amplicon lengths. The final output is a tab-delimited text file detailing the 24-loci MIRU-VNTR pattern of the input sequence.

Validation

The MIRU-profiler was validated on four datasets: complete genomes from NCBI-GenBank (n = 11), complete genomes for locally isolated strains sequenced using PacBio (n = 4), complete genomes for BCG vaccine strains (n = 2) and draft genomes based on 250 bp paired-end Illumina reads (n = 106).

Results

The digital MIRU-VNTR results were identical to the experimental genotyping results for complete genomes of locally isolated strains, BCG vaccine strains and five out of 11 genomes from the NCBI-GenBank. For draft genomes based on short Illumina reads, 21 out of 24 loci were inferred with a high accuracy, while a number of inaccuracies were recorded for three specific loci (ETRA, QUB11b and QUB26). One of the unique features of the MIRU-profiler was its ability to process multiple genomes in a batch. This feature was tested on all complete M. tuberculosis genome (n = 157), for which results were successfully obtained in approximately 14 min.

Conclusion

The MIRU-profiler is a rapid tool for inference of digital MIRU-VNTR profile from the assembled genome sequences. The tool can accurately infer repeat numbers at the standard 24 or 21/24 MIRU-VNTR loci from the complete or draft genomes respectively. Thus, the tool is expected to bridge the communication gap between the laboratories using WGS and those using the conventional MIRU-VNTR typing.

Insight into multidrug-resistant Beijing genotype Mycobacterium tuberculosis isolates in Myanmar

OBJECTIVES

Myanmar is a World Health Organization high tuberculosis (TB) burden country with a high multidrug-resistant (MDR)-TB burden. Of significance, a high prevalence of the Beijing genotype of Mycobacterium tuberculosis (MTB) among MDR-MTB has been reported previously. A detailed genetic characterization of TB clinical isolates was performed in order to explore whether there is an association between the prevalence of the Beijing MTB genotype and MDR-TB in Myanmar.

METHODS

A total of 265 MDR-MTB clinical isolates collected in 2010 and 2012 were subjected to spoligotyping, mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis, single nucleotide polymorphism (SNP) typing, and drug resistance-associated gene sequencing, including rpoC to detect potential compensatory evolution.

RESULTS

Of the total MDR-MTB isolates, 79.2% (210/265) were of the Beijing genotype, the majority of which were the 'modern' subtype. Beijing genotype isolates were differentiated by 15-locus MIRU-VNTR and a high clustering rate (53.0%) was observed in the modern subtype. These MIRU-VNTR patterns were similar to Beijing genotype clones spreading across Russia and Central Asia. A high prevalence of katG Ser315Thr, and genetic evidence of extensive drug resistance (XDR) and pre-XDR and compensatory mutations in rpoC were observed among clustered isolates.

CONCLUSIONS

MDR-MTB strains of the Beijing genotype might be spreading in Myanmar and present a major challenge to TB control in this country.

Mycobacterium tuberculosis Transmission among Elderly Persons, Yamagata Prefecture, Japan, 2009-2015

In many countries with low to moderate tuberculosis (TB) incidence, cases have shifted to elderly persons. It is unclear, however, whether these cases are associated with recent Mycobacterium tuberculosis transmission or represent reactivation of past disease. During 2009–2015, we performed a population-based TB investigation in Yamagata Prefecture, Japan, using in-depth contact tracing and 24-loci variable-number tandem-repeat typing optimized for Beijing family M. tuberculosis strains. We analyzed 494 strains, of which 387 (78.3%) were derived from elderly patients. Recent transmission with an epidemiologic link was confirmed in 22 clusters (70 cases). In 17 (77.3%) clusters, the source patient was elderly; 11 (64.7%) of the 17 clusters occurred in a hospital or nursing home. In this setting, the increase in TB cases was associated with M. tuberculosis transmissions from elderly persons. Prevention of transmission in places where elderly persons gather will be an effective strategy for decreasing TB incidence among predominantly elderly populations.

Major genotype families and epidemic clones of Mycobacterium tuberculosis in Omsk region, Western Siberia, Russia, marked by a high burden of tuberculosis-HIV coinfection

This population-based study characterized Mycobacterium tuberculosis isolates from HIV-positive and HIV-negative tuberculosis (TB) patients in the Omsk region in Western Siberia, Russia. We sought to gain insight into the major genotype families and epidemic and endemic clones of M. tuberculosis in the area with a high burden and adverse trend of TB/HIV coinfection. The study collection included M. tuberculosis isolates from 207 newly-diagnosed patients with pulmonary TB; 55 (26.5%) of patients were HIV-infected. The M. tuberculosis isolates were subjected to drug susceptibility testing and molecular typing based on spoligotyping and analysis of the robust genotype and cluster-specific markers. Patients with disseminated TB disease were more prevalent in the HIV-positive (34.5%) than in the HIV-negative group (4.6%) (P < .001). The Beijing genotype was predominant (62.3% of isolates), and its major subtypes were 94-32-cluster (Central Asian/Russian strain, n = 80) and B0/W148-cluster (successful Russian strain, n = 28). The main non-Beijing families were represented by Latin-American Mediterranean (14.5%), T family (11.1%), Ural (5.8%), and Haarlem (3.9%). Under multivariate logistic regression analysis, MDR was associated with Beijing genotype and not associated with HIV coinfection status (P < .001). Beijing genotype isolates were found more frequently in TB/HIV patients than in TB HIV-negative patients (74.5% versus 57.9%, respectively; P = .031). The non-Beijing genotypes were mainly drug susceptible except for the drug-resistant Ural SIT262 isolates. To summarize, the alarming situation in the Omsk region in Siberia regarding TB/HIV coinfection is seriously influenced by the active circulation of M. tuberculosis isolates of MDR-associated Beijing genotype. Among the non-Beijing families, emergence of the drug-resistant Ural family strains of spoligotype SIT262 warrants attention.

Bi-logistic model for disease dynamics caused by Mycobacterium tuberculosis in Russia

In this work, we explore epidemiological dynamics by the example of tuberculosis in Russian Federation. It has been shown that the epidemiological dynamics correlates linearly with the virulence of Mycobacterium tuberculosis during the period 1987-2012. To construct an appropriate model, we have analysed (using LogLet decomposition method) epidemiological World Health Organization (WHO) data (period 1980-2014) and obtained, as result of their integration, a curve approximated by a bi-logistic function. This fact allows a subdivision of the whole population into parts, each of them satisfies the Verhulst-like models with different constant virulences introduced into each subsystem separately. Such a subdivision could be interconnected with the heterogeneous structure of mycobacterial population that has a high ability of adaptation to the host and strong mutability.

Evolution of Mycobacterium tuberculosis: New Insights into Pathogenicity and Drug Resistance

The tuberculosis agent Mycobacterium tuberculosis has undergone a long and selective evolution toward human infection and represents one of the most widely spread pathogens due to its efficient aerosol-mediated human-to-human transmission. With the availability of more and more genome sequences, the evolutionary trajectory of this obligate pathogen becomes visible, which provides us with new insights into the molecular events governing evolution of the bacterium and its ability to accumulate drug-resistance mutations. In this review, we summarize recent developments in mycobacterial research related to this matter that are important for a better understanding of the current situation and future trends and developments in the global epidemiology of tuberculosis, as well as for possible public health intervention possibilities.

Impact of Genetic Diversity on the Biology of Mycobacterium tuberculosis Complex Strains

Tuberculosis (TB) remains the most deadly bacterial infectious disease worldwide. Its treatment and control are threatened by increasing numbers of multidrug-resistant (MDR) or nearly untreatable extensively drug-resistant (XDR) strains. New concepts are therefore urgently needed to understand the factors driving the TB epidemics and the spread of different strain populations, especially in association with drug resistance. Classical genotyping and, more recently, whole-genome sequencing (WGS) revealed that the world population of tubercle bacilli is more diverse than previously thought. Several major phylogenetic lineages can be distinguished, which are associated with their sympatric host population. Distinct clonal (sub)populations can even coexist within infected patients. WGS is now used as the ultimate approach for differentiating clinical isolates and for linking phenotypic to genomic variation from lineage to strain levels. Multiple lines of evidence indicate that the genetic diversity of TB strains translates into pathobiological consequences, and key molecular mechanisms probably involved in differential pathoadaptation of some main lineages have recently been identified. Evidence also accumulates on molecular mechanisms putatively fostering the emergence and rapid expansion of particular MDR and XDR strain groups in some world regions. However, further integrative studies will be needed for complete elucidation of the mechanisms that allow the pathogen to infect its host, acquire multidrug resistance, and transmit so efficiently. Such knowledge will be key for the development of the most effective new diagnostics, drugs, and vaccination strategies.

Pyrazinamide resistance-conferring mutations in pncA and the transmission of multidrug resistant TB in Georgia

Background

The ongoing epidemic of multidrug-resistant tuberculosis (MDR-TB) in Georgia highlights the need for more effective control strategies. A new regimen to treat MDR-TB that includes pyrazinamide (PZA) is currently being evaluated and PZA resistance status will largely influence the success of current and future treatment strategies. PZA susceptibility testing was not routinely performed at the National Reference Laboratory (NRL) in Tbilisi between 2010 and September 2015. We here provide a first insight into the prevalence of PZA resistant TB in this region.

Methods

Phenotypic susceptibility to PZA was determined in a convenience collection of well-characterised TB patient isolates collected at the NRL in Tbilisi between 2012 and 2013. In addition, the pncA gene was sequenced and whole genome sequencing was performed on two isolates.

Results

Out of 57 isolates tested 33 (57.9%) showed phenotypic drug resistance to PZA and had a single pncA mutation. All of these 33 isolates were MDR-TB strains. pncA mutations were absent in all but one of the 24 PZA susceptible isolate. In total we found 18 polymorphisms in the pncA gene. From the two major MDR-TB clusters represented (94–32 and 100–32), 10 of 15, 67.0% and 13 of 14, 93.0% strains, respectively were PZA resistant. We also identified a member of the potentially highly transmissive clade A strain carrying the characteristic I6L substitution in PncA. Another strain with the same MLVA type as the clade A strain acquired a different mutation in pncA and was genetically more distantly related suggesting that different branches of this particular lineage have been introduced into this region.

Conclusion

In this high MDR-TB setting more than half of the tested MDR-TB isolates were resistant to PZA. As PZA is part of current and planned MDR-TB treatment regimens this is alarming and deserves the attention of health authorities. Based on our typing and sequence analysis results we conclude that PZA resistance is the result of primary transmission as well as acquisition within the patient and recommend prospective genotyping and PZA resistance testing in high MDR-TB settings. This is of utmost importance in order to preserve bacterial susceptibility to PZA to help protect (new) second line drugs in PZA containing regimens.

A Nonsynonymous SNP Catalog of Mycobacterium tuberculosis Virulence Genes and Its Use for Detecting New Potentially Virulent Sublineages

Mycobacterium tuberculosis is divided into several distinct lineages, and various genetic markers such as IS-elements, VNTR, and SNPs are used for lineage identification. We propose an M. tuberculosis classification approach based on functional polymorphisms in virulence genes. An M. tuberculosis virulence genes catalog has been established, including 319 genes from various protein groups, such as proteases, cell wall proteins, fatty acid and lipid metabolism proteins, sigma factors, toxin–antitoxin systems. Another catalog of 1,573 M. tuberculosis isolates of different lineages has been developed. The developed SNP-calling program has identified 3,563 nonsynonymous SNPs. The constructed SNP-based phylogeny reflected the evolutionary relationship between lineages and detected new sublineages. SNP analysis of sublineage F15/LAM4/KZN revealed four lineage-specific mutations in cyp125, mce3B, vapC25, and vapB34. The Ural lineage has been divided into two geographical clusters based on different SNPs in virulence genes. A new sublineage, B0/N-90, was detected inside the Beijing-B0/W-148 by SNPs in irtB, mce3F and vapC46. We have found 27 members of B0/N-90 among the 227 available genomes of the Beijing-B0/W-148 sublineage. Whole-genome sequencing of strain B9741, isolated from an HIV-positive patient, was demonstrated to belong to the new B0/N-90 group. A primer set for PCR detection of B0/N-90 lineage-specific mutations has been developed. The prospective use of mce3 mutant genes as genetically engineered vaccine is discussed.

Positive epistasis of major low-cost drug resistance mutations rpoB531-TTG and katG315-ACC depends on the phylogenetic background of Mycobacterium tuberculosis strains

Mycobacterium tuberculosis Beijing genotype strains increasingly circulate in different world regions, either as historical endemic, e.g. in East Asia, or recently imported, e.g. in South America, and this family is regarded as the most successful lineage of the global tuberculosis (TB) epidemic. Here we analysed the transmission capacity of these strains in the context of their phylogenetic background and drug resistance mutations. The study collection included all multidrug resistant (MDR) strains of Beijing genotype isolated in Beijing Chest Hospital, the largest tertiary TB facility in North China, in 2011-2013 (n = 278). Strains were subjected to NTF/IS6110 and 24-loci MIRU-VNTR analysis. Drug resistance mutations were detected in rpoB, katG, inhA and oxyR-ahpC. A total of 58 and 220 strains were assigned to the ancient and modern Beijing sublineages, respectively. 24-MIRU-VNTR clustering was higher in modern versus ancient Beijing strains (35.9% vs. 12.1%; P <0.001). After taking into consideration the presence of rpoB and katG mutations, clustering decreased to 15.9% in modern and 0% in ancient strains. The most frequent combination of mutations (rpoB531-TTG and katG315-ACC) was more prevalent in clustered versus non-clustered isolates in the modern sublineage (23/35 vs. 47/185; P <0.0001). To conclude, a combination of the known low-fitness-cost rpoB531-TTG and katG315-ACC mutations likely facilitates the increased transmission ability of MDR strains of the modern but not ancient Beijing sublineage. Accordingly, positive epistasis of major low-cost drug resistance-conferring mutations is influenced by the phylogenetic background of M. tuberculosis strains.

Strain-specific estimation of epidemic success provides insights into the transmission dynamics of tuberculosis

The transmission dynamics of tuberculosis involves complex interactions of socio-economic and, possibly, microbiological factors. We describe an analytical framework to infer factors of epidemic success based on the joint analysis of epidemiological, clinical and pathogen genetic data. We derive isolate-specific, genetic distance-based estimates of epidemic success, and we represent success-related time-dependent concepts, namely epidemicity and endemicity, by restricting analysis to specific time scales. The method is applied to analyze a surveillance-based cohort of 1,641 tuberculosis patients with minisatellite-based isolate genotypes. Known predictors of isolate endemicity (older age, native status) and epidemicity (younger age, sputum smear positivity) were identified with high confidence (P < 0.001). Long-term epidemic success also correlated with the ability of Euro-American and Beijing MTBC lineages to cause active pulmonary infection, independent of patient age and country of origin. Our results demonstrate how important insights into the transmission dynamics of tuberculosis can be gained from active surveillance data.

Association between genotype and drug resistance profiles of Mycobacterium tuberculosis strains circulating in China in a national drug resistance survey

We describe the population structure of a representative collection of 3,133 Mycobacterium tuberculosis isolates, collected within the framework of a national resistance survey from 2007 in China. Genotyping data indicate that the epidemic strains in China can be divided into seven major complexes, of which 92% belonged to the East Asian (mainly Beijing strains) or the Euro-American lineage. The epidemic Beijing strains in China are closely related to the Beijing B0/W148 strain earlier described in Russia and a large cluster of these strains has spread national wide. The density of Beijing strains is high in the whole of China (average 70%), but the highest prevalence was found North of the Yellow river. The Euro-American lineage consists of three sublineages (sublineage_1, 2, and 3) and is more prevalent in the South. Beijing lineage showed the highest cluster rate of 48% and a significantly higher level of resistance to rifampicin (14%, p<0.001), ethambutol (9%, p = 0.001), and ofloxacin (5%, p = 0.011). Within the Euro-American Lineage, sublineage_3 revealed the highest cluster rate (28%) and presented a significantly elevated level of resistance to streptomycin (44%, p<0.001). Our findings suggest that standardised treatment in this region may have contributed to the successful spread of certain strains: sublineage_3 in the Euro-American lineage may have thrived when streptomycin was used without rifampicin for treatment, while later under DOTS based treatment, in which rifampicin plays a key role, Beijing lineage appears to be spreading.

The Evolution of Strain Typing in the Mycobacterium tuberculosis Complex

Tuberculosis (TB) is a contagious disease with a complex epidemiology. Therefore, molecular typing (genotyping) of Mycobacterium tuberculosis complex (MTBC) strains is of primary importance to effectively guide outbreak investigations, define transmission dynamics and assist global epidemiological surveillance of the disease. Large-scale genotyping is also needed to get better insights into the biological diversity and the evolution of the pathogen. Thanks to its shorter turnaround and simple numerical nomenclature system, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) typing, based on 24 standardized plus 4 hypervariable loci, optionally combined with spoligotyping, has replaced IS6110 DNA fingerprinting over the last decade as a gold standard among classical strain typing methods for many applications. With the continuous progress and decreasing costs of next-generation sequencing (NGS) technologies, typing based on whole genome sequencing (WGS) is now increasingly performed for near complete exploitation of the available genetic information. However, some important challenges remain such as the lack of standardization of WGS analysis pipelines, the need of databases for sharing WGS data at a global level, and a better understanding of the relevant genomic distances for defining clusters of recent TB transmission in different epidemiological contexts. This chapter provides an overview of the evolution of genotyping methods over the last three decades, which culminated with the development of WGS-based methods. It addresses the relative advantages and limitations of these techniques, indicates current challenges and potential directions for facilitating standardization of WGS-based typing, and provides suggestions on what method to use depending on the specific research question.

Evolutionary History and Ongoing Transmission of Phylogenetic Sublineages of Mycobacterium tuberculosis Beijing Genotype in China

Mycobacterium tuberculosis Beijing genotype originated in China and has undergone a dramatic population growth and global spread in the last century. Here, a collection of M. tuberculosis Beijing family isolates from different provinces across all China was genotyped by high-resolution (24-MIRU-VNTR) and low-resolution, high-rank (modern and ancient sublineages) markers. The molecular profiles and global and local phylogenies were compared to the strain phenotype and patient data. The phylogeographic patterns observed in the studied collection demonstrate that large-scale (but not middle/small-scale) distance remains one of the decisive factors of the genetic divergence of M. tuberculosis populations. Analysis of diversity and network topology of the local collections appears to corroborate a recent intriguing hypothesis about Beijing genotype originating in South China. Placing our results within the Eurasian context suggested that important Russian B0/W148 and Asian/Russian A0/94-32 epidemic clones of the Beijing genotype could trace their origins to the northeastern and northwestern regions of China, respectively. The higher clustering of the modern isolates in children and lack of increased MDR rate in any sublineage suggest that not association with drug resistance but other (e.g., speculatively, virulence-related) properties underlie an enhanced dissemination of the evolutionarily recent, modern sublineage of the Beijing genotype in China.

Urgent Implementation in a Hospital Setting of a Strategy To Rule Out Secondary Cases Caused by Imported Extensively Drug-Resistant Mycobacterium tuberculosis Strains at Diagnosis

Current migratory movements require new strategies for rapidly tracking the transmission of high-risk imported Mycobacterium tuberculosis strains. Whole-genome sequencing (WGS) enables us to identify single-nucleotide polymorphisms (SNPs) and therefore design PCRs to track specific relevant strains. However, fast implementation of these strategies in the hospital setting is difficult because professionals working in diagnostics, molecular epidemiology, and genomics are generally at separate institutions. In this study, we describe the urgent implementation of a system that integrates genomics and molecular tools in a genuine high-risk epidemiological alert involving 2 independent importations of extensively drug resistant (XDR) and pre-XDR Beijing M. tuberculosis strains from Russia into Spain. Both cases involved commercial sex workers with long-standing tuberculosis (TB). The system was based on strain-specific PCRs tailored from WGS data that were transferred to the local node that was managing the epidemiological alert. The optimized tests were available for prospective implementation in the local node 33 working days after receiving the primary cultures of the XDR strains and were applied to all 42 new incident cases. An interpretable result was obtained in each case (directly from sputum for 27 stain-positive cases) and corresponded to the amplification profiles for strains other than the targeted pre-XDR and XDR strains, which made it possible to prospectively rule out transmission of these high-risk strains at diagnosis.