Despite the administration of multiple drugs that are highly effective in vitro, tuberculosis (TB) treatment requires prolonged drug administration and is confounded by the emergence of drug-resistant strains. To understand the mechanisms that limit antibiotic efficacy, we performed a comprehensive genetic study to identify Mycobacterium tuberculosis genes that alter the rate of bacterial clearance in drug-treated mice. Several functionally distinct bacterial genes were found to alter bacterial clearance, and prominent among these was the glpK gene that encodes the glycerol-3-kinase enzyme that is necessary for glycerol catabolism. Growth on glycerol generally increased the sensitivity of M. tuberculosis to antibiotics in vitro, and glpK-deficient bacteria persisted during antibiotic treatment in vivo, particularly during exposure to pyrazinamide-containing regimens. Frameshift mutations in a hypervariable homopolymeric region of the glpK gene were found to be a specific marker of multidrug resistance in clinical M. tuberculosis isolates, and these loss-of-function alleles were also enriched in extensively drug-resistant clones. These data indicate that frequently observed variation in the glpK coding sequence produces a drug-tolerant phenotype that can reduce antibiotic efficacy and may contribute to the evolution of resistance. IMPORTANCE TB control is limited in part by the length of antibiotic treatment needed to prevent recurrent disease. To probe mechanisms underlying survival under antibiotic pressure, we performed a genetic screen for M. tuberculosis mutants with altered susceptibility to treatment using the mouse model of TB. We identified multiple genes involved in a range of functions which alter sensitivity to antibiotics. In particular, we found glycerol catabolism mutants were less susceptible to treatment and that common variation in a homopolymeric region in the glpK gene was associated with drug resistance in clinical isolates. These studies indicate that reversible high-frequency variation in carbon metabolic pathways can produce phenotypically drug-tolerant clones and have a role in the development of resistance.

The outcome of Mycobacterium tuberculosis infection and the immunological response to the bacillus Calmette-Guerin (BCG) vaccine are highly variable in humans. Deciphering the relative importance of host genetics, environment, and vaccine preparation for the efficacy of BCG has proven difficult in natural populations. We developed a model system that captures the breadth of immunological responses observed in outbred individual mice, which can be used to understand the contribution of host genetics to vaccine efficacy. This system employs a panel of highly diverse inbred mouse strains, consisting of the founders and recombinant progeny of the “Collaborative Cross” project. Unlike natural populations, the structure of this panel allows the serial evaluation of genetically identical individuals and the quantification of genotype-specific effects of interventions such as vaccination. When analyzed in the aggregate, our panel resembled natural populations in several important respects: the animals displayed a broad range of susceptibility to M. tuberculosis, differed in their immunological responses to infection, and were not durably protected by BCG vaccination. However, when analyzed at the genotype level, we found that these phenotypic differences were heritable. M. tuberculosis susceptibility varied between lines, from extreme sensitivity to progressive M. tuberculosis clearance. Similarly, only a minority of the genotypes was protected by vaccination. The efficacy of BCG was genetically separable from susceptibility to M. tuberculosis, and the lack of efficacy in the aggregate analysis was driven by nonresponsive lines that mounted a qualitatively distinct response to infection. These observations support an important role for host genetic diversity in determining BCG efficacy and provide a new resource to rationally develop more broadly efficacious vaccines.

Effective tuberculosis treatment requires at least 6 months of combination therapy. Alterations in the physiological state of the bacterium during infection are thought to reduce drug efficacy and prolong the necessary treatment period, but the nature of these adaptations remain incompletely defined. To identify specific bacterial functions that limit drug effects during infection, we employed a comprehensive genetic screening approach to identify mutants with altered susceptibility to the first-line antibiotics in the mouse model. We identified many mutations that increase the rate of bacterial clearance, suggesting new strategies for accelerating therapy. In addition, the drug-specific effects of these mutations suggested that different antibiotics are limited by distinct factors. Rifampin efficacy is inferred to be limited by cellular permeability, whereas isoniazid is preferentially affected by replication rate. Many mutations that altered bacterial clearance in the mouse model did not have an obvious effect on drug susceptibility using in vitro assays, indicating that these chemical-genetic interactions tend to be specific to the in vivo environment. This observation suggested that a wide variety of natural genetic variants could influence drug efficacy in vivo without altering behavior in standard drug-susceptibility tests. Indeed, mutations in a number of the genes identified in our study are enriched in drug-resistant clinical isolates, identifying genetic variants that may influence treatment outcome. Together, these observations suggest new avenues for improving therapy, as well as the mechanisms of genetic adaptations that limit it. IMPORTANCE Understanding how Mycobacterium tuberculosis survives during antibiotic treatment is necessary to rationally devise more effective tuberculosis (TB) chemotherapy regimens. Using genome-wide mutant fitness profiling and the mouse model of TB, we identified genes that alter antibiotic efficacy specifically in the infection environment and associated several of these genes with natural genetic variants found in drug-resistant clinical isolates. These data suggest strategies for synergistic therapies that accelerate bacterial clearance, and they identify mechanisms of adaptation to drug exposure that could influence treatment outcome.

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis.

The outcome of Mycobacterium tuberculosis infection and the immunological response to the bacillus CalmetteGuerin (BCG) vaccine are highly variable in humans. Deciphering the relative importance of host genetics, environment, and vaccine preparation for the efficacy of BCG has proven difficult in natural populations. We developed a model system that captures the breadth of immunological responses observed in outbred individual mice, which can be used to understand the contribution of host genetics to vaccine efficacy. This system employs a panel of highly diverse inbred mouse strains, consisting of the founders and recombinant progeny of the "Collaborative Cross" project. Unlike natural populations, the structure of this panel allows the serial evaluation of genetically identical individuals and the quantification of genotype-specific effects of interventions such as vaccination. When analyzed in the aggregate, our panel resembled natural populations in several important respects: the animals displayed a broad range of susceptibility to M. tuberculosis, differed in their immunological responses to infection, and were not durably protected by BCG vaccination. However, when analyzed at the genotype level, we found that these phenotypic differences were heritable. M. tuberculosis susceptibility varied between lines, from extreme sensitivity to progressive M. tuberculosis clearance. Similarly, only a minority of the genotypes was protected by vaccination. The efficacy of BCG was genetically separable from susceptibility to M. tuberculosis, and the lack of efficacy in the aggregate analysis was driven by nonresponsive lines that mounted a qualitatively distinct response to infection. These observations support an important role for host genetic diversity in determining BCG efficacy and provide a new resource to rationally develop more broadly efficacious vaccines.IMPORTANCE Tuberculosis (TB) remains an urgent global health crisis, and the efficacy of the currently used TB vaccine, M. bovis BCG, is highly variable. The design of more broadly efficacious vaccines depends on understanding the factors that limit the protection imparted by BCG. While these complex factors are difficult to disentangle in natural populations, we used a model population of mice to understand the role of host genetic composition in BCG efficacy. We found that the ability of BCG to protect mice with different genotypes was remarkably variable. The efficacy of BCG did not depend on the intrinsic susceptibility of the animal but, instead, correlated with qualitative differences in the immune responses to the pathogen. These studies suggest that host genetic polymorphism is a critical determinant of vaccine efficacy and provide a model system to develop interventions that will be useful in genetically diverse populations.