An ancestral mycobacterial effector promotes dissemination of infection.

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

Mycobacterial Genetic Technologies for Probing the Host-Pathogen Microenvironment.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is one of the oldest and most successful pathogens in the world. Diverse selective pressures encountered within host cells have directed the evolution of unique phenotypic traits, resulting in the remarkable evolutionary success of this largely obligate pathogen. Despite centuries of study, the genetic repertoire utilized by Mtb to drive virulence and host immune evasion remains to be fully understood. Various genetic approaches have been and continue to be developed to tackle the challenges of functional gene annotation and validation in an intractable organism such as Mtb. In vitro and ex vivo systems remain the primary approaches to generate and confirm hypotheses that drive a general understanding of mycobacteria biology. However, it remains of great importance to characterize genetic requirements for successful infection within a host system as in vitro and ex vivo studies fail to fully replicate the complex microenvironment experienced by Mtb. In this review, we evaluate the employment of the mycobacterial genetic toolkit to probe the host-pathogen interface by surveying the current state of mycobacterial genetic studies within host systems, with a major focus on the murine model. Specifically, we discuss the different ways that these tools have been utilized to examine various aspects of infection, including bacterial survival/virulence, bacterial evasion of host immunity, and development of novel antibacterial/vaccine strategies.

Differential Requirement for IRGM Proteins during Tuberculosis Infection in Mice.

Mycobacterium tuberculosis (Mtb) is a bacterium that exclusively resides in human hosts and remains a dominant cause of morbidity and mortality among infectious diseases worldwide. Host protection against Mtb infection is dependent on the function of immunity-related GTPase clade M (IRGM) proteins. Polymorphisms in human IRGM associate with altered susceptibility to mycobacterial disease, and human IRGM promotes the delivery of Mtb into degradative autolysosomes. Among the three murine IRGM orthologs, Irgm1 has been singled out as essential for host protection during Mtb infections in cultured macrophages and in vivo. However, whether the paralogous murine Irgm genes, Irgm2 and Irgm3, play roles in host defense against Mtb or exhibit functional relationships with Irgm1 during Mtb infection remains undetermined. Here, we report that Irgm1-/- mice are indeed acutely susceptible to aerosol infection with Mtb, yet the additional deletion of the paralogous Irgm3 gene restores protective immunity to Mtb infections in Irgm1-deficient animals. Mice lacking all three Irgm genes (panIrgm-/-) are characterized by shifted lung cytokine profiles at 5 and 24 weeks postinfection, but control disease until the very late stages of the infection, when panIrgm-/- mice display increased mortality compared to wild-type mice. Collectively, our data demonstrate that disruptions in the balance between Irgm isoforms is more detrimental to the Mtb-infected host than total loss of Irgm-mediated host defense, a concept that also needs to be considered in the context of human Mtb susceptibility linked to IRGM polymorphisms.

TMEM41B is a host factor required for the replication of diverse coronaviruses including SARS-CoV-2.

Antiviral therapeutics are a front-line defense against virally induced diseases. Because viruses frequently mutate to escape direct inhibition of viral proteins, there is interest in targeting the host proteins that the virus must co-opt to complete its replication cycle. However, a detailed understanding of the interactions between the virus and the host cell is necessary in order to facilitate development of host-directed therapeutics. As a first step, we performed a genome-wide loss of function screen using the alphacoronavirus HCoV-229E to better define the interactions between coronaviruses and host factors. We report the identification and validation of an ER-resident host protein, TMEM41B, as an essential host factor for not only HCoV-229E but also genetically distinct coronaviruses including the pandemic betacoronavirus SARS-CoV-2. We show that the protein is required at an early, but post-receptor engagement, stage of the viral lifecycle. Further, mechanistic studies revealed that although the protein was not enriched at replication complexes, it likely contributes to viral replication complex formation via mobilization of cholesterol and other lipids to facilitate host membrane expansion and curvature. Continued study of TMEM41B and the development of approaches to prevent its function may lead to broad spectrum anti-coronavirus therapeutics.

Mycobacterium tuberculosis Evasion of Guanylate Binding Protein-Mediated Host Defense in Mice Requires the ESX1 Secretion System.

Cell-intrinsic immune mechanisms control intracellular pathogens that infect eukaryotes. The intracellular pathogen Mycobacterium tuberculosis (Mtb) evolved to withstand cell-autonomous immunity to cause persistent infections and disease. A potent inducer of cell-autonomous immunity is the lymphocyte-derived cytokine IFNγ. While the production of IFNγ by T cells is essential to protect against Mtb, it is not capable of fully eradicating Mtb infection. This suggests that Mtb evades a subset of IFNγ-mediated antimicrobial responses, yet what mechanisms Mtb resists remains unclear. The IFNγ-inducible Guanylate binding proteins (GBPs) are key host defense proteins able to control infections with intracellular pathogens. GBPs were previously shown to directly restrict Mycobacterium bovis BCG yet their role during Mtb infection has remained unknown. Here, we examine the importance of a cluster of five GBPs on mouse chromosome 3 in controlling Mycobacterial infection. While M. bovis BCG is directly restricted by GBPs, we find that the GBPs on chromosome 3 do not contribute to the control of Mtb replication or the associated host response to infection. The differential effects of GBPs during Mtb versus M. bovis BCG infection is at least partially explained by the absence of the ESX1 secretion system from M. bovis BCG, since Mtb mutants lacking the ESX1 secretion system become similarly susceptible to GBP-mediated immune defense. Therefore, this specific genetic interaction between the murine host and Mycobacteria reveals a novel function for the ESX1 virulence system in the evasion of GBP-mediated immunity.

The Phagocyte Oxidase Controls Tolerance to Mycobacterium tuberculosis Infection.

Protection from infectious disease relies on two distinct strategies: antimicrobial resistance directly inhibits pathogen growth, whereas infection tolerance protects from the negative impact of infection on host health. A single immune mediator can differentially contribute to these strategies in distinct contexts, confounding our understanding of protection to different pathogens. For example, the NADPH-dependent phagocyte oxidase (Phox) complex produces antimicrobial superoxide and protects from tuberculosis (TB) in humans. However, Phox-deficient mice display no sustained resistance defects to Mycobacterium tuberculosis, suggesting a more complicated role for NADPH Phox complex than strictly controlling bacterial growth. We examined the mechanisms by which Phox contributes to protection from TB and found that mice lacking the Cybb subunit of Phox suffered from a specific defect in tolerance, which was caused by unregulated Caspase-1 activation, IL-1ß production, and neutrophil influx into the lung. These studies imply that a defect in tolerance alone is sufficient to compromise immunity to M. tuberculosis and highlight a central role for Phox and Caspase-1 in regulating TB disease progression.

Statistical analysis of variability in TnSeq data across conditions using zero-inflated negative binomial regression.

BACKGROUND: Deep sequencing of transposon mutant libraries (or TnSeq) is a powerful method for probing essentiality of genomic loci under different environmental conditions. Various analytical methods have been described for identifying conditionally essential genes whose tolerance for insertions varies between two conditions. However, for large-scale experiments involving many conditions, a method is needed for identifying genes that exhibit significant variability in insertions across multiple conditions. RESULTS: In this paper, we introduce a novel statistical method for identifying genes with significant variability of insertion counts across multiple conditions based on Zero-Inflated Negative Binomial (ZINB) regression. Using likelihood ratio tests, we show that the ZINB distribution fits TnSeq data better than either ANOVA or a Negative Binomial (in a generalized linear model). We use ZINB regression to identify genes required for infection of M. tuberculosis H37Rv in C57BL/6 mice. We also use ZINB to perform a analysis of genes conditionally essential in H37Rv cultures exposed to multiple antibiotics. CONCLUSIONS: Our results show that, not only does ZINB generally identify most of the genes found by pairwise resampling (and vastly out-performs ANOVA), but it also identifies additional genes where variability is detectable only when the magnitudes of insertion counts are treated separately from local differences in saturation, as in the ZINB model.

Statistical analysis of genetic interactions in Tn-Seq data.

Tn-Seq is an experimental method for probing the functions of genes through construction of complex random transposon insertion libraries and quantification of each mutant's abundance using next-generation sequencing. An important emerging application of Tn-Seq is for identifying genetic interactions, which involves comparing Tn mutant libraries generated in different genetic backgrounds (e.g. wild-type strain versus knockout strain). Several analytical methods have been proposed for analyzing Tn-Seq data to identify genetic interactions, including estimating relative fitness ratios and fitting a generalized linear model. However, these have limitations which necessitate an improved approach. We present a hierarchical Bayesian method for identifying genetic interactions through quantifying the statistical significance of changes in enrichment. The analysis involves a four-way comparison of insertion counts across datasets to identify transposon mutants that differentially affect bacterial fitness depending on genetic background. Our approach was applied to Tn-Seq libraries made in isogenic strains of Mycobacterium tuberculosis lacking three different genes of unknown function previously shown to be necessary for optimal fitness during infection. By analyzing the libraries subjected to selection in mice, we were able to distinguish several distinct classes of genetic interactions for each target gene that shed light on their functions and roles during infection.

Red cells from ferrochelatase-deficient erythropoietic protoporphyria patients are resistant to growth of malarial parasites.

Many red cell polymorphisms are a result of selective pressure by the malarial parasite. Here, we add another red cell disease to the panoply of erythrocytic changes that give rise to resistance to malaria. Erythrocytes from individuals with erythropoietic protoporphyria (EPP) have low levels of the final enzyme in the heme biosynthetic pathway, ferrochelatase. Cells from these patients are resistant to the growth of Plasmodium falciparum malarial parasites. This phenomenon is due to the absence of ferrochelatase and not an accumulation of substrate, as demonstrated by the normal growth of P falciparum parasites in the EPP phenocopy, X-linked dominant protoporphyria, which has elevated substrate, and normal ferrochelatase levels. This observation was replicated in a mouse strain with a hypomorphic mutation in the murine ferrochelatase gene. The parasite enzyme is not essential for parasite growth as Plasmodium berghei parasites carrying a complete deletion of the ferrochelatase gene grow normally in erythrocytes, which confirms previous studies. That ferrochelatase is essential to parasite growth was confirmed by showing that inhibition of ferrochelatase using the specific competitive inhibitor, N-methylprotoporphyrin, produced a potent growth inhibition effect against cultures of P falciparum. This raises the possibility of targeting human ferrochelatase in a host-directed antimalarial strategy.

Genome-wide screen identifies host loci that modulate Mycobacterium tuberculosis fitness in immunodivergent mice.

Genetic differences among mammalian hosts and among strains of Mycobacterium tuberculosis (Mtb) are well-established determinants of tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host-pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the highly diverse BXD family of strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb-resistant C57BL/6J (B6 or B) and Mtb-susceptible DBA/2J (D2 or D) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters of "endophenotypes," each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted quantitative trait loci (QTL) mapping of these bacterial fitness endophenotypes and identified 140 host-pathogen QTL (hpQTL). We located a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes: Rv0127 (mak), Rv0359 (rip2), Rv0955 (perM), and Rv3849 (espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Genome-wide screen identifies host loci that modulate M. tuberculosis fitness in immunodivergent mice.

Genetic differences among mammalian hosts and Mycobacterium tuberculosis ( Mtb ) strains determine diverse tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host- pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the BXD family of mouse strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb -resistant C57BL/6J (B6 or B ) and Mtb -susceptible DBA/2J (D2 or D ) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters for "endophenotypes", each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted QTL mapping of these bacterial fitness endophenotypes and identified 140 h ost- p athogen quantitative trait loci ( hp QTL). We identified a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes; Rv0127 ( mak ), Rv0359 ( rip2 ), Rv0955 ( perM ), and Rv3849 ( espR ). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Glycerol contributes to tuberculosis susceptibility in male mice with type 2 diabetes.

Diabetes mellitus increases risk for tuberculosis disease and adverse outcomes. Most people with both conditions have type 2 diabetes, but it is unknown if type 1 and type 2 diabetes have identical effects on tuberculosis susceptibility. Here we show that male mice receiving a high-fat diet and streptozotocin to model type 2 diabetes, have higher mortality, more lung pathology, and higher bacterial burden following Mycobacterium tuberculosis infection compared to mice treated with streptozotocin or high-fat diet alone. Type 2 diabetes model mice have elevated plasma glycerol, which is a preferred carbon source for M. tuberculosis. Infection studies with glycerol kinase mutant M. tuberculosis reveal that glycerol utilization contributes to the susceptibility of the type 2 diabetes mice. Hyperglycemia impairs protective immunity against M. tuberculosis in both forms of diabetes, but our data show that elevated glycerol contributes to an additional adverse effect uniquely relevant to type 2 diabetes.

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice.

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.

Host protein kinases required for SARS-CoV-2 nucleocapsid phosphorylation and viral replication.

Multiple coronaviruses have emerged independently in the past 20 years that cause lethal human diseases. Although vaccine development targeting these viruses has been accelerated substantially, there remain patients requiring treatment who cannot be vaccinated or who experience breakthrough infections. Understanding the common host factors necessary for the life cycles of coronaviruses may reveal conserved therapeutic targets. Here, we used the known substrate specificities of mammalian protein kinases to deconvolute the sequence of phosphorylation events mediated by three host protein kinase families (SRPK, GSK-3, and CK1) that coordinately phosphorylate a cluster of serine and threonine residues in the viral N protein, which is required for viral replication. We also showed that loss or inhibition of SRPK1/2, which we propose initiates the N protein phosphorylation cascade, compromised the viral replication cycle. Because these phosphorylation sites are highly conserved across coronaviruses, inhibitors of these protein kinases not only may have therapeutic potential against COVID-19 but also may be broadly useful against coronavirus-mediated diseases.

These Are the Genes You're Looking For: Finding Host Resistance Genes.

Humanity's ongoing struggle with new, re-emerging and endemic infectious diseases serves as a frequent reminder of the need to understand host-pathogen interactions. Recent advances in genomics have dramatically advanced our understanding of how genetics contributes to host resistance or susceptibility to bacterial infection. Here we discuss current trends in defining host-bacterial interactions at the genome-wide level, including screens that harness CRISPR/Cas9 genome editing, natural genetic variation, proteomics, and transcriptomics. We report on the merits, limitations, and findings of these innovative screens and discuss their complementary nature. Finally, we speculate on future innovation as we continue to progress through the postgenomic era and towards deeper mechanistic insight and clinical applications.

Granulocytes act as a niche for Mycobacterium tuberculosis growth.

Granulocyte recruitment to the pulmonary compartment is a hallmark of progressive tuberculosis (TB). This process is well-documented to promote immunopathology, but can also enhance the replication of the pathogen. Both the specific granulocytes responsible for increasing mycobacterial burden and the underlying mechanisms remain obscure. We report that the known immunomodulatory effects of these cells, such as suppression of protective T-cell responses, play a limited role in altering host control of mycobacterial replication in susceptible mice. Instead, we find that the adaptive immune response preferentially restricts the burden of bacteria within monocytes and macrophages compared to granulocytes. Specifically, mycobacteria within inflammatory lesions are preferentially found within long-lived granulocytes that express intermediate levels of the Ly6G marker and low levels of antimicrobial genes. These cells progressively accumulate in the lung and correlate with bacterial load and disease severity, and the ablation of Ly6G-expressing cells lowers mycobacterial burden. These observations suggest a model in which dysregulated granulocytic influx promotes disease by creating a permissive intracellular niche for mycobacterial growth and persistence.

Distinct Bacterial Pathways Influence the Efficacy of Antibiotics against Mycobacterium tuberculosis.

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 FDA-approved drug Alectinib compromises SARS-CoV-2 nucleocapsid phosphorylation and inhibits viral infection in vitro.

While vaccines are vital for preventing COVID-19 infections, it is critical to develop new therapies to treat patients who become infected. Pharmacological targeting of a host factor required for viral replication can suppress viral spread with a low probability of viral mutation leading to resistance. In particular, host kinases are highly druggable targets and a number of conserved coronavirus proteins, notably the nucleoprotein (N), require phosphorylation for full functionality. In order to understand how targeting kinases could be used to compromise viral replication, we used a combination of phosphoproteomics and bioinformatics as well as genetic and pharmacological kinase inhibition to define the enzymes important for SARS-CoV-2 N protein phosphorylation and viral replication. From these data, we propose a model whereby SRPK1/2 initiates phosphorylation of the N protein, which primes for further phosphorylation by GSK-3a/b and CK1 to achieve extensive phosphorylation of the N protein SR-rich domain. Importantly, we were able to leverage our data to identify an FDA-approved kinase inhibitor, Alectinib, that suppresses N phosphorylation by SRPK1/2 and limits SARS-CoV-2 replication. Together, these data suggest that repurposing or developing novel host-kinase directed therapies may be an efficacious strategy to prevent or treat COVID-19 and other coronavirus-mediated diseases.

Common Variants in the Glycerol Kinase Gene Reduce Tuberculosis Drug Efficacy.

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.