Next Article in Journal
Is Gut Microbiota Dysbiosis a Predictor of Increased Susceptibility to Poor Outcome of COVID-19 Patients? An Update
Next Article in Special Issue
Whole Genome Sequence-Based Prediction of Resistance Determinants in High-Level Multidrug-Resistant Campylobacter jejuni Isolates in Lithuania
Previous Article in Journal
Comparison of Microbiological Characteristics and Genetic Diversity between Burkholderia cepacia Complex Isolates from Vascular Access and Other Clinical Infections
Previous Article in Special Issue
Integrative and Conjugative Element ICETh1 Functions as a Pangenomic DNA Capture Module in Thermus thermophilus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plasmid Identification and Plasmid-Mediated Antimicrobial Gene Detection in Norwegian Isolates

by
Abdolrahman Khezri
1,
Ekaterina Avershina
1,2 and
Rafi Ahmad
1,3,4,*
1
Department of Biotechnology, Inland Norway University of Applied Sciences, 2318 Hamar, Norway
2
Laboratory or Postgenomic Technologies, Izmerov Research Institute of Occupational Health, 105275 Moscow, Russia
3
Division of Medical Services-Clinical Microbiology, Inland Hospital, 2609 Lillehammer, Norway
4
Institute of Clinical Medicine, Faculty of Health Sciences, UiT—The Arctic University of Norway, Hansine Hansens veg 18, 9019 Tromsø, Norway
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(1), 52; https://doi.org/10.3390/microorganisms9010052
Submission received: 23 November 2020 / Revised: 22 December 2020 / Accepted: 23 December 2020 / Published: 27 December 2020
(This article belongs to the Special Issue Bacterial Genomes and Evolution by Horizontal Gene Transfer (HGT))

Abstract

:
Norway is known for being one of the countries with the lowest levels of antimicrobial resistance (AMR). AMR, through acquired genes located on transposons or conjugative plasmids, is the horizontal transmission of genes required for a given bacteria to withstand antibiotics. In this work, bioinformatic analysis of whole-genome sequences and hybrid assembled data from Escherichia coli, and Klebsiella pneumoniae isolates from Norwegian patients was performed. For detection of putative plasmids in isolates, the plasmid assembly mode in SPAdes was used, followed by annotation of resulting contigs using PlasmidFinder and two curated plasmid databases (Brooks and PLSDB). Furthermore, ResFinder and Comprehensive Antibiotic Resistance Database (CARD) were used for the identification of antibiotic resistance genes (ARGs). The IncFIB plasmid was detected as the most prevalent plasmid in both E. coli, and K. pneumoniae isolates. Furthermore, ARGs such as aph(3″)-Ib, aph(6)-Id, sul1, sul2, tet(D), and qnrS1 were identified as the most abundant plasmid-mediated ARGs in Norwegian E. coli and K. pneumoniae isolates, respectively. Using hybrid assembly, we were able to locate plasmids and predict ARGs more confidently. In conclusion, plasmid identification and ARG detection using whole-genome sequencing data are heavily dependent on the database of choice; therefore, it is best to use several tools and/or hybrid assembly for obtaining reliable identification results.

1. Introduction

Antimicrobial resistance (AMR) is the ability of microorganisms to resist antimicrobial treatments, especially antibiotics. Infections due to AMR bacteria are a threat to modern health care and are responsible for an estimated 700,000 and 33,000 deaths/year globally and in Europe, respectively [1]. Recently, the World Health Organization (WHO) published a list of pathogens for which urgent global action is needed [2]. Extended spectrum β-lactamases (ESBL) producing and carbapenem-resistant Enterobacteriaceae are among the priority one critical section of the WHO pathogen list. There has been a global rise in infections caused by multi-drug resistant clones of Enterobacteriaceae, particularly Klebsiella pneumoniae and Escherichia coli [3].
AMR can arise through various mechanisms, including mutations of chromosomal genes and the acquisition of antibiotic resistance genes (ARGs) from other strains in a process termed horizontal gene transfer (HGT). It is the sharing of genes through HGT that has largely contributed to the global dissemination of ARGs [4]. The genomes of E. coli and K. pneumoniae are prone to a mutation in stress, depicting the genetic flexibility to upregulate their natural resistance and acquire foreign determinants through HGT due to mobile genetic elements. These elements, such as plasmids, transposons, integrons, and genomic islands, harbor ARGs [5]. Several plasmids like IncF and IncI1 plasmids are known to carry resistance genes in E. coli, K. pneumoniae, and other Enterobacteriaceae [6]. Additionally, the ColE plasmids, which encode colicins, and have killing activity against different bacteria, are also important plasmids [7]. Broad range resistance plasmids are known to be associated with pathogens; for example, a resistance plasmid from Enterobacteriaceae can be transferred to a wide variety of Gram-negative organisms.
Whole-genome sequencing (WGS) is an effective method of tracking the onward transmission of bacteria or resistance plasmid transfer between bacteria. It has made it possible to determine and evaluate an organism’s whole DNA sequence at low costs in a short period of time. It allows for the identification of antimicrobial resistance and the early detection of outbreaks or their epidemiological investigation [8]. Moreover, plasmid assembly and characterization following WGS is a difficult task. This happens because the plasmids tend to contain repeat sequences with sizes greater than sequences generated by sequencing platforms such as Illumina technology (San Diego, CA, USA) [9]. Therefore, the need for in silico plasmid detection has emerged due to the difficulty of plasmid DNA purification if they are longer than 50 kbp [10]. In addition, regarding the need for an efficient plasmid identification tool, ARG databases with comprehensive and accurate gene records are needed to assess AMR prevalence. Although several ARG databases exist, Comprehensive Antibiotic Resistance Database (CARD), and ResFinder are the most effective and have sustainable curation strategies [11]. Recent studies have shown that the hybrid assemblies, which are a combination of Illumina and long-read sequencing (e.g., Oxford Nanopore Technology’s MinION) data, are better at identifying plasmids and ARGs [12]. However, this requires advanced bioinformatic and machine learning methods for WGS data analysis [13,14,15].
Globally, AMR is unevenly distributed. Recently, Klein et al. investigated the drug resistance index (DRI) for 41 countries [16]. They have compared the reported data on antibiotics’ use and their resistance to the treatment of infections caused by microorganisms from the WHO priority list [2]. Norway is among the countries with the lowest DRIs (third lowest), and has a DRI value around four-fold lower than that of the country with the highest DRI; India. However, there is an increasing trend in AMR cases in Norway. For example, the percentage of E. coli with ESBL, causing septicemia, has increased ten-fold in the last ten years [17]. There has also been a slightly increased prevalence of ESBLs for E. coli (6.6% in 2017 and 6.5% in 2018) and Klebsiella spp. (5.3% in 2017 and 6.6% in 2018) [17].
This research has utilized different tools and databases to identify plasmids and predict plasmid-mediated ARGs in both E. coli and K. pneumoniae isolates. Our results indicate that plasmid identification and ARGs prediction are database/tool dependent. In this regard, a hybrid assembly can be considered an efficient way to identify plasmids and predict plasmid-mediated ARGs.

2. Materials and Methods

2.1. Sample Collection and Characterization

In this study, E. coli and K. pneumoniae isolates were collected from blood and urine specimens of Norwegian patients, in collaboration with Oslo University Hospital. The sample overview is in Table 1.

2.2. Library Preparation and Whole-Genome Sequencing

The WGS data used in this study are from our recent work, which was performed at Oslo University Hospital [18]. In brief, DNA was isolated from bacteria colonies using QIAamp DNA minikit (Qiagen, Hilden, Germany) following the manufacturer’s instructions and was quantified using Qubit fluorometer (Life Technologies, Carlsbad, CA, USA). The libraries were constructed using the Nextera XT kit (Illumina Ltd., San Diego, CA, USA) according to the manufacturer’s recommendations. The libraries were sequenced in pair-end mode (2 × 300 bp) using the Illumina MiSeq platform at the Norwegian Sequencing Center (Oslo, Norway). Furthermore, to make a hybrid assembly, we sequenced three more isolates (E. coli 39, K. pneumoniae 23, and 27) using the nanopore and the Illumina MiSeq sequencing platforms. Details regarding library preparation, sequencing, and hybrid assembly have previously been reported [12]. All bioinformatic analyses for both plasmid and hybrid assemblies were identical and performed as described below.

2.3. Bioinformatic Analyses of Bacterial Genomic

2.3.1. Quality Control and Trimming of Illumina Sequences

Initially, Illumina and Nanopore reads were quality checked using FastQC (v 0.11.8 for Linux) [19]. Then Illumina adapters were removed, and low-quality reads (Phred below 25) were filtered out using Trimmomatic with default parameters [20]. Before downstream analyses, trimmed reads were again quality checked using FastQC software.

2.3.2. Plasmid Assembly and in Silico Plasmid Identification

Putative plasmid sequences were assembled using plasmid flag in SPAdes (v 3.14.1 for Linux) [21]. General statistics of the assembled putative plasmids was assessed using QUAST (v 4.6.0 for Linux) [22]. Putative plasmid sequences were further confirmed using PlasmidFinder (software version: 2.0.1, database version: 2020-07-13) with minimum identity and coverage of 95% and 60%, respectively [23]. In addition to PlasmidFinder, the identification of putative plasmids was performed using two other methods. First, plasmid reference sequences were downloaded from the PLSDB database [24] and a curated database developed by Brooks et al. [25], hereafter referred to as Brooks. Later, assembled putative plasmids were BLAST searched (sequence identity >95% and word size 28) against downloaded reference plasmids databases. Initially, hits (contigs) with coverage between 30 to 100% were extracted and utilized for the next step. Then only hits with qcov ≥90% (PTRUE) were considered for downstream analysis. The qcov is unique query coverage per subject, calculated after considering any alignment overlaps between different fragments aligned with that specific subject in the database.

2.3.3. Identification of Plasmids Mediated Antimicrobial Resistance Genes (ARGs)

To identify ARGs hosted by plasmids, the assembled putative plasmids for each isolate were submitted to Resfinder 4.0 [26] and resistance gene identifier tool from Comprehensive Antibiotic Resistance Database—CARD [27]. In both Resfinder and CARD, only hits showing ≥95% identity and ≥98% length coverage were considered as ARGs. Later, only hits sharing the same contig with PTRUE were regarded as true plasmid-mediated ARGs (ARGPlasmid). Hereafter we refer to ARGPlasmid-PlasmidFinder (meaning ARG and PTRUE from PlasmidFinder were found on the same contig), ARGPlasmid-Brooks (meaning ARG and PTRUE from Brooks were found on the same contig), ARGPlasmid-PLSDB (meaning that ARG and PTRUE from PLSDB were found on the same contig).

3. Results

In the current study, putative plasmid sequences and ARGs were identified in silico in E. coli and K. pneumoniae isolates from Norwegian patients. Furthermore, hybrid assemblies from additional three isolates were also analyzed.

3.1. General Statistics of Assembled Plasmid Sequences and Hybrid Assembled Sequences

The general statistics for assembled plasmid sequences and hybrid assembled sequences are shown in Table 2. We observed a higher number and bigger contig size for E. coli than K. pneumoniae isolates. The GC percentage between E. coli and K. pneumoniae was similar, and the N50 values (i.e., the minimum contig length required to cover 50% of the assembled genome sequence) were higher in E. coli, indicating larger contig size of plasmids which denotes good quality of assembly. Regarding hybrid assembled isolates, generally, bigger contigs and higher N50 values were observed (Table 2). Interestingly GC percentage was higher in both K. pneumoniae isolates compared to the E. coli isolate.

3.2. In Silico Plasmid Validation

Assembled plasmid sequences were further validated with the PlasmidFinder online tool and using BLASTn against PLSDB and Brooks plasmid database. Using the PlasmidFinder tool, we identified plasmid replicons in 39 (67%) and 11 (25%) of E. coli and K. pneumoniae isolates. This corresponds to two to three and one to two plasmid replicons per isolate caring plasmids in E. coli and K. pneumoniae, respectively. (Table S1).
The number of putative plasmids (PTRUE) after BLASTn and removing duplicates hits per isolate for both E. coli and K. pneumoniae is shown in Figure 1A. Overall a higher number of PTRUE was detected for E. coli than K. pneumoniae (Table S1). For E. coli, the majority of PTRUE (122 of 173) detected in the Brooks database were detected using PLSDB as well. In K. pneumoniae, almost all the PTRUE (29 of 30) detected using the Brooks database were also detected using the PLSDB database. Additionally, we detected 13 shared plasmids between E. coli and K. pneumoniae using the PLSDB database. In contrast, only three common plasmids were observed between E. coli and K. pneumoniae by employing the Brooks database. The Neighbor-Joining phylogenetic tree of IncFIB plasmids that were most prevalent for both E. coli and K. pneumoniae is presented in Figure 1B. The tree is based on MAFFT (Multiple Alignment using Fast Fourier Transform) alignment of plasmids conserved regions [28]. There were four IncFIB sequences of K. pneumoniae. Two K. pneumoniae isolates formed a separate branch, whereas two others were clustered together with IncFIB plasmid sequences from E. coli. An overview of the top 20 most abundant putative plasmids (PTRUE) and replicons retrieved from each database for both E. coli and K. pneumoniae is shown in Table 3.
Regarding plasmid detection in hybrid assembled sequences, we managed to retrieve a higher number of PTRUE from PLSDB than Brooks and PlasmidFinder. All plasmids detected in Brooks for E. coli 39 and K. pneumoniae 23 were also detected in the PLSDB database (Table S2). An overview of the top five PTRUE and replicons retrieved from each database for hybrid assembled isolates can be seen in Table 4.

3.3. Identification of Plasmid-Mediated ARGs

Plasmid assembled files were used to explore the plasmid-mediated ARGs using ResFinder and CARD databases. As can be seen from Table 5 and using plasmid data, regardless of whether identified ARGs located on the putative plasmids (PTRUE) or not, we identified more ARGs using the CARD database in E. coli isolates compared to ResFinder. For K. pneumoniae, opposite results were observed. Moreover, several predicted ARGs were different after annotating the results to plasmid databases (for E. coli isolates, ARGPlasmid-PlasmidFinder > ARGPlasmid-PLSDB > ARGPlasmid-Brooks and for K. pneumoniae isolates, ARGPlasmid-PLSDB > ARGPlasmid-PlasmidFinder > ARGPlasmid-Brooks). Moreover, using hybrid assembled data for E. coli 39 isolates, we detected a higher number of ARGs using the CARD database. The same results for K. pneumoniae 37 were also observed (Table S2). Surprisingly, no ARGs were detected using K. pneumoniae plasmid or hybrid data as ARGPlasmid-Brooks.
As is apparent from Table 5, ARGPlasmid-PlasmidFinder in E. coli plasmid data represent a group with the highest detected number of ARGs. Further details about ARGPlasmid-PlasmidFinder can be seen in Figure 2 and Table 5. For E. coli plasmid data, the majority of ARGPlasmid were found on the IncFII plasmid. Furthermore, plasmids such as Col(pHAD28) and IncI1-1(Gamma) hosted the least ARGs. Some ARGs such as aph(3″)-Ib, aph(6)-Id, blaTEM-1B, and sul2 were carried by more than one type of plasmid.
The most abundant ARGPlasmid genes for plasmid data from E. coli and K. pneumoniae isolates can be seen in Table 6. The majority of detected ARGPlasmid hits in E. coli isolates, carried by PTRUE from PlasmidFinder and PLSDB, were beta-lactamase gene-variants blaTEM-1B and TEM-1B. For E. coli, ARGPlasmid genes such as aph(3″)-Ib, aph(6)-Id, sul1, sul2, and tet(D) were flagged as mutual ARGPlasmid, observed in all databases. For K. pneumoniae, no ARGPlasmid gene was detected on PTRUE from Brooks. However, the qnrS1 gene was found as a mutual ARGPlasmid harbored by PTRUE from both PLSDB and PlasmidFinder.
The ARGplasmid prediction using hybrid assembled sequences is presented in Table 7. Overall, ARGplasmid prediction using hybrid assembled sequences was more consistent between databases compared to plasmid assembled data (Table S2). For instance, in the hybrid assembled E. coli 39 isolate, ARGPlasmid such as aac(3)-VIa and aadA1 were hosted by PTRUE from all databases. For K. pneumoniae isolate 37, predicted ARGPlasmid-PlasmidFinder and ARGPlasmid-PLSDB were entirely matched. Regarding K. pneumoniae isolate 23, besides an extra predicted ARGPlasmid-PlasmidFinder, all the predicted ARGPlasmid-PLSDB were covered by ARGPlasmid-PlasmidFinder.

4. Discussion

In the current research, the applicability of three different plasmid databases and two antibiotics resistance gene databases were assessed using E. coli and K. pneumoniae assemblies taken from Norwegian patients.
We identified a total number of 490 and 52 exclusive putative plasmids using PLSDB and Brooks databases, respectively. Observed differences might be explained by the content of databases, as the method used for developing the databases and the date of last revision (October 2018 for Brooks and November 2020 for PLSDB) as well as their file size (11,677 and 13,789 entries in Brooks and PLSDB, respectively) are different. Although a BLASTn search against Brooks and PLSDB databases resulted in a higher number of putative plasmids than PlasmidFinder, the method has its disadvantages. For instance, using the BLASTn search, we have detected multiple hits with similar lengths, alignment coverage, and percentage identity for the same assigned contig. Therefore, assigning the putative plasmids as PTRUE was challenging. Similar challenges following BLAST+ have been previously described for the FindPlasmid package [29]. Using PlasmidFinder, researchers can directly upload raw files from sequencing platforms. Therefore, de novo assembly is not required, and PlasmidFinder can perform de novo assembly automatically, though the assembly results are not presented by the tool. On the other hand, manual de novo assembly is required in advance to BLASTn search when using other databases such as Brooks and PLSDB. However, one of the disadvantages of using PlasmidFinder is that it currently only covers Enterobacteriaceae and a few Gram-positive bacterial species.
It is clinically relevant to perform downstream analyses such as the prediction of plasmid associated ARGs following plasmid identification. In this study, ResFinder performed better than CARD to predict plasmid associated antibacterial resistance genes (AMRPlasmid-PlasmidFinder, AMRPlasmid-Brooks, and AMRPlasmid-PLSDB) for both plasmid and hybrid assembled data. In a study comparing the performance of resistance gene databases, both CARD and ResFinder performed equally when submitting a single gene sequence, but CARD performed slightly better for assembled data [30]. Although CARD only accepts FASTA assembly files up to 20 Mb, but in addition to acquired gene information, it contains chromosomal mutation data too. However, ResFinder takes raw files, and assembly is not required. Furthermore, in ResFinder, users can choose between acquired genes or chromosomal mutations. One of the ResFinder advantages is flagging the hit with the true circular term, which indicates whether the hit is plasmid associated or not. Therefore, current data suggest using PlasmidFinder and its associated ResFinder online tools as the first choice to predict plasmid associated ARGs.
In the current study, the AMRPlasmid gene profiles differed between E. coli strains carrying plasmids of the same type. Similar results have been reported for Salmonella entrica isolates in Ghana [31]. This further highlights the mobility of genetic elements between plasmids, resulting in acquiring or losing the ability for antimicrobial resistance. IncF plasmids are known carriers of a broad spectrum of antibiotic resistance genes in E. coli [32,33,34]. In line with this, IncFII plasmids were strongly associated with various resistance genes in our study. These plasmids carried TEM-1B, aph, sul, tetA, and dfr genes conferring resistance to penicillins, aminoglycosides, sulfonamides, tetracyclines, and trimethoprim [35]. IncFIB were the most prevalent plasmids in our dataset, and they exhibited a low association with antibiotic resistance genes. As such, aad (aminoglycoside), sul (sulfonamides), and tet (tetracycline) were located on IncFIB plasmid contig in two cases, highlighting the low, albeit growing, antibiotic resistance in Norway [16]. However, phylogenetic analysis of IncFIB plasmids revealed that plasmid sequences were shared between E. coli and K. pneumoniae, probably indicating its ability for inter-species transfer, which raises a concern over rising antibiotic resistance in Norway [17]. Additionally, in this work, we documented the co-existence of blaCTX−M genes with other genes corresponding to resistance against sulfonamide, aminoglycoside, trimethoprim, and tetracycline. This agrees with previous reports indicating that plasmids harboring blaCTX−M genes frequently also carry other genes encoding resistance to other antimicrobials [36,37,38].
The high sequence error rate in Oxford Nanopore Technologies and incongruity between short/fragmented reads from MiSeq Illumina platform and large repetitive regions in plasmids often results in the inaccurate prediction of plasmid-mediated ARGs. To overcome this issue, hybrid assembly has been suggested [12,39]. In the present research, we found that the prediction of ARGPlasmids following hybrid assembly was more consistent across different databases. Having a less fragmented assembly where the circular plasmids are apparent makes the prediction of ARGplasmid more accurate. The current conclusion regarding the applicability of hybrid assembly for plasmid-mediated ARGs detection previously has been made [40]. Therefore, future work implanting hybrid assembly to identify ARGs in bacteria is worth investigating.

5. Conclusions

In conclusion, we have demonstrated that plasmid detection and plasmid-mediated ARG prediction are challenging and to obtain a reliable result, one must consider different tools and databases. In the present study, a combination of PlasmidFinder and ResFinder tools showed promising results for both E. coli and K. pneumoniae isolates. Plasmid detection and prediction of plasmid-mediated ARG can be facilitated using hybrid assembly. Although Norway is considered as a country with a low antibiotic resistance frequency, current research provides a reasonable argument to tackle the slightly increasing antibiotic resistance issue in Norway.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/9/1/52/s1, Table S1: an overview of complete results for plasmid assembled data from different databases. Table S2: an overview of complete results for hybrid assembled data from different databases.

Author Contributions

Conceptualization, R.A., A.K., and E.A.; methodology, A.K., and E.A.; software, A.K., and E.A.; formal analysis, A.K.; investigation, A.K., E.A.; data curation, A.K.; writing—original draft preparation, A.K.; writing—review and editing, R.A., A.K., and E.A.; visualization, A.K., and E.A.; supervision, R.A.; project administration, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Norwegian research council, grant number 273609 to AMR-Diag. The APC was funded by the AMR-Diag grant and from support from the Inland Norway University of Applied Sciences.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality agreement related to the AMR-Diag project.

Acknowledgments

The authors would like to thank Arne Michael Taxt for the selection of the clinical isolates and Stephan Frye for performing the WGS. We would also like to thank Ambreen Kauser for the bioinformatics work she performed on a sub-set of samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cassini, A.; Högberg, L.D.; Plachouras, D.; Quattrocchi, A.; Hoxha, A.; Simonsen, G.S.; Colomb-Cotinat, M.; E Kretzschmar, M.; Devleesschauwer, B.; Cecchini, M.; et al. Attributable deaths and disability-adjusted life-years caused by infections with antibiotic-resistant bacteria in the EU and the European Economic Area in 2015: A population-level modelling analysis. Lancet Infect. Dis. 2019, 19, 56–66. [Google Scholar] [CrossRef] [Green Version]
  2. WHO. Global Priority List of Antibiotic-Resistant Bacteria to Guide Research, Discovery, and Development of New Antibiotics; WHO: Geneva, Switzerland, 2017. [Google Scholar]
  3. Dunn, S.J.; Connor, C.; McNally, A. The evolution and transmission of multi-drug resistant Escherichia coli and Klebsiella pneumoniae: The complexity of clones and plasmids. Curr. Opin. Microbiol. 2019, 51, 51–56. [Google Scholar] [CrossRef] [PubMed]
  4. Von Wintersdorff, C.J.; Penders, J.; Van Niekerk, J.M.; Mills, N.D.; Majumder, S.; Van Alphen, L.B.; Savelkoul, P.H.; Wolffs, P.F. Dissemination of Antimicrobial Resistance in Microbial Ecosystems through Horizontal Gene Transfer. Front. Microbiol. 2016, 7, 173. [Google Scholar] [CrossRef] [Green Version]
  5. Bennett, P.M. Plasmid encoded antibiotic resistance: Acquisition and transfer of antibiotic resistance genes in bacteria. Br. J. Pharmacol. 2008, 153, S347–S357. [Google Scholar] [CrossRef] [Green Version]
  6. Kaper, J.B.; Nataro, J.P.; Mobley, H.L.T. Pathogenic Escherichia coli. Nat. Rev. Microbiol. 2004, 2, 123–140. [Google Scholar] [CrossRef]
  7. Hiraga, S.; Sugiyama, T.; Itoh, T. Comparative analysis of the replicon regions of eleven ColE2-related plasmids. J. Bacteriol. 1994, 176, 7233–7243. [Google Scholar] [CrossRef] [Green Version]
  8. Köser, C.U.; Ellington, M.J.; Cartwright, E.J.P.; Gillespie, S.H.; Brown, N.M.; Farrington, M.; Holden, M.T.G.; Dougan, G.; Bentley, S.D.; Parkhill, J.; et al. Routine use of microbial whole genome sequencing in diagnostic and public health microbiology. PLoS Pathog. 2012, 8, e1002824. [Google Scholar] [CrossRef] [Green Version]
  9. Arredondo-Alonso, S.; Rogers, M.R.C.; Braat, J.C.; Verschuuren, T.D.; Top, J.; Corander, J.; Willems, R.J.L.; Schürch, A.C. Mlplasmids: A user-friendly tool to predict plasmid- and chromosome-derived sequences for single species. Microbial. Genom. 2018, 4, e000224. [Google Scholar] [CrossRef] [Green Version]
  10. Tolmasky, M.E.; Alonso, J.C. Plasmids: Biology and Impact in Biotechnology and Discovery; ASM Press: Washington, DC, USA, 2015. [Google Scholar]
  11. Lal Gupta, C.; Kumar Tiwari, R.; Cytryn, E. Platforms for elucidating antibiotic resistance in single genomes and complex metagenomes. Environ. Int. 2020, 138, 105667. [Google Scholar] [CrossRef]
  12. Taxt, A.M.; Avershina, E.; Frye, S.A.; Naseer, U.; Ahmad, R. Rapid identification of pathogens, antibiotic resistance genes and plasmids in blood cultures by nanopore sequencing. Sci. Rep. 2020, 10, 7622. [Google Scholar] [CrossRef]
  13. Harstad, H.; Ahmad, R.; Bredberg, A. Nanopore-based DNA sequencing in clinical microbiology: Preliminary assessment of basic requirements. arXiv 2018, arXiv:bio/382580. [Google Scholar]
  14. Yang, Y.; Stathis, D.; Sharma, P.; Paul, K.; Hemani, A.; Grabherr, M.; Ahmad, R. RiBoSOM: Rapid bacterial genome identification using self-organizing map implemented on the synchoros SiLago platform. In ACM International Conference Proceeding Series; Association for Computing Machinery (ACM): Pythagorion, Greece, 2018; pp. 105–114. [Google Scholar]
  15. Stathis, D.; Yang, Y.; Tewari, S.; Hemani, A.; Paul, K.; Grabherr, M.; Ahmad, R. Approximate Computing Applied to Bacterial Genome Identification using Self-Organizing Maps. In Proceedings of the 2019 IEEE Computer Society Annual Symposium on VLSI (ISVLSI), Miami, FL, USA, 15–17 July 2019; pp. 560–567. [Google Scholar] [CrossRef]
  16. Klein, E.Y.; Tseng, K.K.; Pant, S.; Laxminarayan, R. Tracking global trends in the effectiveness of antibiotic therapy using the Drug Resistance Index. BMJ Glob. Health 2019, 4, e001315. [Google Scholar] [CrossRef] [PubMed]
  17. NORM/NORM-VET. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway; Norwegian Veterinary Institute: Tromsø/Oslo, Norway, 2019. [Google Scholar]
  18. Avershina, E. AMR-Diag: Neural Network based Genotype-to-Phenotype Prediction of Resistance towards Beta-Lactams and Colistin in Escherichia Coli and Klebsiella Pneumoniae. Comput. Struct. Biotechnol. J. Elsevier 2020. under review. [Google Scholar]
  19. Andrews, S. FastQC: A Quality Control. Tool for High. Throughput Sequence Data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 24 September 2020).
  20. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  21. Antipov, D.; Hartwick, N.; Shen, M.; Raiko, M.; Lapidus, A.; Pevzner, P.A. PlasmidSPAdes: Assembling plasmids from whole genome sequencing data. Bioinformatics 2016, 32, 3380–3387. [Google Scholar] [CrossRef] [Green Version]
  22. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  23. Carattoli, A.; Zankari, E.; García-Fernández, A.; Larsen, M.V.; Lund, O.; Villa, L.; Aarestrup, F.M.; Hasman, H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 2014, 58, 3895–3903. [Google Scholar] [CrossRef] [Green Version]
  24. Galata, V.; Fehlmann, T.; Backes, C.; Keller, A. PLSDB: A resource of complete bacterial plasmids. Nucleic Acids Res. 2019, 47, D195–D202. [Google Scholar] [CrossRef]
  25. Brooks, L.; MoSistrom, M. A Curated, Comprehensive Database of Plasmid Sequences. Microbiol. Resour. Announc. 2019, 8, e0132518. [Google Scholar] [CrossRef] [Green Version]
  26. Bortolaia, V.; Ruppe, R.S.; Roberts, E.; Schwarz, M.C.; Cattoir, S.; Philippon, V.; Allesoe, A.; Rebelo, R.L.; Florensa, A.R.; Fagelhauer, A.F.; et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 2020, 75, 3491–3500. [Google Scholar] [CrossRef]
  27. Alcock, B.P.R.; Lau, A.R.; Tsang, T.T.Y.; Bouchard, K.K.; Edalatmand, M.; Huynh, A.; Nguyen, W.; Cheng, A.V.; Liu, A.A.; Min, S.; et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020, 48, D517–D525. [Google Scholar] [CrossRef] [PubMed]
  28. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2017, 20, 1160–1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Hall, B.G. Identifying Plasmids in Bacterial Genome Assemblies. arXiv 2018, arXiv:bio/332049. [Google Scholar]
  30. Xavier, B.B.; Das, A.J.; Cochrane, G.; De Ganck, S.; Kumar-Singh, S.; Aarestrup, F.M.; Goossens, H.; Malhotra-Kumar, S. Consolidating and Exploring Antibiotic Resistance Gene Data Resources. J. Clin. Microbiol. 2016, 54, 851–859. [Google Scholar] [CrossRef] [Green Version]
  31. Kudirkiene, E.; Andoh, L.A.; Ahmed, S.; Herrero-Fresno, A.; Dalsgaard, A.; Obiri-Danso, K.; Olsen, J.E. The Use of a Combined Bioinformatics Approach to Locate Antibiotic Resistance Genes on Plasmids from Whole Genome Sequences of Salmonella enterica Serovars from Humans in Ghana. Front. Microbiol. 2018, 9, 1010. [Google Scholar] [CrossRef]
  32. Lyimo, B.; Buza, J.; Subbiah, M.; Temba, S.; Kipasika, H.; Smith, W.; Call, D.R. IncF Plasmids Are Commonly Carried by Antibiotic Resistant Escherichia coli Isolated from Drinking Water Sources in Northern Tanzania. Int. J. Microbiol. 2016, 2016, 3103672. [Google Scholar] [CrossRef] [Green Version]
  33. Torres-González, P.; Bobadilla-del Valle, M.; Tovar-Calderon, E.; Leal-Vega, F.; Hernandez-Cruz, A.; Martinez-Gamboa, A.; Niembro-Ortega, M.D.; Martinez-Gamboa, A.; Ponce-de-Leon, A. Outbreak Caused by Enterobacteriaceae Harboring NDM-1 Metallo-β-Lactamase Carried in an IncFII Plasmid in a Tertiary Care Hospital in Mexico City. Antimicrob. Agents Chemother. 2015, 59, 7080–7083. [Google Scholar] [CrossRef] [Green Version]
  34. Deng, Y.; Zeng, Z.; Chen, S.; He, L.; Liu, Y.; Wu, C.; Chen, Z.; Yao, Q.; Hou, J.; Yang, T.; et al. Dissemination of IncFII plasmids carrying rmtB and qepA in Escherichia coli from pigs, farm workers and the environment. Clin. Microbiol. Infect. 2011, 17, 1740–1745. [Google Scholar] [CrossRef] [Green Version]
  35. Van Hoek, A.H.; Mevius, D.; Guerra, B.; Mullany, P.; Roberts, A.P.; Aarts, H.J. Acquired Antibiotic Resistance Genes: An. Overview. Front. Microbiol. 2011, 2, 203. [Google Scholar]
  36. Liu, B.T.; Yang, Q.E.; Li, L.; Sun, J.; Liao, X.P.; Fang, L.X.; Yang, S.S.; Deng, H.; Liu, Y.H. Dissemination and characterization of plasmids carrying oqxAB-bla CTX-M genes in Escherichia coli isolates from food-producing animals. PLoS ONE 2013, 8, e73947. [Google Scholar]
  37. Sun, H.; Li, S.; Xie, Z.; Yang, F.; Zhu, Y.; Zhao, X.; Jiang, S. A novel multidrug resistance plasmid isolated from an Escherichia coli strain resistant to aminoglycosides. J. Antimicrob. Chemother. 2012, 67, 1635–1638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Huang, W.; Wang, G.; Sebra, R.; Zhuge, J.; Yin, C.; Aguero-Rosenfeld, M.E.; Schuetz, A.N.; Dimitrova, N.; Fallon, J.T. Emergence and Evolution of Multidrug-Resistant Klebsiella pneumoniae with both bla(KPC) and bla(CTX-M) Integrated in the Chromosome. Antimicrob. Agents Chemother. 2017, 61, e00076-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bertrand, D.; Shaw, J.; Kalathiyappan, M.; Ng, A.H.Q.; Kumar, M.S.; Li, C.; Dvornicic, M.; Soldo, J.P.; Koh, J.Y.; Tong, C.; et al. Hybrid metagenomic assembly enables high-resolution analysis of resistance determinants and mobile elements in human microbiomes. Nat. Biotechnol. 2019, 37, 937–944. [Google Scholar] [CrossRef] [PubMed]
  40. Berbers, B.; Ceyssens, P.-J.; Bogaerts, P.; Vanneste, K.; Roosens, N.H.; Marchal, K.; De Keersmaecker, S.C.J. Development of an NGS-Based Workflow for Improved Monitoring of Circulating Plasmids in Support of Risk Assessment of Antimicrobial Resistance Gene Dissemination. Antibiotics 2020, 9, 503. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) Number of uniquely identified putative plasmids (PTRUE) using different databases after removing duplicates hits per isolate, for both E. coli and K. pneumoniae. (B) phylogenic tree of IncFIB plasmids as the most prevalent plasmids for both E. coli and K. pneumoniae.
Figure 1. (A) Number of uniquely identified putative plasmids (PTRUE) using different databases after removing duplicates hits per isolate, for both E. coli and K. pneumoniae. (B) phylogenic tree of IncFIB plasmids as the most prevalent plasmids for both E. coli and K. pneumoniae.
Microorganisms 09 00052 g001
Figure 2. The co-existence of ARGPlasmid genes and different plasmids detected by PlasmidFinder in plasmid data from E. coli isolates. Numbers inside each cell indicating the number of isolates where ARGs were found on the corresponding plasmid.
Figure 2. The co-existence of ARGPlasmid genes and different plasmids detected by PlasmidFinder in plasmid data from E. coli isolates. Numbers inside each cell indicating the number of isolates where ARGs were found on the corresponding plasmid.
Microorganisms 09 00052 g002
Table 1. An overview of the number of clinical isolates collected from Norwegian patients.
Table 1. An overview of the number of clinical isolates collected from Norwegian patients.
IsolatesSource
BloodUrine
E. coli535
K. pneumoniae404
Table 2. An overview of general statistics (mean ± SD) obtained using the QUAST tool for the Scheme 39. and K. pneumoniae 23, 37).
Table 2. An overview of general statistics (mean ± SD) obtained using the QUAST tool for the Scheme 39. and K. pneumoniae 23, 37).
Number of ContigsLargest Contig (bp)Total Length (bp)GC (%)N50
E. coli87.5 ± 484.847,992 ± 65,592271,880 ± 634,15548.8 ± 2.331,029 ± 62,333
K. pneumoniae7.1 ± 7.935,135 ± 48,18498,095 ± 170,35448.5 ± 3.823,022 ± 29,064
E. coli 39602,757,7345,955,16350.51996,338
K. pneumoniae 23155,305,1065,831,97656.755,305,106
K. pneumoniae 37282,952,4495,558,21357.592,952,449
Table 3. An overview of the top 20 most abundant putative plasmids (PTRUE) in both E. coli and K. pneumoniae retrieved from different databases or tools. (n: number of isolates).
Table 3. An overview of the top 20 most abundant putative plasmids (PTRUE) in both E. coli and K. pneumoniae retrieved from different databases or tools. (n: number of isolates).
PlasmidFinderBrooksPLSDB
Replicon NamenPlasmid NamenPlasmid Namen
E. coliIncFIB (AP001918)23pECO-bc65pMTY17816_OXA487
Col15615pUT1894pMTY17823_OXA487
IncFIA9pUM1464p53_E-OXA487
IncFII (29)9pRS2184pECO-bc65
IncFII8pSF-166-14unnamed 75
Col (BS512)4pSaT0404pUTI894
IncFII (pRSB107)4pZH063-14pEC14_1144
Col (MG828)3pECOS883pUM1464
Col82822pMRY16-002_53pRS2184
IncFIC (FII)2pAPEC-O78-ColV3pSF-166-14
IncFII (pCoo)2pPSUO78_13unnamed 14
IncI1-I (Gamma)2pSF-088-13pSaT0404
IncX42pG749_13pZH063_14
Col (pHAD28)1pECAZ147_13p0.1229_14
FIA (pBK30683)1pECSF12plasmid 24
IncFIB (pB171)1pCE10A2unnamed3
IncM11p3PCN0332pF2_18C_Col3
IncQ11pMVAST0167_12pDB42773
IncX11pKPN-7c32pCERC43
p01111pEC732_62pCERC53
K. pneumoniaeIncFIB (K)4pKpn23412-41pEC-2431
IncFIA (HI1)2pK20441pM206-OXA1811
IncFIB (pKPHS1)2pKCTC22421pM513-OXA1811
IncHI1B (pNDM-MAR)2PCN061p41pM517-OXA1811
repB2plasmid B1pM518-OXA1811
Col (KPHS6)1plasmid A1pTHC11-11
Col440I1pRJF2931pTMTA636311
Col440II1pRJF9991pCA081
ColpVC1pNY9_31pKpvK541
IncFIB (AP001918)1pESBL9311pBH100 alpha1
IncFII1unnamed 11pVir_0951321
IncFII (29)1unnamed1pSF072011
IncFII (K)1unnamed 31pYHCC1
IncX31pKp_Goe_579-61pKS221
pKp_Goe_473-51U25P0021
pKp_Goe_832-51pFAM223211
pKp_Goe_024-51p1836601
pKp_Goe_026-51pKP2442_7c3311
pKp_Goe_021-51pOXA-181_291441
pKp_Goe_304-51pEC-NRS181
Table 4. An overview of the first top-five putative plasmids (PTRUE) in hybrid assembled E. coli and K. pneumoniae isolates retrieved from the different databases.
Table 4. An overview of the first top-five putative plasmids (PTRUE) in hybrid assembled E. coli and K. pneumoniae isolates retrieved from the different databases.
Replicon Name from PlasmidFinderPlasmid Name from BrooksPlasmid Name from PLSDB
E. coli 39IncHI2pCFSAN002069_01pCFSAN002069_01
IncHI2Ap25155p25155
p0111PDM02p2EC1-4
PDM04pSH1148_107
PDM04
K. pneumoniae 23IncFIB(pNDM-Mar)pKp_Goe_579-5p5
IncFIB(pQil)pKp_Goe_832-5pEC25-4
IncFII(K)pKp_Goe_304-5p4_VBA2172
pKp_Goe_021-5pRGI01215
pKp_Goe_026-5p1_040074
K. pneumoniae 37IncFII U25P002
IncFII(pKP91)
IncFIB(K)
IncFIA(HI1)
Table 5. Number of total identified ARGs and ARGPlasmid (sharing the same contig with PTRUE from three different databases), using ResFinder and CARD databases in both E. coli and K. pneumoniae and hybrid assembled isolates. The number inside parentheses shows the percentage of ARGPlasmid identified in PTRUE from different plasmid databases. E. coli 39, K. pneumoniae 23, and K. pneumoniae 37 are hybrid assembled isolates.
Table 5. Number of total identified ARGs and ARGPlasmid (sharing the same contig with PTRUE from three different databases), using ResFinder and CARD databases in both E. coli and K. pneumoniae and hybrid assembled isolates. The number inside parentheses shows the percentage of ARGPlasmid identified in PTRUE from different plasmid databases. E. coli 39, K. pneumoniae 23, and K. pneumoniae 37 are hybrid assembled isolates.
E. coliK. pneumoniaeE. coli 39K. pneumoniae 23K. pneumoniae 37
Res FinderCARDRes FinderCARDRes FinderCARDRes FinderCARDRes FinderCARD
Total number of ARGs6110029159542424917
PlasmidFinder31 (51%)27 (27%)1 (3%)1 (6%)6 (54%)6 (11%)14 (58%)12 (50%)6 (54%)5 (29%)
Brooks10 (16%)10 (10%)--2 (22%)1 (2%)----
PLSDB18 (29%)16 (16%)3 (10%)2 (13%)7 (77%)7 (13%)7 (29%)6 (25%)6 (54%)5 (29%)
Table 6. Gene name and number of isolates with most abundant ARGPlasmid for both E. coli and K. pneumoniae isolates.
Table 6. Gene name and number of isolates with most abundant ARGPlasmid for both E. coli and K. pneumoniae isolates.
E. coliK. pneumoniae
ResFinderCARDResFinderCARD
Gene NamenGene NamenGene NamenGene Namen
PlasmidFinderaadA12aadA1qnrS11qnrS11
aph(3″)-Ib4aph(3″)-Ib3
aph(6)-Id4aph(6)-Id4
blaTEM-1A1TEM-18
blaTEM-1B7TEM-401
blaTEM-1C1
dfrA82dfrA82
qnrB191qnrB191
sul11sul11
sul24sul24
tet(A)3qacEdelta11
tet(D)1tet(D)1
BrooksaadA11aadA1
aph(3″)-Ib2aph(3″)-Ib1
aph(6)-Id2aph(6)-Id2
blaTEM-1B1TEM-11
sul11sul11
sul22sul22
tet(D)1tet(D)1
qacEdelta11
PLSDBaadA11aadA1qnrS11qnrS11
aadA51aadA51floR1SHV-11
aph(3″)-Ib2aph(3″)-Ib1blaSHV-991
aph(6)-Id2aph(6)-Id2
blaTEM-1B3TEM-13
blaTEM-1D1
dfrA171
mph(A)1mph(A)1
sul12sul12
sul22sul22
tet(A)1qacEdelta12
tet(D)1tet(D)1
Table 7. Predicted ARGplasmid using hybrid assembled data for both E. coli and K. pneumoniae isolates.
Table 7. Predicted ARGplasmid using hybrid assembled data for both E. coli and K. pneumoniae isolates.
E. coli 39K. pneumoniae 23K. pneumoniae 37
ResFinderCARDResFinderCARDResFinderCARD
PlasmidFinderaac(3)-VIa aac(3)-IIdaac(3)-IIdaac(3)-IIaaac(3)-IIe
aadA1aadAaac(6′)-Ibaac(6′)-Ib10aac(6′)-Ib3
aac(6′)-Ib-cr aac(6′)-Ib-cr
aadA1aadAblaCTX-M-14CTX-M-14
aph(3′)-Iaaph(3′)-IablaTEM-1BTEM-1
aph(3″)-Ibaph(6)-IdcmlA1cmlA5
aph(6)-Id qacEdelta1
blaCTX-M-15CTX-M-15
blaOXA-9OXA-9
blaSHV-12SHV-134
blaTEM-1ATEM-1
blaTEM-1B
catA1catI
dfrA30qacE
sul2sul2
Brooksaac(3)-VIa
aadA1
PLSDBaac(3)-VIa aac(6′)-Ibaac(6′)-Ib10aac(3)-IIaaac(3)-IIe
aadA1aadAaac(6′)-Ib-cr aac(6′)-Ib3
blaCTX-M-2CTX-M-2aadA1aadAaac(6′)-Ib-cr
blaTEM-1BTEM-1blaCTX-M-15CTX-M-15blaCTX-M-14CTX-M-14
dfrA1dfrA1blaOXA-9OXA-9blaTEM-1BTEM-1
sul1sul1blaSHV-12SHV-134cmlA1cmlA5
tet(A)qacEdelta1blaTEM-1ATEM-1 qacEdelta1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khezri, A.; Avershina, E.; Ahmad, R. Plasmid Identification and Plasmid-Mediated Antimicrobial Gene Detection in Norwegian Isolates. Microorganisms 2021, 9, 52. https://doi.org/10.3390/microorganisms9010052

AMA Style

Khezri A, Avershina E, Ahmad R. Plasmid Identification and Plasmid-Mediated Antimicrobial Gene Detection in Norwegian Isolates. Microorganisms. 2021; 9(1):52. https://doi.org/10.3390/microorganisms9010052

Chicago/Turabian Style

Khezri, Abdolrahman, Ekaterina Avershina, and Rafi Ahmad. 2021. "Plasmid Identification and Plasmid-Mediated Antimicrobial Gene Detection in Norwegian Isolates" Microorganisms 9, no. 1: 52. https://doi.org/10.3390/microorganisms9010052

APA Style

Khezri, A., Avershina, E., & Ahmad, R. (2021). Plasmid Identification and Plasmid-Mediated Antimicrobial Gene Detection in Norwegian Isolates. Microorganisms, 9(1), 52. https://doi.org/10.3390/microorganisms9010052

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop