Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

B cell genomics behind cross-neutralization of SARS-CoV-2 variants and SARS-CoV.

Monoclonal antibodies (mAbs) are a focus in vaccine and therapeutic design to counteract severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants. Here, we combined B cell sorting with single-cell VDJ and RNA sequencing (RNA-seq) and mAb structures to characterize B cell responses against SARS-CoV-2. We show that the SARS-CoV-2-specific B cell repertoire consists of transcriptionally distinct B cell populations with cells producing potently neutralizing antibodies (nAbs) localized in two clusters that resemble memory and activated B cells. Cryo-electron microscopy structures of selected nAbs from these two clusters complexed with SARS-CoV-2 spike trimers show recognition of various receptor-binding domain (RBD) epitopes. One of these mAbs, BG10-19, locks the spike trimer in a closed conformation to potently neutralize SARS-CoV-2, the recently arising mutants B.1.1.7 and B.1.351, and SARS-CoV and cross-reacts with heterologous RBDs. Together, our results characterize transcriptional differences among SARS-CoV-2-specific B cells and uncover cross-neutralizing Ab targets that will inform immunogen and therapeutic design against coronaviruses.

Pathogen Cell-to-Cell Variability Drives Heterogeneity in Host Immune Responses.

Encounters between immune cells and invading bacteria ultimately determine the course of infection. These interactions are usually measured in populations of cells, masking cell-to-cell variation that may be important for infection outcome. To characterize the gene expression variation that underlies distinct infection outcomes and monitor infection phenotypes, we developed an experimental system that combines single-cell RNA-seq with fluorescent markers. Probing the responses of individual macrophages to invading Salmonella, we find that variation between individual infected host cells is determined by the heterogeneous activity of bacterial factors in individual infecting bacteria. We illustrate how variable PhoPQ activity in the population of invading bacteria drives variable host type I IFN responses by modifying LPS in a subset of bacteria. This work demonstrates a causative link between host and bacterial variability, with cell-to-cell variation between different bacteria being sufficient to drive radically different host immune responses. This co-variation has implications for host-pathogen dynamics in vivo.

Perturbation-specific transcriptional mapping for unbiased target elucidation of antibiotics.

The rising prevalence of antibiotic resistance threatens human health. While more sophisticated strategies for antibiotic discovery are being developed, target elucidation of new chemical entities remains challenging. In the postgenomic era, expression profiling can play an important role in mechanism-of-action (MOA) prediction by reporting on the cellular response to perturbation. However, the broad application of transcriptomics has yet to fulfill its promise of transforming target elucidation due to challenges in identifying the most relevant, direct responses to target inhibition. We developed an unbiased strategy for MOA prediction, called perturbation-specific transcriptional mapping (PerSpecTM), in which large-throughput expression profiling of wild-type or hypomorphic mutants, depleted for essential targets, enables a computational strategy to address this challenge. We applied PerSpecTM to perform reference-based MOA prediction based on the principle that similar perturbations, whether chemical or genetic, will elicit similar transcriptional responses. Using this approach, we elucidated the MOAs of three molecules with activity against Pseudomonas aeruginosa by comparing their expression profiles to those of a reference set of antimicrobial compounds with known MOAs. We also show that transcriptional responses to small-molecule inhibition resemble those resulting from genetic depletion of essential targets by clustered regularly interspaced short palindromic repeats interference (CRISPRi) by PerSpecTM, demonstrating proof of concept that correlations between expression profiles of small-molecule and genetic perturbations can facilitate MOA prediction when no chemical entities exist to serve as a reference. Empowered by PerSpecTM, this work lays the foundation for an unbiased, readily scalable, systematic reference-based strategy for MOA elucidation that could transform antibiotic discovery efforts.

3-Hydroxykynurenine targets kainate receptors to promote defense against infection.

Bacterial infection involves a complex interaction between the pathogen and host where the outcome of infection is not solely determined by pathogen eradication. To identify small molecules that promote host survival by altering the host-pathogen dynamic, we conducted an in vivo chemical screen using zebrafish embryos and found that treatment with 3-hydroxykynurenine (3-HK) protects from lethal bacterial infection. 3-HK, a metabolite produced through host tryptophan metabolism, has no direct antibacterial activity but enhances host survival by restricting bacterial expansion in macrophages through a systemic mechanism that targets kainate-sensitive glutamate receptors. These findings reveal a new pathway by which tryptophan metabolism and kainate-sensitive glutamate receptors function and interact to modulate immunity, with important implications for the coordination between the immune and nervous systems in pathological conditions.

Massively multiplexed nucleic acid detection with Cas13.

The great majority of globally circulating pathogens go undetected, undermining patient care and hindering outbreak preparedness and response. To enable routine surveillance and comprehensive diagnostic applications, there is a need for detection technologies that can scale to test many samples1-3 while simultaneously testing for many pathogens4-6. Here, we develop Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (CARMEN), a platform for scalable, multiplexed pathogen detection. In the CARMEN platform, nanolitre droplets containing CRISPR-based nucleic acid detection reagents7 self-organize in a microwell array8 to pair with droplets of amplified samples, testing each sample against each CRISPR RNA (crRNA) in replicate. The combination of CARMEN and Cas13 detection (CARMEN-Cas13) enables robust testing of more than 4,500 crRNA-target pairs on a single array. Using CARMEN-Cas13, we developed a multiplexed assay that simultaneously differentiates all 169 human-associated viruses with at least 10 published genome sequences and rapidly incorporated an additional crRNA to detect the causative agent of the 2020 COVID-19 pandemic. CARMEN-Cas13 further enables comprehensive subtyping of influenza A strains and multiplexed identification of dozens of HIV drug-resistance mutations. The intrinsic multiplexing and throughput capabilities of CARMEN make it practical to scale, as miniaturization decreases reagent cost per test by more than 300-fold. Scalable, highly multiplexed CRISPR-based nucleic acid detection shifts diagnostic and surveillance efforts from targeted testing of high-priority samples to comprehensive testing of large sample sets, greatly benefiting patients and public health9-11.

Large-scale chemical-genetics yields new M. tuberculosis inhibitor classes.

New antibiotics are needed to combat rising levels of resistance, with new Mycobacterium tuberculosis (Mtb) drugs having the highest priority. However, conventional whole-cell and biochemical antibiotic screens have failed. Here we develop a strategy termed PROSPECT (primary screening of strains to prioritize expanded chemistry and targets), in which we screen compounds against pools of strains depleted of essential bacterial targets. We engineered strains that target 474 essential Mtb genes and screened pools of 100-150 strains against activity-enriched and unbiased compound libraries, probing more than 8.5 million chemical-genetic interactions. Primary screens identified over tenfold more hits than screening wild-type Mtb alone, with chemical-genetic interactions providing immediate, direct target insights. We identified over 40 compounds that target DNA gyrase, the cell wall, tryptophan, folate biosynthesis and RNA polymerase, as well as inhibitors that target EfpA. Chemical optimization yielded EfpA inhibitors with potent wild-type activity, thus demonstrating the ability of PROSPECT to yield inhibitors against targets that would have eluded conventional drug discovery.

Chemical-genetic interaction mapping links carbon metabolism and cell wall structure to tuberculosis drug efficacy.

Current chemotherapy against Mycobacterium tuberculosis (Mtb), an important human pathogen, requires a multidrug regimen lasting several months. While efforts have been made to optimize therapy by exploiting drug­drug synergies, testing new drug combinations in relevant host environments remains arduous. In particular, host environments profoundly affect the bacterial metabolic state and drug efficacy, limiting the accuracy of predictions based on in vitro assays alone. In this study, we utilized conditional Mtb knockdown mutants of essential genes as an experimentally tractable surrogate for drug treatment and probe the relationship between Mtb carbon metabolism and chemical­genetic interactions (CGIs). We examined the antitubercular drugs isoniazid, rifampicin, and moxifloxacin and found that CGIs are differentially responsive to the metabolic state, defining both environment-independent and -dependent interactions. Specifically, growth on the in vivo­relevant carbon source, cholesterol, reduced rifampicin efficacy by altering mycobacterial cell surface lipid composition. We report that a variety of perturbations in cell wall synthesis pathways restore rifampicin efficacy during growth on cholesterol, and that both environment-independent and cholesterol-dependent in vitro CGIs could be leveraged to enhance bacterial clearance in the mouse infection model. Our findings present an atlas of chemical­genetic­environmental interactions that can be used to optimize drug­drug interactions, as well as provide a framework for understanding in vitro correlates of in vivo efficacy.

A dual-plasmid CRISPR/Cas9-based method for rapid and efficient genetic disruption in Mycobacterium abscessus.

Mycobacterium abscessus is increasingly recognized for causing infections that are notoriously difficult to treat, owing to its large arsenal of intrinsic antibiotic resistance mechanisms. Tools for the genetic manipulation of the pathogen are critical for enabling a better understanding of M. abscessus biology, pathogenesis, and antibiotic resistance mechanisms. However, existing methods are largely recombination-based, which are relatively inefficient. Meanwhile, CRISPR/Cas9 has revolutionized the field of genome editing including its recent adaptation for use in mycobacteria. In this study, we report a streamlined and efficient method for rapid genetic disruptions in M. abscessus. Harnessing the CRISPR1 loci from Streptococcus thermophilus, we have developed a dual-plasmid workflow that introduces Cas9 and sgRNA cassettes in separate steps but requires no other additional factors to engineer mutations in single genes or multiple genes simultaneously or sequentially using multiple targeting sgRNAs. Importantly, the efficiency of mutant generation is several orders of magnitude higher than reported for homologous recombination-based methods. This work, thus, reports the first application of CRISPR/Cas9 for gene editing in M. abscessus and is an important tool in the arsenal for the genetic manipulation of this human pathogen. IMPORTANCE: Mycobacterium abscessus is an opportunistic pathogen of increasing clinical importance due to its poor clinical outcomes and limited treatment options. Drug discovery and development in this highly antibiotic-resistant species will require further understanding of M. abscessus biology, pathogenesis, and antibiotic resistance mechanisms. However, existing methods for facile genetic engineering are relatively inefficient. This study reports on the first application of CRISPR/Cas9 for gene editing in M. abscessus using a dual-plasmid workflow. We establish that our method is easily programmable, efficient, and versatile for genetic disruptions in M. abscessus. This is a critical advancement to facilitating targeted gene function studies in this emerging pathogen.

Fluorogenic Probes of the Mycobacterial Membrane as Reporters of Antibiotic Action.

The genus Mycobacterium includes species such as Mycobacterium tuberculosis, which can cause deadly human diseases. These bacteria have a protective cell envelope that can be remodeled to facilitate their survival in challenging conditions. Understanding how such conditions affect membrane remodeling can facilitate antibiotic discovery and treatment. To this end, we describe an optimized fluorogenic probe, N-QTF, that reports on mycolyltransferase activity, which is vital for cell division and remodeling. N-QTF is a glycolipid probe that can reveal dynamic changes in the mycobacterial cell envelope in both fast- and slow-growing mycobacterial species. Using this probe to monitor the consequences of antibiotic treatment uncovered distinct cellular phenotypes. Even antibiotics that do not directly inhibit cell envelope biosynthesis cause conspicuous phenotypes. For instance, mycobacteria exposed to the RNA polymerase inhibitor rifampicin release fluorescent extracellular vesicles (EVs). While all mycobacteria release EVs, fluorescent EVs were detected only in the presence of RIF, indicating that exposure to the drug alters EV content. Macrophages exposed to the EVs derived from RIF-treated cells released lower levels of cytokines, suggesting the EVs moderate immune responses. These data suggest that antibiotics can alter EV content to impact immunity. Our ability to see such changes in EV constituents directly results from exploiting these chemical probes.

Evaluation of expanded 2-aminobenzothiazole library as inhibitors of a model histidine kinase and virulence suppressors in Pseudomonas aeruginosa.

Bacterial resistance to antibiotics is a rapidly increasing threat to human health. New strategies to combat resistant organisms are desperately needed. One potential avenue is targeting two-component systems, which are the main bacterial signal transduction pathways used to regulate development, metabolism, virulence, and antibiotic resistance. These systems consist of a homodimeric membrane-bound sensor histidine kinase, and a cognate effector, the response regulator. Histidine kinases play an essential role in the regulation of multiple virulence mechanisms including toxin production, immune evasion, and antibiotic resistance. Targeting virulence, as opposed to development of bactericidal compounds, could reduce evolutionary pressure for acquired resistance. Additionally, compounds targeting the highly conserved catalytic and adenosine triphosphate-binding (CA) domain have the potential to impair multiple two-component systems that regulate virulence in one or more pathogens. We conducted in vitro structure-activity relationship studies of 2-aminobenzothiazole-based inhibitors designed to target the CA domain. We found that these compounds, which inhibit the model histidine kinase, HK853 from Thermotoga maritima, have anti-virulence activities inPseudomonas aeruginosa, reducing motility phenotypes and toxin production associated with the pathogenic functions of this bacterium.

Dual transcriptional analysis reveals adaptation of host and pathogen to intracellular survival of Pseudomonas aeruginosa associated with urinary tract infection.

Long-term survival of bacterial pathogens during persistent bacterial infections can be associated with antibiotic treatment failure and poses a serious public health problem. Infections caused by the Gram-negative pathogen Pseudomonas aeruginosa, which can cause both acute and chronic infections, are particularly challenging due to its high intrinsic resistance to antibiotics. The ineffectiveness of antibiotics is exacerbated when bacteria reside intracellularly within host cells where they can adopt a drug tolerant state. While the early steps of adherence and entry of P. aeruginosa into mammalian cells have been described, the subsequent fate of internalized bacteria, as well as host and bacterial molecular pathways facilitating bacterial long-term survival, are not well defined. In particular, long-term survival within bladder epithelial cells has not been demonstrated and this may have important implications for the understanding and treatment of UTIs caused by P. aeruginosa. Here, we demonstrate and characterize the intracellular survival of wild type (WT) P. aeruginosa inside bladder epithelial cells and a mutant with a disruption in the bacterial two-component regulator AlgR that is unable to survive intracellularly. Using simultaneous dual RNA-seq transcriptional profiling, we define the transcriptional response of intracellular bacteria and their corresponding invaded host cells. The bacterial transcriptional response demonstrates that WT bacteria rapidly adapt to the stress encountered in the intracellular environment in contrast to ΔalgR bacteria. Analysis of the host transcriptional response to invasion suggests that the NF-κB signaling pathway, previously shown to be required for extracellular bacterial clearance, is paradoxically also required for intracellular bacterial survival. Lastly, we demonstrate that intracellular survival is important for pathogenesis of P. aeruginosa in vivo using a model of murine urinary tract infection. We propose that the unappreciated ability of P. aeruginosa to survive intracellularly may play an important role in contributing to the chronicity and recurrence of P. aeruginosa in urinary tract infections.

Early administration of COVID-19 convalescent plasma with high titer antibody content by live viral neutralization assay is associated with modest clinical efficacy.

The efficacy of COVID-19 convalescent plasma (CCP) as a treatment for hospitalized patients with COVID-19 remains somewhat controversial; however, many studies have not evaluated CCP documented to have high neutralizing antibody titer by a highly accurate assay. To evaluate the correlation of the administration of CCP with titer determined by a live viral neutralization assay with 7- and 28-day death rates during hospitalization, a total of 23 118 patients receiving a single unit of CCP were stratified into two groups: those receiving high titer CCP (>250 50% inhibitory dilution, ID50; n = 13 636) or low titer CCP (≤250 ID50; n = 9482). Multivariable Cox regression was performed to assess risk factors. Non-intubated patients who were transfused with high titer CCP showed 1.1% and 1.7% absolute reductions in overall 7- and 28-day death rates, respectively, compared to those non-intubated patients receiving low titer CCP. No benefit of CCP was observed in intubated patients. The relative benefit of high titer CCP was confirmed in multivariable Cox regression. Administration of CCP with high titer antibody content determined by live viral neutralization assay to non-intubated patients is associated with modest clinical efficacy. Although shown to be only of modest clinical benefit, CCP may play a role in the future should viral variants develop that are not neutralized by other available therapeutics.

Defining the core essential genome of Pseudomonas aeruginosa.

Genomics offered the promise of transforming antibiotic discovery by revealing many new essential genes as good targets, but the results fell short of the promise. While numerous factors contributed to the disappointing yield, one factor was that essential genes for a bacterial species were often defined based on a single or limited number of strains grown under a single or limited number of in vitro laboratory conditions. In fact, the essentiality of a gene can depend on both the genetic background and growth condition. We thus developed a strategy for more rigorously defining the core essential genome of a bacterial species by studying many pathogen strains and growth conditions. We assessed how many strains must be examined to converge on a set of core essential genes for a species. We used transposon insertion sequencing (Tn-Seq) to define essential genes in nine strains of Pseudomonas aeruginosa on five different media and developed a statistical model, FiTnEss, to classify genes as essential versus nonessential across all strain-medium combinations. We defined a set of 321 core essential genes, representing 6.6% of the genome. We determined that analysis of four strains was typically sufficient in P. aeruginosa to converge on a set of core essential genes likely to be essential across the species across a wide range of conditions relevant to in vivo infection, and thus to represent attractive targets for novel drug discovery.

Cross-Sectional Assessment of SARS-CoV-2 Viral Load by Symptom Status in Massachusetts Congregate Living Facilities.

Transmission of coronavirus disease 2019 (COVID-19) from people without symptoms confounds societal mitigation strategies. From April to June 2020, we tested nasopharyngeal swabs by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) from 15 514 staff and 16 966 residents of nursing homes and assisted living facilities in Massachusetts. Cycle threshold (Ct) distributions were very similar between populations with (n = 739) and without (n = 2179) symptoms at the time of sampling (mean Ct, 25.7 vs 26.4; ranges 12-38). However, as local cases waned, those without symptoms shifted towards higher Ct. With such similar viral load distributions, existing testing modalities should perform comparably regardless of symptoms, contingent upon time since infection.

SARS-CoV-2 antibody persistence in COVID-19 convalescent plasma donors: Dependency on assay format and applicability to serosurveillance.

BACKGROUND: Antibody response duration following severe acute respiratory syndrome coronavirus 2 infection tends to be variable and depends on severity of disease and method of detection. STUDY DESIGN AND METHODS: COVID-19 convalescent plasma from 18 donors was collected longitudinally for a maximum of 63-129 days following resolution of symptoms. All the samples were initially screened by the Ortho total Ig test to confirm positivity and subsequently tested with seven additional direct sandwich or indirect binding assays (Ortho, Roche, Abbott, Broad Institute) directed against a variety of antigen targets (S1, receptor binding domain, and nucleocapsid [NC]), along with two neutralization assays (Broad Institute live virus PRNT and Vitalant Research Institute [VRI] Pseudovirus reporter viral particle neutralization [RVPN]). RESULTS: The direct detection assays (Ortho total Ig total and Roche total Ig) showed increasing levels of antibodies over the time period, in contrast to the indirect IgG assays that showed a decline. Neutralization assays also demonstrated declining responses; the VRI RVPN pseudovirus had a greater rate of decline than the Broad PRNT live virus assay. DISCUSSION: These data show that in addition to variable individual responses and associations with disease severity, the detection assay chosen contributes to the heterogeneous results in antibody stability over time. Depending on the scope of the research, one assay may be preferable over another. For serosurveillance studies, direct, double Ag-sandwich assays appear to be the best choice due to their stability; in particular, algorithms that include both S1- and NC-based assays can help reduce the rate of false-positivity and discriminate between natural infection and vaccine-derived seroreactivity.

Whole-organism phenotypic screening for anti-infectives promoting host health.

To date, antibiotics have been identified on the basis of their ability to kill bacteria or inhibit their growth rather than directly for their capacity to improve clinical outcomes of infected patients. Although historically successful, this approach has led to the development of an antibiotic armamentarium that suffers from a number of shortcomings, including the inevitable emergence of resistance and, in certain infections, suboptimal efficacy leading to long treatment durations, infection recurrence, or high mortality and morbidity rates despite apparent bacterial sterilization. Conventional antibiotics fail to address the complexities of in vivo bacterial physiology and virulence, as well as the role of the host underlying the complex, dynamic interactions that cause disease. New interventions are needed, aimed at host outcome rather than microbiological cure. Here we review the role of screening models for cellular and whole-organism infection, including worms, flies, zebrafish, and mice, to identify novel therapeutic strategies and discuss their future implications.

Multi-institute analysis of carbapenem resistance reveals remarkable diversity, unexplained mechanisms, and limited clonal outbreaks.

Carbapenem-resistant Enterobacteriaceae (CRE) are among the most severe threats to the antibiotic era. Multiple different species can exhibit resistance due to many different mechanisms, and many different mobile elements are capable of transferring resistance between lineages. We prospectively sampled CRE from hospitalized patients from three Boston-area hospitals, together with a collection of CRE from a single California hospital, to define the frequency and characteristics of outbreaks and determine whether there is evidence for transfer of strains within and between hospitals and the frequency with which resistance is transferred between lineages or species. We found eight species exhibiting resistance, with the majority of our sample being the sequence type 258 (ST258) lineage of Klebsiella pneumoniae There was very little evidence of extensive hospital outbreaks, but a great deal of variation in resistance mechanisms and the genomic backgrounds carrying these mechanisms. Local transmission was evident in clear phylogeographic structure between the samples from the two coasts. The most common resistance mechanisms were KPC (K. pneumoniae carbapenemases) beta-lactamases encoded by blaKPC2, blaKPC3, and blaKPC4, which were transferred between strains and species by seven distinct subgroups of the Tn4401 element. We also found evidence for previously unrecognized resistance mechanisms that produced resistance when transformed into a susceptible genomic background. The extensive variation, together with evidence of transmission beyond limited clonal outbreaks, points to multiple unsampled transmission chains throughout the continuum of care, including asymptomatic carriage and transmission of CRE. This finding suggests that to control this threat, we need an aggressive approach to surveillance and isolation.

Mutations in pmrB Confer Cross-Resistance between the LptD Inhibitor POL7080 and Colistin in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a major bacterial pathogen associated with a rising prevalence of antibiotic resistance. We evaluated the resistance mechanisms of P. aeruginosa against POL7080, a species-specific, first-in-class antibiotic in clinical trials that targets the lipopolysaccharide transport protein LptD. We isolated a series of POL7080-resistant strains with mutations in the two-component sensor gene pmrB Transcriptomic and confocal microscopy studies support a resistance mechanism shared with colistin, involving lipopolysaccharide modifications that mitigate antibiotic cell surface binding.

Systematic, multiparametric analysis of Mycobacterium tuberculosis intracellular infection offers insight into coordinated virulence.

A key to the pathogenic success of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is the capacity to survive within host macrophages. Although several factors required for this survival have been identified, a comprehensive knowledge of such factors and how they work together to manipulate the host environment to benefit bacterial survival are not well understood. To systematically identify Mtb factors required for intracellular growth, we screened an arrayed, non-redundant Mtb transposon mutant library by high-content imaging to characterize the mutant-macrophage interaction. Based on a combination of imaging features, we identified mutants impaired for intracellular survival. We then characterized the phenotype of infection with each mutant by profiling the induced macrophage cytokine response. Taking a systems-level approach to understanding the biology of identified mutants, we performed a multiparametric analysis combining pathogen and host phenotypes to predict functional relationships between mutants based on clustering. Strikingly, mutants defective in two well-known virulence factors, the ESX-1 protein secretion system and the virulence lipid phthiocerol dimycocerosate (PDIM), clustered together. Building upon the shared phenotype of loss of the macrophage type I interferon (IFN) response to infection, we found that PDIM production and export are required for coordinated secretion of ESX-1-substrates, for phagosomal permeabilization, and for downstream induction of the type I IFN response. Multiparametric clustering also identified two novel genes that are required for PDIM production and induction of the type I IFN response. Thus, multiparametric analysis combining host and pathogen infection phenotypes can be used to identify novel functional relationships between genes that play a role in infection.