http://orcid.org/0000-0003-1230-5138
Figures
Abstract
Pseudomonas aeruginosa has been implicated in a wide range of post-operation wound and lung infections. A wide range of acquired resistance and virulence markers indicate surviving strategy of P. aeruginosa. Complete-genome analysis has been identified as efficient approach towards understanding the pathogenicity of this organism. This study was designed to sequence the entire genome of P. aeruginosa UY1PSABAL and UY1PSABAL2; determine drug-resistance profiles and virulence factors of the isolates; assess factors that contribute toward stability of the genomes; and thereafter determine evolutionary relationships between the strains and other isolates from similar sources. The genomes of the MDR P. aeruginosa UY1PSABAL and UY1PSABAL2 were sequenced on the Illumina Miseq platform. The raw sequenced reads were assessed for quality using FastQC v.0.11.5 and filtered for low quality reads and adapter regions using Trimmomatic v.0.36. The de novo genome assembly was made with SPAdes v.3.13 and annotated using Prokka v.2.1.1 annotation pipeline; Rapid Annotation using Subsytems Technology (RAST) server v.2.0; and PATRIC annotation tool v.3.6.2. Antimicrobial resistance genes and virulence determinants were searched through the functional annotation data generated from Prokka, RAST and PATRIC annotation pipelines; In addition to ResFinder and Comprehensive Antibiotic Resistance Database (CARD) which were employed to determine resistance genes. The PHAge Search Tool Enhanced Release (PHASTER) web server was used for the rapid identification and annotation of prophage sequences within bacterial genome. Predictive secondary metabolites were identified with AntiSMASH v.5.0. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and cas genes regions were also investigated with the CRISPRone and CRISPRFinder server. The genome sizes of 7.0 and 6.4 Mb were determined for UY1PSABAL and UY1PSABAL2 strains with G+C contents of 66.1% and 66.48% respectively. β-lactamines resistance genes blaPAO, aminoglycoside phosphorylating enzymes genes aph(3’)-IIb, fosfomycine resistance gene fosA, vancomycin vanW and tetracycline tetA were among identified resistance genes harboured in both isolates. UY1PSABAL bore additional aph(6)-Id, aph(3'')-Ib, ciprofloxacin-modifying enzyme crpP and ribosomal methylation enzyme rmtB. Both isolates were found harbouring virulence markers such as flagella and type IV pili; and also present various type III secretion systems such as exoA, exoS, exoU, exoT. Secondary metabolites such as pyochelin and pyoverdine with iron uptake activity were found within the genomes as well as quorum-sensing systems, and various fragments for prophages and insertion sequences. Only the UY1PSABAL2 contains CRISPR-Cas system. The phylogeny revealed a very close evolutionary relationship between UY1PSABAL and the similar strain isolated from Malaysia; the same trend was observed between UY1PSABAL2 and the strain from Chinese origin. Complete analyses of the entire genomes provide a wide range of information towards understanding pathogenicity of the pathogens in question.
Citation: Madaha EL, Mienie C, Gonsu HK, Bughe RN, Fonkoua MC, Mbacham WF, et al. (2020) Whole-genome sequence of multi-drug resistant Pseudomonas aeruginosa strains UY1PSABAL and UY1PSABAL2 isolated from human broncho-alveolar lavage, Yaoundé, Cameroon. PLoS ONE 15(9): e0238390. https://doi.org/10.1371/journal.pone.0238390
Editor: Yung-Fu Chang, Cornell University, UNITED STATES
Received: June 2, 2020; Accepted: August 14, 2020; Published: September 4, 2020
Copyright: © 2020 Madaha et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession HYPERLINK "VLHS00000000" and HYPERLINK "VLHR00000000" for UY1PSABAL and UY1PSABAL2 respectively. The version described in this paper is the first version. The raw sequencing data sets have also been deposited in the NCBI Sequence Read Archive (SRA) with accession numbers HYPERLINK "SRP214447" and HYPERLINK "SRP214569" respectively.
Funding: This work was supported by the research grants from the Food Safety and Security Niche Area, Faculty of Natural and Agricultural Sciences, North West University, South Africa and the Incentive Funding for Rated Researchers by the National Research Foundation, South Africa provided to C.N.A. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Pseudomonas aeruginosa is an oxidase positive Gram-negative rod non-fermentative aerobic organism, mobile with a monotrichous polar flagellum [1]. It is an opportunistic pathogen isolated mostly in health care associated infections [2]. Like other Gram-negative bacteria, P. aeruginosa can express resistance through numerous mechanisms such as production of antibiotic modifying or degrading enzymes, active efflux pump and target modification [3]. Some strains show resistance to all clinically important antibiotics families. Strains resistant to at least three classes of anti-pseudomonas antibiotics are recognised as a multi-drug resistant Pseudomonas aeruginosa (MDR-PSA) [4]. The spread of MDR-PSA is increasing both within communities and hospital environments leading to severe infection cases and eventual therapeutic impasse [5]. Infections due to multidrug resistant pathogens impose overwhelming burden on the public health through increased hospital stay, cost of treatment, direct economic loss and death of patients from previously treatable disease conditions [6–9].
Persistence of MDR-PSA in the host environment is mainly due to various virulence factors expressed and reinforced by their multiple resistance mechanisms. Understanding the mechanisms involved in bacterial-host interaction can help to find new targets to thwart their spread. Shlaes [5] in his work: the perfect storm, mentioned that if the world strive to understand how resistance spread among microorganisms, then it could stop it and if the mechanisms by which the bacteria become resistant could be well scrutinized, we could find ways around them with discovery of novel antibiotics. Various tools have been developed to this effect and brought relevant knowledge concerning bacterial persistence and resistance.
Whole-genome sequencing (WGS) is one of the latest approaches to study resistance mechanisms in organisms. It gives massive information on the genes present within a pathogen, since this technique is capable of processing multiple DNA sequences in parallel with high throughput at reduced cost and time [10, 11]. It is a powerful protocol, which has revolutionized the study of microbial genomes. Next generation sequence has been used to investigate numerous genetic problems. To the best of our knowledge, there is no published data on P. aeruginosa whole-genome sequencing from clinical samples from Cameroon. Proper understanding of the evolution and physiology of these clinically relevant organisms will be archived by obtaining and comparing their complete genome. This study therefore evaluates the entire genomes of the two selected strains of MDR-PSA isolated from broncho-alveolar-lavage of human patients with lung infections in Yaoundé, Cameroon with a view to obtaining insight into the genes which are involved in the wide range resistance and virulence of the strains.
Methodology
Ethics statement
Ethical clearance was sought from the Central Region Ethical committee of the Central region, Cameroon, and an approval number (0259/CRESHC/2019) was assigned to the study. Administrative authorizations were obtained from the hospitals where bacteria isolates were obtained.
Identification and selection of the bacteria strains
The bacterial strains were isolated from broncho-alveolar-lavage of human patients with lung infections and identified through PCR amplification of Pseudomonas species-specific 16S rRNA [F: GACGGGTGAGTAATGCCTA; R: CACTGGTGTTCCTTCCTATA] and P. aeruginosa specific 16S rRNA [F: GGGGGATCTTCGGACCTCA; R: TCCTTAGAGTGCCCACCCG] regions [2]. Both P. aeruginosa UY1PSABAL and P. aeruginosa UY1PSABAL2 were selected based on their multi-drug resistance profiling and ability to form strong biofilms on polystyrene at 25 and 37 °C from our previous study [12]. The susceptibility test was carried out using disk-diffusion method following the CASFM / EUCAST: Société Française de Microbiologie [13]. The test antibiotic drugs include: piperacillin (30μg); piperacillin-tazobactam (30/6μg); ticarcillin (75μg); ticarcillin-clavulanic acid (75/10μg); cefepime (30μg); ceftazidime (10μg); cefoperaxone (30μg); imipenem (10μg); meropenem (10μg); aztreonam (30μg); ciprofloxacin (5μg); levofloxacin (5μg); amikacin (30μg); gentamicin (10μg); netilmicin (10μg) and tobramycin (10μg).
Complete genome sequence
The sequenced genomes of the MDR P. aeruginosa UY1PSABAL and UY1PSABAL2 were obtained through whole genome sequencing using an Illumina Miseq platform. One ng of the genomic DNA was tagmented with the Nextera XT DNA library prep kit according to the manufacturer’s protocol. The kit reagents fragment the DNA with simultaneous addition of adapter sequences. The libraries were amplified with a limited-cycle PCR program (12 cycles) to add the index 1 (i7) and index 2 (i5) adapters, containing sequences required for cluster generation of the Illumina flow cell. The library was purified using 0.6x Agencourt AMPure XP beads (Beckman Coulter). The quality and sizes of the resulting DNA fragments were evaluated on a 1.5% agarose gel. The libraries were quantified with a fluorometric method (Qubit, Life Technologies) and normalised to4 nM using a standard dilution method. The libraries were pooled, denatured with 0.1 N NaOH and diluted to the final loading concentration of 12 pmol. An identically treated PhiX control was added to a final concentration of 1%. Paired-end sequencing was done on an Illumina MiSeq platform using a MiSeq Reagent Kit V3 600 cycles.
Quality control, trimming, assembling and annotation
The raw sequenced reads were assessed for quality using FastQC v.0.11.5 [14] and filtered for low quality reads and adapter regions using Trimmomatic v.0.36 [15]. The de novo genome assembly was made with SPAdes v.3.13 [16] and annotated using Prokka v.2.1.1 annotation pipeline [17]; Rapid Annotation using Subsytems Technology (RAST) server v.2.0 [18]; and PATRIC annotation tool v.3.6.2 [19].
Identification of the resistome, virulome and Mobile Genetic Elements (MGEs)
Several clinically important antimicrobial resistance genes and virulence determinants were searched through the functional annotation data generated from Prokka, RAST and PATRIC annotation pipelines. In addition to ResFinder [20] and Comprehensive Antibiotic Resistance Database (CARD) [21] that were employed to determine the presence of resistance genes, the PHAge Search Tool Enhanced Release (PHASTER) web server [22] was used for the rapid identification and annotation of prophage sequences within bacterial genomes. The PHAge Search Tool Enhanced Release (PHASTER) web server [22] was used for the rapid identification and annotation of prophage sequences within bacterial genome. Predictive secondary metabolites were identified with AntiSMASH v.5.0 [23]. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated genes (Cas) were also investigated with the CRISPRone [24] and CRISPRFinder server [25].
Phylogeny.
Evolutionary relationship between the strains UY1PSABAL and UY1PSABAL2; and twenty four other strains of randomly selected P. aeruginosa isolated from Homo sapiens between 2014 and 2016 across Africa, America, Asia and Europe was determined using PATRIC Fastree platform [19]. The Newick result file was then assessed using the Figtree v1.4.3 software to view the tree and annotation. The generated tree was labelled according to the codified name, which are in consistency with the strain’s country of collection and the GenBank accession number.
Results
Basic statistical features of the genome assembly
The genome assembly size of the P. aeruginosa UY1PSABAL is 7,029,327 bp distributed into 136 contigs and having a G+C content of 66.1%. The number of protein-coding sequences within the genome is about 6,907, with 63 non-coding RNAs (including tRNAs and rRNAs). The N50 of the genome is 118,711 bp while the L50 is 18. On the other hand, the genome assembly size of the P. aeruginosa UY1PSABAL2 is 6,354,435 bp along 96 contigs with 66.48% G+C contents. The genome consists of about 6,046 protein-coding sequences and 63 non-coding RNAs. The genome’s N50 is 215,321 bp and the L50 is 10. The circular genome representing the entire genome assembly is shown in Fig 1.
The circular genome was generated with PATRIC sever v3.6.2.
Antimicrobial resistance profiling
Phenotypic characterization of the antimicrobial resistance patterns of the P. aeruginosa UY1PSABAL and UY1PSABAL2 revealed that both strains have acquired resistance against a large proportion of clinically important antibiotics. The UY1PSABAL showed resistance against 60% of the test drugs while the UY1PSABAL2 was resistant against 41% but both strains displayed similar resistance patterns against 12% of the drugs tested (Table 1).
On the other hand, a number of resistance genes harboured within the genomes of the UY1PSABAL and the UY1PSABAL2 strains were identified through functional annotation generated from the ResFinder and CARD platforms. Nine predicted regions encoding for levofloxacin resistance genes, in addition to other coding sequences that confer resistance to quinolones such as crpP, gyrA (variant/mutant), parE (variant/mutant) was detected only within the genome of UY1PSABAL while basS was detected only in the UY1PSABAL2. Both strains harboured resistance genes to β-lactams, aminoglycosides, tetracyclines, fosfomycins and glycopeptides (Table 2).
Identification of virulent determinants
The protein-encoding sequences putative for important virulent factors such as flagella protein biosynthesis, adherence motility, endotoxin, Type IV pili Adherence Twitching motility, ion uptake, anti-phagocytosis (Serum resistance), and Type I, Type II, Type III and Type VI secretion systems among others were found common to both isolates (S1 Table).
Identification of the Mobile Genetic Elements (MGE)
Transposable elements.
A plethora of CDS putative for different transposable elements were found in both isolates (S2 Table). A total number of 29 coding sequences that are known to encode for transposable elements were found associated with P. aeruginosa UY1PSABAL, among which is a single phage transposase located between 20617–18833 nucleotide position along the reverse DNA strand. Transposase proteins InsN and InsO belonging to the same insertion sequence family IS911 were found common to both isolates at various positions. The insertion element IS407 coding for transposase proteins commonly associated with Burkholderia multivorans was only found within the genome of UY1PSABAL2.
Prophage and CRISPR arrays.
Prophage search by PHASTER showed the presence of coding sequences (CDS) which are predictive for five intact, three questionable and two incomplete prophages in the UY1PSABAL. On the other hand, only one putative CDS was found each for intact, questionable and incomplete prophages within UY1PSABAL2. A search through functional annotation for prophage related genes revealed 263 phage related genes in UYPSABAL and 83 identified in UYPSABAL2. Analyses of both strains for CRISPR-Cas system revealed three CRISPR arrays with associate cas genes, only within UY1PSABAL2 (Table 3). The CRISPRone software, together with the CRISPR-Cas system, revealed six cas-protein encoding sequences: cas1, cas3, cas8f, cas5f, cas7f and cas6f within UY1PSABAL2 genome.
Secondary metabolites
Determination of secondary metabolites within the two isolates revealed nine different secondary metabolites distributed on 17 gene clusters and 15 gene clusters within the genomes of UY1PSABAL and UY1PSABAL2 respectively. Many of the identified secondary metabolites are shared by both isolates. The secondary metabolites most frequently found in the isolates were peptide synthetase cluster (NRPS)–like betalactone cluster, phenazine, homoserine lactone, bacteriocin, non-ribosomal NRPS, NRPS–like, N-acetylglutaminylglutamine amide (NAGGN), thiopeptide and phenazine. The two NRPS clusters found in both isolates showed 100% similarity, to 2-amino-4-methoxy-trans-3-butenoic acid and pyochelin synthesis clusters. Cyclodipeptide synthetase (CDPS) is the secondary metabolite found only in UY1PSABAL2.
Phylogeny
The phylogenetic tree revealing evolutionary relationship between UY1PSABAL and UY1PSABAL2 as compared with other selected strains of P. aeruginosa isolated from similar sources across continents is shown in Fig 2. A comparative analysis of the phylogeny showed that UY1PSABAL (Cameroon_VLHS00000000) and Malaysia_MPCQ00000000 belong to the same clade and both have about 0.9% divergence compared to the common ancestor. UY1PSABAL2 (Cameroon_VLHR00000000) strain and China_VWQQ01000000 belonged to the same monophyletic group and hence may have originated from a common ancestor. They are closely derived from the common ancestor with less divergence percentage (about 0.3%). Comparing with other strains sampled for this phylogenetic tree building, UY1PSABAL1 and UY1PSABAL2 strains are closely related to Asian strains.
The tree was generated using PATRIC v3.6.2.
Discussion
Genome-based characterisation of isolates has provided solutions to several problems in species identification and pathogenicity. The entire genomes of the two MDR-PSA isolates used in this study were sequenced with a view to accessing their genetic make-up.
Motility and bacteria adherence
Flagella are the most common apparatus used by bacteria for motility. Numerous encoding sequences putative for flagella formation were observed within the current isolates. This could server to improve their motility and consequently pathogenicity. Bacteria equipped with flagella are found to be more resistant to surfactant protein A (SP-A) [26] which is an important lung innate immune protein that kills microbial pathogens through opsonization and membrane permeability [27]. Both UY1PSABAL and UY1PSABAL2 were found equipped with flagella and this could have contributed to their resistance to opsonization by cells in the lung where the strains were isolated. In addition, specific pili such as type IV pili have been reportedly playing significant role in promoting motility and antimicrobial resistance in P. aeruginosa [27]. Piliated P. aeruginosa have been implicated in severe cases of pneumonia, bacteraemia, and increased mortality as compared to non-piliated strains [28]. Presence of this cellular component in both MDR-PSA strains under study could be associated with severity of lung infections [29].
Cellular secretion systems
A number of secretion systems were identified common to both isolates. For instance, type II secretion systems that secrete folded proteins, such as pseudolysin (lasB); phospholipase C (PlcH); or lipase (LipA), from the periplasm into the extracellular milieu [30]. Type III secretion systems with robust virulence factor that allows translocation of specific bacterial effector proteins into the host cell cytosol or the plasma membrane [31, 32]. The presence of a functional T3SS was associated with bacterial persistence in the lungs, higher relapse rates, and increased mortality in patients with acute respiratory infections caused by P. aeruginosa [33]. Type VI (T6SS) is also an identified secretion systems that plays crucial role in competition and pathogenesis [34]. Both isolates contained T1SS, T2SS, T3SS and T6SS which are typical of P. aeruginosa secretion systems. These are channels by which the bacteria uptake or extrude molecules or toxins. Thus paved ways for the UY1PSABAL and UY1PSABAL2 to infect and persist in the host cells.
Alginate production/regulation and anti-phagocytosis
Alginate serves as an extracellular matrix component involved in biofilm formation. It confers resistance to antibiotics and prevents phagocytosis by the immune cells such as macrophages of the host [35]. Synthesis of alginate is a complex process, involving twelve genes which are belong to the algA-algD operon [36]. Alginate production results in P. aeruginosa conversion from the non-mucoid to mucoid phenotype. Mucoid P. aeruginosa strains have been associated with a significant clinical deterioration. Beside protein-coding genes for alginate production, others anti-phagocytosis genes were equally found within the genomes of UY1PSABAL and UY1PSABAL2. These features represent important virulent traits in both isolates. There are novel therapeutic treatments targeting alginate synthesis [37], which can be a possibility for treatment in infections by UY1PSABAL and UY1PSABAL2.
Quorum sensing
Bacteria can display different behaviours depending on cell density, and quorum sensing systems (QS) play a key role in sensing population density [38]. Homoserine lactone is the most studied quorum sensing auto-inducer molecule involved in the las and rhl QS systems. The rhlI, lasI homoserin lactone synthase, pqsE Pseudomonas quinolone system gene belonging to rhl, las and pqs QS system genes and lasR gene encoding for a transcriptional regulator were found in the UY1PSABAL and UY1PSABAL2 genomes. The QS plays a major role in virulence of P. aeruginosa [37]. It allows P. aeruginosa, for example, in the presence of high bacterial numbers to produce two P. aeruginosa toxins found within these isolates: an elastase (protease) and a pyocyanin (syderophore) [39], which will subsequently harm their host.
Hydrolytic enzymes and ion uptake system
Proteases are enzymes with protein or peptide hydrolysis activity. Zinc metalloprotease such as LasA and LasB have elastolytic activity on human tissue, most especially in the lungs. Other hydrolytic enzymes, phospholipase C encoded by plcH and plcN genes can also be active on lung cells’ biomolecules such as Phospholipids contained in surfactants [40]. Among the phospholipases produced by P. aeruginosa, the thermo-stable haemolysin PlcH has the particularity of attenuating the reactive oxygen species (ROS) produced by neutrophils [31]. The genome of UY1PSABAL and UY1PSABAL2 contain the above-mentioned phospholipase C genes as well as zinc metalloprotease virulence factor genes. This further equipped the strains with capability to hydrolyse host cell proteins and avoid ROS produced by host immune cells, hence, could responsible for the bacterial persistence in the host environment.
Moreover, biological metal ions, including Co, Cu, Fe, Mg, Mn, Mo, Ni and Zn ions, are necessary for the survival and the growth of all microorganisms. They can serve as a cofactor in various biological mechanisms and are essential for microbial pathogenicity [41, 42]. Iron (Fe3+) is an essential element for living organisms. It is involved in numerous cellular processes of various organisms [43]. The two siderophores: pyochelin (Pch) and pyoverdine (Pvd) observed within UY1PSABAL and UY1PSABAL2 help in iron assimilation. Pyochelin (Pch) and pyoverdine (Pvd) chelate iron in the extracellular medium; the iron is then transported into the cell via a specific outer membrane transporter; FptA and FpvA respectively [42].
Antibiotic resistance
The P. aeruginosa UY1PSABAL and P. aeruginosa UY1PSABAL2 were characterized as multidrug-resistant (MDR) after the phenotypic profiling revealed that both isolates are resistant to a wide range of clinically important antibiotics. This was further confirmed following the genomic analysis of the resistance genes (Table 2). Several published studies on P. aeruginosa genome characterization have highlighted the prevalence of blaPAO and blaOXA50 which are classified as β-lactamases resistance genes [44, 45]; and this does not make a difference in the P. aeruginosa strains under study. In the study conducted by Girlich and colleagues [32], it was confirmed that insertion of blaOXA-50 into a P. aeruginosa host, decreases susceptibility to ampicillin, ticarcillin, moxalactam and meropenem. Thus, presence of blaOXA-50 within the UY1PSABAL2 genome could explain its resistance to meropenem. The genes blaOXA-395 and blaOXA-486 which are present respectively within UY1PASBAL and UY1PSABAL2, are commonly acquired resistance genes found within P. aeruginosa genome [45, 46].
Furthermore, quinolones resistance genes were also found common to both isolates. UY1PSABAL genome sequence showed nine predicted regions for levofloxacin and plasmid encoded ciprofloxacin resistance gene crpP which has been demonstrated to induce resistance against ciprofloxacin [47]. A gyrA variant/mutant and parE variant/mutant were observed within the same strain. The UY1PSABAL bearing a wide range of quinolone resistance markers explains its resistance to whole anti-pseudomonas quinolone (ciprofloxacin and levofloxacin) tested as compared to UY1PSABAL2, which did not carry any of these resistance genes or resistant variants, and eventually showed no resistance to quinolone antibiotics.
Additionally, rmtB is a 16S rRNA methylase which confers resistance to almost all aminoglycosides except streptomycin [48]. This ability was verified in UY1PSABAL and was confirmed harbouring the rmtB gene and consequently displayed resistance to all tested aminoglycosides (amikacin, gentamicin, gentamicin and tobramycin) represented by AMKR GENR NETR TOBR. Apart from the rmtB gene, aminoglycosides modifying enzymes (AMEs) were found within UY1PASBAL and UY1PSABAL2 genomes, and it was exclusively aminoglycoses phosphotransferases (APHs) that mostly confer resistance to amikacin and isepamicin [49]. The aph(3’)-IIb variant found in both UY1PASBAL and UY1PSABAL2 genome has been observed among MDR P. aeruginosa from several studies [44, 45, 50].
Identification of Class B and C tetracyclin; fosfomycin fosA; vancomycin B-type (vanW); and phenicols catB7 resistance genes within the UY1PASBAL and UY1PSABAL2 genome suggest that these strains can express resistance to above-mentioned antibiotic families. The presence of vancomycin resistance gene within these genomes is worrisome because of the possibility of interspecies transfer [51] of this gene to Staphylococcus aureus for which vancomycin is the front line treatment antibiotic [52]. Efflux pump systems belong to five different classes. In agreement with several scientific reports in the literature [3, 53], Resistance-Nodulation-Division family (RND) was the most represented also within our isolate. And the efflux pump gene determinants were, as well, found in high numbers within both strains. It is important to note that efflux pump resistance determinants were found predominant within the UY1PSABAL2 strain.
Mobile genetic elements and CRISPR-Cas system
Several coding sequences putative for mobile genetic elements are located within the genome of both isolates; these include, transposable elements and prophages. Transposable elements such as insertion sequences (IS) play a major role in the prokaryote genomes plasticity and are generally small mobile elements that typically carry one and sometimes two transposase (tnp) genes [54]. Transposase InsO, for IS911 and InsN transposase for IS407 (Burkholderia multivorans) were found within UY1PSABAL and UY1PSABAL2. Insertion sequence IS911, which is often located upstream of several β-lactamases [55], were also identified in both UY1PSABAL and UY1PSABAL2 strains. The IS407 has been demonstrated to promote the expression of silent lac operon in Transposons in Burkholderia cepacia. Transposons are DNA segments that can move from one genetic location to another. Transposon 4652 (Tn4652) identified within the genome of UY1PSABAL is a derivative of the toluene degradation transposon Tn4651 that belongs to the Tn3 which is a non-composite transposon bearing the β-lactamase resistance gene family of transposons [56]. The transposable elements found within the current isolates could have introduced the resistance genes that serve as the assets to resist antimicrobial agents, most especially β-lactams.
Clustered regularly interspaced short palindromic repeats (CRISPR) with associated cas genes is an adaptive immunity against mobile genetic elements, especially phages [57]. The combined results of PHASTER and RAST annotation showed a high number of intact prophage related genes in UY1PSABAL, indicating that various bacteriophage residues have infected this strain. The infection of UY1PSABAL was probably been made easier as a result of its lack of the “immune system” against bacteriophages and insertion sequences due to the absence of CRISPR-cas system. On the contrary, UY1PSABAL2 contained fewer prophages and IS related coding sequences; the presence of CRISPR-cas system might have protected the strain from invasion by bacteriophages.
Secondary metabolites
Secondary metabolites are naturally produced compounds that serve as defensive molecules for the organisms producing them [58]. Various gene clusters coding for non-ribosomal peptide synthetases (NRPSs), non-proteinogenic peptides which are often toxic for prokaryotes and eukaryotes [59] were identified within the two isolates. Two NRPS clusters appeared to match with 100% similarity to clusters encoding 2-amino-4-methoxy-trans-3-butenoic acid and pyochelin within both UY1PSABAL and UY1PSABAL2 genomes. Some of the NRSPs were similar to pyoverdine cluster with less than 100% similarity. Pyoverdine is one of the common siderophores found within P. aeruginosa species and could represent a novel drug or vaccine target [37]. The P. aeruginosa toxin l,2-amino-4-methoxy-trans-3-butenoic acid (AMB) is a non-proteinogenic amino acid, which is toxic to prokaryotes and eukaryotes [60].
N-acetylglutaminyl-glutamine amide gene (NAGGN) cluster which plays a key role in cell osmoprotection of the bacterial cells were also observed. It is the most abundant solute at high salt concentration within P. aeruginosa [61]. β-lactone is another class of natural product found within these isolates. They have a unique reactivity among natural products and very suitable for the inhibition of hydrolases, transferases, ligases, and oxidoreductases unlike β-lactams which primarily target penicillin-binding proteins [62]. Phenazine is another secondary metabolite produced by P. aeruginosa and equally found within the genome of UY1PSABAL and UY1PSABAL2. Phenazine acts by altering the expression of immunomodulatory proteins by human airway epithelial cells [63].
Additionally, bacteriocin, the heat-stable ribosomally synthesized antibacterial peptides that are active by either killing or inhibiting the growth of other neighbouring bacteria [64] is also found within the genomes. This could serve as signalling peptide in the quorum sensing system [65]. Thiopeptide cluster is a biosynthetic gene cluster for a class of antibiotics widely distributed in genomes and metagenomes of the human microbial ecosystem [66]. Thiopeptides are known to be inactive against Gram-negative bacteria, nevertheless, thiostrepton; a thiopeptide antibiotic was found to stimulate P. aeruginosa biofilm formation [67]. This might have contributed to the production of biofilm exhibited by both isolates.
Phylogeny
Analysis of the phylogenetic relationship between the UY1PSABAL and UY1PSABAL2 with other strains revealed that the strains under study were closely related to Asian (Malaysian and Chinese) strains. This could be attributed to the consequence of population dynamic. In fact, the last two decades were characterized with increased economic exchange between Asia and Africa; leading to increase in the cross-border travel between the two continents [68]. When populations move, they carry with them pathogens and cross-exchange of genetic material could occur. The phylogenetic proximity of our isolates to Malaysian and Chinese strains can be explained by the improved economic migration and social interaction between the two regions.
Conclusion
This work focused on the characterisation of two MDR P. aeruginosa strains isolated from broncho-alveolar-lavage of lung-infected patients. The WGS approach was used to obtain a deeper insight into the complete DNA sequences of UY1PSABAL and UY1PSABAL2 strains. Numerous resistance and virulent genes observed in the genome of these isolates indicate their potential public implications, the severity of infections they could cause in humans and the difficulties in the management of these complications. The diverse MGEs identified in both strains demonstrate how dynamic each strain was during its evolutionary history. UY1PSABAL and UY1PSABAL2 strains are well equipped with secondary metabolites which serve as important asset for them to invade and persist in the host environment. Further studies are required to have a wider understanding of genetic features of P. aeruginosa isolates circulating within communities in Cameroon. This could ultimately bring about informed government policy and adequate community awareness toward curbing the transmission and continuous emergence of such dangerous pathogens.
Supporting information
S1 Table. Important virulence determinants located withinUY1PSABAL and UY1PSABAL2 genomes.
https://doi.org/10.1371/journal.pone.0238390.s001
(DOCX)
S2 Table. Locations and sizes of identified transposable elements in P. aeruginosa UY1PSABAL and UY1PSABAL2.
https://doi.org/10.1371/journal.pone.0238390.s002
(DOCX)
References
- 1.
Wisplinghoff H, Pseudomonas spp., Acinetobacter spp. and Miscellaneous Gram-Negative Bacilli in: Cohen, J., Powderly, W., Opal, S. (Eds.), Infecious deseases, 2017; pp. 1579–1599.
- 2. Spilker T, Coenye T, Vandamme P, LiPuma JJ. PCR-based assay for differentiation of Pseudomonas aeruginosa from other Pseudomonas species recovered from cystic fibrosis patients, J Clin Microbiol, 2004; 42: 2074–2079. pmid:15131172
- 3. Breidenstein EB, de la Fuente-Núñez C, Hancock RE. Pseudomonas aeruginosa: all roads lead to resistance, Trends microbiol, 2011; 19: 419–426. pmid:21664819
- 4. Magiorakos AP, Srinivasan A, Carey R, Carmeli Y, Falagas M, Giske C, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance, Clin Microbiol infec, 2012; 18: 268–281.
- 5.
Shlaes DM. Antibiotics: the perfect storm. Springer Science & Business Media, CT, USA, 2010.
- 6. De Angelis G, Fiori B, Menchinelli G, D’Inzeo T, Liotti FM, Morandotti GA, et al. Incidence and antimicrobial resistance trends in bloodstream infections caused by ESKAPE and Escherichia coli at a large teaching hospital in Rome, a 9-year analysis (2007–2015), Eur J Clin Microbiol Infec Dis, 2018; 37: 1627–1636.
- 7. Bassetti M, Vena A, Croxatto A, Righi E, Guery BJD. How to manage Pseudomonas aeruginosa infections, Drugs in Context, 2018b; 7: 212527
- 8. Bassetti M, Peghin M, Vena A, Giacobbe DRJFiM. Treatment of infections due to MDR Gram-negative bacteria, Front Med, 2019; 6: 74.
- 9. Alayande KA, Aiyegoro OA, Ateba CN. Probiotics in Animal Husbandry: Applicability and Associated Risk Factors, Sustainability, 2020; 12: 1087.
- 10. Zhou X, Ren L, Li Y, Zhang M, Yu Y, Yu J. The next-generation sequencing technology: a technology review and future perspective, Sci China Life Sci, 2010; 53: 44–57. pmid:20596955
- 11. Schwarze K, Buchanan J, Taylor JC, Wordsworth S. Are whole-exome and whole-genome sequencing approaches cost-effective? A systematic review of the literature, Genet Med, 2018; 20: 1122–1130. pmid:29446766
- 12. Madaha EL, Gonsu HK, Bughe RN, Fonkoua MC, Ateba CN, Mbacham WF. Occurrence of blaTEM and blaCTXMgenes and biofilm- forming ability among clinical isolates of Pseudomonas aeruginosa and Acinetobacter baumannii in Yaoundé, Cameroon. Microorganisms 2020; 8 (708): 1–22.
- 13.
CASFM / EUCAST: Société Française de Microbiologie, recommandations, 2016
- 14.
Andrews S. FastQC: a quality control tool for high throughput sequence data, 2010.
- 15. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data, Bioinformatics, 2014; 30: 2114–2120. pmid:24695404
- 16. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing, J. Comput. Biol, 2012; 19: 455–477. pmid:22506599
- 17. Seemann T. Prokka: rapid prokaryotic genome annotation, Bioinformatics, 2014; 30: 2068–2069. pmid:24642063
- 18. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST), Nucleic Acids Res, 2013; 42: D206–D214. pmid:24293654
- 19. Wattam AR, Davis JJ, Assaf R, Boisvert S, Brettin T, Bun C, et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center, Nucleic Acids Res, 2016; 45: D535–D542. pmid:27899627
- 20. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes, J Antimicrob Chemother, 2012; 67: 2640–2644. pmid:22782487
- 21. Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database, Nucleic Acids Res., gkw1004, 2016.
- 22. Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool, Nucleic Acids Res, 2016; 44: W16–W21. pmid:27141966
- 23. Blin K, Shaw S, Steinke K, Villebro R, Ziemert N, Lee SY, et al. antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline, Nucleic Acids Res., 2019.
- 24. Zhang Q, Ye Y. Not all predicted CRISPR–Cas systems are equal: isolated cas genes and classes of CRISPR like elements, BMC bioinform, 2017; 18: 92.
- 25. Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats, Nucleic Acids Res, 2007; 35: W52–W57. pmid:17537822
- 26. Zhang S, McCormack FX, Levesque RC, O’Toole GA, Lau GW. The flagellum of Pseudomonas aeruginosa is required for resistance to clearance by surfactant protein A, PloS one, 2007; 2: 1–11.
- 27. Tan RM, Kuang Z, Hao Y, Lau GW. Type IV Pilus of Pseudomonasaeruginosa Confers Resistance to Antimicrobial Activities of the Pulmonary Surfactant Protein-A, J innate immun, 2014; 6: 227–239. pmid:24080545
- 28. Tang H, Kays M, Prince A. Role of Pseudomonas aeruginosa pili in acute pulmonary infection, Infect Immun, 1995; 63: 1278–1285. pmid:7890385
- 29. Tan RM, Kuang Z, Hao Y, Lee F, Lee T, Lee RJ, et al. Type IV pilus glycosylation mediates resistance of Pseudomonas aeruginosa to opsonic activities of the pulmonary surfactant protein A. Infect Immun. 2014; 83: 1339–1346.
- 30. Costa TR, Felisberto-Rodrigues C, Meir A, Prevost MS, Redzej A, Trokter M, et al. Secretion systems in Gram-negative bacteria: structural and mechanistic insights, Nat Rev Microbiol, 2015; 13: 343–359. pmid:25978706
- 31. de Bentzmann S, Plésiat P. Pseudomonas aeruginosa: une virulence complexe, Revue Francophone des Laboratoires 2011; 73–81.
- 32. Girlich D, Naas T, Nordmann P. Biochemical characterization of the naturally occurring oxacillinase OXA-50 of Pseudomonas aeruginosa, Antimicrob Agents chemother, 2004; 48: 2043–2048. pmid:15155197
- 33.
Anantharajah A. Contribution to the understanding of Pseudomonas aeruginosa virulence in search of innovative therapeutic approaches: Focus on the type III secretion system, Secteur des sciences de la santé. Université catholique de Louvain, Belgium, p. 231, 2017.
- 34. Wood TE, Howard SA, Förster A, Nolan LM, Manoli E, Bullen NP, et al. The Pseudomonas aeruginosa T6SS Delivers a Periplasmic Toxin that Disrupts Bacterial Cell Morphology, Cell Rep, 2019; 29: 187–201. pmid:31577948
- 35. Rehman ZU, Rehm BH. Dual roles of Pseudomonas aeruginosa AlgE in secretion of the virulence factor alginate and formation of the secretion complex, Appl Environ Microbiol, 2013; 79: 2002–2011. pmid:23335756
- 36. Pulcrano G, Iula DV, Raia V, Rossano F, Catania MR. Different mutations in mucA gene of Pseudomonas aeruginosa mucoid strains in cystic fibrosis patients and their effect on algU gene expression, New Microbiol, 2012; 35: 295–305. pmid:22842599
- 37. Fothergill JL, Winstanley C, James CE. Novel therapeutic strategies to counter Pseudomonas aeruginosa infections, Expert Rev anti-infect ther, 2012; 10: 219–235. pmid:22339195
- 38. Veesenmeyer JL, Hauser AR, Lisboa T, Rello J. Pseudomonas aeruginosa virulence and therapy: evolving translational strategies, CritCare Med, 2009; 37: 1777.
- 39. Hauser AR. Pseudomonas aeruginosa: so many virulence factors, so little time, Crit Care Med, 2011; 39: 1–3.
- 40. Fazeli N, Momtaz H. Virulence gene profiles of multidrug-resistant Pseudomonas aeruginosa isolated from Iranian hospital infections, Iran Red Crescent Med J, 2014; 16: 1–10.
- 41. Braud A, Hannauer M, Mislin GL, Schalk IJ. The Pseudomonas aeruginosa pyochelin-iron uptake pathway and its metal specificity, J bacteriol, 2009; 191, 3517–3525. pmid:19329644
- 42. Schalk IJ, Cunrath O. An overview of the biological metal uptake pathways in Pseudomonas aeruginosa, Environ microbiol, 2016; 18: 3227–3246. pmid:27632589
- 43. Hofmann M, Heine T, Schulz V, Hofmann S, Tischler D. Draft genomes and initial characterization of siderophore producing pseudomonads isolated from mine dump and mine drainage, Biotechnol Rep, 2020; 25: 1–8.
- 44. Grandjean T, Le Guern R, Duployez C, Faure K, Kipnis E, Dessein R., Draft Genome Sequences of Two Pseudomonas aeruginosa Multidrug-Resistant Clinical Isolates, PAL0. 1 and PAL1. 1, Microbiol Resour Announc, 2018; 7: 1–2.
- 45. Hussain M, Suliman M, Ahmed A, Altayb H, Elneima E. Draft genome sequence of a multidrug-resistant Pseudomonas aeruginosa strain isolated from a patient with a urinary tract infection in Khartoum, Sudan, Genome Announc, 2017; 5: 1–2.
- 46. Taiaroa G, Samuelsen Ø, Kristensen T, Økstad OAL, Heikal A. Complete Genome Sequence of Pseudomonas aeruginosa K34-7, a Carbapenem-Resistant Isolate of the High-Risk Sequence Type 233, Microbiol Resour Announc, 2018; 7: 1–2.
- 47. Chávez-Jacobo VM, Hernández-Ramírez KC, Romo-Rodríguez P, Pérez-Gallardo RV, Campos-García J, Gutiérrez-Corona JF, et al. CrpP is a novel ciprofloxacin-modifying enzyme encoded by the Pseudomonas aeruginosa pUM505 plasmid, Antimicrob Agents Chemother, 2018; 62: 1–11.
- 48. Nie L, Lv Y, Yuan M, Hu X, Nie T, Yang X, et al. Genetic basis of high level aminoglycoside resistance in Acinetobacter baumanniifrom Beijing, China, Acta Pharm Sin B, 2014; B4: 295–300.
- 49. Poole K. Aminoglycoside resistance in Pseudomonas aeruginosa, Antimicrob Agents Chemother, 2005; 49: 479–487. pmid:15673721
- 50. Subedi D, Vijay AK, Kohli GS, Rice SA, Willcox MJS. Comparative genomics of clinical strains of Pseudomonas aeruginosa strains isolated from different geographic sites, Sci Rep, 2018b; 8:15668.
- 51.
Finch RG, Greenwood D, Whitley RJ, Norrby SR. Antibiotic and chemotherapy e-book. Elsevier Health Sciences, 2010.
- 52. Amina AH, Amin M. Horizontal gene transfer events among different species of bacteria, J Am Sci, 2010; 6: 534–544.
- 53. Bassetti M., Vena A., Croxatto A., Righi E., Guery B., How to manage Pseudomonas aeruginosa infections, Drugs in context, 2018a; 7: 1–18.
- 54.
Fong IW, Shlaes D, Drlica K. Antimicrobial Resistance in the 21st Century. Springer, 2018.
- 55. Vandecraen J, Chandler M, Aertsen A, Van Houdt R. The impact of insertion sequences on bacterial genome plasticity and adaptability, CritRev Microbiol, 2017; 43: 709–730.
- 56. Hõrak R, Kivisaar M. Expression of the Transposase Gene tnpAof Tn4652 Is Positively Affected by Integration Host Factor, J bacteriol, 1998; 180: 2822–2829. pmid:9603867
- 57.
Alkhnbashi O, Meier T, Mitrofanov A, Backofen R, Voß BJM. CRISPR-Cas Bioinformatics, 2019;
- 58. James K.-D. Animal Metabolites: From Amphibians, Reptiles, Aves/Birds, and Invertebrates, Pharmacognosy, 2017; pp. 401–411.
- 59. Rojas MN, Lee X, Waridel P, Maspoli A, Imker HJ, Chai T, et al. The Pseudomonas aeruginosa antimetabolite L -2-amino-4-methoxy-trans-3-butenoic acid (AMB) is made from glutamate and two alanine residues via a thiotemplate-linked tripeptide precursor, Front Microbiol, 2015a; 6: 1–13.
- 60. Rojas MN, Lee X, Waridel P, Maspoli A, Imker HJ, Chai T, et al. The Pseudomonas aeruginosa antimetabolite L-2-amino-4-methoxy-trans-3-butenoic acid (AMB) is made from glutamate and two alanine residues via a thiotemplate-linked tripeptide precursor, Front Microbiol, 2015; 6: 170. pmid:25814981
- 61. Sagot B, Gaysinski M, Mehiri M, Guigonis J-M, Le Rudulier D, Alloing G. Osmotically induced synthesis of the dipeptide N-acetylglutaminylglutamine amide is mediated by a new pathway conserved among bacteria, Proc Natl Acad Sci, 2010; 107: 12652–12657. pmid:20571117
- 62. Kreitler DF, Gemmell EM, Schaffer JE, Wencewicz TA, Gulick AM. The structural basis of N-acyl-α-amino-β-lactone formation catalyzed by a nonribosomal peptide synthetase, Nat Commun, 2019; 10: 1–13. pmid:30602773
- 63. Pierson LS, Pierson EA, biotechnology, Metabolism and function of phenazines in bacteria: impacts on the behavior of bacteria in the environment and biotechnological processes, Appl Microbiol, 2010; 86: 1659–1670.
- 64.
Kastin A. Handbook of biologically active peptides. Academic press, United State of America, 2013.
- 65. Perez RH, Zendo T, Sonomoto K. Novel bacteriocins from lactic acid bacteria (LAB): various structures and applications, Microb Cell fact, 2014; 13: S3. pmid:25186038
- 66. Donia MS, Cimermancic P, Schulze CJ, Brown LCW, Martin J, Mitreva M, et al. et al. A systematic analysis of biosynthetic gene clusters in the human microbiome reveals a common family of antibiotics, Cell Press, 2014; 158, 1402–1414.
- 67. Ranieri MR, Chan DC, Yaeger LN, Rudolph M, Kerabalas-Pittman S, Abdo H, et al. Thiostrepton hijacks pyoverdine receptors to inhibit growth of Pseudomonas aeruginosa, Antimicrob Agents Chemother, 2019; 63: 1–13.
- 68. Sullivan J, Cheng J. Contextualising Chinese Migration to Africa, JAsian Afr Stud, 2018; 53: 1173–1187.