Defining the core essential genome of Pseudomonas aeruginosa

Towards new antibiotics: P. aeruginosa FabF ligands discovered by crystallographic fragment screening followed by hit expansion

There is an urgent need for new antibiotics. FabF (3-oxoacyl-[acyl-carrier-protein] synthase 2), which catalyses the rate limiting condensation reaction in the fatty acid synthesis II pathway, is an attractive target. Very few inhibitors of FabF are known and most are derived from natural products. In an effort to further explore the chemical space of FabF ligands, we have carried out fragment screening by X-ray crystallography against an intermediated state-mimicking variant of P. aeruginosa FabF (PaFabF C164Q). This screen has resulted in 48 hits out of which 16 bind in or close to the malonyl-CoA or fatty acid binding site or an adjacent dimer interface. None of the closer investigated fragments were active in a binding assay, but the same was the case for fragments derived from a potent FabF inhibitor. For hit optimization, we focused on the two fragments binding close to the catalytic residues of FabF. Different strategies were followed in the optimization process: exploration of commercially available analogues, fragment merging, virtual screening of a combinatorial make-on-demand space, and design and in-house synthesis of analogues. In total, more than 90 analogues of the hit compounds were explored, and for 10 of those co-crystal structures could be determined. The most potent ligand was discovered using manual structure-based design and has a binding affinity of 65 μM. This data package forms a strong foundation for the development of more potent and diverse FabF inhibitors.

Quinoline ATP Synthase Inhibitors with Activity Against Multidrug Resistant Acinetobacter baumannii and Pseudomonas aeruginosa

The Gram-negative, pathogenic bacteria Acinetobacter baumannii (AB) and Pseudomonas aeruginosa (PA) have been identified as a particular threat due to rising multidrug resistance, and antibiotics with novel mechanisms of action are needed. Bacterial bioenergetics is a promising but underdeveloped drug target since the complexes of oxidative phosphorylation are critical to cell survival in these organisms. Building from our previous work using quinoline derivatives to inhibit the ATP synthase of PA, we report a new set of 14 quinoline derivatives that demonstrates potent inhibition of the AB ATP synthase, with the best inhibitor having an IC50 of 230 ng/mL in vitro, expands the quinoline structure-activity relationship against the PA enzyme, and establishes molecular strategies for achieving selectivity between PA and AB. Furthermore, several compounds demonstrated potent antibacterial activity against multidrug resistant strains of AB and PA indicating ATP synthase as a promising new area for broad spectrum antibiotic development in AB.

The Iron Metalloproteome of Pseudomonas aeruginosa Under Oxic and Anoxic Conditions

Pseudomonas aeruginosa is a major contributor to human infections and is widely distributed in the environment. Its ability for growth under aerobic and anaerobic conditions provides adaptability to environmental changes and in confronting immune responses. We applied native 2-dimensional metalloproteomics to P. aeruginosa to examine how use of iron within the metallome responds to oxic and anoxic conditions. Analyses revealed four iron peaks comprised of metalloproteins with synergistic functions, including: 1) respiratory and metabolic enzymes, 2) oxidative stress response enzymes, 3) DNA synthesis and nitrogen assimilation enzymes, and 4) denitrification enzymes and related copper enzymes. Fe peaks were larger under anoxic conditions, consistent with increased iron demand due to anaerobic metabolism and with the denitrification peak absent under oxic conditions. Three ferritins co-eluted with the first and third iron peaks, localizing iron storage with these functions. Several enzymes were more abundant at low oxygen, including alkylhydroperoxide reductase C that deactivates organic radicals produced by denitrification, all three classes of ribonucleotide reductases (including monomers and oligomer forms), ferritin (increasing in ratio relative to bacterioferritin), and denitrification enzymes. Superoxide dismutase and homogentisate 1,2-dioxygenase were more abundant at high oxygen. Several Fe peaks contained iron metalloproteins that co-eluted earlier than their predicted size, implying additional protein-protein interactions and suggestive of cellular organization that contributes to iron prioritization in Pseudomonas with its large genome and flexible metabolism. This study characterized the iron metalloproteome of one of the more complex prokaryotic microorganisms, attributing enhanced iron use under anaerobic denitrifying metabolism to its specific metalloprotein constituents.

Throwing a spotlight on genomic dark matter: the power and potential of transposon-insertion sequencing

Linking genotype to phenotype is a central goal in biology. In the microbiological field, transposon mutagenesis is a technique that has been widely used since the 1970's to facilitate this connection. The development of modern 'omics approaches and next-generation sequencing, have allowed high-throughput association between genes and their putative function. In 2009, four different variations of modern transposon-insertion sequencing (TIS) approaches were published, being referred to as transposon-directed insertion-site sequencing (TraDIS), transposon sequencing (Tn-seq), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). These approaches exploit a similar concept to allow estimation of the essentiality or contribution to fitness of each gene in any bacterial genome. The main rationale is to perform a comparative analysis of the abundance of specific transposon mutants under one or more selective conditions. The approaches themselves only vary in the transposon used for mutagenesis, and in the methodology used for sequencing library preparation. In this review, we discuss how TIS approaches have been used to facilitate a major shift in our fundamental understanding of bacterial biology in a range of areas. We focus on several aspects including pathogenesis, biofilm development, polymicrobial interactions in various ecosystems, and antimicrobial resistance. These studies have provided new insight into bacterial physiology and revealed predicted functions for hundreds of genes previously representing genomic 'dark matter'. We also discuss how TIS approaches have been used to understand complex bacterial systems and interactions and how future developments of TIS could continue to accelerate and enrich our understanding of bacterial biology.

Architecture of Pseudomonas aeruginosa glutamyl-tRNA synthetase defines a subfamily of dimeric class Ib aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (AaRSs) are an ancient family of structurally diverse enzymes that are divided into two major classes. The functionalities of most AaRSs are inextricably linked to their oligomeric states. While GluRSs were previously classified as monomers, the current investigation reveals that the form expressed in Pseudomonas aeruginosa is a rotationally pseudosymmetrical homodimer featuring intersubunit tRNA binding sites. Both subunits display a highly bent, "pipe strap" conformation, with the anticodon binding domain directed toward the active site. The tRNA binding sites are similar in shape to those of the monomeric GluRSs, but are formed through an approximately 180-degree rotation of the anticodon binding domains and dimerization via the anticodon and D-arm binding domains. As a result, each anticodon binding domain is poised to recognize the anticodon loop of a tRNA bound to the adjacent protomer. Additionally, the anticodon binding domain has an α-helical C-terminal extension containing a conserved lysine-rich consensus motif positioned near the predicted location of the acceptor arm, suggesting dual functions in tRNA recognition. The unique architecture of PaGluRS broadens the structural diversity of the GluRS family, and member synthetases of all bacterial AaRS subclasses have now been identified that exhibit oligomerization.

Carbapenem-Resistant Pseudomonas aeruginosa's Resistome: Pan-Genomic Plasticity, the Impact of Transposable Elements and Jumping Genes

Pseudomonas aeruginosa, a Gram-negative, motile bacterium, may cause significant infections in both community and hospital settings, leading to substantial morbidity and mortality. This opportunistic pathogen can thrive in various environments, making it a public health concern worldwide. P. aeruginosa's genomic pool is highly dynamic and diverse, with a pan-genome size ranging from 5.5 to 7.76 Mbp. This versatility arises from its ability to acquire genes through horizontal gene transfer (HGT) via different genetic elements (GEs), such as mobile genetic elements (MGEs). These MGEs, collectively known as the mobilome, facilitate the spread of genes encoding resistance to antimicrobials (ARGs), resistance to heavy metals (HMRGs), virulence (VGs), and metabolic functions (MGs). Of particular concern are the acquired carbapenemase genes (ACGs) and other β-lactamase genes, such as classes A, B [metallo-β-lactamases (MBLs)], and D carbapenemases, which can lead to increased antimicrobial resistance. This review emphasizes the importance of the mobilome in understanding antimicrobial resistance in P. aeruginosa.

Discovery of a Pseudomonas aeruginosa-specific small molecule targeting outer membrane protein OprH-LPS interaction by a multiplexed screen

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, 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 pathogen biology.

Anaerobiosis and Mutations Can Reduce Susceptibility of Pseudomonas aeruginosa to Tobramycin Without Reducing the Cellular Concentration of the Antibiotic

Chronic infections of Pseudomonas aeruginosa are commonly treated with tobramycin. During infections, the bacteria can exist under conditions of oxygen deprivation that render them less susceptible to this antibiotic. The aims of this research were to investigate the genetic basis of tobramycin resistance under anaerobic conditions, and to investigate the effects of anaerobiosis and mutations on the cellular concentration of tobramycin. Ten mutants with lowered susceptibility to tobramycin than wild-type bacteria were evolved from a laboratory reference strain under anaerobic conditions. Mutations were identified by genome sequencing. Mutations had arisen most frequently in the fusA1 gene that encodes elongation factor EF-G1A and in genes involved in twitching motility. Cellular concentrations of tobramycin were then measured. Mutations in fusA1 or absence of the MexXY efflux pump that is associated with tobramycin resistance did not alter the cellular tobramycin concentration under either anaerobic or aerobic conditions. Anaerobic growth reduced the cellular concentration of tobramycin, relative to aerobically grown bacteria, in some but not all of five tested P. aeruginosa isolates. Overall, our findings indicate that anaerobiosis and mutations that reduce aminoglycoside effectiveness do not lower the cellular concentration of antibiotic but instead reduce susceptibility through other mechanisms.

Micromix: web infrastructure for visualizing and remixing microbial 'omics data

Micromix is a flexible web platform for sharing and integrating microbial omics data, including RNA sequencing and transposon-insertion sequencing. Currently, the lack of solutions for making data web-accessible results in omics data being fragmented across supplementary spreadsheets or languishing as raw read data in public repositories. Micromix solves this problem and can be easily deployed on a standard web server or using cloud services. It is organism-agnostic, accommodates data and annotations from various sources, and allows filtering based on KEGG pathways, Gene Ontology terms, and curated gene sets. Visualizations are provided through a plug-in system that integrates existing visualization services and allows rapid development of new services, with available plug-ins currently supporting interactive heatmap and clustering functions. Users can upload their own data in a variety of formats to perform integrative analyses in the context of existing datasets. To support collaborative research, Micromix allows sharing of interactive sessions that maintain defined filtering and/or visualization options. We demonstrate the utility of Micromix with case studies focusing on the SPI-2 pathogenicity island in Salmonella enterica and polysaccharide utilization loci in Bacteroides thetaiotaomicron, showcasing the platform's capabilities for integrating, filtering, and visualizing diverse functional genomic datasets. Micromix is available at http://micromix.systems.

Pangenome graphs and their applications in biodiversity genomics

Complete datasets of genetic variants are key to biodiversity genomic studies. Long-read sequencing technologies allow the routine assembly of highly contiguous, haplotype-resolved reference genomes. However, even when complete, reference genomes from a single individual may bias downstream analyses and fail to adequately represent genetic diversity within a population or species. Pangenome graphs assembled from aligned collections of high-quality genomes can overcome representation bias by integrating sequence information from multiple genomes from the same population, species or genus into a single reference. Here, we review the available tools and data structures to build, visualize and manipulate pangenome graphs while providing practical examples and discussing their applications in biodiversity and conservation genomics across the tree of life.

Urinary tract infections and catheter-associated urinary tract infections caused by Pseudomonas aeruginosa

SUMMARYUrinary tract infection (UTI) is one of the most common infections in otherwise healthy individuals. UTI is also common in healthcare settings where patients often require urinary catheters to alleviate urinary retention. The placement of a urinary catheter often leads to catheter-associated urinary tract infection (CAUTI) caused by a broad range of opportunistic pathogens, commonly referred to as ESKAPE (Enterococcus, Staphylococcus, Klebsiella, Acinetobacter, Pseudomonas, and Enterobacter) pathogens. Our understanding of CAUTI is complicated by the differences in pathogens, in initial microbial load, changes that occur due to the duration of catheterization, and the relationship between infection (colonization) and disease symptoms. To advance our understanding of CAUTI, we reviewed UTI and CAUTI caused by Pseudomonas aeruginosa which is unique in that it is not commonly found associated with human microbiomes. For this reason, the ability of P. aeruginosa to cause UTI and CAUTI requires the introduction of the bacteria to the bladder from catheterization. Once in the host, the virulence factors used by P. aeruginosa in these infections remain an area of ongoing research. In this review, we will discuss studies that focus on P. aeruginosa UTI and CAUTI to better understand the infection dynamics and outcome in clinical settings, virulence factors associated with P. aeruginosa isolated from the urinary tract, and animal studies to test which bacterial factors are required for this infection. Understanding how P. aeruginosa can cause UTI and CAUTI can provide an understanding of how these infections initiate and progress and may provide possible strategies to limit these infections.

The identification of essential cellular genes is critical for validating drug targets

Accurately identifying biological targets is crucial for advancing treatment options. Essential genes, vital for cell or organism survival, hold promise as potential drug targets in disease treatment. Although many studies have sought to identify essential genes as therapeutic targets in medicine and bioinformatics, systematic reviews on their relationship with drug targets are relatively rare. This work presents a comprehensive analysis to aid in identifying essential genes as potential targets for drug discovery, encompassing their relevance, identification methods, successful case studies, and challenges. This work will facilitate the identification of essential genes as therapeutic targets, thereby boosting new drug development.

A molecular comparison of [Fe-S] cluster-based homeostasis in Escherichia coli and Pseudomonas aeruginosa

Iron-sulfur [Fe-S] clusters are essential protein cofactors allowing bacteria to perceive environmental redox modification and to adapt to iron limitation. Escherichia coli, which served as a bacterial model, contains two [Fe-S] cluster biogenesis systems, ISC and SUF, which ensure [Fe-S] cluster synthesis under balanced and stress conditions, respectively. However, our recent phylogenomic analyses revealed that most bacteria possess only one [Fe-S] cluster biogenesis system, most often SUF. The opportunist human pathogen Pseudomonas aeruginosa is atypical as it harbors only ISC. Here, we confirmed the essentiality of ISC in P. aeruginosa under both normal and stress conditions. Moreover, P. aeruginosa ISC restored viability, under balanced growth conditions, to an E. coli strain lacking both ISC and SUF. Reciprocally, the E. coli SUF system sustained growth and [Fe-S] cluster-dependent enzyme activities of ISC-deficient P. aeruginosa. Surprisingly, an ISC-deficient P. aeruginosa strain expressing E. coli SUF showed defects in resistance to H2O2 stress and paraquat, a superoxide generator. Similarly, the P. aeruginosa ISC system did not confer stress resistance to a SUF-deficient E. coli mutant. A survey of 120 Pseudomonadales genomes confirmed that all but five species have selected ISC over SUF. While highlighting the great versatility of bacterial [Fe-S] cluster biogenesis systems, this study emphasizes that their contribution to cellular homeostasis must be assessed in the context of each species and its own repertoire of stress adaptation functions. As a matter of fact, despite having only one ISC system, P. aeruginosa shows higher fitness in the face of ROS and iron limitation than E. coli.

IMPORTANCE

ISC and SUF molecular systems build and transfer Fe-S cluster to cellular apo protein clients. The model Escherichia coli has both ISC and SUF and study of the interplay between the two systems established that the ISC system is the house-keeping one and SUF the stress-responding one. Unexpectedly, our recent phylogenomic analysis revealed that in contrast to E. coli (and related enterobacteria such as Salmonella), most bacteria have only one system, and, in most cases, it is SUF. Pseudomonas aeruginosa fits the general rule of having only one system but stands against the rule by having ISC. This study aims at engineering P. aeruginosa harboring E. coli systems and vice versa. Comparison of the recombinants allowed to assess the functional versatility of each system while appreciating their contribution to cellular homeostasis in different species context.

High-Throughput Discovery of Synthetic Siderophores for Trojan Horse Antibiotics

To cause infection, bacterial pathogens must overcome host immune factors and barriers to nutrient acquisition. Reproducing these aspects of host physiology in vitro has shown great promise for antibacterial drug discovery. When used as a bacterial growth medium, human serum replicates several aspects of the host environment, including innate immunity and iron limitation. We previously reported that a high-throughput chemical screen using serum as the growth medium enabled the discovery of novel growth inhibitors overlooked by conventional screens. Here, we report that a subset of compounds from this high-throughput serum screen display an unexpected growth enhancing phenotype and are enriched for synthetic siderophores. We selected 35 compounds of diverse chemical structure and quantified their ability to enhance bacterial growth in human serum. We show that many of these compounds chelate iron, suggesting they were acting as siderophores and providing iron to the bacteria. For two different pharmacophores represented among these synthetic siderophores, conjugation to the β-lactam antibiotic ampicillin imparted iron-dependent enhancement in antibacterial activity. Conjugation of the most potent growth-enhancing synthetic siderophore with the monobactam aztreonam produced MLEB-22043, a broad-spectrum antibiotic with significantly improved activity against Klebsiella pneumoniae, Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa. This synthetic siderophore-monobactam conjugate uses multiple TonB-dependent transporters for uptake into P. aeruginosa. Like aztreonam, MLEB-22043 demonstrated activity against metallo-β-lactamase expressing bacteria, and, when combined with the β-lactamase inhibitor avibactam, was active against clinical strains coexpressing the NDM-1 metallo-β-lactamase and serine β-lactamases. Our work shows that human serum is an effective bacterial growth medium for the high-throughput discovery of synthetic siderophores, enabling the development of novel Trojan Horse antibiotics.

High-throughput transposon mutagenesis in the family Enterobacteriaceae reveals core essential genes and rapid turnover of essentiality

The Enterobacteriaceae are a scientifically and medically important clade of bacteria, containing the model organism Escherichia coli, as well as major human pathogens including Salmonella enterica and Klebsiella pneumoniae. Essential gene sets have been determined for several members of the Enterobacteriaceae, with the Keio E. coli single-gene deletion library often regarded as a gold standard. However, it remains unclear how gene essentiality varies between related strains and species. To investigate this, we have assembled a collection of 13 sequenced high-density transposon mutant libraries from five genera within the Enterobacteriaceae. We first assess several gene essentiality prediction approaches, investigate the effects of transposon density on essentiality prediction, and identify biases in transposon insertion sequencing data. Based on these investigations, we develop a new classifier for gene essentiality. Using this new classifier, we define a core essential genome in the Enterobacteriaceae of 201 universally essential genes. Despite the presence of a large cohort of variably essential genes, we find an absence of evidence for genus-specific essential genes. A clear example of this sporadic essentiality is given by the set of genes regulating the σE extracytoplasmic stress response, which appears to have independently acquired essentiality multiple times in the Enterobacteriaceae. Finally, we compare our essential gene sets to the natural experiment of gene loss in obligate insect endosymbionts that have emerged from within the Enterobacteriaceae. This isolates a remarkably small set of genes absolutely required for survival and identifies several instances of essential stress responses masked by redundancy in free-living bacteria.IMPORTANCEThe essential genome, that is the set of genes absolutely required to sustain life, is a core concept in genetics. Essential genes in bacteria serve as drug targets, put constraints on the engineering of biological chassis for technological or industrial purposes, and are key to constructing synthetic life. Despite decades of study, relatively little is known about how gene essentiality varies across related bacteria. In this study, we have collected gene essentiality data for 13 bacteria related to the model organism Escherichia coli, including several human pathogens, and investigated the conservation of essentiality. We find that approximately a third of the genes essential in any particular strain are non-essential in another related strain. Surprisingly, we do not find evidence for essential genes unique to specific genera; rather it appears a substantial fraction of the essential genome rapidly gains or loses essentiality during evolution. This suggests that essentiality is not an immutable characteristic but depends crucially on the genomic context. We illustrate this through a comparison of our essential genes in free-living bacteria to genes conserved in 34 insect endosymbionts with naturally reduced genomes, finding several cases where genes generally regarded as being important for specific stress responses appear to have become essential in endosymbionts due to a loss of functional redundancy in the genome.

Genome-wide fitness analysis of Salmonella enterica reveals aroA mutants are attenuated due to iron restriction in vitro

Salmonella enterica is a globally disseminated pathogen that is the cause of over 100 million infections per year. The resulting diseases are dependent upon host susceptibility and the infecting serovar. As S. enterica serovar Typhimurium induces a typhoid-like disease in mice, this model has been used extensively to illuminate various aspects of Salmonella infection and host responses. Due to the severity of infection in this model, researchers often use strains of mice resistant to infection or attenuated Salmonella. Despite decades of research, many aspects of Salmonella infection and fundamental biology remain poorly understood. Here, we use a transposon insertion sequencing technique to interrogate the essential genomes of widely used isogenic wild-type and attenuated S. Typhimurium strains. We reveal differential essential pathways between strains in vitro and provide a direct link between iron starvation, DNA synthesis, and bacterial membrane integrity.IMPORTANCESalmonella enterica is an important clinical pathogen that causes a high number of deaths and is increasingly resistant to antibiotics. Importantly, S. enterica is used widely as a model to understand host responses to infection. Understanding how Salmonella survives in vivo is important for the design of new vaccines to combat this pathogen. Live attenuated vaccines have been used clinically for decades. A widely used mutation, aroA, is thought to attenuate Salmonella by restricting the ability of the bacterium to access particular amino acids. Here we show that this mutation limits the ability of Salmonella to acquire iron. These observations have implications for the interpretation of many previous studies and for the use of aroA in vaccine development.

Bismuth-based drugs sensitize Pseudomonas aeruginosa to multiple antibiotics by disrupting iron homeostasis

Pseudomonas aeruginosa infections are difficult to treat due to rapid development of antibiotic drug resistance. The synergistic combination of already-in-use drugs is an alternative to developing new antibiotics to combat antibiotic-resistant bacteria. Here we demonstrate that bismuth-based drugs (bismuth subsalicylate, colloidal bismuth subcitrate) in combination with different classes of antibiotics (tetracyclines, macrolides, quinolones, rifamycins and so on) can eliminate multidrug-resistant P. aeruginosa and do not induce development of antibiotic resistance. Bismuth disrupts iron homeostasis by binding to P. aeruginosa siderophores. Inside cells, bismuth inhibits the electron transport chain, dissipates the proton motive force and impairs efflux pump activity by disrupting iron-sulfur cluster-containing enzymes, including respiration complexes. As a result, bismuth facilitates antibiotic accumulation inside bacteria, enhancing their efficacy. The combination therapy shows potent antibacterial efficacy and low toxicity in an ex vivo bacteraemia model and increases the survival rate of mice in in vivo mouse lung-infection models. Our findings highlight the potential of bismuth-based drugs to be repurposed to combat P. aeruginosa infections in combination with clinically used antibiotics.

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.

'Wild Type'

In this opinion piece, we consider the meaning of the term 'wild type' in the context of microbiology. This is especially pertinent in the post-genomic era, where we have a greater awareness of species diversity than ever before. Genomic heterogeneity, in vitro evolution/selection pressures, definition of 'the wild', the size and importance of the pan-genome, gene-gene interactions (epistasis), and the nature of the 'wild-type gene' are all discussed. We conclude that wild type is an outdated and even misleading phrase that should be gradually phased out.

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.

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.

Genome-wide fitness analysis identifies genes required for in vitro growth and macrophage infection by African and global epidemic pathovariants of Salmonella enterica Enteritidis

Salmonella enterica Enteritidis is the second most common serovar associated with invasive non-typhoidal Salmonella (iNTS) disease in sub-Saharan Africa. Previously, genomic and phylogenetic characterization of S . enterica Enteritidis isolates from the human bloodstream led to the discovery of the Central/Eastern African clade (CEAC) and West African clade, which were distinct from the gastroenteritis-associated global epidemic clade (GEC). The African S . enterica Enteritidis clades have unique genetic signatures that include genomic degradation, novel prophage repertoires and multi-drug resistance, but the molecular basis for the enhanced propensity of African S . enterica Enteritidis to cause bloodstream infection is poorly understood. We used transposon insertion sequencing (TIS) to identify the genetic determinants of the GEC representative strain P125109 and the CEAC representative strain D7795 for growth in three in vitro conditions (LB or minimal NonSPI2 and InSPI2 growth media), and for survival and replication in RAW 264.7 murine macrophages. We identified 207 in vitro -required genes that were common to both S . enterica Enteritidis strains and also required by S . enterica Typhimurium, S . enterica Typhi and Escherichia coli , and 63 genes that were only required by individual S . enterica Enteritidis strains. Similar types of genes were required by both P125109 and D7795 for optimal growth in particular media. Screening the transposon libraries during macrophage infection identified 177 P125109 and 201 D7795 genes that contribute to bacterial survival and replication in mammalian cells. The majority of these genes have proven roles in Salmonella virulence. Our analysis uncovered candidate strain-specific macrophage fitness genes that could encode novel Salmonella virulence factors.

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.