Search

Background: Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Methods: We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains. Findings: Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08–2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28–1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13–2·12), independent of other patient and bacterial factors. Interpretation: Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented., The Lancet Microbe, 4 (7), ISSN:2666-5247

In settings with high tuberculosis (TB) endemicity, distinct genotypes of the Mycobacterium tuberculosis complex (MTBC) often differ in prevalence. However, the factors leading to these differences remain poorly understood. Here we studied the MTBC population in Dar es Salaam, Tanzania over a six-year period, using 1,082 unique patient-derived MTBC whole-genome sequences (WGS) and associated clinical data. We show that the TB epidemic in Dar es Salaam is dominated by multiple MTBC genotypes introduced to Tanzania from different parts of the world during the last 300 years. The most common MTBC genotypes deriving from these introductions exhibited differences in transmission rates and in the duration of the infectious period, but little differences in overall fitness, as measured by the effective reproductive number. Moreover, measures of disease severity and bacterial load indicated no differences in virulence between these genotypes during active TB. Instead, the combination of an early introduction and a high transmission rate accounted for the high prevalence of L3.1.1, the most dominant MTBC genotype in this setting. Yet, a longer co-existence with the host population did not always result in a higher transmission rate, suggesting that distinct life-history traits have evolved in the different MTBC genotypes. Taken together, our results point to bacterial factors as important determinants of the TB epidemic in Dar es Salaam., PLoS Pathogens, 19 (4), ISSN:1553-7374, ISSN:1553-7366

Genome-wide association studies rely on the statistical inference of untyped variants, called imputation, to increase the coverage of genotyping arrays. However, the results are often suboptimal in populations underrepresented in existing reference panels and array designs, since the selected single nucleotide polymorphisms (SNPs) may fail to capture population-specific haplotype structures, hence the full extent of common genetic variation. Here, we propose to sequence the full genomes of a small subset of an underrepresented study cohort to inform the selection of population-specific add-on tag SNPs and to generate an internal population-specific imputation reference panel, such that the remaining array-genotyped cohort could be more accurately imputed. Using a Tanzania-based cohort as a proof-of-concept, we demonstrate the validity of our approach by showing improvements in imputation accuracy after the addition of our designed add-on tags to the base H3Africa array., Author summary Genome-wide association studies, which study the association between genetic variants and various phenotypes, typically rely on genotyping arrays. Only a small proportion of genetic variants within the genome are typed on genotyping arrays. Untyped variants are statistically inferred through a process known as genotype imputation, where correlations between variants (haplotypes) observed in external reference panels are leveraged to infer untyped variants in the study population. However, for study populations that are underrepresented in existing reference panels, the quality of imputation is often sub-optimal. This is because typed variants incorporated on existing genotyping arrays can be unsuitable for the study population, and haplotype structures can be different between the reference and the study population. Here, we illustrate an approach to select a custom set of population-specific typed variants to improve genotype imputation in such underrepresented populations.

Background The bacteria that compose the Mycobacterium tuberculosis complex (MTBC) cause tuberculosis (TB) in humans and in different animals, including livestock. Much progress has been made in understanding the population structure of the human-adapted members of the MTBC by combining phylogenetics with genomics. Accompanying the discovery of new genetic diversity, a body of operational nomenclature has evolved to assist comparative and molecular epidemiological studies of human TB. By contrast, for the livestock-associated MTBC members, Mycobacterium bovis, M. caprae and M. orygis, there has been a lack of comprehensive nomenclature to accommodate new genetic diversity uncovered by emerging phylogenomic studies. We propose to fill this gap by putting forward a new nomenclature covering the main phylogenetic groups within M. bovis, M. caprae and M. orygis. Methods We gathered a total of 8,736 whole-genome sequences (WGS) from public sources and 39 newly sequenced strains, and selected a subset of 829 WGS, representative of the worldwide diversity of M. bovis, M. caprae and M. orygis. We used phylogenetics and genetic diversity patterns inferred from WGS to define groups. Results We propose to divide M. bovis, M. caprae and M. orygis in three main phylogenetic lineages, which we named La1, La2 and La3, respectively. Within La1, we identified several monophyletic groups, which we propose to classify into eight sublineages (La1.1-La1.8). These sublineages differed in geographic distribution, with some being geographically restricted and others globally widespread, suggesting different expansion abilities. To ease molecular characterization of these MTBC groups by the community, we provide phylogenetically informed, single nucleotide polymorphisms that can be used as barcodes for genotyping. These markers were implemented in KvarQ and TB-Profiler, which are platform-independent, open-source tools. Conclusions Our results contribute to an improved classification of the genetic diversity within the livestock-associated MTBC, which will benefit future molecular epidemiological and evolutionary studies.

22 páginas, 5 figuras, 2 tablas., Each day, approximately 27,000 people become ill with tuberculosis (TB), and 4,000 die from this disease. Pulmonary TB is the main clinical form of TB, and affects the lungs with a considerably heterogeneous manifestation among patients. Immunomodulation by an interplay of host-, environment-, and pathogen-associated factors partially explains such heterogeneity. Microbial communities residing in the host's airways have immunomodulatory effects, but it is unclear if the inter-individual variability of these microbial communities is associated with the heterogeneity of pulmonary TB. Here, we investigated this possibility by characterizing the microbial composition in the sputum of 334 TB patients from Tanzania, and by assessing its association with three aspects of disease manifestations: sputum mycobacterial load, severe clinical findings, and chest x-ray (CXR) findings. Compositional data analysis of taxonomic profiles based on 16S-rRNA gene amplicon sequencing and on whole metagenome shotgun sequencing, and graph-based inference of microbial associations revealed that the airway microbiome of TB patients was shaped by inverse relationships between Streptococcus and two anaerobes: Selenomonas and Fusobacterium. Specifically, the strength of these microbial associations was negatively correlated with Faith's phylogenetic diversity (PD) and with the accumulation of transient genera. Furthermore, low body mass index (BMI) determined the association between abnormal CXRs and community diversity and composition. These associations were mediated by increased abundance of Selenomonas and Fusobacterium, relative to the abundance of Streptococcus, in underweight patients with lung parenchymal infiltrates and in comparison to those with normal chest x-rays. And last, the detection of herpesviruses and anelloviruses in sputum microbial assemblage was linked to co-infection with HIV. Given the anaerobic metabolism of Selenomonas and Fusobacterium, and the hypoxic environment of lung infiltrates, our results suggest that in underweight TB patients, lung tissue remodeling toward anaerobic conditions favors the growth of Selenomonas and Fusobacterium at the expense of Streptococcus. These new insights into the interplay among particular members of the airway microbiome, BMI, and lung parenchymal lesions in TB patients, add a new dimension to the long-known association between low BMI and pulmonary TB. Our results also drive attention to the airways virome in the context of HIV-TB coinfection., This work was supported by grants from the Bernese Lung Association, the Swiss National Science Foundation (grant nos. 310030_188888, IZRJZ3_164171, IZLSZ3_170834, and CRSII5_177163), and the European Research Council (309540-EVODRTB). The Tuberculosis Cohort Study in Dar es Salaam, Tanzania, is funded by the Rudolf Geigy Foundation, Basel, Switzerland. MT is supported by a Ph.D. fellowship from the Swiss Institute of Bioinformatics (funded by the Rudolf Geigy Foundation).

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.

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.

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this Review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient, and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations, and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

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

The Lineage 2–Beijing (L2–Beijing) sub-lineage of Mycobacterium tuberculosis has received much attention due to its high virulence, fast disease progression, and association with antibiotic resistance. Despite several reports of the recent emergence of L2–Beijing in Africa, no study has investigated the evolutionary history of this sub-lineage on the continent. In this study, we used whole genome sequences of 781 L2 clinical strains from 14 geographical regions globally distributed to investigate the origins and onward spread of this lineage in Africa. Our results reveal multiple introductions of L2–Beijing into Africa linked to independent bacterial populations from East- and Southeast Asia. Bayesian analyses further indicate that these introductions occurred during the past 300 years, with most of these events pre-dating the antibiotic era. Hence, the success of L2–Beijing in Africa is most likely due to its hypervirulence and high transmissibility rather than drug resistance., Frontiers in Ecology and Evolution, 7, ISSN:2296-701X

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.

Figure S3. Plot of 100 RTL decays (-r =10) for the Influenza dataset (black dots). All decays are similar, but there is a larger variability compared to the TB dataset. Red dots indicate the range of RTL among the different runs for the corresponding iteration, the variability among runs increases steadily and it reaches a maximum of 2.5% when the tree is reduced to 40% of the RTL. (PDF 56 kb)

Figure S2. Plot of 100 RTL decays (-r =100) for the TB dataset (black dots). All decays are similar. Red dots indicates the range of RTL among the different runs for the corresponding iteration, the variability among runs increases slowly at first, it reaches a maximum of 1% when the tree is reduced to 20% of the RTL. (PDF 25 kb)

Figure S1. Root-to-tip regression for the complete influenza dataset (2063 leaves, black dots and black regression line, in gray the 99% confidence interval) and four reduced trees (orange: 99% or RTL , red: 90% of RTL, green: 75% of RTL, blue: 50% of RTL. All trees were re-rooted with the best-fit method implemented in TempEst. (PDF 12 kb)

Antibiotic resistance is one of the biggest threats to human health globally. Alarmingly, multidrug-resistant and extensively drug-resistant Mycobacterium tuberculosis have now spread worldwide. Some key antituberculosis antibiotics are prodrugs, for which resistance mechanisms are mainly driven by mutations in the bacterial enzymatic pathway required for their bioactivation. We have developed drug-like molecules that activate a cryptic alternative bioactivation pathway of ethionamide in M. tuberculosis, circumventing the classic activation pathway in which resistance mutations have now been observed. The first-of-its-kind molecule, named SMARt-420 (Small Molecule Aborting Resistance), not only fully reverses ethionamide-acquired resistance and clears ethionamide-resistant infection in mice, it also increases the basal sensitivity of bacteria to ethionamide.

OBJECTIVES For tuberculosis (TB) transmission to occur, an uninfected individual must inhale the previously infected breath. Our objective was to identify potential TB transmission hotspots in metropolitan city of Dar es Salaam, Tanzania and to model the annual risk of TB transmission in different locations of public importance. METHODS We collected indoor carbon dioxide (CO2) data from markets, prisons, night clubs, public transportation, religious and social halls, and from schools. Study volunteers recorded social contacts at each of the locations. We then estimated the annual risks of TB transmission using a modified Wells-Riley equation for different locations. RESULTS The annual risks of TB transmission were highest among prison inmates (41.6%) and drivers (20.3%) in public transport. Lower transmission risks were found in central markets (4.8% for traders, but 0.5% for their customers), passengers on public transport (2.4%), public schools (4.0%), nightclubs (1.7%), religious (0.13%), and social halls (0.12%). CONCLUSION For the first time in a country representative of sub-Saharan Africa, we modelled the risk of TB transmission in important public locations by using a novel approach of studying airborne transmission. This approach can guide identification of TB transmission hotspots and targeted interventions to reach WHO's ambitious End TB targets.

Zoonotic tuberculosis is a risk for human health, especially when animals are in close contact with humans. Mycobacterium tuberculosis was cultured from several organs, including lung tissue and gastric mucosa, of three captive elephants euthanized in a Swiss zoo. The elephants presented weight loss, weakness and exercise intolerance. Molecular characterization of the M. tuberculosis isolates by spoligotyping revealed an identical profile, suggesting a single source of infection. Multilocus variable-number of tandem-repeat analysis (MLVA) elucidated two divergent populations of bacteria and mixed infection in one elephant, suggesting either different transmission chains or prolonged infection over time. A total of eight M. tuberculosis isolates were subjected to whole-genome sequence (WGS) analysis, confirming a single source of infection and indicating the route of transmission between the three animals. Our findings also show that the methods currently used for epidemiological investigations of M. tuberculosis infections should be carefully applied on isolates from elephants. Moreover the importance of multiple sampling and analysis of within-host mycobacterial clonal populations for investigations of transmission is demonstrated.

Treatment of multidrug-resistant Mycobacterium tuberculosis is a challenge. This letter describes the emergence of resistance to new therapies, bedaquiline and delamanid.

Pyrazinamide (PZA) plays the important role in shortening the tuberculosis treatment period and in treating MDR-TB. Phenotypic PZA susceptibility methods are limited because they require specialized acidified media, which increases costs and complexity. In this study we developed a genotypic high resolution melt (HRM) analysis technique to detect pncA mutations associated with PZA resistant Mycobacterium tuberculosis. Seven overlapping primer pairs were designed to cover the entire pncA gene and upstream regions. Each gene segment was individually amplified by real-time PCR followed by HRM analysis. The assay was evaluated on 98 clinical M. tuberculosis isolates (41 PZA susceptible by MGIT method, 55 PZA resistant, 2 undetermined). HRM was 94% concordant to full-length sequencing results, with most discrepancies attributable to mixed populations per HRM or transversions. Sequencing and HRM yielded 82% and 84% concordance, respectively, to phenotypic PZA susceptibilities by MGIT, with most discrepancies attributable to isolates with wild-type pncA but phenotypic PZA resistance. This HRM technique is a simple and high-throughput method for screening clinical M. tuberculosis samples for PZA resistance.

There is considerable variability in the outcome of Mycobacterium tuberculosis infection. We hypothesized that Mycobacterium africanum was less likely than M. tuberculosis to transmit and progress to tuberculosis disease.In a cohort study of patients with tuberculosis and their household contacts in The Gambia, we categorized 1808 HIV-negative tuberculosis contacts according to exposure to M. tuberculosis or M. africanum. Positive skin test results indicated transmission, and development of tuberculosis during 2 years of follow-up indicated progression to disease.Transmission rates were similar, but rates of progression to disease were significantly lower in contacts exposed to M. africanum than in those exposed to M. tuberculosis (1.0% vs. 2.9%; hazard ratio [HR], 3.1 [95% confidence interval {CI}, 1.1-8.7]). Within M. tuberculosis sensu stricto, contacts exposed to a Beijing family strain were most likely to progress to disease (5.6%; HR relative to M. africanum, 6.7 [95% CI, 2.0-22]).M. africanum and M. tuberculosis transmit equally well to household contacts, but contacts exposed to M. africanum are less likely to progress to tuberculosis disease than those exposed to M. tuberculosis. The variable rate of progression by lineage suggests that tuberculosis variability matters in clinical settings and should be accounted for in studies evaluating tuberculosis vaccines and treatment regimens for latent tuberculosis infection.

Mycobacterium tuberculosis remains a major cause of morbidity and mortality worldwide. Studies have reported human pathogens to have geographically structured population genetics, some of which have been linked to ancient human migrations. However, no study has addressed the potential evolutionary consequences of such longstanding human-pathogen associations. Here, we demonstrate that the global population structure of M. tuberculosis is defined by six phylogeographical lineages, each associated with specific, sympatric human populations. In an urban cosmopolitan environment, mycobacterial lineages were much more likely to spread in sympatric than in allopatric patient populations. Tuberculosis cases that did occur in allopatric hosts disproportionately involved high-risk individuals with impaired host resistance. These observations suggest that mycobacterial lineages are adapted to particular human populations. If confirmed, our findings have important implications for tuberculosis control and vaccine development.

Table S1. Characteristics of the TB patients with and without follow-up. (DOCX 15 kb)

Whole genome sequencing allows rapid detection of drug-resistant isolates. However, the availability of high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data has thus far been limited.We determined drug resistance profiles of 176 genetically diverse clinical isolates 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 DST results with predicted drug resistance profiles inferred by whole genome sequencing.Classification of strains by the two phenotypic DST methods into resistotype/wild type populations was concordant in 73-99 % of cases, depending on the drug. Our data suggests that the established critical concentration (5 mg/L) for ethambutol resistance (MGIT 960 system) is too high and may misclassify strains as susceptible, compared to 7H10 agar dilution. Increased minimal inhibitory concentrations were explained by mutations identified by whole genome sequencing. Using whole genome sequences, we were able to predict quantitative drug resistance levels for the majority of drug resistance mutations. Predicting quantitative levels of drug resistance by whole genome sequencing was partially limited due to incompletely understood drug resistance mechanisms. The overall sensitivity and specificity of whole genome-based DST were 86.8 % and 94.5 %, respectively.Despite some limitations, whole genome sequencing has the potential to infer resistance profiles without the need for time-consuming phenotypic methods.

We analysed 312 drug-resistant genomes of Mycobacterium tuberculosis (Mtb) collected from HIV coinfected and HIV negative TB patients from nine countries with a high tuberculosis burden. We found that rifampicin-resistant Mtb 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.

LSP-PCR and LRPS sample analyses. (DOC 43 kb)

Supplementary text, figures and tables. In this document we elaborate on the efficacy of treatment (Section 1), expand on the reasoning behind our choice of variant calling approach, describe the outcome of population simulations providing the basis for our demographic expectations, illustrate the phylogeny of our isolates, and finally describe the data contained in the other additional files. The display items in the document include the following tables and figures. Table S1 Phenotypic drug susceptibility testing for first line drugs. Table S2 Treatment combinations given to patients. Table S3 Average sequencing depth of individual samples. Table S4 Variations in the number of detected f-SNPs in patients at each time point is consistent with v-SNP dynamics. Table S5 Transition probabilities for the detection of v-SNPs. Figure S1 Treatment efficacy in patient groups. Figure S2 Error profiles of simulated reads. Figure S3 Culture-induced minor variants in MTBC. Figure S4 Retained minor variants in our patient population. Figure S5 Demographic trends of MTBC populations within the host are consistent with bacterial killing. Figure S6 Sequencing depth does not bias v-SNP detection. Figure S7 The number of detected v-SNPs for each patient varies over time. Figure S8 Cumulative distributions of allele frequencies. Figure S9 Frequency trajectories of recurrent alleles vary, showing a dynamic and heterogeneous population. Figure S10 Sweeps of v-SNPs occur independently of the type of v-SNP affected the association of v-SNPs with drug resistance. Figure S11 Most patients are infected with strains of the Beijing family. (PDF 532 kb)

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

Fluorescence applification curve readings: Specific Rv004c and Rv0129c primer/probe set containing lineage specific SNP were used to analyze 70 and 33 MTBC isolates respectively. Amplification fluorescence curves for Rv004c and Rv0129 target region are shown in panel (a) and panel (b) respectively. (DOC 123 kb)

LRPS is more sensitive than LSP-PCR typing: H37Rv genomic DNA was serially diluted and analyzed with RD 724 (N = 10) or RD 750 (N = 10) deletion primers: lane M = molecular weight markers, lane 1 positive control (MTB Uganda family or lineage 3 strain), lane 2 negative control, lane 3- 12 H37Rv DNA diluted from 100-10 ng; the minimum dilution of genomic DNA that can be amplified was in lane 10 (30ng) (Plate a & c). Plate b & d are the corresponding LRPS using Rv004c/probe or Rv0129c/probe set with H37Rv genomic DNA diluted from (10-1) ng per assay, the minimum dilution of DNA that can be amplified was 3ng (See arrow on the amplification curves plate b & d). (DOC 179 kb)

Primer/probe used in genotyping MTB. (DOCX 21 kb)

Mutations in the genes of the F420 signaling pathway, including dnn, fgd1, fbiA, fbiB, fbiC, and fbiD, of Mycobacterium tuberculosis (Mtb) complex can lead to delamanid resistance. We searched for such mutations among 129 Mtb 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 ten strains with mutations, we determined the minimum inhibitory concentration (MIC) of delamanid. We found one strain from a delamanid-naïve patient carrying the natural polymorphism Tyr29del (ddn) that was associated with a critical MIC to delamanid.

BACKGROUND: Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen. METHODS: Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes. RESULTS: Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p

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., PLoS Biology, 17 (5), ISSN:1544-9173, ISSN:1545-7885

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.

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

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