Defining the core essential genome of Pseudomonas aeruginosa

A deoxyviolacein-based transposon insertion vector for pigmented tracer studies

Pigments provide a simple means to rapidly visually ascertain the quantities or presence of specific microbes in a complex community. The selection of pigment-producing colonies that are simple to differentiate from common colony phenotypes provides a high degree of certainty for the identity of pigment-tagged strains. Successful employment of pigment production is dependent on various intrinsic factors related to proper levels of gene expression and pigment production that are not always easy to predict and vary within each microbe. We have constructed a simple transposon system that incorporates the genes for the production of deoxyviolacein, a pigment produced from intracellular reserves of the amino acid tryptophan, to randomly insert these genes throughout the genome. This tool allows the user to select from many thousands of potential sites throughout a bacterial genome for an ideal location to generate the desired amount of pigment. We have applied this system to a small selection of endophytes and other model bacteria to differentiate these strains from complex communities and confirm their presence after several weeks in natural environments. We provide two examples of applications using the pigments to trace strains following introduction into plant tissues or to produce a reporter strain for extracellular nitrogen compound sensing. We recognize that this tool could have far broader utility in other applications and microbes, and describe the methodology for use by the greater scientific community.

Pangenome evaluation of gene essentiality in Streptococcus pyogenes

Bacterial species often consist of strains with variable gene content, collectively referred to as the pangenome. Variations in the genetic makeup of strains can alter bacterial physiology and fitness. To define biologically relevant genes of a genome, genome-wide transposon mutant libraries have been used to identify genes essential for survival or virulence in a given strain. Such phenotypic studies have been conducted in four different genotypes of the human pathogen Streptococcus pyogenes, yet challenges exist in comparing results across studies conducted in different genetic backgrounds and conditions. To advance genotype to phenotype inferences across different S. pyogenes strains, we built a pangenome database of 249 S. pyogenes reference genomes. We systematically re-analyzed publicly available transposon sequencing datasets from S. pyogenes using a transposon sequencing-specific analysis pipeline, Transit. Across four genetic backgrounds and nine phenotypic conditions, 355 genes were essential for survival, corresponding to ~24% of the core genome. Clusters of Orthologous Genes (COG) categories related to coenzyme and lipid transport and growth functions were overrepresented as essential. Finally, essential operons across S. pyogenes genotypes were defined, with an increased number of essential operons detected under in vivo conditions. This study provides an extendible database to which new studies can be added, and a searchable html-based resource to direct future investigations into S. pyogenes biology.IMPORTANCEStreptococcus pyogenes is a human-adapted pathogen occupying restricted ecological niches. Understanding the essentiality of genes across different strains and experimental conditions is important to direct research questions and efforts to prevent the large burden of disease caused by S. pyogenes. To this end we systematically reanalyzed transposon sequencing studies in S. pyogenes using transposon sequencing-specific methods, integrating them into an extendible meta-analysis framework. This provides a repository of gene essentiality in S. pyogenes which was used to highlight specific genes of interest and for the community to guide future phenotypic studies.

Integrating bacterial molecular genetics with chemical biology for renewed antibacterial drug discovery

The application of dyes to understanding the aetiology of infection inspired antimicrobial chemotherapy and the first wave of antibacterial drugs. The second wave of antibacterial drug discovery was driven by rapid discovery of natural products, now making up 69% of current antibacterial drugs. But now with the most prevalent natural products already discovered, ∼107 new soil-dwelling bacterial species must be screened to discover one new class of natural product. Therefore, instead of a third wave of antibacterial drug discovery, there is now a discovery bottleneck. Unlike natural products which are curated by billions of years of microbial antagonism, the vast synthetic chemical space still requires artificial curation through the therapeutics science of antibacterial drugs - a systematic understanding of how small molecules interact with bacterial physiology, effect desired phenotypes, and benefit the host. Bacterial molecular genetics can elucidate pathogen biology relevant to therapeutics development, but it can also be applied directly to understanding mechanisms and liabilities of new chemical agents with new mechanisms of action. Therefore, the next phase of antibacterial drug discovery could be enabled by integrating chemical expertise with systematic dissection of bacterial infection biology. Facing the ambitious endeavour to find new molecules from nature or new-to-nature which cure bacterial infections, the capabilities furnished by modern chemical biology and molecular genetics can be applied to prospecting for chemical modulators of new targets which circumvent prevalent resistance mechanisms.

CRISPRi functional genomics in bacteria and its application to medical and industrial research

SUMMARYFunctional genomics is the use of systematic gene perturbation approaches to determine the contributions of genes under conditions of interest. Although functional genomic strategies have been used in bacteria for decades, recent studies have taken advantage of CRISPR (clustered regularly interspaced short palindromic repeats) technologies, such as CRISPRi (CRISPR interference), that are capable of precisely modulating expression of all genes in the genome. Here, we discuss and review the use of CRISPRi and related technologies for bacterial functional genomics. We discuss the strengths and weaknesses of CRISPRi as well as design considerations for CRISPRi genetic screens. We also review examples of how CRISPRi screens have defined relevant genetic targets for medical and industrial applications. Finally, we outline a few of the many possible directions that could be pursued using CRISPR-based functional genomics in bacteria. Our view is that the most exciting screens and discoveries are yet to come.

Cofitness network connectivity determines a fuzzy essential zone in open bacterial pangenome

Most in silico evolutionary studies commonly assumed that core genes are essential for cellular function, while accessory genes are dispensable, particularly in nutrient-rich environments. However, this assumption is seldom tested genetically within the pangenome context. In this study, we conducted a robust pangenomic Tn-seq analysis of fitness genes in a nutrient-rich medium for Sinorhizobium strains with a canonical open pangenome. To evaluate the robustness of fitness category assignment, Tn-seq data for three independent mutant libraries per strain were analyzed by three methods, which indicates that the Hidden Markov Model (HMM)-based method is most robust to variations between mutant libraries and not sensitive to data size, outperforming the Bayesian and Monte Carlo simulation-based methods. Consequently, the HMM method was used to classify the fitness category. Fitness genes, categorized as essential (ES), advantage (GA), and disadvantage (GD) genes for growth, are enriched in core genes, while nonessential genes (NE) are over-represented in accessory genes. Accessory ES/GA genes showed a lower fitness effect than core ES/GA genes. Connectivity degrees in the cofitness network decrease in the order of ES, GD, and GA/NE. In addition to accessory genes, 1599 out of 3284 core genes display differential essentiality across test strains. Within the pangenome core, both shared quasi-essential (ES and GA) and strain-dependent fitness genes are enriched in similar functional categories. Our analysis demonstrates a considerable fuzzy essential zone determined by cofitness connectivity degrees in Sinorhizobium pangenome and highlights the power of the cofitness network in understanding the genetic basis of ever-increasing prokaryotic pangenome data.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

A role for the stringent response in ciprofloxacin resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a major cause of nosocomial infections and the leading cause of chronic lung infections in cystic fibrosis and chronic obstructive pulmonary disease patients. Antibiotic treatment remains challenging because P. aeruginosa is resistant to high concentrations of antibiotics and has a remarkable ability to acquire mutations conferring resistance to multiple groups of antimicrobial agents. Here we report that when P. aeruginosa is plated on ciprofloxacin (cipro) plates, the majority of cipro-resistant (ciproR) colonies observed at and after 48 h of incubation carry mutations in genes related to the Stringent Response (SR). Mutations in one of the major SR components, spoT, were present in approximately 40% of the ciproR isolates. Compared to the wild-type strain, most of these isolates had decreased growth rate, longer lag phase and altered intracellular ppGpp content. Also, 75% of all sequenced mutations were insertions and deletions, with short deletions being the most frequently occurring mutation type. We present evidence that most of the observed mutations are induced on the selective plates in a subpopulation of cells that are not instantly killed by cipro. Our results suggests that the SR may be an important contributor to antibiotic resistance acquisition in P. aeruginosa.

"Multiplexed screen identifies a Pseudomonas aeruginosa -specific small molecule targeting the outer membrane protein OprH and its interaction with LPS"

The surge of antimicrobial resistance threatens efficacy of current antibiotics, particularly against Pseudomonas aeruginosa , a highly resistant gram-negative pathogen. The asymmetric outer membrane (OM) of P. aeruginosa combined with its array of efflux pumps provide a barrier to xenobiotic accumulation, thus making antibiotic discovery challenging. We adapted PROSPECT 1 , a target-based, whole-cell screening strategy, to discover small molecule probes that kill P. aeruginosa mutants depleted for essential proteins localized at the OM. We identified BRD1401, a small molecule that has specific activity against a P. aeruginosa mutant depleted for the essential lipoprotein, OprL. Genetic and chemical biological studies identified that BRD1401 acts by targeting the OM β-barrel protein OprH to disrupt its interaction with LPS and increase membrane fluidity. Studies with BRD1401 also revealed an interaction between OprL and OprH, directly linking the OM with peptidoglycan. Thus, a whole-cell, multiplexed screen can identify species-specific chemical probes to reveal novel pathogen biology.

Siderophores promote cooperative interspecies and intraspecies cross-protection against antibiotics in vitro

The antibiotic cefiderocol hijacks iron transporters to facilitate its uptake and resists β-lactamase degradation. While effective, resistance has been detected clinically with unknown mechanisms. Here, using experimental evolution, we identified cefiderocol resistance mutations in Pseudomonas aeruginosa. Resistance was multifactorial in host-mimicking growth media, led to multidrug resistance and paid fitness costs in cefiderocol-free environments. However, kin selection drove some resistant populations to cross-protect susceptible individuals from killing by increasing pyoverdine secretion via a two-component sensor mutation. While pyochelin sensitized P. aeruginosa to cefiderocol killing, pyoverdine and the enterobacteria siderophore enterobactin displaced iron from cefiderocol, preventing uptake by susceptible cells. Among 113 P. aeruginosa intensive care unit clinical isolates, pyoverdine production directly correlated with cefiderocol tolerance, and high pyoverdine producing isolates cross-protected susceptible P. aeruginosa and other Gram-negative bacteria. These in vitro data show that antibiotic cross-protection can occur via degradation-independent mechanisms and siderophores can serve unexpected protective cooperative roles in polymicrobial communities.

Identification of complex III, NQR, and SDH as primary bioenergetic enzymes during the stationary phase of Pseudomonas aeruginosa cultured in urine-like conditions

Pseudomonas aeruginosa is a common cause of urinary tract infections by strains that are often multidrug resistant, representing a major challenge to the world's health care system. This microorganism has a highly adaptable metabolism that allows it to colonize many environments, including the urinary tract. In this work, we have characterized the metabolic strategies used by stationary phase P. aeruginosa cells cultivated in urine-like media to understand the adaptations used by this microorganism to survive and produce disease. Our proteomics results show that cells rely on the Entner-Duodoroff pathway, pentose phosphate pathway, the Krebs cycle/ glyoxylate shunt and the aerobic oxidative phosphorylation to survive in urine-like media and other conditions. A deep characterization of the oxidative phosphorylation showed that the respiratory rate of stationary phase cells is increased 3-4 times compared to cells in the logarithmic phase of growth, indicating that the aerobic metabolism plays critical roles in the stationary phase of cells grown in urine like media. Moreover, the data show that respiratory complex III, succinate dehydrogenase and the NADH dehydrogenase NQR have important functions and could be used as targets to develop new antibiotics against this bacterium.

Genomic and metabolic versatility of Pseudomonas aeruginosa contributes to its inter-kingdom transmission and survival

Pseudomonas aeruginosa is one of the most versatile bacteria with renowned pathogenicity and extensive drug resistance. The diverse habitats of this bacterium include fresh, saline and drainage waters, soil, moist surfaces, taps, showerheads, pipelines, medical implants, nematodes, insects, plants, animals, birds and humans. The arsenal of virulence factors produced by P. aeruginosa includes pyocyanin, rhamnolipids, siderophores, lytic enzymes, toxins and polysaccharides. All these virulent elements coupled with intrinsic, adaptive and acquired antibiotic resistance facilitate persistent colonization and lethal infections in different hosts. To date, treating pulmonary diseases remains complicated due to the chronic secondary infections triggered by hospital-acquired P. aeruginosa. On the contrary, this bacterium can improve plant growth by suppressing phytopathogens and insects. Notably, P. aeruginosa is one of the very few bacteria capable of trans-kingdom transmission and infection. Transfer of P. aeruginosa strains from plant materials to hospital wards, animals to humans, and humans to their pets occurs relatively often. Recently, we have identified that plant-associated P. aeruginosa strains could be pathologically similar to clinical isolates. In this review, we have highlighted the genomic and metabolic factors that facilitate the dominance of P. aeruginosa across different biological kingdoms and the varying roles of this bacterium in plant and human health.

Structural assembly of the bacterial essential interactome

The study of protein interactions in living organisms is fundamental for understanding biological processes and central metabolic pathways. Yet, our knowledge of the bacterial interactome remains limited. Here, we combined gene deletion mutant analysis with deep-learning protein folding using AlphaFold2 to predict the core bacterial essential interactome. We predicted and modeled 1402 interactions between essential proteins in bacteria and generated 146 high-accuracy models. Our analysis reveals previously unknown details about the assembly mechanisms of these complexes, highlighting the importance of specific structural features in their stability and function. Our work provides a framework for predicting the essential interactomes of bacteria and highlight the potential of deep-learning algorithms in advancing our understanding of the complex biology of living organisms. Also, the results presented here offer a promising approach to identify novel antibiotic targets.

Translational PK/PD for the Development of Novel Antibiotics-A Drug Developer's Perspective

Antibiotic development traditionally involved large Phase 3 programs, preceded by Phase 2 studies. Recognizing the high unmet medical need for new antibiotics and, in some cases, challenges to conducting large clinical trials, regulators created a streamlined clinical development pathway in which a lean clinical efficacy dataset is complemented by nonclinical data as supportive evidence of efficacy. In this context, translational Pharmacokinetic/Pharmacodynamic (PK/PD) plays a key role and is a major contributor to a "robust" nonclinical package. The classical PK/PD index approach, proven successful for established classes of antibiotics, is at the core of recent antibiotic approvals and the current antibacterial PK/PD guidelines by regulators. Nevertheless, in the case of novel antibiotics with a novel Mechanism of Action (MoA), there is no prior experience with the PK/PD index approach as the basis for translating nonclinical efficacy to clinical outcome, and additional nonclinical studies and PK/PD analyses might be considered to increase confidence. In this review, we discuss the value and limitations of the classical PK/PD approach and present potential risk mitigation activities, including the introduction of a semi-mechanism-based PK/PD modeling approach. We propose a general nonclinical PK/PD package from which drug developers might choose the studies most relevant for each individual candidate in order to build up a "robust" nonclinical PK/PD understanding.

Evaluating the In Vivo Virulence of Environmental Pseudomonas aeruginosa Using Microinjection Model of Zebrafish (Danio rerio)

(1) Background: Microinjection of zebrafish (Danio rerio) embryos offers a promising model for studying the virulence and potential environmental risks associated with Pseudomonas aeruginosa. (2) Methods: This work aimed to develop a P. aeruginosa infection model using two parallel exposition pathways on zebrafish larvae with microinjection into the yolk and the perivitelline space to simultaneously detect the invasive and cytotoxic features of the examined strains. The microinjection infection model was validated with 15 environmental and clinical strains of P. aeruginosa of various origins, antibiotic resistance profiles, genotypes and phenotypes: both exposition pathways were optimized with a series of bacterial dilutions, different drop sizes (injection volumes) and incubation periods. Besides mortality, sublethal symptoms of the treated embryos were detected and analyzed. (3) Results: According to the statistical evaluation of our results, the optimal parameters (dilution, drop size and incubation period) were determined. (4) Conclusions: The tested zebrafish embryo microinjection infection model is now ready for use to determine the in vivo virulence and ecological risk of environmental P. aeruginosa.

Machine learning identification of Pseudomonas aeruginosa strains from colony image data

When grown on agar surfaces, microbes can produce distinct multicellular spatial structures called colonies, which contain characteristic sizes, shapes, edges, textures, and degrees of opacity and color. For over one hundred years, researchers have used these morphology cues to classify bacteria and guide more targeted treatment of pathogens. Advances in genome sequencing technology have revolutionized our ability to classify bacterial isolates and while genomic methods are in the ascendancy, morphological characterization of bacterial species has made a resurgence due to increased computing capacities and widespread application of machine learning tools. In this paper, we revisit the topic of colony morphotype on the within-species scale and apply concepts from image processing, computer vision, and deep learning to a dataset of 69 environmental and clinical Pseudomonas aeruginosa strains. We find that colony morphology and complexity under common laboratory conditions is a robust, repeatable phenotype on the level of individual strains, and therefore forms a potential basis for strain classification. We then use a deep convolutional neural network approach with a combination of data augmentation and transfer learning to overcome the typical data starvation problem in biological applications of deep learning. Using a train/validation/test split, our results achieve an average validation accuracy of 92.9% and an average test accuracy of 90.7% for the classification of individual strains. These results indicate that bacterial strains have characteristic visual 'fingerprints' that can serve as the basis of classification on a sub-species level. Our work illustrates the potential of image-based classification of bacterial pathogens and highlights the potential to use similar approaches to predict medically relevant strain characteristics like antibiotic resistance and virulence from colony data.

An expanded CRISPR-Cas9-assisted recombineering toolkit for engineering genetically intractable Pseudomonas aeruginosa isolates

Much of our current understanding of microbiology is based on the application of genetic engineering procedures. Since their inception (more than 30 years ago), methods based largely on allelic exchange and two-step selection processes have become a cornerstone of contemporary bacterial genetics. While these tools are established for adapted laboratory strains, they have limited applicability in clinical or environmental isolates displaying a large and unknown genetic repertoire that are recalcitrant to genetic modifications. Hence, new tools allowing genetic engineering of intractable bacteria must be developed to gain a comprehensive understanding of them in the context of their biological niche. Herein, we present a method for precise, efficient and rapid engineering of the opportunistic pathogen Pseudomonas aeruginosa. This procedure relies on recombination of short single-stranded DNA facilitated by targeted double-strand DNA breaks mediated by a synthetic Cas9 coupled with the efficient Ssr recombinase. Possible applications include introducing single-nucleotide polymorphisms, short or long deletions, and short DNA insertions using synthetic single-stranded DNA templates, drastically reducing the need of PCR and cloning steps. Our toolkit is encoded on two plasmids, harboring an array of different antibiotic resistance cassettes; hence, this approach can be successfully applied to isolates displaying natural antibiotic resistances. Overall, this toolkit substantially reduces the time required to introduce a range of genetic manipulations to a minimum of five experimental days, and enables a variety of research and biotechnological applications in both laboratory strains and difficult-to-manipulate P. aeruginosa isolates.

Microbial Primer: Transposon directed insertion site sequencing (TraDIS): A high throughput method for linking genotype to phenotype

Genetic screens are a key tool for linking phenotype and genotype. Transposon mutagenesis was one of the first genetic methodologies to associate genetic loci with phenotypes. The advent of next-generation sequencing transformed the use of this technique allowing rapid interrogation of whole genomes for genes that correlate with phenotype. One method is transposon directed insertion-site sequencing (TraDIS). Here we describe the method, recent developments in technology, and the advantages and disadvantages of this method compared to other genetic screening tools.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multi-locus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

In vivo gene expression profile of Haemophilus influenzae during human pneumonia

Haemophilus influenzae is a major cause of community-acquired pneumonia. While studied extensively in various laboratory models, less is known about the cell function while inside the human lung. We present the first analysis of the global gene expression of H. influenzae while the bacteria are in the lung during pneumonia (in vivo conditions) and contrast it with bacterial isolates that have been cultured under standard laboratory conditions (in vitro conditions). Patients with pneumonia were recruited from emergency departments and intensive care units during 2018-2020 (n = 102). Lower respiratory samples were collected for bacterial culture and RNA extraction. Patient samples with H. influenzae (n = 8) and colonies from bacterial cultures (n = 6) underwent RNA sequencing. The reads were then pseudo-aligned to core and pan genomes created from 15 reference strains. While bacteria cultured in vitro clustered tightly by principal component analysis of core genome (n = 1067) gene expression, bacteria in the patient samples had more diverse transcriptomic signatures and did not group with their lab-cultured counterparts. In total, 328 core genes were significantly differentially expressed between in vitro and in vivo conditions. The most highly upregulated genes in vivo included tbpA and fbpA, which are involved in the acquisition of iron from transferrin, and the stress response gene msrAB. The biosynthesis of nucleotides/purines and molybdopterin-scavenging processes were also significantly enriched in vivo. In contrast, major metabolic pathways and iron-sequestering genes were downregulated under this condition. In conclusion, extensive transcriptomic differences were found between bacteria while in the human lung and bacteria that were cultured in vitro. IMPORTANCE The human-specific pathogen Haemophilus influenzae is generally not well suited for studying in animal models, and most laboratory models are unlikely to approximate the diverse environments encountered by bacteria in the human airways accurately. Thus, we have examined the global gene expression of H. influenzae during pneumonia. Extensive differences in the global gene expression profiles were found in H. influenzae while in the human lung compared to bacteria that were grown in the laboratory. In contrast, the gene expression profiles of isolates collected from different patients were found to cluster together when grown under the same laboratory conditions. Interesting observations were made of how H. influenzae acquires and uses iron and molybdate, endures oxidative stress, and regulates central metabolism while in the lung. Our results indicate important processes during infection and can guide future research on genes and pathways that are relevant in the pathogenesis of H. influenzae pneumonia.

Reconstruction of the last bacterial common ancestor from 183 pangenomes reveals a versatile ancient core genome

BACKGROUND

Cumulative sequencing efforts have yielded enough genomes to construct pangenomes for dozens of bacterial species and elucidate intraspecies gene conservation. Given the diversity of organisms for which this is achievable, similar analyses for ancestral species are feasible through the integration of pangenomics and phylogenetics, promising deeper insights into the nature of ancient life.

RESULTS

We construct pangenomes for 183 bacterial species from 54,085 genomes and identify their core genomes using a novel statistical model to estimate genome-specific error rates and underlying gene frequencies. The core genomes are then integrated into a phylogenetic tree to reconstruct the core genome of the last bacterial common ancestor (LBCA), yielding three main results: First, the gene content of modern and ancestral core genomes are diverse at the level of individual genes but are similarly distributed by functional category and share several poorly characterized genes. Second, the LBCA core genome is distinct from any individual modern core genome but has many fundamental biological systems intact, especially those involving translation machinery and biosynthetic pathways to all major nucleotides and amino acids. Third, despite this metabolic versatility, the LBCA core genome likely requires additional non-core genes for viability, based on comparisons with the minimal organism, JCVI-Syn3A.

CONCLUSIONS

These results suggest that many cellular systems commonly conserved in modern bacteria were not just present in ancient bacteria but were nearly immutable with respect to short-term intraspecies variation. Extending this analysis to other domains of life will likely provide similar insights into more distant ancestral species.

Analysis of the Plasmid-Based ts-Mutant ΔfabA/pTS-fabA Reveals Its Lethality under Aerobic Growth Conditions That Is Suppressed by Mild Overexpression of desA at a Restrictive Temperature in Pseudomonas aeruginosa

It is uncertain whether PA1610|fabA is essential or dispensable for growth on LB-agar plates under aerobic conditions in Pseudomonas aeruginosa PAO1. To examine its essentiality, we disrupted fabA in the presence of a native promoter-controlled complementary copy on ts-plasmid. In this analysis, we showed that the plasmid-based ts-mutant ΔfabA/pTS-fabA failed to grow at a restrictive temperature, consistent with the observation by Hoang and Schweizer (T. T. Hoang, H. P. Schweizer, J Bacteriol 179:5326-5332, 1997, https://doi.org/10.1128/jb.179.17.5326-5332.1997), and expanded on this by showing that ΔfabA exhibited curved cell morphology. On the other hand, strong induction of fabA-OE or PA3645|fabZ-OE impeded the growth of cells displaying oval morphology. Suppressor analysis revealed a mutant sup gene that suppressed a growth defect but not cell morphology of ΔfabA. Genome resequencing and transcriptomic profiling of sup identified PA0286|desA, whose promoter carried a single-nucleotide polymorphism (SNP), and transcription was significantly upregulated (level increase of >2-fold, P < 0.05). By integration of the SNP-bearing promoter-controlled desA gene into the chromosome of ΔfabA/pTS-fabA, we showed that the SNP is sufficient for ΔfabA to phenocopy the sup mutant. Furthermore, mild induction of the araC-PBAD-controlled desA gene but not desB rescued ΔfabA. These results validated that mild overexpression of desA fully suppressed the lethality but not the curved cell morphology of ΔfabA. Similarly, Zhu et al. (Zhu K, Choi K-H, Schweizer HP, Rock CO, Zhang Y-M, Mol Microbiol 60:260-273, 2006, https://doi.org/10.1111/j.1365-2958.2006.05088.x) showed that multicopy desA partially alleviated the slow growth phenotype of ΔfabA, the difference in which was that ΔfabA was viable. Taken together, our results demonstrate that fabA is essential for aerobic growth. We propose that the plasmid-based ts-allele is useful for exploring the genetic suppression interaction of essential genes of interest in P. aeruginosa. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen whose multidrug resistance demands new drug development. Fatty acids are essential for viability, and essential genes are ideal drug targets. However, the growth defect of essential gene mutants can be suppressed. Suppressors tend to be accumulated during the construction of essential gene deletion mutants, hampering the genetic analysis. To circumvent this issue, we constructed a deletion allele of fabA in the presence of a native promoter-controlled complementary copy in the ts-plasmid. In this analysis, we showed that ΔfabA/pTS-fabA failed to grow at a restrictive temperature, supporting its essentiality. Suppressor analysis revealed desA, whose promoter carried a SNP and whose transcription was upregulated. We validated that both the SNP-bearing promoter-controlled and regulable PBAD promoter-controlled desA suppressed the lethality of ΔfabA. Together, our results demonstrate that fabA is essential for aerobic growth. We propose that plasmid-based ts-alleles are suitable for genetic analysis of essential genes of interest.

Core defense hotspots within Pseudomonas aeruginosa are a consistent and rich source of anti-phage defense systems

Bacteria use a diverse arsenal of anti-phage immune systems, including CRISPR-Cas and restriction enzymes. Recent advances in anti-phage system discovery and annotation tools have unearthed many unique systems, often encoded in horizontally transferred defense islands, which can be horizontally transferred. Here, we developed Hidden Markov Models (HMMs) for defense systems and queried microbial genomes on the NCBI database. Out of the 30 species with >200 completely sequenced genomes, our analysis found Pseudomonas aeruginosa exhibits the greatest diversity of anti-phage systems, as measured by Shannon entropy. Using network analysis to identify the common neighbors of anti-phage systems, we identified two core defense hotspot loci (cDHS1 and cDHS2). cDHS1 is up to 224 kb (median: 26 kb) with varied arrangements of more than 30 distinct immune systems across isolates, while cDHS2 has 24 distinct systems (median: 6 kb). Both cDHS regions are occupied in a majority of P. aeruginosa isolates. Most cDHS genes are of unknown function potentially representing new anti-phage systems, which we validated by identifying a novel anti-phage system (Shango) commonly encoded in cDHS1. Identifying core genes flanking immune islands could simplify immune system discovery and may represent popular landing spots for diverse MGEs carrying anti-phage systems.

Outsmarting Pathogens with Antibody Engineering

There is growing interest in identifying antibodies that protect against infectious diseases, especially for high-risk individuals and pathogens for which no vaccine is yet available. However, pathogens that manifest as opportunistic or latent infections express complex arrays of virulence-associated proteins and are adept at avoiding immune responses. Some pathogens have developed strategies to selectively destroy antibodies, whereas others create decoy epitopes that trick the host immune system into generating antibodies that are at best nonprotective and at worst enhance pathogenesis. Antibody engineering strategies can thwart these efforts by accessing conserved neutralizing epitopes, generating Fc domains that resist capture or degradation and even accessing pathogens hidden inside cells. Design of pathogen-resistant antibodies can enhance protection and guide development of vaccine immunogens against these complex pathogens. Here, we discuss general strategies for design of antibodies resistant to specific pathogen defense mechanisms.

Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches

The Pseudomonas aeruginosa genome can change to adapt to different ecological niches. We compared four genomes from a Mexican hospital and 59 genomes from GenBank from different niches, such as urine, sputum, and environmental. The ST analysis showed that high-risk STs (ST235, ST773, and ST27) were present in the genomes of the three niches from GenBank, and the STs of Mexican genomes (ST167, ST2731, and ST549) differed from the GenBank genomes. Phylogenetic analysis showed that the genomes were clustering according to their ST and not their niche. When analyzing the genomic content, we observed that environmental genomes had genes involved in adapting to the environment not found in the clinics and that their mechanisms of resistance were mutations in antibiotic resistance-related genes. In contrast, clinical genomes from GenBank had resistance genes, in mobile/mobilizable genetic elements in the chromosome, except for the Mexican genomes that carried them mostly in plasmids. This was related to the presence of CRISPR-Cas and anti-CRISPR; however, Mexican strains only had plasmids and CRISPR-Cas. bla_OXA-488 (a variant of bla_OXA50) with higher activity against carbapenems was more prevalent in sputum genomes. The virulome analysis showed that exoS was most prevalent in the genomes of urinary samples and exoU and pldA in sputum samples. This study provides evidence regarding the genetic variability among P. aeruginosa isolated from different niches.

Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Pseudomonas aeruginosa is an opportunistic pathogen that causes difficult-to-treat infections. Two well-studied divergent P. aeruginosa strain types, PAO1 and PA14, have significant genomic heterogeneity, including diverse accessory genes present in only some strains. Genome content comparisons find core genes that are conserved across both PAO1 and PA14 strains and accessory genes that are present in only a subset of PAO1 and PA14 strains. Here, we use recently assembled transcriptome compendia of publicly available P. aeruginosa RNA sequencing (RNA-seq) samples to create two smaller compendia consisting of only strain PAO1 or strain PA14 samples with each aligned to their cognate reference genome. We confirmed strain annotations and identified other samples for inclusion by assessing each sample's median expression of PAO1-only or PA14-only accessory genes. We then compared the patterns of core gene expression in each strain. To do so, we developed a method by which we analyzed genes in terms of which genes showed similar expression patterns across strain types. We found that some core genes had consistent correlated expression patterns across both compendia, while others were less stable in an interstrain comparison. For each accessory gene, we also determined core genes with correlated expression patterns. We found that stable core genes had fewer coexpressed neighbors that were accessory genes. Overall, this approach for analyzing expression patterns across strain types can be extended to other groups of genes, like phage genes, or applied for analyzing patterns beyond groups of strains, such as samples with different traits, to reveal a deeper understanding of regulation. IMPORTANCE Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. Understanding how these different strain-level traits manifest is important for identifying targets that regulate different traits of interest. With the availability of thousands of PAO1 and PA14 samples, we created two strain-specific RNA-seq compendia where each one contains hundreds of samples from PAO1 or PA14 strains and used them to compare the expression patterns of core genes that are conserved in both strain types and to determine which core genes have expression patterns that are similar to those of accessory genes that are unique to one strain or the other using an approach that we developed. We found a subset of core genes with different transcriptional patterns across PAO1 and PA14 strains and identified those core genes with expression patterns similar to those of strain-specific accessory genes.

Pan-Genomics of Escherichia albertii for Antibiotic Resistance Profiling in Different Genome Fractions and Natural Product Mediated Intervention: In Silico Approach

Escherichia albertii is an emerging, enteric pathogen of significance. It was first isolated in 2003 from a pediatric diarrheal sample from Bangladesh. In this study, a comprehensive in silico strategy was followed to first list out antibiotic-resistant genes from core, accessory and unique genome fractions of 95 available genomes of E. albertii. Then, 56 drug targets were identified from the core essential genome. Finally, ZipA, an essential cell division protein that stabilizes the FtsZ protofilaments by cross-linking them and serves as a cytoplasmic membrane anchor for the Z ring, was selected for further downstream processing. It was computationally modeled using a threading approach, followed by virtual screening of two phytochemical libraries, Ayurvedic (n = 2103 compounds) and Traditional Chinese Medicine (n = 36,043 compounds). ADMET profiling, followed by PBPK modeling in the central body compartment, in a population of 250 non-diseased, 250 cirrhotic and 250 renally impaired people was attempted. ZINC85624912 from Chinese medicinal library showed the highest bioavailability and plasma retention. This is the first attempt to simulate the fate of natural products in the body through PBPK. Dynamics simulation of 20 ns for the top three compounds from both libraries was also performed to validate the stability of the compounds. The obtained information from the current study could aid wet-lab scientists to work on the scaffold of screened drug-like compounds from natural resources and could be useful in our quest for therapy against antibiotic-resistant E. albertii.

Diverse Partners of the Partitioning ParB Protein in Pseudomonas aeruginosa

In the majority of bacterial species, the tripartite ParAB-parS system, composed of an ATPase (ParA), a DNA-binding protein (ParB), and its target parS sequence(s), assists in the chromosome partitioning. ParB forms large nucleoprotein complexes at parS(s), located in the vicinity of origin of chromosomal replication (oriC), which after replication are subsequently positioned by ParA in cell poles. Remarkably, ParA and ParB participate not only in the chromosome segregation but through interactions with various cellular partners they are also involved in other cell cycle-related processes, in a species-specific manner. In this work, we characterized Pseudomonas aeruginosa ParB interactions with the cognate ParA, showing that the N-terminal motif of ParB is required for these interactions, and demonstrated that ParAB-parS-mediated rapid segregation of newly replicated ori domains prevented structural maintenance of chromosome (SMC)-mediated cohesion of sister chromosomes. Furthermore, using proteome-wide techniques, we have identified other ParB partners in P. aeruginosa, which encompass a number of proteins, including the nucleoid-associated proteins NdpA(PA3849) and NdpA2, MinE (PA3245) of Min system, and transcriptional regulators and various enzymes, e.g., CTP synthetase (PA3637). Among them are also NTPases PA4465, PA5028, PA3481, and FleN (PA1454), three of them displaying polar localization in bacterial cells. Overall, this work presents the spectrum of P. aeruginosa ParB partners and implicates the role of this protein in the cross-talk between chromosome segregation and other cellular processes. IMPORTANCE In Pseudomonas aeruginosa, a Gram-negative pathogen causing life-threatening infections in immunocompromised patients, the ParAB-parS system is involved in the precise separation of newly replicated bacterial chromosomes. In this work, we identified and characterized proteins interacting with partitioning protein ParB. We mapped the domain of interactions with its cognate ParA partner and showed that ParB-ParA interactions are crucial for the chromosome segregation and for proper SMC action on DNA. We also demonstrated ParB interactions with other DNA binding proteins, metabolic enzymes, and NTPases displaying polar localization in the cells. Overall, this study uncovers novel players cooperating with the chromosome partition system in P. aeruginosa, supporting its important regulatory role in the bacterial cell cycle.

Genome-scale model of Pseudomonas aeruginosa metabolism unveils virulence and drug potentiation

Pseudomonas aeruginosa is one of the leading causes of hospital-acquired infections. To decipher the metabolic mechanisms associated with virulence and antibiotic resistance, we have developed an updated genome-scale model (GEM) of P. aeruginosa. The model (iSD1509) is an extensively curated, three-compartment, and mass-and-charge balanced BiGG model containing 1509 genes, the largest gene content for any P. aeruginosa GEM to date. It is the most accurate with prediction accuracies as high as 92.4% (gene essentiality) and 93.5% (substrate utilization). In iSD1509, we newly added a recently discovered pathway for ubiquinone-9 biosynthesis which is required for anaerobic growth. We used a modified iSD1509 to demonstrate the role of virulence factor (phenazines) in the pathogen survival within biofilm/oxygen-limited condition. Further, the model can mechanistically explain the overproduction of a drug susceptibility biomarker in the P. aeruginosa mutants. Finally, we use iSD1509 to demonstrate the drug potentiation by metabolite supplementation, and elucidate the mechanisms behind the phenotype, which agree with experimental results.

Bacterial envelope stress responses: Essential adaptors and attractive targets

Millions of deaths a year across the globe are linked to antimicrobial resistant infections. The need to develop new treatments and repurpose of existing antibiotics grows more pressing as the growing antimicrobial resistance pandemic advances. In this review article, we propose that envelope stress responses, the signaling pathways bacteria use to recognize and adapt to damage to the most vulnerable outer compartments of the microbial cell, are attractive targets. Envelope stress responses (ESRs) support colonization and infection by responding to a plethora of toxic envelope stresses encountered throughout the body; they have been co-opted into virulence networks where they work like global positioning systems to coordinate adhesion, invasion, microbial warfare, and biofilm formation. We highlight progress in the development of therapeutic strategies that target ESR signaling proteins and adaptive networks and posit that further characterization of the molecular mechanisms governing these essential niche adaptation machineries will be important for sparking new therapeutic approaches aimed at short-circuiting bacterial adaptation.

Nutrient stress is a target for new antibiotics

Novel approaches are required to address the looming threat of pan-resistant Gram-negative pathogens and forestall the rise of untreatable infections. Unconventional targets that are uniquely important during infection and tractable to high-throughput drug discovery methods hold high potential for innovation in antibiotic discovery programs. In this context, inhibitors of bacterial nutrient stress are particularly exciting candidates for future antibiotic development. Amino acid, nucleotide, and vitamin biosynthesis pathways are critical for bacterial growth in nutrient-limiting conditions in the laboratory and the host. Although historically dismissed as dispensable for pathogens, a wealth of transposon mutagenesis and single-mutant studies have emerged which demonstrate that several such pathways are critical for infection. Indeed, high-throughput screens of diverse synthetic compounds and natural products have uncovered inhibitors of nutrient biosynthesis. Herein, we review bacterial nutrient biosynthesis and its role during host infection. Further, we explore screening platforms developed to search for inhibitors of these targets and highlight successes among these. Finally, we feature important and sometimes surprising connections between bacterial nutrient biosynthesis, antibiotic activity, and antibiotic resistance.

Genome-wide screen in human plasma identifies multifaceted complement evasion of Pseudomonas aeruginosa

Pseudomonas aeruginosa, an opportunistic Gram-negative pathogen, is a leading cause of bacteremia with a high mortality rate. We recently reported that P. aeruginosa forms a persister-like sub-population of evaders in human plasma. Here, using a gain-of-function transposon sequencing (Tn-seq) screen in plasma, we identified and validated previously unknown factors affecting bacterial persistence in plasma. Among them, we identified a small periplasmic protein, named SrgA, whose expression leads to up to a 100-fold increase in resistance to killing. Additionally, mutants in pur and bio genes displayed higher tolerance and persistence, respectively. Analysis of several steps of the complement cascade and exposure to an outer-membrane-impermeable drug, nisin, suggested that the mutants impede membrane attack complex (MAC) activity per se. Electron microscopy combined with energy-dispersive X-ray spectroscopy (EDX) revealed the formation of polyphosphate (polyP) granules upon incubation in plasma of different size in purD and wild-type strains, implying the bacterial response to a stress signal. Indeed, inactivation of ppk genes encoding polyP-generating enzymes lead to significant elimination of persisting bacteria from plasma. Through this study, we shed light on a complex P. aeruginosa response to the plasma conditions and discovered the multifactorial origin of bacterial resilience to MAC-induced killing.

The Pseudomonas aeruginosa RpoH (σ32) Regulon and Its Role in Essential Cellular Functions, Starvation Survival, and Antibiotic Tolerance

The bacterial heat-shock response is regulated by the alternative sigma factor, σ³² (RpoH), which responds to misfolded protein stress and directs the RNA polymerase to the promoters for genes required for protein refolding or degradation. In P. aeruginosa, RpoH is essential for viability under laboratory growth conditions. Here, we used a transcriptomics approach to identify the genes of the RpoH regulon, including RpoH-regulated genes that are essential for P. aeruginosa. We placed the rpoH gene under control of the arabinose-inducible P_BAD promoter, then deleted the chromosomal rpoH allele. This allowed transcriptomic analysis of the RpoH (σ³²) regulon following a short up-shift in the cellular concentration of RpoH by arabinose addition, in the absence of a sudden change in temperature. The P. aeruginosa ∆rpoH (P_BAD-rpoH) strain grew in the absence of arabinose, indicating that some rpoH expression occurred without arabinose induction. When arabinose was added, the rpoH mRNA abundance of P. aeruginosa ∆rpoH (P_BAD-rpoH) measured by RT-qPCR increased five-fold within 15 min of arabinose addition. Transcriptome results showed that P. aeruginosa genes required for protein repair or degradation are induced by increased RpoH levels, and that many genes essential for P. aeruginosa growth are induced by RpoH. Other stress response genes induced by RpoH are involved in damaged nucleic acid repair and in amino acid metabolism. Annotation of the hypothetical proteins under RpoH control included proteins that may play a role in antibiotic resistances and in non-ribosomal peptide synthesis. Phenotypic analysis of P. aeruginosa ∆rpoH (P_BAD-rpoH) showed that it is impaired in its ability to survive during starvation compared to the wild-type strain. P. aeruginosa ∆rpoH (P_BAD-rpoH) also had increased sensitivity to aminoglycoside antibiotics, but not to other classes of antibiotics, whether cultured planktonically or in biofilms. The enhanced aminoglycoside sensitivity of the mutant strain may be due to indirect effects, such as the build-up of toxic misfolded proteins, or to the direct effect of genes, such as aminoglycoside acetyl transferases, that are regulated by RpoH. Overall, the results demonstrate that RpoH regulates genes that are essential for viability of P. aeruginosa, that it protects P. aeruginosa from damage from aminoglycoside antibiotics, and that it is required for survival during nutrient-limiting conditions.

SurA-like and Skp-like Proteins as Important Virulence Determinants of the Gram Negative Bacterial Pathogens

In the Gram-negative bacteria, many important virulence factors reach their destination via two-step export systems, and they must traverse the periplasmic space before reaching the outer membrane. Since these proteins must be maintained in a structure competent for transport into or across the membrane, they frequently require the assistance of chaperones. Based on the results obtained for the model bacterium Escherichia coli and related species, it is assumed that in the biogenesis of the outer membrane proteins and the periplasmic transit of secretory proteins, the SurA peptidyl-prolyl isomerase/chaperone plays a leading role, while the Skp chaperone is rather of secondary importance. However, detailed studies carried out on several other Gram-negative pathogens indicate that the importance of individual chaperones in the folding and transport processes depends on the properties of client proteins and is species-specific. Taking into account the importance of SurA functions in bacterial virulence and severity of phenotypes due to surA mutations, this folding factor is considered as a putative therapeutic target to combat microbial infections. In this review, we present recent findings regarding SurA and Skp proteins: their mechanisms of action, involvement in processes related to virulence, and perspectives to use them as therapeutic targets.

Elucidation of Essential Genes and Mutant Fitness during Adaptation toward Nitrogen Fixation Conditions in the Endophyte Azoarcus olearius BH72 Revealed by Tn-Seq

Azoarcus olearius BH72 is a diazotrophic model endophyte that contributes fixed nitrogen to its host plant, Kallar grass, and expresses nitrogenase genes endophytically. Despite extensive studies on biological nitrogen fixation (BNF) of diazotrophic endophytes, little is known about global genetic players involved in survival under respective physiological conditions. Here, we report a global genomic screen for putatively essential genes of A. olearius employing Tn5 transposon mutagenesis with a modified transposon combined with high-throughput sequencing (Tn-Seq). A large Tn5 master library of ~6 × 105 insertion mutants of strain BH72 was obtained. Next-generation sequencing identified 183,437 unique insertion sites into the 4,376,040-bp genome, displaying one insertion every 24 bp on average. Applying stringent criteria, we describe 616 genes as putatively essential for growth on rich medium. COG (Clusters of Orthologous Groups) assignment of the 564 identified protein-coding genes revealed enrichment of genes related to core cellular functions and cell viability. To mimic gradual adaptations toward BNF conditions, the Tn5 mutant library was grown aerobically in synthetic medium or microaerobically on either combined or atmospheric nitrogen. Enrichment and depletion analysis of Tn5 mutants not only demonstrated the role of BNF- and metabolism-related proteins but also revealed that, strikingly, many genes relevant for plant-microbe interactions decrease bacterial competitiveness in pure culture, such type IV pilus- and bacterial envelope-associated genes. IMPORTANCE A constantly growing world population and the daunting challenge of climate change demand new strategies in agricultural crop production. Intensive usage of chemical fertilizers, overloading the world's fields with organic input, threaten terrestrial and marine ecosystems as well as human health. Long overlooked, the beneficial interaction of endophytic bacteria and grasses has attracted ever-growing interest in research in the last decade. Capable of biological nitrogen fixation, diazotrophic endophytes not only provide a valuable source of combined nitrogen but also are known for diverse plant growth-promoting effects, thereby contributing to plant productivity. Elucidation of an essential gene set for a prominent model endophyte such as A. olearius BH72 provides us with powerful insights into its basic lifestyle. Knowledge about genes detrimental or advantageous under defined physiological conditions may point out a way of manipulating key steps in the bacterium's lifestyle and plant interaction toward a more sustainable agriculture.

Genome-Wide Analysis of Gene Expression Noise Brought About by Transcriptional Regulation in Pseudomonas aeruginosa

The part of expression noise that is brought about by transcriptional regulation (represented here as NTR) is an important criterion for estimating the regulatory mode of a gene. However, characterization of NTR is an under-explored area, and there is little knowledge regarding the genome-wide NTR in the model pathogen Pseudomonas aeruginosa. Here, with a library of dual-color transcriptional reporters, we estimated the NTR for over 90% of the promoters in P. aeruginosa. Most promoters exhibit low NTR, while 42 and 115 promoters with high NTR were screened out in the exponential and the stationary growth phases, respectively. Specifically, a rearrangement of NTR was found in promoters involved in amino acid metabolism when bacteria enter the exponential phase. In addition, during the stationary phase, high NTR was found in a wide range of iron-related promoters involving siderophore synthesis and heme uptake, ExsA-regulated promoters involving bacterial virulence, and FleQ-regulated promoters involving biofilm development. We also found a large-scale negative dependence of transcriptional regulation between high-NTR promoters belonging to different functional categories. Our findings offer a global view of transcriptional heterogeneity in P. aeruginosa. IMPORTANCE The phenotypic diversity of Pseudomonas aeruginosa is frequently observed in research, suggesting that bacteria adopt strategies such as bet-hedging to survive ever-changing environments. Gene expression noise (GEN) is the major source of phenotypic diversity. Large GEN from transcriptional regulation (represented as NTR) represent an evolutionary necessity to maintain the copy number diversity of certain proteins in the population. Here, we provide a system-wide view of NTR in P. aeruginosa under nutrient-rich and stressed conditions. High NTR was found in genes involved in flagella biosynthesis and amino acid metabolism under both conditions. Specially, iron acquisition genes exhibited high NTR in the stressed condition, suggesting a great diversity of iron physiology in P. aeruginosa. We further revealed a global negative dependence of transcriptional regulation between those high-NTR genes under the stressed condition, suggesting a mutually exclusive relationship between different bacterial survival strategies.

Systems-Wide Dissection of Organic Acid Assimilation in Pseudomonas aeruginosa Reveals a Novel Path To Underground Metabolism

The human pathogen Pseudomonas aeruginosa (Pa) is one of the most frequent and severe causes of nosocomial infection. This organism is also a major cause of airway infections in people with cystic fibrosis (CF). Pa is known to have a remarkable metabolic plasticity, allowing it to thrive under diverse environmental conditions and ecological niches; yet, little is known about the central metabolic pathways that sustain its growth during infection or precisely how these pathways operate. In this work, we used a combination of 'omics approaches (transcriptomics, proteomics, metabolomics, and ¹³C-fluxomics) and reverse genetics to provide systems-level insight into how the infection-relevant organic acids succinate and propionate are metabolized by Pa. Moreover, through structural and kinetic analysis of the 2-methylcitrate synthase (2-MCS; PrpC) and its paralogue citrate (CIT) synthase (GltA), we show how these two crucial enzymatic steps are interconnected in Pa organic acid assimilation. We found that Pa can rapidly adapt to the loss of GltA function by acquiring mutations in a transcriptional repressor, which then derepresses prpC expression. Our findings provide a clear example of how "underground metabolism," facilitated by enzyme substrate promiscuity, "rewires" Pa metabolism, allowing it to overcome the loss of a crucial enzyme. This pathogen-specific knowledge is critical for the advancement of a model-driven framework to target bacterial central metabolism. IMPORTANCE Pseudomonas aeruginosa is an opportunistic human pathogen that, due to its unrivalled resistance to antibiotics, ubiquity in the built environment, and aggressiveness in infection scenarios, has acquired the somewhat dubious accolade of being designated a "critical priority pathogen" by the WHO. In this work, we uncover the pathways and mechanisms used by P. aeruginosa to grow on a substrate that is abundant at many infection sites: propionate. We found that if the organism is prevented from metabolizing propionate, the substrate turns from being a convenient nutrient source into a potent poison, preventing bacterial growth. We further show that one of the enzymes involved in these reactions, 2-methylcitrate synthase (PrpC), is promiscuous and can moonlight for another essential enzyme in the cell (citrate synthase). Indeed, mutations that abolish citrate synthase activity (which would normally prevent the cell from growing) can be readily overcome if the cell acquires additional mutations that increase the expression of PrpC. This is a nice example of the evolutionary utility of so-called "underground metabolism."

Gluconacetobacter diazotrophicus Gene Fitness during Diazotrophic Growth

Plant growth-promoting (PGP) bacteria are important to the development of sustainable agricultural systems. PGP microbes that fix atmospheric nitrogen (diazotrophs) could minimize the application of industrially derived fertilizers and function as a biofertilizer. The bacterium Gluconacetobacter diazotrophicus is a nitrogen-fixing PGP microbe originally discovered in association with sugarcane plants, where it functions as an endophyte. It also forms endophyte associations with a range of other agriculturally relevant crop plants. G. diazotrophicus requires microaerobic conditions for diazotrophic growth. We generated a transposon library for G. diazotrophicus and cultured the library under various growth conditions and culture medium compositions to measure fitness defects associated with individual transposon inserts (transposon insertion sequencing [Tn-seq]). Using this library, we probed more than 3,200 genes and ascertained the importance of various genes for diazotrophic growth of this microaerobic endophyte. We also identified a set of essential genes. IMPORTANCE Our results demonstrate a succinct set of genes involved in diazotrophic growth for G. diazotrophicus, with a lower degree of redundancy than what is found in other model diazotrophs. The results will serve as a valuable resource for those interested in biological nitrogen fixation and will establish a baseline data set for plant free growth, which could complement future studies related to the endophyte relationship.

Surviving the host: Microbial metabolic genes required for growth of Pseudomonas aeruginosa in physiologically-relevant conditions

Pseudomonas aeruginosa, like other pathogens, adapts to the limiting nutritional environment of the host by altering patterns of gene expression and utilizing alternative pathways required for survival. Understanding the genes essential for survival in the host gives insight into pathways that this organism requires during infection and has the potential to identify better ways to treat infections. Here, we used a saturated transposon insertion mutant pool of P. aeruginosa strain PAO1 and transposon insertion sequencing (Tn-Seq), to identify genes conditionally important for survival under conditions mimicking the environment of a nosocomial infection. Conditions tested included tissue culture medium with and without human serum, a murine abscess model, and a human skin organoid model. Genes known to be upregulated during infections, as well as those involved in nucleotide metabolism, and cobalamin (vitamin B₁₂) biosynthesis, etc., were required for survival in vivo- and in host mimicking conditions, but not in nutrient rich lab medium, Mueller Hinton broth (MHB). Correspondingly, mutants in genes encoding proteins of nucleotide and cobalamin metabolism pathways were shown to have growth defects under physiologically-relevant media conditions, in vivo, and in vivo-like models, and were downregulated in expression under these conditions, when compared to MHB. This study provides evidence for the relevance of studying P. aeruginosa fitness in physiologically-relevant host mimicking conditions and identified metabolic pathways that represent potential novel targets for alternative therapies.

Recombinant Actifensin and Defensin-d2 Induce Critical Changes in the Proteomes of Multidrug-Resistant Pseudomonas aeruginosa and Candida albicans

Drug-resistant strains of Pseudomonas aeruginosa and Candida albicans pose serious threats to human health because of their propensity to cause fatal infections. Defensin and defensin-like antimicrobial peptides (AMPs) are being explored as new lines of antimicrobials, due to their broad range of activity, low toxicity, and low pathogen resistance. Defensin-d2 and actifensin are AMPs from spinach and Actinomyces ruminicola, respectively, whose mechanisms of action are yet to be clearly elucidated. This study investigated the mechanisms of action of the recombinant AMPs through label-free quantitative proteomics. The data are available at PRIDE with accession number PXD034169. A total of 28 and 9 differentially expressed proteins (DEPs) were identified in the treated P. aeruginosa and C. albicans, respectively, with a 2-fold change threshold and P values of <0.05. Functional analysis revealed that the DEPs were involved in DNA replication and repair, translation, and membrane transport in P. aeruginosa, while they were related mainly to oxidative phosphorylation, RNA degradation, and energy metabolism in C. albicans. Protein-protein interactions showed that the DEPs formed linear or interdependent complexes with one another, indicative of functional interaction. Subcellular localization indicated that the majority of DEPs were cytoplasmic proteins in P. aeruginosa, while they were of nuclear or mitochondrial origin in C. albicans. These results show that recombinant defensin-d2 and actifensin can elicit complex multiple organism responses that cause cell death in P. aeruginosa and C. albicans. IMPORTANCE AMPs are considered essential alternatives to conventional antimicrobials because of their broad-spectrum efficacy and low potential for resistance by target cells. In this study, we established that the recombinant AMPs defensin-d2 and actifensin exert proteomic changes in P. aeruginosa and C. albicans within 1 h after treatment. We also found that the DEPs in peptide-treated P. aeruginosa are related to ion transport and homeostasis, molecular functions including nucleic and amino acid metabolism, and structural biogenesis and activity, while the DEPs in treated C. albicans are mainly involved in membrane synthesis and mitochondrial metabolism. Our results also highlight ATP synthase as a potential drug target for multidrug-resistant P. aeruginosa and C. albicans.

Pseudomonas aeruginosa reference strains PAO1 and PA14: A genomic, phenotypic, and therapeutic review

Pseudomonas aeruginosa is a ubiquitous, motile, gram-negative bacterium that has been recently identified as a multi-drug resistant pathogen in critical need of novel therapeutics. Of the approximately 5,000 strains, PAO1 and PA14 are common laboratory reference strains, modeling moderately and hyper-virulent phenotypes, respectively. PAO1 and PA14 have been instrumental in facilitating the discovery of novel drug targets, testing novel therapeutics, and supplying critical genomic information on the bacterium. While the two strains have contributed to a wide breadth of knowledge on the natural behaviors and therapeutic susceptibilities of P. aeruginosa, they have demonstrated significant deviations from observations in human infections. Many of these deviations are related to experimental inconsistencies in laboratory strain environment that complicate and, at times, terminate translation from laboratory results to clinical applications. This review aims to provide a comparative analysis of the two strains and potential methods to improve their clinical relevance.

BosR: A novel biofilm-specific regulator in Pseudomonas aeruginosa

Biofilms are the most common cause of bacterial infections in humans and notoriously hard to treat due to their ability to withstand antibiotics and host immune defenses. To overcome the current lack of effective antibiofilm therapies and guide future design, the identification of novel biofilm-specific gene targets is crucial. In this regard, transcriptional regulators have been proposed as promising targets for antimicrobial drug design. Therefore, a Transposon insertion sequencing approach was employed to systematically identify regulators phenotypically affecting biofilm growth in Pseudomonas aeruginosa PA14 using the TnSeq analysis tools Bio-TraDIS and TRANSIT. A screen of a pool of 300,000 transposon insertion mutants identified 349 genes involved in biofilm growth on hydroxyapatite, including 47 regulators. Detection of 19 regulatory genes participating in well-established biofilm pathways validated the results. An additional 28 novel prospective biofilm regulators suggested the requirement for multiple one-component transcriptional regulators. Biofilm-defective phenotypes were confirmed for five one-component transcriptional regulators and a protein kinase, which did not affect motility phenotypes. The one-component transcriptional regulator bosR displayed a conserved role in P. aeruginosa biofilm growth since its ortholog in P. aeruginosa strain PAO1 was also required for biofilm growth. Microscopic analysis of a chromosomal deletion mutant of bosR confirmed the role of this regulator in biofilm growth. Overall, our results highlighted that the gene network driving biofilm growth is complex and involves regulators beyond the primarily studied groups of two-component systems and cyclic diguanylate signaling proteins. Furthermore, biofilm-specific regulators, such as bosR, might constitute prospective new drug targets to overcome biofilm infections.

A bacterial pan-genome makes gene essentiality strain-dependent and evolvable

Many bacterial species are represented by a pan-genome, whose genetic repertoire far outstrips that of any single bacterial genome. Here we investigate how a bacterial pan-genome might influence gene essentiality and whether essential genes that are initially critical for the survival of an organism can evolve to become non-essential. By using Transposon insertion sequencing (Tn-seq), whole-genome sequencing and RNA-seq on a set of 36 clinical Streptococcus pneumoniae strains representative of >68% of the species' pan-genome, we identify a species-wide 'essentialome' that can be subdivided into universal, core strain-specific and accessory essential genes. By employing 'forced-evolution experiments', we show that specific genetic changes allow bacteria to bypass essentiality. Moreover, by untangling several genetic mechanisms, we show that gene essentiality can be highly influenced by and/or be dependent on: (1) the composition of the accessory genome, (2) the accumulation of toxic intermediates, (3) functional redundancy, (4) efficient recycling of critical metabolites and (5) pathway rewiring. While this functional characterization underscores the evolvability potential of many essential genes, we also show that genes with differential essentiality remain important antimicrobial drug target candidates, as their inactivation almost always has a severe fitness cost in vivo.

Experimental and Analytical Approaches for Improving the Resolution of Randomly Barcoded Transposon Insertion Sequencing (RB-TnSeq) Studies

Randomly barcoded transposon insertion sequencing (RB-TnSeq) is an efficient, multiplexed method to determine microbial gene function during growth under a selection condition of interest. This technique applies to growth, tolerance, and persistence studies in a variety of hosts, but the wealth of data generated can complicate the identification of the most critical gene targets. Experimental and analytical methods for improving the resolution of RB-TnSeq are proposed, using Pseudomonas putida KT2440 as an example organism. Several key parameters, such as baseline media selection, substantially influence the determination of gene fitness. We also present options to increase statistical confidence in gene fitness, including increasing the number of biological replicates and passaging the baseline culture in parallel with selection conditions. These considerations provide practitioners with several options to identify genes of importance in TnSeq data sets, thereby streamlining metabolic characterization.

The Impact of Chromate on Pseudomonas aeruginosa Molybdenum Homeostasis

Acquisition of the trace-element molybdenum via the high-affinity ATP-binding cassette permease ModABC is essential for Pseudomonas aeruginosa respiration in anaerobic and microaerophilic environments. This study determined the X-ray crystal structures of the molybdenum-recruiting solute-binding protein ModA from P. aeruginosa PAO1 in the metal-free state and bound to the group 6 metal oxyanions molybdate, tungstate, and chromate. Pseudomonas aeruginosa PAO1 ModA has a non-contiguous dual-hinged bilobal structure with a single metal-binding site positioned between the two domains. Metal binding results in a 22° relative rotation of the two lobes with the oxyanions coordinated by four residues, that contribute six hydrogen bonds, distinct from ModA orthologues that feature an additional oxyanion-binding residue. Analysis of 485 Pseudomonas ModA sequences revealed conservation of the metal-binding residues and β-sheet structural elements, highlighting their contribution to protein structure and function. Despite the capacity of ModA to bind chromate, deletion of modA did not affect P. aeruginosa PAO1 sensitivity to chromate toxicity nor impact cellular accumulation of chromate. Exposure to sub-inhibitory concentrations of chromate broadly perturbed P. aeruginosa metal homeostasis and, unexpectedly, was associated with an increase in ModA-mediated molybdenum uptake. Elemental analyses of the proteome from anaerobically grown P. aeruginosa revealed that, despite the increase in cellular molybdenum upon chromate exposure, distribution of the metal within the proteome was substantially perturbed. This suggested that molybdoprotein cofactor acquisition may be disrupted, consistent with the potent toxicity of chromate under anaerobic conditions. Collectively, these data reveal a complex relationship between chromate toxicity, molybdenum homeostasis and anaerobic respiration.

Integrated genomics and chemical biology herald an era of sophisticated antibacterial discovery, from defining essential genes to target elucidation

The golden age of antibiotic discovery in the 1940s-1960s saw the development and deployment of many different classes of antibiotics, revolutionizing the field of medicine. Since that time, our ability to discover antibiotics of novel structural classes or mechanisms has not kept pace with the ever-growing threat of antibiotic resistance. Recently, advances at the intersection of genomics and chemical biology have enabled efforts to better define the vulnerabilities of essential gene targets, to develop sophisticated whole-cell chemical screening methods that reveal target biology early, and to elucidate small molecule targets and modes of action more effectively. These new technologies have the potential to expand the chemical diversity of antibiotic candidates, as well as the breadth of targets. We illustrate how the latest tools of genomics and chemical biology are being integrated to better understand pathogen vulnerabilities and antibiotic mechanisms in order to inform a new era of antibiotic discovery.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

A curated collection of Klebsiella metabolic models reveals variable substrate usage and gene essentiality

The Klebsiella pneumoniae species complex (KpSC) is a set of seven Klebsiella taxa that are found in a variety of niches and are an important cause of opportunistic health care-associated infections in humans. Because of increasing rates of multi-drug resistance within the KpSC, there is a growing interest in better understanding the biology and metabolism of these organisms to inform novel control strategies. We collated 37 sequenced KpSC isolates isolated from a variety of niches, representing all seven taxa. We generated strain-specific genome-scale metabolic models (GEMs) for all 37 isolates and simulated growth phenotypes on 511 distinct carbon, nitrogen, sulfur, and phosphorus substrates. Models were curated and their accuracy was assessed using matched phenotypic growth data for 94 substrates (median accuracy of 96%). We explored species-specific growth capabilities and examined the impact of all possible single gene deletions using growth simulations in 145 core carbon substrates. These analyses revealed multiple strain-specific differences, within and between species, and highlight the importance of selecting a diverse range of strains when exploring KpSC metabolism. This diverse set of highly accurate GEMs could be used to inform novel drug design, enhance genomic analyses, and identify novel virulence and resistance determinants. We envisage that these 37 curated strain-specific GEMs, covering all seven taxa of the KpSC, provide a valuable resource to the Klebsiella research community.