Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis.

Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen.

The Genetic Background Modulates the Evolution of Fluoroquinolone-Resistance in Mycobacterium tuberculosis.

Fluoroquinolones (FQ) form the backbone in experimental treatment regimens against drug-susceptible tuberculosis. However, little is known on whether the genetic variation present in natural populations of Mycobacterium tuberculosis (Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used nine genetically distinct drug-susceptible clinical isolates of Mtb and measured their frequency of resistance to the FQ ofloxacin (OFX) in vitro. We found that the Mtb genetic background led to differences in the frequency of OFX-resistance (OFX-R) that spanned two orders of magnitude and substantially modulated the observed mutational profiles for OFX-R. Further, in vitro assays showed that the genetic background also influenced the minimum inhibitory concentration and the fitness effect conferred by a given OFX-R mutation. To test the clinical relevance of our in vitro work, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinical Mtb isolates, and found substantial Mtb lineage-dependent variability. Comparison of the clinical and the in vitro mutational profiles for FQ-R showed that 51% and 39% of the variability in the clinical frequency of FQ-R gyrA mutation events in Lineage 2 and Lineage 4 strains, respectively, can be attributed to how Mtb evolves FQ-R in vitro. As the Mtb genetic background strongly influenced the evolution of FQ-R in vitro, we conclude that the genetic background of Mtb also impacts the evolution of FQ-R in the clinic.

HIV Coinfection Is Associated with Low-Fitness rpoB Variants in Rifampicin-Resistant Mycobacterium tuberculosis.

We analyzed 312 drug-resistant genomes of Mycobacterium tuberculosis isolates collected from HIV-coinfected and HIV-negative TB patients from nine countries with a high tuberculosis burden. We found that rifampicin-resistant M. tuberculosis strains isolated from HIV-coinfected patients carried disproportionally more resistance-conferring mutations in rpoB that are associated with a low fitness in the absence of the drug, suggesting these low-fitness rpoB variants can thrive in the context of reduced host immunity.

Natural Polymorphisms in Mycobacterium tuberculosis Conferring Resistance to Delamanid in Drug-Naive Patients.

Mutations in the genes of the F420 signaling pathway of Mycobacterium tuberculosis complex, including dnn, fgd1, fbiA, fbiB, fbiC, and fbiD, can lead to delamanid resistance. We searched for such mutations among 129 M. tuberculosis strains from Asia, South America, and Africa using whole-genome sequencing; 70 (54%) strains had at least one mutation in one of the genes. For 10 strains with mutations, we determined the MIC of delamanid. We found one strain from a delamanid-naive patient carrying the natural polymorphism Tyr29del (ddn) that was associated with a critical delamanid MIC.

Whole-Genome Sequencing for Drug Resistance Profile Prediction in Mycobacterium tuberculosis.

Whole-genome sequencing allows rapid detection of drug-resistant Mycobacterium tuberculosis isolates. However, the availability of high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data have thus far been limited. We determined drug resistance profiles of 176 genetically diverse clinical M. tuberculosis isolates from the Democratic Republic of the Congo, Ivory Coast, Peru, Thailand, and Switzerland by quantitative phenotypic DST for 11 antituberculous drugs using the BD Bactec MGIT 960 system and 7H10 agar dilution to generate a cross-validated phenotypic DST readout. We compared DST results with predicted drug resistance profiles inferred by whole-genome sequencing. Classification of strains by the two phenotypic DST methods into resistotype/wild-type populations was concordant in 73 to 99% of cases, depending on the drug. Our data suggest that the established critical concentration (5 mg/liter) for ethambutol resistance (MGIT 960 system) is too high and misclassifies strains as susceptible, unlike 7H10 agar dilution. Increased minimal inhibitory concentrations were explained by mutations identified by whole-genome sequencing. Using whole-genome sequences, we were able to predict quantitative drug resistance levels for the majority of drug resistance mutations. Predicting quantitative levels of drug resistance by whole-genome sequencing was partially limited due to incompletely understood drug resistance mechanisms. The overall sensitivity and specificity of whole-genome-based DST were 86.8% and 94.5%, respectively. Despite some limitations, whole-genome sequencing has the potential to infer resistance profiles without the need for time-consuming phenotypic methods.

Revised Interpretation of the Hain Lifescience GenoType MTBC To Differentiate Mycobacterium canettii and Members of the Mycobacterium tuberculosis Complex.

Using 894 phylogenetically diverse genomes of the Mycobacterium tuberculosis complex (MTBC), we simulated in silico the ability of the Hain Lifescience GenoType MTBC assay to differentiate the causative agents of tuberculosis. Here, we propose a revised interpretation of this assay to reflect its strengths (e.g., it can distinguish some strains of Mycobacterium canettii and variants of Mycobacterium bovis that are not intrinsically resistant to pyrazinamide) and limitations (e.g., Mycobacterium orygis cannot be differentiated from Mycobacterium africanum).

Treemmer: a tool to reduce large phylogenetic datasets with minimal loss of diversity.

BACKGROUND: Large sequence datasets are difficult to visualize and handle. Additionally, they often do not represent a random subset of the natural diversity, but the result of uncoordinated and convenience sampling. Consequently, they can suffer from redundancy and sampling biases. RESULTS: Here we present Treemmer, a simple tool to evaluate the redundancy of phylogenetic trees and reduce their complexity by eliminating leaves that contribute the least to the tree diversity. CONCLUSIONS: Treemmer can reduce the size of datasets with different phylogenetic structures and levels of redundancy while maintaining a sub-sample that is representative of the original diversity. Additionally, it is possible to fine-tune the behavior of Treemmer including any kind of meta-information, making Treemmer particularly useful for empirical studies.

Transmission of Drug-Resistant Leprosy in Guinea-Conakry Detected Using Molecular Epidemiological Approaches.

Molecular drug susceptibility testing was performed on skin biopsies from 24 leprosy patients from Guinea-Conakry for the first time. We identified primary drug resistance in 4 cases and a dapsone-resistant cluster caused by the same strain. Primary transmission of drug-resistant Mycobacterium leprae, including a rifampicin-resistant strain, is reported.

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.

Prisons as ecological drivers of fitness-compensated multidrug-resistant Mycobacterium tuberculosis.

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance1. Drug resistance-conferring mutations frequently cause fitness costs in bacteria2-5. Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations6-10. However, the clinical relevance of compensatory evolution remains poorly understood. Here we show that, in the country of Georgia, during a 6-year nationwide study, 63% of MDR-TB was due to patient-to-patient transmission. Compensatory mutations and patient incarceration were independently associated with transmission. Furthermore, compensatory mutations were overrepresented among isolates from incarcerated individuals that also frequently spilled over into the non-incarcerated population. As a result, up to 31% of MDR-TB in Georgia was directly or indirectly linked to prisons. We conclude that prisons fuel the epidemic of MDR-TB in Georgia by acting as ecological drivers of fitness-compensated strains with high transmission potential.

Genomic epidemiological analysis identifies high relapse among individuals with recurring tuberculosis and provides evidence of recent household-related transmission of tuberculosis in Ghana.

OBJECTIVE: To retrospectively investigate the cause of recurring tuberculosis (rcTB) among participants with pulmonary TB recruited from a prospective population-based study conducted between July 2012 and December 2015. METHODS: Mycobacterium tuberculosis complex isolates obtained from rcTB cases were characterized by standard mycobacterial genotyping tools, whole-genome sequencing, and phylogenetic analysis carried out to assess strain relatedness. RESULTS: The majority (58.3%, 21/36) of study participants with rcTB episodes had TB recurrence within 12 months post treatment. TB strains with isoniazid (INH) resistance were found in 19.4% (7/36) of participants at the primary episode, of which 29% (2/7) were also rifampicin-resistant. On TB recurrence, an INH-resistant strain was found in a larger proportion of participants, 27.8% (10/36), of which 40% (4/10) were MDR-TB strains. rcTB was attributed to relapse (same strain) in 75.0% (27/36) of participants and 25.0% (9/36) to re-infection. CONCLUSION: Our findings indicate that previous unresolved infectiondue to inadequate treatment, may be the major cause of rcTB.

Mortality from drug-resistant tuberculosis in high-burden countries comparing routine drug susceptibility testing with whole-genome sequencing: a multicentre cohort study.

BACKGROUND: Drug resistance threatens global tuberculosis control. We aimed to examine mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and whole-genome sequencing (WGS). METHODS: In this multicentre cohort study, we collected pulmonary Mycobacterium tuberculosis isolates and clinical data from individuals with tuberculosis from antiretroviral therapy programmes and tuberculosis clinics in Côte d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, Peru, South Africa, and Thailand, stratified by HIV status and drug resistance. Sites tested drug susceptibility using routinely available methods. WGS was done on Illumina HiSeq 2500 in the USA and Switzerland, and TBprofiler was used to analyse the genomes. We included individuals aged 16 years or older with pulmonary tuberculosis (bacteriologically confirmed or clinically diagnosed). We analysed mortality in multivariable logistic regression models adjusted for sex, age, HIV status, history of tuberculosis, and sputum positivity. FINDINGS: Between Sept 1, 2014, and July 4, 2016, of 634 patients included in our previous analysis, we included 582 patients with tuberculosis (median age 33 years [IQR 27-43], 225 [39%] women, and 247 [42%] HIV-positive). Based on WGS, 339 (58%) isolates were pan-susceptible, 35 (6%) monoresistant, 146 (25%) multidrug-resistant, and 24 (4%) pre-extensively drug-resistant (pre-XDR) or XDR. The analysis of mortality was based on 530 patients; 63 (12%) died and 77 (15%) patients received inappropriate treatment. Mortality ranged from 6% (18 of 310) in patients with pan-susceptible tuberculosis to 39% (nine of 23) in patients with pre-XDR or XDR tuberculosis. The adjusted odds ratio for mortality was 4·92 (95% CI 2·47-9·78) among undertreated patients, compared with appropriately treated patients. INTERPRETATION: In seven countries with a high burden of tuberculosis, we observed discrepancies between drug resistance patterns obtained locally and WGS. The underdiagnosis of drug resistance resulted in inappropriate treatment and higher mortality. WGS can provide accurate and detailed drug resistance information required to improve the outcomes of drug-resistant tuberculosis in high-burden settings. Our results support WHO's call for point-of-care tests based on WGS. FUNDING: National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, and Swiss National Center for Mycobacteria.

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.

Local adaptation in populations of Mycobacterium tuberculosis endemic to the Indian Ocean Rim.

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world's new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics. For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.

An African origin for Mycobacterium bovis.

BACKGROUND AND OBJECTIVES: Mycobacterium bovis and Mycobacterium caprae are two of the most important agents of tuberculosis in livestock and the most important causes of zoonotic tuberculosis in humans. However, little is known about the global population structure, phylogeography and evolutionary history of these pathogens. METHODOLOGY: We compiled a global collection of 3364 whole-genome sequences from M.bovis and M.caprae originating from 35 countries and inferred their phylogenetic relationships, geographic origins and age. RESULTS: Our results resolved the phylogenetic relationship among the four previously defined clonal complexes of M.bovis, and another eight newly described here. Our phylogeographic analysis showed that M.bovis likely originated in East Africa. While some groups remained restricted to East and West Africa, others have subsequently dispersed to different parts of the world. CONCLUSIONS AND IMPLICATIONS: Our results allow a better understanding of the global population structure of M.bovis and its evolutionary history. This knowledge can be used to define better molecular markers for epidemiological investigations of M.bovis in settings where whole-genome sequencing cannot easily be implemented. LAY SUMMARY: During the last few years, analyses of large globally representative collections of whole-genome sequences (WGS) from the human-adapted Mycobacterium tuberculosis complex (MTBC) lineages have enhanced our understanding of the global population structure, phylogeography and evolutionary history of these pathogens. In contrast, little corresponding data exists for M. bovis, the most important agent of tuberculosis in livestock. Using whole-genome sequences of globally distributed M. bovis isolates, we inferred the genetic relationships among different M. bovis genotypes distributed around the world. The most likely origin of M. bovis is East Africa according to our inferences. While some M. bovis groups remained restricted to East and West Africa, others have subsequently dispersed to different parts of the world driven by cattle movements.

Whole Genome Sequencing and Spatial Analysis Identifies Recent Tuberculosis Transmission Hotspots in Ghana.

Whole genome sequencing (WGS) is progressively being used to investigate the transmission dynamics of Mycobacterium tuberculosis complex (MTBC). We used WGS analysis to resolve traditional genotype clusters and explored the spatial distribution of confirmed recent transmission clusters. Bacterial genomes from a total of 452 MTBC isolates belonging to large traditional clusters from a population-based study spanning July 2012 and December 2015 were obtained through short read next-generation sequencing using the illumina HiSeq2500 platform. We performed clustering and spatial analysis using specified R packages and ArcGIS. Of the 452 traditional genotype clustered genomes, 314 (69.5%) were confirmed clusters with a median cluster size of 7.5 genomes and an interquartile range of 4-12. Recent tuberculosis (TB) transmission was estimated as 24.7%. We confirmed the wide spread of a Cameroon sub-lineage clone with a cluster size of 78 genomes predominantly from the Ablekuma sub-district of Accra metropolis. More importantly, we identified a recent transmission cluster associated with isoniazid resistance belonging to the Ghana sub-lineage of lineage 4. WGS was useful in detecting unsuspected outbreaks; hence, we recommend its use not only as a research tool but as a surveillance tool to aid in providing the necessary guided steps to track, monitor, and control TB.

A sister lineage of the Mycobacterium tuberculosis complex discovered in the African Great Lakes region.

The human- and animal-adapted lineages of the Mycobacterium tuberculosis complex (MTBC) are thought to have expanded from a common progenitor in Africa. However, the molecular events that accompanied this emergence remain largely unknown. Here, we describe two MTBC strains isolated from patients with multidrug resistant tuberculosis, representing an as-yet-unknown lineage, named Lineage 8 (L8), seemingly restricted to the African Great Lakes region. Using genome-based phylogenetic reconstruction, we show that L8 is a sister clade to the known MTBC lineages. Comparison with other complete mycobacterial genomes indicate that the divergence of L8 preceded the loss of the cobF genome region - involved in the cobalamin/vitamin B12 synthesis - and gene interruptions in a subsequent common ancestor shared by all other known MTBC lineages. This discovery further supports an East African origin for the MTBC and provides additional molecular clues on the ancestral genome reduction associated with adaptation to a pathogenic lifestyle.

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

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

Whole genome sequencing of Mycobacterium tuberculosis: current standards and open issues.

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has rapidly progressed from a research tool to a clinical application for the diagnosis and management of tuberculosis and in public health surveillance. This development has been facilitated by drastic drops in cost, advances in technology and concerted efforts to translate sequencing data into actionable information. There is, however, a risk that, in the absence of a consensus and international standards, the widespread use of WGS technology may result in data and processes that lack harmonization, comparability and validation. In this Review, we outline the current landscape of WGS pipelines and applications, and set out best practices for M. tuberculosis WGS, including standards for bioinformatics pipelines, curated repositories of resistance-causing variants, phylogenetic analyses, quality control and standardized reporting.

A New Phylogenetic Framework for the Animal-Adapted Mycobacterium tuberculosis Complex.

Tuberculosis (TB) affects humans and other animals and is caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Previous studies have shown that there are at least nine members of the MTBC infecting animals other than humans; these have also been referred to as ecotypes. However, the ecology and the evolution of these animal-adapted MTBC ecotypes are poorly understood. Here we screened 12,886 publicly available MTBC genomes and newly sequenced 17 animal-adapted MTBC strains, gathering a total of 529 genomes of animal-adapted MTBC strains. Phylogenomic and comparative analyses confirm that the animal-adapted MTBC members are paraphyletic with some members more closely related to the human-adapted Mycobacterium africanum Lineage 6 than to other animal-adapted strains. Furthermore, we identified four main animal-adapted MTBC clades that might correspond to four main host shifts; two of these clades are hypothesized to reflect independent cattle domestication events. Contrary to what would be expected from an obligate pathogen, MTBC nucleotide diversity was not positively correlated with host phylogenetic distances, suggesting that host tropism in the animal-adapted MTBC seems to be driven by contact rates and demographic aspects of the host population rather by than host relatedness. By combining phylogenomics with ecological data, we propose an evolutionary scenario in which the ancestor of Lineage 6 and all animal-adapted MTBC ecotypes was a generalist pathogen that subsequently adapted to different host species. This study provides a new phylogenetic framework to better understand the evolution of the different ecotypes of the MTBC and guide future work aimed at elucidating the molecular mechanisms underlying host range.