Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Accessory genome of the multi-drug resistant ocular isolate of Pseudomonas aeruginosa PA34

  • Dinesh Subedi ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    d.subedi@unsw.edu.au

    Affiliation School of Optometry and Vision Science, University of New South Wales, Sydney, Australia

  • Gurjeet Singh Kohli,

    Roles Data curation, Formal analysis, Investigation

    Affiliation The Singapore Centre for Environmental Life Sciences Engineering, Nanyang Technological University, Singapore, Singapore

  • Ajay Kumar Vijay,

    Roles Conceptualization, Supervision, Writing – original draft, Writing – review & editing

    Affiliation School of Optometry and Vision Science, University of New South Wales, Sydney, Australia

  • Mark Willcox,

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliation School of Optometry and Vision Science, University of New South Wales, Sydney, Australia

  • Scott A. Rice

    Roles Conceptualization, Resources, Supervision, Writing – review & editing

    Affiliations The Singapore Centre for Environmental Life Sciences Engineering, Nanyang Technological University, Singapore, Singapore, The School of Biological Sciences, Nanyang Technological University, Singapore, Singapore, The ithree Institute, The University of Technology Sydney, Sydney, New South Wales, Australia

Abstract

Bacteria can acquire an accessory genome through the horizontal transfer of genetic elements from non-parental lineages. This leads to rapid genetic evolution allowing traits such as antibiotic resistance and virulence to spread through bacterial communities. The study of complete genomes of bacterial strains helps to understand the genomic traits associated with virulence and antibiotic resistance. We aimed to investigate the complete accessory genome of an ocular isolate of Pseudomonas aeruginosa strain PA34. We obtained the complete genome of PA34 utilising genome sequence reads from Illumina and Oxford Nanopore Technology followed by PCR to close any identified gaps. In-depth genomic analysis was performed using various bioinformatics tools. The susceptibility to heavy metals and cytotoxicity was determined to confirm expression of certain traits. The complete genome of PA34 includes a chromosome of 6.8 Mbp and two plasmids of 95.4 Kbp (pMKPA34-1) and 26.8 Kbp (pMKPA34-2). PA34 had a large accessory genome of 1,213 genes and had 543 unique genes not present in other strains. These exclusive genes encoded features related to metal and antibiotic resistance, phage integrase and transposons. At least 24 genomic islands (GIs) were predicated in the complete chromosome, of which two were integrated into novel sites. Eleven GIs carried virulence factors or replaced pathogenic genes. A bacteriophage carried the aminoglycoside resistance gene (AAC(3)-IId). The two plasmids carried other six antibiotic resistance genes. The large accessory genome of this ocular isolate plays a large role in shaping its virulence and antibiotic resistance.

Introduction

Pseudomonas aeruginosa is associated with many opportunistic and nosocomial human infections such as pneumonia, septicaemia, corneal ulcers (microbial keratitis) and chronic infections in cystic fibrotic lungs [1, 2]. Antibiotic resistance in this bacterium is alarmingly on the rise and P. aeruginosa has been included in one of the top three priority pathogens urgently requiring new antimicrobial therapies for treatment by the World Health Organization [3]. Resistance to different antimicrobials in P. aeruginosa is due to its inherent capacity to oppose the action of antibiotics and its capacity to acquire genetic elements that often carry antibiotic resistance genes [4].

The study of the mechanisms used by different strains of P. aeruginosa to become resistant to antibiotics or acquire virulence have benefited from genome sequencing and these have identified a number of mobile genetic elements (MGEs) [57]. However, many of the genomes investigated are not completely closed and, in this context, accessory genes associated with resistance or virulence in addition to mutations may be missed. This may limit the understanding whether resistance and virulence genes are present on MGEs or “standard” chromosomal regions as well as their prevalence in general [8]. For example, less than 5% of P. aeruginosa draft genomes are complete; furthermore, only 37 plasmids of the total (3003 as of 15/10/2018) in the NCBI P. aeruginosa database are complete. Additionally, many of the isolates sequenced are derived from cystic fibrosis infections, while relatively few ocular and other isolates are represented. This is an important source of genetic information for P. aeruginosa as specific subpopulations are thought to be associated with microbial keratitis [9]. Notably, those subpopulations are predominantly characterised by possession of genes associated with type IV pili twitching motility, cytotoxicity (exoU), and certain genomic islands [9]. Indeed, to date, the complete genome of only one P. aeruginosa keratitis isolate has been published [10]. Very little information is available about the accessory genomes of ocular isolates.

P. aeruginosa strain PA34 (referred to as PA34 hereafter) is a multi-drug resistant microbial keratitis isolate which is resistant to gentamicin, imipenem, ciprofloxacin and moxifloxacin. Analysis of its draft genome revealed that PA34 belongs to sequence type 1284 based on multi-locus sequence typing and carried at least 12 acquired resistance genes and possessed the exoU gene [11], an important effector gene in the pathogenesis of microbial keratitis [12]. A class I integron (In1427) that carries two antibiotic resistance genes has been shown to be integrated into a Tn3-like transposon carried by PA34 [13]. However, the genetic context of these resistance and pathogenic genes has not been completely elucidated. Hence, in this study we sought to analyse the complete genome of PA34, to examine the genetic structure of accessory genome and to compare the complete genome of PA34 with other genomes of P. aeruginosa from the public database.

Materials and methods

Bacterial strain and microbiology

The strain Pseudomonas aeruginosa PA34 was isolated in 1997 from the cornea of a microbial keratitis patient in a tertiary eye care centre L.V. Prasad Eye Institute, Hyderabad, India. The cause of microbial keratitis was recorded as trauma induced by a stone in the right eye of a 21-year old male. Three complete genomes of P. aeruginosa were used for comparison: I) P. aeruginosa PAO1, which has the most complete and curated annotations [14, 15], II) P. aeruginosa PA14, an exoU positive strain [16] and III) P. aeruginosa VRFPA04, an Indian microbial keratitis isolate [10].

The minimum inhibitory concentration (MIC) to three heavy metals (mercury, copper and cobalt) was determined by broth dilution methods as described elsewhere with minor modification [17, 18]. PA34 and PAO1 (as the control strain) were examined against Hg++, Cu++ and Co++ by using twofold serial dilutions from 16mM– 0.31mM. The metal salt solutions (HgCl2, CuCl2 and CoCl2) were prepared in deionised water, sterilized by membrane filtration and added to Mueller-Hinton broth (Oxoid Ltd, Hampshire, UK) to obtain the required concentrations. The experiments were performed in triplicate and repeated three times. Wilcoxon matched-pairs signed rank test was used to confirm significant difference between metal tolerance and the strain. A p-value of <0.05 was taken as significant.

Cytotoxicity was examined by the trypan blue dye exclusion assay [19]. Briefly, human corneal epithelial cells (HCEC) were grown in 96 well plates in the presence of SHEM (DMEM (Gibco, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (Gibco), 1.05mM CaCl2, 0.5% DMSO, 2ng/mL epidermal growth factor (Gibco), 1% ITS-X (Gibco)). 200μL of 5x105 CFU/mL P. aeruginosa in SHEM was used to expose confluent HCEC. Following 3 h incubation at 37°C, the cells were stained with 0.4% trypan blue. The stained cells were observed by microscopy and photographs were taken for the quantitative determination of cytotoxicity. P. aeruginosa PAO1 (a non-cytotoxic strain) were used as a control. The experiment was performed in triplicate and repeated three times.

MinION sequencing and complete genome assembly

Genomic DNA was extracted using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA) as per the manufacturer’s protocol. The amount of extracted DNA was determined using a Qubit fluorometer (Life Technologies, Carlsbad, CA, USA). The library was prepared using the Ligation Sequencing Kit 1D R9 Version as per the manufacturer’s instructions. The sequencing library was loaded into the Flow Cell Mk 1 Spot-ON of the MinION system using Library Loading Bead Kit R9 Version according to the manufacturer's instructions. The raw reads were obtained using the MinKNOW v1.7.14 in a 24-h run experiment and base calling was performed using Albacore v2.0.2. A total of 50,831 reads ranging from 174bp to 58,887bp (mean = 3,292bp) were obtained. The sequence reads from Illumina (which was obtained using MiSeq (Illumina, San Diego, CA) generating 300bp paired-end reads [11]) and MinION data were assembled using three hybrid assembly pipe-lines; Unicycler v0.4.3 [20], Japsa v1.8.0 [21] and Spades v3.12.0 [22]. Spades resulted in the best assembly in terms of contig numbers (Table 1). The genome coverage was calculated using BBmap v35.82 [23].

All contigs were reordered with reference to the three complete genomes of P. aeruginosa. (PAO1, PA14, and VRFPA04) using MUMmer v3.0 [24]. Gaps between the contigs were determined by ABACAS v1.3.1 [25] using the “circular reference genomes” flag. Pseudogenomes of different sizes were observed after comparison with different reference strains and in all comparisons, contigs 6 and 7 were observed as unused contigs indicating they were not part of the chromosome. This was further verified by performing a BLAST of unused contigs against reference genomes. The fewest gaps between contigs in the PA34 genome (99 bp) were observed when PA14 was used as the scaffold compared to the other strains. Further, the size of the PA34 pseudogenome was slightly less than the size of the draft genome (~6.8 Mbp) [11] when using PA14 as the scaffold genome. A gap of 99 bp is an indication of a joining point between two contigs [25]. Therefore, this comparison was considered as the best and was further used to complete the genome of PA34. Primers were designed for each gap (S1 and S2 Tables) using Primer3web v4.1.0 [26] and the gap size was verified by PCR. Using this approach, we were able to close the gaps between contigs to generate a single contiguous chromosome based on contigs 1–5. Contigs 6 and 7 had a lower G+C content than rest of the contigs and showed similarity with different plasmids based BLASTn searches against NCBI database. Contigs 6 and 7 were therefore estimated to represent two different plasmids and were both circularised, referred to here as pMKPA34-1 and pMKPA34-2 (Plasmid 1 and 2 of Microbial Keratitis strain PA34).

Bioinformatics analysis

The complete chromosomal genome was annotated using NCBI Prokaryotic Annotation Pipeline [27]. Plasmids were annotated by Prokka v1.7 [28] followed by manual examination and curation using information from ISsaga [29], the Rapid Annotations using Subsystems Technology (RAST v2.0) [30] and NCBI BLASTn searches. Genomic islands were predicted on the basis of MAUVE [31] whole genome alignment against P. aeruginosa strains PAO1, PA14 and VRFPA04. Through this comparison, DNA blocks of four contiguous open reading frames were predicted as genomic islands of PA34 if these blocks were present in PA34 but absent in anyone of the three reference genomes used for comparison [5]. Pangenome analysis was conducted using Roary v3.12.0 [32], and all four genomes (PAO1, PA14, VRFPA04 and PA34) were annotated using Prokka to avoid annotation biases. Other software used for analysis were CRISPRCasFinder database [33] CGview [34], Easyfig [35], SnapGene Viewer v4.3.2 [36], VENNY v2.1.0 [37] and Resfinder V3.0 [38]. The nucleotide sequence of the complete chromosome and two plasmids were made available in the NCBI database under GenBank accession numbers CP032552, MH547560 and MH547561.

Results and discussion

General features of P. aeruginosa PA34 genome

The statistics of the complete genome of P. aeruginosa PA34 (Table 2) were broadly matched with other published complete genomes of P. aeruginosa [39]. A total of 6,462 coding sequences (CDS) were predicted in the 6.8 Mbp genome of PA34 using NCBI Prokaryotic Genome Annotation Pipeline (PGAP) [40]. Of the predicted CDS, 6,314 were predicted to form functional proteins and 148 were pseudogenes.

In the genomic comparison with three reference genomes (PAO1, PA14 and VRFPA04) and the genome of PA34, it was estimated there were 7,643 orthologs in the pangenome. Of those orthologs, 5078 were common amongst all four strains and are suggested to represent the size of the core genome. PA34 had the largest accessory genome compared to these reference strains (having 1,213 genes [~20% of total genes]) that had 543 genes unique to PA34. These exclusive genes were determined to encode features related to metal and antibiotic resistance, phage integrase and transposons. Furthermore, in the context of individual reference genomes, PA34 has 886, 737, 946 genes with no orthologs in PAO1, PA14, and VRFPA04, respectively (Fig 1) (See S1 File for full annotations of exclusive orthologs). A high proportion of strain-specific accessory genomic elements are common in P. aeruginosa strains and this represents the diversity of its genome [41, 42]. However, in this study, we observed a lower range (3–10%) of exclusive orthologues than previously reported range (10–20%) [41]. The difference observed may be related to the use of the complete genome for the analysis over the use of draft genomes in previous studies. Since draft genomes are likely to be mis-assembled due to the presence of a high content of repetitive sequences in bacteria [43], analysis based on the complete genome is expected to be more relevant. PA34 shared 124 orthologs with VRFPA04, an eye isolate, and those genes were associated with chloramphenicol resistance, pathogenesis (type IV secretion system) and phages.

thumbnail
Fig 1. Graphical representation of genomic comparison of P. aeruginosa strains.

Venn diagram showing the number of common and exclusive orthologs between and amongst four strains of P. aeruginosa. The total number of genes examined of each strain is shown in parentheses. The number were predicted by Roary pangenome analysis and figure was created using VENNY v2.1.

https://doi.org/10.1371/journal.pone.0215038.g001

The diversity of core and accessory genomes observed in this study contradicts findings on environmental species P. fluorescens, where comparative genomics of three isolates demonstrated that they shared a core of 61% genes (vs 66% in this study) and 22–27% of unique genes to each isolate [44]. Furthermore, comparative genomics of three strains of the plant pathogen P. syringae showed a smaller core genome of 56% [45] than that of P. aeruginosa observed in this and another study [46]. The overall low diversity amongst P. aeruginosa strains compared to other Pseudomonas spp. may result from the selection of closely related strains (i.e. those containing exoU) in the present study.

No significant CRISPR and associated Cas genes were observed in PA34 based on comparison with the CRISPRCasFinder database [33]. CRISPR-Cas system may be negatively correlated with the size of the accessory genome and acquired antibiotic resistance genes because this system restricts the invasion and incorporation of MGEs. Thus, strains without a CRISPR-Cas system are expected to have a larger genome size [47].

Genomic islands

Genomic islands (GIs) refer to blocks of horizontally acquired DNA that integrate into certain loci in the core genome [48]. Such loci are known as Regions of Genomic Plasticity (RGPs). GIs should have a minimum of four contiguous open reading frames, which are not present in any of the set of genomes to which they are compared [5]. More than 80 RGPs have been identified in the genome of P. aeruginosa [42]. The position of RGPs in the genome are indicated by homologous flanking loci in PAO1 [5]. RGPs constitute a major portion of the accessory genome and are essential for the adaption of P. aeruginosa in diverse habitats. In this study, 24 RGPs were observed in the complete genome of PA34 (Table 3 and Fig 2) when compared with genomes of P. aeruginosa strains PAO1, PA14 and VRFPA04. A study has reported 27 to 37 RGPs in individual genomes of five P. aeruginosa strains [5]. The difference in number of observed RGPs may be due to a difference in the number of genomes used for comparison. Nevertheless, out of 24 GIs in the current study, two GIs were integrated into loci which have not been reported previously [42] and were deemed as GIs MKPA34-GI1 and MKPA34-GI2 in this study (Table 3). The GI MKPA34-GI1 was 68.9 Kbp, was integrated between PA2858 and PA2859 homologues of PAO1 and carried genes for mercury and chromate resistance. The GI MKPA34-GI2 had a size of 39.9 Kbp, encodes genes associated with phages and was integrated into loci flanking 4856/4857 of PAO1. In addition, three GIs (RGP23, RGP56 and RGP84) carried phage related genes; phages can be associated with virulence [49]. Our analysis showed that at least five GIs were associated with replacement or insertion of pathogenic genes and one GI (RGP23) carried an acquired aminoglycoside resistance gene (AAC(3)-IId). Three GIs (RGP2, RGP50, and RGP57) were associated with the type IV secretion system.

thumbnail
Table 3. Regions of genomic plasticity (RGP) [5, 42] with associated features in PA34.

https://doi.org/10.1371/journal.pone.0215038.t003

thumbnail
Fig 2. The chromosomal map of P. aeruginosa PA34 and position of genomic islands.

Circles from outer to inner represent; circle 1. positions of Coding Sequences (CDS) in plus strand (blue), 2. positions of CDS in minus strand (blue), 3. Location of region of Genomic Plasticity (RGPs) (red) (in the same order as listed in Table 3), 4. BLASTn comparison against PAO1 (green) 5. BLASTn comparison against PA14 (orange), 6. BLASTn comparison against VRFPA04 (purple). 7. G+C content and deviation from the average, 8. G+C skew in green (+) and purple 9. scale in kbp.

https://doi.org/10.1371/journal.pone.0215038.g002

Integrative conjugative elements (ICEs) are chromosomally integrated self-transmissible MGEs that can exist as circular extrachromosomal elements and are antecedents of various GIs [48]. Our analysis revealed three ICEs in the complete genome of PA34. The ICEs observed in RGP5 and RGP29 were related to clc-like ICE elements (Genbank Accession number AJ617740). The parental clc element contains genes for chlorocatechol (clc) degradation [50]. However, these genes were lost from both RGP5 and RGP29 of PA34. Nevertheless, genes required for integration and conjugation were identified in both of the ICEs in RGP5 and RGP29 (Fig 3A) indicating that these elements may have the capacity to transfer within and between species. This may be the reason for observing clc-like GIs in two different loci in this strain. In addition, RGP5 carried mercury resistance genes, which are not present in the parental clc elements. This finding suggests that the clc-like elements may incorporate antibiotic resistance genes during genetic recombination.

thumbnail
Fig 3. Graphical representation of BLASTn comparison of integrative conjugative elements of PA34.

Protein-coding regions are represented by the orange arrows and key features/associated genes are shown. The gradient grey shading represents regions of nucleotide sequence identity (100% to 65%) determined by BLASTn analysis. Figures are drawn to scale using Easyfig [35]. (A) A comparison of two genomic islands (RGP 29 and RGP5) of PA34 with a clc-like integrative conjugative element (Genbank Accession number AJ617740) (parA = plasmid partition protein, traD = conjugative factor, intB13 = phage integrase, merA/B/R = mercury resistance operon, clc = chlorocatechol (clc) degradation). (B) A comparison of the genomic island RGP41 of PA34 with pKLC102 (GenBank Accession number AY257538) and PAPI1(GenBank Accession number AY273869) (soj = chromosome partitioning protein, dnaB = replicative DNA helicase, ssb = single strand binding protein, pilP/pilV = type IV pilus biogenesis and transfer, traD = conjugative factor, cupD/rcsC/pvrR = cell surface fimbriae operon, parE = plasmid stabilisation protein, xerC = integrase).

https://doi.org/10.1371/journal.pone.0215038.g003

The third ICE observed in RGP41 of PA34 was related to the pKLC102 family (GenBank Accession number AY257538) and was similar to P. aeruginosa pathogenicity island PAPI-1 (GenBank Accession number AY273869) first identified in P. aeruginosa clone C (Fig 3B). PAPI-1 is self-transmissible and carries virulence factors such as genes for type IV pilus biogenesis and the cupD gene clusters The cupD genes are essential for the formation of cell surface fimbriae that are involved in biofilm formation [51]. Interestingly, the cupD gene cluster appears to have been lost from the RGP41 of PA34, although the genes for type IV pilus biogenesis are still present.

In addition, various members of pKCL102 ICE family are found to be associated with the carriage of exoU/spcU genes [48]. These GIs have frequently been referred as exoU islands such as exoU island A, exoU island B, exoU island C and PAPI-2 [52, 53]. Possession of the exoU gene, that encodes the type III secretion system effector cytotoxin ExoU, markedly enhances the virulence of P. aeruginosa [54]. An exoU-island of 7.5 Kbp carrying five open reading frames was observed in PA34 (RGP7). In contrast to other exoU-islands that contain genes associated with integration, transposition, conjugations and sometimes antibiotic resistance [55], the exoU island observed here contains only three genes that have homology to a site-specific integrase (Fig 4A). The presence of both PAPI-1 and PAPI-2 like elements that act synergistically in virulence [56] may indicate an enhanced virulence of PA34.

thumbnail
Fig 4. Graphical representation of BLASTn comparison of various genomic islands of PA34.

Protein-coding regions are represented by the orange arrows and key features/associated genes are shown. The gradient grey shading represents regions of nucleotide sequence identity (100% to 79%) determined by BLASTn analysis. Figures are drawn to scale using Easyfig [35]. (A) A comparison of exoU-island (RGP7) of PA34 with four exoU-islands of different strains (exoU island A, exoU island B, exoU island C and PAPI2) (xerC = site specific integrase, exoU = type III secretion system effector protein ExoU, spcU = ExoU chaperon protein, PA0975/PA0989 = franking loci in PAO1) (B) A comparison of a genomic island (MKPA34-GI2) of PA34 with Pseudomonas phage MP38 (GenBank Accession number NC_011611) (traA = transposase A, hni = host nuclease inhibitor, vip = virion protein, tcp = tail component protein).

https://doi.org/10.1371/journal.pone.0215038.g004

To determine whether the exoU gene was functional, cytotoxicity of PA34 was examined in a human corneal epithelium cell line and compared with P. aeruginosa PAO1 which is a non-cytotoxic invasive strain [57]. Microscopic examination after staining with trypan blue indicated that PA34 was highly cytotoxic (Fig 5) showing a large number of dead cells. Based on the percentage of dead cells, maximum cytotoxicity (of score 4) was observed in PA34 [19]. Previous studies have shown that possession of exoU is associated with in vitro and in vivo cytotoxicty, for example by using knockout mutants of P. aeruginosa strains such as PA103 and PAO1 [58, 59]. The results in the current study suggest that exoU in RGP7 of PA34 was transcribed and translated into a functional protein that retained cytotoxicity. The finding also demonstrates that the size of the exoU island was not be associated with activity. However, direct confirmation of this result by producing exoU knockout of P. aeruginosa PA34 is required to confirm this result.

thumbnail
Fig 5. Microscopic photograph of cytotoxic effect of (A) Pseudomonas aeruginosa PAO1 (the control strain) and (B) Pseudomonas aeruginosa PA34 for human corneal epithelium cells strained by trypan blue.

Dead cells are strained with trypan blue.

https://doi.org/10.1371/journal.pone.0215038.g005

A mercury resistance operon was also observed in the MKPA34-GI1 RGP that also carried chromate resistance genes and features related to transposition and conjugation. A BLASTn search against the plasmid database in NCBI did not show significant sequence similarity with MKPA34-GI1. However, 51% of the MKPA34-GI1 sequence was loosely identical (205 matches with an average identity of 97%) to a the complete genome of Pseudomonas stutzeri strain 273 (GenBank accession number CP015641.1). Pseudomonas stutzeri is an environmental organism and has the capacity to degrade pollutants such as toxic metals [60]. These observations indicate that the MKPA34-GI1 may have been acquired by horizontal transfer from organisms in mercury polluted environments. The ability of MKPA34-GI1 to transfer between strains needs to be confirmed by further experiments that will also help to confirm the MKPA34-GI1 RGP loci is mobile.

In order to compare the phenotypic susceptibility to heavy metals, MIC to Hg++, Cu++ and Co++ were examined. High mercury tolerance was observed in PA34 (p<0.05) in comparison to PAO1 (Fig 6), possibly due to the presence of mercury resistance operon in two different GIs (RGP5 and MKPA34-GI1) in PA34. Copper tolerance was slightly higher in PA34 than PAO1 (p>0.05). Although bacteria require copper as a cofactor for many metabolic processes and can tolerate low concentration of Cu++, decreasing sensitivity to copper can be associated with acquisition of copper resistance genes [61]. The presence of a copper resistance operon in RGP23 in PA34 may be associated with higher MIC to Cu++. On the other hand, cobalt tolerance, for which no acquired genes were observed in this study was low in both PA34 and PAO1 (p<0.05). P. aeruginosa can tolerate Cu++ and other heavy metals including Co++ mainly due to P-type ATPase transporter and resistance nodulation cell division efflux pump [62]. In PAO1, the P-type ATPase encoded by PA3920 was found to be important for copper tolerance [63]. In addition, bacterial plasmids also carry metal resistance determinants [64, 65]. However, plasmid associated metal resistance determinates were not detected in PA34.

thumbnail
Fig 6. Minimum inhibitory concentration (MIC) of P. aeruginosa PA34 and P. aeruginosa PAO1 to mercury, copper and cobalt (Asterisk (*) represents statistically significant).

https://doi.org/10.1371/journal.pone.0215038.g006

The Pseudomonas phage MP38 was identified in the MKPA34-GI2 (flanking loci homologous to 4856/4857 of PAO1) (Fig 4B). MP38 is D3112-like phage which is a transposable and has been isolated from many clinical isolates of P. aeruginosa [66]. We examined the presence of MP38 in genomes of different P. aeruginosa and found that a similar phage was integrated into VRFPA04 near the PA4204 homologue of PAO1. This confirms the result of another study [67] that has shown that transposable phages may have variable integration sites, indicating that the phage MP38 does not have any specific integration site in the genome of P. aeruginosa. Therefore, the MKPA34-GI2 integration site may not be a constant RGP for P. aeruginosa.

Four RGPs (RGP31, RGP73, RGP9 and RGP60) of PA34 were associated with replacement of pathogenic genes and were related to lipo-polysaccharide O-antigen, pyoverdine (pvdE), flagella (fliC) and pilus (pilA) synthesis. Despite the replacement islands contains the same genes and being integrated into the same loci in the core genome, they highly diverse between strains [68]. These components are important for interaction of a bacterium with external environments including other species or hosts and hence are under continuous selection pressure [6871]. The pvdE and fliC orthologs have been shown to vary greatly between different strains of P. aeruginosa [72]. The replacement islands may contribute to higher virulence in PA34 although this needs to be tested for swimming or twitching motility and for expression of pvdE under low iron condition.

The largest GI (RGP23) observed in PA34 was 125.1 Kbp and carried resistance genes for aminoglycosides (AAC(3)-IId) and tunicamycin. An insertion sequence of the family IS1182 was observed in this GI, which may be responsible for carriage of these resistance genes. A BLASTn search revealed that the antibiotic resistance genes were best matched with those of Acinetobacter sp. WCHA45 plasmid pNDM1_010045 (GenBank accession number NZ_CP028560.1) and were similar to many from other members of the Enterobacteriaceae where IS1182 is also present. In addition, a phage tail protein gp37, which was first identified in Enterobacter phage T4 [73] was found inserted into this island. This suggests that this resistance island may be derived from phages.

Plasmid features

The Illumina and ONT reads of the whole genome sequence of PA34 were assembled into seven contigs (> 500 bp) using hybrid strategies in SPAdes. Out of seven contigs, contigs 6 and 7 had a G+C content that less than that of other contigs and showed a significant match with plasmids in the NCBI database. This suggests that contigs 6 and 7 are likely to represent two different plasmids carried by PA34 and were named here as pMKPA34-1 and pMKPA34-2. pMKPA34-1 is 95.4 kbp with 57% G+C content and pMKPA34-2 is 26.8 kbp with 61% G+C content. Automatic annotation followed by manual confirmation with BLASTn revealed 98 CDS in pMKPA34-1 and 33 CDS in pMKPA34-2. Out of 98 predicted genes in pMKPA34-1, 46 were predicted to encode proteins with unknown functions whilst all 33 CDS in pMKPA34-2 were predicted to encode known proteins.

General features of pMKPA34-1.

The putative plasmid pMKPA34-1 contains a replication gene (repE), a chromosome partition gene (smc), a plasmid stabilization gene (parB) and a plasmid conjugal transfer mating pair stabilization protein (traN). These genes may help replication, maintenance and transfer of the plasmid [74, 75] (Fig 7). Based on BLASTn searches against the NCBI plasmid database, pMKPA34-1 had the best match (31% query cover with 94% identity) with plasmid pLIB119 of P. stutzeri strain 19SMN4 (GenBank accession number CP007510), followed by Citrobacter freundii strain 18–1 plasmid pBKPC18-1 (23% query cover with 99% identity) (Fig 7). The parA, smc and repE genes display high similarity to the corresponding genes of pLIB119. The presence of the traN gene, which encodes an outer membrane protein and is essential for F-mediated bacterial conjugation [75], indicates the potential exchangeability of pMKPA34-1. This transfer system is common amongst bacterial plasmids notably in Enterobacteriaceae [76] and its presence in pMKPA34-1 may indicate that the plasmid could be exchanged with members of the Enterobacteriaceae.

thumbnail
Fig 7. The map of P. aeruginosa PA34 plasmid pMKPA34-1.

Circles from outer to inner represent; circle 1 and 4: positions of the matched regions representing mobile elements (orange), similarity with other plasmids (green), plasmid functions (marron), and resistance genes (red). Circle 2. Coding Sequence (CDS) in plus strand (purple), Circle 3. positions of CDS in minus strand (purple), Circle 4. G+C content and deviation from the average, Circle 5. G+C skew in green (+) and purple (–) and Circle 6. scale in kbp.

https://doi.org/10.1371/journal.pone.0215038.g007

Transposase, integrase and recombinase genes formed a major portion of the mobile genetic elements in pMKPA34-1. ISsage analysis identified two putative prophage integrases, and site-specific recombinases (xerC and xerD) of the tyrosine family. The XerC-XerD system, which is found on both the chromosome and on plasmids helps separation and segregation of newly replicated bacterial dimeric chromosome by resolving it to monomers [77]. This system is essential in the segregation of multicopy plasmids and contributes to the stable inheritance of multicopy plasmids [78, 79].

Plasmid pMKPA34-1 carries at least six antibiotic resistance genes. Five resistance genes (trimethoprim; dfrA15, chloramphenicol; cmlA1, aminoglycosides; APH(3”)-Ib/APH(6)-Id, and beta-lactam; bla NPS-1) are located in the Tn3-like transposon which also carries a class I integron (In1427) possessing dfrA15 and cmlA1 [13]. Additionally, pMKPA34-1 carries a multi-drug efflux gene (acrB). The acrB gene has high similarity to the corresponding gene of pBKPC18-1 (GenBank accession number CP022275), a resistance plasmid isolated from Citrobacter freundii strain 18–1. This efflux pump is similar to acrA and acrB of E. coli and responsible for resistance to hydrophilic compounds that include disinfectants [80]. The results suggest these movable resistance genes may be related to enteric bacteria.

General features of pMKPA34-2.

The plasmid pMKPA34-2 lacks identifiable replication genes. However, the genome coverage statistics (109x vs. 50x) (Table 1) suggests it may have a higher copy number than the pMKPA34-1. pMKPA34-2 may use alternative mechanisms for replication. In the BLASTn search against the NCBI plasmid database, pMKPA34-2 was best matched (29% query cover with 88% identity) with plasmid pSSE-ATCC-43845 of Salmonella enterica subsp. enterica serovar Senftenberg strain 775WP (GenBank accession number CP016838), followed by plasmid tig00000727 of Klebsiella pneumoniae strain AR_0158 (GenBank accession number CP021699) (14% query cover with 99% identity) (Fig 8). Moreover, pMKPA34-2 carried the series of genes (tnsA, tnsB, tnsC, tnsD and tnsE) that are associated with transposition and are related to Tn7 transposons, a phage integrase and a putative resolvase. Upstream of these mobile genetic elements, a multi drug export protein gene (mepA) was identified. MepA is multi drug efflux transporter of MATE (Multidrug And Toxic Compound Extrusion) family whose role in P. aeruginosa is not well understood.[81]. However, MepA has been associated with resistance to many antibiotics and microbicidal dyes (crystal violet and ethidium bromide) in Staphylococcus aureus [82]. Indeed, the presence of the mepA gene in this plasmid and its association with mobile genetic elements suggest that this gene may have been acquired through horizontal gene transfer.

thumbnail
Fig 8. The map of P. aeruginosa PA34 plasmid pMKPA34-2.

Circles from outer to inner represent; circle 1 and 4: positions of the matched regions representing mobile elements (orange), similarity with other plasmids (green) and resistance genes (red). Circle 2. CDS in plus strand(purple), Circle 3. positions of CDS in minus strand (purple), Circle 4. G+C content and deviation from the average, Circle 5. G+C skew in green (+) and purple (–) and Circle 6. scale in kbp.

https://doi.org/10.1371/journal.pone.0215038.g008

Conclusions

The large accessory genome of P. aeruginosa strain PA34 indicates that this strain has a diverse genomic structure. The strain harbours twenty-four genomic islands and two plasmids that carry various metal and antibiotic resistance genes as well as several genes associated with virulence. This may be associated with the observance of higher antibiotic resistance, mercury tolerance and in-vitro cytotoxicity of PA34. The in-silico analysis showed that six antibiotic resistance genes were present in two different plasmids and one antibiotic resistance gene plus various mercury resistance genes were integrated into different GIs and these resistance genes have sequence similarities with that of either other environmental or enteric bacteria. Furthermore, the genome of the PA34 has been integrated by phage element (gp37) that has its origin in enteric bacteria and carried aminoglycoside resistance gene (AAC(3)-IId). These findings suggest that resistance and virulence in PA34 may have evolved due to environmental selection pressure where organisms acquire traits to survive predation by other inhabitants. Thus acquired traits enhance the pathogenesis process in human. Given the eye is susceptible to being infected by environmental strains of P. aeruginosa, examination of a larger number of eye isolates may be necessary to uncover any additional acquired genetic features associated with microbial keratitis. This will help our understanding of different aspects of Pseudomonas keratitis.

Supporting information

S1 Table. Gaps in pseudogenome of PA34 with comparison with PA14.

https://doi.org/10.1371/journal.pone.0215038.s001

(DOCX)

S2 Table. Primer proposed for the gap validation.

https://doi.org/10.1371/journal.pone.0215038.s002

(DOCX)

S1 File. Annotation of exclusive orthologs of strains.

https://doi.org/10.1371/journal.pone.0215038.s003

(XLS)

Acknowledgments

The authors would like to acknowledge Associate Professor Federico Lauro of the Asian School of the Environment and the Singapore Centre for Environmental Life Sciences Engineering (SCELSE) for providing Oxford Nanopore Technology facility for DNA sequencing. We are also thankful to UNSW high performance computing facility KATANA for providing us cluster time for data analysis. Thanks to Sidra Sarwat UNSW for her help with cytotoxicity assay.

References

  1. 1. Lyczak JB, Cannon CL, Pier GB. Establishment of Pseudomonas aeruginosa infection: lessons from a versatile opportunist. Microbes Infect. 2000;2(9):1051–60. pmid:10967285.
  2. 2. Stapleton F, Dart JK, Seal DV, Matheson M. Epidemiology of Pseudomonas aeruginosa keratitis in contact lens wearers. Epidemiol Infect. 1995;114(3):395–402. Epub 1995/06/01. pmid:7781727.
  3. 3. WHO priority pathogens list for R&D of new antibiotics [Internet]. 2017.
  4. 4. Livermore DM. Multiple mechanisms of antimicrobial resistance in Pseudomonas aeruginosa: our worst nightmare? Clin Infect Dis. 2002;34(5):634–40. pmid:11823954.
  5. 5. Mathee K, Narasimhan G, Valdes C, Qiu X, Matewish JM, Koehrsen M, et al. Dynamics of Pseudomonas aeruginosa genome evolution. Proc Natl Acad Sci U S A. 2008;105(8):3100–5. pmid:18287045.
  6. 6. Jeukens J, Boyle B, Kukavica-Ibrulj I, Ouellet MM, Aaron SD, Charette SJ, et al. Comparative genomics of isolates of a Pseudomonas aeruginosa epidemic strain associated with chronic lung infections of cystic fibrosis patients. PLoS One. 2014;9(2):e87611. Epub 2014/02/08. pmid:24505294.
  7. 7. Stewart L, Ford A, Sangal V, Jeukens J, Boyle B, Kukavica-Ibrulj I, et al. Draft genomes of 12 host-adapted and environmental isolates of Pseudomonas aeruginosa and their positions in the core genome phylogeny. Pathog Dis. 2014;71(1):20–5. pmid:24167005.
  8. 8. Beatson SA, Walker MJ. Microbiology. Tracking antibiotic resistance. Science. 2014;345(6203):1454–5. pmid:25237090.
  9. 9. Stewart RM, Wiehlmann L, Ashelford KE, Preston SJ, Frimmersdorf E, Campbell BJ, et al. Genetic characterization indicates that a specific subpopulation of Pseudomonas aeruginosa is associated with keratitis infections. J Clin Microbiol. 2011;49(3):993–1003. pmid:21227987.
  10. 10. Murugan N, Malathi J, Umashankar V, Madhavan HN. Unraveling genomic and phenotypic nature of multidrug-resistant (MDR) Pseudomonas aeruginosa VRFPA04 isolated from keratitis patient. Microbiol Res. 2016;193:140–9. Epub 2016/11/09. pmid:27825482.
  11. 11. Subedi D, Vijay AK, Kohli GS, Rice SA, Willcox M. Comparative genomics of clinical strains of Pseudomonas aeruginosa strains isolated from different geographic sites. Sci Rep. 2018;8(1):15668. Epub 2018/10/26. pmid:30353070.
  12. 12. Fleiszig SM, Lee EJ, Wu C, Andika RC, Vallas V, Portoles M, et al. Cytotoxic strains of Pseudomonas aeruginosa can damage the intact corneal surface in vitro. Clao j. 1998;24(1):41–7. pmid:9474453
  13. 13. Subedi D, Vijay AK, Kohli GS, Rice SA, Willcox M. Nucleotide sequence analysis of NPS-1 beta-lactamase and a novel integron (In1427)-carrying transposon in an MDR Pseudomonas aeruginosa keratitis strain. J Antimicrob Chemother. 2018. Epub 2018/03/17. pmid:29547943.
  14. 14. Klockgether J, Munder A, Neugebauer J, Davenport CF, Stanke F, Larbig KD, et al. Genome diversity of Pseudomonas aeruginosa PAO1 laboratory strains. J Bacteriol. 2010;192(4):1113–21. pmid:20023018.
  15. 15. Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, Hickey MJ, et al. Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen. Nature. 2000;406(6799):959–64. pmid:10984043.
  16. 16. Lee DG, Urbach JM, Wu G, Liberati NT, Feinbaum RL, Miyata S, et al. Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial. Genome Biol. 2006;7(10):R90. Epub 2006/10/14. pmid:17038190.
  17. 17. Clinical and Laboratory Standards Institute (CLSI). Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically; approved standard—Ninth Edition. CLSI. 2012;32(M07-A9).
  18. 18. Teitzel GM, Parsek MR. Heavy metal resistance of biofilm and planktonic Pseudomonas aeruginosa. Appl Environ Microb. 2003;69(4):2313–20. pmid:12676715.
  19. 19. Fleiszig SM, Zaidi TS, Preston MJ, Grout M, Evans DJ, Pier GB. Relationship between cytotoxicity and corneal epithelial cell invasion by clinical isolates of Pseudomonas aeruginosa. Infect Immun. 1996;64(6):2288–94. Epub 1996/06/01. pmid:8675339.
  20. 20. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS computational biology. 2017;13(6):e1005595. pmid:28594827
  21. 21. Nguyen SH, Duarte TP, Coin LJ, Cao MD. Real-time demultiplexing Nanopore barcoded sequencing data with npBarcode. Bioinformatics. 2017;33(24):3988–90. pmid:28961965
  22. 22. 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(5):455–77. pmid:22506599.
  23. 23. Bushnell B. BBMap short read aligner. University of California, Berkeley, California URL http://sourceforgenet/projects/bbmap. 2016.
  24. 24. Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biology. 2004;5(2):R12. pmid:14759262
  25. 25. Assefa S, Keane TM, Otto TD, Newbold C, Berriman M. ABACAS: algorithm-based automatic contiguation of assembled sequences. Bioinformatics. 2009;25(15):1968–9. pmid:19497936
  26. 26. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115–e. pmid:22730293
  27. 27. Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Research. 2016;44(14):6614–24. pmid:27342282
  28. 28. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9. pmid:24642063.
  29. 29. Varani AM, Siguier P, Gourbeyre E, Charneau V, Chandler M. ISsaga is an ensemble of web-based methods for high throughput identification and semi-automatic annotation of insertion sequences in prokaryotic genomes. Genome Biol. 2011;12(3):R30. pmid:21443786
  30. 30. 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. 2014;42(D1):D206–D14. pmid:24293654
  31. 31. Darling AE, Mau B, Perna NT. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5(6):e11147. pmid:20593022.
  32. 32. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(22):3691–3. pmid:26198102
  33. 33. Couvin D, Bernheim A, Toffano-Nioche C, Touchon M, Michalik J, Néron B, et al. CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res. 2018;46(W1):W246–W51. pmid:29790974
  34. 34. Grant JR, Arantes AS, Stothard P. Comparing thousands of circular genomes using the CGView Comparison Tool. BMC Genomics. 2012;13:202. pmid:22621371.
  35. 35. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. Epub 2011/02/01. pmid:21278367.
  36. 36. Biotech G. SnapGene Viewer. Glick B, editor. 3(3).
  37. 37. Venny. An interactive tool for comparing lists with Venn's diagrams [Internet]. 2007–2015. http://bioinfogp.cnb.csic.es/tools/venny/index.html.
  38. 38. 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(11):2640–4. pmid:22782487.
  39. 39. Winsor GL, Griffiths EJ, Lo R, Dhillon BK, Shay JA, Brinkman FS. Enhanced annotations and features for comparing thousands of Pseudomonas genomes in the Pseudomonas genome database. Nucleic Acids Res. 2016;44(D1):D646–53. Epub 2015/11/19. pmid:26578582.
  40. 40. Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44(14):6614–24. Epub 2016/06/28. pmid:27342282.
  41. 41. Pohl S, Klockgether J, Eckweiler D, Khaledi A, Schniederjans M, Chouvarine P, et al. The extensive set of accessory Pseudomonas aeruginosa genomic components. FEMS Microbiol Lett. 2014;356(2):235–41. pmid:24766399.
  42. 42. Klockgether J, Cramer N, Wiehlmann L, Davenport CF, Tummler B. Pseudomonas aeruginosa genomic structure and diversity. Front Microbiol. 2011;2:150. pmid:21808635.
  43. 43. Ricker N, Qian H, Fulthorpe RR. The limitations of draft assemblies for understanding prokaryotic adaptation and evolution. Genomics. 2012;100(3):167–75. https://doi.org/10.1016/j.ygeno.2012.06.009. pmid:22750556
  44. 44. Silby MW, Cerdeno-Tarraga AM, Vernikos GS, Giddens SR, Jackson RW, Preston GM, et al. Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens. Genome Biol. 2009;10(5):R51. pmid:19432983.
  45. 45. Almeida NF, Yan S, Lindeberg M, Studholme DJ, Schneider DJ, Condon B, et al. A draft genome sequence of Pseudomonas syringae pv. tomato T1 reveals a type III effector repertoire significantly divergent from that of Pseudomonas syringae pv. tomato DC3000. Mol Plant Microbe Interact. 2009;22(1):52–62. Epub 2008/12/09. pmid:19061402.
  46. 46. Ozer EA, Allen JP, Hauser AR. Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt. BMC Genomics. 2014;15(1):737. pmid:25168460
  47. 47. van Belkum A, Soriaga LB, LaFave MC, Akella S, Veyrieras JB, Barbu EM, et al. Phylogenetic Distribution of CRISPR-Cas Systems in Antibiotic-Resistant Pseudomonas aeruginosa. MBio. 2015;6(6):e01796–15. Epub 2015/11/26. pmid:26604259.
  48. 48. Kung VL, Ozer EA, Hauser AR. The accessory genome of Pseudomonas aeruginosa. Microbiol Mol Biol Rev. 2010;74(4):621–41. pmid:21119020.
  49. 49. Winstanley C, Langille MG, Fothergill JL, Kukavica-Ibrulj I, Paradis-Bleau C, Sanschagrin F, et al. Newly introduced genomic prophage islands are critical determinants of in vivo competitiveness in the Liverpool Epidemic Strain of Pseudomonas aeruginosa. Genome Res. 2009;19(1):12–23. Epub 2008/12/03. pmid:19047519.
  50. 50. Ravatn R, Studer S, Springael D, Zehnder AJ, van der Meer JR. Chromosomal integration, tandem amplification, and deamplification in Pseudomonas putida F1 of a 105-kilobase genetic element containing the chlorocatechol degradative genes from Pseudomonas sp. Strain B13. J Bacteriol. 1998;180(17):4360–9. pmid:9721270.
  51. 51. Mikkelsen H, Ball G, Giraud C, Filloux A. Expression of Pseudomonas aeruginosa CupD fimbrial genes is antagonistically controlled by RcsB and the EAL-containing PvrR response regulators. PLoS One. 2009;4(6):e6018. Epub 2009/06/24. pmid:19547710.
  52. 52. Kulasekara BR, Kulasekara HD, Wolfgang MC, Stevens L, Frank DW, Lory S. Acquisition and evolution of the exoU locus in Pseudomonas aeruginosa. J Bacteriol. 2006;188(11):4037–50. pmid:16707695.
  53. 53. He J, Baldini RL, Deziel E, Saucier M, Zhang Q, Liberati NT, et al. The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes. Proc Natl Acad Sci U S A. 2004;101(8):2530–5. pmid:14983043.
  54. 54. Hauser AR, Cobb E, Bodi M, Mariscal D, Valles J, Engel JN, et al. Type III protein secretion is associated with poor clinical outcomes in patients with ventilator-associated pneumonia caused by Pseudomonas aeruginosa. Crit Care Med. 2002;30(3):521–8. pmid:11990909.
  55. 55. Roy Chowdhury P, Scott M, Worden P, Huntington P, Hudson B, Karagiannis T, et al. Genomic islands 1 and 2 play key roles in the evolution of extensively drug-resistant ST235 isolates of Pseudomonas aeruginosa. Open Biol. 2016;6(3). pmid:26962050.
  56. 56. Harrison EM, Carter ME, Luck S, Ou HY, He X, Deng Z, et al. Pathogenicity islands PAPI-1 and PAPI-2 contribute individually and synergistically to the virulence of Pseudomonas aeruginosa strain PA14. Infect Immun. 2010;78(4):1437–46. Epub 2010/02/04. pmid:20123716.
  57. 57. Evans DJ, Frank DW, Finck-Barbancon V, Wu C, Fleiszig SM. Pseudomonas aeruginosa invasion and cytotoxicity are independent events, both of which involve protein tyrosine kinase activity. Infect Immun. 1998;66(4):1453–9. pmid:9529067.
  58. 58. Finck-Barbancon V, Goranson J, Zhu L, Sawa T, Wiener-Kronish JP, Fleiszig SM, et al. ExoU expression by Pseudomonas aeruginosa correlates with acute cytotoxicity and epithelial injury. Mol Microbiol. 1997;25(3):547–57. pmid:9302017.
  59. 59. Allewelt M, Coleman FT, Grout M, Priebe GP, Pier GB. Acquisition of expression of the Pseudomonas aeruginosa ExoU cytotoxin leads to increased bacterial virulence in a murine model of acute pneumonia and systemic spread. Infect Immun. 2000;68(7):3998. pmid:10858214
  60. 60. Lalucat J, Bennasar A, Bosch R, García-Valdés E, Palleroni NJ. Biology of Pseudomonas stutzeri. Microbiol Mol Biol Rev. 2006;70(2):510–47. pmid:16760312.
  61. 61. Elguindi J, Wagner J, Rensing C. Genes involved in copper resistance influence survival of Pseudomonas aeruginosa on copper surfaces. J Appl Microbiol. 2009;106(5):1448–55. Epub 02/23. pmid:19239551.
  62. 62. Nies DH. Efflux-mediated heavy metal resistance in prokaryotes. FEMS Microbiol Rev. 2003;27(2–3):313–39. pmid:12829273.
  63. 63. Teitzel GM, Geddie A, De Long SK, Kirisits MJ, Whiteley M, Parsek MR. Survival and growth in the presence of elevated copper: transcriptional profiling of copper-stressed Pseudomonas aeruginosa. J Bacteriol. 2006;188(20):7242–56. Epub 2006/10/04. pmid:17015663.
  64. 64. Ohtake H, Cervantes C, Silver S. Decreased chromate uptake in Pseudomonas fluorescens carrying a chromate resistance plasmid. J Bacteriol. 1987;169(8):3853–6. Epub 1987/08/01. pmid:3112130.
  65. 65. Brown NL, Misra TK, Winnie JN, Schmidt A, Seiff M, Silver S. The nucleotide sequence of the mercuric resistance operons of plasmid R100 and transposon Tn501: further evidence for mer genes which enhance the activity of the mercuric ion detoxification system. Mol Gen Genet. 1986;202(1):143–51. Epub 1986/01/01. pmid:3007931.
  66. 66. Wang PW, Chu L, Guttman DS. Complete sequence and evolutionary genomic analysis of the Pseudomonas aeruginosa transposable bacteriophage d3112. J Bacteriol. 2004;186(2):400–10. pmid:14702309.
  67. 67. Vergnaud G, Midoux C, Blouin Y, Bourkaltseva M, Krylov V, Pourcel C. Transposition behavior revealed by high-resolution description of Pseudomonas aeruginosa Saltovirus integration sites. Viruses. 2018;10(5). pmid:29735891.
  68. 68. Smith EE, Sims EH, Spencer DH, Kaul R, Olson MV. Evidence for diversifying selection at the pyoverdine locus of Pseudomonas aeruginosa. J Bacteriol. 2005;187(6):2138–47. Epub 2005/03/04. pmid:15743962.
  69. 69. Chibeu A, Ceyssens P-J, Hertveldt K, Volckaert G, Cornelis P, Matthijs S, et al. The adsorption of Pseudomonas aeruginosa bacteriophage φKMV is dependent on expression regulation of type IV pili genes. FEMS Microbiol Lett. 2009;296(2):210–8. pmid:19459952
  70. 70. Kuzio J, Kropinski AM. O-antigen conversion in Pseudomonas aeruginosa PAO1 by bacteriophage D3. J Bacteriol. 1983;155(1):203–12. pmid:6190794
  71. 71. Pier GB. Pseudomonas aeruginosa lipopolysaccharide: a major virulence factor, initiator of inflammation and target for effective immunity. Int J Med Microbiol. 2007;297(5):277–95. pmid:17466590
  72. 72. Subedi D, Vijay AK, Kohli GS, Rice SA, Willcox M. Comparative genomics of clinical strains of Pseudomonas aeruginosa strains isolated from different geographic sites. Sci Rep. 2018;8(1):15668. pmid:30353070
  73. 73. Cerritelli ME, Wall JS, Simon MN, Conway JF, Steven AC. Stoichiometry and domainal organization of the long tail-fiber of bacteriophage T4: a hinged viral adhesin. J Mol Biol. 1996;260(5):767–80. Epub 1996/08/02. pmid:8709154.
  74. 74. del Solar G, Giraldo R, Ruiz-Echevarría MJ, Espinosa M, Díaz-Orejas R. Replication and Control of Circular Bacterial Plasmids. Microbiol Mol Biol Rev. 1998;62(2):434–64. pmid:9618448.
  75. 75. Klimke WA, Rypien CD, Klinger B, Kennedy RA, Rodriguez-Maillard JM, Frost LS. The mating pair stabilization protein, TraN, of the F plasmid is an outer-membrane protein with two regions that are important for its function in conjugation. Microbiology. 2005;151(Pt 11):3527–40. pmid:16272376.
  76. 76. Lawley TD, Klimke WA, Gubbins MJ, Frost LS. F factor conjugation is a true type IV secretion system. FEMS Microbiol Lett. 2003;224(1):1–15. pmid:12855161.
  77. 77. Castillo F, Benmohamed A, Szatmari G. Xer site specific recombination: Double and single recombinase systems. Front Microbiol. 2017;8:453. Epub 2017/04/05. pmid:28373867.
  78. 78. Arciszewska LK, Sherratt DJ. Xer site-specific recombination in vitro. EMBO J. 1995;14(9):2112–20. pmid:7744017.
  79. 79. Blakely GW, Davidson AO, Sherratt DJ. Sequential strand exchange by XerC and XerD during site-specific recombination at dif. J Biol Chem. 2000;275(14):9930–6. Epub 2000/04/01. pmid:10744667.
  80. 80. Ma D, Cook DN, Alberti M, Pon NG, Nikaido H, Hearst JE. Genes acrA and acrB encode a stress-induced efflux system of Escherichia coli. Mol Microbiol. 1995;16(1):45–55. pmid:7651136.
  81. 81. Kuroda T, Tsuchiya T. Multidrug efflux transporters in the MATE family. Biochimica et Biophysica Acta (BBA)—Proteins and Proteomics. 2009;1794(5):763–8. https://doi.org/10.1016/j.bbapap.2008.11.012.
  82. 82. Kaatz GW, McAleese F, Seo SM. Multidrug resistance in Staphylococcus aureus due to overexpression of a novel Multidrug and Toxin Extrusion (MATE) transport protein. Antimicrob Agents Chemother. 2005;49(5):1857–64. pmid:15855507