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

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.

Transmission of drug-resistant Mycobacterium tuberculosis isolates between Finnish- and foreign-born cases, 2014-2021: A molecular epidemiological study

BACKGROUND

Data on the molecular epidemiology and transmission of drug-resistant Mycobacterium tuberculosis (MTB) in low-incidence settings with immigration from high-incidence settings is limited.

METHOD

We included 115 drug-resistant (DR) MTB isolates with whole-genome sequencing data isolated in Finland between 2014 and 2021. Potential transmission clusters were identified using a threshold of 12 single-nucleotide polymorphisms (SNPs). Highly related clusters were identified using a threshold of 5 SNPs.

RESULT

Of the 115 DR MTB isolates, 31 (27.0%) isolates were from Finnish-born cases and 84 (73.0%) were from foreign-born cases. The proportion of multidrug-resistant (MDR) MTB isolates (30/84, 35.7%) from foreign-born cases was higher than that of MDR MTB isolates from Finnish-born cases (8/31, 25.8%). Lineage 2 (40/115, 34.8%) and lineage 4 (40/115, 34.8%) were the most prevalent lineages. A total of 25 (21.7%) isolates were classified into eight potential transmission clusters (≤12 SNPs). Furthermore, five highly related clusters (≤5 SNPs) were identified, including three DR MTB isolates from Finnish-born cases and 14 DR isolates from foreign-born cases.

CONCLUSION

The risk of DR MTB transmission between Finnish- and foreign-born persons is not negligible. Further research on clustering analysis in drug-susceptible MTB is worth to inform tuberculosis management and control in low-incidence settings with increasing immigration.

The Recent Transmission and Associated Risk Factor of Mycobacterium tuberculosis in Golmud City, China

Background

Tuberculosis (TB) remains a severe public health problem globally, and it is essential to comprehend the transmission pattern to control tuberculosis. Herein, we evaluated the drug-resistant characteristics, recent transmission, and associated risk factors of TB in Golmud, Qinghai, China.

Methods

In this study, we performed a population-based study of patients diagnosed with TB in Golmud from 2013 to 2018. Drug-susceptibility testing and whole-genome sequencing were performed on 133 Mycobacterium tuberculosis strains. The genomic clustering rate was calculated to evaluate the level of recent transmission. Risk factors were identified by logistic regression analysis.

Results

Our results showed that 46.97% (62/132) of strains were phylogenetically clustered and formed into 23 transmission clusters, suggesting a high recent transmission of TB in the area. 12.78% (17/133) strains were multidrug-resistant/rifampicin tuberculosis (MDR/RR-TB), with a high drug-resistant burden. Based on drug resistance gene analysis, we found 23 strains belonging to genotype MDR/RR-TB, where some strains may have borderline mutations. Among these strains, 65.2% (15/23) were found within putative transmission clusters. Additionally, risk factor analysis showed that recent transmission of TB happened more in patients with Tibetan nationality or older age.

Conclusion

Overall our study indicates that the recent transmissions of MTB strains, especially genotypic MDR/RR strains, drive the tuberculosis epidemic in Golmud, which could contribute to developing effective TB prevention and control strategies.

Evolution and transmission of antibiotic resistance is driven by Beijing lineage Mycobacterium tuberculosis in Vietnam

IMPORTANCE

Drug-resistant tuberculosis (TB) infection is a growing and potent concern, and combating it will be necessary to achieve the WHO's goal of a 95% reduction in TB deaths by 2035. While prior studies have explored the evolution and spread of drug resistance, we still lack a clear understanding of the fitness costs (if any) imposed by resistance-conferring mutations and the role that Mtb genetic lineage plays in determining the likelihood of resistance evolution. This study offers insight into these questions by assessing the dynamics of resistance evolution in a high-burden Southeast Asian setting with a diverse lineage composition. It demonstrates that there are clear lineage-specific differences in the dynamics of resistance acquisition and transmission and shows that different lineages evolve resistance via characteristic mutational pathways.

Towards the reconstruction of a global TB history using a new pipeline "TB-Annotator"

Mycobacterium tuberculosis complex (MTBC) has a population structure consisting of 9 human and animal lineages. The genomic diversity within these lineages is a pathogenesis factor that affects virulence, transmissibility, host response, and antibiotic resistance. Hence it is important to develop improved information systems for tracking and understanding the spreading and evolution of genomes. We present results obtained thanks to a new informatics platform for computational biology of MTBC, that uses a convenience sample from public/private SRAs, designated as TB-Annotator. Version 1 was a first interactive graphic-based web tool based on 15,901 representative genomes. Version 2, still interactive, is a more sophisticated database, developed using the Snakemake Workflow Management System (WMS) that allows an unsupervised global and scalable analysis of the content of the USA National Center for Biotechnology Information Short Read Archives database. This platform analyzes nucleotide variants, the presence/absence of genes, known regions of difference and detect new deletions, the insertion sites of mobile genetic elements, and allows phylogenetic trees to be built, imported in a graphical interface and interactively analyzed between the data and the tree. The objective of TB-Annotator is triple: detect recent epidemiological links, reconstruct distant phylogeographical histories as well as perform more complex phenotypic/genotypic Genome-Wide Association Studies (GWAS). In this paper, we compare the various taxonomic SNPs-based labels and hierarchies previously described in recent reference papers for L1, and present a comparative analysis that allows identification of alias and thus provides the basis of a future unifying naming scheme for L1 sublineages. We present a global phylogenetic tree built with RAxML-NG, and one on L2; at the time of writing, we characterized about 200 sublineages, with many new ones; a detail tree for Modern L2 and a hierarchical scheme allowing to facilitate L2 lineage assignment are also presented.

A revised SNP-based barcoding scheme for typing Mycobacterium tuberculosis complex isolates

The development of whole-genome sequencing technologies is gradually leading to a more detailed description of the population structure of the Mycobacterium tuberculosis complex (MTBC). In this study, we correlated previously published classifications on a collection of more than 10,000 genomes and proposed a new, comprehensive nomenclature that unifies the existing ones. In total, we identified 169 lineages and sublineages of M. tuberculosis/M. africanum and 9 animal-adapted species. For the purpose of organizing these genotypes in a more streamlined manner, we stratified them into five hierarchical levels. To represent the classification and compare it with the reference, we compiled a confirmatory data set of 670 high-quality isolates, which includes all genotypes and species of MTBC, and this confirmatory data set can serve as a basis for further studies. We proposed a set of 213 robust barcoding single-nucleotide polymorphisms and a suitable workflow for reliable differentiation of genotypes and species within the complex. This work integrates the results of all the major systematized studies to date to provide an understanding of the global diversity of the MTBC population structure. The results of this work may ultimately help to reliably determine the pathogen genotype and associate it with traits that reflect its prevalence, virulence, vaccination, and treatment efficiency, as well as to reliably find natural features revealed during its spread. IMPORTANCE Through years of research into the Mycobacterium tuberculosis complex (MTBC), a number of ambiguous phylogenetic classifications have emerged, which often overlap with one another. In the present study, we have combined all major studies on MTBC classification and inferred a unified, most complete to date classification and accompanying SNP barcodes.

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.

Whole Genomic Analysis Revealed High Genetic Diversity and Drug-Resistant Characteristics of Mycobacterium tuberculosis in Guangxi, China

Background

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is a major public health issue in China. Nevertheless, the prevalence and drug resistance characteristics of isolates vary in different regions and provinces. In this study, we investigated the population structure, transmission dynamics and drug-resistant profiles of Mtb in Guangxi, located on the border of China.

Methods

From February 2016 to April 2017, 462 clinical M. tuberculosis isolates were selected from 5 locations in Guangxi. Drug-susceptibility testing was performed using 6 common anti-tuberculosis drugs. The genotypic drug resistance and transmission dynamics were analyzed by the whole genome sequence.

Results

Our data showed that the Mtb in Guangxi has high genetic diversity including Lineage 1 to Lineage 4, and mostly belong to Lineage 2 and Lineage 4. Novelty, 9.6% of Lineage 2 isolates were proto-Beijing genotype (L2.1), which is rare in China. About 12.6% of isolates were phylogenetically clustered and formed into 28 transmission clusters. We observed that the isolates with the high resistant rate of isoniazid (INH, 21.2%), followed by rifampicin (RIF, 13.2%), and 6.7%, 12.1%, 6.7% and 1.9% isolates were resistant to ethambutol (EMB), streptomycin (SM), ofloxacin (OFL) and kanamycin (KAN), respectively. Among these, 6.5% and 3.3% of isolates belong to MDR-TB and Pre-XDR, respectively, with a high drug-resistant burden. Genetic analysis identified the most frequently encountered mutations of INH, RIF, EMB, SM, OFL and KAN were katG_Ser315Thr (62.2%), rpoB_Ser450Leu (42.6%), embB_Met306Vol (45.2%), rpsL_Lys43Arg (53.6%), gyrA_Asp94Gly (29.0%) and rrs_A1401G (66.7%), respectively. Additionally, we discovered that isolates from border cities are more likely to be drug-resistant than isolates from non-border cities.

Conclusion

Our findings provide a deep analysis of the genomic population characteristics and drug-resistant of M. tuberculosis in Guangxi, which could contribute to developing effective TB prevention and control strategies.

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

BACKGROUND

. The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern.

RESULTS

. The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes).

CONCLUSIONS

. We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

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

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

The relative transmission fitness of multidrug-resistant Mycobacterium tuberculosis in a drug resistance hotspot

Multidrug-resistant tuberculosis (MDR-TB) is among the most frequent causes of death due to antimicrobial resistance. Although only 3% of global TB cases are MDR, geographical hotspots with up to 40% of MDR-TB have been observed in countries of the former Soviet Union. While the quality of TB control and patient-related factors are known contributors to such hotspots, the role of the pathogen remains unclear. Here we show that in the country of Georgia, a known hotspot of MDR-TB, MDR Mycobacterium tuberculosis strains of lineage 4 (L4) transmit less than their drug-susceptible counterparts, whereas most MDR strains of L2 suffer no such defect. Our findings further indicate that the high transmission fitness of these L2 strains results from epistatic interactions between the rifampicin resistance-conferring mutation RpoB S450L, compensatory mutations in the RNA polymerase, and other pre-existing genetic features of L2/Beijing clones that circulate in Georgia. We conclude that the transmission fitness of MDR M. tuberculosis strains is heterogeneous, but can be as high as drug-susceptible forms, and that such highly drug-resistant and transmissible strains contribute to the emergence and maintenance of hotspots of MDR-TB. As these strains successfully overcome the metabolic burden of drug resistance, and given the ongoing rollout of new treatment regimens against MDR-TB, proper surveillance should be implemented to prevent these strains from acquiring resistance to the additional drugs.

Phylogenetic relationships of Mycobacterium tuberculosis isolates in Poland: The emergence of Beijing genotype among multidrug-resistant cases

Introduction

The epidemiological situation of tuberculosis (TB) in Poland urges for its continuous and scrupulous monitoring. The objective of this study was to explore the genetic diversity of multidrug-resistant (MDR) and drug-susceptible (DS) Mycobacterium tuberculosis isolates from Poland with a combination of spoligotyping and high-resolution mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis. The results were placed in the Northern and Eastern Europe context.

Methods

The study included 89 (39 MDR and 50 DS) M. tuberculosis isolates collected from as many patients between 2018 and 2021 in Poland. The analysis was done using spoligotyping, and MIRU-VNTR typing at 24 standard loci. The data were compared to those available on Poland and neighbors and global M. tuberculosis datasets.

Results

The main identified families were Beijing (28.1%) and Haarlem (16.8%) while 34.8% of isolates were in the heterogeneous L4-unclassified group. Although the Beijing family was the most prevalent (61.5%) among MDR-TB cases, it accounted for only 2% of DS isolates. Among foreign-born patients, a higher ratio of MDR isolates were observed when compared with those who Poland-born (64.3% vs. 40%). Furthermore, all patients from the Former Soviet Union (FSU) countries were infected with MDR-TB.

Discussion

Whereas DS M. tuberculosis population in Poland is dominated by L4 isolates, MDR isolates are mostly of the Beijing genotype. The rise in the prevalence of the Beijing isolates in Poland, coupled with high proportion of the Beijing genotype among foreign-born TB patients may reflect an ongoing transmission of this family, imported to Poland mainly from FSU countries.

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

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.

Genotypic and phenotypic characteristics of <i>Mycobacterium tuberculosis</i> drug resistance in TB children

Background. Russian Federation is included in the list of 30 countries with the highest burden of tuberculosis, including MDR tuberculosis. The most important part of this problem is the primary MDR/XDR TB in children.The aim: a comparative analysis of the phenotypic and genotypic profile of drug resistance to anti-tuberculosis drugs (ATP) according to whole genome sequencing of M. tuberculosis strains from children.Materials and methods. Whole genome sequencing (WGS) results of 61 M. tuberculosis isolates from children with tuberculosis in 2006–2020 in the Russian Federation were analyzed for anti-TB drug resistance mutations, according to the WHO catalog and were compared with the results of phenotypic drug sensitivity.Results. The M. tuberculosis belonged to two genetic groups: Beijing genotype – 82 % (50/61) dominant Central Asian Russian (31/50) and B0/W148 (16/50) subtypes, and non-Beijing (Ural, S, LAM) – 18 % (11/61). Three isolates belonged to Asian Ancestral subtype (3/50). Of the 61 isolates, only 14.7 % (9/61) were sensitive to antiTB drugs, 49.2 % (30/61) were MDR and 14.7 % (9/61) were pre-XDR. Comparison of the resistance profile (MDR/pre-XDR) with genotype revealed an upward shift for Beijing isolates, in particular Beijing B0/W148 (15/16) subline compared to other Beijing (19/34) (Chi-square with Yates correction = 5.535; p < 0.05) and nonBeijing (5/12) (Chi-square with Yates correction = 6.741; p < 0.05) subtypes. Discrepancies between genotypic and phenotypic drug resistance profiles were found in 11.5 % (7/61) of cases.Conclusions. Based on the analysis of WGS data, the genotypic characteristics of M. tuberculosis and the most complete set of drug resistance mutations were obtained, indicating a significant prevalence in MDR and pre-XDR TB of cases caused by epidemic subtypes of Beijing (B0/W148 and Central Asian Russian). The molecular mechanisms of adaptation of M. tuberculosis to the treatment of anti-TB drugs are not unique for the child population but reflect the general processes of the spread of MDR/XDR in Russia.

Investigating the Diversity of Tuberculosis Spoligotypes with Dimensionality Reduction and Graph Theory

The spoligotype is a graphical description of the CRISPR locus present in Mycobacterium tuberculosis, which has the particularity of having only 68 possible spacers. This spoligotype, which can be easily obtained either in vitro or in silico, allows to have a summary information of lineage or even antibiotic resistance (when known to be associated to a particular cluster) at a lower cost. The objective of this article is to show that this representation is richer than it seems, and that it is under-exploited until now. We first recall an original way to represent these spoligotypes as points in the plane, allowing to highlight possible sub-lineages, particularities in the animal strains, etc. This graphical representation shows clusters and a skeleton in the form of a graph, which led us to see these spoligotypes as vertices of an unconnected directed graph. In this paper, we therefore propose to exploit in detail the description of the variety of spoligotypes using a graph, and we show to what extent such a description can be informative.

Whole Genome Sequencing Assessing Impact of Diabetes Mellitus on Tuberculosis Mutations and Type of Recurrence in India

BACKGROUND

Evidence describing the impact of diabetes mellitus (DM) on the recurrence and mutation rate of Mycobacterium tuberculosis (Mtb) is limited.

METHODS

This study was nested in 3 cohort studies of tuberculosis (TB) patients with and without DM in India. Paired Mtb isolates recovered at baseline and treatment failure/recurrence underwent whole genome sequencing. We compared acquisition of single-nucleotide polymorphisms (SNPs), TB drug resistance mutations, and type of recurrence (endogenous reactivation [<8 SNPs] or exogenous reinfection [≥8 SNPs]) by DM status.

RESULTS

Of 1633 enrolled in the 3 parent cohorts, 236 (14.5%) had microbiologically confirmed TB treatment failure/recurrence; 76 Mtb isolate pairs were available for sequencing (22 in TB-DM and 54 in TB-only). The SNP acquisition rate was overall was 0.43 (95% confidence interval [CI], .25-.64) per 1 person-year (PY); 0.77 (95% CI, .40-1.35) per 1 PY, and 0.44 (95% CI, .19-.86) per 1 PY at treatment failure and recurrence, respectively. Significant difference in SNP rates by DM status was seen at recurrence (0.21 [95% CI, .04-.61]) per 1 PY for TB-only vs 1.28 (95% CI, .41-2.98) per 1 PY for TB-DM; P = .02). No significant difference in SNP rates by DM status was observed at treatment failure. Acquired TB drug resistance was seen in 4 of 18 (22%) in TB-DM vs 4 of 45 (9%) in TB-only (P = .21). Thirteen (17%) participants had exogenous reinfection; the reinfection rate at recurrence was 25% (3/12) for TB-DM vs 17% (4/24) in TB-only (P = .66).

CONCLUSIONS

Considerable intrahost Mtb mutation rates were present at recurrence among patients with DM in India. One-fourth of patients with DM had exogenous reinfection at recurrence.

Identification of region of difference and H37Rv-related deletion in Mycobacterium tuberculosis complex by structural variant detection and genome assembly

Mycobacterium tuberculosis complex (MTBC), the main cause of TB in humans and animals, is an extreme example of genetic homogeneity, whereas it is still nevertheless separated into various lineages by numerous typing methods, which differ in phenotype, virulence, geographic distribution, and host preference. The large sequence polymorphism (LSP), incorporating region of difference (RD) and H37Rv-related deletion (RvD), is considered to be a powerful means of constructing phylogenetic relationships within MTBC. Although there have been many studies on LSP already, focusing on the distribution of RDs in MTBC and their impact on MTB phenotypes, a crumb of new lineages or sub-lineages have been excluded and RvDs have received less attention. We, therefore, sampled a dataset of 1,495 strains, containing 113 lineages from the laboratory collection, to screen for RDs and RvDs by structural variant detection and genome assembly, and examined the distribution of RvDs in MTBC, including RvD2, RvD5, and cobF region. Consistent with genealogical delineation by single nucleotide polymorphism (SNP), we identified 125 RDs and 5 RvDs at the species, lineage, or sub-lineage levels. The specificities of RDs and RvDs were further investigated in the remaining 10,218 strains, suggesting that most of them were highly specific to distinct phylogenetic groups, could be used as stable genetic markers in genotyping. More importantly, we identified 34 new lineage or evolutionary branch specific RDs and 2 RvDs, also demonstrated the distribution of known RDs and RvDs in MTBC. This study provides novel details about deletion events that have occurred in distinct phylogenetic groups and may help to understand the genealogical differentiation.

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.

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

Close Related Drug-Resistance Beijing Isolates of Mycobacterium tuberculosis Reveal a Different Transcriptomic Signature in a Murine Disease Progression Model

Mycobacterium tuberculosis (MTB) lineage 2/Beijing is associated with high virulence and drug resistance worldwide. In Colombia, the Beijing genotype has circulated since 1997, predominantly on the pacific coast, with the Beijing-Like SIT-190 being more prevalent. This genotype conforms to a drug-resistant cluster and shows a fatal outcome in patients. To better understand virulence determinants, we performed a transcriptomic analysis with a Beijing-Like SIT-190 isolate (BL-323), and Beijing-Classic SIT-1 isolate (BC-391) in progressive tuberculosis (TB) murine model. Bacterial RNA was extracted from mice lungs on days 3, 14, 28, and 60. On average, 0.6% of the total reads mapped against MTB genomes and of those, 90% against coding genes. The strains were independently associated as determined by hierarchical cluster and multidimensional scaling analysis. Gene ontology showed that in strain BL-323 enriched functions were related to host immune response and hypoxia, while proteolysis and protein folding were enriched in the BC-391 strain. Altogether, our results suggested a differential bacterial transcriptional program when evaluating these two closely related strains. The data presented here could potentially impact the control of this emerging, highly virulent, and drug-resistant genotype.

Emergence of multi-drug-resistant Mycobacterium tuberculosis in Niger: A snapshot based on whole-genome sequencing

Background

Among other West African countries experiencing the high endemicity of deadly tuberculosis, the situation in Niger is poorly evidenced by microbiological investigations.

Methodology/Principal findings

The study of 42 isolates of Mycobacterium tuberculosis from Niger by whole genome sequencing using Illumina iSeq technology yielded four M. tuberculosis lineages: Indo-Oceanic L1 (n = 1) (2.3%), East-Asian (n = 1) (2.3%), East-African Indian L3 (n = 2) (4.7%) and Euro-American L4 (n = 38) (90.4%). The sub-lineage L4.1.3 comprising 18 isolates (47.3%) was predominant, followed by the L4.6.2.2 sub-lineage (Cameroon genotype, n = 13 isolates) (34.2%). Investigating drug resistance profile for 12 antibiotics found 8/42 (19%) pan-susceptible isolates and 34/42 (81%) resistant isolates; with 40/42 (95.2%) isolates being susceptible to clofazimine-bedaquiline.

Conclusions/Significance

These unprecedented data from Niger highlight the dynamics of tuberculosis transmission and drug resistance in Niger and may assist tuberculosis control in this country which continues to support a high burden of tuberculosis.

Tuberculosis is a major public health problem in Niger, in West Africa. Niger has an emerging problem with multidrug-resistant tuberculosis (MDR-TB). Whole genome sequencing was used to understand the epidemiology of tuberculosis and genetics of multi-drug resistance among patients from the regions in Niger. In this study, most isolates of M. tuberculosis from this dataset belonged to the L4.6.2.2 sub-lineage and L4.1.3 sub-lineage within the Euro-American lineage. Thirty-four out of 42 (81%) isolates were detected as resistant isolates. Our study highlights the need for epidemiological surveillance and more concerted efforts to ensure that patients are put through treatment.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

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.

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.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Whole-genome comparative analysis at the lineage/sublineage level discloses relationships between Mycobacterium tuberculosis genotype and clinical phenotype

Background

Human tuberculosis (TB) caused by members of the Mycobacterium tuberculosis complex (MTBC) is the main cause of death among infectious diseases worldwide. Pulmonary TB (PTB) is the most common clinical phenotype of the disease, but some patients develop an extrapulmonary (EPTB) phenotype in which any organ or tissue can be affected. MTBC species include nine phylogenetic lineages, with some appearing globally and others being geographically restricted. EPTB can or not have pulmonary involvement, challenging its diagnosis when lungs are not implicated, thus causing an inadequate treatment. Finding evidence of a specific M. tuberculosis genetic background associated with EPTB is epidemiologically relevant due to the virulent and multidrug-resistant strains isolated from such cases. Until now, the studies conducted to establish associations between M. tuberculosis lineages and PTB/EPTB phenotypes have shown inconsistent results, which are attributed to the strain predominance from specific M. tuberculosis lineages/sublineages in the samples analyzed and the use of low-resolution phylogenetic tools that have impaired sublineage discrimination abilities. The present work elucidates the relationships between the MTBC strain lineages/sublineages and the clinical phenotypes of the disease as well as the antibiotic resistance of the strains.

Methods

To avoid biases, we retrieved the raw genomic reads (RGRs) of all (n = 245) the M. tuberculosis strains worldwide causing EPTB available in databases and an equally representative sample of the RGRs (n = 245) of PTB strains. A multiple alignment was constructed, and a robust maximum likelihood phylogeny based on single-nucleotide polymorphisms was generated, allowing effective strain lineage/sublineage assignment.

Results

A significant Odds Ratio (OR range: 1.8–8.1) association was found between EPTB and the 1.1.1, 1.2.1, 4.1.2.1 and ancestral Beijing sublineages. Additionally, a significant association between PTB with 4.3.1, 4.3.3, and 4.5 and Asian African 2 and Europe/Russia B0/W148 modern Beijing sublineages was found. We also observed a significant association of Lineage 3 strains with multidrug resistance (OR 3.8; 95% CI [1.1–13.6]), as well as between modern Beijing sublineages and antibiotic resistance (OR 4.3; 3.8–8.6). In this work, it was found that intralineage diversity can drive differences in the immune response that triggers the PTB/EPTB phenotype.

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.

A genome epidemiological study of mycobacterium tuberculosis in subpopulations with high and low incidence rate in Guangxi, South China

BACKGROUND

Tuberculosis (TB) is caused by a bacterium called Mycobacterium tuberculosis (Mtb). China is the third in top 8 high TB burden countries and Guangxi is one of the high incidence areas in South China. Determine bacterial factors that affected TB incidence rate is a step toward Ending the TB epidemic.

RESULTS

Genomes of M. tuberculosis cultures from a relatively high and low incidence region in Guangxi have been sequenced. 347 of 358(96.9%) were identified as M. tuberculosis. All the strains belong to Lineage 2 and Lineage 4, except for one in Lineage 1. We found that the genetic structure of the M. tuberculosis population in each county varies enormously. Low incidence rate regions have a lower prevalence of Beijing genotypes than other regions. Four isolates which harbored mutT4-48 also had mutT2-58 mutations. It is suggested that strains from the ancestors of modern Beijing lineage is circulating in Guangxi. Strains of modern Beijing lineage (OR=2.04) were more likely to acquire drug resistances than Lineage 4. Most of the lineage differentiation SNPs are related to cell wall biosynthetic pathways.

CONCLUSIONS

These results provided a higher resolution to better understand the history of transmission of M. tuberculosis from/to South China. And the incidence rate of tuberculosis might be affected by bacterial population structure shaped by demographic history. Our findings also support the hypothesis that Modern Beijing lineage originated in South China.

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.

Drug-Resistant Tuberculosis in Pet Ring-Tailed Lemur, Madagascar

We diagnosed tuberculosis in an illegally wild-captured pet ring-tailed lemur manifesting lethargy, anorexia, and cervical lymphadenopathy. Whole-genome sequencing confirmed the Mycobacterium tuberculosis isolate belonged to lineage 3 and harbored streptomycin resistance. We recommend reverse zoonosis prevention and determination of whether lemurs are able to maintain M. tuberculosis infection.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Estimation of the global burden of Mycobacterium tuberculosis lineage 1

Tuberculosis is still problematic as it affects large numbers of people globally. Mycobacterium tuberculosis Lineage 1 (L1) or Indo Oceanic Lineage, one of widespread major lineages, has a specific geographic distribution and high mortality. It is highly diverse and endemic in several high burden countries. However, studies on the global burden of L1 and its sublineages remain limited. This may lead to the underestimation of the importance of its variance in developing and applying tuberculosis control measures. This study aimed to estimate the number of patients infected with M. tuberculosis L1 and its sublineages worldwide. The proportion of L1 among tuberculosis patients was searched in published reports from countries around the world and the number of patients was calculated based on a WHO report on country incidences and populations. The numbers of patients infected with the five major sublineages, namely L1.1.1, L1.1.2, L1.1.3, L1.2.1, and L1.2.2 were estimated where information was available. It was found that L1 accounted for 28% of global tuberculosis cases in 2012 and 2018. Over 80% of the L1 global burden was in India, the Philippines, Indonesia and Bangladesh, which are also among the countries with highest absolute numbers of tuberculosis patients in the world. Globally, the estimated number of patients infected with M. tuberculosis L1.2.1 and L1.1.2 was over 1.1 million and of patients infected with L1.1.1 was about 200,000. This study demonstrated that L1 contributes significantly to the global burden of tuberculosis. To achieve the End TB Strategy, more attention needs to be paid to the responses of M. tuberculosis L1 to various control measures.

CRISPRbuilder-TB: "CRISPR-builder for tuberculosis". Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA

Mycobacterium tuberculosis complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains.

In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in M. tuberculosis. We first show that the process we set up, that we coin as “CRISPRbuilder-TB” (https://github.com/cguyeux/CRISPRbuilder-TB), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental in vitro approach and the basic in silico spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS6110 insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of M. tuberculosis CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.

In this article, we tackle the bioinformatical issue of the reconstruction of the Mycobacterium tuberculosis complex CRISPR locus using short read sequences without requiring genome assembly. We first show that many complete genomes, as found in public databases and often reconstructed by de novo assemblies, often contain errors on this locus as well as on other repeated sequences. We provide an in-depth description of our new method, designated as ‘CRISPRbuilder-TB’, and we show that our method provides much more exhaustive and reliable information (on DR variants, spacer diversity, global structure) than Crass and CRISPR_detector. The new and unsuspected genomic diversity we detected is described in a companion paper. Scripts are available to adapt the tool to other species.

Proteogenomic Approach for Mycobacterium tuberculosis Investigation

Recent advances in MS/MS technology have made it possible to use proteomic data to predict protein-coding sequences. This approach is called proteogenomics, and it allows to correctly translate start and stop sites and to reveal new open reading frames. Here, we focus on using proteogenomics to improve the annotation of Mycobacteriumtuberculosis strains. We also describe detail procedures of the extraction of proteins and their further preparation for LC-MS/MS analysis and outline the main steps of data analysis.

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.

Genomic evidence supporting the clonal expansion of extensively drug-resistant tuberculosis bacteria belonging to a rare protoBeijing genotype

Tuberculosis disease (TB), caused by Mycobacterium tuberculosis, is a major public health issue in Thailand. The high prevalence of modern Beijing (Lineage 2.2.1) strains has been associated with multi- and extensively drug-resistant infections (MDR-, XDR-TB), complicating disease control. The impact of rarer proto-Beijing (L2.1) strains is less clear. In our study of thirty-seven L2.1 clinical isolates spanning thirteen years, we found a high prevalence of XDR-TB cases (32.4%). With ≤ 12 pairwise SNP distances, 43.2% of L2.1 patients belong to MDR-TB or XDR-TB transmission clusters suggesting a high level of clonal expansion across four Thai provinces. All XDR-TB (100%) were likely due to transmission rather than inadequate treatment. We found a 47 mutation signature and a partial deletion of the fadD14 gene in the circulating XDR-TB cluster, which can be used for surveillance of this rare and resilient M. tuberculosis strain-type that is causing increasing health burden. We also detected three novel deletion positions, a deletion of 1285 bp within desA3 (Rv3230c), large deletions in the plcB, plcA, and ppe38 gene which may play a role in the virulence, pathogenesis or evolution of the L2.1 strain-type.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

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.

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.

Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia

The Mycobacterium tuberculosis Beijing genotype is a clinically and epidemiologically important lineage that is subdivided into ancient/ancestral and modern strains. In our previous study in western Siberia, we identified variable number of tandem repeats (VNTR)-based clusters within the early ancient sublineage of the Beijing genotype characterized by an unexpectedly high rate of extensive drug resistance (XDR). In the current study, next generation sequencing data were analysed to gain insight into genomic signatures underlying drug resistance of these strains. A total of 184 genomes of the Beijing early ancient sublineage from Russia (16), China (15), Japan (36), Korea (25), Vietnam (18), Thailand (73), and the USA (1) were used for phylogenetic analysis. The drug-resistant profile was deduced genotypically. The Russian isolates were distributed into two clusters and were all drug resistant, mainly pre-XDR and XDR. The largest of these clusters included only Russian isolates from remote locations in both Asian and European parts of the country. All its isolates had a quadruple drug resistance (to isoniazid, rifampin, ethambutol and streptomycin) due to the 6-mutation signature (KatG Ser315Thr, KatG Ile335Val, RpoB Ser450Leu, RpoC Asp485Asn, EmbB Gln497Arg, and RpsL Lys43Arg). In most samples, it was complemented with additional and different pncA, gyrA and rrs mutations leading to the pre-XDR/XDR genotype. Phylogenomic analysis indicates a distant origin of this Russian resistant cluster in the early 1970s but location and circumstances are yet to be clarified.

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.

Genetic profiling of Mycobacterium tuberculosis revealed "modern" Beijing strains linked to MDR-TB from Southwestern Colombia

Beijing strains of Mycobacterium tuberculosis (lineage 2) have been associated with drug-resistance and transmission of tuberculosis worldwide. Most of the Beijing strains identified in the Colombian Pacific coast have exhibited a multidrug resistant (MDR) phenotype. We sought to evaluate the clonality and sublineage of Beijing strains circulating in Southwestern Colombia. Thirty-seven Beijing strains were identified through spoligotyping out of 311 clinical isolates collected in 9 years from 2002-2010. Further analysis by MIRU-VNTR 24 loci was conducted for the Beijing strains. For sublineage classification, deletions of RD105, RD207, and RD131 and point mutations at fbpB, mutT2, and acs were evaluated. Drug-resistance associated mutations to first- and second-line anti-TB drugs were also evaluated. Additionally, two Beijing strains were Illumina-whole genome sequenced (one MDR and one drug-susceptible). Among the 37 Beijing strains characterized, 36 belonged to the SIT190 type from which 28 were MDR, four pre-extensively drug resistant (XDR) TB, and four XDR-TB. The remaining strain was SIT1 and drug susceptible. MIRU-VNTR analysis allowed the identification of three Beijing clusters and two unique strains. Beijing strains were confirmed as "modern" sublineage. The mutations rpoB S531L and katG S315T were the most common among MDR strains. Moreover, the two strains evaluated by whole genome sequencing (WGS) shared most of the genetic features with the sublineage 2.2.1 "modern" Beijing previously characterized from Asian strains. WGS analysis of the MDR strain revealed the presence of eight SNPs previously reported in other MDR "Beijing-like" strains from Colombia. The presence of "modern" Beijing strains in Southwestern Colombia, most of them with MDR phenotype, suggests a different origin of this M. tuberculosis sublineage compared to other Beijing strains found in neighboring South American countries. This work may serve as a genetic baseline to study the evolution and spread of M. tuberculosis Beijing strains in Colombia, which play an important role in the propagation of MDR-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.

The epidemic of multidrug resistant tuberculosis in China in historical and phylogenetic perspectives

OBJECTIVES

For the past decade, the epidemic of multidrug resistance tuberculosis (MDR-TB) stays high in China. We investigated the possible driving forces behind the epidemics from phylogenetic and historical perspectives.

METHODS

420 representative strains were selected from the first national drug resistance survey based on their genotypes, drug susceptibility patterns and geographic information. We reconstructed the phylogeny by whole genome sequencing and compared it to the global phylogeny including MDR outbreaks reported in other settings. We estimated the historical trajectory of population dynamics by Bayesian Skygrid plot for all strains and MDR-TB alone. Integrating geographic information and mutations in drug resistance related genes, we investigated the spatial scale of transmission, recent selection of drug resistant mutant, and mechanism for fitness restoration.

RESULTS

Three new subgroups within Beijing clade are described for the first time, but none of the MDR-TB outbreak strains reported in other high MDR-TB burden settings is identified. The overall epidemics experienced two successive phases of expansion at different rates between 1660s and 1950s, followed by a sharp decline till today. Four fifths of the clustered MDR-TB strains suggest transmission of DR strains and nearly half suggest recent selection of (additional) mutations in rpoB. Among all identified transmission events, about one fifth occurred between far distant locations. Possible intergenic and intragenic compensatory mutations both presented in our dataset at comparable frequencies.

CONCLUSIONS

MDR-TB epidemic in China is not yet driven by the spread of a few highly successful clonal expansions but by repeated emergence of smaller and currently less successful clusters. However, internal migration and undertreatment could escalate MDR-TB epidemic. To prevent generating of drug resistance and restoration of fitness as well as to stop transmission of MDR-TB at early stage, national TB control program needs to strengthen management of floating populations and promote universal drug susceptibility testing in China.