Identification of bacterial determinants of tuberculosis infection and treatment outcomes: a phenogenomic analysis of clinical strains.

BACKGROUND: Bacterial diversity could contribute to the diversity of tuberculosis infection and treatment outcomes observed clinically, but the biological basis of this association is poorly understood. The aim of this study was to identify associations between phenogenomic variation in Mycobacterium tuberculosis and tuberculosis clinical features. METHODS: We developed a high-throughput platform to define phenotype-genotype relationships in M tuberculosis clinical isolates, which we tested on a set of 158 drug-sensitive M tuberculosis strains sampled from a large tuberculosis clinical study in Ho Chi Minh City, Viet Nam. We tagged the strains with unique genetic barcodes in multiplicate, allowing us to pool the strains for in-vitro competitive fitness assays across 16 host-relevant antibiotic and metabolic conditions. Relative fitness was quantified by deep sequencing, enumerating output barcode read counts relative to input normalised values. We performed a genome-wide association study to identify phylogenetically linked and monogenic mutations associated with the in-vitro fitness phenotypes. These genetic determinants were further associated with relevant clinical outcomes (cavitary disease and treatment failure) by calculating odds ratios (ORs) with binomial logistic regressions. We also assessed the population-level transmission of strains associated with cavitary disease and treatment failure using terminal branch length analysis of the phylogenetic data. FINDINGS: M tuberculosis clinical strains had diverse growth characteristics in host-like metabolic and drug conditions. These fitness phenotypes were highly heritable, and we identified monogenic and phylogenetically linked variants associated with the fitness phenotypes. These data enabled us to define two genetic features that were associated with clinical outcomes. First, mutations in Rv1339, a phosphodiesterase, which were associated with slow growth in glycerol, were further associated with treatment failure (OR 5·34, 95% CI 1·21-23·58, p=0·027). Second, we identified a phenotypically distinct slow-growing subclade of lineage 1 strains (L1.1.1.1) that was associated with cavitary disease (OR 2·49, 1·11-5·59, p=0·027) and treatment failure (OR 4·76, 1·53-14·78, p=0·0069), and which had shorter terminal branch lengths on the phylogenetic tree, suggesting increased transmission. INTERPRETATION: Slow growth under various antibiotic and metabolic conditions served as in-vitro intermediate phenotypes underlying the association between M tuberculosis monogenic and phylogenetically linked mutations and outcomes such as cavitary disease, treatment failure, and transmission potential. These data suggest that M tuberculosis growth regulation is an adaptive advantage for bacterial success in human populations, at least in some circumstances. These data further suggest markers for the underlying bacterial processes that contribute to these clinical outcomes. FUNDING: National Health and Medical Research Council/A∗STAR, National Institutes of Allergy and Infectious Diseases, National Institute of Child Health and Human Development, and the Wellcome Trust Fellowship in Public Health and Tropical Medicine.

Rifampicin tolerance and growth fitness among isoniazid-resistant clinical Mycobacterium tuberculosis isolates: an in-vitro longitudinal study.

Antibiotic tolerance in Mycobacterium tuberculosis leads to less effective bacterial killing, poor treatment responses and resistant emergence. There is limited understanding of antibiotic tolerance in clinical isolates of M. tuberculosis. Therefore, we investigated the rifampicin tolerance of M. tuberculosis isolates, with or without pre-existing isoniazid-resistance. In-vitro rifampicin survival fractions determined by minimum duration of killing assay in isoniazid susceptible (n=119) and resistant (n=84) M. tuberculosis isolates. Rifampicin tolerance was correlated with bacterial growth, rifampicin minimum inhibitory concentrations (MICs) and isoniazid-resistant mutations. The longitudinal isoniazid-resistant isolates were analyzed for rifampicin tolerance based on collection time from patients and associated emergence of genetic variants. The median duration of rifampicin exposure reducing the M. tuberculosis surviving fraction by 90% (minimum duration of killing-MDK90) increased from 1.23 (95%CI 1.11; 1.37) and 1.31 (95%CI 1.14; 1.48) to 2.55 (95%CI 2.04; 2.97) and 1.98 (95%CI 1.69; 2.56) days, for IS and IR respectively, during 15 to 60 days of incubation respectively. Increase in MDK90 time indicated the presence of fast and slow growing tolerant sub-populations. A range of 6 log10-fold survival fraction enabled classification of tolerance as low, medium or high and revealed isoniazid-resistance association with increased tolerance with faster growth (OR=2.68 for low vs. medium, OR=4.42 for low vs. high, P-trend=0.0003). The high tolerance in longitudinal isoniazid-resistant isolates was specific to those collected during rifampicin treatment in patients and associated with bacterial genetic microvariants. Our study identifies a range of rifampicin tolerance and reveals that isoniazid resistance is associated with higher tolerance with growth fitness. Furthermore, rifampicin treatment may select isoniazid-resistant isolate microvariants with higher rifampicin tolerance, with survival potential similar to multi-drug resistant isolates. These findings suggest that isoniazid-resistant tuberculosis needs to be evaluated for rifampicin tolerance or needs further improvement in treatment regimen. It is made available under a CC-BY 4.0 International license.

High-throughput phenogenotyping of Mycobacteria tuberculosis clinical strains reveals bacterial determinants of treatment outcomes.

Background: Combatting the tuberculosis (TB) epidemic caused by Mycobacterium tuberculosis ( Mtb ) necessitates a better understanding of the factors contributing to patient clinical outcomes and transmission. While host and environmental factors have been evaluated, the impact of Mtb genetic background and phenotypic diversity is underexplored. Previous work has made associations between Mtb genetic lineages and some clinical and epidemiological features, but the bacterial traits underlying these connections are largely unknown. Methods: We developed a high-throughput functional genomics platform for defining genotype-phenotype relationships across a panel of Mtb clinical isolates. These phenotypic fitness profiles function as intermediate traits which can be linked to Mtb genetic variants and associated with clinical and epidemiological outcomes. We applied this approach to a collection of 158 Mtb strains from a study of Mtb transmission in Ho Chi Minh City, Vietnam. Mtb strains were genetically tagged in multiplicate, which allowed us to pool the strains and assess in vitro competitive fitness using deep sequencing across a set of 14 host-relevant antibiotic and metabolic conditions. Phylogenetic and monogenic associations with these intermediate traits were identified and then associated with clinical outcomes. Findings: Mtb clinical strains have a broad range of growth and drug response dynamics that can be clustered by their phylogenetic relationships. We identified novel monogenic associations with Mtb fitness in various metabolic and antibiotic conditions. Among these, we find that mutations in Rv1339 , a phosphodiesterase, which were identified through their association with slow growth in glycerol, are further associated with treatment failure. We also identify a previously uncharacterized subclade of Lineage 1 strains (L1.1.1.1) that is phenotypically distinguished by slow growth under most antibiotic and metabolic stress conditions in vitro . This clade is associated with cavitary disease, treatment failure, and demonstrates increased transmission potential. Interpretation: High-throughput phenogenotyping of Mtb clinical strains enabled bacterial intermediate trait identification that can provide a mechanistic link between Mtb genetic variation and patient clinical outcomes. Mtb strains associated with cavitary disease, treatment failure, and transmission potential display intermediate phenotypes distinguished by slow growth under various antibiotic and metabolic conditions. These data suggest that Mtb growth regulation is an adaptive advantage for host bacterial success in human populations, in at least some circumstances. These data further suggest markers for the underlying bacterial processes that govern these clinical outcomes. Funding: National Institutes of Allergy and Infectious Diseases: P01 AI132130 (SS, SMF); P01 AI143575 (XW, SMF); U19 AI142793 (QL, SMF); 5T32AI132120-03 (SS); 5T32AI132120-04 (SS); 5T32AI049928-17 (SS) Wellcome Trust Fellowship in Public Health and Tropical Medicine: 097124/Z/11/Z (NTTT) National Health and Medical Research Council (NHMRC)/A*STAR joint call: APP1056689 (SJD) The funding sources had no involvement in study methodology, data collection, analysis, and interpretation nor in the writing or submission of the manuscript. Research in context: Evidence before this study: We used different combinations of the words mycobacterium tuberculosis, tuberculosis, clinical strains, intermediate phenotypes, genetic barcoding, phenogenomics, cavitary disease, treatment failure, and transmission to search the PubMed database for all studies published up until January 20 th , 2022. We only considered English language publications, which biases our search. Previous work linking Mtb determinants to clinical or epidemiological data has made associations between bacterial lineage, or less frequently, genetic polymorphisms to in vitro or in vivo models of pathogenesis, transmission, and clinical outcomes such as cavitary disease, treatment failure, delayed culture conversion, and severity. Many of these studies focus on the global pandemic Lineage 2 and Lineage 4 Mtb strains due in part to a deletion in a polyketide synthase implicated in host-pathogen interactions. There are a number of Mtb GWAS studies that have led to novel genetic determinants of in vitro drug resistance and tolerance. Previous Mtb GWAS analyses with clinical outcomes did not experimentally test any predicted phenotypes of the clinical strains. Published laboratory-based studies of Mtb clinical strains involve relatively small numbers of strains, do not identify the genetic basis of relevant phenotypes, or link findings to the corresponding clinical outcomes. There are two recent studies of other pathogens that describe phenogenomic analyses. One study of 331 M. abscessus clinical strains performed one-by-one phenotyping to identify bacterial features associated with clearance of infection and another details a competition experiment utilizing three barcoded Plasmodium falciparum clinical isolates to assay antimalarial fitness and resistance. Added value of this study: We developed a functional genomics platform to perform high-throughput phenotyping of Mtb clinical strains. We then used these phenotypes as intermediate traits to identify novel bacterial genetic features associated with clinical outcomes. We leveraged this platform with a sample of 158 Mtb clinical strains from a cross sectional study of Mtb transmission in Ho Chi Minh City, Vietnam. To enable high-throughput phenotyping of large numbers of Mtb clinical isolates, we applied a DNA barcoding approach that has not been previously utilized for the high-throughput analysis of Mtb clinical strains. This approach allowed us to perform pooled competitive fitness assays, tracking strain fitness using deep sequencing. We measured the replicative fitness of the clinical strains in multiplicate under 14 metabolic and antibiotic stress condition. To our knowledge, this is the largest phenotypic screen of Mtb clinical isolates to date. We performed bacterial GWAS to delineate the Mtb genetic variants associated with each fitness phenotype, identifying monogenic associations with several conditions. We then defined Mtb phenotypic and genetic features associated with clinical outcomes. We find that a subclade of Mtb strains, defined by variants largely involved in fatty acid metabolic pathways, share a universal slow growth phenotype that is associated with cavitary disease, treatment failure and increased transmission potential in Vietnam. We also find that mutations in Rv1339 , a poorly characterized phosphodiesterase, also associate with slow growth in vitro and with treatment failure in patients. Implications of all the available evidence: Phenogenomic profiling demonstrates that Mtb strains exhibit distinct growth characteristics under metabolic and antibiotic stress conditions. These fitness profiles can serve as intermediate traits for GWAS and association with clinical outcomes. Intermediate phenotyping allows us to examine potential processes by which bacterial strain differences contribute to clinical outcomes. Our study identifies clinical strains with slow growth phenotypes under in vitro models of antibiotic and host-like metabolic conditions that are associated with adverse clinical outcomes. It is possible that the bacterial intermediate phenotypes we identified are directly related to the mechanisms of these outcomes, or they may serve as markers for the causal yet unidentified bacterial determinants. Via the intermediate phenotyping, we also discovered a surprising diversity in Mtb responses to the new anti-mycobacterial drugs that target central metabolic processes, which will be important in considering roll-out of these new agents. Our study and others that have identified Mtb determinants of TB clinical and epidemiological phenotypes should inform efforts to improve diagnostics and drug regimen design.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation.

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam.

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Bacterial risk factors for treatment failure and relapse among patients with isoniazid resistant tuberculosis.

BACKGROUND: Drug resistant tuberculosis (TB) is increasing in prevalence worldwide. Treatment failure and relapse is known to be high for patients with isoniazid resistant TB treated with standard first line regimens. However, risk factors for unfavourable outcomes and the optimal treatment regimen for isoniazid resistant TB are unknown. This cohort study was conducted when Vietnam used the eight month first line treatment regimen and examined risk factors for failure/relapse among patients with isoniazid resistant TB. METHODS: Between December 2008 and June 2011 2090 consecutive HIV-negative adults (≥18 years of age) with new smear positive pulmonary TB presenting at participating district TB units in Ho Chi Minh City were recruited. Participants with isoniazid resistant TB identified by Microscopic Observation Drug Susceptibility (MODS) had extended follow-up for 2 years with mycobacterial culture to test for relapse. MGIT drug susceptibility testing confirmed 239 participants with isoniazid resistant, rifampicin susceptible TB. Bacterial and demographic factors were analysed for association with treatment failure and relapse. RESULTS: Using only routine programmatic sputum smear microscopy for assessment, (months 2, 5 and 8) 30/239 (12.6%) had an unfavourable outcome by WHO criteria. Thirty-nine patients were additionally detected with unfavourable outcomes during 2 year follow up, giving a total of 69/239 (28.9%) of isoniazid (INH) resistant cases with unfavourable outcome by 2 years of follow-up. Beijing lineage was the only factor significantly associated with unfavourable outcome among INH-resistant TB cases during 2 years of follow-up. (adjusted OR = 3.16 [1.54-6.47], P = 0.002). CONCLUSION: One third of isoniazid resistant TB cases suffered failure/relapse within 2 years under the old eight month regimen. Over half of these cases were not identified by standard WHO recommended treatment monitoring. Intensified research on early identification and optimal regimens for isoniazid resistant TB is needed. Infection with Beijing genotype of TB is a significant risk factor for bacterial persistence on treatment resulting in failure/relapse within 2 years. The underlying mechanism of increased tolerance for standard drug regimens in Beijing genotype strains remains unknown.

Timing of initiation of antiretroviral therapy in human immunodeficiency virus (HIV)--associated tuberculous meningitis.

BACKGROUND: The optimal time to initiate antiretroviral therapy (ART) in human immunodeficiency virus (HIV)-associated tuberculous meningitis is unknown. METHODS: We conducted a randomized, double-blind, placebo-controlled trial of immediate versus deferred ART in patients with HIV-associated tuberculous meningitis to determine whether immediate ART reduced the risk of death. Antiretroviral drugs (zidovudine, lamivudine, and efavirenz) were started either at study entry or 2 months after randomization. All patients were treated with standard antituberculosis treatment, adjunctive dexamethasone, and prophylactic co-trimoxazole and were followed up for 12 months. We conducted intention-to-treat, per-protocol, and prespecified subgroup analyses. RESULTS: A total of 253 patients were randomized, 127 in the immediate ART group and 126 in the deferred ART group; 76 and 70 patients died within 9 months in the immediate and deferred ART groups, respectively. Immediate ART was not significantly associated with 9-month mortality (hazard ratio [HR], 1.12; 95% confidence interval [CI], .81-1.55; P = .50) or the time to new AIDS events or death (HR, 1.16; 95% CI, .87-1.55; P = .31). The percentage of patients with severe (grade 3 or 4) adverse events was high in both arms (90% in the immediate ART group and 89% in the deferred ART group; P = .84), but there were significantly more grade 4 adverse events in the immediate ART arm (102 in the immediate ART group vs 87 in the deferred ART group; P = .04). CONCLUSIONS: Immediate ART initiation does not improve outcome in patients presenting with HIV-associated tuberculous meningitis. There were significantly more grade 4 adverse events in the immediate ART arm, supporting delayed initiation of ART in HIV-associated tuberculous meningitis. Clinical Trials Registration. ISRCTN63659091.