Genomic and transcriptomic characterization of Pseudomonas aeruginosa small colony variants derived from a chronic infection model

Genomic Comparative of Pseudomonas aeruginosa Small Colony Variant, Mucoid and Non-mucoid Phenotypes Obtained from a Patient with Cystic Fibrosis During Respiratory Exacerbations

Pseudomonas aeruginosa, the most prevalent opportunistic pathogen in chronic obstructive pulmonary disease, associated with high morbidity and mortality in patients with cystic fibrosis (CF), is practically impossible to be eradicated from the airways in chronicity. Its extraordinary genomic plasticity is possibly associated with high antimicrobial resistance, virulence factors, and its phenotypic diversity. The occurrence of P. aeruginosa isolates promoting airway infection, showing mucoid, non-mucoid, and small colony variant (SCV) phenotypes, was observed simultaneously, in the present study, in sputum cultures obtained from a male CF young patient with chronic pulmonary infection for over a decade. The isolates belonged to a new ST (2744) were obtained in two moments of exacerbation of the respiratory disease, in which he was hospitalized. Genetic background and phenotypic analysis indicated that the isolates exhibited multi- and pan-antimicrobial resistant profiles, as well as non-susceptible to polymyxin and predominantly hypermutable (HPM) phenotypes. Whole genome sequencing showed variations in genome sizes, coding sequences and their determinants of resistance and virulence. The annotated genomes were compared for antimicrobial resistance, hypermutability, and SCV characteristics. We highlight the lack of reported genetic determinants of SCV emergence and HPM phenotypes, which can be explained in part due to the very short time between collections of isolates. To the best of our knowledge, this is the first report of genome sequencing of P. aeruginosa SCV from a CF patient in Brazil.

Structural genome variants of Pseudomonas aeruginosa clone C and PA14 strains

Plasticity of Pseudomonas aeruginosa chromosomes is mainly driven by an extended accessory genome that is shaped by insertion and deletion events. Further modification of the genome composition can be induced by chromosomal inversion events which lead to relocation of genes in the affected genomic DNA segments, modify the otherwise highly conserved core genome synteny and could even alter the location of the replication terminus. Although the genome of the first sequenced strain, PAO1, displayed such a large genomic inversion, knowledge on such recombination events in the P. aeruginosa population is limited. Several large inversions had been discovered in the late 1990s in cystic fibrosis isolates of the major clonal lineage C by physical genome mapping, and subsequent work on these examples led to the characterization of the DNA at the recombination breakpoints and a presumed recombination mechanism. Since then, the topic was barely addressed in spite of the compilation of thousands of P. aeruginosa genome sequences that are deposited in databases. Due to the use of second-generation sequencing, genome contig assembly had usually followed synteny blueprints provided by the existing reference genome sequences. Inversion detection was not feasible by these approaches, as the respective read lengths did not allow reliable resolution of sequence repeats that are typically found at the borders of inverted segments. In this study, we applied PacBio and MinION long-read sequencing to isolates of the mentioned clone C collection. Confirmation of inversions predicted from the physical mapping data demonstrated that unbiased sequence assembly of such read datasets allows the detection of genomic inversions and the resolution of the recombination breakpoint regions. Additional long-read sequencing of representatives of the other major clonal lineage, PA14, revealed large inversions in several isolates, from cystic fibrosis origin as well as from other sources. These findings indicated that inversion events are not restricted to strains from chronic infection background, but could be widespread in the P. aeruginosa population and contribute to genome plasticity. Moreover, the monitored examples emphasized the role of small mobile DNA units, such as IS elements or transposons, and accessory DNA elements in the inversion-related recombination processes.

Ribosomal RNA operons define a central functional compartment in the Streptomyces chromosome

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. These bacteria possess a large linear chromosome genetically compartmentalized: core genes are grouped in the central part, while terminal regions are populated by poorly conserved genes. In exponentially growing cells, chromosome conformation capture unveiled sharp boundaries formed by ribosomal RNA (rrn) operons that segment the chromosome into multiple domains. Here we further explore the link between the genetic distribution of rrn operons and Streptomyces genetic compartmentalization. A large panel of genomes of species representative of the genus diversity revealed that rrn operons and core genes form a central skeleton, the former being identifiable from their core gene environment. We implemented a new nomenclature for Streptomyces genomes and trace their rrn-based evolutionary history. Remarkably, rrn operons are close to pericentric inversions. Moreover, the central compartment delimited by rrn operons has a very dense, nearly invariant core gene content. Finally, this compartment harbors genes with the highest expression levels, regardless of gene persistence and distance to the origin of replication. Our results highlight that rrn operons are structural boundaries of a central functional compartment prone to transcription in Streptomyces.

Pseudomonas aeruginosa Strains from Both Clinical and Environmental Origins Readily Adopt a Stable Small-Colony-Variant Phenotype Resulting from Single Mutations in c-di-GMP Pathways

A subpopulation of small-colony variants (SCVs) is a frequently observed feature of Pseudomonas aeruginosa isolates obtained from colonized cystic fibrosis lungs. Since most SCVs have until now been isolated from clinical samples, it remains unclear how widespread the ability of P. aeruginosa strains to develop this phenotype is and what the genetic mechanism(s) behind the emergence of SCVs are according to the origin of the isolate. In the present work, we investigated the ability of 22 P. aeruginosa isolates from various environmental origins to spontaneously adopt an SCV-like smaller alternative morphotype distinguishable from that of the ancestral parent strain under laboratory culture conditions. We found that all the P. aeruginosa strains tested could adopt an SCV phenotype, regardless of their origin. Whole-genome sequencing of SCVs obtained from clinical and environmental sources revealed single mutations exclusively in two distinct c-di-GMP signaling pathways, the Wsp and YfiBNR pathways. We conclude that the ability to switch to an SCV phenotype is a conserved feature of P. aeruginosa and results from the acquisition of a stable genetic mutation, regardless of the origin of the strain. IMPORTANCE P. aeruginosa is an opportunistic pathogen that thrives in many environments. It poses a significant health concern, notably because this bacterium is the most prevalent pathogen found in the lungs of people with cystic fibrosis. In infected hosts, its persistence is considered related to the emergence of an alternative small-colony-variant (SCV) phenotype. By reporting the distribution of P. aeruginosa SCVs in various nonclinical environments and the involvement of c-di-GMP in SCV emergence from both clinical and environmental strains, this work contributes to understanding a conserved adaptation mechanism used by P. aeruginosa to adapt readily in all environments. Hindering this adaptation strategy could help control persistent infection by P. aeruginosa.

Mutation to ispA Produces Stable Small-Colony Variants of Pseudomonas aeruginosa That Have Enhanced Aminoglycoside Resistance

Pseudomonas aeruginosa is a major pathogen in burn wound infections. We present one of the first reports of small-colony variant (SCV) emergence of P. aeruginosa, taken from a patient under aminoglycosides for a persistent burn wound infection. We confirm the causative role of a single ispA mutation in SCV emergence and increased aminoglycoside resistance. IspA is involved in the synthesis of ubiquinone, providing a possible link between electron transport and SCV formation in P. aeruginosa.

Real-time monitoring of Pseudomonas aeruginosa biofilm growth dynamics and persister cells' eradication

Biofilm formation and the appearance of persister cells with low metabolic rates are key factors affecting conventional treatment failure and antibiotic resistance. Using impedance-based measurements, crystal violet staining and traditional culture we have studied the biofilm growth dynamics of 13 Pseudomonas aeruginosa strains under the effect of seven conventional antibiotics. Real-time growth quantifications revealed that the exposure of established P. aeruginosa biofilms to certain concentrations of ciprofloxacin, ceftazidime and tobramycin induced the emergence of persister cells, that showed different morphology and pigmentation, as well increased antibiotic resistance. Whole-genome sequencing of wildtype and persister cells identified several SNPs, a genomic inversion and a genomic duplication in one of the strains. However, these mutations were not uniquely associated with persisters, suggesting that the persistent phenotype may be related to metabolic and transcriptional changes. Given that mannitol has been proposed to activate bacterial metabolism, the synergistic combination of mannitol and ciprofloxacin was evaluated on clinical 48 h P. aeruginosa biofilms. When administered at doses ≥320 mg/L, mannitol was capable of preventing persister cell formation by efficiently activating dormant bacteria and making them susceptible to the antibiotic. These results were confirmed using viable colony counting. As the tested ciprofloxacin-mannitol combination appeared to fully eradicate mature biofilms, we conclude that impedance-based biofilm diagnostics, which permits antibiotic susceptibility testing and the identification of persister cells, is of great potential for the clinical practice and could aid in establishing treatment breakpoints for emerging biofilm-related infections.

PaReBrick: PArallel REarrangements and BReaks identification toolkit

MOTIVATION

High plasticity of bacterial genomes is provided by numerous mechanisms including horizontal gene transfer and recombination via numerous flanking repeats. Genome rearrangements such as inversions, deletions, insertions and duplications may independently occur in different strains, providing parallel adaptation or phenotypic diversity. Specifically, such rearrangements might be responsible for virulence, antibiotic resistance and antigenic variation. However, identification of such events requires laborious manual inspection and verification of phyletic pattern consistency.

RESULTS

Here, we define the term 'parallel rearrangements' as events that occur independently in phylogenetically distant bacterial strains and present a formalization of the problem of parallel rearrangements calling. We implement an algorithmic solution for the identification of parallel rearrangements in bacterial populations as a tool PaReBrick. The tool takes a collection of strains represented as a sequence of oriented synteny blocks and a phylogenetic tree as input data. It identifies rearrangements, tests them for consistency with a tree, and sorts the events by their parallelism score. The tool provides diagrams of the neighbors for each block of interest, allowing the detection of horizontally transferred blocks or their extra copies and the inversions in which copied blocks are involved. We demonstrated PaReBrick's efficiency and accuracy and showed its potential to detect genome rearrangements responsible for pathogenicity and adaptation in bacterial genomes.

AVAILABILITY AND IMPLEMENTATION

PaReBrick is written in Python and is available on GitHub: https://github.com/ctlab/parallel-rearrangements.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Depletion forces drive reversible capture of live bacteria on non-adhesive surfaces

Because bacterial adhesion to surfaces is associated with infections and biofilm growth, it has been a longstanding goal to develop coatings that minimize biomolecular adsorption and eliminate bacteria adhesion. We demonstrate that, even on carefully-engineered non-bioadhesive coatings such as polyethylene glycol (PEG) layers that prevent biomolecule adsorption and cell adhesion, depletion interactions from non-adsorbing polymer in solution (such as 10 K PEG or 100 K PEO) can cause adhesion and retention of Escherichia coli cells, defeating the antifouling functionality of the coating. The cells are immobilized and remain viable on the timescale of the study, at least up to 45 minutes. When the polymer solution is replaced by buffer, cells rapidly escape from the surface, consistent with expectations for the reversibility of depletion attractions. The dissolved polymer additionally causes cells to aggregate in solution and aggregates rapidly dissociate to singlets upon tenfold dilution in buffer, also consistent with depletion. Hydrodynamic forces can substantially reduce the adhesion of aggregates on surfaces in conditions where single cells adhere via depletion. The findings reported here suggest that because bacteria thrive in polymer-rich environments both in vivo and in situ, depletion interactions may make it impossible to avoid bacterial retention on surfaces.

Extracellular DNA released by glycine-auxotrophic Staphylococcus epidermidis small colony variant facilitates catheter-related infections

Though a definitive link between small colony variants (SCVs) and implant-related staphylococcal infections has been well-established, the specific underlying mechanism remains an ill-explored field. The present study analyzes the role SCVs play in catheter infection by performing genomic and metabolic analyses, as well as analyzing biofilm formation and impacts of glycine on growth and peptidoglycan-linking rate, on a clinically typical Staphylococcus epidermidis case harboring stable SCV, normal counterpart (NC) and nonstable SCV. Our findings reveal that S. epidermidis stable SCV carries mutations involved in various metabolic processes. Metabolome analyses demonstrate that two biosynthetic pathways are apparently disturbed in SCV. One is glycine biosynthesis, which contributes to remarkable glycine shortage, and supplementation of glycine restores growth and peptidoglycan-linking rate of SCV. The other is overflow of pyruvic acid and acetyl-CoA, leading to excessive acetate. SCV demonstrates higher biofilm-forming ability due to rapid autolysis and subsequent eDNA release. Despite a remarkable decline in cell viability, SCV can facilitate in vitro biofilm formation and in vivo survival of NC when co-infected with its normal counterparts. This work illustrates an intriguing strategy utilized by a glycine-auxotrophic clinical S. epidermidis SCV isolate to facilitate biofilm-related infections, and casts a new light on the role of SCV in persistent infections.

Connection Between Chromosomal Location and Function of CtrA Phosphorelay Genes in Alphaproteobacteria

Most bacterial chromosomes are circular, with replication starting at one origin (ori) and proceeding on both replichores toward the terminus (ter). Several studies have shown that the location of genes relative to ori and ter can have profound effects on regulatory networks and physiological processes. The CtrA phosphorelay is a gene regulatory system conserved in most alphaproteobacteria. It was first discovered in Caulobacter crescentus where it controls replication and division into a stalked and a motile cell in coordination with other factors. The locations of the ctrA gene and targets of this response regulator on the chromosome affect their expression through replication-induced DNA hemi-methylation and specific positioning along a CtrA activity gradient in the dividing cell, respectively. Here we asked to what extent the location of CtrA regulatory network genes might be conserved in the alphaproteobacteria. We determined the locations of the CtrA phosphorelay and associated genes in closed genomes with unambiguously identifiable ori from members of five alphaproteobacterial orders. The location of the phosphorelay genes was the least conserved in the Rhodospirillales followed by the Sphingomonadales. In the Rhizobiales a trend toward certain chromosomal positions could be observed. Compared to the other orders, the CtrA phosphorelay genes were conserved closer to ori in the Caulobacterales. In contrast, the genes were highly conserved closer to ter in the Rhodobacterales. Our data suggest selection pressure results in differential positioning of CtrA phosphorelay and associated genes in alphaproteobacteria, particularly in the orders Rhodobacterales, Caulobacterales and Rhizobiales that is worth deeper investigation.

Genetic characteristics and phylogenetic analysis of Brazilian clinical strains of Pseudomonas aeruginosa harboring CRISPR/Cas systems

The CRISPR-Cas are adaptive immune systems found in archaea and bacteria, responsible for providing sequence-specific resistance against foreign DNA. Strains of Pseudomonas aeruginosa may carry CRISPR/Cas system types I-F, I-E and/or I-C; however, several aspects related to the epidemiology and functionality of these systems have not yet been revealed. Here, we report 13 genomes of clinical strains of P. aeruginosa from Brazil that were positive for CRISPR/Cas system types I-F and I-E, a rare feature in this species. The phylogenetic tree, which was constructed with 161 other publicly available genomes, suggested no direct relationship between positive strains, and the various types of CRISPR/Cas systems were spread throughout the tree. Comparative analysis of the genetic locations of type I-F and a specific orphan CRISPR array (without cas genes), named the LES locus, showed sequence similarities between this orphan locus and type I-F, but these LES loci were inserted in a different genomic location. We also report the presence of prophages, the presence of anti-CRISPR genes, and possibly the presence of self-targeting spacers. Here, we conclude that CRISPR/Cas is highly associated with certain lineages and is spread throughout the phylogenetic tree, showing no clear pattern of evolutionary distribution. Moreover, the LES locus might be a CRISPR1 locus related to type I-F that may have been misplaced and maintained over time. Furthermore, strains carrying I-F and I-E are rare, and not all of them are closely related. Further functional work is needed to better comprehend if aspects reported in this study are functional, including the LES locus, self-targeting spacers, anti-CRISPR protection, and I-F/I-E-carrying lineages.