Search

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) that vary in virulence and transmissibility. While genome-wide association studies have uncovered several mutations conferring drug resistance, much less is known about the factors underlying other bacterial phenotypes. Variation in the outcome of tuberculosis infection and diseases has been attributed primarily to patient and environmental factors, but recent evidence indicates an additional role for the genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the effect of genetic variation on the strain-specific metabolic adaptive capacity and vulnerability. To define the functionality of single-nucleotide polymorphisms (SNPs) systematically, we developed a constraint-based approach that integrates metabolomic and genomic data. Our model-based predictions correctly classify SNP effects in pyruvate kinase and suggest a genetic basis for strain-specific inherent baseline susceptibility to the antibiotic para-aminosalicylic acid. Our method is broadly applicable across microbial life, opening possibilities for the development of more selective treatment strategies.

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.

Coscolla et al., 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., M.C. is supported by the Ramón y Cajal programme from the Ministerio de Ciencia, Innovación y Universidades. This work was supported by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) (research award to M.C.), Ministerio de Ciencia, Innovación y Universidades (grant number RTI2018-094399-A-I00 to M.C.) and Consellería de Educació de la Generalitat Valenciana (grant number SEJI/2019/011 to M.C.), the Swiss National Science Foundation (grants 310030_188888, IZRJZ3_164171, IZLSZ3_170834 and CRSII5_177163 to S. G.), the European Research Council (883582-ECOEVODRTB to S. G.) and Wellcome (grant number 097134/Z/11/Z to D. Y.-M).

24 páginas, 3 figuras, 1 tabla. The sequence data generated by this study has been deposited on SRA (https://www.ncbi.nlm.nih.gov/sra) under the accession number PRJNA670836. Extended data is available here: https://github.com/fmenardo/MTBC_L1_L3. DOI: https://doi.org/10.5281/zenodo.4609804 (Menardo, 2021)., 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., This work was supported by the Swiss National Science Foundation (grants 310030_188888, CRSII5_177163, IZRJZ3_164171 and IZLSZ3_170834) and the European Research Council (309540 EVODRTB and 883582-ECOEVODRTB).

Lineage 1 (L1) and 3 (L3) are two lineages of theMycobacterium tuberculosiscomplex (MTBC), causing tuberculosis (TB) in humans. L1 and L3 are endemic to 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. Here we analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We show that South Asia played a central role in the dispersion of these two lineages to neighboring regions. Moreover, we found that L1 exhibits signatures of local adaptation at theesxHlocus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. 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.

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

SummaryThe Mycobacterium tuberculosis (Mtb) complex comprises seven phylogenetically distinct human-adapted lineages exhibiting different geographical distribution and degrees of pathogenicity. Among these, Lineage 1 (L1) has been associated with low virulence whereas Lineage 2 (L2) has been linked to hyper-virulence, enhanced transmission and drug resistance. Here, we conducted multi-layer comparative analyses using whole genome sequencing data combined with quantitative transcriptomic and proteomic profiling of a set of L1 and L2 clinical strains, each grown under two different conditions in vitro. Our data revealed different degrees of correlation between transcript and protein abundances across clinical strains and functional gene categories, indicating variable levels of post-transcriptional regulation in the tested lineages. Contrasting genomic and gene expression data showed that the magnitude of the transcriptional and translational changes was proportional to the phylogenetic distance between strains, with one out of three single nucleotide polymorphisms leading to a transcriptional and/or translational change on average. We devised a new genome-scale transcriptional regulatory model and identified several master transcription factors, strongly linked to the sigma factor network, whose targets were differentially regulated between the two lineages. These differences resulted in a higher basal expression of DosR proteins and a stronger response to nitric oxide (NO) exposure in L2 compared to L1. These patterns are most likely responsible for the shorter NO-induced growth arrest in L2 observed. Given the limited genetic variation between strains, it appears that phenotypic differences in Mtb are substantially driven by differences in the regulation of biochemical networks through master transcriptional regulators.HighlightsProteomic and transcriptomic characterization of fully sequenced diverse L1 and L2 clinical isolates of Mtb.Post-transcriptional control mechanisms for regulatory and virulence genes are mitigated in Mtb L2.By applying a genome-scale transcriptional framework, DosR, Rv1985c, Lsr2 and Rv0691c are identified as master transcription factors responsible for differential target gene expression in L2 strains compared to L1.L1 and L2 DosR proteins respond differently to nitric oxide stress, thus determining a relevant phenotype.

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.

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.

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., Some types of mutations occur more frequently than expected. This study shows that such bias —an excess of transitions over transversions—influences the evolution of antibiotic resistance in a key global pathogen, Mycobacterium tuberculosis.

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.

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.

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.

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.

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.

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.

Background Monitoring drug resistance in Mycobacterium tuberculosis is essential to curb the spread of tuberculosis (TB). Unfortunately, drug susceptibility testing is currently not available in Papua New Guinea (PNG) and that impairs TB control in this country. We report for the first time M. tuberculosis mutations associated with resistance to first and second-line anti-TB drugs in Madang, PNG. A molecular cluster analysis was performed to identify M. tuberculosis transmission in that region. Results Phenotypic drug susceptibility tests showed 15.7% resistance to at least one drug and 5.2% multidrug resistant (MDR) TB. Rifampicin resistant strains had the rpoB mutations D516F, D516Y or S531L; Isoniazid resistant strains had the mutations katG S315T or inhA promoter C15T; Streptomycin resistant strains had the mutations rpsL K43R, K88Q, K88R), rrs A514C or gidB V77G. The molecular cluster analysis indicated evidence for transmission of resistant strain. Conclusions We observed a substantial rate of MDR-TB in the Madang area of PNG associated with mutations in specific genes. A close monitoring of drug resistance is therefore urgently required, particularly in the presence of drug-resistant M. tuberculosis transmission. In the absence of phenotypic drug susceptibility testing in PNG, molecular assays for drug resistance monitoring would be of advantage.