Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis

The first insight into Mycobacterium tuberculosis complex isolates in the lower northern region in Thailand

BACKGROUND

Tuberculosis (TB) remains an important infectious disease and different genotypes have been reported. This study aimed to investigate the genetic diversity and molecular epidemiology of TB in the lower northern region of Thailand, where genotyping data are limited.

METHODS

A total of 159 Mycobacterium tuberculosis complex (MTBC) isolates from this region were genotyped by spoligotyping and the major spoligotypes were further subdivided by the mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) method.

RESULTS

Spoligotyping identified 34 types and classified them into 14 clusters. East African-Indian (EAI) groups were the most frequent (44.7%), followed by Beijing (36.5%), with a higher prevalence of drug resistance. By 15-loci MIRU-VNTR typing, the major groups of the Beijing and EAI2_NTB were further differentiated into 44 and 21 subtypes forming 9 and 5 subclusters with cluster rates of 0.26 and 0.44, respectively. The Hunter-Gaston Discriminatory Index among the Beijing and EAI2_NTB groups were 0.987 and 0.931, respectively, indicating high diversity.

CONCLUSIONS

This is the first look at the MTBC genotypes in the lower northern region of Thailand, which could aid in understanding the distribution and potential spread of MTBC and Mycobacterium bovis in the target region to support TB control in Thailand.

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.

High-throughput phenogenotyping of Mycobacteria tuberculosis clinical strains reveals bacterial determinants of treatment outcomes

Background

Combatting the tuberculosis (TB) epidemic caused by Mycobacterium tuberculosis ( Mtb ) necessitates a better understanding of the factors contributing to patient clinical outcomes and transmission. While host and environmental factors have been evaluated, the impact of Mtb genetic background and phenotypic diversity is underexplored. Previous work has made associations between Mtb genetic lineages and some clinical and epidemiological features, but the bacterial traits underlying these connections are largely unknown.

Methods

We developed a high-throughput functional genomics platform for defining genotype-phenotype relationships across a panel of Mtb clinical isolates. These phenotypic fitness profiles function as intermediate traits which can be linked to Mtb genetic variants and associated with clinical and epidemiological outcomes. We applied this approach to a collection of 158 Mtb strains from a study of Mtb transmission in Ho Chi Minh City, Vietnam. Mtb strains were genetically tagged in multiplicate, which allowed us to pool the strains and assess in vitro competitive fitness using deep sequencing across a set of 14 host-relevant antibiotic and metabolic conditions. Phylogenetic and monogenic associations with these intermediate traits were identified and then associated with clinical outcomes.

Findings

Mtb clinical strains have a broad range of growth and drug response dynamics that can be clustered by their phylogenetic relationships. We identified novel monogenic associations with Mtb fitness in various metabolic and antibiotic conditions. Among these, we find that mutations in Rv1339 , a phosphodiesterase, which were identified through their association with slow growth in glycerol, are further associated with treatment failure. We also identify a previously uncharacterized subclade of Lineage 1 strains (L1.1.1.1) that is phenotypically distinguished by slow growth under most antibiotic and metabolic stress conditions in vitro . This clade is associated with cavitary disease, treatment failure, and demonstrates increased transmission potential.

Interpretation

High-throughput phenogenotyping of Mtb clinical strains enabled bacterial intermediate trait identification that can provide a mechanistic link between Mtb genetic variation and patient clinical outcomes. Mtb strains associated with cavitary disease, treatment failure, and transmission potential display intermediate phenotypes distinguished by slow growth under various antibiotic and metabolic conditions. These data suggest that Mtb growth regulation is an adaptive advantage for host bacterial success in human populations, in at least some circumstances. These data further suggest markers for the underlying bacterial processes that govern these clinical outcomes.

Funding

National Institutes of Allergy and Infectious Diseases: P01 AI132130 (SS, SMF); P01 AI143575 (XW, SMF); U19 AI142793 (QL, SMF); 5T32AI132120-03 (SS); 5T32AI132120-04 (SS); 5T32AI049928-17 (SS) Wellcome Trust Fellowship in Public Health and Tropical Medicine: 097124/Z/11/Z (NTTT) National Health and Medical Research Council (NHMRC)/A*STAR joint call: APP1056689 (SJD) The funding sources had no involvement in study methodology, data collection, analysis, and interpretation nor in the writing or submission of the manuscript.

Research in context

Evidence before this study: We used different combinations of the words mycobacterium tuberculosis, tuberculosis, clinical strains, intermediate phenotypes, genetic barcoding, phenogenomics, cavitary disease, treatment failure, and transmission to search the PubMed database for all studies published up until January 20 th , 2022. We only considered English language publications, which biases our search. Previous work linking Mtb determinants to clinical or epidemiological data has made associations between bacterial lineage, or less frequently, genetic polymorphisms to in vitro or in vivo models of pathogenesis, transmission, and clinical outcomes such as cavitary disease, treatment failure, delayed culture conversion, and severity. Many of these studies focus on the global pandemic Lineage 2 and Lineage 4 Mtb strains due in part to a deletion in a polyketide synthase implicated in host-pathogen interactions. There are a number of Mtb GWAS studies that have led to novel genetic determinants of in vitro drug resistance and tolerance. Previous Mtb GWAS analyses with clinical outcomes did not experimentally test any predicted phenotypes of the clinical strains. Published laboratory-based studies of Mtb clinical strains involve relatively small numbers of strains, do not identify the genetic basis of relevant phenotypes, or link findings to the corresponding clinical outcomes. There are two recent studies of other pathogens that describe phenogenomic analyses. One study of 331 M. abscessus clinical strains performed one-by-one phenotyping to identify bacterial features associated with clearance of infection and another details a competition experiment utilizing three barcoded Plasmodium falciparum clinical isolates to assay antimalarial fitness and resistance. Added value of this study: We developed a functional genomics platform to perform high-throughput phenotyping of Mtb clinical strains. We then used these phenotypes as intermediate traits to identify novel bacterial genetic features associated with clinical outcomes. We leveraged this platform with a sample of 158 Mtb clinical strains from a cross sectional study of Mtb transmission in Ho Chi Minh City, Vietnam. To enable high-throughput phenotyping of large numbers of Mtb clinical isolates, we applied a DNA barcoding approach that has not been previously utilized for the high-throughput analysis of Mtb clinical strains. This approach allowed us to perform pooled competitive fitness assays, tracking strain fitness using deep sequencing. We measured the replicative fitness of the clinical strains in multiplicate under 14 metabolic and antibiotic stress condition. To our knowledge, this is the largest phenotypic screen of Mtb clinical isolates to date. We performed bacterial GWAS to delineate the Mtb genetic variants associated with each fitness phenotype, identifying monogenic associations with several conditions. We then defined Mtb phenotypic and genetic features associated with clinical outcomes. We find that a subclade of Mtb strains, defined by variants largely involved in fatty acid metabolic pathways, share a universal slow growth phenotype that is associated with cavitary disease, treatment failure and increased transmission potential in Vietnam. We also find that mutations in Rv1339 , a poorly characterized phosphodiesterase, also associate with slow growth in vitro and with treatment failure in patients. Implications of all the available evidence: Phenogenomic profiling demonstrates that Mtb strains exhibit distinct growth characteristics under metabolic and antibiotic stress conditions. These fitness profiles can serve as intermediate traits for GWAS and association with clinical outcomes. Intermediate phenotyping allows us to examine potential processes by which bacterial strain differences contribute to clinical outcomes. Our study identifies clinical strains with slow growth phenotypes under in vitro models of antibiotic and host-like metabolic conditions that are associated with adverse clinical outcomes. It is possible that the bacterial intermediate phenotypes we identified are directly related to the mechanisms of these outcomes, or they may serve as markers for the causal yet unidentified bacterial determinants. Via the intermediate phenotyping, we also discovered a surprising diversity in Mtb responses to the new anti-mycobacterial drugs that target central metabolic processes, which will be important in considering roll-out of these new agents. Our study and others that have identified Mtb determinants of TB clinical and epidemiological phenotypes should inform efforts to improve diagnostics and drug regimen design.

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.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

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.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

Molecular epidemiology and drug susceptibility profiles of Mycobacterium tuberculosis complex isolates from Northern Ghana

OBJECTIVE

We conducted a cross-sectional study in the five administrative regions of Northern Ghana to determine the diversity of Mycobacterium tuberculosis complex (MTBC) sub/lineages and their susceptibility to isoniazid (INH) and rifampicin (RIF).

METHODS

Sputum specimens were collected and cultured from 566 pulmonary tuberculosis patients reporting to 17 health facilities from 2015 to 2019. Mycobacterial isolates obtained from solid cultures were confirmed as members of the MTBC by PCR amplification of IS6110 and rpoß and assigned lineages and sub-lineages using spoligotyping.

RESULTS

Of 294 mycobacterial isolates recovered, MTBC species identified were: M. tuberculosis sensu stricto (Mtbss) 241 (82.0%), M. africanum 41 (13.9%) and M. bovis four (1.4%) with eight (2.7%) unidentified. The human-adapted lineages (L) identified (N=279) were L1 (8/279, 2.9%), L2 (15/279, 5.4%), L3 (7/279, 2.5%), L4 (208/279, 74.5%), L5 (13/279, 4.7%) and L6 (28/279, 10.0%) with three unidentified lineages. Among the 208 L4, the dominant sub-lineages in the region were the Cameroon 120/208 (57.7%) and Ghana 50/208 (24.0%). We found 4.4% (13/294) and 0.7% (2/294) of the patients infected with MTBC isolates resistant to INH only and RIF only, respectively, with 2.4% (7/294) being infected with MDR strains. Whereas L6 was associated with the elderly, we identified that the Ghana sub-lineage of L4 was associated with both INH and MDR (p<0.05), making them important TB pathogens in Northern Ghana and a growing public health concern.

Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Genotyping and Molecular Characterization of Fluoroquinolone's Resistance Among Multidrug-Resistant Mycobacterium tuberculosis in Southwest of China

Although fluoroquinolones (FQs) are the backbone drugs for the treatment of multidrug-resistant tuberculosis (MDR-TB), the knowledge about the resistance pattern and molecular characterization of new-generation FQs in Chongqing is limited. This study aimed to investigate the resistance rate and mutation types of later-generation FQs against MDR-TB in Chongqing, and further to explore the relationship between different genotypes and phenotypes. A total of 967 clinical strains were characterized using multilocus sequence typing and drug susceptibility testing, followed by analysis of genotype/phenotype association. The 229 (23.7%, 229/967) isolates were identified as MDR-TB. The most effective agent against MDR-TB was gatifloxacin (GFX) (20.1%, 46/229), and the highest resistant rate was observed in ofloxacin (OFX) (41.0%, 94/229). Of the 190 strains (83.0%) identified as Beijing genotype, 111 isolates were modern Beijing genotype (58.4%) and 79 isolates were ancient Beijing genotype (41.6%). By analyzing 94 OFX-resistant isolates, 13 isolates were clustered with the cumulative clustering rate of 13.8% (13/94). Of the 91 isolates (39.7%, 91/229) with a mutation in gyrA gene, mutation in codon 94 was the most prevalent. Only 15 isolates (6.6%, 15/229) harbored a mutation in gyrB gene. There was no significant difference in the mutation rate of gyrA gene between Beijing and non-Beijing genotype, clustered isolates, and nonclustered isolates (p > 0.05).

Does the BCG vaccine have different effects on strains of tuberculosis?

Several explanations have been suggested concerning the variety in bacille Calmette-Guerin (BCG) vaccine efficacy on strains of Mycobacterium tuberculosis (Mtb). This study aimed to compare the effect of BCG vaccination history in the prevention of the occurrence of Mtb-Beijing and non-Beijing strains. In this cross-sectional study, 64 patients with pulmonary tuberculosis (TB) were recruited from the Iranian border provinces (North West and West). Isolates were subjected to restriction fragment length polymorphism (RFLP) analysis, using the insertion sequence IS6110 as a probe (IS6110 RFLP) and drug susceptibility testing using the proportion method. Samples were analyzed with Gel Compare II 6.6 and spss version 18. The mean age [standard deviation (SD)] of the patients was 54·4 (SD = 17·0). Overall, 49 cases (76·56%) had no BCG vaccination scar. The prevalence of Beijing strains was 9·38% and drug resistance proportion among the isolates was 14·1% (nine cases). There was a significant relationship between Beijing strains and tuberculosis (TB)-drug resistance in isolates (χ²  = 26·29, P < 0·001). There was also a strong association between vaccination history and Beijing strains (χ²  = 13·23, P = 0·002). Also, a statistical relationship was observed between Beijing strains and drug-resistant TB among patients with a history of vaccination (χ²  = 7·47, P = 0·002). This association was not maintained in the unvaccinated group (P = 0·102). These findings confirm the claim that the vaccine has different effects on different subspecies of tuberculosis. The cause of the high probability of drug resistance in patients with Beijing-TB and vaccination history requires further investigation with a higher sample size.

Genetic Variation Putatively Associated with Mycobacterium tuberculosis Resistance to Perchlozone, a New Thiosemicarbazone: Clues from Whole Genome Sequencing and Implications for Treatment of Multidrug-Resistant Tuberculosis

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.

Two Novel katG Mutations Conferring Isoniazid Resistance in Mycobacterium tuberculosis

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, is among the top 10 leading causes of death worldwide. The treatment course for TB is challenging; it requires antibiotic administration for at least 6 months, and bacterial drug resistance makes treatment even more difficult. Understanding the mechanisms of resistance is important for improving treatment. To investigate new mechanisms of isoniazid (INH) resistance, we obtained three INH-resistant (INH-R) M. tuberculosis clinical isolates collected by the Taiwan Centers for Disease Control (TCDC) and sequenced genes known to harbor INH resistance-conferring mutations. Then, the relationship between the mutations and INH resistance of these three INH-R isolates was investigated. Sequencing of the INH-R isolates identified three novel katG mutations resulting in R146P, W341R, and L398P KatG proteins, respectively. To investigate the correlation between the observed INH-R phenotypes of the clinical isolates and these katG mutations, wild-type katG from H37Rv was expressed on a plasmid (pMN437-katG) in the isolates, and their susceptibilities to INH were determined. The plasmid expressing H37Rv katG restored INH susceptibility in the two INH-R isolates encoding the W341R KatG and L398P KatG proteins. In contrast, no phenotypic change was observed in the KatG R146P isolate harboring pMN437-katG. H37Rv isogenic mutant with W341R KatG or L398P KatG was further generated. Both showed resistant to INH. In conclusion, W341R KatG and L398P KatG conferred resistance to INH in M. tuberculosis, whereas R146P KatG did not affect the INH susceptibility of M. tuberculosis.

Mycobacterial Virulence Factors: Surface-Exposed Lipids and Secreted Proteins

The clinically important Mycobacterium tuberculosis (M. tb) and related mycobacterial pathogens use various virulence mechanisms to survive and cause disease in their hosts. Several well-established virulence factors include the surface-exposed lipids in the mycobacterial outer membrane, as well as the Esx family proteins and the Pro-Glu (PE)/ Pro-Pro-Glu (PPE) family proteins secreted by type VII secretion systems (T7SS). Five ESX T7SS exist in M. tb and three—EsxA secretion system-1 (ESX-1), ESX-3, and ESX-5—have been implicated in virulence, yet only the structures of ESX-3 and ESX-5 have been solved to date. Here, we summarize the current research on three outer membrane lipids—phthiocerol dimycocerosates, phenolic glycolipids, and sulfolipids—as well as the secretion machinery and substrates of three mycobacterial T7SS—ESX-1, ESX-3, and ESX-5. We propose a structural model of the M. tb ESX-1 system based on the latest structural findings of the ESX-3 and ESX-5 secretion apparatuses to gain insight into the transport mechanism of ESX-associated virulence factors.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Lineage-specific differences in lipid metabolism and its impact on clinical strains of Mycobacterium tuberculosis

Mycobacterium tuberculosis (M. tb) is the causative agent of TB and its incidences has been on the rise since 1993. Lipid metabolism is an imperative metabolic process, which grants M. tb the ability to utilize host-derived lipids as a secondary source of nutrition during infection. In addition to degrading host lipids, M. tb is proficient at using lipids, such as cholesterol, to facilitate its entry into macrophages. Mycolic acids, constituents of the mycobacterial cell wall, offer protection and aid in persistence of the bacterium. These are effectively synthesized using a complex fatty acid synthase system. Many pathogenesis studies have reported differences in lipid-metabolism of clinical strains of M. tb that belongs to diverse lineages of the Mycobacterium tuberculosis complex (MTBC). East-Asian and Euro-American lineages possess "unique" cell wall-associated lipids compared to the less transmissible Ethiopian lineage, which may offer these lineages a competitive advantage. Therefore, it is crucial to comprehend the complexities among the MTBC lineages with lipid metabolism and their impact on virulence, transmissibility and pathogenesis. Thus, this review provides an insight into lipid metabolism in various lineages of the MTBC and their impact on virulence and persistence during infection, as this may provide critical insight into developing novel therapeutics to combat TB.

Natural and trained innate immunity against Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb) infection, remains a major global health emergency. It is estimated that one third of global population are affected, predominantly with latent granuloma form of the disease. Mtb co-evolved with humans, for its obligatory intra-macrophage phagosome habitat and slow replication, balanced against unique mycobacterial innate immunity, which appears to be highly complex. TB is transmitted via cough aerosol Mtb inhalation. Bovine TB attenuated Bacillus Calmette Guerin (BCG) live vaccine has been in practice for protection of young children from severe disseminated Mtb infection, but not sufficiently for their lungs, as obtained by trials in TB endemic community. To augment BCG vaccine-driven innate and adaptive immunity for neonates and better protection against adult pulmonary TB, a number of BCG pre-vaccination based, subset vaccine candidates have been tested via animal preclinical, followed by safe clinical trials. BCG also enhances innate macrophage trained immunity and memory, through primordial intracellular Toll-like receptors (TLRs) 7 and 9, which recognise distinct mycobacterial molecular pattern signature. This signature is transmitted by TLR signalling via nuclear factor-κB, for activating innate immune transcription and expression of gene profiling in a mycobacterial signature-specific manner. These are epigenetically imprinted in reprogramming of distinct chromatin areas for innate immune memory, to be recalled following lung reinfection. Unique TB innate immunity and its trained memory are considered independent from adaptive immune B and T cells. On the other hand, adaptive immunity is crucial in Mtb containment in granulomatous latency, supported by innate immune cell infiltration. In nearly 5-10 % of susceptible people, latent TB may be activated due to immune evasion by Mtb from intracellular phagosome within macrophage, perpetrating TB. However, BCG and new recombinant BCG vaccines have the capacity, as indicated in pre- and clinical trials, to overcome such Mtb evasion. Various strategies include pro-inflammatory-bactericidal type 1 polarisation (M1) phenotype of the infected macrophage, involving thrombospondin-TLR pathway. Saprophytic M. smegmatis-based recombinant vaccines are also promising candidates against TB. BCG vaccination of neonates/infants in TB endemic countries also reduced their pneumonia caused by various microbes independent of TB immunity. Here, we discuss host immune response against Mtb, its immune evasion strategies, and the important role innate immunity plays in the development of protection against TB.

A landscape of genomic alterations at the root of a near-untreatable tuberculosis epidemic

BACKGROUND

Atypical Beijing genotype Mycobacterium tuberculosis strains are widespread in South Africa and have acquired resistance to up to 13 drugs on multiple occasions. It is puzzling that these strains have retained fitness and transmissibility despite the potential fitness cost associated with drug resistance mutations.

METHODS

We conducted Illumina sequencing of 211 Beijing genotype M. tuberculosis isolates to facilitate the detection of genomic features that may promote acquisition of drug resistance and restore fitness in highly resistant atypical Beijing forms. Phylogenetic and comparative genomic analysis was done to determine changes that are unique to the resistant strains that also transmit well. Minimum inhibitory concentration (MIC) determination for streptomycin and bedaquiline was done for a limited number of isolates to demonstrate a difference in MIC between isolates with and without certain variants.

RESULTS

Phylogenetic analysis confirmed that two clades of atypical Beijing strains have independently developed resistance to virtually all the potent drugs included in standard (pre-bedaquiline) drug-resistant TB treatment regimens. We show that undetected drug resistance in a progenitor strain was likely instrumental in this resistance acquisition. In this cohort, ethionamide (ethA A381P) resistance would be missed in first-line drug-susceptible isolates, and streptomycin (gidB L79S) resistance may be missed due to an MIC close to the critical concentration. Subsequent inadequate treatment historically led to amplification of resistance and facilitated spread of the strains. Bedaquiline resistance was found in a small number of isolates, despite lack of exposure to the drug. The highly resistant clades also carry inhA promoter mutations, which arose after ethA and katG mutations. In these isolates, inhA promoter mutations do not alter drug resistance, suggesting a possible alternative role.

CONCLUSION

The presence of the ethA mutation in otherwise susceptible isolates from ethionamide-naïve patients demonstrates that known exposure is not an adequate indicator of drug susceptibility. Similarly, it is demonstrated that bedaquiline resistance can occur without exposure to the drug. Inappropriate treatment regimens, due to missed resistance, leads to amplification of resistance, and transmission. We put these results into the context of current WHO treatment regimens, underscoring the risks of treatment without knowledge of the full drug resistance profile.

Circulation of M. tuberculosis Beijing genotype in Latin America and the Caribbean

Lineage 2 (East Asian), which includes the Beijing genotype, is one of the most prevalent lineages of Mycobacterium tuberculosis (Mtb) throughout the world. The Beijing family is associated to hypervirulence and drug-resistant tuberculosis. The study of this genotype's circulation in Latin America is crucial for achieving total control of TB, the goal established by the World Health Organization, for the American sub-continent, before 2035. In this sense, the present work presents an overview of the status of the Beijing genotype for this region, with a bibliographical review, and data analysis of MIRU-VNTRs for available Beijing isolates. Certain countries present a prevalent trend of <5%, suggesting low transmissibility for the region, with the exception of Cuba (17.2%), Perú (16%) and Colombia (5%). Minimum Spanning Tree analysis, obtained from MIRU-VNTR data, shows distribution of specific clonal complex strains in each country. From this data, in most countries, we found that molecular epidemiology has not been a tool used for the control of TB, suggesting that the Beijing genotype may be underestimated in Latin America. It is recommended that countries with the highest incidence of the Beijing genotype use effective control strategies and increased care, as a requirement for public health systems.

Dispersal of Mycobacterium tuberculosis Driven by Historical European Trade in the South Pacific

Mycobacterium tuberculosis (Mtb) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background

The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining.

Results

We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination.

Conclusions

The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

Reference set of Mycobacterium tuberculosis clinical strains: A tool for research and product development

The Mycobacterium tuberculosis complex (MTBC) causes tuberculosis (TB) in humans and various other mammals. The human-adapted members of the MTBC comprise seven phylogenetic lineages that differ in their geographical distribution. There is growing evidence that this phylogeographic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC laboratory strains H37Rv and Erdman, both of which belong to Lineage 4. Relying on only a few laboratory-adapted strains can be misleading as study results might not be directly transferrable to clinical settings where patients are infected with a diverse array of strains, including drug-resistant variants. Here, we argue for the need to expand TB research and development by incorporating the phylogenetic diversity of the MTBC. To facilitate such work, we have assembled a group of 20 genetically well-characterized clinical strains representing the seven known human-adapted MTBC lineages. With the "MTBC clinical strains reference set" we aim to provide a standardized resource for the TB community. We hope it will enable more direct comparisons between studies that explore the physiology of MTBC beyond the laboratory strains used thus far. We anticipate that detailed phenotypic analyses of this reference strain set will increase our understanding of TB biology and assist in the development of new control tools that are broadly effective.

A deletion in the RD105 region confers resistance to multiple drugs in Mycobacterium tuberculosis

BACKGROUND

The emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), especially those that are multidrug resistant poses a serious threat to global tuberculosis control. However, the mechanism underlying the occurrence of drug resistance against more than one drug is poorly understood. Given that the Beijing/W strains are associated with outbreaks and multidrug resistance, they may harbor a genetic advantage and provide useful insight into the disease. One marker found in all Beijing/W Mtb strains is a deletion of RD105 region that results in a gene fusion, Rv0071/74, with a variable number (3-9 m) of VDP (V: Val, D: Asp; P: Pro) repeats (coded by gtggacccg repeat sequences) at the N-terminal. Here, we report that this variable number of VDP repeats in Rv0071/74 regulates the development of multidrug resistance.

RESULTS

We collected and analyzed 1255 Beijing/W clinical strains. The results showed that the number of VDP repeats in Rv0071/74 was related to the development of multidrug resistance, and the deletion of Rv0071/74-9 m from Beijing/W clinical strain restored drug susceptibility. Rv0071/74-9 m also increased resistance to multiple drugs when transferred to different mycobacterial strains. Cell-free assays indicate that the domain carrying 4-9 VDP repeats (4-9 m) showed a variable binding affinity with peptidoglycan and Rv0071/74 cleaves peptidoglycan. Furthermore, Rv0071/74-9 m increased cell wall thickness and reduced the intracellular concentration of antibiotics.

CONCLUSIONS

These findings not only identify Rv0071/74 with VDP repeats as a newly identified multidrug resistance gene but also provide a new model for the development of multiple drug resistance.

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.

Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates

Tuberculosis presents a global health challenge. Mycobacterium tuberculosis is divided into several lineages, each with a different geographical distribution. M. tuberculosis lineage 1 (L1) is common in the high-burden areas in East Africa and Southeast Asia. Although the founder effect contributes significantly to the phylogeographic profile, co-evolution between the host and M. tuberculosis may also play a role. Here, we reported the genomic analysis of 480 L1 isolates from patients in northern Thailand. The studied bacterial population was genetically diverse, allowing the identification of a total of 18 sublineages distributed into three major clades. The majority of isolates belonged to L1.1 followed by L1.2.1 and L1.2.2. Comparison of the single nucleotide variant (SNV) phylogenetic tree and the clades defined by spoligotyping revealed some monophyletic clades representing EAI2_MNL, EAI2_NTM and EAI6_BGD1 spoligotypes. Our work demonstrates that ambiguity in spoligotype assignment could be partially resolved if the entire DR region is investigated. Using the information to map L1 diversity across Southeast Asia highlighted differences in the dominant strain-types in each individual country, despite extensive interactions between populations over time. This finding supported the hypothesis that there is co-evolution between the bacteria and the host, and have implications for tuberculosis disease control.

A complete high-quality MinION nanopore assembly of an extensively drug-resistant Mycobacterium tuberculosis Beijing lineage strain identifies novel variation in repetitive PE/PPE gene regions

A better understanding of the genomic changes that facilitate the emergence and spread of drug-resistant Mycobacterium tuberculosis strains is currently required. Here, we report the use of the MinION nanopore sequencer (Oxford Nanopore Technologies) to sequence and assemble an extensively drug-resistant (XDR) isolate, which is part of a modern Beijing sub-lineage strain, prevalent in Western Province, Papua New Guinea. Using 238-fold coverage obtained from a single flow-cell, de novo assembly of nanopore reads resulted into one contiguous assembly with 99.92 % assembly accuracy. Incorporation of complementary short read sequences (Illumina) as part of consensus error correction resulted in a 4 404 064 bp genome with 99.98 % assembly accuracy. This assembly had an average nucleotide identity of 99.7 % relative to the reference genome, H37Rv. We assembled nearly all GC-rich repetitive PE/PPE family genes (166/168) and identified variants within these genes. With an estimated genotypic error rate of 5.3 % from MinION data, we demonstrated identification of variants to include the conventional drug resistance mutations, and those that contribute to the resistance phenotype (efflux pumps/transporter) and virulence. Reference-based alignment of the assembly allowed detection of deletions and insertions. MinION sequencing provided a fully annotated assembly of a transmissible XDR strain from an endemic setting and showed its utility to provide further understanding of genomic processes within Mycobacterium tuberculosis.

Importance of differential identification of Mycobacterium tuberculosis strains for understanding differences in their prevalence, treatment efficacy, and vaccine development

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a serious global health problem in the 21^st century because of its high mortality. Mtb is an extremely successful human-adapted pathogen that displays a multifactorial ability to control the host immune response and to evade killing by drugs, resulting in the breakdown of BCG vaccine-conferred anti-TB immunity and development of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mtb. Although genetic components of the genomes of the Mtb complex strains are highly conserved, showing over 99% similarity to other bacterial genera, recently accumulated evidence suggests that the genetic diversity of the Mtb complex strains has implications for treatment outcomes, development of MDR/XDR Mtb, BCG vaccine efficacy, transmissibility, and epidemiological outbreaks. Thus, new insights into the pathophysiological features of the Mtb complex strains are required for development of novel vaccines and for control of MDR/XDR Mtb infection, eventually leading to refinement of treatment regimens and the health care system. Many studies have focused on the differential identification of Mtb complex strains belonging to different lineages because of differences in their virulence and geographical dominance. In this review, we discuss the impact of differing genetic characteristics among Mtb complex strains on vaccine efficacy, treatment outcome, development of MDR/XDR Mtb strains, and epidemiological outbreaks by focusing on the best-adapted human Mtb lineages. We further explore the rationale for differential identification of Mtb strains for more effective control of TB in clinical and laboratory settings by scrutinizing current diagnostic methods.

The truncated Rv2820c of Mycobacterium tuberculosis Beijing family augments intracellular survival of M. smegmatis by altering cytokine profile and inhibiting NO generation

Genetic variations among genes of Mycobacterium tuberculosis may be associated with antigenic variation and immune evasion, which complicates the pathogenesis of M. tuberculosis. The hyper-virulent M. tuberculosis Beijing strains harbored several large sequence deletions, among which RD207 attributed to the deletion of CRISPR loci and several Cas genes. RD207 also gave rise to a truncated gene Rv2820c-Bj with 60% deletion in length at the 3'-end and a new 3'-end of five amino acid mutations. It has been reported that Rv2820c-Bj correlated with enhanced intracellular survival of M. smegmatis in macrophages when compared to its full-length counterpart Rv2820c in M. tuberculosis, however, the respective contribution of the truncation and the new 3'-end of Rv2820c-Bj to this enhancement was unclear. Here, by infecting THP-1 macrophages with Ms_Rv2820c-Bj, Ms_Rv2820c and MS_Rv2820c-Tr (expressing the truncated Rv2820c without five amino acid mutations at 3'-end), we found only Ms_Rv2820c-Bj was responsible for the enhancement of survival of M. smegmatis in macrophages. Furthermore, we detected that Ms_Rv2820c-Tr and Ms_Rv2820c-Bj induced similar cytokine profile and NO production after infection of macrophages, which was distinctly different from Ms_Rv2820c. However, Ms_Rv2820c-Bj evoked higher levels of interleukin-10 (IL-10) and lower levels of interleukin- 6 (IL-6), interleukin-1β (IL-1β) and interleukin-12 (IL-12) in infected THP-1 macrophages than Ms_Rv2820c-Tr. Accordingly, we concluded that the new 3'-end of Rv2820c-Bj was important to dampen host defense and enhance the intracellular survival of M. smegmatis.

Influence of Stress and Antibiotic Resistance on Cell-Length Distribution in Mycobacterium tuberculosis Clinical Isolates

Mycobacterial cellular variations in growth and division increase heterogeneity in cell length, possibly contributing to cell-to-cell variation in host and antibiotic stress tolerance. This may be one of the factors influencing Mycobacterium tuberculosis persistence to antibiotics. Tuberculosis (TB) is a major public health problem in developing countries, antibiotic persistence, and emergence of antibiotic resistance further complicates this problem. We wanted to investigate the factors influencing cell-length distribution in clinical M. tuberculosis strains. In parallel we examined M. tuberculosis cell-length distribution in a large set of clinical strains (n = 158) from ex vivo sputum samples, in vitro macrophage models, and in vitro cultures. Our aim was to understand the influence of clinically relevant factors such as host stresses, M. tuberculosis lineages, antibiotic resistance, antibiotic concentrations, and disease severity on the cell size distribution in clinical M. tuberculosis strains. Increased cell size and cell-to-cell variation in cell length were associated with bacteria in sputum and infected macrophages rather than liquid culture. Multidrug-resistant (MDR) strains displayed increased cell length heterogeneity compared to sensitive strains in infected macrophages and also during growth under rifampicin (RIF) treatment. Importantly, increased cell length was also associated with pulmonary TB disease severity. Supporting these findings, individual host stresses, such as oxidative stress and iron deficiency, increased cell-length heterogeneity of M. tuberculosis strains. In addition we also observed synergism between host stress and RIF treatment in increasing cell length in MDR-TB strains. This study has identified some clinical factors contributing to cell-length heterogeneity in clinical M. tuberculosis strains. The role of these cellular adaptations to host and antibiotic tolerance needs further investigation.

Cmr is a redox-responsive regulator of DosR that contributes to M. tuberculosis virulence

Mycobacterium tuberculosis (MTb) is the causative agent of pulmonary tuberculosis (TB). MTb colonizes the human lung, often entering a non-replicating state before progressing to life-threatening active infections. Transcriptional reprogramming is essential for TB pathogenesis. In vitro, Cmr (a member of the CRP/FNR super-family of transcription regulators) bound at a single DNA site to act as a dual regulator of cmr transcription and an activator of the divergent rv1676 gene. Transcriptional profiling and DNA-binding assays suggested that Cmr directly represses dosR expression. The DosR regulon is thought to be involved in establishing latent tuberculosis infections in response to hypoxia and nitric oxide. Accordingly, DNA-binding by Cmr was severely impaired by nitrosation. A cmr mutant was better able to survive a nitrosative stress challenge but was attenuated in a mouse aerosol infection model. The complemented mutant exhibited a ∼2-fold increase in cmr expression, which led to increased sensitivity to nitrosative stress. This, and the inability to restore wild-type behaviour in the infection model, suggests that precise regulation of the cmr locus, which is associated with Region of Difference 150 in hypervirulent Beijing strains of Mtb, is important for TB pathogenesis.

Evolutionary pathway analysis and unified classification of East Asian lineage of Mycobacterium tuberculosis

Due to its rapid spread and association with the numerous outbreaks, the global spread of East Asian lineage of Mycobacterium tuberculosis strains presents a global concern. Although there were many attempts to describe its population structure, no consensus has been reached yet. To define unbiased classification that will facilitate future studies of this lineage, we analyzed the performance and congruence of eight different genotyping schemes based on phylogenetic analysis of 1,398 strains from 32 countries using whole-genome sequencing (WGS) data. We confirm that East Asian lineage comprises two major clades, designated proto-Beijing, which harbors unusual 43-signal spoligoprofile, and Beijing, with well-known spoligoprofile (deleted signals from 1 to 34). We show that different genotyping methods give high consistency results in description of ancient Beijing strains while the classification of modern Beijing strains is significantly divergent due to star-shaped phylogeny. Using WGS data we intersect different studies and for the first time provide balanced classification with well-defined major groups and their genetic markers. Our reconstructed phylogenetic tree can also be used for further analysis of epidemiologically important clusters and their ancestors as well as white spots of unclassified strains, which are prospective areas of research.

Comparative genomic analysis of Mycobacterium tuberculosis Beijing-like strains revealed specific genetic variations associated with virulence and drug resistance

Isolates of the Mycobacterium tuberculosis lineage 2/East-Asian are considered one of the most successful strains due to their increased pathogenicity, hyper-virulence associated with drug resistance, and high transmission. Recent studies in Colombia have shown that the Beijing-like genotype is associated with multidrug-resistance and high prevalence in the southwest of the country, but the genetic basis of its success in dissemination is unknown. In contribution to this matter, we obtained the whole sequences of six genomes of clinical isolates assigned to the Beijing-like genotype. The genomes were compared with the reference genome of M. tuberculosis H37Rv and 53 previously published M. tuberculosis genomes. We found that the six Beijing-like isolates belong to a modern Beijing sub-lineage and share specific genomic variants: i.e. deletion in the PPE8 gene, in Rv3806c (ubiA) responsible of high ethambutol resistance and in Rv3862c (whiB6) which is involved in granuloma formation and virulence, are some of them. Moreover, each isolated has exclusively single nucleotide polymorphisms (SNPs) in genes related with cell wall processes and cell metabolism. We identified polymorphisms in genes related to drug resistance that could explain the drug-resistant phenotypes found in the six isolates from Colombia. We hypothesize that changes due to these genetic variations contribute to the success of these strains. Finally, we analyzed the IS6110 insertion sequences finding very low variance between them, suggesting that SNPs is the major cause of variability found in Beijing-like strains circulating in Colombia.

First Insight into the Molecular Epidemiology of Mycobacterium tuberculosis Isolates from the Minority Enclaves of Southwestern China

Yunnan is a province located in southwestern China. The aim of this study was to investigate the genetic diversity of Mycobacterium tuberculosis (MTB) strains circulating in Yunnan Province. We used spoligotyping and a 12-locus high-resolution VNTR set to identify a total of 271 MTB isolates collected from six prefectures' Yunnan. All the 271 patients were classified as 11 different ethnic groups, including 133 Han patients (49.1%) and 138 minority patients (50.9%). Spoligotyping analyses revealed that the largest two spoligotypes were SIT1 (Beijing family, n = 136) and SIT53 (T family, n = 35). Statistical analysis indicated that the proportion of Beijing genotype in Qujing was significantly higher than that in the ethnic enclaves (P ≤ 0.01). Compared with the percentage of Beijing-family isolates from patients under 25 years of age (72.7%), those from patients aged 45~64 years (44.6%, P < 0.01) and older than 64 years (48.1%) (P = 0.04) were significantly lower. Beijing genotype strains (23.8%, 36/151) showed higher clustering rate than non-Beijing genotype strains (16/120, 13.3%, P = 0.029). In conclusion, our data demonstrated that the Beijing genotype was the predominant genotype in Yunnan Province. The distribution of Beijing genotype strains showed geographic diversity. In addition, Beijing genotype was associated with recent transmission rather than drug resistance.

Rapid and simple identification of Beijing genotype strain of Mycobacterium tuberculosis using a loop-mediated isothermal amplification assay

Beijing genotype strains of Mycobacterium tuberculosis are geographically widespread and pose a notorious public health problem, these strains causing outbreaks of multidrug-resistant tuberculosis (TB); some studies have reported an association with drug resistance. Because the prevalence of Beijing strain has a substantial impact on TB control programs, the availability of a rapid and reliable method for detecting these strains is important for epidemiological monitoring of their circulation. The main methods currently used to identify Beijing genotype strains are IS6110 DNA fingerprinting, spoligotyping and PCR to detect specific deletions such as region of difference (RD)207. More recently, multiplex PCR assay using a Beijing-specific single nucleotide polymorphism (SNP) has been developed for detecting Beijing lineage strains. However, these methods are time-consuming and technically demanding. In the present study, a loop-mediated isothermal amplification (LAMP) assay that allows specific identification of Beijing genotype strain was developed. This Beijing genotype strain-identifying LAMP assay was performed 214 clinical isolates and the results compared with those of conventional PCR that targeted RD207 and Rv0679c-targreting multiplex PCR for Beijing lineage identification. LAMP assay showed 100% sensitivity and specificity compared with RD207-PCR. Furthermore, the sensitivity and specificity were 99.3% and 100%, respectively, compared with Rv0679c-multiplex PCR. This LAMP assay could be used routinely in local laboratories to monitor the prevalence of the Beijing genotype strain and thereby used to help control the spread of these potentially highly virulent and drug resistant strains.

Structure and variation of CRISPR and CRISPR-flanking regions in deleted-direct repeat region Mycobacterium tuberculosis complex strains

Background

CRISPR and CRISPR-flanking genomic regions are important for molecular epidemiology of Mycobacterium tuberculosis complex (MTBC) strains, and potentially for adaptive immunity to phage and plasmid DNA, and endogenous roles in the bacterium. Genotyping in the Israel National Mycobacterium Reference Center Tel-Aviv of over 1500 MTBC strains from 2008–2013 showed three strains with validated negative 43-spacer spoligotypes, that is, with putatively deleted direct repeat regions (deleted-DR/CRISPR regions). Two isolates of each of three negative spoligotype MTBC (a total of 6 isolates) were subjected to Next Generation Sequencing (NGS). As positive controls, NGS was performed for three intact-DR isolates belonging to T3_Eth, the largest multiple-drug-resistant (MDR)-containing African-origin cluster in Israel. Other controls consisted of NGS reads and complete whole genome sequences from GenBank for 20 intact-DR MTBC and for 1 deleted-DR MTBC strain recognized as CAS by its defining RD deletion.

Results

NGS reads from negative spoligotype MTBC mapped to reference H37Rv NC_000962.3 suggested that the DR/CRISPR regions were completely deleted except for retention of the middle IS6110 mobile element. Clonally specific deletion of CRISPR-flanking genes also was observed, including deletion of at least cas2 and cas1 genes. Genomic RD deletions defined lineages corresponding to the major spoligotype families Beijing, EAI, and Haarlem, consistent with 24 loci MIRU-VNTR profiles. Analysis of NGS reads, and analysis of contigs obtained by manual PCR confirmed that all 43 gold standard DR/CRISPR spacers were missing in the deleted-DR genomes.

Conclusions

Although many negative spoligotype strains are recorded as spoligotype-international-type (SIT) 2669 in the SITVIT international database, this is the first time to our knowledge that it has been shown that negative spoligotype strains are found in at least 4 different 24 loci MIRU-VNTR and RD deletion families. We report for the first time negative spoligotype-associated total loss of CRISPR region spacers and repeats, with accompanying clonally specific loss of flanking genes, including at least CRISPR-associated genes cas2 and cas1. Since cas1 deleted E.coli shows increased sensitivity to DNA damage and impaired chromosomal segregation, we discussed the possibility of a similar phenotype in the deleted-DR strains and Beijing family strains as both lack the cas1 gene.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3560-6) contains supplementary material, which is available to authorized users.

Unique Regulation of the DosR Regulon in the Beijing Lineage of Mycobacterium tuberculosis

The DosR regulon, a set of 48 genes normally expressed in Mycobacterium tuberculosis under conditions that inhibit aerobic respiration, is controlled via the DosR-DosS/DosT two-component system. While the regulon requires induction in most M. tuberculosis isolates, for members of the Beijing lineage, its expression is uncoupled from the need for signaling. In our attempts to understand the mechanistic basis for this uncoupling in the Beijing background, we previously reported the identification of two synonymous single-nucleotide polymorphisms (SNPs) within the adjacent Rv3134c gene. In the present study, we have interrogated the impact of these SNPs on dosR expression in wild-type strains, as well as a range of dosR-dosS-dosT mutants, for both Beijing and non-Beijing M. tuberculosis backgrounds. In this manner, we have unequivocally determined that the C601T dosR promoter SNP is the sole requirement for the dramatic shift in the pattern of DosR regulon expression seen in this globally important lineage. Interestingly, we also show that DosT is completely nonfunctional within these strains. Thus, a complex series of evolutionary steps has led to the present-day Beijing DosR phenotype that, in turn, potentially confers a fitness advantage in the face of some form of host-associated selective pressure.

IMPORTANCE

Mycobacterium tuberculosis strains of the Beijing lineage have been described as being of enhanced virulence compared to other lineages, and in certain regions, they are associated with the dramatic spread of multidrug-resistant tuberculosis (TB). In terms of trying to understand the functional basis for these broad epidemiological phenomena, it is interesting that, in contrast to the other major lineages, the Beijing strains all constitutively overexpress members of the DosR regulon. Here, we identify the mutational events that led to the evolution of this unique phenotype. In addition, our work highlights the fact that important phenotypic differences exist between distinct M. tuberculosis lineages, with the potential to impact the efficacy of diagnosis, vaccination, and treatment programs.

Risk factors for pulmonary cavitation in tuberculosis patients from China

Pulmonary cavitation is one of the most frequently observed clinical characteristics in tuberculosis (TB). The objective of this study was to investigate the potential risk factors associated with cavitary TB in China. A total of 385 smear-positive patients were enrolled in the study, including 192 (49.9%) patients with cavitation as determined by radiographic findings. Statistical analysis revealed that the distribution of patients with diabetes in the cavitary group was significantly higher than that in the non-cavitary group (adjusted odds ratio (OR) (95% confidence interval (CI)):12.08 (5.75–25.35), P<0.001). Similarly, we also found that the proportion of individuals with multidrug-resistant TB in the cavitary group was also higher than that in the non-cavitary group (adjusted OR (95% CI): 2.48 (1.52–4.07), P<0.001). Of the 385 Mycobacterium tuberculosis strains, 330 strains (85.7%) were classified as the Beijing genotype, which included 260 strains that belonged to the modern Beijing sublineage and 70 to the ancient Beijing sublineage. In addition, there were 80 and 31 strains belonging to large and small clusters, respectively. Statistical analysis revealed that cavitary disease was observed more frequently among the large clusters than the small clusters (P=0.037). In conclusion, our findings demonstrate that diabetes and multidrug resistance are risk factors associated with cavitary TB. In addition, there was no significant difference in the cavitary presentation between patients infected with the Beijing genotype strains and those infected with the non-Beijing genotype strains.

Parallel reaction monitoring of clinical Mycobacterium tuberculosis lineages reveals pre-existent markers of rifampicin tolerance in the emerging Beijing lineage

The spread of multidrug resistant Mycobacterium tuberculosis is one of the major challenges in tuberculosis control. In Eurasia, the spread of multidrug resistant tuberculosis is driven by the M. tuberculosis Beijing genotype. In this study, we examined whether selective advantages are present in the proteome of Beijing isolates that contribute to the emergence of this genotype. To this end, we compared the proteome of M. tuberculosis Beijing to that of M. tuberculosis H37Rv, both in the presence and absence of the first-line antibiotic rifampicin. During rifampicin exposure, both M. tuberculosis genotypes express proteins belonging to the DosR dormancy regulon, which induces a metabolically hypoactive-, drug tolerant phenotype. However, these markers of rifampicin tolerance were already more abundant in the M. tuberculosis Beijing isolate prior to drug exposure. To determine whether the a priori high abundance of specific proteins contribute to the formation of antibiotic resistance in M. tuberculosis Beijing, we quantified the abundance of 33 selected proteins in 27 clinical isolates from the five most common M. tuberculosis lineages using parallel reaction monitoring. The observed pre-existing high abundance of dormancy proteins in Beijing strains provides an evolutionary advantage that allows these strains to persist for prolonged periods during rifampicin treatment.

SIGNIFICANCE

M. tuberculosis is the leading cause of death by a bacterial infection worldwide. Treatment-regimen to eradicate this pathogen make use of the first-line antibiotic rifampicin, which is considered to be the cornerstone of modern day anti-tuberculosis treatment. Despite the potency of rifampicin, there is an increasing occurrence of rifampicin resistant mutants in a specific cluster of M. tuberculosis, the Beijing genotype. Using both a data dependent acquisition and a targeted proteomic approach we identified markers of rifampicin tolerance to be high abundant in members of the M. tuberculosis Beijing genotype, already prior drug exposure. The identification of this M. tuberculosis Beijing specific trait will contribute to improved diagnostics and treatment of M. tuberculosis.

Variable transcriptional adaptation between the laboratory (H37Rv) and clinical strains (S7 and S10) of Mycobacterium tuberculosis under hypoxia

Tuberculosis continues to be a major public health problem in many parts of the world, despite intensified efforts taken to control the disease. The remarkable success of M. tuberculosis as a pathogen is largely due to its ability to persist within the host for long periods. To develop the effective intervention strategies, understanding the biology of persistence is highly required. Accumulating evidences showed oxygen deprivation (hypoxia) as a potential stimulus for triggering the transition of M. tuberculosis to a non-replicating persistent state analogous to latency in vivo. To date, in vitro hypoxia experimental models used the laboratory adapted isolate H37Rv and very little is known about the behavior of clinical isolates that are involved during disease outbreaks. Hence, we compared the transcription profiles of H37Rv and two south Indian clinical isolates (S7 and S10) under hypoxia to find differences in gene expression pattern. The main objective of this current work is to find "differentially regulated genes" (genes that are down regulated in H37Rv but upregulated in both the clinical isolates) under hypoxia. Microarray results showed, a total of 502 genes were down regulated in H37Rv under hypoxia and 10 out of 502 genes were upregulated in both the clinical isolates. Thus, giving less importance to down regulated genes based on H37Rv model strain might exclude the true representative gene candidates in clinical isolates. Our study suggests the use of most prevalent clinical isolates for in vitro experimental model to minimize the variation in understanding the adaptation mechanisms of the strains.

A rapid PCR assay to characterize the intact pks15/1 gene, a virulence marker in Mycobacterium tuberculosis

Intact pks15/1 is involved in the biosynthesis of phenolic glycolipids and proposed as a marker for virulence and phylogeny in tuberculosis. Identification of intact condition is achieved mainly by DNA sequencing. For this reason the aim of this study was to develop a reproducible endpoint PCR-assay to characterize it.

Host and Microbial Predictors of Childhood Extrathoracic Tuberculosis and Tuberculosis Meningitis

BACKGROUND

Childhood tuberculosis (TB) is a major infectious disease causing morbidity and mortality in children and yet has been largely ignored until recently. This study is the first study to characterize childhood TB in China incorporating both Mycobacterium tuberculosis genetic characteristics and patient data.

METHODS

We analyzed a total of 331 culture-confirmed childhood TB cases and 158 M. tuberculosis isolates from a subset of the study sample. Univariate and multivariate logistic regression analyses were performed to identify host and microbial predictors for having extrathoracic TB alone, concurrent thoracic and extrathoracic TB and TB meningitis (TBM), respectively.

RESULTS

Fifty-eight percent (192/331) of the study subjects had extrathoracic TB, and 139 (42.0%) cases had TBM. Both age of less than 5 years [adjusted odds ratio (OR): 4.52; 95% confidence interval (CI): 1.27-16.16] and female sex (adjusted OR: 2.72; 95% CI: 1.03-7.18) were significantly associated with extrathoracic TB alone, whereas living in rural area (adjusted OR: 2.35; 95% CI: 1.06-5.18) was significantly associated with thoracic-extrathoracic TB. Age of less than 5 years was also strongly associated with TBM (adjusted OR: 3.63; 95% CI: 1.64-8.05). Sixty-four percent (101/158) of the study isolates were Beijing lineage strains. Infection with Beijing lineage strains was significantly associated with thoracic-extrathoracic TB (adjusted OR: 2.39; 95% CI: 1.11-5.15) and TBM (adjusted OR: 2.25; 95% CI: 1.10-4.60).

CONCLUSIONS

Both microbial and host factors can affect the outcome of M. tuberculosis infection in children. Future studies incorporating host and pathogen data from different populations are warranted to develop new strategies for childhood TB control.

Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese

The Beijing family is the most successful genotype of Mycobacterium tuberculosis and responsible for more than a quarter of the global tuberculosis epidemic. As the predominant genotype in East Asia, the Beijing family has been emerging in various areas of the world and is often associated with disease outbreaks and antibiotic resistance. Revealing the origin and historical dissemination of this strain family is important for understanding its current global success. Here we characterized the global diversity of this family based on whole-genome sequences of 358 Beijing strains. We show that the Beijing strains endemic in East Asia are genetically diverse, whereas the globally emerging strains mostly belong to a more homogenous subtype known as "modern" Beijing. Phylogeographic and coalescent analyses indicate that the Beijing family most likely emerged around 30,000 y ago in southern East Asia, and accompanied the early colonization by modern humans in this area. By combining the genomic data and genotyping result of 1,793 strains from across China, we found the "modern" Beijing sublineage experienced massive expansions in northern China during the Neolithic era and subsequently spread to other regions following the migration of Han Chinese. Our results support a parallel evolution of the Beijing family and modern humans in East Asia. The dominance of the "modern" Beijing sublineage in East Asia and its recent global emergence are most likely driven by its hypervirulence, which might reflect adaption to increased human population densities linked to the agricultural transition in northern China.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.