BACTOME-a reference database to explore the sequence- and gene expression-variation landscape of Pseudomonas aeruginosa clinical isolates

Extracellular vesicles released by host epithelial cells during Pseudomonas aeruginosa infection function as homing beacons for neutrophils

BACKGROUND

Pseudomonas aeruginosa (PA) is an opportunistic pathogen that can cause sight threatening infections in the eye and fatal infections in the cystic fibrosis airway. Extracellular vesicles (EVs) are released by host cells during infection and by the bacteria themselves; however, there are no studies on the composition and functional role of host-derived EVs during PA infection of the eye or lung. Here we investigated the composition and capacity of EVs released by PA infected epithelial cells to modulate innate immune responses in host cells.

METHODS

Human telomerase immortalized corneal epithelial cells (hTCEpi) cells and human telomerase immortalized bronchial epithelial cells (HBECs) were treated with a standard invasive test strain of Pseudomonas aeruginosa, PAO1, for 6 h. Host derived EVs were isolated by qEV size exclusion chromatography. EV proteomic profiles during infection were compared using mass spectrometry and functional studies were carried out using hTCEpi cells, HBECs, differentiated neutrophil-like HL-60 cells, and primary human neutrophils isolated from peripheral blood.

RESULTS

EVs released from PA infected corneal epithelial cells increased pro-inflammatory cytokine production in naïve corneal epithelial cells and induced neutrophil chemotaxis independent of cytokine production. The EVs released from PA infected bronchial epithelial cells were also chemotactic although they failed to induce cytokine secretion from naïve HBECs. At the proteomic level, EVs derived from PA infected corneal epithelial cells exhibited lower complexity compared to bronchial epithelial cells, with the latter having reduced protein expression compared to the non-infected control.

CONCLUSIONS

This is the first study to comprehensively profile EVs released by corneal and bronchial epithelial cells during Pseudomonas infection. Together, these findings show that EVs released by PA infected corneal and bronchial epithelial cells function as potent mediators of neutrophil migration, contributing to the exuberant neutrophil response that occurs during infection in these tissues.

tRNA epitranscriptome determines pathogenicity of the opportunistic pathogen Pseudomonas aeruginosa

The success of bacterial pathogens depends on the coordinated expression of virulence determinants. Regulatory circuits that drive pathogenesis are complex, multilayered, and incompletely understood. Here, we reveal that alterations in tRNA modifications define pathogenic phenotypes in the opportunistic pathogen Pseudomonas aeruginosa. We demonstrate that the enzymatic activity of GidA leads to the introduction of a carboxymethylaminomethyl modification in selected tRNAs. Modifications at the wobble uridine base (cmnm5U34) of the anticodon drives translation of transcripts containing rare codons. Specifically, in P. aeruginosa the presence of GidA-dependent tRNA modifications modulates expression of genes encoding virulence regulators, leading to a cellular proteomic shift toward pathogenic and well-adapted physiological states. Our approach of profiling the consequences of chemical tRNA modifications is general in concept. It provides a paradigm of how environmentally driven tRNA modifications govern gene expression programs and regulate phenotypic outcomes responsible for bacterial adaption to challenging habitats prevailing in the host niche.

A comprehensive review of genomics, transcriptomics, proteomics, and metabolomic insights into the differentiation of Pseudomonas aeruginosa from the planktonic to biofilm state: A multi-omics approach

Biofilm formation by Pseudomonas aeruginosa is primarily responsible for chronic wound and lung infections in humans. These infections are persistent owing to the biofilm's high tolerance to antimicrobials and constantly changing environmental factors. Understanding the mechanism governing biofilm formation can help to develop therapeutics explicitly directed against the molecular markers responsible for this process. After numerous years of research, many genes responsible for both in vitro and in vivo biofilm development remain unidentified. However, there is no "all in one" complete in vivo or in vitro biofilm model. Recent findings imply that the shift from planktonic bacteria to biofilms is a complicated and interrelated differentiation process. Research on the applications of omics technologies in P. aeruginosa biofilm development is ongoing, and these approaches hold great promise for expanding our knowledge of the mechanisms of biofilm formation. This review discusses the different factors that affect biofilm formation and compares P. aeruginosa biofilm formation using the omics approaches targeting essential biological macromolecules, such as DNA, RNA, Protein, and metabolome. Furthermore, we have outlined the application of currently available omics tools, such as genomics, proteomics, metabolomics, transcriptomics, and integrated multi-omics methodologies, to understand the differential gene expression (biofilm vs. planktonic bacteria) of P. aeruginosa biofilms.

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.

Challenges in using transcriptome data to study the c-di-GMP signaling network in Pseudomonas aeruginosa clinical isolates

In the Pseudomonas aeruginosa type strain PA14, 40 genes are known to encode for diguanylate cyclases (DGCs) and/or phosphodiesterases (PDEs), which modulate the intracellular pool of the nucleotide second messenger c-di-GMP. While in general, high levels of c-di-GMP drive the switch from highly motile phenotypes towards a sessile lifestyle, the different c-di-GMP modulating enzymes are responsible for smaller and in parts nonoverlapping phenotypes. In this study, we sought to utilize previously recorded P. aeruginosa gene expression datasets on 414 clinical isolates to uncover transcriptional changes as a result of a high expression of genes encoding DGCs. This approach might provide a unique opportunity to bypass the problem that for many c-di-GMP modulating enzymes it is not known under which conditions their expression is activated. However, while we demonstrate that the selection of subgroups of clinical isolates with high versus low expression of sigma factor encoding genes served the identification of their downstream regulons, we were unable to confirm the predicted DGC regulons, because the high c-di-GMP associated phenotypes were rapidly lost in the clinical isolates,. Further studies are needed to determine the specific mechanisms underlying the loss of cyclase activity upon prolonged cultivation of clinical P. aeruginosa isolates.

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.

Computationally Efficient Assembly of Pseudomonas aeruginosa Gene Expression Compendia

Thousands of Pseudomonas aeruginosa RNA sequencing (RNA-seq) gene expression profiles are publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). In this work, the transcriptional profiles from hundreds of studies performed by over 75 research groups were reanalyzed in aggregate to create a powerful tool for hypothesis generation and testing. Raw sequence data were uniformly processed using the Salmon pseudoaligner, and this read mapping method was validated by comparison to a direct alignment method. We developed filtering criteria to exclude samples with aberrant levels of housekeeping gene expression or an unexpected number of genes with no reported values and normalized the filtered compendia using the ratio-of-medians method. The filtering and normalization steps greatly improved gene expression correlations for genes within the same operon or regulon across the 2,333 samples. Since the RNA-seq data were generated using diverse strains, we report the effects of mapping samples to noncognate reference genomes by separately analyzing all samples mapped to cDNA reference genomes for strains PAO1 and PA14, two divergent strains that were used to generate most of the samples. Finally, we developed an algorithm to incorporate new data as they are deposited into the SRA. Our processing and quality control methods provide a scalable framework for taking advantage of the troves of biological information hibernating in the depths of microbial gene expression data and yield useful tools for P. aeruginosa RNA-seq data to be leveraged for diverse research goals. IMPORTANCE Pseudomonas aeruginosa is a causative agent of a wide range of infections, including chronic infections associated with cystic fibrosis. These P. aeruginosa infections are difficult to treat and often have negative outcomes. To aid in the study of this problematic pathogen, we mapped, filtered for quality, and normalized thousands of P. aeruginosa RNA-seq gene expression profiles that were publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The resulting compendia facilitate analyses across experiments, strains, and conditions. Ultimately, the workflow that we present could be applied to analyses of other microbial species.

Transcriptome Dynamics of Pseudomonas aeruginosa during Transition from Overlapping To Non-Overlapping Cell Cycles

Bacteria either duplicate their chromosome once per cell division or a new round of replication is initiated before the cells divide, thus cell cycles overlap. Here, we show that the opportunistic pathogen Pseudomonas aeruginosa switches from fast growth with overlapping cell cycles to sustained slow growth with only one replication round per cell division when cultivated under standard laboratory conditions. The transition was characterized by fast-paced, sequential changes in transcriptional activity along the ori-ter axis of the chromosome reflecting adaptation to the metabolic needs during both growth phases. Quorum sensing (QS) activity was highest at the onset of the slow growth phase with non-overlapping cell cycles. RNA sequencing of subpopulations of these cultures sorted based on their DNA content, revealed a strong gene dosage effect as well as specific expression patterns for replicating and nonreplicating cells. Expression of flagella and mexE, involved in multidrug efflux was restricted to cells that did not replicate, while those that did showed a high activity of the cell division locus and recombination genes. A possible role of QS in the formation of these subpopulations upon switching to non-overlapping cell cycles could be a subject of further research. IMPORTANCE The coordination of gene expression with the cell cycle has so far been studied only in a few bacteria, the bottleneck being the need for synchronized cultures. Here, we determined replication-associated effects on transcription by comparing Pseudomonas aeruginosa cultures that differ in their growth mode and number of replicating chromosomes. We further show that cell cycle-specific gene regulation can be principally identified by RNA sequencing of subpopulations from cultures that replicate only once per cell division and that are sorted according to their DNA content. Our approach opens the possibility to study asynchronously growing bacteria from a wide phylogenetic range and thereby enhance our understanding of the evolution of cell cycle control on the transcriptional level.

Repetitive Exposure to Bacteriophage Cocktails against Pseudomonas aeruginosa or Escherichia coli Provokes Marginal Humoral Immunity in Naïve Mice

Phage therapy of ventilator-associated pneumonia (VAP) is of great interest due to the rising incidence of multidrug-resistant bacterial pathogens. However, natural or therapy-induced immunity against therapeutic phages remains a potential concern. In this study, we investigated the innate and adaptive immune responses to two different phage cocktails targeting either Pseudomonas aeruginosa or Escherichia coli-two VAP-associated pathogens-in naïve mice without the confounding effects of a bacterial infection. Active or UV-inactivated phage cocktails or buffers were injected intraperitoneally daily for 7 days in C57BL/6J wild-type mice. Blood cell analysis, flow cytometry analysis, assessment of phage distribution and histopathological analysis of spleens were performed at 6 h, 10 days and 21 days after treatment start. Phages reached the lungs and although the phage cocktails were slightly immunogenic, phage injections were well tolerated without obvious adverse effects. No signs of activation of innate or adaptive immune cells were observed; however, both active phage cocktails elicited a minimal humoral response with secretion of phage-specific antibodies. Our findings show that even repetitive injections lead only to a minimal innate and adaptive immune response in naïve mice and suggest that systemic phage treatment is thus potentially suitable for treating bacterial lung infections.

Distinct Long- and Short-Term Adaptive Mechanisms in Pseudomonas aeruginosa

Heterogeneous environments such as the chronically infected cystic fibrosis lung drive the diversification of Pseudomonas aeruginosa populations into, e.g., mucoid, alginate-overproducing isolates or small-colony variants (SCVs). In this study, we performed extensive genome and transcriptome profiling on a clinical SCV isolate that exhibited high cyclic diguanylate (c-di-GMP) levels and a mucoid phenotype. We observed a delayed, stepwise decrease of the high levels of c-di-GMP as well as alginate gene expression upon passaging the SCV under noninducing, rich medium growth conditions over 7 days. Upon prolonged passaging, this lagging reduction of the high c-di-GMP levels under noninducing planktonic conditions (reminiscent of a hysteretic response) was followed by a phenotypic switch to a large-colony morphology, which could be linked to mutations in the Gac/Rsm signaling pathway. Complementation of the Gac/Rsm signaling-negative large-colony variants with a functional GacSA system restored the SCV colony morphotype but was not able to restore the high c-di-GMP levels of the SCV. Our data thus suggest that expression of the SCV colony morphotype and modulation of c-di-GMP levels are genetically separable and follow different evolutionary paths. The delayed switching of c-di-GMP levels in response to fluctuating environmental conditions might provide a unique opportunity to include a time dimension to close the gap between short-term phenotypic and long-term genetic adaptation to biofilm-associated growth conditions. IMPORTANCE Extreme environments, such as those encountered during an infection process in the human host, make effective bacterial adaptation inevitable. While bacteria adapt individually by activating stress responses, long-term adaptation of bacterial communities to challenging conditions can be achieved via genetic fixation of favorable traits. In this study, we describe a two-pronged bacterial stress resistance strategy in the opportunistic pathogen Pseudomonas aeruginosa. We show that the production of adjusted elevated c-di-GMP levels, which drive protected biofilm-associated phenotypes in vivo, resembles a stable hysteretic response which prevents unwanted frequent switching. Cellular hysteresis might provide a link between individual adaptability and evolutionary adaptation to ensure the evolutionary persistence of host-adapted stress response strategies.

Effects of the Sex Steroid Hormone Estradiol on Biofilm Growth of Cystic Fibrosis Pseudomonas aeruginosa Isolates

Women with cystic fibrosis (CF) have a significantly lower life expectancy compared to men, which is indicated by an earlier impairment of lung function due to chronic colonization with biofilm formed by Pseudomonas aeruginosa. There is growing evidence that blood serum concentrations of the steroid sex hormone estradiol (E₂) correlate with the occurrence of pulmonary exacerbations in CF but also play a role in the mucoid switch of P. aeruginosa. This study aims to shed light on possible microbiological reasons for sexual dimorphism in CF by investigating the influence of E₂ on biofilm formation of P. aeruginosa CF isolates. For this purpose, 10 CF isolates of the respiratory tract derived from different CF patients have been treated with E₂ in a microtiter plate biofilm model. Biofilms have been examined by crystal violet assays, field emission scanning electron microscopy (FE-SEM), 3D laser scanning microscopy (LSM), and quorum sensing (QS) reporter assays of the supernatants taken from biofilms. This allowed us to simultaneously investigate the effects of E₂ on attached biofilm mass, biofilm ultrastructure, and QS activity. Upon E₂ treatment, six out of 10 investigated CF isolates showed an increase of attached biofilm mass, whereas biofilms from two tested non-CF laboratory strains (PAO1 and ATCC19660) did not. Moreover, FE-SEM and 3D LSM analyses of the E₂ responsive CF biofilms revealed ultrastructural remodeling of biofilm structure at different scales with increased formation of prominent biofilm spots, enhanced coverage with extracellular polymeric substance (EPS), and extended average surface roughness. QS activity measurements performed in biofilm supernatants via luminescence acyl homoserine lactone (AHL) reporter assays further showed that E₂ treatment may also modulate QS signaling, as shown in an E₂ sensitive CF isolate. Together, our results suggest the biofilm modulating effects of E₂ on various clinical CF isolates that are documented by both biomass and ultrastructural changes of biofilms. The gained new insight into the influence of steroid hormones on P. aeruginosa biofilm phenotypes might pave the way for novel future approaches in personalized medicine based on the patients’ sex and hormonal status.

Evolution of Habitat-Dependent Antibiotic Resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that usually causes difficult-to-treat infections due to its low intrinsic antibiotic susceptibility and outstanding capacity for becoming resistant to antibiotics. In addition, it has a remarkable metabolic versatility, being able to grow in different habitats, from natural niches to different and changing inpatient environments. Study of the environmental conditions that shape genetic and phenotypic changes of P. aeruginosa toward antibiotic resistance supposes a novelty, since experimental evolution assays are usually performed with well-defined antibiotics in regular laboratory growth media. Therefore, in this work we address the extent to which the nutrients’ availability may constrain the evolution of antibiotic resistance. We determined that P. aeruginosa genetic trajectories toward resistance to tobramycin, ceftazidime, and ceftazidime-avibactam are different when evolving in laboratory rich medium, urine, or synthetic sputum. Furthermore, our study, linking genotype with phenotype, showed a clear impact of each analyzed environment on both the fitness and resistance level associated with particular resistance mutations. This indicates that the phenotype associated with specific resistance mutations is variable and dependent on the bacterial metabolic state in each particular habitat. Our results support that the design of evolution-based strategies to tackle P. aeruginosa infections should be based on robust patterns of evolution identified within each particular infection and body location.

IMPORTANCE Predicting evolution toward antibiotic resistance (AR) and its associated trade-offs, such as collateral sensitivity, is important to design evolution-based strategies to tackle AR. However, the effect of nutrients' availability on such evolution, particularly those that can be found under in vivo infection conditions, has been barely addressed. We analyzed the evolutionary patterns of P. aeruginosa in the presence of antibiotics in different media, including urine and synthetic sputum, whose compositions are similar to the ones in infections, finding that AR evolution differs, depending on growth conditions. Furthermore, the representative mutants isolated under each condition tested render different AR levels and fitness costs, depending on nutrients’ availability, supporting the idea that environmental constraints shape the phenotypes associated with specific AR mutations. Consequently, the selection of AR mutations that render similar phenotypes is environment dependent. The analysis of evolution patterns toward AR requires studying growth conditions mimicking those that bacteria face during in vivo evolution.

Pseudomonas aeruginosa Is More Tolerant Under Biofilm Than Under Planktonic Growth Conditions: A Multi-Isolate Survey

Biofilm-associated bacteria exhibit profound changes in bacterial physiology. They thrive in the environment but also in the human host in protected sessile communities. Antimicrobial therapy usually fails, despite the absence of genotypic resistance, and it is commonly accepted that biofilm-grown bacteria are up to 1,000-fold more resistant than planktonic cells. We are only at the beginning to understand the reasons for biofilm recalcitrance, and systematic approaches to describe biofilm-induced tolerance phenotypes are lacking. In this study, we investigated a large and highly diverse collection of 352 clinical Pseudomonas aeruginosa isolates for their antimicrobial susceptibility profiles under biofilm growth conditions towards the antibiotics ciprofloxacin, tobramycin, and colistin. We discovered characteristic patterns of drug-specific killing activity and detected conditional tolerance levels far lower (in the range of the minimal inhibitory concentration (MIC)), but also far higher (up to 16,000-fold increase compared to planktonic cells) than generally believed. This extremely broad distribution of biofilm-induced tolerance phenotypes across the clinical isolates was greatly influenced by the choice of the antibiotic. We furthermore describe cross-tolerance against ciprofloxacin and tobramycin, but not colistin, and observed an additive activity between biofilm-induced tolerance and genetically determined resistance. This became less evident when the biofilm-grown cells were exposed to very high antibiotic concentrations. Although much more remains to be learned on the molecular mechanisms underlying biofilm-induced tolerance, our data on intra-species variations in tolerance profiles provide valuable new insights. Furthermore, our observation that colistin appears to act independently of the tolerance mechanisms of individual clinical strains could make colistin a valuable therapeutic option in chronic biofilm-associated infections characterized by the presence of particularly tolerant strains.

Evolution of biofilm-adapted gene expression profiles in lasR-deficient clinical Pseudomonas aeruginosa isolates

The overall success of a pathogenic microbe depends on its ability to efficiently adapt to challenging conditions in the human host. Long-term evolution experiments track and predict adaptive trajectories and have contributed significantly to our understanding of the driving forces of bacterial adaptation. In this study, we conducted a cross-sectional study instead of long-term longitudinal evolution experiments. We analyzed the transcriptional profiles as well as genomic sequence variations of a large number of clinical Pseudomonas aeruginosa isolates that have been recovered from different infected human sites. Convergent changes in gene expression patterns were found in different groups of clinical isolates. The majority of repeatedly observed expression patterns could be attributed to a defective lasR gene, which encodes the major quorum-sensing regulator LasR. Strikingly, the gene expression pattern of the lasR-defective strains appeared to reflect a transcriptional response that evolves in a direction consistent with growth within a biofilm. In a process of genetic assimilation, lasR-deficient P. aeruginosa isolates appear to constitutively express a biofilm-adapted transcriptional profile and no longer require a respective environmental trigger. Our results demonstrate that profiling the functional consequences of pathoadaptive mutations in clinical isolates reveals long-term evolutionary pathways and may explain the success of lasR mutants in the opportunistic pathogen P. aeruginosa in a clinical context.

The Pseudomonas aeruginosa whole genome sequence: A 20th anniversary celebration

Toward the end of August 2000, the 6.3 Mbp whole genome sequence of Pseudomonas aeruginosa strain PAO1 was published. With 5570 open reading frames (ORFs), PAO1 had the largest microbial genome sequenced up to that point in time-including a large proportion of metabolic, transport and antimicrobial resistance genes supporting its ability to colonize diverse environments. A remarkable 9% of its ORFs were predicted to encode proteins with regulatory functions, providing new insight into bacterial network complexity as a function of network size. In this celebratory article, we fast forward 20 years, and examine how access to this resource has transformed our understanding of P. aeruginosa. What follows is more than a simple review or commentary; we have specifically asked some of the leaders in the field to provide personal reflections on how the PAO1 genome sequence, along with the Pseudomonas Community Annotation Project (PseudoCAP) and Pseudomonas Genome Database (pseudomonas.com), have contributed to the many exciting discoveries in this field. In addition to bringing us all up to date with the latest developments, we also ask our contributors to speculate on how the next 20 years of Pseudomonas research might pan out.

Single-Nucleotide Polymorphism-Based Genetic Diversity Analysis of Clinical Pseudomonas aeruginosa Isolates

Extensive use of next-generation sequencing has the potential to transform our knowledge on how genomic variation within bacterial species impacts phenotypic versatility. Because different environments have unique selection pressures, they drive divergent evolution. However, there is also parallel or convergent evolution of traits in independent bacterial isolates inhabiting similar environments. The application of tools to describe population-wide genomic diversity provides an opportunity to measure the predictability of genetic changes underlying adaptation. Here, we describe patterns of sequence variations in the core genome among 99 individual Pseudomonas aeruginosa clinical isolates and identified single-nucleotide polymorphisms that are the basis for branching of the phylogenetic tree. We also identified single-nucleotide polymorphisms that were acquired independently, in separate lineages, and not through inheritance from a common ancestor. Although our results demonstrate that the Pseudomonas aeruginosa core genome is highly conserved and in general, not subject to adaptive evolution, instances of parallel evolution will provide an opportunity to uncover genetic changes that underlie phenotypic diversity.

Targeting Bacterial Gyrase with Cystobactamid, Fluoroquinolone, and Aminocoumarin Antibiotics Induces Distinct Molecular Signatures in Pseudomonas aeruginosa

The design of novel antibiotics relies on a profound understanding of their mechanism of action. While it has been shown that cellular effects of antibiotics cluster according to their molecular targets, we investigated whether compounds binding to different sites of the same target can be differentiated by their transcriptome or metabolome signatures. The effects of three fluoroquinolones, two aminocoumarins, and two cystobactamids, all inhibiting bacterial gyrase, on Pseudomonas aeruginosa at subinhibitory concentrations could be distinguished clearly by RNA sequencing as well as metabolomics. We observed a strong (2.8- to 212-fold) induction of autolysis-triggering pyocins in all gyrase inhibitors, which correlated with extracellular DNA (eDNA) release. Gyrase B-binding aminocoumarins induced the most pronounced changes, including a strong downregulation of phenazine and rhamnolipid virulence factors. Cystobactamids led to a downregulation of a glucose catabolism pathway. The study implies that clustering cellular mechanisms of action according to the primary target needs to take class-dependent variances into account.

IMPORTANCE Novel antibiotics are urgently needed to tackle the growing worldwide problem of antimicrobial resistance. Bacterial pathogens possess few privileged targets for a successful therapy: the majority of existing antibiotics as well as current candidates in development target the complex bacterial machinery for cell wall synthesis, protein synthesis, or DNA replication. An important mechanistic question addressed by this study is whether inhibiting such a complex target at different sites with different compounds has similar or differentiated cellular consequences. Using transcriptomics and metabolomics, we demonstrate that three different classes of gyrase inhibitors can be distinguished by their molecular signatures in P. aeruginosa. We describe the cellular effects of a promising, recently identified gyrase inhibitor class, the cystobactamids, in comparison to those of the established gyrase A-binding fluoroquinolones and the gyrase B-binding aminocoumarins. The study results have implications for mode-of-action discovery approaches based on target-specific reference compounds, as they highlight the intraclass variability of cellular compound effects.

RNA atlas of human bacterial pathogens uncovers stress dynamics linked to infection

Bacterial processes necessary for adaption to stressful host environments are potential targets for new antimicrobials. Here, we report large-scale transcriptomic analyses of 32 human bacterial pathogens grown under 11 stress conditions mimicking human host environments. The potential relevance of the in vitro stress conditions and responses is supported by comparisons with available in vivo transcriptomes of clinically important pathogens. Calculation of a probability score enables comparative cross-microbial analyses of the stress responses, revealing common and unique regulatory responses to different stresses, as well as overlapping processes participating in different stress responses. We identify conserved and species-specific 'universal stress responders', that is, genes showing altered expression in multiple stress conditions. Non-coding RNAs are involved in a substantial proportion of the responses. The data are collected in a freely available, interactive online resource (PATHOgenex).

Compensatory evolution of Pseudomonas aeruginosa's slow growth phenotype suggests mechanisms of adaptation in cystic fibrosis

Long-term infection of the airways of cystic fibrosis patients with Pseudomonas aeruginosa is often accompanied by a reduction in bacterial growth rate. This reduction has been hypothesised to increase within-patient fitness and overall persistence of the pathogen. Here, we apply adaptive laboratory evolution to revert the slow growth phenotype of P. aeruginosa clinical strains back to a high growth rate. We identify several evolutionary trajectories and mechanisms leading to fast growth caused by transcriptional and mutational changes, which depend on the stage of adaptation of the strain. Return to high growth rate increases antibiotic susceptibility, which is only partially dependent on reversion of mutations or changes in the transcriptional profile of genes known to be linked to antibiotic resistance. We propose that similar mechanisms and evolutionary trajectories, in reverse direction, may be involved in pathogen adaptation and the establishment of chronic infections in the antibiotic-treated airways of cystic fibrosis patients.

The Genomic Basis of Rapid Adaptation to Antibiotic Combination Therapy in Pseudomonas aeruginosa

Combination therapy is a common antibiotic treatment strategy that aims at minimizing the risk of resistance evolution in several infectious diseases. Nonetheless, evidence supporting its efficacy against the nosocomial opportunistic pathogen Pseudomonas aeruginosa remains elusive. Identification of the possible evolutionary paths to resistance in multidrug environments can help to explain treatment outcome. For this purpose, we here performed whole-genome sequencing of 127 previously evolved populations of P. aeruginosa adapted to sublethal doses of distinct antibiotic combinations and corresponding single-drug treatments, and experimentally characterized several of the identified variants. We found that alterations in the regulation of efflux pumps are the most favored mechanism of resistance, regardless of the environment. Unexpectedly, we repeatedly identified intergenic variants in the adapted populations, often with no additional mutations and usually associated with genes involved in efflux pump expression, possibly indicating a regulatory function of the intergenic regions. The experimental analysis of these variants demonstrated that the intergenic changes caused similar increases in resistance against single and multidrug treatments as those seen for efflux regulatory gene mutants. Surprisingly, we could find no substantial fitness costs for a majority of these variants, most likely enhancing their competitiveness toward sensitive cells, even in antibiotic-free environments. We conclude that the regulation of efflux is a central target of antibiotic-mediated selection in P. aeruginosa and that, importantly, changes in intergenic regions may represent a usually neglected alternative process underlying bacterial resistance evolution, which clearly deserves further attention in the future.

Phenotypic and Transcriptomic Analyses of Seven Clinical Stenotrophomonas maltophilia Isolates Identify a Small Set of Shared and Commonly Regulated Genes Involved in the Biofilm Lifestyle

Stenotrophomonas maltophilia is one of the most frequently isolated multidrug-resistant nosocomial opportunistic pathogens. It contributes to disease progression in cystic fibrosis (CF) patients and is frequently isolated from wounds, infected tissues, and catheter surfaces. On these diverse surfaces S. maltophilia lives in single-species or multispecies biofilms. Since very little is known about common processes in biofilms of different S. maltophilia isolates, we analyzed the biofilm profiles of 300 clinical and environmental isolates from Europe of the recently identified main lineages Sgn3, Sgn4, and Sm2 to Sm18. The analysis of the biofilm architecture of 40 clinical isolates revealed the presence of multicellular structures and high phenotypic variability at a strain-specific level. Further, transcriptome analyses of biofilm cells of seven clinical isolates identified a set of 106 shared strongly expressed genes and 33 strain-specifically expressed genes. Surprisingly, the transcriptome profiles of biofilm versus planktonic cells revealed that just 9.43% ± 1.36% of all genes were differentially regulated. This implies that just a small set of shared and commonly regulated genes is involved in the biofilm lifestyle. Strikingly, iron uptake appears to be a key factor involved in this metabolic shift. Further, metabolic analyses implied that S. maltophilia employs a mostly fermentative growth mode under biofilm conditions. The transcriptome data of this study together with the phenotypic and metabolic analyses represent so far the largest data set on S. maltophilia biofilm versus planktonic cells. This study will lay the foundation for the identification of strategies for fighting S. maltophilia biofilms in clinical and industrial settings.IMPORTANCE Microorganisms living in a biofilm are much more tolerant to antibiotics and antimicrobial substances than planktonic cells are. Thus, the treatment of infections caused by microorganisms living in biofilms is extremely difficult. Nosocomial infections (among others) caused by S. maltophilia, particularly lung infection among CF patients, have increased in prevalence in recent years. The intrinsic multidrug resistance of S. maltophilia and the increased tolerance to antimicrobial agents of its biofilm cells make the treatment of S. maltophilia infection difficult. The significance of our research is based on understanding the common mechanisms involved in biofilm formation of different S. maltophilia isolates, understanding the diversity of biofilm architectures among strains of this species, and identifying the differently regulated processes in biofilm versus planktonic cells. These results will lay the foundation for the treatment of S. maltophilia biofilms.

Organism-specific depletion of highly abundant RNA species from bacterial total RNA

High-throughput sequencing has become a standard tool for transcriptome analysis. The depletion of overrepresented RNA species from sequencing libraries plays a key role in establishing potent and cost-efficient RNA-seq routines. Commercially available kits are known to obtain good results for the reduction of ribosomal RNA (rRNA). However, we found that the transfer-messenger RNA (tmRNA) was frequently highly abundant in rRNA-depleted samples of Pseudomonas aeruginosa, consuming up to 25 % of the obtained reads. The tmRNA fraction was particularly high in samples taken from stationary cultures. This suggests that overrepresentation of this RNA species reduces the mRNA fraction when cells are grown under challenging conditions. Here, we present an RNase-H-based depletion protocol that targets the tmRNA in addition to ribosomal RNAs. We were able to increase the mRNA fraction to 93–99% and therefore outperform not only the commercially Ribo-off kit (Vazyme) operating by the same principle but also the formerly widely used Ribo-Zero kit (Illumina). Maximizing the read share of scientifically interesting RNA species enhances the discriminatory potential of next-generation RNA-seq experiments and, therefore, can contribute to a better understanding of the transcriptomic landscape of bacterial pathogens and their used mechanisms in host infection.

Host-induced spermidine production in motile Pseudomonas aeruginosa triggers phagocytic uptake

Exploring the complexity of host-pathogen communication is vital to understand why microbes persist within a host, while others are cleared. Here, we employed a dual-sequencing approach to unravel conversational turn-taking of dynamic host-pathogen communications. We demonstrate that upon hitting a host cell, motile Pseudomonas aeruginosa induce a specific gene expression program. This results in the expression of spermidine on the surface, which specifically activates the PIP3-pathway to induce phagocytic uptake into primary or immortalized murine cells. Non-motile bacteria are more immunogenic due to a lower expression of arnT upon host-cell contact, but do not produce spermidine and are phagocytosed less. We demonstrate that not only the presence of pathogen inherent molecular patterns induces immune responses, but that bacterial motility is linked to a host-cell-induced expression of additional immune modulators. Our results emphasize on the value of integrating microbiological and immunological findings to unravel complex and dynamic host-pathogen interactions.

Untargeted LC-MS Metabolomics Differentiates Between Virulent and Avirulent Clinical Strains of Pseudomonas aeruginosa

Pseudomonas aeruginosa is a facultative pathogen that can cause, inter alia, acute or chronic pneumonia in predisposed individuals. The gram-negative bacterium displays considerable genomic and phenotypic diversity that is also shaped by small molecule secondary metabolites. The discrimination of virulence phenotypes is highly relevant to the diagnosis and prognosis of P. aeruginosa infections. In order to discover small molecule metabolites that distinguish different virulence phenotypes of P. aeruginosa, 35 clinical strains were cultivated under standard conditions, characterized in terms of virulence and biofilm phenotype, and their metabolomes were investigated by untargeted liquid chromatography-mass spectrometry. The data was both mined for individual candidate markers as well as used to construct statistical models to infer the virulence phenotype from metabolomics data. We found that clinical strains that differed in their virulence and biofilm phenotype also had pronounced divergence in their metabolomes, as underlined by 332 features that were significantly differentially abundant with fold changes greater than 1.5 in both directions. Important virulence-associated secondary metabolites like rhamnolipids, alkyl quinolones or phenazines were found to be strongly upregulated in virulent strains. In contrast, we observed little change in primary metabolism. A hitherto novel cationic metabolite with a sum formula of C₁₂H₁₅N₂ could be identified as a candidate biomarker. A random forest model was able to classify strains according to their virulence and biofilm phenotype with an area under the Receiver Operation Characteristics curve of 0.84. These findings demonstrate that untargeted metabolomics is a valuable tool to characterize P. aeruginosa virulence, and to explore interrelations between clinically important phenotypic traits and the bacterial metabolome.

There is no hiding if you Seq: recent breakthroughs in Pseudomonas aeruginosa research revealed by genomic and transcriptomic next-generation sequencing

The advent of next-generation sequencing technology has revolutionized the field of prokaryotic genetics and genomics by allowing interrogation of entire genomes, transcriptomes and global transcription factor binding profiles. As more studies employing these techniques have been performed, the state of the art regarding prokaryotic gene regulation has developed from the level of individual genes to genetic regulatory networks and systems biology. When applied to bacterial pathogens, particularly valuable insights have been gained into the regulation of virulence-associated genes, their relative importance to bacterial survival in planktonic, biofilm or host infection scenarios, antimicrobial resistance and the molecular details of host-pathogen interactions. This review outlines some of the latest developments and applications of next-generation sequencing techniques that have used primarily Pseudomonas aeruginosa as a model system. We focus particularly on insights into Pseudomonas virulence and infection that have been gained from these approaches and the future directions in which this field could develop.

Parallel evolutionary paths to produce more than one Pseudomonas aeruginosa biofilm phenotype

Studying parallel evolution of similar traits in independent within-species lineages provides an opportunity to address evolutionary predictability of molecular changes underlying adaptation. In this study, we monitored biofilm forming capabilities, motility, and virulence phenotypes of a plethora of phylogenetically diverse clinical isolates of the opportunistic pathogen Pseudomonas aeruginosa. We also recorded biofilm-specific and planktonic transcriptional responses. We found that P. aeruginosa isolates could be stratified based on the production of distinct organismal traits. Three major biofilm phenotypes, which shared motility and virulence phenotypes, were produced repeatedly in several isolates, indicating that the phenotypes evolved via parallel or convergent evolution. Of note, while we found a restricted general response to the biofilm environment, the individual groups of biofilm phenotypes reproduced biofilm transcriptional profiles that included the expression of well-known biofilm features, such as surface adhesive structures and extracellular matrix components. Our results provide insights into distinct ways to make a biofilm and indicate that genetic adaptations can modulate multiple pathways for biofilm development that are followed by several independent clinical isolates. Uncovering core regulatory pathways that drive biofilm-associated growth and tolerance towards environmental stressors promises to give clues to host and environmental interactions and could provide useful targets for new clinical interventions.

Adapting to the Airways: Metabolic Requirements of Pseudomonas aeruginosa during the Infection of Cystic Fibrosis Patients

Pseudomonas aeruginosa is one of the major causes of morbidity and mortality of cystic fibrosis patients. During the infection, the bacteria colonize the nutritional rich lung mucus, which is present in the airway secretions in the patients, and they adapt their phenotype accordingly to the lung environment. In the airways, P. aeruginosa undergoes a broad metabolic rewiring as a consequence of the nutritional and stressful complexity of the lungs. However, the role of such metabolic rewiring on the infection outcome is poorly understood. Here, we review the metabolic evolution of clinical strains of P. aeruginosa during a cystic fibrosis lung infection and the metabolic functions operating in vivo under patho-physiological conditions. Finally, we discuss the perspective of modeling the cystic fibrosis environment using genome scale metabolic models of P. aeruginosa. Understanding the physiological changes occurring during the infection may pave the way to a more effective treatment for P. aeruginosa lung infections.

The Core Proteome of Biofilm-Grown Clinical Pseudomonas aeruginosa Isolates

Comparative genomics has greatly facilitated the identification of shared as well as unique features among individual cells or tissues, and thus offers the potential to find disease markers. While proteomics is recognized for its potential to generate quantitative maps of protein expression, comparative proteomics in bacteria has been largely restricted to the comparison of single cell lines or mutant strains. In this study, we used a data independent acquisition (DIA) technique, which enables global protein quantification of large sample cohorts, to record the proteome profiles of overall 27 whole genome sequenced and transcriptionally profiled clinical isolates of the opportunistic pathogen Pseudomonas aeruginosa. Analysis of the proteome profiles across the 27 clinical isolates grown under planktonic and biofilm growth conditions led to the identification of a core biofilm-associated protein profile. Furthermore, we found that protein-to-mRNA ratios between different P. aeruginosa strains are well correlated, indicating conserved patterns of post-transcriptional regulation. Uncovering core regulatory pathways, which drive biofilm formation and associated antibiotic tolerance in bacterial pathogens, promise to give clues to interactions between bacterial species and their environment and could provide useful targets for new clinical interventions to combat biofilm-associated infections.

Establishment of an induced memory response in Pseudomonas aeruginosa during infection of a eukaryotic host

In a given habitat, bacterial cells often experience recurrent exposures to the same environmental stimulus. The ability to memorize the past event and to adjust current behaviors can lead to efficient adaptation to the recurring stimulus. Here we demonstrate that the versatile bacterium Pseudomonas aeruginosa adopts a virulence phenotype after serial passage in the invertebrate model host Galleria mellonella. The virulence phenotype was not linked to the acquisition of genetic variations and was sustained for several generations, despite cultivation of the ex vivo virulence-adapted P. aeruginosa cells under rich medium conditions in vitro. Transcriptional reprogramming seemed to be induced by a host-specific food source, as reprogramming was also observed upon cultivation of P. aeruginosa in rich medium supplemented with polyunsaturated long-chain fatty acids. The establishment of induced memory responses adds a time dimension and seems to fill the gap between long-term evolutionary genotypic adaptation and short-term induced individual responses. Efforts to unravel the fundamental mechanisms that underlie the carry-over effect to induce such memory responses will continue to be of importance as hysteretic behavior can serve survival of bacterial populations in changing and challenging habitats.

Antibiotic resistance in Pseudomonas aeruginosa - Mechanisms, epidemiology and evolution

Antibiotics are powerful drugs used in the treatment of bacterial infections. The inappropriate use of these medicines has driven the dissemination of antibiotic resistance (AR) in most bacteria. Pseudomonas aeruginosa is an opportunistic pathogen commonly involved in environmental- and difficult-to-treat hospital-acquired infections. This species is frequently resistant to several antibiotics, being in the "critical" category of the WHO's priority pathogens list for research and development of new antibiotics. In addition to a remarkable intrinsic resistance to several antibiotics, P. aeruginosa can acquire resistance through chromosomal mutations and acquisition of AR genes. P. aeruginosa has one of the largest bacterial genomes and possesses a significant assortment of genes acquired by horizontal gene transfer (HGT), which are frequently localized within integrons and mobile genetic elements (MGEs), such as transposons, insertion sequences, genomic islands, phages, plasmids and integrative and conjugative elements (ICEs). This genomic diversity results in a non-clonal population structure, punctuated by specific clones that are associated with significant morbidity and mortality worldwide, the so-called high-risk clones. Acquisition of MGEs produces a fitness cost in the host, that can be eased over time by compensatory mutations during MGE-host coevolution. Even though plasmids and ICEs are important drivers of AR, the underlying evolutionary traits that promote this dissemination are poorly understood. In this review, we provide a comprehensive description of the main strategies involved in AR in P. aeruginosa and the leading drivers of HGT in this species. The most recently developed genomic tools that allowed a better understanding of the features contributing for the success of P. aeruginosa are discussed.

Non-invasive, ratiometric determination of intracellular pH in Pseudomonas species using a novel genetically encoded indicator

The ability of Pseudomonas species to thrive in all major natural environments (i.e. terrestrial, freshwater and marine) is based on its exceptional capability to adapt to physicochemical changes. Thus, environmental bacteria have to tightly control the maintenance of numerous physiological traits across different conditions. The intracellular pH (pHi ) homoeostasis is a particularly important feature, since the pHi influences a large portion of the biochemical processes in the cell. Despite its importance, relatively few reliable, easy-to-implement tools have been designed for quantifying in vivo pHi changes in Gram-negative bacteria with minimal manipulations. Here we describe a convenient, non-invasive protocol for the quantification of the pHi in bacteria, which is based on the ratiometric fluorescent indicator protein PHP (pH indicator for Pseudomonas). The DNA sequence encoding PHP was thoroughly adapted to guarantee optimal transcription and translation of the indicator in Pseudomonas species. Our PHP-based quantification method demonstrated that pHi is tightly regulated over a narrow range of pH values not only in Pseudomonas, but also in other Gram-negative bacterial species such as Escherichia coli. The maintenance of the cytoplasmic pH homoeostasis in vivo could also be observed upon internal (e.g. redirection of glucose consumption pathways in P. putida) and external (e.g. antibiotic exposure in P. aeruginosa) perturbations, and the PHP indicator was also used to follow dynamic changes in the pHi upon external pH shifts. In summary, our work describes a reliable method for measuring pHi in Pseudomonas, allowing for the detailed investigation of bacterial pHi homoeostasis and its regulation.

Core regulon of the global anaerobic regulator Anr targets central metabolism functions in Pseudomonas species

A comparative genome analysis of the global anaerobic regulator Anr regulon in five species of Pseudomonas with different life style was performed. Expression of this regulator was detected in all analyzed Pseudomonas. The predicted Anr regulon (pan-regulon) consisted of 253 genes. However, only 11 Anr-boxes located upstream of qor/hemF, hemN, cioA/PA3931, azu, rpsL, gltP, orthologous to PA2867, cspD, tyrZ, slyD and oprG, were common to all species. Whole genome in silico prediction of metabolic pathways identified genes belonging to heme biosynthesis, cytochromes and Entner-Doudoroff pathway as members of Anr regulon in all strains. Extending genome analysis to 28 Pseudomonas spp. spanning all phylogenetic groups showed Anr-boxes conservation in genes related to these functions. When present, genes related to anaerobic metabolism were predicted to hold Anr-boxes. Focused on the genomes of eight P. aeruginosa isolates of diverse origins, we observed a conserved regulon, sharing nearly 80% of the genes, indicating its key role in this opportunistic pathogen. The results suggest that the core Anr regulon comprises genes involved in central metabolism and aerobic electron transport chain, whereas those genes related to anaerobic metabolism and other functions constitute the accessory Anr-regulon, thereby differentially contributing to the ecological fitness of each Pseudomonas species.

Biofilms 2018: A diversity of microbes and mechanisms

The 8th ASM Conference on Biofilms was held in Washington D.C. on October 7-11, 2018. This very highly subscribed meeting represented a wide breadth of current research in biofilms, and included over 500 attendees, 12 sessions with 64 oral presentations, and four poster sessions with about 400 posters.