Mutations in dnaA and a cryptic interaction site increase drug resistance in Mycobacterium tuberculosis.

Genomic dissection of antibiotic resistance in bacterial pathogens has largely focused on genetic changes conferring growth above a single critical concentration of drug. However, reduced susceptibility to antibiotics-even below this breakpoint-is associated with poor treatment outcomes in the clinic, including in tuberculosis. Clinical strains of Mycobacterium tuberculosis exhibit extensive quantitative variation in antibiotic susceptibility but the genetic basis behind this spectrum of drug susceptibility remains ill-defined. Through a genome wide association study, we show that non-synonymous mutations in dnaA, which encodes an essential and highly conserved regulator of DNA replication, are associated with drug resistance in clinical M. tuberculosis strains. We demonstrate that these dnaA mutations specifically enhance M. tuberculosis survival during isoniazid treatment via reduced expression of katG, the activator of isoniazid. To identify DnaA interactors relevant to this phenotype, we perform the first genome-wide biochemical mapping of DnaA binding sites in mycobacteria which reveals a DnaA interaction site that is the target of recurrent mutation in clinical strains. Reconstructing clinically prevalent mutations in this DnaA interaction site reproduces the phenotypes of dnaA mutants, suggesting that clinical strains of M. tuberculosis have evolved mutations in a previously uncharacterized DnaA pathway that quantitatively increases resistance to the key first-line antibiotic isoniazid. Discovering genetic mechanisms that reduce drug susceptibility and support the evolution of high-level drug resistance will guide development of biomarkers capable of prospectively identifying patients at risk of treatment failure in the clinic.

Genetic analysis of Vibrio parahaemolyticus intestinal colonization.

Vibrio parahaemolyticus is the most common cause of seafood-borne gastroenteritis worldwide and a blight on global aquaculture. This organism requires a horizontally acquired type III secretion system (T3SS2) to infect the small intestine, but knowledge of additional factors that underlie V. parahaemolyticus pathogenicity is limited. We used transposon-insertion sequencing to screen for genes that contribute to viability of V. parahaemolyticus in vitro and in the mammalian intestine. Our analysis enumerated and controlled for the host infection bottleneck, enabling robust assessment of genetic contributions to in vivo fitness. We identified genes that contribute to V. parahaemolyticus colonization of the intestine independent of known virulence mechanisms in addition to uncharacterized components of T3SS2. Our study revealed that toxR, an ancestral locus in Vibrio species, is required for V. parahaemolyticus fitness in vivo and for induction of T3SS2 gene expression. The regulatory mechanism by which V. parahaemolyticus ToxR activates expression of T3SS2 resembles Vibrio cholerae ToxR regulation of distinct virulence elements acquired via lateral gene transfer. Thus, disparate horizontally acquired virulence systems have been placed under the control of this ancestral transcription factor across independently evolved human pathogens.

A Cytosine Methyltransferase Modulates the Cell Envelope Stress Response in the Cholera Pathogen [corrected].

DNA methylation is a key epigenetic regulator in all domains of life, yet the effects of most bacterial DNA methyltransferases on cellular processes are largely undefined. Here, we used diverse techniques, including bisulfite sequencing, transcriptomics, and transposon insertion site sequencing to extensively characterize a 5-methylcytosine (5mC) methyltransferase, VchM, in the cholera pathogen, Vibrio cholerae. We have comprehensively defined VchM's DNA targets, its genetic interactions and the gene networks that it regulates. Although VchM is a relatively new component of the V. cholerae genome, it is required for optimal V. cholerae growth in vitro and during infection. Unexpectedly, the usually essential σE cell envelope stress pathway is dispensable in ∆vchM V. cholerae, likely due to its lower activation in this mutant and the capacity for VchM methylation to limit expression of some cell envelope modifying genes. Our work illuminates how an acquired DNA methyltransferase can become integrated within complex cell circuits to control critical housekeeping processes.

Correction: A Cytosine Methytransferase Modulates the Cell Envelope Stress Response in the Cholera Pathogen.

[This corrects the article DOI: 10.1371/journal.pgen.1005666.].

Peptidoglycan synthesis in Mycobacterium tuberculosis is organized into networks with varying drug susceptibility.

Peptidoglycan (PG), a complex polymer composed of saccharide chains cross-linked by short peptides, is a critical component of the bacterial cell wall. PG synthesis has been extensively studied in model organisms but remains poorly understood in mycobacteria, a genus that includes the important human pathogen Mycobacterium tuberculosis (Mtb). The principle PG synthetic enzymes have similar and, at times, overlapping functions. To determine how these are functionally organized, we carried out whole-genome transposon mutagenesis screens in Mtb strains deleted for ponA1, ponA2, and ldtB, major PG synthetic enzymes. We identified distinct factors required to sustain bacterial growth in the absence of each of these enzymes. We find that even the homologs PonA1 and PonA2 have unique sets of genetic interactions, suggesting there are distinct PG synthesis pathways in Mtb. Either PonA1 or PonA2 is required for growth of Mtb, but both genetically interact with LdtB, which has its own distinct genetic network. We further provide evidence that each interaction network is differentially susceptible to antibiotics. Thus, Mtb uses alternative pathways to produce PG, each with its own biochemical characteristics and vulnerabilities.

Phosphorylation of the Peptidoglycan Synthase PonA1 Governs the Rate of Polar Elongation in Mycobacteria.

Cell growth and division are required for the progression of bacterial infections. Most rod-shaped bacteria grow by inserting new cell wall along their mid-section. However, mycobacteria, including the human pathogen Mycobacterium tuberculosis, produce new cell wall material at their poles. How mycobacteria control this different mode of growth is incompletely understood. Here we find that PonA1, a penicillin binding protein (PBP) capable of transglycosylation and transpeptidation of cell wall peptidoglycan (PG), is a major governor of polar growth in mycobacteria. PonA1 is required for growth of Mycobacterium smegmatis and is critical for M. tuberculosis during infection. In both cases, PonA1's catalytic activities are both required for normal cell length, though loss of transglycosylase activity has a more pronounced effect than transpeptidation. Mutations that alter the amount or the activity of PonA1 result in abnormal formation of cell poles and changes in cell length. Moreover, altered PonA1 activity results in dramatic differences in antibiotic susceptibility, suggesting that a balance between the two enzymatic activities of PonA1 is critical for survival. We also find that phosphorylation of a cytoplasmic region of PonA1 is required for normal activity. Mutations in a critical phosphorylated residue affect transglycosylase activity and result in abnormal rates of cell elongation. Together, our data indicate that PonA1 is a central determinant of polar growth in mycobacteria, and its governance of cell elongation is required for robust cell fitness during both host-induced and antibiotic stress.

High-resolution genetic analysis of the requirements for horizontal transmission of the ESBL plasmid from Escherichia coli O104:H4.

Horizontal dissemination of the genes encoding extended spectrum beta-lactamases (ESBLs) via conjugative plasmids is facilitating the increasingly widespread resistance of pathogens to beta-lactam antibiotics. However, there is relatively little known about the regulatory factors and mechanisms that govern the spread of these plasmids. Here, we carried out a high-throughput, transposon insertion site sequencing analysis (TnSeq) to identify genes that enable the maintenance and transmission of pESBL, an R64 (IncI1)-related resistance plasmid that was isolated from Escherichia coli O104:H4 linked to a recent large outbreak of gastroenteritis. With a few exceptions, the majority of the genes identified as required for maintenance and transmission of pESBL matched those of their previously defined R64 counterparts. However, our analyses of the high-density transposon insertion library in pESBL also revealed two very short and linked regions that constitute a previously unrecognized regulatory system controlling spread of IncI1 plasmids. In addition, we investigated the function of the pESBL-encoded M.EcoGIX methyltransferase, which is also encoded by many other IncI1 and IncF plasmids. This enzyme proved to protect pESBL from restriction in new hosts, suggesting it aids in expanding the plasmid's host range. Collectively, our work illustrates the power of the TnSeq approach to enable rapid and comprehensive analyses of plasmid genes and sequences that facilitate the dissemination of determinants of antibiotic resistance.

ARTIST: high-resolution genome-wide assessment of fitness using transposon-insertion sequencing.

Transposon-insertion sequencing (TIS) is a powerful approach for deciphering genetic requirements for bacterial growth in different conditions, as it enables simultaneous genome-wide analysis of the fitness of thousands of mutants. However, current methods for comparative analysis of TIS data do not adjust for stochastic experimental variation between datasets and are limited to interrogation of annotated genomic elements. Here, we present ARTIST, an accessible TIS analysis pipeline for identifying essential regions that are required for growth under optimal conditions as well as conditionally essential loci that participate in survival only under specific conditions. ARTIST uses simulation-based normalization to model and compensate for experimental noise, and thereby enhances the statistical power in conditional TIS analyses. ARTIST also employs a novel adaptation of the hidden Markov model to generate statistically robust, high-resolution, annotation-independent maps of fitness-linked loci across the entire genome. Using ARTIST, we sensitively and comprehensively define Mycobacterium tuberculosis and Vibrio cholerae loci required for host infection while limiting inclusion of false positive loci. ARTIST is applicable to a broad range of organisms and will facilitate TIS-based dissection of pathways required for microbial growth and survival under a multitude of conditions.

A Genome-Wide Screen Reveals that the Vibrio cholerae Phosphoenolpyruvate Phosphotransferase System Modulates Virulence Gene Expression.

Diverse environmental stimuli and a complex network of regulatory factors are known to modulate expression of Vibrio cholerae's principal virulence factors. However, there is relatively little known about how metabolic factors impinge upon the pathogen's well-characterized cascade of transcription factors that induce expression of cholera toxin and the toxin-coregulated pilus (TCP). Here, we used a transposon insertion site (TIS) sequencing-based strategy to identify new factors required for expression of tcpA, which encodes the major subunit of TCP, the organism's chief intestinal colonization factor. Besides identifying most of the genes known to modulate tcpA expression, the screen yielded ptsI and ptsH, which encode the enzyme I (EI) and Hpr components of the V. cholerae phosphoenolpyruvate phosphotransferase system (PTS). In addition to reduced expression of TcpA, strains lacking EI, Hpr, or the associated EIIA(Glc) protein produced less cholera toxin (CT) and had a diminished capacity to colonize the infant mouse intestine. The PTS modulates virulence gene expression by regulating expression of tcpPH and aphAB, which themselves control expression of toxT, the central activator of virulence gene expression. One mechanism by which PTS promotes virulence gene expression appears to be by modulating the amounts of intracellular cyclic AMP (cAMP). Our findings reveal that the V. cholerae PTS is an additional modulator of the ToxT regulon and demonstrate the potency of loss-of-function TIS sequencing screens for defining regulatory networks.

Protein complexes and proteolytic activation of the cell wall hydrolase RipA regulate septal resolution in mycobacteria.

Peptidoglycan hydrolases are a double-edged sword. They are required for normal cell division, but when dysregulated can become autolysins lethal to bacteria. How bacteria ensure that peptidoglycan hydrolases function only in the correct spatial and temporal context remains largely unknown. Here, we demonstrate that dysregulation converts the essential mycobacterial peptidoglycan hydrolase RipA to an autolysin that compromises cellular structural integrity. We find that mycobacteria control RipA activity through two interconnected levels of regulation in vivo-protein interactions coordinate PG hydrolysis, while proteolysis is necessary for RipA enzymatic activity. Dysregulation of RipA protein complexes by treatment with a peptidoglycan synthase inhibitor leads to excessive RipA activity and impairment of correct morphology. Furthermore, expression of a RipA dominant negative mutant or of differentially processed RipA homologues reveals that RipA is produced as a zymogen, requiring proteolytic processing for activity. The amount of RipA processing differs between fast-growing and slow-growing mycobacteria and correlates with the requirement for peptidoglycan hydrolase activity in these species. Together, the complex picture of RipA regulation is a part of a growing paradigm for careful control of cell wall hydrolysis by bacteria during growth, and may represent a novel target for chemotherapy development.

Letting sleeping dos lie: does dormancy play a role in tuberculosis?

Mycobacterium tuberculosis, which causes tuberculosis, remains a major human public health threat. This is largely due to a sizeable reservoir of latently infected individuals, who may relapse into active disease decades after first acquiring the infection. Furthermore, patients have a very slow response to treatment of active disease. Latency and antibiotic tolerance are commonly taken as a proxy for dormancy, a stable nonreplicative state. However, latency is a clinical term that is solely defined by a lack of disease indicators. The actual state of the bacterium in human latency is not well understood. Here we evaluate the results of several in vitro models of dormancy and consider the applicability of various animal models for studying aspects of human latency and resistance to killing by antibiotics. Furthermore, we propose a model for the initiation of dormancy and resuscitation during infection.

High-resolution definition of the Vibrio cholerae essential gene set with hidden Markov model-based analyses of transposon-insertion sequencing data.

The coupling of high-density transposon mutagenesis to high-throughput DNA sequencing (transposon-insertion sequencing) enables simultaneous and genome-wide assessment of the contributions of individual loci to bacterial growth and survival. We have refined analysis of transposon-insertion sequencing data by normalizing for the effect of DNA replication on sequencing output and using a hidden Markov model (HMM)-based filter to exploit heretofore unappreciated information inherent in all transposon-insertion sequencing data sets. The HMM can smooth variations in read abundance and thereby reduce the effects of read noise, as well as permit fine scale mapping that is independent of genomic annotation and enable classification of loci into several functional categories (e.g. essential, domain essential or 'sick'). We generated a high-resolution map of genomic loci (encompassing both intra- and intergenic sequences) that are required or beneficial for in vitro growth of the cholera pathogen, Vibrio cholerae. This work uncovered new metabolic and physiologic requirements for V. cholerae survival, and by combining transposon-insertion sequencing and transcriptomic data sets, we also identified several novel noncoding RNA species that contribute to V. cholerae growth. Our findings suggest that HMM-based approaches will enhance extraction of biological meaning from transposon-insertion sequencing genomic data.

Cell separation in Vibrio cholerae is mediated by a single amidase whose action is modulated by two nonredundant activators.

Synthesis and hydrolysis of septal peptidoglycan (PG) are critical processes at the conclusion of cell division that enable separation of daughter cells. Cleavage of septal PG is mediated by PG amidases, hydrolytic enzymes that release peptide side chains from the glycan strand. Most gammaproteobacteria, including Escherichia coli, encode several functionally redundant periplasmic amidases. However, members of the Vibrio genus, including the enteric pathogen Vibrio cholerae, encode only a single PG amidase, AmiB. Here, we show that V. cholerae AmiB is crucial for cell division and growth. Genetic and biochemical analyses indicated that AmiB is regulated by two activators, EnvC and NlpD, at least one of which is required for AmiB's localization to the cell division site. Localization of the activators (and thus of AmiB) is dependent upon the cell division protein FtsN. These factors mediate septal PG cleavage in E. coli as well; however, their precise roles vary between the two organisms in a number of ways. Notably, even though V. cholerae EnvC and NlpD appear to be functionally redundant under most growth conditions tested, NlpD is specifically required for intestinal colonization in the infant mouse model of cholera and for V. cholerae resistance against bile salts, perhaps due to environmental regulation of AmiB or its activators. Collectively, our findings reveal that although the cellular components that enable cleavage of septal PG appear to be generally conserved between E. coli and V. cholerae, they can be combined into diverse functional regulatory networks.

Differential requirement for PBP1a and PBP1b in in vivo and in vitro fitness of Vibrio cholerae.

We investigated the roles of the Vibrio cholerae high-molecular-weight bifunctional penicillin binding proteins, PBP1a and PBP1b, in the fitness of this enteric pathogen. Using a screen for synthetic lethality, we found that the V. cholerae PBP1a and PBP1b proteins, like their Escherichia coli homologues, are each essential in the absence of the other and in the absence of the other's putative activator, the outer membrane lipoproteins LpoA and LpoB, respectively. Comparative analyses of V. cholerae mutants suggest that PBP1a/LpoA of V. cholerae play a more prominent role in generating and/or maintaining the pathogen's cell wall than PBP1b/LpoB. V. cholerae lacking PBP1b or LpoB exhibited wild-type growth under all conditions tested. In contrast, V. cholerae lacking PBP1a or LpoA exhibited growth deficiencies in minimal medium, in the presence of deoxycholate and bile, and in competition assays with wild-type cells both in vitro and in the infant mouse small intestine. PBP1a pathway mutants are particularly impaired in stationary phase, which renders them sensitive to a product(s) present in supernatants from stationary-phase wild-type cells. The marked competitive defect of the PBP1a pathway mutants in vivo was largely absent when exponential-phase cells rather than stationary-phase cells were used to inoculate suckling mice. Thus, at least for V. cholerae PBP1a pathway mutants, the growth phase of the inoculum is a key modulator of infectivity.

Targeted protein degradation in mycobacteria uncovers antibacterial effects and potentiates antibiotic efficacy.

Proteolysis-targeting chimeras (PROTACs) represent a new therapeutic modality involving selectively directing disease-causing proteins for degradation through proteolytic systems. Our ability to exploit targeted protein degradation (TPD) for antibiotic development remains nascent due to our limited understanding of which bacterial proteins are amenable to a TPD strategy. Here, we use a genetic system to model chemically-induced proximity and degradation to screen essential proteins in Mycobacterium smegmatis (Msm), a model for the human pathogen M. tuberculosis (Mtb). By integrating experimental screening of 72 protein candidates and machine learning, we find that drug-induced proximity to the bacterial ClpC1P1P2 proteolytic complex leads to the degradation of many endogenous proteins, especially those with disordered termini. Additionally, TPD of essential Msm proteins inhibits bacterial growth and potentiates the effects of existing antimicrobial compounds. Together, our results provide biological principles to select and evaluate attractive targets for future Mtb PROTAC development, as both standalone antibiotics and potentiators of existing antibiotic efficacy.

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.

CD4 T cells and CD8α+ lymphocytes are necessary for intravenous BCG-induced protection against tuberculosis in macaques.

Tuberculosis (TB) is a major cause of morbidity and mortality worldwide despite widespread intradermal (ID) BCG vaccination in newborns. We previously demonstrated that changing the route and dose of BCG vaccination from 5×105 CFU ID to 5×107 CFU intravenous (IV) resulted in prevention of infection and disease in a rigorous, highly susceptible non-human primate model of TB. Identifying the immune mechanisms of protection for IV BCG will facilitate development of more effective vaccines against TB. Here, we depleted select lymphocyte subsets in IV BCG vaccinated macaques prior to Mtb challenge to determine the cell types necessary for that protection. Depletion of CD4 T cells or all CD8α expressing lymphoycytes (both innate and adaptive) resulted in loss of protection in most macaques, concomitant with increased bacterial burdens (~4-5 log10 thoracic CFU) and dissemination of infection. In contrast, depletion of only adaptive CD8αß+ T cells did not significantly reduce protection against disease. Our results demonstrate that CD4 T cells and innate CD8α+ lymphocytes are critical for IV BCG-induced protection, supporting investigation of how eliciting these cells and their functions can improve future TB vaccines.

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.

CD8+ lymphocytes are critical for early control of tuberculosis in macaques.

The functional role of CD8+ lymphocytes in tuberculosis remains poorly understood. We depleted innate and/or adaptive CD8+ lymphocytes in macaques and showed that loss of all CD8α+ cells (using anti-CD8α antibody) significantly impaired early control of Mycobacterium tuberculosis (Mtb) infection, leading to increased granulomas, lung inflammation, and bacterial burden. Analysis of barcoded Mtb from infected macaques demonstrated that depletion of all CD8+ lymphocytes allowed increased establishment of Mtb in lungs and dissemination within lungs and to lymph nodes, while depletion of only adaptive CD8+ T cells (with anti-CD8ß antibody) worsened bacterial control in lymph nodes. Flow cytometry and single-cell RNA sequencing revealed polyfunctional cytotoxic CD8+ lymphocytes in control granulomas, while CD8-depleted animals were unexpectedly enriched in CD4 and γδ T cells adopting incomplete cytotoxic signatures. Ligand-receptor analyses identified IL-15 signaling in granulomas as a driver of cytotoxic T cells. These data support that CD8+ lymphocytes are required for early protection against Mtb and suggest polyfunctional cytotoxic responses as a vaccine target.

Interaction and modulation of two antagonistic cell wall enzymes of mycobacteria.

Bacterial cell growth and division require coordinated cell wall hydrolysis and synthesis, allowing for the removal and expansion of cell wall material. Without proper coordination, unchecked hydrolysis can result in cell lysis. How these opposing activities are simultaneously regulated is poorly understood. In Mycobacterium tuberculosis, the resuscitation-promoting factor B (RpfB), a lytic transglycosylase, interacts and synergizes with Rpf-interacting protein A (RipA), an endopeptidase, to hydrolyze peptidoglycan. However, it remains unclear what governs this synergy and how it is coordinated with cell wall synthesis. Here we identify the bifunctional peptidoglycan-synthesizing enzyme, penicillin binding protein 1 (PBP1), as a RipA-interacting protein. PBP1, like RipA, localizes both at the poles and septa of dividing cells. Depletion of the ponA1 gene, encoding PBP1 in M. smegmatis, results in a severe growth defect and abnormally shaped cells, indicating that PBP1 is necessary for viability and cell wall stability. Finally, PBP1 inhibits the synergistic hydrolysis of peptidoglycan by the RipA-RpfB complex in vitro. These data reveal a post-translational mechanism for regulating cell wall hydrolysis and synthesis through protein-protein interactions between enzymes with antagonistic functions.