Gene-Gene Interactions Dictate Ciprofloxacin Resistance in Pseudomonas aeruginosa and Facilitate Prediction of Resistance Phenotype from Genome Sequence Data

Keeping up with the pathogens: improved antimicrobial resistance detection and prediction from Pseudomonas aeruginosa genomes

BACKGROUND

Antimicrobial resistance (AMR) is an intensifying threat that requires urgent mitigation to avoid a post-antibiotic era. Pseudomonas aeruginosa represents one of the greatest AMR concerns due to increasing multi- and pan-drug resistance rates. Shotgun sequencing is gaining traction for in silico AMR profiling due to its unambiguity and transferability; however, accurate and comprehensive AMR prediction from P. aeruginosa genomes remains an unsolved problem.

METHODS

We first curated the most comprehensive database yet of known P. aeruginosa AMR variants. Next, we performed comparative genomics and microbial genome-wide association study analysis across a Global isolate Dataset (n = 1877) with paired antimicrobial phenotype and genomic data to identify novel AMR variants. Finally, the performance of our P. aeruginosa AMR database, implemented in our AMR detection and prediction tool, ARDaP, was compared with three previously published in silico AMR gene detection or phenotype prediction tools-abritAMR, AMRFinderPlus, ResFinder-across both the Global Dataset and an analysis-naïve Validation Dataset (n = 102).

RESULTS

Our AMR database comprises 3639 mobile AMR genes and 728 chromosomal variants, including 75 previously unreported chromosomal AMR variants, 10 variants associated with unusual antimicrobial susceptibility, and 281 chromosomal variants that we show are unlikely to confer AMR. Our pipeline achieved a genotype-phenotype balanced accuracy (bACC) of 85% and 81% across 10 clinically relevant antibiotics when tested against the Global and Validation Datasets, respectively, vs. just 56% and 54% with abritAMR, 58% and 54% with AMRFinderPlus, and 60% and 53% with ResFinder. ARDaP's superior performance was predominantly due to the inclusion of chromosomal AMR variants, which are generally not identified with most AMR identification tools.

CONCLUSIONS

Our ARDaP software and associated AMR variant database provides an accurate tool for predicting AMR phenotypes in P. aeruginosa, far surpassing the performance of current tools. Implementation of ARDaP for routine AMR prediction from P. aeruginosa genomes and metagenomes will improve AMR identification, addressing a critical facet in combatting this treatment-refractory pathogen. However, knowledge gaps remain in our understanding of the P. aeruginosa resistome, particularly the basis of colistin AMR.

Quorum-sensing regulation of phenazine production heightens Pseudomonas aeruginosa resistance to ciprofloxacin

Quorum sensing is a type of cell-cell communication that modulates various biological activities of bacteria. Previous studies indicate that quorum sensing contributes to the evolution of bacterial resistance to antibiotics, but the underlying mechanisms are not fully understood. In this study, we grew Pseudomonas aeruginosa in the presence of sub-lethal concentrations of ciprofloxacin, resulting in a large increase in ciprofloxacin minimal inhibitory concentration. We discovered that quorum sensing-mediated phenazine biosynthesis was significantly enhanced in the resistant isolates, where the quinolone circuit was the predominant contributor to this phenomenon. We found that production of pyocyanin changed carbon flux and showed that the effect can be partially inhibited by the addition of pyruvate to cultures. This study illustrates the role of quorum sensing-mediated phenotypic resistance and suggests a strategy for its prevention.

Binding Study of Antibacterial Drug Ciprofloxacin with Imidazolium-Based Ionic Liquids Having Different Halide Anions: A Spectroscopic and Density Functional Theory Analysis

Herein, we have shown the interaction of an antibiotic drug ciprofloxacin (CIP) with three surface-active ionic liquids (ILs), having the same cation and different anions, namely, 1-decyl-3-methylimidazoliumtetrafluoroborate [C10mim][BF4], 1-decyl-3-methylimidazolium bromide [C10mim][Br], and 1-decyl-3-methylimidazolium chloride [C10mim][Cl]. This study has been performed by exploiting various spectroscopic techniques such as steady-state fluorescence, time-resolved fluorescence, and UV-visible spectroscopy. The fluorescence emission study of CIP with ILs was performed at various concentrations of all the three ILs. The emission spectra of CIP decreased in the presence of ILs, suggesting complex formation between CIP-IL. The effect of different concentrations of ILs on the emission spectra of CIP was exploited in terms of quenching and binding parameters. Further, fluorescence emission study was validated by the time-resolved fluorescence technique by measuring the average lifetime (τavg) of CIP in the presence of all the three ILs. The τavg value of the drug changed with the addition of ILs, which suggests complex formation between the drug and ILs. This complex formation was also confirmed by UV-visible spectroscopy results of CIP with all the three ILs. Further, for evaluating the thermodynamic parameters of the CIP-IL interactions, isothermal titration calorimetry (ITC) was performed. The ITC experiment yielded the thermodynamic parameters, ΔH (change in the enthalpy of association), ΔG (Gibbs free energy change), ΔS (entropy change), and binding constant (Ka). The binding parameters driven by ITC revealed that CIP-IL interactions are spontaneous in nature and enthalpy-driven, involving hydrophobic forces. Further, the classical density functional theory (DFT) calculations were performed, which provided deep insight for CIP-IL complex formation.

The PitA protein contributes to colistin susceptibility in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic pathogen that causes a wide range of problematic infections in individuals with predisposing conditions. Infections can be treated with colistin but some isolates are resistant to this antibiotic. To better understand the genetic basis of resistance, we experimentally evolved 19 independent resistant mutants from the susceptible laboratory strain PAO1. Whole genome sequencing identified mutations in multiple genes including phoQ and pmrB that have previously been associated with resistance, pitA that encodes a phosphate transporter, and carB and eno that encode enzymes of metabolism. Individual mutations were engineered into the genome of strain PAO1. Mutations in pitA, pmrB and phoQ increased the minimum inhibitory concentration (MIC) for colistin 8-fold, making the bacteria resistant. Engineered pitA/phoQ and pitA/pmrB double mutants had higher MICs than single mutants, demonstrating additive effects on colistin susceptibility. Single carB and eno mutations did not increase the MIC suggesting that their effect is dependent on the presence of other mutations. Many of the resistant mutants had increased susceptibility to β-lactams and lower growth rates than the parental strain demonstrating that colistin resistance can impose a fitness cost. Two hundred and fourteen P. aeruginosa isolates from a range of sources were tested and 18 (7.8%) were colistin resistant. Sequence variants in genes identified by experimental evolution were present in the 18 resistant isolates and may contribute to resistance. Overall our results identify pitA mutations as novel contributors to colistin resistance and demonstrate that resistance can reduce fitness of the bacteria.

Gene-Gene Interactions Reduce Aminoglycoside Susceptibility of Pseudomonas aeruginosa through Efflux Pump-Dependent and -Independent Mechanisms

Pseudomonas aeruginosa causes a wide range of acute and chronic infections. Aminoglycosides are a cornerstone of treatment, but isolates are often resistant. The purpose of this research was to better understand the genetic basis of aminoglycoside resistance in P. aeruginosa. Bioinformatic approaches identified mutations in resistance-associated genes in the clinical isolates of P. aeruginosa. The common mutations were then engineered into the genome of P. aeruginosa reference strain PAO1. Mutations in the elongation factor gene fusA1 caused the biggest reduction in aminoglycoside susceptibility, with mutations in the two-component regulator gene amgS and the efflux pump regulator gene mexZ having less impact. This susceptibility was further reduced by combinations of mutations. Mutations in fusA1, amgS and mexZ all increased the expression of the mexXY efflux pump that is strongly associated with aminoglycoside resistance. Furthermore, the fusA1 amgS mexZ triple mutant had the highest efflux pump gene expression. Engineering fusA1 and amgS mutants lacking this efflux pump showed that fusA1 and amgS also reduce aminoglycoside susceptibility through additional mechanisms. The fusA1 and amgS mutations reduced bacterial growth, showing that these mutations have a fitness cost. Our findings demonstrate the complex interplay between mutations, efflux pump expression and other mechanisms for reducing the susceptibility of P. aeruginosa to aminoglycosides.

Ceftazidime resistance in Pseudomonas aeruginosa is multigenic and complex

Pseudomonas aeruginosa causes a wide range of severe infections. Ceftazidime, a cephalosporin, is a key antibiotic for treating infections but a significant proportion of isolates are ceftazidime-resistant. The aim of this research was to identify mutations that contribute to resistance, and to quantify the impacts of individual mutations and mutation combinations. Thirty-five mutants with reduced susceptibility to ceftazidime were evolved from two antibiotic-sensitive P. aeruginosa reference strains PAO1 and PA14. Mutations were identified by whole genome sequencing. The evolved mutants tolerated ceftazidime at concentrations between 4 and 1000 times that of the parental bacteria, with most mutants being ceftazidime resistant (minimum inhibitory concentration [MIC] ≥ 32 mg/L). Many mutants were also resistant to meropenem, a carbapenem antibiotic. Twenty-eight genes were mutated in multiple mutants, with dacB and mpl being the most frequently mutated. Mutations in six key genes were engineered into the genome of strain PAO1 individually and in combinations. A dacB mutation by itself increased the ceftazidime MIC by 16-fold although the mutant bacteria remained ceftazidime sensitive (MIC < 32 mg/L). Mutations in ampC, mexR, nalC or nalD increased the MIC by 2- to 4-fold. The MIC of a dacB mutant was increased when combined with a mutation in ampC, rendering the bacteria resistant, whereas other mutation combinations did not increase the MIC above those of single mutants. To determine the clinical relevance of mutations identified through experimental evolution, 173 ceftazidime-resistant and 166 sensitive clinical isolates were analysed for the presence of sequence variants that likely alter function of resistance-associated genes. dacB and ampC sequence variants occur most frequently in both resistant and sensitive clinical isolates. Our findings quantify the individual and combinatorial effects of mutations in different genes on ceftazidime susceptibility and demonstrate that the genetic basis of ceftazidime resistance is complex and multifactorial.

A fast lasso-based method for inferring higher-order interactions

Large-scale genotype-phenotype screens provide a wealth of data for identifying molecular alterations associated with a phenotype. Epistatic effects play an important role in such association studies. For example, siRNA perturbation screens can be used to identify combinatorial gene-silencing effects. In bacteria, epistasis has practical consequences in determining antimicrobial resistance as the genetic background of a strain plays an important role in determining resistance. Recently developed tools scale to human exome-wide screens for pairwise interactions, but none to date have included the possibility of three-way interactions. Expanding upon recent state-of-the-art methods, we make a number of improvements to the performance on large-scale data, making consideration of three-way interactions possible. We demonstrate our proposed method, Pint, on both simulated and real data sets, including antibiotic resistance testing and siRNA perturbation screens. Pint outperforms known methods in simulated data, and identifies a number of biologically plausible gene effects in both the antibiotic and siRNA models. For example, we have identified a combination of known tumour suppressor genes that is predicted (using Pint) to cause a significant increase in cell proliferation.

Deciphering the genetic network and programmed regulation of antimicrobial resistance in bacterial pathogens

Antimicrobial resistance (AMR) in bacteria is an important global health problem affecting humans, animals, and the environment. AMR is considered as one of the major components in the "global one health". Misuse/overuse of antibiotics in any one of the segments can impact the integrity of the others. In the presence of antibiotic selective pressure, bacteria tend to develop several defense mechanisms, which include structural changes of the bacterial outer membrane, enzymatic processes, gene upregulation, mutations, adaptive resistance, and biofilm formation. Several components of mobile genetic elements (MGEs) play an important role in the dissemination of AMR. Each one of these components has a specific function that lasts long, irrespective of any antibiotic pressure. Integrative and conjugative elements (ICEs), insertion sequence elements (ISs), and transposons carry the antimicrobial resistance genes (ARGs) on different genetic backbones. Successful transfer of ARGs depends on the class of plasmids, regulons, ISs proximity, and type of recombination systems. Additionally, phage-bacterial networks play a major role in the transmission of ARGs, especially in bacteria from the environment and foods of animal origin. Several other functional attributes of bacteria also get successfully modified to acquire ARGs. These include efflux pumps, toxin-antitoxin systems, regulatory small RNAs, guanosine pentaphosphate signaling, quorum sensing, two-component system, and clustered regularly interspaced short palindromic repeats (CRISPR) systems. The metabolic and virulence state of bacteria is also associated with a range of genetic and phenotypic resistance mechanisms. In spite of the availability of a considerable information on AMR, the network associations between selection pressures and several of the components mentioned above are poorly understood. Understanding how a pathogen resists and regulates the ARGs in response to antimicrobials can help in controlling the development of resistance. Here, we provide an overview of the importance of genetic network and regulation of AMR in bacterial pathogens.

Rapid fluoroquinolone resistance detection in Pseudomonas aeruginosa using mismatch amplification mutation assay-based real-time PCR

Background. Antimicrobial resistance (AMR) is an ever-increasing global health concern. One crucial facet in tackling the AMR epidemic is earlier and more accurate AMR diagnosis, particularly in the dangerous and highly multi-drug-resistant ESKAPE pathogen, Pseudomonas aeruginosa .

Objectives. We aimed to develop two SYBR Green-based mismatch amplification mutation assays (SYBR-MAMAs) targeting GyrA T83I (gyrA248) and GyrA D87N, D87Y and D87H (gyrA259). Together, these variants cause the majority of fluoroquinolone (FQ) AMR in P. aeruginosa .

Methods. Following assay validation, the gyrA248 and gyrA259 SYBR-MAMAs were tested on 84 Australian clinical P. aeruginosa isolates, 46 of which demonstrated intermediate/full ciprofloxacin resistance according to antimicrobial susceptibility testing.

Results. Our two SYBR-MAMAs correctly predicted an AMR phenotype in the majority (83%) of isolates with intermediate/full FQ resistance. All FQ-sensitive strains were predicted to have a sensitive phenotype. Whole-genome sequencing confirmed 100 % concordance with SYBR-MAMA genotypes.

Conclusions. Our GyrA SYBR-MAMAs provide a rapid and cost-effective method for same-day identification of FQ AMR in P. aeruginosa . An additional SYBR-MAMA targeting the GyrB S466Y/S466F variants would increase FQ AMR prediction to 91 %. Clinical implementation of our assays will permit more timely treatment alterations in cases where decreased FQ susceptibility is identified, leading to improved patient outcomes and antimicrobial stewardship.

Resistance to Fluoroquinolones in Pseudomonas aeruginosa from Human, Animal, Food and Environmental Origin: The Role of CrpP and Mobilizable ICEs

Fluoroquinolone resistance and the associated genetic mechanisms were assessed by antimicrobial susceptibility and whole genome sequencing in 56 Pseudomonas aeruginosa strains from human, animal, food and environmental origins. P. aeruginosa PAO1, PA7 and PA14 reference strains were also included in the study. Twenty-two strains (37%) were resistant to, at least, one fluoroquinolone agent. Correlation between the number of changes in GyrA and ParC proteins and the level of fluoroquinolone resistance was observed. Mutations or absence of genes, such as mexZ, mvaT and nalD encoding efflux pumps regulators, were also found in resistant strains. The crpP gene was detected in 43 strains (72.9%; 17 of them non-clinical strains), and coded seven different CrpP variants, including a novel one (CrpP-7). The crpP gene was located in 23 different chromosomal mobile integrative and conjugative elements (ICEs), inserted in two tRNAs integration sites. A great variety of structures was detected in the crpP-ICEs elements, e.g., the fimbriae related cup clusters, the mercury resistance mer operon, the pyocin S5 or S8 bacteriocin encoding genes, and mobilization genes. The location of crpP-like genes in mobilizable ICEs and linked to heavy metal resistance and virulence factors is of significant concern in P. aeruginosa. This work provides a genetic explanation of the fluoroquinolone resistance and crpP-associated pathogenesis of P. aeruginosa from a One-Health approach.

The First Saudi Report of Novel and Common Mutations in the gyrA and parC Genes Among Pseudomonas Spp. Clinical Isolates Recovered from Taif Area

Background and Aims

Reports examine quinolone resistance mechanisms among Pseudomonas spp. are sporadic in the Kingdom of Saudi Arabia (KSA). We previously examined the genetic bases of plasmid-mediated quinolone resistance among Pseudomonas spp. clinical isolates. This study investigated chromosomally mediated quinolone resistance mechanisms via investigation of the mutations in the gyrA and parC genes.

Methods

The minimum inhibitory concentration (MIC) to different quinolones was determined. Twenty-nine quinolone resistant Pseudomonas spp. clinical isolates were included. The gyrA and parC genes were sequenced by Sanger capillary electrophoresis. Multiple sequence alignment for the translated gyrA and parC genes was performed to identify mutation sites.

Results

Of the 29 isolates, 27 isolates were P. aeruginosa and two were P. putida. The cluster analysis of the quinolone susceptibility pattern revealed seven susceptibility phenotypes (A-G) based on susceptibility patterns rather than the MIC values. Also, 22 different susceptibility phenotypes were detected based on MIC values. All isolates exhibited a missense mutation at position 83 (S83I/T/F) of the gyrA gene in addition to six missense mutations at positions outside the QRDR of this gene. In addition, 82.8% (24/29) of the isolates harbored a missense mutation in the parC gene at position 87 (S87L), along with six novel mutations outside the QRDR of the parC gene. Haplotyping of the gyrA, parC, and the overall QRDR revealed six, 10, and 13 different haplotypes, respectively.

Conclusion

This study documents the incidence of the commonly reported mutations in the gyrA and parC genes in addition to novel mutations in these genes among Pseudomonas spp. clinical isolates recovered from KSA. Together with our previous findings, these data provide an insight into the genetic background of quinolone resistance among Pseudomonas spp. clinical isolates in KSA.

Aminoglycoside-Modifying Enzymes Are Sufficient to Make Pseudomonas aeruginosa Clinically Resistant to Key Antibiotics

Aminoglycosides are widely used to treat infections of Pseudomonas aeruginosa. Genes encoding aminoglycoside-modifying enzymes (AMEs), acquired by horizontal gene transfer, are commonly associated with aminoglycoside resistance, but their effects have not been quantified. The aim of this research was to determine the extent to which AMEs increase the antibiotic tolerance of P. aeruginosa. Bioinformatics analysis identified AME-encoding genes in 48 out of 619 clinical isolates of P. aeruginosa, with ant(2′)-Ia and aac(6′)-Ib3, which are associated with tobramcyin and gentamicin resistance, being the most common. These genes and aph(3′)-VIa (amikacin resistance) were deleted from antibiotic-resistant strains. Antibiotic minimum inhibitory concentrations (MICs) were reduced by up to 64-fold, making the mutated bacteria antibiotic-sensitive in several cases. Introduction of the same genes into four antibiotic-susceptible P. aeruginosa strains increased the MIC by up to 128-fold, making the bacteria antibiotic-resistant in all cases. The cloned genes also increased the MIC in mutants lacking the MexXY-OprM efflux pump, which is an important contributor to aminoglycoside resistance, demonstrating that AMEs and this efflux pump act independently in determining levels of aminoglycoside tolerance. Quantification of the effects of AMEs on antibiotic susceptibility demonstrates the large effect that these enzymes have on antibiotic resistance.

Express Yourself: Quantitative Real-Time PCR Assays for Rapid Chromosomal Antimicrobial Resistance Detection in Pseudomonas aeruginosa

The rise of antimicrobial-resistant (AMR) bacteria is a global health emergency. One critical facet of tackling this epidemic is more rapid AMR diagnosis in serious multidrug-resistant pathogens like Pseudomonas aeruginosa. Here, we designed and then validated two multiplex quantitative real-time PCR (qPCR) assays to simultaneously detect differential expression of the resistance-nodulation-division efflux pumps MexAB-OprM, MexCD-OprJ, MexEF-OprN, and MexXY-OprM, the AmpC β-lactamase, and the porin OprD, which are commonly associated with chromosomally encoded AMR. Next, qPCRs were tested on 15 sputa from 11 participants with P. aeruginosa respiratory infections to determine AMR profiles in vivo. We confirmed multiplex qPCR testing feasibility directly on sputa, representing a key advancement in in vivo AMR diagnosis. Notably, comparison of sputa with their derived isolates grown in Luria-Bertani broth (±2.5% NaCl) or a 5-antibiotic cocktail showed marked expression differences, illustrating the difficulty in replicating in vivo expression profiles in vitro. Cystic fibrosis sputa showed significantly reduced mexE and mexY expression compared with chronic obstructive pulmonary disease sputa, despite harboring fluoroquinolone- and aminoglycoside-resistant strains, indicating that these loci do not contribute to AMR in vivo. oprD was also significantly downregulated in cystic fibrosis sputa, even in the absence of contemporaneous carbapenem use, suggesting a common adaptive trait in chronic infections that may affect carbapenem efficacy. Sputum ampC expression was highest in participants receiving carbapenems (6.7 to 15×), some of whom were simultaneously receiving cephalosporins, the latter of which would be rendered ineffective by the upregulated ampC. Our qPCR assays provide valuable insights into the P. aeruginosa resistome, and their use on clinical specimens will permit timely treatment alterations that will improve patient outcomes and antimicrobial stewardship measures.