The effect of antibiotic selection on collateral effects and evolvability of uropathogenic Escherichia coli.

Trimethoprim is recommended as a first-line treatment of urinary tract infections (UTIs) in the UK. In 2018, 31.4% of Escherichia coli isolated from UTIs in England were trimethoprim-resistant, leading to overreliance on other first and second-line antibiotics. Here, we assessed whether, in principle, prior selection with trimethoprim results in collateral effects to other antibiotics recommended for the treatment of UTIs. As collateral effects, we considered changes in susceptibility, mutation-selection window and population establishment probability. We selected 10 trimethoprim-resistant derivatives from three clinical isolates of uropathogenic Escherichia coli. We found that mutations conferring trimethoprim resistance did not have any collateral effects on fosfomycin. In contrast, resistance to trimethoprim resulted in decreased susceptibility (collateral resistance) to nitrofurantoin, below the clinical breakpoint and narrowed the mutation-selection window, thereby reducing the maximum concentration for selection of nitrofurantoin resistance mutations. Our analyses demonstrate that multiple collateral responses should be accounted for when predicting and optimising antibiotic use, limiting future antimicrobial resistance emergence.

Phylogenomic and genomic analysis reveals unique and shared genetic signatures of Mycobacterium kansasii complex species.

Species belonging to the Mycobacterium kansasii complex (MKC) are frequently isolated from humans and the environment and can cause serious diseases. The most common MKC infections are caused by the species M. kansasii (sensu stricto), leading to tuberculosis-like disease. However, a broad spectrum of virulence, antimicrobial resistance and pathogenicity of these non-tuberculous mycobacteria (NTM) are observed across the MKC. Many genomic aspects of the MKC that relate to these broad phenotypes are not well elucidated. Here, we performed genomic analyses from a collection of 665 MKC strains, isolated from environmental, animal and human sources. We inferred the MKC pangenome, mobilome, resistome, virulome and defence systems and show that the MKC species harbours unique and shared genomic signatures. High frequency of presence of prophages and different types of defence systems were observed. We found that the M. kansasii species splits into four lineages, of which three are lowly represented and mainly in Brazil, while one lineage is dominant and globally spread. Moreover, we show that four sub-lineages of this most distributed M. kansasii lineage emerged during the twentieth century. Further analysis of the M. kansasii genomes revealed almost 300 regions of difference contributing to genomic diversity, as well as fixed mutations that may explain the M. kansasii's increased virulence and drug resistance.

Single nucleotide variation catalog from clinical isolates mapped on tertiary and quaternary structures of ESX-1-related proteins reveals critical regions as putative Mtb therapeutic targets.

Proteins encoded by the ESX-1 genes of interest are essential for full virulence in all Mycobacterium tuberculosis complex (Mtbc) lineages, the pathogens causing the highest mortality worldwide. Identifying critical regions in these ESX-1-related proteins could provide preventive or therapeutic targets for Mtb infection, the game changer needed for tuberculosis control. We analyzed a compendium of whole genome sequences of clinical Mtb isolates from all lineages from >32,000 patients and identified single nucleotide polymorphisms. When mutations corresponding to all non-synonymous single nucleotide polymorphisms were mapped on structural models of the ESX-1 proteins, fully conserved regions emerged. Some could be assigned to known quaternary structures, whereas others could be predicted to be involved in yet-to-be-discovered interactions. Some mutants had clonally expanded (found in >1% of the isolates); these mutants were mostly located at the surface of globular domains, remote from known intra- and inter-molecular protein-protein interactions. Fully conserved intrinsically disordered regions of proteins were found, suggesting that these regions are crucial for the pathogenicity of the Mtbc. Altogether, our findings highlight fully conserved regions of proteins as attractive vaccine antigens and drug targets to control Mtb virulence. Extending this approach to the whole Mtb genome as well as other microorganisms will enhance vaccine development for various pathogens. IMPORTANCE: We mapped all non-synonymous single nucleotide polymorphisms onto each of the experimental and predicted ESX-1 proteins' structural models and inspected their placement. Varying sizes of conserved regions were found. Next, we analyzed predicted intrinsically disordered regions within our set of proteins, finding two putative long stretches that are fully conserved, and discussed their potential essential role in immunological recognition. Combined, our findings highlight new targets for interfering with Mycobacterium tuberculosis complex virulence.

Inferring diet, disease and antibiotic resistance from ancient human oral microbiomes.

The interaction between a host and its microbiome is an area of intense study. For the human host, it is known that the various body-site-associated microbiomes impact heavily on health and disease states. For instance, the oral microbiome is a source of various pathogens and potential antibiotic resistance gene pools. The effect of historical changes to the human host and environment to the associated microbiome, however, has been less well explored. In this review, we characterize several historical and prehistoric events which are considered to have impacted the oral environment and therefore the bacterial communities residing within it. The link between evolutionary changes to the oral microbiota and the significant societal and behavioural changes occurring during the pre-Neolithic, Agricultural Revolution, Industrial Revolution and Antibiotic Era is outlined. While previous studies suggest the functional profile of these communities may have shifted over the centuries, there is currently a gap in knowledge that needs to be filled. Biomolecular archaeological evidence of innate antimicrobial resistance within the oral microbiome shows an increase in the abundance of antimicrobial resistance genes since the advent and widespread use of antibiotics in the modern era. Nevertheless, a lack of research into the prevalence and evolution of antimicrobial resistance within the oral microbiome throughout history hinders our ability to combat antimicrobial resistance in the modern era.

Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches.

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

Genomic transmission clusters and circulating lineages of Mycobacterium tuberculosis among refugees residing in refugee camps in Ethiopia.

BACKGROUND: Understanding the transmission dynamics of Mycobacterium tuberculosis (Mtb) could benefit the design of tuberculosis (TB) prevention and control strategies for refugee populations. Whole Genome Sequencing (WGS) has not yet been used to document the Mtb transmission dynamics among refugees in Ethiopia. We applied WGS to accurately identify transmission clusters and Mtb lineages among TB cases in refugee camps in Ethiopia. METHOD AND DESIGN: We conducted a cross-sectional study of 610 refugees in refugee camps in Ethiopia presenting with symptoms of TB. WGS data of 67 isolates was analyzed using the Maximum Accessible Genome for Mtb Analysis (MAGMA) pipeline; iTol and FigTree were used to visualize phylogenetic trees, lineages, and the presence of transmission clusters. RESULTS: Mtb culture-positive refugees originated from South Sudan (52/67, 77.6%), Somalia (9/67, 13.4%). Eritrea (4/67, 6%), and Sudan (2/67, 3%). The majority (52, 77.6%) of the isolates belonged to Mtb lineage (L) 3, and one L9 was identified from a Somalian refugee. The vast majority (82%) of the isolates were pan-susceptible Mtb, and none were multi-drug-resistant (MDR)-TB. Based on the 5-single nucleotide polymorphisms cutoff, we identified eight potential transmission clusters containing 23.9% of the isolates. Contact investigation confirmed epidemiological links with either family or social interaction within the refugee camps or with neighboring refugee camps. CONCLUSION: Four lineages (L1, L3, L4, and L9) were identified, with the majority of strains being L3, reflecting the Mtb L3 dominance in South Sudan, where the majority of refugees originated from. Recent transmission among refugees was relatively low (24%), likely due to the short study period. The improved understanding of the Mtb transmission dynamics using WGS in refugee camps could assist in designing effective TB control programs for refugees.

Borderline rpoB mutations transmit at the same rate as common rpoB mutations in a tuberculosis cohort in Bangladesh.

The spread of multidrug-resistant tuberculosis (MDR-TB) is a growing problem in many countries worldwide. Resistance to one of the primary first-line drugs, rifampicin, is caused by mutations in the Mycobacterium tuberculosis rpoB gene. So-called borderline rpoB mutations confer low-level resistance, in contrast to more common rpoB mutations which confer high-level resistance. While some borderline mutations show lower fitness in vitro than common mutations, their in vivo fitness is currently unknown. We used a dataset of 394 whole genome sequenced MDR-TB isolates from Bangladesh, representing around 44 % of notified MDR-TB cases over 6 years, to look at differences in transmission clustering between isolates with borderline rpoB mutations and those with common rpoB mutations. We found a relatively low percentage of transmission clustering in the dataset (34.8 %) but no difference in clustering between different types of rpoB mutations. Compensatory mutations in rpoA, rpoB, and rpoC were associated with higher levels of transmission clustering as were lineages two, three, and four relative to lineage one. Young people as well as patients with high sputum smear positive TB were more likely to be in a transmission cluster. Our findings show that although borderline rpoB mutations have lower in vitro growth potential this does not translate into lower transmission potential or in vivo fitness. Proper detection of these mutations is crucial to ensure they do not go unnoticed and spread MDR-TB within communities.

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.

Transmission, distribution and drug resistance-conferring mutations of extensively drug-resistant tuberculosis in the Western Cape Province, South Africa.

Extensively drug-resistant tuberculosis (XDR-TB), defined as resistance to at least isoniazid (INH), rifampicin (RIF), a fluoroquinolone (FQ) and a second-line injectable drug (SLID), is difficult to treat and poses a major threat to TB control. The transmission dynamics and distribution of XDR Mycobacterium tuberculosis (Mtb) strains have not been thoroughly investigated. Using whole genome sequencing data on 461 XDR-Mtb strains, we aimed to investigate the geographical distribution of XDR-Mtb strains in the Western Cape Province of South Africa over a 10 year period (2006-2017) and assess the association between Mtb sub-lineage, age, gender, geographical patient location and membership or size of XDR-TB clusters. First, we identified transmission clusters by excluding drug resistance-conferring mutations and using the 5 SNP cutoff, followed by merging clusters based on their most recent common ancestor. We then consecutively included variants conferring resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA), SLIDs and FQs in the cluster definition. Cluster sizes were classified as small (2-4 isolates), medium (5-20 isolates), large (21-100 isolates) or very large (>100 isolates) to reflect the success of individual strains. We found that most XDR-TB strains were clustered and that including variants conferring resistance to INH, RIF, EMB, PZA and SLIDs in the cluster definition did not significantly reduce the proportion of clustered isolates (85.5-82.2 %) but increased the number of patients belonging to small clusters (4.3-12.4 %, P=0.56). Inclusion of FQ resistance-conferring variants had the greatest effect, with 11 clustered isolates reclassified as unique while the number of clusters increased from 17 to 37. Lineage 2 strains (lineage 2.2.1 typical Beijing or lineage 2.2.2 atypical Beijing) showed the large clusters which were spread across all health districts of the Western Cape Province. We identified a significant association between residence in the Cape Town metropole and cluster membership (P=0.016) but no association between gender, age and cluster membership or cluster size (P=0.39). Our data suggest that the XDR-TB epidemic in South Africa probably has its origin in the endemic spread of MDR Mtb and pre-XDR Mtb strains followed by acquisition of FQ resistance, with more limited transmission of XDR Mtb strains. This only became apparent with the inclusion of drug resistance-conferring variants in the definition of a cluster. In addition to the prevention of amplification of resistance, rapid diagnosis of MDR, pre-XDR and XDR-TB and timely initiation of appropriate treatment is needed to reduce transmission of difficult-to-treat TB.

Multidrug-resistant tuberculosis control in Rwanda overcomes a successful clone that causes most disease over a quarter century.

SUMMARY BACKGROUND: Multidrug-resistant (MDR) tuberculosis (TB) poses an important challenge in TB management and control. Rifampicin resistance (RR) is a solid surrogate marker of MDR-TB. We investigated the RR-TB clustering rates, bacterial population dynamics to infer transmission dynamics, and the impact of changes to patient management on these dynamics over 27 years in Rwanda. METHODS: We analysed whole genome sequences of a longitudinal collection of nationwide RR-TB isolates. The collection covered three important periods: before programmatic management of MDR-TB (PMDT; 1991-2005), the early PMDT phase (2006-2013), in which rifampicin drug-susceptibility testing (DST) was offered to retreatment patients only, and the consolidated phase (2014-2018), in which all bacteriologically confirmed TB patients had rifampicin DST done mostly via Xpert MTB/RIF assay. We constructed clusters based on a 5 SNP cut-off and resistance conferring SNPs. We used Bayesian modelling for dating and population size estimations, TransPhylo to estimate the number of secondary cases infected by each patient, and multivariable logistic regression to assess predictors of being infected by the dominant clone. RESULTS: Of 308 baseline RR-TB isolates considered for transmission analysis, the clustering analysis grouped 259 (84.1%) isolates into 13 clusters. Within these clusters, a single dominant clone was discovered containing 213 isolates (82.2% of clustered and 69.1% of all RR-TB), which we named the "Rwanda Rifampicin-Resistant clone" (R3clone). R3clone isolates belonged to Ugandan sub-lineage 4.6.1.2 and its rifampicin and isoniazid resistance were conferred by the Ser450Leu mutation in rpoB and Ser315Thr in katG genes, respectively. All R3clone isolates had Pro481Thr, a putative compensatory mutation in the rpoC gene that likely restored its fitness. The R3clone was estimated to first arise in 1987 and its population size increased exponentially through the 1990s', reaching maximum size (∼84%) in early 2000 s', with a declining trend since 2014. Indeed, the highest proportion of R3clone (129/157; 82·2%, 95%CI: 75·3-87·8%) occurred between 2000 and 13, declining to 64·4% (95%CI: 55·1-73·0%) from 2014 onward. We showed that patients with R3clone detected after an unsuccessful category 2 treatment were more likely to generate secondary cases than patients with R3clone detected after an unsuccessful category 1 treatment regimen. CONCLUSIONS: RR-TB in Rwanda is largely transmitted. Xpert MTB/RIF assay as first diagnostic test avoids unnecessary rounds of rifampicin-based TB treatment, thus preventing ongoing transmission of the dominant R3clone. As PMDT was intensified and all TB patients accessed rifampicin-resistance testing, the nationwide R3clone burden declined. To our knowledge, our findings provide the first evidence supporting the impact of universal DST on the transmission of RR-TB.

Nanopore Sequencing for Mycobacterium tuberculosis: a Critical Review of the Literature, New Developments, and Future Opportunities.

The next-generation, short-read sequencing technologies that generate comprehensive, whole-genome data with single nucleotide resolution have already advanced tuberculosis diagnosis, treatment, surveillance, and source investigation. Their high costs, tedious and lengthy processes, and large equipment remain major hurdles for research use in high tuberculosis burden countries and implementation into routine care. The portable next-generation sequencing devices developed by Oxford Nanopore Technologies (ONT) are attractive alternatives due to their long-read sequence capability, compact low-cost hardware, and continued improvements in accuracy and throughput. A systematic review of the published literature demonstrated limited uptake of ONT sequencing in tuberculosis research and clinical care. Of the 12 eligible articles presenting ONT sequencing data on at least one Mycobacterium tuberculosis sample, four addressed software development for long-read ONT sequencing data with potential applications for M. tuberculosis. Only eight studies presented results of ONT sequencing of M. tuberculosis, of which five performed whole-genome and three did targeted sequencing. Based on these findings, we summarize the standard processes, reflect on the current limitations of ONT sequencing technology, and the research needed to overcome the main hurdles. The low capital cost, portable nature and continued improvement in the performance of ONT sequencing make it an attractive option for sequencing for research and clinical care, but limited data are available on its application in the tuberculosis field. Important research investment is needed to unleash the full potential of ONT sequencing for tuberculosis research and care.

Reconstituting the genus Mycobacterium.

The definition of a genus has wide-ranging implications both in terms of binomial species names and also evolutionary relationships. In recent years, the definition of the genus Mycobacterium has been debated due to the proposed split of this genus into five new genera (Mycolicibacterium, Mycolicibacter, Mycolicibacillus, Mycobacteroides and an emended Mycobacterium). Since this group of species contains many important obligate and opportunistic pathogens, it is important that any renaming of species does not cause confusion in clinical treatment as outlined by the nomen periculosum rule (56a) of the Prokaryotic Code. In this study, we evaluated the proposed and original genus boundaries for the mycobacteria, to determine if the split into five genera was warranted. By combining multiple approaches for defining genus boundaries (16S rRNA gene similarity, amino acid identity index, average nucleotide identity, alignment fraction and percentage of conserved proteins) we show that the original genus Mycobacterium is strongly supported over the proposed five-way split. Thus, we propose that the original genus label be reapplied to all species within this group, with the proposed five genera potentially used as sub-genus complex names.

Insights into the ancestry evolution of the Mycobacterium tuberculosis complex from analysis of Mycobacterium riyadhense.

Current evolutionary scenarios posit the emergence of Mycobacterium tuberculosis from an environmental saprophyte through a cumulative process of genome adaptation. Mycobacterium riyadhense, a related bacillus, is being increasingly isolated from human clinical cases with tuberculosis-like symptoms in various parts of the world. To elucidate the evolutionary relationship between M. riyadhense and other mycobacterial species, including members of the M. tuberculosis complex (MTBC), eight clinical isolates of M. riyadhense were sequenced and analyzed. We show, among other features, that M. riyadhense shares a large number of conserved orthologs with M. tuberculosis and shows the expansion of toxin/antitoxin pairs, PE/PPE family proteins compared with other non-tuberculous mycobacteria. We observed M. riyadhense lacks wecE gene which may result in the absence of lipooligosaccharides (LOS) IV. Comparative transcriptomic analysis of infected macrophages reveals genes encoding inducers of Type I IFN responses, such as cytosolic DNA sensors, were relatively less expressed by macrophages infected with M. riyadhense or M. kansasii, compared to BCG or M. tuberculosis. Overall, our work sheds new light on the evolution of M. riyadhense, its relationship to the MTBC, and its potential as a system for the study of mycobacterial virulence and pathogenesis.

Quantifying transmission fitness costs of multi-drug resistant tuberculosis.

As multi-drug resistant tuberculosis (MDR-TB) continues to spread, investigating the transmission potential of different drug-resistant strains becomes an ever more pressing topic in public health. While phylogenetic and transmission tree inferences provide valuable insight into possible transmission chains, phylodynamic inference combines evolutionary and epidemiological analyses to estimate the parameters of the underlying epidemiological processes, allowing us to describe the overall dynamics of disease spread in the population. In this study, we introduce an approach to Mycobacterium tuberculosis (M. tuberculosis) phylodynamic analysis employing an existing computationally efficient model to quantify the transmission fitness costs of drug resistance with respect to drug-sensitive strains. To determine the accuracy and precision of our approach, we first perform a simulation study, mimicking the simultaneous spread of drug-sensitive and drug-resistant tuberculosis (TB) strains. We analyse the simulated transmission trees using the phylodynamic multi-type birth-death model (MTBD, (Kühnert et al., 2016)) within the BEAST2 framework and show that this model can estimate the parameters of the epidemic well, despite the simplifying assumptions that MTBD makes compared to the complex TB transmission dynamics used for simulation. We then apply the MTBD model to an M. tuberculosis lineage 4 dataset that primarily consists of MDR sequences. Some of the MDR strains additionally exhibit resistance to pyrazinamide - an important first-line anti-tuberculosis drug. Our results support the previously proposed hypothesis that pyrazinamide resistance confers a transmission fitness cost to the bacterium, which we quantify for the given dataset. Importantly, our sensitivity analyses show that the estimates are robust to different prior distributions on the resistance acquisition rate, but are affected by the size of the dataset - i.e. we estimate a higher fitness cost when using fewer sequences for analysis. Overall, we propose that MTBD can be used to quantify the transmission fitness cost for a wide range of pathogens where the strains can be appropriately divided into two or more categories with distinct properties.

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

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

Identification and Characterization of Mycobacterial Species Using Whole-Genome Sequences.

Whole-genome sequencing (WGS) has shown immense value in enabling identification and characterization of bacterial taxa. This is particularly true for mycobacteria, where culture-based characterization becomes delayed by the inherently slow growth rate of these organisms. This chapter reviews the general techniques behind WGS and their optimization, existing techniques for species-level identification and the advantages of WGS for this purpose, and a variety of useful tools for the genomic characterization of mycobacterial strains.

Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history.

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum, and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.