Search

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.

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

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

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

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

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

INTRODUCTION:Recurrent tuberculosis (TB) is one of the main challenges in TB control. Genotyping based on Mycobacterial Interspersed Repetitive Units-Variable Tandem Repeats (MIRU-VNTR) has been widely used to differentiate between relapse and reinfection, which are the two main causes of recurrent TB. There is a lack of data regarding the causes of TB recurrence in Georgia, and while differentiating between relapse and reinfection plays a key role in defining appropriate interventions, the required genotyping methodologies have not been implemented. The objective of this study was to implement MIRU-VNTR genotyping at the National Center for Tuberculosis and Lung Diseases (NCTBLD) and differentiate between relapse and reinfection in multidrug resistant (MDR-) TB patients from Tbilisi, Georgia. METHODS:Recurrent MDR tuberculosis cases from 2014-2016 diagnosed at NCTLD were included in the study when bacterial samples from both episodes were available. Genotyping based on the MIRU-VNTR 24 loci was implemented and used for differentiating between relapse and reinfection. Paired samples showing the same MIRU-VNTR pattern or one locus difference were classified as relapse, while two and more loci differences were treated as reinfection. Exact logistic regression was used to identify predictors of recurrence. RESULTS:Thirty two MDR-TB patients (64 samples) were included and MIRU-VNTR 24 typing was performed on the corresponding paired samples. Of the 32 patients, 25 (83.3%) were identified as relapse while 5 (16.7%) were due to re-infection. Patients with a history of incarceration were significantly associated with TB reinfection (p< 0.05). CONCLUSION:Recurrent TB in MDR patients in Georgia are mainly caused by relapse, raising concerns on the efficacy of the TB control program. An association between incarceration and reinfection likely reflects high levels of ongoing TB transmission in prisons, indicating the need for better TB infection control measures in these settings. Our results add to the rationale for implementing genotypic surveillance of TB more broadly to support TB control in Georgia.

BackgroundHuman tuberculosis (TB) is caused by seven phylogenetic lineages of the Mycobacterium tuberculosis complex (MTBC), Lineage 1-7. Recent advances in rapid genotyping of MTBC based on single nucleotide polymorphisms (SNP), allow for phylogenetically robust strain classification, paving the way for defining genotype-phenotype relationships in clinical settings. Such studies have revealed that, in addition to host and environmental factors, strain variation in the MTBC influences the outcome of TB infection and disease. In Tanzania, such molecular epidemiological studies of TB however are scarce in spite of a high TB burden.Methods and findingsHere we used SNP-typing to characterize a nationwide collection of 2,039 MTBC clinical isolates representative of 1.6% of all new and retreatment TB cases notified in Tanzania during 2012 and 2013. Four lineages, namely Lineage 1-4 were identified within the study population. The distribution and frequency of these lineages varied across regions but overall, Lineage 4 was the most frequent (n = 866, 42.5%), followed by Lineage 3 (n = 681, 33.4%) and 1 (n = 336, 16.5%), with Lineage 2 being the least frequent (n = 92, 4.5%). We found Lineage 2 to be independently associated with female sex (adjusted odds ratio [aOR] 2.14; 95% confidence interval [95% CI] 1.31 - 3.50, p = 0.002) and retreatment cases (aOR 1.67; 95% CI 0.95 - 2.84, p = 0. 065) in the study population. We found no associations between MTBC lineage and patient age or HIV status. Our sublineage typing based on spacer oligotyping on a subset of Lineage 1, 3 and 4 strains revealed the presence of mainly EAI, CAS and LAM families. Finally, we detected low levels of multidrug resistant isolates among a subset of 144 retreatment cases.ConclusionsThis study provides novel insights into the MTBC lineages and the possible influence of pathogen-related factors on the TB epidemic in Tanzania.

BackgroundDiabetes Mellitus (DM) is a known risk factor for tuberculosis (TB) but little is known on TB-Diabetes Mellitus (TBDM) co-morbidity in Sub-Saharan Africa.MethodsConsecutive TB cases registered at a tertiary facility in Ghana were recruited from September 2012 to April 2016 and screened for DM using random blood glucose and glycated hemoglobin (HbA1c) level. TB patients were tested for other clinical parameters including HIV co-infection and TB lesion location. Mycobacterial isolates obtained from collected sputum samples were characterized by standard methods. Associations between TBDM patients' epidemiological as well as microbiological variables were assessed.ResultsThe prevalence of DM at time of diagnosis among 2990 enrolled TB cases was 9.4% (282/2990). TBDM cases were significantly associated with weight loss, poor appetite, night sweat and fatigue (pConclusionOur findings confirms the importance of screening for diabetes during TB diagnosis and highlights the association between genetic diversity and diabetes. in Ghana.

BackgroundLimited data existed exclusively describing Mycobacterium tuberculosis lineage 3 (MTB-L3), sub-lineages, and clinical manifestations in Kampala, Uganda. This study sought to elucidate the circulating MTB-L3 sub-lineages and their corresponding clinical phenotypes.MethodA total of 141 M. tuberculosis isolates were identified as M. tuberculosis lineage 3 using Single nucleotide polymorphism (SNP) marker analysis method. To ascertain the sub-lineages/sub-strains within the M. tuberculosis lineage 3, the direct repeat (DR) loci for all the isolates was examined for sub-lineage specific signatures as described in the SITVIT2 database. The infecting sub-strains were matched with patients' clinical and demographic characteristics to identify any possible association.ResultThe data showed 3 sub-lineages circulating with CAS 1 Delhi accounting for 55% (77/141), followed by CAS 1-Kili 16% (22/141) and CAS 2/CAS 8% (12/141). Remaining isolates 21% (30/141) were unclassifiable. To explore whether the sub-lineages differ in their ability to cause increased severe disease, we used extent of lung involvement as a proxy for severe disease. Multivariable analysis showed no association between M. tuberculosis lineage 3 sub-lineages with severe disease. The risk factors associated with severe disease include having a positive smear (OR = 9.384; CI 95% = 2.603-33.835), HIV (OR = 0.316; CI 95% = 0.114-0.876), lymphadenitis (OR = 0. 171; CI 95% = 0.034-0.856) and a BCG scar (OR = 0.295; CI 95% = 0.102-0.854).ConclusionIn Kampala, Uganda, there are three sub-lineages of M. tuberculosis lineage 3 that cause disease of comparable severity with CAS-Dehli as the most prevalent. Having HIV, lymphadenitis, a BCG scar and a smear negative status is associated with reduced severe disease.

Bacteria belonging to the genus Mycobacterium are predominantly responsible for pulmonary diseases; most notably Mycobacterium tuberculosis causes granulomatous pulmonary infections. Here we describe a novel slow growing mycobacterial species isolated from respiratory samples from five patients, four with underlying pulmonary disease. The isolates were characterized by biochemical and molecular techniques, including whole genome sequencing. Biochemical characteristics generally match those of M. marinum and M. ulcerans; however, the most striking difference of the new species is its ability to grow at 37°C. The new species was found to grow in human macrophages, but not amoebae, suggesting a pathogenic rather than an environmental lifestyle. Phylogenetic analysis reveals a deep-rooting relationship to M. marinum and M. ulcerans. A complete genome sequence was obtained through combining short and long-read sequencing, providing a genome of 5.6 Mb. The genome appears to be highly intact, syntenic with that of M. marinum, with very few insertion sequences. A vast array of virulence factors includes 283 PE/PPE surface-associated proteins, making up 10% of the coding capacity, and 22 non-ribosomal peptide synthase clusters. A comparison of six clinical isolates from the five patients shows that they differ by up to two single nucleotide polymorphisms, suggesting a common source of infection. Our findings are in accordance with the recognition of a new taxonomic entity. We propose the name M. basiliense, as all isolates were found in patients from the Basel area of Switzerland.

Abstract Background Vitamin D level is inversely associated with tuberculosis (TB) and diabetes (DM). Vitamin D could be a mediator in the association between TB and DM. We examined the associations between vitamin D, TB and DM. Methods Consecutive adults with TB and sex- and age-matched volunteers were included in a case-control study in Dar es Salaam, Tanzania. Glycemia and total vitamin D (25(OH)D) were measured at enrolment and after TB treatment in cases. The association between low 25(OH)D (

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

Bacterial factors favoring the unprecedented multidrug-resistant tuberculosis (MDR-TB) epidemic in the former Soviet Union remain unclear. We utilized whole genome sequencing and Bayesian statistics to analyze the evolutionary history, temporal emergence of resistance and transmission networks of MDR Mycobacterium tuberculosis complex isolates from Karakalpakstan, Uzbekistan (2001–2006). One clade (termed Central Asian outbreak, CAO) dating back to 1974 (95% HPD 1969–1982) subsequently acquired resistance mediating mutations to eight anti-TB drugs. Introduction of standardized WHO-endorsed directly observed treatment, short-course in Karakalpakstan in 1998 likely selected for CAO-strains, comprising 75% of sampled MDR-TB isolates in 2005/2006. CAO-isolates were also identified in a published cohort from Russia (2008–2010). Similarly, the presence of mutations supposed to compensate bacterial fitness deficits was associated with transmission success and higher drug resistance rates. The genetic make-up of these MDR-strains threatens the success of both empirical and standardized MDR-TB therapies, including the newly WHO-endorsed short MDR-TB regimen in Uzbekistan.