Search

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 ofMycobacterium tuberculosis(Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used a set ofMtbstrains that included nine genetically distinct drug-susceptible clinical isolates, and measured their frequency of resistance to the FQ ofloxacin (OFX)in vitro. We found that theMtbgenetic 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. Furtherin vitroassays 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 ourin vitrowork, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinicalMtbisolates, and found substantialMtblineage-dependent variability. Comparison of the clinical and thein vitromutational profiles for FQ-R showed that 45% and 19% of the variability in the clinical frequency of FQ-RgyrAmutations in Lineage 2 and Lineage 4 strains, respectively, can be attributed to howMtbevolves FQ-Rin vitro. As theMtbgenetic background strongly influenced the evolution of FQ-Rin vitro, we conclude that the genetic background ofMtbalso impacts the evolution of FQ-R in the clinic.SignificanceNewer generations of fluoroquinolones form the backbone in many experimental treatment regimens againstM. tuberculosis(Mtb). While the genetic variation in natural populations ofMtbcan influence resistance evolution to multiple different antibiotics, it is unclear whether it modulates fluoroquinolone-resistance evolution as well. Using a combination ofin vitroassays coupled with genomic analysis of clinical isolates, we provide the first evidence illustrating theMtbgenetic background’s substantial role in fluoroquinolone-resistance evolution, and highlight the importance of bacterial genetics when studying the prevalence of fluoroquinolone-resistance inMtb. Our work may provide insights into how to maximize the timespan in which fluoroquinolones remain effective in clinical settings, whether as part of current standardized regimens, or in new regimens againstMtb.

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) and presents variable disease outcomes. The variation has primarily been attributed to host and environmental factors, but recent evidence indicates an additional role of genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the potential role of genetic variations in conferring strain-specific adaptive capacity and vulnerability. To systematically identify functionality of single nucleotide polymorphisms (SNPs), we developed a constraint-based approach that integrates metabolomic and genomic data. Model-based predictions were systematically tested against independent metabolome data; they correctly classified SNP effects in pyruvate kinase and suggested a genetic basis for strain-specific sensitivity to the antibiotic para-aminosalicylic acid. Our method is broadly applicable to mutations in enzyme-encoding genes across microbial life, opening new possibilities for identifying strain-specific metabolic vulnerabilities that could lead to more selective treatment strategies.

Whole genome sequencing allows rapid detection of drug-resistantM. tuberculosisisolates. However, high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data have thus far been lacking.We determined drug resistance profiles of 176 genetically diverse clinicalM. tuberculosisisolates from 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 phenotypic drug susceptibility results with predicted drug resistance profiles inferred by whole genome sequencing.Both phenotypic DST methods identically classified the strains into resistant/susceptible in 73-99% of the cases, depending on the drug. Changes in minimal inhibitory concentrations were readily explained by mutations identified by whole genome sequencing. Using the whole genome sequences we were able to predict quantitative drug resistance levels where wild type and mutant MIC distributions did not overlap. The utility of genome sequences to predict quantitative levels of drug resistance was partially limited due to incompletely understood mechanisms influencing the expression of phenotypic drug resistance. The overall sensitivity and specificity of whole genome-based DST were 86.8% and 94.5%, respectively.Despite some limitations, whole genome sequencing has high predictive power to infer resistance profiles without the need for time-consuming phenotypic methods.One sentence summaryWhole genome sequencing of clinicalM. tuberculosisisolates accurately predicts drug resistance profiles and may replace culture-based drug susceptibility testing in the future.

Summary Background Drug resistance is a challenge for the global control of tuberculosis. We examined mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and in a reference laboratory. Methods This multicentre cohort study was done in Cote d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, South Africa, Peru, and Thailand. We collected Mycobacterium tuberculosis isolates and clinical data from adult patients aged 16 years or older. Patients were stratified by HIV status and tuberculosis drug resistance. Molecular or phenotypic drug susceptibility testing was done locally and at the Swiss National Center for Mycobacteria, Zurich, Switzerland. We examined mortality during treatment according to drug susceptibility test results and treatment adequacy in multivariable logistic regression models adjusting for sex, age, sputum microscopy, and HIV status. Findings We obtained M tuberculosis isolates from 871 patients diagnosed between 2013 and 2016. After exclusion of 237 patients, 634 patients with tuberculosis were included in this analysis; the median age was 33·2 years (IQR 26·9–42·5), 239 (38%) were women, 272 (43%) were HIV-positive, and 69 (11%) patients died. Based on the reference laboratory drug susceptibility test, 394 (62%) strains were pan-susceptible, 45 (7%) monoresistant, 163 (26%) multidrug-resistant (MDR), and 30 (5%) had pre-extensively or extensively drug resistant (pre-XDR or XDR) tuberculosis. Results of reference and local laboratories were concordant for 513 (81%) of 634 patients and discordant for 121 (19%) of 634. Overall, sensitivity to detect any resistance was 90·8% (95% CI 86·5–94·2) and specificity 84·3% (80·3–87·7). Mortality ranged from 6% (20 of 336) in patients with pan-susceptible tuberculosis treated according to WHO guidelines to 57% (eight of 14) in patients with resistant strains who were under-treated. In logistic regression models, compared with concordant drug susceptibility test results, the adjusted odds ratio of death was 7·33 (95% CI 2·70–19·95) for patients with discordant results potentially leading to under-treatment. Interpretation Inaccurate drug susceptibility testing by comparison with a reference standard leads to under-treatment of drug-resistant tuberculosis and increased mortality. Rapid molecular drug susceptibility test of first-line and second-line drugs at diagnosis is required to improve outcomes in patients with MDR tuberculosis and pre-XDR or XDR tuberculosis. Funding National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, Swiss National Center for Mycobacteria.

TheMycobacterium tuberculosiscomplex (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 phylogenetic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC reference 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 lab 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 universally effective.

SummaryAntimicrobial resistance poses a threat to global health and the economy. It is widely accepted that, in the absence of antibiotics, drug resistance mutations carry a fitness cost. In the case of rifampicin resistance in fast-growing bacteria, this cost stems from a reduced transcription rate of the RNA polymerase resulting in slower ribosome biosynthesis. However, this relationship does not apply in the slow-growingMycobacterium tuberculosis, where the true mechanism of fitness cost of rifampicin resistance as well as the impact of compensatory evolution remain unknown. Here we show, using global transcriptomic and proteomic profiling of selectedM. tuberculosismutants and clinical strains, that the fitness cost of rifampicin resistance inM. tuberculosisis the result of the physiological burden caused by aberrant gene expression. We further show that the perceived burden can be increased, effectively suppressing the emergence of drug resistance.

Whole-genome sequencing allows rapid detection of drug-resistant

BackgroundHuman tuberculosis (TB) is caused by seven phylogenetic lineages of theMycobacterium tuberculosiscomplex (MTBC), Lineage 1–7. Recent advances in rapid genotyping of MTBC based on single nucleotide polymorphisms (SNP), allow for rapid and phylogenetically robust strain classification, paving the way for defining genotype-phenotype relationships in clinical settings. Such studies have revealed that, in addition to host and environmental factors, different strains of the MTBC influence the outcome of TB infection and disease. In Tanzania, such molecular epidemiological studies of TB however are scarce in spite of a high TB burden.Methods and FindingsHere we used a SNP-typing method to genotype a nationwide collection of 2,039 MTBC clinical isolates obtained from new and retreatment TB cases diagnosed in 2012 and 2013. Four lineages, namely Lineage 1–4 were identified. The distribution and frequency of these lineages varied across the regions but overall, Lineage 4 was the most frequent (n=866, 42.5%), followed by Lineage 3 (n=681, 33.4%) and 1 (n=336, 16.5%), with Lineage 2 being the least frequent (n=92, 4.5%). A total of 64 (3.1%) isolates could not be assigned to any lineage. We found Lineage 2 to be associated with female sex (adjusted odds ratio [aOR] 2.25; 95% confidence interval [95% CI] 1.38 – 3.70, pConclusionsThis study provides novel insights into the influence of pathogen-related factors on the TB epidemic in Tanzania.

Immigrants from regions with a high incidence of tuberculosis (TB) are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000 to 2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-locus mycobacterial interspersed repetitive-unit–variable-number tandem-repeat [MIRU-VNTR] typing and spoligotyping). Here, we used whole-genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR typing clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of the clusters involving only foreign-born patients. The overall clustering proportion was 17% (90 patients; 95% confidence interval [CI], 14 to 21%) by standard genotyping but only 8% (43 patients; 95% CI, 6 to 11%) by WGS. The clustering proportion was 17% (67/401; 95% CI, 13 to 21%) by standard genotyping and 7% (26/401; 95% CI, 4 to 9%) by WGS among foreign-born patients and 19% (23/119; 95% CI, 13 to 28%) and 14% (17/119; 95% CI, 9 to 22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence of an association between birth origin and transmission (adjusted odds ratio of 2.2 and 95% CI of 0.9 to 5.5 comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland compared to WGS, particularly among immigrants from regions with a high TB incidence, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in settings with a low incidence of TB.

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., This work was supported by the Swiss National Science Foundation (Grants 310030_166687, IZRJZ3_164171, IZLSZ3_170834, and CRSII5_177163), the European Research Council (309540-EVODRTB) and SystemsX.ch.

BackgroundDrug resistance and HIV co-infection are challenges for the global control of tuberculosis.MethodsWe collected Mycobacterium tuberculosis isolates from adult patients in Côte d’Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, South Africa, Peru, and Thailand, stratified by HIV status and tuberculosis drug resistance. Molecular or phenotypic drug susceptibility testing (DST) was done locally and at the Swiss tuberculosis reference laboratory. We examined mortality during treatment according to DST results and treatment adequacy in logistic regression models adjusting for sex, age, sputum microscopy and HIV status.Findings634 tuberculosis patients were included; median age was 33.2 years, 239 (37.7%) were female, 272 (42.9%) HIV-positive and 69 (10.9%) patients died. Based on the reference laboratory DST, 394 (62.2%) strains were pan-susceptible, 45 (7.1%) mono-resistant, 163 (25.7%) multidrug-resistant (MDR-TB), and 30 (4.7%) had pre-extensive or extensive drug resistance (pre-XDR/XDR-TB). Results of reference and local laboratories were discordant in 121 (19.1%) cases, corresponding to a sensitivity of 84.3% and a specificity of 90.8%. In patients with drug-resistant tuberculosis, discordant results were associated with increased mortality (risk ratio 1.81; 95% CI 1.07-3.07). In logistic regression, compared to adequately treated patients with pan-susceptible strains, the adjusted odds ratio for death was 4.23 (95% CI 2.16-8.29) for adequately treated patients with drug-resistant strains and 21.54 (95% CI 3.36-138.1) for inadequately treated patients with drug-resistant strains. HIV status was not associated with mortality.InterpretationUsing a reference laboratory standard, inaccurate DST leading to inappropriate treatment of drug-resistant tuberculosis, but not HIV infection, contributed to mortality.

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., BMC Bioinformatics, 19, ISSN:1471-2105

Rapidly spreading antibiotic resistance and the low discovery rate of new antimicrobial compounds demand more effective strategies for early drug discovery. One bottleneck in the drug discovery pipeline is the identification of the modes of action (MoAs) of new compounds. We have developed a rapid systematic metabolome profiling strategy to classify the MoAs of bioactive compounds. The method predicted MoA-specific metabolic responses in the nonpathogenic bacterium Mycobacterium smegmatis after treatment with 62 reference compounds with known MoAs and different metabolic and nonmetabolic targets. We then analyzed a library of 212 new antimycobacterial compounds with unknown MoAs from a drug discovery effort by the pharmaceutical company GlaxoSmithKline (GSK). More than 70% of these new compounds induced metabolic responses in M. smegmatis indicative of known MoAs, seven of which were experimentally validated. Only 8% (16) of the compounds appeared to target unconventional cellular processes, illustrating the difficulty in discovering new antibiotics with different MoAs among compounds used as monotherapies. For six of the GSK compounds with potentially new MoAs, the metabolome profiles suggested their ability to interfere with trehalose and lipid metabolism. This was supported by whole-genome sequencing of spontaneous drug-resistant mutants of the pathogen Mycobacterium tuberculosis and in vitro compound-proteome interaction analysis for one of these compounds. Our compendium of drug-metabolome profiles can be used to rapidly query the MoAs of uncharacterized antimicrobial compounds and should be a useful resource for the drug discovery community.

Background: Drug resistance and HIV co-infection are challenges for the global control of tuberculosis. Methods: We collected Mycobacterium tuberculosis isolates from adult patients in Cote d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, South Africa, Peru, and Thailand, stratified by HIV status and tuberculosis drug resistance. Molecular or phenotypic drug susceptibility testing (DST) was done locally and at the Swiss tuberculosis reference laboratory. We examined mortality during treatment according to DST results and treatment adequacy in logistic regression models adjusting for sex, age, sputum microscopy and HIV status. Findings: 634 tuberculosis patients were included; median age was 33.2 years, 239 (37.7%) were female, 272 (42.9%) HIV-positive and 69 (10.9%) patients died. Based on the reference laboratory DST, 394 (62.2%) strains were pan-susceptible, 45 (7.1%) mono-resistant, 163 (25.7%) multidrug-resistant (MDR-TB), and 30 (4.7%) had pre-extensive or extensive drug resistance (pre-XDR/XDR-TB). Results of reference and local laboratories were discordant in 121 (19.1%) cases, corresponding to a sensitivity of 84.3% and a specificity of 90.8%. In patients with drug-resistant tuberculosis, discordant results were associated with increased mortality (risk ratio 1.81; 95% CI 1.07-3.07). In logistic regression, compared to adequately treated patients with pan-susceptible strains, the adjusted odds ratio for death was 4.23 (95% CI 2.16-8.29) for adequately treated patients with drug-resistant strains and 21.54 (95% CI 3.36-138.1) for inadequately treated patients with drug-resistant strains. HIV status was not associated with mortality. Interpretation: Using a reference laboratory standard, inaccurate DST leading to inappropriate treatment of drug-resistant tuberculosis, but not HIV infection, contributed to mortality. Funding Statement: This research was supported by the Swiss National Science Foundation (grant numbers 153442, 310030_166687 and 174281), the National Institutes of Allergy and Infectious Diseases (NIAID) under award numbers U01 AI096299, U01 AI069919, U01 AI069924, U01 AI069911, U01 AI069907, U01 AI096186, and U01 AI069923, and Swiss National Center for Mycobacteria, University of Zurich, Switzerland. Declaration of Interests: AA has received honoraria fees from Jensen-Cilag, Gilead and Bristol-Myers Squibb. All other authors have no conflicts of interest to declare. Ethics Approval Statement:Local institutional review boards or ethics committees approved the study at all participating sites. Informed consent was obtained where requested per local regulations. The study was also approved by the Cantonal Ethics Committee in Bern, Switzerland.

Drug resistance is best thought of as an ongoing biological process. Resistant bacteria must emerge, become established and ultimately transmit in order to be relevant to human health. In this context, genetic diversity can influence the rate and likelihood of resistance emerging; it can also modulate the net physiological impact of resistance and the propensity of an organism to improve any defects that arise from it. Combined, these effects can have an impact on a larger scale, with highly transmissible drug-resistant bacterial strains posing a formidable threat to global health. These considerations are pertinent to the future of tuberculosis control as well. In this chapter, we review our current understanding of the impact of genetic diversity in the broadest sense on the evolution of drug-resistant members of the Mycobacterium tuberculosis complex.

BackgroundUnderstanding transmission dynamics is useful for tuberculosis (TB) control. We conducted a population-based molecular epidemiological study to understand TB transmission in Ghana.MethodsMycobacterium tuberculosiscomplex (MTBC) isolates obtained from prospectively-sampled pulmonary TB patients between July, 2012 and December, 2015 were confirmed as MTBC using IS6110PCR. MTBC lineages were identified by large sequence polymorphism and single nucleotide polymorphism assays and further characterized using spoligotyping and standard 15-loci MIRU-VNTR typing. We used the n-1 method to estimate recent TB transmission and identified associated risk factors using logistic regression analysis.FindingsOut of 2,309 MTBC isolates, we identified 1,082 (46·9%) single cases with 1,227 (53·1%) isolates belonging to one of 276 clustered cases (clustering range; 2-35). Recent TB transmission rate was estimated to be 41·2%. While we see no significant difference in the recent transmission rates between lineages ofMycobacterium africanum(lineage-5 (31·8%); lineage-6 (24·7%), p=0·118), we found that lineage-4 belonging to theM. tuberculosistransmitted significantly higher (44·9%, pInterpretationsOur findings indicate high recent TB transmission suggesting occurrences of unsuspected outbreaks. The observed reduced transmission rate ofM. africanumsuggests other factor(s) may be responsible for its continuous presence in West Africa.FundingWellcome Trust Intermediate Fellowship Grant 097134/Z/11/Z to Dorothy Yeboah-Manu.

Mycobacterium africanum(Maf) causes up to half of human tuberculosis in West Africa, but little is known on this pathogen. We compared the genomes of 253Mafclinical isolates from Ghana, including both L5 and L6. We found that the genomic diversity of L6 was higher than in L5, and the selection pressures differed between both groups. Regulatory proteins appeared to evolve neutrally in L5 but under purifying selection in L6. Conversely, human T cell epitopes were under purifying selection in L5, but under positive selection in L6. Although only 10% of the T cell epitopes were variable, mutations were mostly lineage-specific. Our findings indicate thatMafL5 and L6 are genomically distinct, possibly reflecting different ecological niches.

Background. Whole-genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single-nucleotide polymorphism (SNP) typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. Methods. On the basis of genome sequences of 3 historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1642 patient isolates and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. Results. We identified 68 patients associated with the outbreak strain. Most received a tuberculosis diagnosis in 1991-1995, but cases were observed until 2011. Two thirds were homeless and/or substance abusers. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into 3 subclusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. Conclusions. Strain-specific SNP genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real time and at high resolution

Background Bovine tuberculosis (BTB) is caused by Mycobacterium bovis, which belongs to the Mycobacterium tuberculosis complex. Mycobacterium bovis have been described to be responsible of most cases of bovine tuberculosis. Although M. tuberculosis, M. africanum and non-complex mycobacteria were isolated from cattle. In Morocco, so far, no molecular studies were conducted to characterize the strains responsible of BTB. The present study aims to characterize M. bovis in Morocco. The present study was conducted in slaughterhouses in Rabat and El Jadida. Samples were collected from 327 slaughtered animals with visible lesions suggesting BTB. Results A total of 225 isolates yielded cultures, 95% (n = 215) of them were acid-fast (AF). Sixty eight per cent of the AF positive samples were confirmed as tuberculous mycobacteria (n = 147), 99% of these (n = 146) having RD9 and among the latter, 98% (n = 143) positive while 2% (n = 3) negative for RD4 A total of 134 samples were analyzed by spoligotyping of which 14 were in cluster and with 41 different spoligotypes, ten of them were new patterns (23%). The most prevalent spoligotypes were SB0121, SB0265, and SB0120, and were already identified in many other countries, such as Algeria, Spain, Tunisia, the United States and Argentina. Conclusion The shared borders between Algeria and Morocco, in addition to the previous importation of cattle from Europe and the US could explain the similarities found in M. bovis spoligotypes. On the other hand, the desert of Morocco could be considered as an efficient barrier preventing the introduction of BTB to Morocco from West Central and East Africa. Our findings suggest a low level endemic transmission of BTB similar to other African countries. However, more research is needed for further knowledge about the transmission patterns of BTB in Morocco.

Drug resistance is best thought of as an ongoing biological process. Resistant bacteria must emerge, become established and ultimately transmit in order to be relevant to human health. In this context, genetic diversity can influence the rate and likelihood of resistance emerging; it can also modulate the net physiological impact of resistance and the propensity of an organism to improve any defects that arise from it. Combined, these effects can have an impact on a larger scale, with highly transmissible drug-resistant bacterial strains posing a formidable threat to global health. These considerations are pertinent to the future of tuberculosis control as well. In this chapter, we review our current understanding of the impact of genetic diversity in the broadest sense on the evolution of drug-resistant members of the Mycobacterium tuberculosis complex.

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration., We thank S. Lecher, S. Li and J. Zallet for technical support. Calculations were performed at the sciCORE scientific computing core facility at the University of Basel. This work was supported by the Swiss National Science Foundation (grants 310030_166687 (S.G.) and 320030_153442 (M.E.) and Swiss HIV Cohort Study grant 740 to L.F.), the European Research Council (309540-EVODRTB to S.G.), TB-PAN-NET (FP7-223681 to S.N.), PathoNgenTrace projects (FP7-278864-2 to S.N.), SystemsX.ch (S.G.), the German Center for Infection Research (DZIF; S.N.), the Novartis Foundation (S.G.), the Natural Science Foundation of China (91631301 to Q.G.), and the National Institute of Allergy and Infectious Diseases (5U01-AI069924-05) of the US National Institutes of Health (M.E.).

Recent surveillance data of multidrug-resistant tuberculosis (MDR-TB) reported the highest rates of resistance ever documented. As further amplification of resistance in MDR strains of Mycobacterium tuberculosis occurs, extensively drug-resistant (XDR) and totally drug resistant (TDR) TB are beginning to emerge. Whilst for the most part, the epidemiological factors involved in the spread of MDR-TB are understood, insights into the bacterial drivers of MDR-TB have been gained only recently, largely owing to novel technologies and research in other organisms. Herein, we review recent findings on how bacterial factors such as persistence, hypermutation, the complex interrelationship between drug resistance and fitness, compensatory evolution, and epistasis affect the evolution of multidrug resistance in M. tuberculosis. Improved knowledge of these factors will help better predict the future trajectory of MDR-TB, and contribute to the development of new tools and strategies to combat this growing public health threat.

Rifampin resistance in clinical isolates of Mycobacterium tuberculosis arises primarily through the selection of bacterial variants harboring mutations in the 81-bp rifampin resistance-determining region of the rpoB gene. While these mutations were shown to infer a fitness cost in the absence of antibiotic pressure, compensatory mutations in rpoA and rpoC were identified which restore the fitness of rifampin-resistant bacteria carrying mutations in rpoB . To investigate the epidemiological relevance of these compensatory mutations, we analyzed 286 drug-resistant and 54 drug-susceptible clinical M. tuberculosis isolates from the Western Cape, South Africa, a high-incidence setting of multidrug-resistant tuberculosis. Sequencing of a portion of the RpoA-RpoC interaction region of the rpoC gene revealed that 23.5% of all rifampin-resistant isolates tested carried a nonsynonymous mutation in this region. These putative compensatory mutations in rpoC were associated with transmission, as 30.8% of all rifampin-resistant isolates with an IS 6110 restriction fragment length polymorphism (RFLP) pattern belonging to a recognized RFLP cluster harbored putative rpoC mutations. Such mutations were present in only 9.4% of rifampin-resistant isolates with unique RFLP patterns ( P < 0.01). Moreover, these putative compensatory mutations were associated with specific strain genotypes and the rpoB S531L rifampin resistance mutation. Among isolates harboring this rpoB mutation, 44.1% also harbored rpoC mutations, while only 4.1% of the isolates with other rpoB mutations exhibited mutations in rpoC ( P < 0.001). Our study supports a role for rpoC mutations in the transmission of multidrug-resistant tuberculosis and illustrates how epistatic interactions between drug resistance-conferring mutations, compensatory mutations, and different strain genetic backgrounds might influence compensatory evolution in drug-resistant M. tuberculosis .

Background Mycobacterium africanum comprises two phylogenetic lineages within the M. tuberculosis complex (MTBC) and is an important cause of human tuberculosis (TB) in West Africa. The reasons for this geographic restriction of M. africanum remain unclear. Here, we performed a prospective study to explore associations between the characteristics of TB patients and the MTBC lineages circulating in Ghana. Method We genotyped 1,211 MTBC isolates recovered from pulmonary TB patients recruited between 2012 and 2014 using single nucleotide polymorphism typing and spoligotyping. Associations between patient and pathogen variables were assessed using univariate and multivariate logistic regression. Results Of the 1,211 MTBC isolates analysed, 71.9 % (871) belonged to Lineage 4; 12.6 % (152) to Lineage 5 (also known as M. africanum West-Africa 1), 9.2 % (112) to Lineage 6 (also known as M. africanum West-Africa 2) and 0.6 % (7) to Mycobacterium bovis. Univariate analysis revealed that Lineage 6 strains were less likely to be isoniazid resistant compared to other strains (odds ratio = 0.25, 95 % confidence interval (CI): 0.05–0.77, P

SETTING: Madang and surroundings, Papua New Guinea (PNG). OBJECTIVE: To characterise the genetic diversity and drug susceptibility of Mycobacterium tuberculosis isolates collected in Madang and surroundings. DESIGN: M. tuberculosis was isolated from sputum samples from active pulmonary tuberculosis cases. Drug resistance profiles were obtained by drug susceptibility testing. M. tuberculosis lineages were identified by single nucleotide polymorphisms and sub-typing was performed by spoligotyping. Spoligotyping and 24 locus mycobacterial interspersed repetitive units-variable number of tandem repeats were combined to identify clustered isolates. RESULTS: The 173 M. tuberculosis isolates collected belonged predominantly to the Euro-American lineage (Lineage 4) and the East-Asian lineage (Lineage 2). Multidrug-resistant M. tuberculosis were observed in 5.2% of isolates. Lineage 2 M. tuberculosis, which includes the 'Beijing' genotype, was significantly associated with any drug resistance (OR 5.2, 95%CI 1.8-15.1). Cluster analyses showed 44% molecularly clustered isolates, suggesting transmission of M. tuberculosis in the community, including transmission of primary drug-resistant M. tuberculosis. CONCLUSION: These data provide the first insight into the molecular characteristics of M. tuberculosis in the Madang area of PNG, and indicate substantial drug resistance with evidence of ongoing transmission

Epidemics of drug-resistant bacteria emerge worldwide, even as resistant strains frequently have reduced fitness compared to their drug-susceptible counterparts. Data from model systems suggest that the fitness cost of antimicrobial resistance can be reduced by compensatory mutations; however, there is limited evidence that compensatory evolution has any significant role in the success of drug-resistant bacteria in human populations. Here we describe a set of compensatory mutations in the RNA polymerase genes of rifampicin-resistant M. tuberculosis, the etiologic agent of human tuberculosis (TB). M. tuberculosis strains harboring these compensatory mutations showed a high competitive fitness in vitro. Moreover, these mutations were associated with high fitness in vivo, as determined by examining their relative clinical frequency across patient populations. Of note, in countries with the world's highest incidence of multidrug-resistant (MDR) TB, more than 30% of MDR clinical isolates had this form of mutation. Our findings support a role for compensatory evolution in the global epidemics of MDR TB.

Mycobacterium tuberculosis harbours little DNA sequence diversity compared with other bacteria. However, there is mounting evidence that strain-to-strain variation in this organism has been underestimated. We review our current understanding of the genetic diversity among M. tuberculosis clinical strains and discuss the relevance of this diversity for the ongoing global epidemics of drug-resistant tuberculosis. Based on findings in other bacteria, we propose that epistatic interactions between pre-existing differences in strain genetic background, acquired drug-resistance-conferring mutations and compensatory changes could play a role in the emergence and spread of drug-resistant M. tuberculosis.

Summary In recent years, the percentage of tuberculosis (TB) cases in immigrants in Spain has increased markedly. In this context, discrimination between cases with potentially imported TB and cases likely to have acquired it by recent transmission after arrival is a basic issue. In this study, we evaluated molecular markers to obtain information on the geographic origin of Mycobacterium tuberculosis (MTB) strains and to determine whether they could help to identify imported TB cases. We analyzed the pks 15/1 gene, spoligotype, and MIRU locus 24 in 168 MTB isolates from immigrants of 32 nationalities and from a selection of autochthonous cases. In non-Asian immigrants we could not detect major differences with the autochthonous cases. For the Asian cases, we found some specific features for these markers. pks 15/1 was intact in 69% of the Asian patients, but this gene had a 7-bp deletion in all non-Asian and Spanish (non-Beijing) cases. The spoligotype-defined lineages EAI and CAS, and the allele with two repetitions in MIRU locus 24 was found exclusively among Asian immigrants. The analysis of certain bacterial molecular markers could help to discriminate between potentially imported TB cases and those more probably acquired in the host country.

Aims and objectives Resistant-, multidrug resistant- (MDR) and even extensively drug resistant- (XDR) Mycobacterium tuberculosis complex (MTBC) strains have emerged worldwide and represent a serious challenge for global tuberculosis (TB) control. Interestingly, high rates of MDR TB have been associated with particular phylogenetic lineages of the MTBC such as the Beijing lineage and the strong clonal expansion of particular MDR strains in Eastern Europe. However, a systematic investigation of the evolution and cross border transmission of particular highly transmissible MDR strains in Eastern Europe has not been performed so far. Methods To address this question, an in-depth analysis was performed of the population structure (based on 24 loci MIRU-VNTR-typing and spoligotyping) of strains from different settings in Eastern Europe, such as Uzbekistan, Turkmenistan, and Kazakhstan. Results The results showed that two major clones of the Beijing lineage are the drivers of the MDR epidemic in Eastern Europe. In addition, retrospective investigation of an MDR strain collection from Germany using MIRU-VNTR and IS6110 RFLP patterns revealed that both clones were circulating for more than 15 years. In-depth analysis of genome characteristics confirmed the longitudinal transmission of the two major MDR clones and allowed the identification of particular beneficial combinations of mutations potentially compensating for the fitness defect of the MDR/XDR strains. Moreover, these allele combinations coincide with an enhanced transmissibility compared with other circulating MDR/XDR strain types and promote the selection of extremely drug resistant but also highly virulent MTBC strain variants. Conclusion Taken together, these data argue strongly for an expansion of particular MDR Beijing clones in Eastern Europe that started at least 15 years ago. Evolution and selection over the past years has resulted in “hyper transmissible” MDR-TB strains that have the potential to spread rapidly and, thus, further complicate TB control.

BackgroundDifferent combinations of variable number of tandem repeat (VNTR) loci have been proposed for genotyping Mycobacterium tuberculosis complex (MTBC). Existing VNTR schemes show different discriminatory capacity among the six human MTBC lineages. Here, we evaluated the discriminatory power of a "customized MIRU12" loci format proposed previously by Comas et al. based on the standard 24 loci defined by Supply et al. for VNTR-typing of MTBC in Ghana.MethodOne hundred and fifty-eight MTBC isolates classified into Lineage 4 and Lineage 5 were used to compare a customized lineage-specific panel of 12 MIRU-VNTR loci ("customized MIRU-12") to the standard MIRU-15 genotyping scheme. The resolution power of each typing method was determined based on the Hunter-Gaston- Discriminatory Index (HGDI). A minimal set of customized MIRU-VNTR loci for typing Lineages 4 (Euro-American) and 5 (M. africanum West African 1) strains from Ghana was defined based on the cumulative HGDI.Results and conclusionAmong the 106 Lineage 4 strains, the customized MIRU-12 identified a total of 104 distinct genotypes consisting of 2 clusters of 2 isolates each (clustering rate 1.8%), and 102 unique strains while standard MIRU-15 yielded a total of 105 different genotypes, including 1 cluster of 2 isolates (clustering rate: 0.9%) and 104 singletons. Among, 52 Lineage 5 isolates, customized MIRU-12 genotyping defined 51 patterns with 1 cluster of 2 isolates (clustering rate: 0.9%) and 50 unique strains whereas MIRU-15 classified all 52 strains as unique. Cumulative HGDI values for customized MIRU-12 for Lineages 4 and 5 were 0.98 respectively whilst that of standard MIRU-15 was 0.99. A union of loci from the customised MIRU-12 and standard MIRU-15 revealed a set of customized eight highly discriminatory loci: 4052, 2163B, 40, 4165, 2165, 10,16 and 26 with a cumulative HGDI of 0.99 for genotyping Lineage 4 and 5 strains from Ghana.

The authors show that some mycobacteria carrying mutations conferring resistance to two antibiotics have a higher competitive fitness than corresponding strains carrying only one of these mutations. Moreover, the double-resistant strains exhibiting the highest competitive fitness in the laboratory are overrepresented in clinical settings with a high burden of extensively drug-resistant tuberculosis., Background and objectives: Multidrug resistant (MDR) bacteria are a growing threat to global health. Studies focusing on single antibiotics have shown that drug resistance is often associated with a fitness cost in the absence of drug. However, little is known about the fitness cost associated with resistance to multiple antibiotics. Methodology: We used Mycobacterium smegmatis as a model for human tuberculosis (TB) and an in vitro competitive fitness assay to explore the combined fitness effects and interaction between mutations conferring resistance to rifampicin (RIF) and ofloxacin (OFX); two of the most important first- and second-line anti-TB drugs, respectively. Results: We found that 4 out of 17 M. smegmatis mutants (24%) resistant to RIF and OFX showed a statistically significantly higher or lower competitive fitness than expected when assuming a multiplicative model of fitness effects of each individual mutation. Moreover, 6 out of the 17 double drug-resistant mutants (35%) had a significantly higher fitness than at least one of the corresponding single drug-resistant mutants. The particular combinations of resistance mutations associated with no fitness deficit in M. smegmatis were the most frequent among 151 clinical isolates of MDR and extensively drug-resistant (XDR) Mycobacterium tuberculosis from South Africa. Conclusions and implications: Our results suggest that epistasis between drug resistance mutations in mycobacteria can lead to MDR strains with no fitness deficit, and that these strains are positively selected in settings with a high burden of drug-resistant TB. Taken together, our findings support a role for epistasis in the evolution and epidemiology of MDR- and XDR-TB.

Tuberculosis caused 20% of all human deaths in the Western world between the seventeenth and nineteenth centuries and remains a cause of high mortality in developing countries. In analogy to other crowd diseases, the origin of human tuberculosis has been associated with the Neolithic Demographic Transition, but recent studies point to a much earlier origin. We analyzed the whole genomes of 259 M. tuberculosis complex (MTBC) strains and used this data set to characterize global diversity and to reconstruct the evolutionary history of this pathogen. Coalescent analyses indicate that MTBC emerged about 70,000 years ago, accompanied migrations of anatomically modern humans out of Africa and expanded as a consequence of increases in human population density during the Neolithic period. This long coevolutionary history is consistent with MTBC displaying characteristics indicative of adaptation to both low and high host densities.

Bacterial factors may contribute to the global emergence and spread of drug-resistant tuberculosis (TB). Only a few studies have reported on the interactions between different bacterial factors. We studied drug-resistant Mycobacterium tuberculosis isolates from a nationwide study conducted from 2000 to 2008 in Switzerland. We determined quantitative drug resistance levels of first-line drugs by using Bactec MGIT-960 and drug resistance genotypes by sequencing the hot-spot regions of the relevant genes. We determined recent transmission by molecular methods and collected clinical data. Overall, we analyzed 158 isolates that were resistant to isoniazid, rifampin, or ethambutol, 48 (30.4%) of which were multidrug resistant. Among 154 isoniazid-resistant strains, katG mutations were associated with high-level and inhA promoter mutations with low-level drug resistance. Only katG ( S315T ) (65.6% of all isoniazid-resistant strains) and inhA promoter −15C/T (22.7%) were found in molecular clusters. M. tuberculosis lineage 2 (includes Beijing genotype) was associated with any drug resistance (adjusted odds ratio [OR], 3.0; 95% confidence interval [CI], 1.7 to 5.6; P < 0.0001). Lineage 1 was associated with inhA promoter −15C/T mutations (OR, 6.4; 95% CI, 2.0 to 20.7; P = 0.002). We found that the genetic strain background influences the level of isoniazid resistance conveyed by particular mutations (interaction tests of drug resistance mutations across all lineages; P < 0.0001). In conclusion, M. tuberculosis drug resistance mutations were associated with various levels of drug resistance and transmission, and M. tuberculosis lineages were associated with particular drug resistance-conferring mutations and phenotypic drug resistance. Our study also supports a role for epistatic interactions between different drug resistance mutations and strain genetic backgrounds in M. tuberculosis drug resistance.

There is increasing evidence that strain variation in Mycobacterium tuberculosis complex (MTBC) might influence the outcome of tuberculosis infection and disease. To assess genotype-phenotype associations, phylogenetically robust molecular markers and appropriate genotyping tools are required. Most current genotyping methods for MTBC are based on mobile or repetitive DNA elements. Because these elements are prone to convergent evolution, the corresponding genotyping techniques are suboptimal for phylogenetic studies and strain classification. By contrast, single nucleotide polymorphisms (SNP) are ideal markers for classifying MTBC into phylogenetic lineages, as they exhibit very low degrees of homoplasy. In this study, we developed two complementary SNP-based genotyping methods to classify strains into the six main human-associated lineages of MTBC, the "Beijing" sublineage, and the clade comprising Mycobacterium bovis and Mycobacterium caprae. Phylogenetically informative SNPs were obtained from 22 MTBC whole-genome sequences. The first assay, referred to as MOL-PCR, is a ligation-dependent PCR with signal detection by fluorescent microspheres and a Luminex flow cytometer, which simultaneously interrogates eight SNPs. The second assay is based on six individual TaqMan real-time PCR assays for singleplex SNP-typing. We compared MOL-PCR and TaqMan results in two panels of clinical MTBC isolates. Both methods agreed fully when assigning 36 well-characterized strains into the main phylogenetic lineages. The sensitivity in allele-calling was 98.6% and 98.8% for MOL-PCR and TaqMan, respectively. Typing of an additional panel of 78 unknown clinical isolates revealed 99.2% and 100% sensitivity in allele-calling, respectively, and 100% agreement in lineage assignment between both methods. While MOL-PCR and TaqMan are both highly sensitive and specific, MOL-PCR is ideal for classification of isolates with no previous information, whereas TaqMan is faster for confirmation. Furthermore, both methods are rapid, flexible and comparably inexpensive.

1 Swiss Tropical and Public Health Institute, Basel, Switzerland, 2 University of Basel, Basel, Switzerland, 3 National Committee for Leprosy, Buruli Ulcer, Yaws and Leishmaniasis Control, Department of Disease Control, Ministry of Public Health, Yaounde, Cameroon, 4 Bankim District Hospital, Bankim, Cameroon, 5 Universitat Heidelberg, Heidelberg, Germany, 6 FAIRMED Africa Regional Office, Yaounde, Cameroon

BackgroundMycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as "Beijing" has repeatedly been associated with drug resistance and has been emerging in some parts of the world. "Beijing" strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical "Beijing" spoligotyping pattern.Methods and findingsMTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this "Pseudo-Beijing" spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical "Beijing" strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC.ConclusionsWe provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of "Beijing" strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.

IS6110, a Mycobacterium tuberculosis complex species-specific insertion element, is targeted primarily for DNA fingerprinting of M. tuberculosis strains. The number and chromosomal positions of copies of this element have been found to be highly variable between unrelated strains and have been exploited for molecular epidemiological purpose but the utility of IS6110 as an informative marker of strain phylogeny has yet to be demonstrated. In the current study, a recently proposed IS6110-targetting PCR based typing methodology, fluorescent amplified fragment length polymorphism (fAFLP) was applied to a global panel of 166 of the clinically more predominant ‘modern’ strains characterised by spoligotype and, where available, Variable Number Tandem Repeat (VNTR) to identify potentially evolutionarily informative common fragments that could define strains as belonging to established genetic lineages. These common fragments are hereby proposed to be ancestral insertion sites present in common ancestors of these strains rather than preferential insertion sites or ‘hot spots’. Results indicate that the exact same spoligotype and VNTR-defined lineages are reflected in the fragment patterns but with greater resolution and are able to clearly define the very distinct Haarlem, LAM, X and also the currently ill-defined T and S and lineages and spoligotypes designated as U, or unknown, without ambiguity. The biogeographical patterns generated reflect the migration of mankind across the globe and indicate that only four successful clones (or individual bacteria) gave rise to virtually all of the tuberculosis in Europe and the Americas. Potential lies in the application of the data to determine IS6110 evolutionary events that have occurred during the evolution of these lineages.

In this multicenter study, the reliability of two nonradiometric, fully automated systems, the MB/BacT and BACTEC MGIT 960 systems, for testing the susceptibilities of 82 Mycobacterium tuberculosis strains to isoniazid, rifampin, ethambutol, and streptomycin was evaluated in comparison with the radiometric BACTEC 460TB system. The arbitration of discrepant results was done by the reanalysis of the strain, the determination of the MIC, and the molecular characterization of some resistance determinants. The overall level of agreement with BACTEC 460TB results was 96% with the MB/BacT test and 97.2% with the BACTEC MGIT 960 system. With both methods, the level of agreement with BACTEC 460TB results was 96.3% for isoniazid, 98.8% for rifampin, and 98.8% for ethambutol. The level of agreement for streptomycin was 90.2% with MB/BacT and 97.5% with BACTEC MGIT 960. Overall, there were 11 very major errors and 2 major errors with the MB/BacT method and 5 very major errors and 2 major errors with the BACTEC MGIT 960 system. In general, the MB/BacT and BACTEC MGIT 960 systems showed good performance for susceptibility testing with first-line antituberculosis drugs.

Mycobacterium africanum is a member of the Mycobacterium tuberculosis complex (MTBC) and an important cause of human tuberculosis in West Africa that is rarely observed elsewhere. Here we genotyped 613 MTBC clinical isolates from Ghana, and searched for associations between the different phylogenetic lineages of MTBC and patient variables. We found that 17.1% (105/613) of the MTBC isolates belonged to M. africanum, with the remaining belonging to M. tuberculosis sensu stricto. No M. bovis was identified in this sample. M. africanum was significantly more common in tuberculosis patients belonging to the Ewe ethnic group (adjusted odds ratio: 3.02; 95% confidence interval: 1.67–5.47, p, Author Summary Tuberculosis remains one of the main global public health problems. Human tuberculosis is caused by bacteria known as the Mycobacterium tuberculosis complex (MTBC). The MTBC includes a variant called Mycobacterium africanum, which causes up to half of all tuberculosis cases in West Africa. For reasons unknown, M. africanum does not occur in other parts of the world. To explore the possible reasons for this geographic restriction of M. africanum, we analysed a large collection of bacterial strains isolated from tuberculosis patients in Ghana. We genetically characterized these bacterial isolates and collected relevant socio-demographic and epidemiological data. We found tuberculosis patients infected with M. africanum were more likely to belong to the Ewe ethnic group, compared to patients carrying other MTBC bacteria. The Ewes are indigenous inhabitants of coastal regions in West Africa that have previously been shown to harbour a high prevalence of M. africanum. Our findings support the hypothesis that different variants of MTBC have adapted to different human populations, and offer a possible explanation for the geographical restriction of M. africanum to West Africa.

Background High-throughput DNA sequencing produces vast amounts of data, with millions of short reads that usually have to be mapped to a reference genome or newly assembled. Both reference-based mapping and de novo assembly are computationally intensive, generating large intermediary data files, and thus require bioinformatics skills that are often lacking in the laboratories producing the data. Moreover, many research and practical applications in microbiology require only a small fraction of the whole genome data. Results We developed KvarQ, a new tool that directly scans fastq files of bacterial genome sequences for known variants, such as single nucleotide polymorphisms (SNP), bypassing the need of mapping all sequencing reads to a reference genome and de novo assembly. Instead, KvarQ loads “testsuites” that define specific SNPs or short regions of interest in a reference genome, and directly synthesizes the relevant results based on the occurrence of these markers in the fastq files. KvarQ has a versatile command line interface and a graphical user interface. KvarQ currently ships with two “testsuites” for Mycobacterium tuberculosis, but new “testsuites” for other organisms can easily be created and distributed. In this article, we demonstrate how KvarQ can be used to successfully detect all main drug resistance mutations and phylogenetic markers in 880 bacterial whole genome sequences. The average scanning time per genome sequence was two minutes. The variant calls of a subset of these genomes were validated with a standard bioinformatics pipeline and revealed >99% congruency. Conclusion KvarQ is a user-friendly tool that directly extracts relevant information from fastq files. This enables researchers and laboratory technicians with limited bioinformatics expertise to scan and analyze raw sequencing data in a matter of minutes. KvarQ is open-source, and pre-compiled packages with a graphical user interface are available at http://www.swisstph.ch/kvarq. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-881) contains supplementary material, which is available to authorized users.

The phylogeographic population structure of Mycobacterium tuberculosis suggests local adaptation to sympatric human populations. We hypothesized that HIV infection, which induces immunodeficiency, will alter the sympatric relationship between M. tuberculosis and its human host. To test this hypothesis, we performed a nine-year nation-wide molecular-epidemiological study of HIV–infected and HIV–negative patients with tuberculosis (TB) between 2000 and 2008 in Switzerland. We analyzed 518 TB patients of whom 112 (21.6%) were HIV–infected and 233 (45.0%) were born in Europe. We found that among European-born TB patients, recent transmission was more likely to occur in sympatric compared to allopatric host–pathogen combinations (adjusted odds ratio [OR] 7.5, 95% confidence interval [95% CI] 1.21–infinity, p = 0.03). HIV infection was significantly associated with TB caused by an allopatric (as opposed to sympatric) M. tuberculosis lineage (OR 7.0, 95% CI 2.5–19.1, p, Author Summary Human tuberculosis (TB) caused by Mycobacterium tuberculosis kills 1.5 million people each year. M. tuberculosis has been affecting humans for millennia, suggesting that different strain lineages may be adapted to specific human populations. The combination of a particular strain lineage and its corresponding patient population can be classified as sympatric (e.g. Euro-American lineage in Europeans) or allopatric (e.g. East-Asian lineage in Europeans). We hypothesized that infection with the human immunodeficiency virus (HIV), which impairs the human immune system, will interfere with this host–pathogen relationship. We performed a nation-wide molecular-epidemiological study of HIV–infected and HIV–negative TB patients between 2000 and 2008 in Switzerland. We found that HIV infection was associated with the less adapted allopatric lineages among patients born in Europe, and this was not explained by social or other patient factors such as increased social mixing in HIV–infected individuals. Strikingly, the association between HIV infection and less adapted M. tuberculosis lineages was stronger in patients with more pronounced immunodeficiency. Our observation was replicated in a second independent panel of M. tuberculosis strains collected during a population-based study in the Canton of Bern. In summary, our study provides evidence that the sympatric host–pathogen relationship in TB is disrupted by HIV infection.

Background Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown. Methods and Findings We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42–4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43–5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status. Conclusions We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian region.

Drug-resistant tuberculosis is a growing threat to global public health. Recent efforts to understand the evolution of drug resistance have shown that changes in drug–target interactions are only the first step in a longer adaptive process. The emergence of transmissible drug-resistant Mycobacterium tuberculosis is the result of a multitude of additional genetic mutations, many of which interact, a phenomenon known as epistasis. The varied effects of these epistatic interactions include compensating for the reduction of the biological cost associated with the development of drug resistance, increasing the level of resistance, and possibly accommodating broader changes in the physiology of resistant bacteria. Knowledge of these processes and our ability to detect them as they happen informs the development of diagnostic tools and better control strategies. In particular, the use of whole genome sequencing combined with surveillance efforts in the field could provide a powerful instrument to prevent future epidemics of drug-resistant tuberculosis.

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance(1). Drug resistance-conferring mutations frequently cause fitness costs in bacteria(2–5). Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations(6–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.

Summary 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 Cote 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.