Newly Identified Mycobacterium africanum Lineage 10, Central Africa.

Analysis of genome sequencing data from >100,000 genomes of Mycobacterium tuberculosis complex using TB-Annotator software revealed a previously unknown lineage, proposed name L10, in central Africa. Phylogenetic reconstruction suggests L10 could represent a missing link in the evolutionary and geographic migration histories of M. africanum.

Towards a better understanding of the long-lasting evolutionary history of Mycobacterium tuberculosis.

The daily increasing sequencing of Mycobacterium tuberculosis has made it possible to establish an advanced phylogeny of this bacterium. It currently includes 9 lineages mainly affecting humans, completed by animal lineages, which form the Mycobacterium tuberculosis complex. Inherited from various historical approaches, this phylogeny is now based on Single Nucleotide Polymorphisms (SNPs), of which updates are frequently proposed. We present here evidence that the task needs refinements: some lineages have currently suboptimal defining SNPs, and many sublineages still need to be named and characterized. These findings are based on a new tool specifically designed to index the entire existing sequencing data. In this article, we focus on lineages 4.5, 4.7, 6 and 7. We take the opportunity to present some evidence that TB-annotator shows strong relevance, identifying well supported sublineages, as well as good global agreement with previous findings.

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.

Comparison of in silico predicted Mycobacterium tuberculosis spoligotypes and lineages from whole genome sequencing data.

Bacterial strain-types in the Mycobacterium tuberculosis complex underlie tuberculosis disease, and have been associated with drug resistance, transmissibility, virulence, and host-pathogen interactions. Spoligotyping was developed as a molecular genotyping technique used to determine strain-types, though recent advances in whole genome sequencing (WGS) technology have led to their characterization using SNP-based sub-lineage nomenclature. Notwithstanding, spoligotyping remains an important tool and there is a need to study the congruence between spoligotyping-based and SNP-based sub-lineage assignation. To achieve this, an in silico spoligotype prediction method ("Spolpred2") was developed and integrated into TB-Profiler. Lineage and spoligotype predictions were generated for > 28 k isolates and the overlap between strain-types was characterized. Major spoligotype families detected were Beijing (25.6%), T (18.6%), LAM (13.1%), CAS (9.4%), and EAI (8.3%), and these broadly followed known geographic distributions. Most spoligotypes were perfectly correlated with the main MTBC lineages (L1-L7, plus animal). Conversely, at lower levels of the sub-lineage system, the relationship breaks down, with only 65% of spoligotypes being perfectly associated with a sub-lineage at the second or subsequent levels of the hierarchy. Our work supports the use of spoligotyping (membrane or WGS-based) for low-resolution surveillance, and WGS or SNP-based systems for higher-resolution studies.

The future of CRISPR in Mycobacterium tuberculosis infection.

Clustered Regularly Interspaced Short Palindromic repeats (CRISPR)-Cas systems rapidly raised from a bacterial genetic curiosity to the most popular tool for genetic modifications which revolutionized the study of microbial physiology. Due to the highly conserved nature of the CRISPR locus in Mycobacterium tuberculosis, the etiological agent of one of the deadliest infectious diseases globally, initially, little attention was paid to its CRISPR locus, other than as a phylogenetic marker. Recent research shows that M. tuberculosis has a partially functional Type III CRISPR, which provides a defense mechanism against foreign genetic elements mediated by the ancillary RNAse Csm6. With the advent of CRISPR-Cas based gene edition technologies, our possibilities to explore the biology of M. tuberculosis and its interaction with the host immune system are boosted. CRISPR-based diagnostic methods can lower the detection threshold to femtomolar levels, which could contribute to the diagnosis of the still elusive paucibacillary and extrapulmonary tuberculosis cases. In addition, one-pot and point-of-care tests are under development, and future challenges are discussed. We present in this literature review the potential and actual impact of CRISPR-Cas research on human tuberculosis understanding and management. Altogether, the CRISPR-revolution will revitalize the fight against tuberculosis with more research and technological developments.

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.

Mycobacterium tuberculosis genetic features associated with pulmonary tuberculosis severity.

OBJECTIVES: Mycobacterium tuberculosis (Mtb) infections result in a wide spectrum of clinical presentations but without proven Mtb genetic determinants. Herein, we hypothesized that the genetic features of Mtb clinical isolates, such as specific polymorphisms or microdiversity, may be linked to tuberculosis (TB) severity. METHODS: A total of 234 patients with pulmonary TB (including 193 drug-susceptible and 14 monoresistant cases diagnosed between 2017 and 2020 and 27 multidrug-resistant cases diagnosed between 2010 and 2020) were stratified according to TB disease severity, and Mtb genetic features were explored using whole genome sequencing, including heterologous single-nucleotide polymorphism (SNP), calling to explore microdiversity. Finally, we performed a structural equation modeling analysis to relate TB severity to Mtb genetic features. RESULTS: The clinical isolates from patients with mild TB carried mutations in genes associated with host-pathogen interaction, whereas those from patients with moderate/severe TB carried mutations associated with regulatory mechanisms. Genome-wide association study identified an SNP in the promoter of the gene coding for the virulence regulator espR, statistically associated with moderate/severe disease. Structural equation modeling and model comparisons indicated that TB severity was associated with the detection of Mtb microdiversity within clinical isolates and to the espR SNP. CONCLUSION: Taken together, these results provide a new insight to better understand TB pathophysiology and could provide a new prognosis tool for pulmonary TB severity.

Mycobacterium tuberculosis Diversity Exploration: A Way to Serve the Three Main Weapons against Epidemics, Hygiene, Vaccine Development and Chemotherapy.

As highlighted by the SARS pandemic which is still ongoing, the battle against pathogens relies on three main "weapons": hygiene, vaccine development and chemotherapy strategies [...].

Consistency of Mycobacterium tuberculosis Complex Spoligotyping between the Membrane-Based Method and In Silico Approach.

To tackle the spread of tuberculosis (TB), epidemiological studies are undertaken worldwide to investigate TB transmission chains. Clustered regulatory interspaced short palindromic repeats (CRISPR) locus diversity, also called spoligotyping, is a widely used genotyping assay for the characterization of Mycobacterium tuberculosis complex (MTBC). We compared herein the spoligotyping of MTBC clinical isolates using a membrane-based method (following an initial PCR step) and whole-genome sequencing (WGS)-based method (i.e., in silico spoligotyping). All MTBC strains isolated at the Lyon University Hospital, France, between November 2016 and December 2020 were included (n = 597). Spoligotyping profiles were also used for species identification among the MTBC. Outputs of both methods were analyzed, and discrepant results were investigated thanks to CRISPRbuilder-TB. The overall agreement was 85.7%. Spacer discrepancies observed between the methods were due to the insertion of IS6110 within the direct repeat (DR) sequence upstream or downstream of spacers, mutated DR sequences, or truncated spacers. Discrepancies did not impact species identification. Although spoligotyping-based species identification was inconclusive for 29 isolates, SNP-based phylogeny conducted after WGS allowed the identification of 23 M. tuberculosis (Mtb), 2 M. canettii, and 4 mixed MTBC infections. WGS yielded very few discrepancies compared to membrane-based spoligotyping. Overall agreement was significantly improved (92.4%) by the CRISPR locus reconstruction using CRISPRbuilder-TB for the MTBC isolates with the shared international type 53 in silico spoligotyping. A smooth transition from the membrane-based to the in silico-based genotyping of M. tuberculosis isolates is, therefore, possible for TB diagnosis and epidemiologic survey. IMPORTANCE Whole-genome sequencing (WGS) has profoundly transformed the perspectives of tuberculosis (TB) diagnosis, providing a better discriminatory power to determine relatedness between Mycobacterium tuberculosis complex (MTBC) isolates. Previous genotyping approaches, such as spoligotyping consisting of an initial PCR step followed by reverse dot hybridization, are currently being replaced by WGS. Several pipelines have been developed to extract a spoligotype from WGS data (in silico spoligotyping) allowing for the continuity of MTBC molecular surveys before and after WGS implementation. The present study found very good overall agreement between hybridization to membrane-based spoligotyping and in silico spoligotyping, indicating the possibility of a smooth transition from the traditional to the in silico-based genotyping of MTBC isolates for TB diagnosis and epidemiological survey.

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.

Molecular epidemiology of Mycobacterium tuberculosis in Brazil before the whole genome sequencing era: a literature review.

Molecular-typing can help in unraveling epidemiological scenarios and improvement for disease control strategies. A literature review of Mycobacterium tuberculosis transmission in Brazil through genotyping on 56 studies published from 1996-2019 was performed. The clustering rate for mycobacterial interspersed repetitive units - variable tandem repeats (MIRU-VNTR) of 1,613 isolates were: 73%, 33% and 28% based on 12, 15 and 24-loci, respectively; while for RFLP-IS6110 were: 84% among prison population in Rio de Janeiro, 69% among multidrug-resistant isolates in Rio Grande do Sul, and 56.2% in general population in São Paulo. These findings could improve tuberculosis (TB) surveillance and set up a solid basis to build a database of Mycobacterium genomes.

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.

Spoligotyping of Mycobacterium tuberculosis isolates using Luminex®-based method in Lebanon.

INTRODUCTION: Data about the genotypes of circulating Mycobacterium tuberculosis isolates (MTB) in Lebanon are scarce. This study was undertaken to reveal the spoligotypes of MTB isolates recovered from patients in Lebanon. METHODOLOGY: MTB isolates from 49 patients living in Lebanon were recovered and identified. The samples were heat killed and subjected to DNA extraction. Spoligotyping was performed using microbeads from TB-SPOL Kit and the fluorescence intensity was measured using Luminex 200®. Generated patterns were assigned to families using the SITVIT2 international database of the Pasteur Institute of Guadeloupe and compared. RESULTS: The spoligotyping of the 49 MTB isolates revealed that 31 isolates belonged to Lineage 4 (Euro-American, 63.3%), 12 to Lineage 3 (East- African Indian, 24.5%), 3 to Lineage 2 (East Asian, 6%) and 2 were unknown. Over half of the genotypes (16 of 30) harbored SIT127 supposed to belong to the L4.5 sublineage. One isolate belonging to the rare Manu-Ancestor SIT523 was recovered for the first time in Lebanon, being associated with highly virulent extensively drug-resistant (XDR) MTB phenotype. CONCLUSION: The application of the Spoligotyping Multiplex Luminex® method is an efficient, discriminatory and rapid method to use for first-lane genotyping of MTB isolates. Though humble numbers were tested, this study is one of the first to describe the genomic diversity and epidemiology of MTB isolates of Lebanon, and suggests an increasing prevalence of SIT127 in the country.

Prediction of the hidden genotype of mixed infection strains in Iranian tuberculosis patients.

BACKGROUND: Patients with mixed-strain Mycobacterium tuberculosis infections may be at a high risk of poor treatment outcomes. However, the mechanisms through which mixed infections affect the clinical manifestations are not well recognized. Evidence suggests that failure to detect the pathogen diversity within the host can influence the clinical results. We aimed to investigate the effects of different genotypes in mixed infections and determine their relationship with heteroresistance in the treatment of Iranian tuberculosis patients. METHODS: One of the genotypes was identified in the culture and another genotype pattern in the mixed infection was predicted by comparing the pattern of MIRU-VNTR between the clinical specimens and their respective cultures in each patient. For all patients, the drug susceptibility testing was carried out on three single colonies from each clinical sample. The follow-up of patients was carried out during six months of treatment. RESULTS: Based on MIRU-VNTR profiles of clinical samples, we showed that 55.6% (25/45) of the Iranian patients included in the study had mixed infections. Patients with mixed infections had a higher rate of treatment failure, compared to others (P=0.03). By comparing clinical sample profiles to profiles obtained after culture, we were able to distinguish between major and hidden genotypes. Among hidden genotypes, Haarlem (L4.1.2) and Beijing (L2) were associated to treatment failure (6/8 patients). CONCLUSIONS: To conclude, we propose a procedure using the MIRU-VNTR method to identify the different genotypes in mixed infections. The present findings suggest that genotypes with potentially higher pathogenicity may not be detected when performing experimental culture in patients with mixed infections.

NTF-RINT, a new method for the epidemiological surveillance of MDR Mycobacterium tuberculosis L2/Beijing strains.

The most widely discussed antibiotic-resistant tuberculosis strains ("W" and "B0/W148", "CAO") belong to L2/Beijing Lineage and are characterized by IS6110 insertion sequences at the NTF locus. We present a high-throughput, microbead-based method, called NTF-RINT for detection of IS in NTF and Rifampicin and Isoniazid Typing. This method provides tuberculosis diagnostic confirmation, screens for the so-called modern L2/Beijing sublineage and detects mutations involved in resistance to Rifampicin (RIF) and Isoniazid (INH).

Diversity of Mycobacterium tuberculosis in the Middle Fly District of Western Province, Papua New Guinea: microbead-based spoligotyping using DNA from Ziehl-Neelsen-stained microscopy preparations.

Tuberculosis remains the world's leading cause of death from an infectious agent, and is a serious health problem in Papua New Guinea (PNG) with an estimated 36,000 new cases each year. This study describes the genetic diversity of Mycobacterium tuberculosis among tuberculosis patients in the Balimo/Bamu region in the Middle Fly District of Western Province in PNG, and investigates rifampicin resistance-associated mutations. Archived Ziehl-Neelsen-stained sputum smears were used to conduct microbead-based spoligotyping and assess genotypic resistance. Among the 162 samples included, 80 (49.4%) generated spoligotyping patterns (n = 23), belonging predominantly to the L2 Lineage (44%) and the L4 Lineage (30%). This is consistent with what has been found in other PNG regions geographically distant from Middle Fly District of Western Province, but is different from neighbouring South-East Asian countries. Rifampicin resistance was identified in 7.8% of the successfully sequenced samples, with all resistant samples belonging to the L2/Beijing Lineage. A high prevalence of mixed L2/L4 profiles was suggestive of polyclonal infection in the region, although this would need to be confirmed. The method described here could be a game-changer in resource-limited countries where large numbers of archived smear slides could be used for retrospective (and prospective) studies of M. tuberculosis genetic epidemiology.

Mycobacterium tuberculosis lineage 1 genetic diversity in Pará, Brazil, suggests common ancestry with east-African isolates potentially linked to historical slave trade.

Lineage 1 (L1) is one of seven Mycobacterium tuberculosis complex (MTBC) lineages. The objective of this study was to improve the complex taxonomy of L1 using phylogenetic SNPs, and to look for the origin of the main L1 sublineage prevalent in Para, Brazil. We developed a high-throughput SNPs-typing assay based on 12-L1-specific SNPs. This assay allowed us to experimentally retrieve SNP patterns on nine of these twelve SNPs in 277 isolates previously tentatively assigned to L1 spoligotyping-based sublineages. Three collections were used: Pará-Brazil (71); RIVM, the Netherlands (102), Madagascar (104). One-hundred more results were generated in Silico using the PolyTB database. Based on the final SNPs combination, the samples were classified into 11 clusters (C1-C11). Most isolates within a SNP-based cluster shared a mutual spoligotyping-defined lineage. However, L1/EAI1-SOM (SIT48) and L1/EAI6-BGD1 (SIT591) showed a poor correlation with SNP data and are not monophyletic. L1/EAI8-MDG and L1/EAI3-IND belonged to C5; this result suggests that they share a common ancestor. L1.1.3/SIT129, a spoligotype pattern found in SNPs-cluster C6, was found to be shared between Pará/Brazil and Malawi. SIT129 was independently found to be highly prevalent in Mozambique, which suggests a migration history from East-Africa to Brazil during the 16th-18th slave trade period to Northern Brazil.

TB-EFI, a novel 18-Plex microbead-based method for prediction of second-line drugs and ethambutol resistance in Mycobacterium tuberculosis complex.

Several diagnostic tests are being developed to detect drug resistance in tuberculosis. In line with previous developments detecting rifampicin and isoniazid resistance using microbead-based systems (spoligoriftyping and TB-SPRINT), we present here an assay called TB-EFI detecting mutations involved in resistance to ethambutol, fluoroquinolones and the three classical injectable drugs (kanamycin, amikacin and capreomycin) in Mycobacterium tuberculosis. The proposed test includes both wild-type and mutant probes for each targeted locus. Basic analysis can be performed manually. An upgraded interpretation is made available in Excel 2016®. Using a reference set of 61 DNA extracts, we show that TB-EFI provides perfect concordance with pyrosequencing. Concordance between genotypic resistance and phenotypic DST was relatively good (72 to 98% concordance), with lower efficiency for fluoroquinolones and ethambutol due to some untargeted mutations. When compared to phenotypical resistance, performances were similar to those obtained with Hain MTBDRsl assay, possibly thanks to the use of automatized processing of data although some mutations involved in fluoroquinolone resistance could not be included. When applied on three uncharacterized sets, phenotype could be predicted for 51% to 98% depending on the setting and the drug investigated, detecting one extensively drug-resistant isolate in each of a Pakistan and a Brazilian set of 91 samples, and 9 XDR among 43 multi-resistant Kazakhstan samples. By allowing high-throughput detection of second-line drugs resistance and of resistance to ethambutol that is often combined to second-line treatments, TB-EFI is a cost-effective assay for large-scale worldwide surveillance of resistant tuberculosis and XDR-TB control.

Reverse line probe assay for cheap detection of Single Nucleotide Polymorphisms in Mycobacterium tuberculosis.

More and more Single Nucleotide Polymosrphisms of interest among pathogenic organisms are described with the advent of Whole Genome Sequencing but WGS approach is still too expensive, time consuming, and relying on bioinformatical means that are not available in many developing countries. This study presents a low-cost reverse hybridization line probe technique for detecting SNPs in Mycobacterium tuberculosis. The proposed test is able to detect mutations in the RRDR of rpoB gene in M. tuberculosis with specificity and sensitivity of 98% and 100%, respectively and for an average cost of less than €3 per sample. The technique proved efficient not only on pure DNA samples extracted from culture isolates but also on crude extracts from clinical samples. The flexibility of the platform allows to get it transformed to any kind of test detection, hence, building a bridge between rich countries performing SNP discovery and countries with high burden that can target these SNPs on the collected samples.