Search

Abstract The Mycobacterium tuberculosis complex (MTBC) comprises nine human-adapted lineages that differ in their geographical distribution. Local adaptation of specific MTBC genotypes to the respective human host population has been invoked in this context. We aimed to assess if bacterial genetics governs MTBC pathogenesis or if local co-adaptation translates into differential susceptibility of human macrophages to infection by different MTBC genotypes. We generated macrophages from cryopreserved blood mononuclear cells of Tanzanian tuberculosis patients, from which the infecting MTBC strains had previously been phylogenetically characterized. We infected these macrophages ex vivo with a phylogenetically similar MTBC strain (“matched infection”) or with strains representative of other MTBC lineages (“mismatched infection”). We found that L1 infections resulted in a significantly lower bacterial burden and that the intra-cellular replication rate of L2 strains was significantly higher compared the other MTBC lineages, irrespective of the MTBC lineage originally infecting the patients. Moreover, L4-infected macrophages released significantly greater amounts of TNF-α, IL-6, IL-10, MIP-1β, and IL-1β compared to macrophages infected by all other strains. While our results revealed no measurable effect of local adaptation, they further highlight the strong impact of MTBC phylogenetic diversity on the variable outcome of the host–pathogen interaction in human tuberculosis.

ABSTRACT Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) remains one of the deadliest infectious diseases globally. Timely diagnosis is a key step in the management of TB patients and in the prevention of further transmission events. Current diagnostic tools are limited in these regards. There is an urgent need for new accurate non-sputum-based diagnostic tools for the detection of symptomatic as well as subclinical TB. In this study, we recruited 52 symptomatic TB patients (sputum Xpert MTB/RIF positive) and 58 household contacts to assess the accuracy of a sequence-specific hybridization assay that detects the presence of Mtb cell-free DNA in urine. Using sputum Xpert MTB/RIF as a reference test, the magnetic bead-capture assay could discriminate active TB from healthy household contacts with an overall sensitivity of 72.1% [confidence interval (CI) 0.59–0.86] and specificity of 95.5% (CI 0.90–1.02) with a positive predictive value of 93.9% and negative predictive value of 78.2%. The detection of Mtb-specific DNA in urine suggested four asymptomatic TB infection cases that were confirmed in all instances either by concomitant Xpert MTB/RIF sputum testing or by follow-up investigation raising the specificity of the index test to 100%. We conclude that sequence-specific hybridization assays on urine specimens hold promise as non-invasive tests for the detection of subclinical TB.IMPORTANCEThere is an urgent need for a non-sputum-based diagnostic tool allowing sensitive and specific detection of all forms of tuberculosis (TB) infections. In that context, we performed a case-control study to assess the accuracy of a molecular detection method enabling the identification of cell-free DNA from Mycobacterium tuberculosis that is shed in the urine of tuberculosis patients. We present accuracy data that would fulfill the target product profile for a non-sputum test. In addition, recent epidemiological data suggested that up to 50% of individuals secreting live bacilli do not present with symptoms at the time of screening. We report, here, that the investigated index test could also detect instances of asymptomatic TB infections among household contacts.

Thanks to its expertise in clinical research, epidemiology, infectious diseases, microbiology, parasitology, public health, translational research and tropical medicine, coupled with deeply rooted partnerships with institutions in low- and middle-income countries (LMICs), the Swiss Tropical and Public Health Institute (Swiss TPH) has been a key contributor in many drug research and development consortia involving academia, pharma and product development partnerships. Our know-how of the maintenance of parasites and their life-cycles in the laboratory, plus our strong ties to research centres and disease control programme managers in LMICs with access to field sites and laboratories, have enabled systems for drug efficacy testing in vitro and in vivo, clinical research, and modelling to support the experimental approaches. Thus, Swiss TPH has made fundamental contributions towards the development of new drugs – and the better use of old drugs – for neglected tropical diseases and infectious diseases of poverty, such as Buruli ulcer, Chagas disease, food-borne trematodiasis (e.g. clonorchiasis, fascioliasis and opisthorchiasis), human African trypanosomiasis, leishmaniasis, malaria, schistosomiasis, soil-transmitted helminthiasis and tuberculosis. In this article, we show case the success stories of molecules to which Swiss TPH has made a substantial contribution regarding their use as anti-infective compounds with the ultimate aim to improve people’s health and well-being.

T cell activation markers (TAM) expressed by antigen-specific T cells constitute promising candidates to attest the presence of an active infection by Mycobacterium tuberculosis (Mtb). Reciprocally, their modulation may be used to assess antibiotic treatment efficacy and eventually attest disease resolution. We hypothesized that the phenotype of Mtb-specific T cells may be quantitatively impacted by the load of bacteria present in a patient. We recruited 105 Tanzanian adult tuberculosis (TB) patients and obtained blood before and after 5 months of antibiotic treatment. We studied relationships between patients' clinical characteristics of disease severity and microbiological as well as molecular proxies of bacterial load in sputum at the time of diagnosis. Besides, we measured by flow cytometry the expression of CD38 or CD27 on CD4+ T cells producing interferon gamma (IFN-γ) and/or tumor necrosis factor alpha (TNF-α) in response to a synthetic peptide pool covering the sequences of Mtb antigens ESAT-6, CFP-10, and TB10.4. Reflecting the difficulty to extrapolate bacterial burden from a single end-point read-out, we observed statistically significant but weak correlations between Xpert MTB/RIF, molecular bacterial load assay and time to culture positivity. Unlike CD27, the resolution of CD38 expression by antigen-specific T cells was observed readily following 5 months of antibiotic therapy. However, the intensity of CD38-TAM signals measured at diagnosis did not significantly correlate with Mtb 16S RNA or rpoB DNA detected in patients' sputa. Altogether, our data support CD38-TAM as an accurate marker of infection resolution independently of sputum bacterial load.

Abstract Background Differences in rural and urban settings could account for distinct characteristics in the epidemiology of tuberculosis (TB). We comparatively studied epidemiological features of TB and helminth co-infections in adult patients from rural and urban settings of Tanzania. Methods Adult patients (≥ 18 years) with microbiologically confirmed pulmonary TB were consecutively enrolled into two cohorts in Dar es Salaam, with ~ 4.4 million inhabitants (urban), and Ifakara in the sparsely populated Kilombero District with ~ 400 000 inhabitants (rural). Clinical data were obtained at recruitment. Stool and urine samples were subjected to diagnose helminthiases using Kato-Katz, Baermann, urine filtration, and circulating cathodic antigen tests. Differences between groups were assessed by χ 2, Fisher’s exact, and Wilcoxon rank sum tests. Logistic regression models were used to determine associations. Results Between August 2015 and February 2017, 668 patients were enrolled, 460 (68.9%) at the urban and 208 (31.1%) at the rural site. Median patient age was 35 years (interquartile range [IQR]: 27–41.5 years), and 454 (68%) were males. Patients from the rural setting were older (median age 37 years vs. 34 years, P = 0.003), had a lower median body mass index (17.5 kg/m2 vs. 18.5 kg/m2, P

BACKGROUND:Helminth infections can negatively affect the immunologic host control, which may increase the risk of progression from latent Mycobacterium tuberculosis infection to tuberculosis (TB) disease and alter the clinical presentation of TB. We assessed the prevalence and determined the clinical relevance of helminth co-infection among TB patients and household contact controls in urban Tanzania. METHODOLOGY:Between November 2013 and October 2015, we enrolled adult (≥18 years) sputum smear-positive TB patients and household contact controls without TB during an ongoing TB cohort study in Dar es Salaam, Tanzania. We used Baermann, FLOTAC, Kato-Katz, point-of-care circulating cathodic antigen, and urine filtration to diagnose helminth infections. Multivariable logistic regression models with and without random effects for households were used to assess for associations between helminth infection and TB. PRINCIPAL FINDINGS:A total of 597 TB patients and 375 household contact controls were included. The median age was 33 years and 60.2% (585/972) were men. The prevalence of any helminth infection among TB patients was 31.8% (190/597) and 25.9% (97/375) among controls. Strongyloides stercoralis was the predominant helminth species (16.6%, 161), followed by hookworm (9.0%, 87) and Schistosoma mansoni (5.7%, 55). An infection with any helminth was not associated with TB (adjusted odds ratio (aOR) 1.26, 95% confidence interval (CI): 0.88-1.80, p = 0.22), but S. mansoni infection was (aOR 2.15, 95% CI: 1.03-4.45, p = 0.040). Moreover, S. mansoni infection was associated with lower sputum bacterial load (aOR 2.63, 95% CI: 1.38-5.26, p = 0.004) and tended to have fewer lung cavitations (aOR 0.41, 95% CI: 0.12-1.16, p = 0.088). CONCLUSIONS/SIGNIFICANCE:S. mansoni infection was an independent risk factor for active TB and altered the clinical presentation in TB patients. These findings suggest a role for schistosomiasis in modulating the pathogenesis of human TB. Treatment of helminths should be considered in clinical management of TB and TB control programs.

Decentralization of Directly Observed Treatment (DOT) for tuberculosis (TB) to the community (home-based DOT) has improved the coverage of TB treatment and reduced the burden to the health care facilities (facility-based DOT). We aimed to compare TB treatment outcomes in home-based and facility-based DOT under programmatic conditions in an urban setting with a high TB burden.A retrospective analysis of a cohort of adult TB patients (≥15 years) routinely notified between 2010 and 2013 in two representative TB sub-districts in the Temeke district, Dar es Salaam, Tanzania. We assessed differences in treatment outcomes by calculating Risk Ratios (RRs). We used logistic regression to assess the association between DOT and treatment outcomes.Data of 4,835 adult TB patients were analyzed, with a median age of 35 years, 2,943 (60.9%) were men and TB/HIV co-infection prevalence of 39.9%. A total of 3,593 (74.3%) patients were treated under home-based DOT. Patients on home-based DOT were more likely to die compared to patients on facility-based DOT (RR 2.04, 95% Confidence Interval [95% CI]: 1.52-2.73), and more likely to complete TB treatment (RR 1.14, 95% CI: 1.06-1.23), but less likely to have a successful treatment outcome (RR 0.94, 95% CI: 0.92-0.97). Home-based DOT was preferred by women (adjusted Odds Ratio [aOR] 1.55, 95% CI: 1.34-1.80, p

The genetic diversity of 134 serogroup X Neisseria meningitis isolates from Africa, Europe, and North America was analyzed by multilocus sequence typing and pulsed-field gel electrophoresis. Although most European and American isolates were highly diverse, one clonal grouping was identified in sporadic disease and carrier strains isolated over the last 2 decades in the United Kingdom, the Netherlands, Germany, and the United States. In contrast to the diversity in the European and American isolates, most carrier and disease isolates recovered during the last 30 years in countries in the African meningitis belt belonged to a second clonal grouping. During the last decade, these bacteria have caused meningitis outbreaks in Niger and Ghana. These results support the development of a comprehensive conjugate vaccine that would include serogroup X polysaccharide.

The aim of the present study was to determine whether there is a correlation between phylogenetic relationship and inflammatory response amongst a panel of clinical isolates representative of the global diversity of the human Mycobacterium tuberculosis Complex (MTBC). Measurement of cytokines from infected human peripheral blood monocyte-derived macrophages revealed a wide variation in the response to different strains. The same pattern of high or low response to individual strains was observed for different pro-inflammatory cytokines and chemokines, and was conserved across multiple human donors. Although each major phylogenetic lineage of MTBC included strains inducing a range of cytokine responses, we found that overall inflammatory phenotypes differed significantly across lineages. In particular, comparison of evolutionarily modern lineages demonstrated a significant skewing towards lower early inflammatory response. The differential response to ancient and modern lineages observed using GM-CSF derived macrophages was also observed in autologous monocyte-derived dendritic cells and murine bone marrow-derived macrophages, but not in human unfractionated peripheral blood mononuclear cells. We hypothesize that the reduced immune responses to modern lineages contribute to more rapid disease progression and transmission, which might be a selective advantage in the context of expanding human populations. In addition to the lineage effects, the large strain-to-strain variation in innate immune responses elicited by MTBC will need to be considered in tuberculosis vaccine development.

The Kassena-Nankana District of northern Ghana lies in the African "meningitis belt" where epidemics of meningococcal meningitis have been reoccurring every eight to 12 years for the last 100 years. The dynamics of meningococcal colonisation and disease are incompletely understood, and hence we embarked on a long-term study to determine how levels of colonisation with different bacterial serogroups change over time, and how the patterns of disease relate to such changes.Between February 1998 and November 2005, pharyngeal carriage of Neisseria meningitidis in the Kassena-Nankana District was studied by twice-yearly colonisation surveys. Meningococcal disease was monitored throughout the eight-year study period, and patient isolates were compared to the colonisation isolates. The overall meningococcal colonisation rate of the study population was 6.0%. All culture-confirmed patient isolates and the majority of carriage isolates were associated with three sequential waves of colonisation with encapsulated (A ST5, X ST751, and A ST7) meningococci. Compared to industrialised countries, the colonising meningococcal population was less constant in genotype composition over time and was genetically less diverse during the peaks of the colonisation waves, and a smaller proportion of the isolates was nonserogroupable. We observed a broad age range in the healthy carriers, resembling that of meningitis patients during large disease epidemics.The observed lack of a temporally stable and genetically diverse resident pharyngeal flora of meningococci might contribute to the susceptibility to meningococcal disease epidemics of residents in the African meningitis belt. Because capsular conjugate vaccines are known to impact meningococcal carriage, effects on herd immunity and potential serogroup replacement should be monitored following the introduction of such vaccines.

A previous study in Pará, Northern Brazil, described a strain of Mycobacterium tuberculosis with a unique genotype (SIT2517/T1) associated with multidrug-resistant tuberculosis (MDR-TB). To improve our understanding of MDR-TB transmission dynamics of these strains within this region, we performed phenotypic and genotypic drug susceptibility testing (pDST/gDST), 24-loci mycobacterial interspersed repetitive units (MIRU-VNTR) genotyping, whole-genome sequencing (WGS) and geo-epidemiology analysis. Of the 28 SIT2517/T1 isolates, 19 (67.9%) could be genotyped by 24-loci MIRU-VNTR and 15 by WGS. All belonged to sublineage 4.1.1.3, distinct from other representative Lineage 4 isolates identified in Brazil. The MDR phenotype determined by pDST was confirmed by gDST, the latter also demonstrating the presence of additional mutations conferring pre-extensively drug-resistance (pre-XDR). Discrepancies between gDST and pDST were observed for pyrazinamide and fluoroquinolones. Thirteen out of 15 isolates analyzed by WGS were clustered when applying a 12 single nucleotide polymorphisms (SNPs) cutoff. The SIT2517/T1 isolates were distributed across the metropolitan regions of Belém and Collares municipalities, showing no geographic clustering. WGS-transmission network analysis revealed a high likelihood of direct transmission and the formation of two closely linked transmission chains. This study highlights the need to implement TB genomic surveillance in the Brazilian Amazon region.

Persons living with HIV are known to be at increased risk of developing tuberculosis (TB) disease upon infection with Mycobacterium tuberculosis (Mtb). However, it has remained unclear how HIV co-infection affects subsequent Mtb transmission from these patients. Here, we customized a Bayesian phylodynamic framework to estimate the effects of HIV co-infection on the Mtb transmission dynamics from sequence data. We applied our model to four Mtb genomic datasets collected in sub-Saharan African countries with a generalized HIV epidemic. Our results confirm that HIV co-infection is a strong risk factor for developing active TB. Additionally, we demonstrate that HIV co-infection is associated with a reduced effective reproductive number for TB. Stratifying the population by CD4+ T-cell count yielded similar results, suggesting that, in this context, CD4+ T-cell count is not a better predictor of Mtb transmissibility than HIV infection status alone. Together, our genome-based analyses complement observational household contact studies, and more firmly establish the negative association between HIV co-infection and Mtb transmissibility.

Geographical hotspots with high frequency of multi-drug resistant tuberculosis (MDR-TB) have been observed in several locations, such as the country of Georgia. Here, the authors analyse genomic sequences from tuberculosis bacteria isolated from Georgia to show that the transmission fitness of MDR-TB strains is heterogeneous, and highly drug-resistant and transmissible isolates contribute to the emergence and maintenance of MDR-TB hotspots.

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

BackgroundThe epidemiology of Mycobacterium tuberculosis complex (MTBC) lineage 5 (L5) infections in Ghana revealed a significantly increased prevalence in Ewes compared to other self-reported ethnic groups. In that context, we sought to investigate the early phase of tuberculosis (TB) infection using ex vivo infection of macrophages derived from the blood of Ewe and Akan ethnic group volunteers with MTBC L4 and L5 strains.MethodsThe study participants consisted of 16 controls, among which self-reported Akan and Ewe ethnicity was equally represented, as well as 20 cured TB cases consisting of 11 Akans and 9 Ewes. Peripheral blood mononuclear cells were isolated from both healthy controls and cured TB cases. CD14+ monocytes were isolated and differentiated into monocyte-derived macrophages (MDMs) before infection with L4 or L5 endemic strains. The bacterial load was assessed after 2 hours (uptake) as well as 3 and 7 days post-infection.ResultsWe observed a higher capacity of MDMs from Ewes to phagocytose L4 strains (p < 0.001), translating into a higher bacillary load on day 7 (p < 0.001) compared to L5, despite the higher replication rate of L5 in Ewe MDMs (fold change: 1.4 vs. 1.2, p = 0.03) among the controls. On the contrary, within macrophages from Akans, we observed a significantly higher phagocytic uptake of L5 (p < 0.001) compared to L4, also translating into a higher load on day 7 (p = 0.04). However, the replication rate of L4 in Akan MDMs was higher than that of L5 (fold change: L4 = 1.2, L4 = 1.1, p = 0.04). Although there was no significant difference in the uptake of L4 and L5 among cured TB cases, there was a higher bacterial load of both L4 (p = 0.02) and L5 (p = 0.02) on day 7 in Ewe MDMs.ConclusionOur results suggest that host ethnicity (driven by host genetic diversity), MTBC genetic diversity, and individual TB infection history are all acting together to modulate the outcome of macrophage infections by MTBC.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

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.

BackgroundResistance to tuberculosis (TB) drugs has become a major threat to global control efforts. Early case detection and drug susceptibility profiling of the infecting bacteria are essential for appropriate case management. The objective of this study was to determine the drug susceptibility profiles of difficult-to-treat (DTT) TB patients in Ghana.MethodsSputum samples obtained from DTT-TB cases from health facilities across Ghana were processed for rapid diagnosis and detection of drug resistance using the Genotype MTBDRplus and Genotype MTBDRsl.v2 from Hain Life science.ResultsA total of 298 (90%) out of 331 sputum samples processed gave interpretable bands out of which 175 (58.7%) were resistant to at least one drug (ANYr); 16.8% (50/298) were isoniazid-mono-resistant (INHr), 16.8% (50/298) were rifampicin-mono-resistant (RIFr), and 25.2% (75/298) were MDR. 24 (13.7%) of the ANYr were additionally resistant to at least one second line drug: 7.4% (2 RIFr, 1 INHr, and 10 MDR samples) resistant to only FQs and 2.3% (2 RIFr, 1 INHr, and 1 MDR samples) resistant to AMG drugs kanamycin (KAN), amikacin (AMK), capreomycin (CAP), and viomycin (VIO). Additionally, there were 4.0% (5 RIFr and 2 MDR samples) resistant to both FQs and AMGs. 81 (65.6%) out of 125 INH-resistant samples including INHr and MDR had katG-mutations (MT) whereas 15 (12%) had inhApro-MT. The remaining 28 (22.4%) had both katG and inhA MT. All the 19 FQ-resistant samples were gyrA mutants whereas the 10 AMGs were rrs (3), eis (3) as well as rrs, and eis co-mutants (4). Except for the seven pre-XDR samples, no sample had eis MT.ConclusionThe detection of several pre-XDR TB cases in Ghana calls for intensified drug resistance surveillance and monitoring of TB patients to, respectively, ensure early diagnosis and treatment compliance.

Nine-banded armadillos (Dasypus novemcinctus) are naturally infected with Mycobacterium leprae and are implicated in the zoonotic transmission of leprosy in the United States. In Mexico, the existence of such a reservoir remains to be characterized. We describe a wild armadillo infected by M. leprae in the state of Nuevo León, Mexico.

Abstract CD4 T cell phenotyping-based blood assays have the potential to meet WHO target product profiles (TPP) of non-sputum-biomarker-based tests to diagnose tuberculosis (TB). Yet, substantial refinements are required to allow their implementation in clinical settings. This study assessed the real time performance of a simplified T cell activation marker (TAM)-TB assay to detect TB in adults from one millilitre of blood with a 24 h turnaround time. We recruited 479 GeneXpert positive cases and 108 symptomatic but GeneXpert negative controls from presumptive adult TB patients in the Temeke District of Dar-es-Salaam, Tanzania. TAM-TB assay accuracy was assessed by comparison with a composite reference standard comprising GeneXpert and solid culture. A single millilitre of fresh blood was processed to measure expression of CD38 or CD27 by CD4 T cells producing IFN-γ and/or TNF-α in response to a synthetic peptide pool covering the sequences of Mycobacterium tuberculosis (Mtb) ESAT-6, CFP-10 and TB10.4 antigens on a 4-color FACSCalibur apparatus. Significantly superior to CD27 in accurately diagnosing TB, the CD38-based TAM-TB assay specificity reached 93.4% for a sensitivity of 82.2% with an area under the receiver operating characteristics curve of 0.87 (95% CI 0.84–0.91). The assay performance was not significantly affected by HIV status. To conclude, we successfully implemented TAM-TB immunoassay routine testing with a 24 h turnaround time at district level in a resource limited setting. Starting from one millilitre of fresh blood and being not influenced by HIV status, TAM-TB assay format and performance appears closely compatible with the optimal TPP accuracy criteria defined by WHO for a non-sputum confirmatory TB test.

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

Polyclonal infections occur when at least two unrelated strains of the same pathogen are detected in an individual. Here, Moreno-Molina et al. analyse sputum and surgical resections from tuberculosis patients, showing that the magnitude of polyclonal infections can be underestimated when only testing sputum samples.

Summary: Background: Molecular diagnostics are considered the most promising route to achievement of rapid, universal drug susceptibility testing for Mycobacterium tuberculosis complex (MTBC). We aimed to generate a WHO-endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction. Methods: In this systematic analysis, we used a candidate gene approach to identify mutations associated with resistance or consistent with susceptibility for 13 WHO-endorsed antituberculosis drugs. We collected existing worldwide MTBC whole-genome sequencing data and phenotypic data from academic groups and consortia, reference laboratories, public health organisations, and published literature. We categorised phenotypes as follows: methods and critical concentrations currently endorsed by WHO (category 1); critical concentrations previously endorsed by WHO for those methods (category 2); methods or critical concentrations not currently endorsed by WHO (category 3). For each mutation, we used a contingency table of binary phenotypes and presence or absence of the mutation to compute positive predictive value, and we used Fisher's exact tests to generate odds ratios and Benjamini-Hochberg corrected p values. Mutations were graded as associated with resistance if present in at least five isolates, if the odds ratio was more than 1 with a statistically significant corrected p value, and if the lower bound of the 95% CI on the positive predictive value for phenotypic resistance was greater than 25%. A series of expert rules were applied for final confidence grading of each mutation. Findings: We analysed 41 137 MTBC isolates with phenotypic and whole-genome sequencing data from 45 countries. 38 215 MTBC isolates passed quality control steps and were included in the final analysis. 15 667 associations were computed for 13 211 unique mutations linked to one or more drugs. 1149 (7·3%) of 15 667 mutations were classified as associated with phenotypic resistance and 107 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was more than 80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were identified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs. Interpretation: We present the first WHO-endorsed catalogue of molecular targets for MTBC drug susceptibility testing, which is intended to provide a global standard for resistance interpretation. The existence of this catalogue should encourage the implementation of molecular diagnostics by national tuberculosis programmes. Funding: Unitaid, Wellcome Trust, UK Medical Research Council, and Bill and Melinda Gates Foundation.

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.

Some strains of Mycobacterium tuberculosis seem to be able to avoid host defense systems. Here the authors stratify patients by severity of tuberculosis and find correlations with the level of IL-1β production by macrophages exposed to these isolates.

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.

Introduction Mycobacterial Interspersed Repetitive Units–Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia. Methods We used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method. Results For the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each. Conclusion Our study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

IntroductionMycobacterial Interspersed Repetitive Units-Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia.MethodsWe used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method.ResultsFor the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each.ConclusionOur study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

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.

Summary: Background: South Africa has a high burden of rifampicin-resistant tuberculosis (including multidrug-resistant [MDR] tuberculosis), with increasing rifampicin-monoresistant (RMR) tuberculosis over time. Resistance acquisition during first-line tuberculosis treatment could be a key contributor to this burden, and HIV might increase the risk of acquiring rifampicin resistance. We assessed whether HIV during previous treatment was associated with RMR tuberculosis and resistance acquisition among a retrospective cohort of patients with MDR or rifampicin-resistant tuberculosis. Methods: In this retrospective cohort study, we included all patients routinely diagnosed with MDR or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa, between Jan 1, 2008, and Dec 31, 2017. Patient-level data were obtained from a prospective database, complemented by data on previous tuberculosis treatment and HIV from a provincial health data exchange. Stored MDR or rifampicin-resistant tuberculosis isolates from patients underwent whole-genome sequencing (WGS). WGS data were used to infer resistance acquisition versus transmission, by identifying genomically unique isolates (single nucleotide polymorphism threshold of five). Logistic regression analyses were used to assess factors associated with RMR tuberculosis and genomic uniqueness. Findings: The cohort included 2041 patients diagnosed with MDR or rifampicin-resistant tuberculosis between Jan 1, 2008, and Dec 31, 2017; of those, 463 (22·7%) with RMR tuberculosis and 1354 (66·3%) with previous tuberculosis treatment. In previously treated patients, HIV positivity during previous tuberculosis treatment versus HIV negativity (adjusted odds ratio [OR] 2·07, 95% CI 1·35–3·18), and three or more previous tuberculosis treatment episodes versus one (1·96, 1·21–3·17) were associated with RMR tuberculosis. WGS data showing MDR or rifampicin-resistant tuberculosis were available for 1169 patients; 360 (30·8%) isolates were identified as unique. In previously treated patients, RMR tuberculosis versus MDR tuberculosis (adjusted OR 4·96, 3·40–7·23), HIV positivity during previous tuberculosis treatment (1·71, 1·03–2·84), and diagnosis in 2013–17 (1·42, 1·02–1·99) versus 2008–12, were associated with uniqueness. In previously treated patients with RMR tuberculosis, HIV positivity during previous treatment (adjusted OR 5·13, 1·61–16·32) was associated with uniqueness as was female sex (2·50 [1·18–5·26]). Interpretation: These data suggest that HIV contributes to rifampicin-resistance acquisition during first-line tuberculosis treatment and that this might be driving increasing RMR tuberculosis over time. Large-scale prospective cohort studies are required to further quantify this risk. Funding: Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust.

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

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.

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

A variety of biological differences exist between strains and phylogenetic lineages in the Mycobacterium tuberculosis complex (MTBC). Here, the authors perform RNAseq and methylome analysis of the main clades of the MTBC and show variation in transcriptional profiles is mainly due to mutations in transcriptional start sites.

The fast evolution of pathogenic viruses has allowed for the development of phylodynamic approaches that extract information about the epidemiological characteristics of viral genomes. Thanks to advances in whole genome sequencing, they can be applied to slowly evolving bacterial pathogens like Mycobacterium tuberculosis. In this study, we investigate and compare the epidemiological dynamics underlying two M. tuberculosis outbreaks using phylodynamic methods. Specifically, we (i) test if the outbreak data sets contain enough genetic variation to estimate short-term evolutionary rates and (ii) reconstruct epidemiological parameters such as the effective reproduction number.The first outbreak occurred in the Swiss city of Bern (1987–2012) and was caused by a drug-susceptible strain belonging to the phylogenetic M. tuberculosis Lineage 4. The second outbreak was caused by a multidrug-resistant (MDR) strain of Lineage 2, imported from the Wat Tham Krabok (WTK) refugee camp in Thailand into California.There is little temporal signal in the Bern data set and moderate temporal signal in the WTK data set. Thanks to its high sampling proportion (90%) the Bern outbreak allows robust estimation of epidemiological parameters despite the poor temporal signal. Conversely, there is much uncertainty in the epidemiological estimates concerning the sparsely sampled (9%) WTK outbreak. Our results suggest that both outbreaks peaked around 1990, although they were only recognized as outbreaks in 1993 (Bern) and 2004 (WTK). Furthermore, individuals were infected for a significantly longer period (around 9 years) in the WTK outbreak than in the Bern outbreak (4–5 years).Our work highlights both the limitations and opportunities of phylodynamic analysis of outbreaks involving slowly evolving pathogens: (i) estimation of the evolutionary rate is difficult on outbreak time scales and (ii) a high sampling proportion allows quantification of the age of the outbreak based on the sampling times, and thus allows for robust estimation of epidemiological parameters. Keywords: Phylodynamic analysis, Tuberculosis outbreak, Epidemic dynamics, Transmission dynamics

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.

ObjectiveTo determine the prevalence of tuberculosis (TB) disease and infection as well as incident TB disease among people who use drugs (PWUD) attending Medication Assisted Treatment (MAT) clinics in Dar-es-Salaam, Tanzania.MethodsIn this prospective cohort study, a total of 901 consenting participants were enrolled from November 2016 to February 2017 and a structured questionnaire administered to them through the open data kit application on android tablets. Twenty-two months later, we revisited the MAT clinics and reviewed 823 of the 901 enrolled participant's medical records in search for documentation on TB disease diagnosis and treatment. Medical records reviewed included those of participants whom at enrolment were asymptomatic, not on TB disease treatment, not on TB preventive therapy and those who had a documented tuberculin skin test (TST) result.ResultsOf the 823 medical records reviewed 22 months after enrolment, 42 had documentation of being diagnosed with TB disease and initiated on TB treatment. This is equivalent to a TB disease incidence rate of 2,925.2 patients per 100,000 person years with a total follow up time of 1,440 person-years. At enrolment the prevalence of TB disease and TB infection was 2.6% and 54% respectively and the HIV prevalence was 44% and 16% among females and males respectively.ConclusionPWUD attending MAT clinics bear an extremely high burden of TB and HIV and are known to have driven TB epidemics in a number of countries. Our reported TB disease incidence is 12 times that of the general Tanzanian incidence of 237 per 100,000 further emphasizing that this group should be prioritized for TB screening, testing and treatment. Gender specific approaches should also be developed as female PWUDs are markedly more affected with HIV and TB disease than male PWUDs.

BACKGROUND The evaluation of procedures for drug susceptibility prediction of Mycobacterium tuberculosis based on genomic data against the conventional reference method test based on culture is realistic considering the scenario of growing number of tools proposals based on whole-genome sequences (WGS). OBJECTIVES This study aimed to evaluate drug susceptibility testing (DST) outcome based on WGS tools and the phenotypic methods performed on isolates of M. tuberculosis Lineage 1 from the state of Pará, Brazil, generally associated with low levels of drug resistance. METHODOLOGY Culture based DST was performed using the Proportion Method in Löwenstein-Jensen medium on 71 isolates that had been submitted to WGS. We analysed the seven main genome sequence-based tools for resistance and lineage prediction applied to M. tuberculosis and for comparison evaluation we have used the Kappa concordance test. FINDINGS When comparing the WGS-based tools against the DST, we observed the highest level of agreement using TB-profiler. Among the tools, TB-profiler, KvarQ and Mykrobe were those which identified the largest number of TB-MDR cases. Comparing the four most sensitive tools regarding resistance prediction, agreement was observed for 43 genomes. MAIN CONCLUSIONS Drug resistance profiling using next-generation sequencing offers rapid assessment of resistance-associated mutations, therefore facilitating rapid access to effective treatment.

Findings from previous comparative genomics studies of the Mycobacterium tuberculosis complex (MTBC) suggest genomic variation among the genotypes may have phenotypic implications. We investigated the diversity in the phenotypic profiles of the main prevalent MTBC genotypes in West Africa. Thirty-six whole genome sequenced drug susceptible MTBC isolates belonging to lineages 4, 5 and 6 were included in this study. The isolates were phenotypically characterized for urease activity, tween hydrolysis, Thiophen-2-Carboxylic Acid Hydrazide (TCH) susceptibility, nitric oxide production, and growth rate in both liquid (7H9) and solid media (7H11 and Löwenstein-Jensen (L-J)). Lineage 4 isolates showed the highest growth rate in both liquid (p = 0.0003) and on solid (L-J) media supplemented with glycerol (p

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

Genetic studies of both the human host and Mycobacterium tuberculosis (MTB) demonstrate independent association with tuberculosis (TB) risk. However, neither explains a large portion of disease risk or severity. Based on studies in other infectious diseases and animal models of TB, we hypothesized that the genomes of the two interact to modulate risk of developing active TB or increasing the severity of disease, when present. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which there were 3 MTB lineages of which L4-Ugandan (L4.6) is the most recent. TB severity, measured using the Bandim TBscore, was modeled as a function of host SNP genotype, MTB lineage, and their interaction, within two independent cohorts of TB cases, N = 113 and 121. No association was found between lineage and severity, but association between multiple polymorphisms in IL12B and TBscore was replicated in two independent cohorts (most significant rs3212227, combined p = 0.0006), supporting previous associations of IL12B with TB susceptibility. We also observed significant interaction between a single nucleotide polymorphism (SNP) in SLC11A1 and the L4-Ugandan lineage in both cohorts (rs17235409, meta p = 0.0002). Interestingly, the presence of the L4-Uganda lineage in the presence of the ancestral human allele associated with more severe disease. These findings demonstrate that IL12B is associated with severity of TB in addition to susceptibility, and that the association between TB severity and human genetics can be due to an interaction between genes in the two species, consistent with host-pathogen coevolution in TB.

Mycobacterium tuberculosis (Mtb) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb.

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

Abstract Background Culture contamination with environmental bacteria is a major challenge in tuberculosis (TB) laboratories in hot and humid climate zones. We studied the effect of cetylpyridinium chloride (CPC) preservation on culture results and performance of Xpert MTB/RIF. Methods Consecutive sputum samples from microscopy smear-positive TB patients were collected. Two-hundred samples were equally split in two aliquots, one aliquot was treated with CPC and stored at ambient temperature for 7 days. The second aliquot was immediately processed. Samples were decontaminated for 20, 15 or 10 min, and subsequently cultured on Löwenstein-Jensen medium. Furthermore, 50 samples were stored for 7, 14 and 21 days, and 100 CPC-pretreated samples tested by Xpert MTB/RIF. Results CPC pretreated samples showed a higher culture yield compared to non-treated sputum samples across all decontamination times: 94% vs. 73% at 10 min (p = 0.01), 94% vs. 64% at 15 min (p = 0.004), and 90% vs. 52% at 20 min (p

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

Abstract Background Combination therapy is one of the most effective tools for limiting the emergence of drug resistance in pathogens. Despite the widespread adoption of combination therapy across diseases, drug resistance rates continue to rise, leading to failing treatment regimens. The mechanisms underlying treatment failure are well studied, but the processes governing successful combination therapy are poorly understood. We address this question by studying the population dynamics of Mycobacterium tuberculosis within tuberculosis patients undergoing treatment with different combinations of antibiotics. Results By combining very deep whole genome sequencing (~1000-fold genome-wide coverage) with sequential sputum sampling, we were able to detect transient genetic diversity driven by the apparently continuous turnover of minor alleles, which could serve as the source of drug-resistant bacteria. However, we report that treatment efficacy has a clear impact on the population dynamics: sufficient drug pressure bears a clear signature of purifying selection leading to apparent genetic stability. In contrast, M. tuberculosis populations subject to less drug pressure show markedly different dynamics, including cases of acquisition of additional drug resistance. Conclusions Our findings show that for a pathogen like M. tuberculosis, which is well adapted to the human host, purifying selection constrains the evolutionary trajectory to resistance in effectively treated individuals. Nonetheless, we also report a continuous turnover of minor variants, which could give rise to the emergence of drug resistance in cases of drug pressure weakening. Monitoring bacterial population dynamics could therefore provide an informative metric for assessing the efficacy of novel drug combinations.

The molecular clock and its phylogenetic applications to genomic data have changed how we study and understand one of the major human pathogens, Mycobacterium tuberculosis (MTB), the etiologic agent of tuberculosis. Genome sequences of MTB strains sampled at different times are increasingly used to infer when a particular outbreak begun, when a drug-resistant clone appeared and expanded, or when a strain was introduced into a specific region. Despite the growing importance of the molecular clock in tuberculosis research, there is a lack of consensus as to whether MTB displays a clocklike behavior and about its rate of evolution. Here we performed a systematic study of the molecular clock of MTB on a large genomic data set (6,285 strains), covering different epidemiological settings and most of the known global diversity. We found that sampling times below 15-20 years were often insufficient to calibrate the clock of MTB. For data sets where such calibration was possible, we obtained a clock rate between 1x10-8 and 5x10-7 nucleotide changes per-site-per-year (0.04-2.2 SNPs per-genome-per-year), with substantial differences between clades. These estimates were not strongly dependent on the time of the calibration points as they changed only marginally when we used epidemiological isolates (sampled in the last 40 years) or three ancient DNA samples (about 1,000 years old) to calibrate the tree. Additionally, the uncertainty and the discrepancies in the results of different methods were sometimes large, highlighting the importance of using different methods, and of considering carefully their assumptions and limitations.

[This corrects the article DOI: 10.1371/journal.ppat.1006939.].