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

Whole-genome analyses of extended-spectrum or AmpC β-lactamase-producing Escherichia coli isolates from companion dogs in Japan

  • Mayo Yasugi ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing

    shishimaru@vet.osakafu-u.ac.jp

    Affiliations Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan, Asian Health Science Research Institiute, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Shingo Hatoya,

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

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Daisuke Motooka,

    Roles Data curation, Formal analysis, Funding acquisition, Writing – original draft, Writing – review & editing

    Affiliation Research Institute for Microbial Diseases, Osaka University, Suita, Osaka, Japan

  • Yuki Matsumoto,

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

    Affiliation Research Institute for Microbial Diseases, Osaka University, Suita, Osaka, Japan

  • Shunsuke Shimamura,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Hiroyuki Tani,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Masaru Furuya,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Keiichiro Mie,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Masami Miyake,

    Roles Writing – review & editing

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan

  • Shota Nakamura,

    Roles Data curation, Formal analysis, Funding acquisition, Writing – original draft, Writing – review & editing

    Affiliation Research Institute for Microbial Diseases, Osaka University, Suita, Osaka, Japan

  • Terumasa Shimada

    Roles Funding acquisition, Resources, Supervision, Writing – review & editing

    Affiliations Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Izumisano, Osaka, Japan, Asian Health Science Research Institiute, Osaka Prefecture University, Izumisano, Osaka, Japan

Abstract

The emergence and global spread of extended-spectrum or AmpC β-lactamase (ESBL/AmpC)-producing Enterobacteriaceae in companion animals have led to the hypothesis that companion animals might be reservoirs for cross-species transmission because of their close contact with humans. However, current knowledge in this field is limited; therefore, the role of companion animals in cross-species transmission remains to be elucidated. Herein, we studied ESBL/AmpC-producing Enterobacteriaceae, Escherichia coli in particular, isolated from extraintestinal sites and feces of companion dogs. Whole-genome sequencing analysis revealed that (i) extraintestinal E. coli isolates were most closely related to those isolated from feces from the same dog, (ii) chromosomal sequences in the ST131/C1-M27 clade isolated from companion dogs were highly similar to those in the ST131/C1-M27 clade of human origin, (iii) certain plasmids, such as IncFII/pMLST F1:A2:B20/blaCTX-M-27, IncI1/pMLST16/blaCTX-M-15, or IncI1/blaCMY-2 from dog-derived E. coli isolates, shared high homology with those from several human-derived Enterobacteriaceae, (iv) chromosomal blaCTX-M-14 was identified in the ST38 isolate from a companion dog, and (v) eight out of 14 tested ESBL/AmpC-producing E. coli isolates (i.e., ST131, ST68, ST405, and ST998) belonged to the human extraintestinal pathogenic E. coli (ExPEC) group. All of the bla-coding plasmids that were sequenced genome-wide were capable of horizontal transfer. These results suggest that companion dogs can spread ESBL/AmpC-producing ExPEC via their feces. Furthermore, at least some ESBL/AmpC-producing ExPECs and bla-coding plasmids can be transmitted between humans and companion dogs. Thus, companion dogs can act as an important reservoir for ESBL/AmpC-producing E. coli in the community.

Introduction

The World Health Organization published a report in 2014 showing very high rates of resistance in bacteria that cause common health-care associated and community-acquired infections worldwide [1]. In 2016, it was estimated that unless action was taken, the number of deaths resulting from antimicrobial resistance (AMR) could balloon to ten million per year by 2050 [2]. AMR bacteria are now one of the most serious issues worldwide. Of these, third-generation cephalosporin-resistant (3GCR) Enterobacteriaceae are of great concern because of an increasing trend in their prevalence globally [3]. 3GCR Enterobacteriaceae produce extended-spectrum or AmpC β-lactamases (ESBL/AmpC) that hydrolyze the β-lactam ring of third-generation cephalosporins [4]. ESBL-producing Enterobacteriaceae emerged during the 1980s. Most of the early ESBLs described were of the TEM or SHV type, mainly found in Klebsiella pneumoniae, and were usually associated with nosocomial outbreaks [5]. Since the late 1990s, the incidence of ESBLs of the CTX-M type, particularly in Escherichia coli, has increased dramatically, and not only in nosocomial environments, but also in community settings [6,7]. The CTX-M type is currently the most prevalent ESBL worldwide. AmpC genes are located on the chromosome in several Enterobacteriaceae, such as Enterobacter spp. and Citrobacter spp [8]. Chromosomal AmpC genes were mobilized onto plasmids by insertion sequences, and plasmid-borne AmpC emerged in K. pneumoniae in the late 1980s [9]. Transmissible plasmids encoding AmpC gene have spread in bacteria that lack or poorly express chromosomal AmpC genes, such as E. coli, K. pneumoniae, and Proteus mirabilis [8,9]; therefore, plasmid-borne AmpC-producing Enterobacteriaceae are now highly prevalent in humans.

Global concerns have also been raised about the dissemination of ESBL/AmpC-producing Enterobacteriaceae in animals. In livestock animals, poultry and poultry products are considered reservoirs of ESBL/AmpC-producing Salmonella enterica and E. coli because of the high incidence of CMY-2-producing E. coli in poultry [10]. Many studies suggested that ESBL-producing E. coli are clonally transferred from livestock to humans, based on typing methods [1113]. More recently, whole-genome analysis of ESBL/AmpC-producing E. coli isolated from humans, livestock animals, and meat products suggested that clonal transmission of the bla gene is a rare event, while the horizontal transfer of temporally stable bla-carrying plasmids is more likely the dominant form of cross-species transmission [11,14]. ESBL/AmpC-producing Enterobacteriaceae are also detected worldwide in companion animals [1517]. Surveillance in Japan showed that approximately 40% of tested diseased dogs and cats carried 3GCR E. coli [18,19]. The emergence and spread of ESBL/AmpC-producing Enterobacteriaceae among companion animals pose a risk of the animals being reservoirs for cross-species transmission, because of their close contact with humans [20,21]. Clonal transfer of ESBL-producing Enterobacteriaceae between companion animals and humans in the same household was previously suggested by typing methods [17,22], and more recently by whole-genome analysis [2325]. However, current knowledge of cross-species transmission is still limited; thus, the precise role of companion animals in cross-species transmission of ESBL/AmpC-producing Enterobacteriaceae is not fully understood. To shed light on this issue, here we studied ESBL/AmpC-producing Enterobacteriaceae, in particular E. coli, isolated from companion dogs and compared the chromosomal and plasmid sequences with those from humans at the genome-wide scale.

Materials and methods

Ethics

Animal samples were collected using protocols approved by the Institutional Review Boards of the Veterinary Medical Center, Osaka Prefecture University. Written informed consent was obtained from the owners. All samples were anonymized before use.

Bacterial isolates

All clinical specimens characterized in this study were collected at the Veterinary Medical Center, Osaka Prefecture University in 2018. All extraintestinal bacterial isolates were isolated, identified, and tested for antimicrobial susceptibility as part of a clinical diagnosis by Japan Clinical Laboratories (Kyoto, Japan). The isolates that were suspected to be ESBL/AmpC-producing Enterobacteriaceae were provided by Japan Clinical Laboratories, cultured in brain heart infusion (BHI) broth at 37°C overnight, and stored at −80°C with glycerol (final concentration, 20%) for future use. Fecal samples were obtained from 12 dogs. The feces were streaked on MacConkey agar containing 4 μg/ml cefotaxime, and cultured at 37°C overnight. After colonies were subcloned, five colonies underwent further biochemical characterization for species identification. Finally, the identified Enterobacteriaceae isolates were cultured in BHI broth supplemented with 4 μg/ml cefotaxime, and stored at −80°C with glycerol for future use.

Antimicrobial susceptibility testing

Minimal inhibitory concentration (MIC) assays for ampicillin, cefazolin, cefmetazole, cefotaxime, ceftazidime, ceftriaxone, cefpirome, imipenem, meropenem, levofloxacin, and ciprofloxacin were performed on the isolates according to the broth microdilution method of the Clinical and Laboratory Standards Institute (CLSI) by following instructions in CLSI document M07-A10 [26]. The results were interpreted using CLSI M100-S26 guidelines [27]. The 3GCR isolates underwent further phenotypic testing for the production of ESBL, AmpC, and/or metallo-β-lactamase as follows: ESBL production was confirmed by a disc diffusion assay using a mixture of antibiotics (cefotaxime or ceftazidime) and clavulanic acid; AmpC production was confirmed by a disc diffusion assay using a mixture of antibiotics (cefmetazole, cefotaxime, or ceftazidime) and aminophenylboronic acid; metallo-β-lactamase production was confirmed by a double-disc synergy assay using antibiotics (imipenem or ceftazidime) and sodium mercaptoacetate. The results were interpreted as described previously [2729].

Genotypes of the ESBL/AmpC-producing Enterobacteriaceae were determined by a direct sequencing method. The isolates were cultured in tryptic soy broth (TSB) supplemented with 4 μg/ml cefotaxime at 37°C overnight. After centrifugation, culture pellets were suspended in distilled water and boiled for 10 min. After another centrifugation step, supernatants were subjected to PCR using GoTaq polymerase (Promega, Madison, WI) with genogroup-specific primers (Table 1). After genogroups were identified, the same samples were subjected to PCR using PrimeSTAR GXL DNA polymerase (Takara, Shiga, Japan) with the genotype-specific primers described in Table 1. After electrophoresis, amplified DNA was extracted with the QIAquick Gel Extraction kit (Qiagen, Hilden, Germany) and sequenced (Macrogen, Kyoto, Japan). Blastp (https://blast.ncbi.nlm.nih.gov/Blast.cgi) was used to find sequence similarities, and to determine the CTX-M and CMY-2 genotype. Single-nucleotide polymorphisms (SNPs) were used to determine the TEM genotype [30,31].

Multilocus sequencing typing (MLST)

The extracted DNA described above under ‘Antimicrobial susceptibility testing’ was amplified by PCR using PrimeSTAR GXL DNA polymerase (Takara) with gene-specific primers (Table 1), and sequenced (Macrogen). MLST analyses for E. coli isolates was performed according to the protocol described in the MLST database for E. coli (https://enterobase.warwick.ac.uk/species/ecoli/allele_st_search) [40].

Pulsed-field gel electrophoresis (PFGE)

PFGE was performed as described previously [41], with slight modifications. Briefly, the isolates were cultured in TSB supplemented with 1 μg/ml cefotaxime at 37°C overnight. After centrifugation, PFGE plugs were prepared by mixing culture pellets with 1% SeaKem Gold Agarose (Lonza, Basel, Switzerland). The plugs were digested with 30 U XbaI at 37°C overnight. The plugs were electrophoresed with a CHEF-DR III system (Biorad, Hercules, CA) in a 1% SeaKem Gold Agarose gel in 0.5× TBE buffer at 14°C and 6 V/cm for 16 h. Switching times were ramped from 2.2 to 54.2 s. Dendrographic analysis of the DNA fragments was performed using FPQuest software (Bio-rad). The Salmonella serotype Braenderup H9812 strain provided by Dr. Umeda, Osaka Institute of Public Health, was used as a DNA size standard.

Whole-genome sequencing (WGS) and data analysis

E. coli isolates were subjected to WGS using the MinION device (Oxford Nanopore Technologies, Oxford, UK) and the MiSeq system (Illumina, San Diego, CA) as described previously [42]. The bacterial isolates were cultured overnight in BHI broth supplemented with 1.28 μg/ml cefotaxime. Genomic DNA was extracted from the cell culture using the DNeasy PowerSoil kit (Qiagen). To prepare MinION sequencing libraries, 1.5–2.5 μg of DNA was sheared using a g-TUBE (Covaris, Woburn, MA) to obtain 8000-bp fragments. One microgram of the sheared DNA was converted to each library using the Ligation Sequencing kit (SQK-LSK109; Oxford Nanopore Technologies) and the Native Barcoding kit (EXP-NBD104; Oxford Nanopore Technologies). To prepare the MiSeq sequencing libraries, Illumina libraries was constructed using the Nextera DNA Flex Library Prep kit (Illumina). Paired-end sequencing was carried out on the MiSeq instrument (Illumina) with the MiSeq Reagent kit version 2 (500 cycles) (Illumina). To construct genomes, a hybrid assembly was conducted using the Unicycler assembler version 0.4.4 [43] and the Flye assembler version 2.6 [44], with sequence reads obtained from the MinION device and the MiSeq instrument. For Flye assembly, assembled sequences were polished using BWA version 0.7.17-r1188 [45] and Pilon version 1.23 [46], with MiSeq sequence reads at three times. Annotation was performed using both the RAST server [47] and the DFAST server [48], followed by manual curation using blastp (https://blast.ncbi.nlm.nih.gov/Blast.cgi). BacWGSTdb 2.0 [49] was used to track the source of whole-genome sequenced isolates. Plasmid replicon typing, identification of resistance and virulence genes and O, H, and FimH typing were carried out using PlasmidFinder 1.3 [50], pMLST 2.0 [50], ResFinder 2.1 [51], VirulenceFinder 2.0 [52], SerotypeFinder [53], and FimTyper [54], respectively. Genetic structures were compared in Easyfig [55].

SNP calling and phylogenetic analysis

Mobile genetic elements (miniature inverted repeats, insertion sequences, composite transposons, pseudo-composite transposons, unit transposons, integrative conjugative elements, integrative mobilizable elements, and cis-mobilizable elements), prophage, and plasmid regions of the genomic sequences of the reference isolate (described below under each result) were predicted by MGEFinder 1.0.2 [56], PHASTER [57], and PlasmidFinder 1.3, respectively, followed by removal of these regions from the sequences. CSI Phylogeny 1.4 [58] was used to call SNPs and to construct a maximum likelihood phylogenetic tree of E. coli isolates inferred from the concatenated alignment of detected SNPs. CSI Phylogeny pipeline was run with default parameters: a minimal depth at SNP positions of 10 reads, a minimal relative depth at SNP positions of 10%, a minimal distance between SNPs of 10bp, a minimal Z-score of 1.96, a minimal SNP quality of 30, and a minimal read mapping quality of 25 [59]. The sequence data of an isolate belonging to a different ST were used as the outgroup.

Conjugal transfer of plasmids

Bacterial conjugation was carried out as described previously [60] with slight modifications. The recipient was a rifampicin-resistant E. coli K-12 DH5α strain [61] provided by Dr. Umeda, Osaka Institute of Public Health. Overnight cultures of E. coli isolates in TSB supplemented with 4 μg/ml cefotaxime were mixed with the cultures of the recipient in TSB supplemented with 50 μg/ml rifampicin at a 1:10 ratio. The bacterial mixture was then cultured in tenfold dilutions of TSB overnight without agitation. Transconjugants were selected on TSB agar supplemented with 4 μg/ml cefotaxime and 50 μg/ml rifampicin. The transformation efficiency was calculated as the ratio of the number of transconjugant to recipient colonies. Transfer of the target plasmid was confirmed by testing antimicrobial susceptibility, and by detecting resistance genes in the transconjugants.

Accession numbers

WGS data were deposited at DDBJ/EMBL/GenBank under the accession numbers BNIP00000000, BNIQ00000000, BNIR00000000, BNIS00000000, BNIT00000000, BNIU00000000, BNIV00000000, BNIW00000000, BNIX00000000, BNIY00000000, BNIZ00000000, BNJA00000000, BNJB00000000, and BNJC00000000.

Results and discussion

Prevalence of 3GCR Enterobacteriaceae in companion dogs

The prevalence of concerning AMR Enterobacteriaceae, from a human public-health perspective (carbapenem-resistant, 3GCR, and fluoroquinolone-resistant (FQR) Enterobacteriaceae) [18,19,62], was investigated in 188 extraintestinal bacteria isolated from 138 companion animals with suspected bacterial infectious disease at the Veterinary Medical Center of Osaka Prefecture University. Three, seven, and nine of the 65 Enterobacteriaceae isolates were suspected to be 3GCR, FQR, and both 3GCR and FQR, respectively, by antimicrobial susceptibility tests performed by the outsourcing service. Carbapenem-resistant Enterobacteriaceae were not detected in this study. These results agree with surveillance reports [18,19], which describe the prevalence of 3GCR and FQR Enterobacteriaceae in companion dogs and cats in Japan. The 12 Enterobacteriaceae isolates suspected of being 3GCR were isolated from extraintestinal sites of companion dogs: urinary bladder, skin abscess, abdominal cavity, nasal cavity, or vagina. At the time, these isolates were not defined as extraintestinal pathogenic Enterobacteriaceae, because commensal Enterobacteriaceae can cause extraintestinal infection when the host is compromised [63]. These 12 isolates were obtained from the outsourcing service, and confirmed to be 3GCR by an MIC test (Table 2). Nine of 12 3GCR Enterobacteriaceae isolates were also confirmed to be resistant to fluoroquinolone (IDs: 001Sk, 019S, 019Sk, 024S, 026S, 056S, 066S, 082S, and 123S in Table 2). Fluoroquinolone-resistant E. coli is very widespread in humans, and FQR is often accompanied by 3GCR, which is mainly caused by the production of ESBLs [64]. For instance, E. coli ST131 evolved from clade B (FQ sensitive) to clade C (FQR) by chromosomal point mutations in the late 1980s, and then acquired plasmids encoding ESBLs, which led to global domination by this clade [65,66]. Enterobacteriaceae with both FQR and 3GCR were also recently reported in companion animals [67,68]. Given that 75% of tested Enterobacteriaceae isolates (9 of 12) in this study were both FQR and 3GCR, this type of Enterobacteriaceae may have become dominant among 3GCR Enterobacteriaceae in companion animals in Osaka, Japan. The 12 isolates were further analyzed for β-lactamase production by disc diffusion and synergy assays (Table 3). Five isolates (001Sk, 019Sk, 056S, 082S, and 123S) produced ESBLs, and six isolates (019S, 024S, 026S, 066S, 112S, and 128S) were AmpC-producing Enterobacteriaceae. The Proteus vulgaris isolate (048Sp) was susceptible to cefotaxime and ceftazidime in the disc diffusion assay, and thus was not characterized for β-lactamase production. All isolates were tested negative for metallo-β-lactamase production.

thumbnail
Table 2. MICs of 3GCR Enterobacteriaceae isolated from extraintestinal sites and feces.

https://doi.org/10.1371/journal.pone.0246482.t002

thumbnail
Table 3. Disc diffusion and synergy test results for 3GCR Enterobateriaecae isolated from extraintestinal sites and feces.

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

Dog feces were then examined for the presence of the ESBL/AmpC-producing Enterobacteriaceae. Fourteen 3GCR Enterobacteriaceae isolates were isolated from 12 dogs (Table 2). Nine of 14 isolates (001Fk, 019Fk, 024F, 026F1, 026F2, 048F, 049F, 082F, and 123F) were ESBL-producing E. coli or K. pneumoniae. Three isolates (019F, 026F3, and 066F) were AmpC-producing E. coli. Two isolates (066Fc and 128F) were both ESBL- and AmpC-producing Enterobacteriaceae (Table 3). These results showed that 81.8% (9 of 11) of the tested companion dogs harbored ESBL/AmpC-producing bacteria in their intestinal tracts.

WGS-based characterization of ESBL/AmpC-producing E. coli isolates from companion dogs

We focused on 19 E. coli isolates from the extraintestinal sites and feces of 11 dogs (S1 Fig). These isolates belonged to seven STs: ST131 (n = 6), ST162 (n = 5), ST68 (n = 2), ST405 (n = 2), ST998 (n = 2), ST10 (n = 1), and ST38 (n = 1). In addition, the isolates harbored one, two, or three copies of ESBL/AmpC genes encoding TEM-1A (n = 4), TEM-1B (n = 4), TEM-1C (n = 1), CTX-M-14 (n = 1), CTX-M-15 (n = 3), CTX-M-27 (n = 6), and CMY-2 (n = 10). To address whether E. coli isolates from companion dogs can infect humans, 14 E. coli isolates were whole-genome sequenced (Tables 4, 5 and S1), and their genomic similarities with strains isolated from humans were evaluated.

thumbnail
Table 4. WGS for 14 strains, with particular respect to resistance genes1.

https://doi.org/10.1371/journal.pone.0246482.t004

thumbnail
Table 5. WGS for 14 strains, with particular respect to virulence genes1.

https://doi.org/10.1371/journal.pone.0246482.t005

First, we investigated the genetic relatedness between multiple isolates from companion dogs belonging to the same ST. Two previous studies revealed that maximum SNP counts were four [69] and 15 [70] for the same outbreak of E. coli foodborne illness. To determine if multiple isolates from companion dogs with the same ST arise from the same source (belong to the same outbreak event), pairwise SNP differences were measured using core genome regions from whole-genome sequences (S2 Table). Six and eight SNP differences were identified for ST68 isolates (026S and 026F3) and ST998 isolates (128S and 128F), respectively. In both ST isolates, one isolate (026S or 128S) was isolated from the extraintestinal site and the other (026F3 or 128F) was from the feces of the same dog. These results strongly suggest that these two isolates in ST68 and ST998, respectively, arose from the same source. PFGE and the resulting dendrogram also revealed that each E. coli isolate from the extraintestinal site had the closest genetic relationship with the isolate from feces of the same dog (>85% pattern similarity) (S1 Fig). Intestinal colonization by human extraintestinal pathogenic E. coli (ExPEC) was demonstrated in companion dogs [71]. Furthermore, closely related or indistinguishable ExPEC strains (by typing methods) have been recovered from humans and their companion animals [72,73]. Taken together, these results suggest that canine extraintestinal ESBL/AmpC-producing E. coli can spread via feces, which may lead to cross-transmission to humans. Four ST131 isolates carried 77–178 SNPs and were not likely to belong to the same outbreak event (S2 Table). In ST162 isolates, three SNP differences were identified between the isolates 019S and 066S, which were derived from different companion dogs (S2 Table). These results strongly suggest that these two isolates arose from the same source (possibly from the same outbreak).

The closely related isolates were then predicted using the genome-based SNP strategy of the single genome analysis tool on the BacWGST database to track the source of the whole-genome sequenced 14 isolates from companion dogs (S3 Table). Isolates 026F2, 048F, 049F, 056S, 082S, 112S, and 123S had the closely related isolates on the BacWGST database. SNP differences ranged from 454–947 between the isolates (049F, 112S, and 123S) and their closest isolates. The presence of many (hundreds or more) SNPs indicates that the isolates were distantly related, implying that they did not originate from the same reservoir population [74]. Therefore, the sources of isolates 049F, 112S, and 123S could not be tracked or identified. Each ST131 isolate (026F2, 048F, 056S, and 082S) differed from its closest isolate by 80–247 SNPs. All the closest isolates were isolated from infection sites in human and were likely to be human ExPECs. These results suggested that ST131 isolates from companion dogs arose from human ExPECs although we were unable to identify the recent ancestor of ST131 isolates because of the high SNP differences.

Next, we focused on AMR-associated genes. Four isolates (026F2, 048F, 056S, and 082S) harbored blaCTX-M-27 in IncF incompatibility group plasmids with the pMLST type F1:A2:B20. These isolates all belonged to the ST131 and O25:H4:FimH30 types. Two isolates (123S and 128F) harbored blaCTX-M-15 in plasmids. The ST38 isolate, 049F, possessed two blaCTX-M-14 genes on its chromosome. Seven isolates (019S, 024S, 026F3, 066S, 112S, 128S, and 128F) harbored blaCMY-2, six of which harbored the gene in plasmids. We could not determine whether the blaCMY-2 gene of isolate 024S was located on its chromosome or on a plasmid, because we did not successfully determine the whole chromosomal sequence. No β-lactamase genes were detected from isolate 026S in our analysis. All the isolates, except 112S, had point mutations in quinolone resistance-determining regions of the protein. Resistance to fluoroquinolone requires the concurrent presence of at least three mutations in the target proteins encoded on the chromosome; the most common are amino acid changes in positions S83 and D87 of GyrA and S80 of ParC [75]. Two isolates, 128S and 128F, showed only one point mutation (S83L in GyrA) and also lacked the plasmid-borne quinolone resistance-associated genes qnr and aac(6)-Ib-cr, resulting in quinolone susceptibility. In contrast, the remaining isolates had at least three mutations as described above, leading to quinolone resistance (Tables 2 and 4).

(i) blaCTX-M-27 in the ST131 lineage.

Four isolates (026F2, 048F, 056S, and 082S) were of the ST131 and O25:H4:FimH30 type. These isolates all harbored blaCTX-M-27 in multireplicon plasmid depicting IncFIA/FIB/FII replicons with the pMLST type F1:A2:B20, and acquired fluoroquinolone resistance by point mutations in parC (S80I and E84V), parE (I529L), and gyrA (S83L and D87N). These results indicate that these four isolates belong to the C1-M27 clade within the C1/H30R clade [65,66]. In 2006, C1-M27 isolates of the ST131 lineage were detected in humans in Japan, and the numbers of this clade escalated in the 2010s [76]. Recently, this clade emerged and spread in France and Germany [77,78]. These reports showed that the ST131/C1-M27 clade is responsible for the substantial increase in ESBL-producing ExPEC in humans. We considered whether C1-M27 isolates could be transmitted between dogs and humans, and therefore evaluated genomic similarities between dog and human isolates. The chromosomal sequences of four C1-M27 isolates from companion dogs were compared with those from the human-derived ST131 lineage C1-M27, the clade C1/H30R harboring blaCTX-M-14, and the clade C2/H30Rx harboring blaCTX-M-15 [78,79], extracted from the NCBI database. A phylogenetic analysis based on genome-wide SNPs revealed that isolates in this study (underlined in Fig 1A) were highly similar to the human-derived C1-M27 clade (colored red in Fig 1A). In particular, E. coli str81009 (a human ExPEC isolated from human urine in UAE in 2009) carried 23–36 SNP differences against the isolates 026F2, 048F, and 056S from companion dogs (S4 Table). These SNP differences (23–36) between E. coli str81009 and the isolates 026F2, 048F, and 056S were lower than those among companion dogs (48–61 SNPs). Taken together, the results suggest that each of the three isolates 026F2, 048F, and 056S arose from human ExPEC and then evolved independently. We next focused on the IncF/pMLST F1:A2:B20 plasmid harboring blaCTX-M-27. The plasmid from 026F2 showed 99.2% nucleotide identity and 100% query coverage with the human-derived pH105 plasmid by blastn search (Fig 1B). pH105 originated from a human-derived E. coli H105 strain that was isolated in Germany in 2010, the chromosome of which shared high similarity with that of isolates 026F2, 048F, and 056S in Fig 1A and S4 Table. Two of the four plasmids (plasmid 2 of 026F2: 026F2p2 and 048Fp2) had genes conferring resistance to aminoglycosides (aadA5, aph(3′′)-1b, and aph(6)-1d), macrolides (mph(A)), tetracycline (tetA), sulfonamide (sul1 and sul2), and trimethoprim (dfrA17), which were identical to those in pH105 (Fig 1C). All the plasmids had an IS26-ΔISEcp1-blaCTX-M-27-ΔIS903D/IS903D-IS26-like structure [8082] (Fig 1D). The structures shared 100% identity in 026F2p2, 048Fp2, and 056Sp3, and were common in the human C1-M27 clade [76,83]. By contrast, the 082Sp2 plasmid contained an intact IS903D element between blaCTX-M-27 and IS26 that appeared to be unique, because there is no report describing an IS26-ΔISEcp1-blaCTX-M-27-IS903D-IS26 structure to the best of our knowledge. Conjugative transfer assays revealed that these four F1:A2:B20 plasmids harboring blaCTX-M-27 could be transmitted horizontally (Table 6). These results strongly suggest that the C1-M27 clade of the ST131 lineage is transmitted between companion dogs and humans, and also that the F1:A2:B20 plasmid with blaCTX-M-27 can be transferred between strains colonizing the two host species.

thumbnail
Fig 1. Characteristics of ESBL-producing E. coli ST131 isolates.

(A) A maximum likelihood chromosomal phylogenetic tree based on core-genome wide SNPs depicts four ESBL-producing E. coli isolates of the ST131 lineage C1-M27 from companion dogs (underlined) along with those from the human-derived ST131 lineage C1-M27 (red), clade C1/H30R harboring blaCTX-M-14 (green), and clade C2/H30Rx harboring blaCTX-M-15 (blue), which were extracted from the NCBI database. E. coli JJ1886 was used as a reference strain. The chromosomal sequence of the ST998 isolate 128S in this study was used as the outgroup. (B) The blaCTX-M-27-coding plasmids isolated in this study were compared with a highly similar plasmid, pH105, from E. coli strain 105. Homologous regions are shaded gray. The orange bars are enlarged in Fig 1C. (C) An enlargement of partial sequences of the blaCTX-M-27-coding plasmid depicted in Fig 1B. Genetic elements are indicated by colors as follows: blaCTX-M-27, red; other antimicrobial resistance genes, orange; IS26 tnpA, blue; IS903D tnpA, green; other mobile genetic elements, black. Genes are indicated by numbers as follows: dfrA17, 1; aadA5, 2; sul1, 3; mph(A), 4; sul2, 5; aph(3′′)-Ib, 6; aph(6)-Id, 7; tetA, 8. (D) Details of the blaCTX-M-27 transposon elements of 026F2p2, 048Fp2, 056Sp3 (upper panel), and 082Sp2 (lower panel). IRL and IRR are inverted repeats of each insertion sequence. Nucleotide sequences of each IR are as follows: IS26-IRL, GGCACTGTTGCAAA; IS26-IRR, TTTGCAACAGTGCC; IS903D-IRL, GGTTTTGTTGAATAAATC; IS903D-IRR, GATTTATTCAACAAAGCC; ISEcp1-IRR, ACTGTTAATTTAGG. Accession numbers of human-derived ST131 isolates and plasmids are as follows: JJ1886, PRJNA218163; TO148, PRJEB27474; EC958, PRJNA283246; B36, PRJEB29930; AZ162, SAMN06187728; AR_0055, PRJNA292904; TO143, PRJEB27473; TO60, PRJEB27477; str81009, PRJNA383781; H105, PRJNA387731; BRG120, SAMD00044968; KUN3594, SAMD00044940; KSEC29, SAMD00044934; BRG274, SAMD00044971; KN94, SAMD00044963; KS46, SAMD00044955; ECpH105, CP021871.1.

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

(ii) Plasmids harboring the blaCTX-M-15 gene.

Two isolates (123S and 128F) harbored the plasmid coding for blaCTX-M-15. We failed to determine the whole plasmid sequence in 128F, but analysis of contigs c5 and c6 of 128F (128Fc5c6) revealed that it belonged to the incompatibility group IncI1/pMLST16, and carried two ESBL genes, blaCTX-M-15 and blaTEM-1B. This plasmid showed high similarity (97.9% identity, 95.0% query coverage) to the previously published plasmid pEC_Bactec from a horse-derived E. coli strain that was isolated in 2010 [84] (Fig 2A). In comparison with the reference plasmid for the IncI1 group pR64 from Salmonella enterica, pEC_Bactec arose by transposition of Tn2 carrying blaTEM and ISEcp1-blaCTX-M-15 [84]. 128Fc5c6 had a conserved transposon–ESBL gene element (Fig 2B) [8587]. The blaTEM-1B gene was located in a Tn2-like transposon that possessed intact 38-bp inverted repeats (IR, colored in green in Fig 2B). The Tn2-like transposon was flanked by 5-bp direct repeats (DR1 in Fig 2B). The tnpA gene of the Tn2 transposon was disrupted between nucleotides 209 and 210 by transposition of the ISEcp1-blaCTX-M-15 element. ISEcp1-mediated transposition created a 5-bp duplication of the target sequence (DR2 in Fig 2B). This mechanism involved the left IR (IRL) of ISEcp1, and a right alternative IR (IRR’) that resembled the IRR of ISEcp1 (IR colored in blue in Fig 2B). The plasmid shared more than 95% homology with pEC545 from a human-derived E. coli isolate (from Vietnam in 2017) and pKHSB1 from human-derived Shigella sonnei, which was prevalent among humans in Vietnam and Korea [88,89] (Fig 2A). IncI1/pMLST16 harboring blaCTX-M-15 was found in various bacterial species, continents, and host species; for instance, human-derived E. coli isolates in Ireland [90], cattle-derived E. coli isolates in the UK [90], human-derived S. enterica isolates [91], and the human-derived S. sonnei isolates in Vietnam and Korea described above [89]. Given that the plasmid containing 128Fc5c6 was capable of horizontal transfer (Table 6), we strongly suggest that this type of plasmid is shared between humans and companion dogs. The plasmids have spread in several continents, hosted by both zoonotic pathogens and commensal bacterial species, circulating among both humans and animals. Companion dogs might contribute to the spread of this type of plasmid among humans.

thumbnail
Fig 2. Characteristics of blaCTX-M-15-coding plasmids.

(A) IncI1/pMLST16 contig 128Fc5c6 coding blaCTX-M-15 was compared with a prototype or the highly homologous plasmids pR64 from S. enterica, pEC_Bactec from E. coli, pEC545 from E. coli, and pKHSB1 from S. sonnei. Homologous regions are shaded gray. Genetic elements are indicated by colors as follows: blaCTX-M-15, red; blaTEM-1, orange; ISEcp1 tnpA, blue; Tn2 tnpA, green; other mobile genetic elements, black. (B) Details of the blaCTX-M-15 transposon element of 128Fc5c6 in Fig 2A. The basepair numbering corresponds to segments of the Tn2 tnpA gene that were separated by the transposon. DR1 and DR2 indicate the 5-bp direct repeats for Tn2 (GAAAA) and ISEcp1 (TCATA), respectively. IRs in green and blue indicate the inverted repeats for Tn2 and ISEcp1, respectively. (C) The blaCTX-M-15-coding IncFIIK/K2:A−:B− plasmid 123Sp5 was compared with pKpQIL, pKF3-94, and pIMP1572 from K. pneumoniae isolates that shared highly homology with 123Sp5. Genes are indicated by colors as follows: blaCTX-M-15, red; blaTEM-1, orange; blaIMP or blaKPC, brown; other antimicrobial resistance genes, pink; ISEcp1 tnpA, blue; Tn3 tnpA, green; other mobile genetic elements, black. Genes are indicated by numbers as follows: tetA, 1; aac(6′)-Ib-cr, 2; ARR-3, 3; dfrA27, 4; aadA16, 5; sul1, 6; mph(A), 7. (D) Details of the blaCTX-M-15 transposon element of 123Sp5 in Fig 2C. The basepair numbering corresponds to segments of the Tn3 tnpA gene that were separated by the transposon. DR1 and DR2 indicate the 5-bp direct repeats for Tn3 (GTTAA) and ISEcp1 (TCATA), respectively. L, R, and R’ for the mobile elements are the left, right, and alternative right IRs, respectively. Nucleotide sequences of each IR are as follows: Tn2 and Tn3-IRL, GGGGTCTGACGCTCAGTGGAACGAAAACTCACGTTAAG; Tn2 and Tn3-IRR, CTTAACGTGAGTTTTCGTTCCACTGAGCGTCAGACCCC; ISEcp1-IRL, CCTAGATTCTACGT; ISEcp1-IRR, ACGTGGAATTTAGG; ISEcp1-IRR’, ACGTAGGTCCCAGG; Tn1721-IRL, GGGGAGCCCGCAGAATTCGGAAAAAATCGTACGCTAAG; IS26-IRL, GGCACTGTTGCAAA; IS26-IRR, TTTGCAACAGTGCC. Accession numbers of database-derived plasmids are as follows: pR64, AP005147; pEC_Bactec, GU371927; pEC545, CP018975; pKHSB1, HF572032; pKpQIL, GU595196; pKF3-94, FJ876826; pIMP1572, MH464586.

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

The other plasmid, 123Sp5, belonged to the incompatibility group IncFIIK with the pMLST type K2:A−:B−. This replicon plasmid type originated from K. pneumoniae [92,93] and was found in several E. coli isolates [94], indicating that IncFIIK/K2:A−:B− can transfer from K. pneumoniae to other bacterial species. Conjugal assays in this study revealed that 123Sp5 was capable of horizontal transfer (Table 6), suggesting that 123Sp5 can be transmitted between these two bacterial species. Indeed, 123Sp5 shared high homology with plasmids from K. pneumoniae: IncFIIK early plasmid pKF3-94 [95] and pIMP1572 [96] (nucleotide identity, 99.97%; query coverage, 79% and 98%, respectively), and all the plasmids harbored blaCTX-M-15 (Fig 2C). pKpQIL, an IncFIIK2 (K2:A−:B−) group plasmid, is one of the most common blaKPC-harboring plasmids [94,97]. 123Sp5 shared less homology (96.81% identity and 51% query coverage) with pKpQIL than it did with pKF3-94 and pIMP1572, and did not possess the carbapenem-resistance gene. 123Sp5 was closely related to pKF3-94, but possessed genes conferring multidrug resistance to aminoglycosides (aadA16), trimethoprim (dfrA27), fluoroquinolones (aac(6cl-Ib-cr), tetracycline (tetA), macrolides (mph(A)), and sulfonamide (sul1), in addition to the EBSL genes. The blaTEM-1C gene was located in a Tn3-like transposon. The Tn3 tnpA gene was disrupted between nucleotides 372 and 373 by the ISEcp1-blaCTX-M-15 element [8487] (Fig 2D). This structure was highly similar to that in 128Fc5c6 (Fig 2B). Furthermore, Tn1721 tnpA composing the tetA mobile element was disrupted by the Tn3 transposon element carrying blaTEM-1C and ISEcp1-blaCTX-M-15 (Fig 2D). This indicates that blaTEM-1C and blaCTX-M-15 were inserted after the tetA element was translocated; and therefore, pKF3-94 did not appear to be the direct ancestor of 123Sp5. Another closely related plasmid, pIMP1572, additionally possessed blaIMP, potentially leading to carbapenem resistance. In this study, it was not clear whether the 123Sp5-like plasmid or E. coli carrying a 123Sp5-like plasmid was transmitted between humans and companion dogs. However, the prevalence of this type of plasmid should be longitudinally monitored, because mobile elements may be easily inserted or deleted on the locus coding for multidrug-resistance genes, and the plasmid can be transmitted horizontally outside of the bacterial species.

(iii) Chromosomal insertion of blaCTX-M-14.

In several bacterial species, chromosomal blaCTX-M (e.g., blaCTX-M-2, blaCTX-M-14, and blaCTX-M-15) were identified likely because of the transposition of blaCTX-M from a plasmid into the chromosome [98100]. The isolate 049F (ST38) possessed two chromosomal copies of blaCTX-M-14 (Fig 3A). The most closely related E. coli strain, 545, possessed the plasmid pEC545, in which the transposon element, ISEcp1-blaCTX-M-14-lamB, shared 100% homology with that in strain 049F (nucleotide positions 89517–93993 in Fig 3B). These results suggested that the blaCTX-M-14 gene in strain 049F was derived from a plasmid-borne transposon element, although the target site duplication (DR in Fig 3B) could not be identified. The other transposon element harboring blaCTX-M-14 in strain 049F (nucleotide positions 31540–35673 in Fig 3B) was identical to, but 343 bp shorter than, the larger element (nucleotide positions 89517–93993). These results suggested that the larger element was partially transferred (e.g., nucleotide positions 89517–93650 in Fig 3B) into the chromosome by the recurrent transposition that was observed in E. coli isolates with chromosomal blaCTX-M-14 [101]. The larger transposon element was inserted adjacent to yicI. Hamamoto et al. [101] found five genes adjacent to the insertion of a blaCTX-M-14 transposition element: tetC, ompN, yhjQ, sraG, and dacD. This study newly identified the yicI gene adjacent to the ISEcp1-blaCTX-M-14-lamB transposon element. Chromosomal blaCTX-M-14 was detected in various STs of human-derived E. coli isolates [98]. To the best of our knowledge, this was the first report of ST38 carrying chromosomal blaCTX-M-14, although ST38 possessing a blaCTX-M-14-coding plasmid was previously reported [64]; and therefore, the stability of chromosomal blaCTX-M-14, and its ability for cross-species transmission remain unclear. The chromosomal location of blaCTX-M-14 is considered to be one factor contributing to the high prevalence of these E. coli isolates in humans [98]. Continual surveillance in companion dogs is needed.

thumbnail
Fig 3. Characteristics of a chromosomal segment coding for blaCTX-M-14 from an E. coli isolate.

The partial chromosomal sequence from strain 049F was compared with that from the most closely related E. coli strain 545 (accession no. NZ_CP018976). Homologous regions are shaded gray. Genetic elements are indicated by colors as follows: blaCTX-M-14, red; lamB, pink; yicI, green; ISEcp1 tnpA, blue; other mobile genetic elements, black. Regions (a) and (b) are enlarged in Fig 3B. (B) Details of the blaCTX-M-14 transposon elements of strain 049F in Fig 3A. Numbers indicate chromosomal nucleotide positions. DRs and IRs indicate direct repeats and inverted repeats, respectively. Nucleotide positions 89517–93650 are identical to positions 31540–35673. Nucleotide sequences of each IR are as follows: ISEcp1-IRL, CCTAGATTCTACGT; ISEcp1-IRR, CGTGGAATTTAGGG. The accession number for pEC545 is CP018973.

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

(iv) Plasmid harboring blaCMY-2.

The most commonly reported AmpC β-lactamase is CMY-2. Plasmid-borne CMY-2 is distributed worldwide particularly in E. coli and Salmonella spp. from humans and animals [14,102]. Livestock animals, poultry in particular, are considered reservoirs for these bacteria and are a possible source of blaCMY-2 in humans [10]. This study showed that 42.9% of strains (6 of 14) possessed the plasmid harboring blaCMY-2 (Table 4). We did not determine whether the blaCMY-2 gene in the isolate 024S was located on the chromosome or a plasmid, because we failed to determine the whole chromosomal sequence. These results suggest that plasmid-borne CMY-2 is disseminated among companion dogs in Osaka, Japan. Various plasmid types, including IncA/C, IncI1, IncI2, and IncFII are associated with the presence of blaCMY-2 [90]. IncA/C and IncI1 plasmids, in particular, are often reported as the most common carriers of blaCMY-2 [90,103]. Of the six plasmids in this study, two (026F3c3 and 112Sc2) belonged to the IncI1 group. Two other plasmids (019Sp3 and 066Sp3) belonged to the IncFII group. The remaining two plasmids (128Sp3 and 128Fc7) belonged to the IncC (formerly IncA/C2) group.

We first focused on the IncI1 incompatibility group (plasmids 026F3c3 and 112Sc2). IncI1 has become one of the most common plasmid families in contemporary Enterobacteriaceae from both human and animal sources, and was confirmed as a transmitter of AmpC genes, particularly the blaCMY-2 gene, in isolates from livestock animals [104]. The 112Sc2 plasmid belonged to the pMLST12 which, in the IncI1 incompatibility group, has especially contributed to the wide spread of blaCMY-2 [104], and now features an epidemic plasmid in poultry and poultry products [10,105]. In the USA, the IncI1/pMLST12 plasmid harboring blaCMY-2 was detected in several companion dogs [106]. The blaCMY-2 gene of 112Sc2 was within a transposon-like structure consisting of ISEcp1-blaCMY-2-blc-sugE (Fig 4A), which is well conserved in IncI1/pMLST12 plasmids harboring blaCMY-2 in E. coli and S. enterica strains [107,108]. The transposon element was inserted in the yagA gene of the IncI1 prototype plasmid pColIb-P9 [109] from S. sonnei, flanked by the 5-bp DR (Fig 4A). The 112Sc2 plasmid exhibited 97.94% identity and 87% query coverage with pColIb-P9 (pMLST11), and was closely related to several pMLST12 plasmids: p19C93.1 from avian-derived E. coli in Japan in 2007 (99.91% identity and 94% query coverage), pJB10 [105] from avian-derived E. coli in Brazil in 2015 (99.99% identity and 88% query coverage), and pCVM29188 [110] from avian-derived S. enterica in the USA in 2003 (99.99% identity and 97% query coverage) (Fig 4B). The 112Sc2 plasmid was also highly similarity (99.81% identity and 84% query coverage) to the pMLST23 plasmid, pCVM22462 [106], from an S. enterica isolate from a dog. The other IncI1 plasmid, 026F3c3, was untypable, likely because the whole plasmid sequence was not determined. The transposon-like element ISEcp1-blaCMY2-blc-sugE was conserved, as in 112Sc2, but the element was inserted in the finQ gene, flanked by a 5-bp DR (Fig 4C). The disruption of the finQ gene by the blaCMY-2 transposon element was first reported in the S. enterica serotype Choleraesuis SCB67 strain in 2004 [111]. In addition, Su et al. [107] revealed that the transposon element carrying blaCMY-2 was inserted in the finQ gene of the transfer region of pColIb-P9 in most of the tested Enterobacteriaceae isolates, including K. pneumoniae, E. coli, S. enterica, Citrobactor spp., S. sonnei, and M. morganii. The 026F3c3 plasmid did not share high homology with the plasmid from S. Choleraesuis SCB67 (99.69% identity and 47% query coverage), suggesting that 026F3c3 was not directly related to the plasmid from S. Choleraesuis SCB67. The finQ gene appears to be a hot spot for insertion of the blaCMY-2 transposon element. Comparison of whole plasmid sequences revealed that 026F3c3 shared 99.23% identity and 88% query coverage with pColIb-P9 (Fig 4D). One of the plasmids most closely related to 026F3c3 was psg_wt5, from avian-derived S. enterica in Singapore in 2016 (99.53% identity and 93% query coverage), which lacked a blaCMY-2 transposon element. Plasmid 026F3c3 also shared high homology (98.67% identity and 96% query coverage) with pCMY-2 from human-derived pathogenic E. coli in Taiwan in 2013, which can be transformed into other bacterial species, such as K. pneumoniae [112]. We could not analyze 024S, because it was unclear whether blaCMY-2 was located on the chromosome or a plasmid, although the partial sequence was highly homologous to that in the 026F3p3 plasmid. For both 112Sc2 and 026F3c3, highly homologous plasmids were identified in unrelated bacterial species isolated in different countries and eras, suggesting that the plasmids have been spreading in several continents over many years, hosted by both zoonotic pathogens and commensal bacterial species, circulating in both humans and animals, including companion dogs.

thumbnail
Fig 4. Characteristics of blaCMY-2-coding plasmids.

(A) (C) (E) (G) Details of the blaCMY-2 transposon elements of the plasmids 112Sc2 (A), 026F3c3 (C), 019Sp3 and 066Sp3 (E), and 128Sp3 and 128Fc7 (G). DRs indicate 4–5-bp direct repeats: TGGGT (A), GATAA (C), GTTC (E), and ATTTC (G). IRs indicate inverted repeats. oriIS and terIS in Fig 4E are oriented and terminated insertion sequences, respectively. (B) (D) (F) (H) The plasmids encoding blaCMY-2 were compared with a prototype and highly homologous plasmids from E. coli (EC), S. enterica (SE), and S. soneii (SS). Genetic elements are indicated by colors as follows: blaCMY-2, red; other resistance genes, orange; blc and sugE, pink; genes disrupted by the ISEcp1-blaCMY-2 cassette, green; tnpA genes for ISEcp1 and IS1294, blue; other mobile genetic elements, black. Genes are indicated by numbers as follows: tetA, 1; aph(6)-Id, 2; aph(3′′)-Ib, 3; sul1, 4; mph(A), 5; floR, 6. Nucleotide sequences of each IR are as follows: ISEcp1-IRL, CCTAGATTCTACGT; ISEcp1-IRR, ACGTGGAATTTAGG; ISEcp1-IRR’ in Fig 4A, GAAGCGTTTGCAGG; ISEcp1-IRR’ in Fig 4C, AACCAGAAAGTCGA; IS26-IRL, GGCACTGTTGCAAA; IS26-IRR, TTTGCAACAGTGCC; IS1294-oriIS, ACGTATAGGAAATTGAAAAAC; IS1294-terIS, TTGCAGCCGCCAGGCTGCCGTG. Accession numbers of database-derived plasmids are as follows: pColIb-P9, AB021078; p19C93.1, LC501481; pJB10, KX452392; pCVM29188, NC011077; pCVM22462, CP009566; psg_wt5, CP037994; pCMY2, LC019731; p74, CP023389; pAR-0429-1, CP044142; p23C50-1, LC501563; p3, CP031549.

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

We next focused on the IncFII/pMLST F40:A−:B− plasmids 019Sp3 and 066Sp3 from different companion dogs. The transposon element IS1294-blaCMY-2-blc-sugE was conserved in the two plasmids. Furthermore, this element was inserted between the genes encoding a hypothetical protein and colicin, flanked by a 4-bp DR in both plasmids [113] (Fig 4E). The two plasmids shared 100% homology; moreover, the core-genome SNP count was three between 019S and 066S (S2 Table). These results indicate that the ST162 E. coli isolates found carrying the IncFII/pMLST F40:A−:B− plasmid with blaCMY-2 can spread among companion dogs. These isolates may be prevalent (or outbreak) among companion dogs in Osaka, Japan. Blastn search showed that one of the most closely related plasmids, p74, from E. coli strain 1105 (97.96% identity and 86% query coverage) was derived from a dog that lacked the blaCMY-2 transposon element (Fig 4F). The most closely related plasmid from the human-derived E. coli strain 2014C-3741 shared less homology (95.10% identity and 81% query coverage) with that from the dog-derived strain. The plasmids from other bacterial species also showed low homologies: 93.92% identity and 74% coverage with pEG356 from S. sonnei, 96.63% identity and 66% coverage with pKP12226 from K. pneumoniae, 98.57% identity and 80% coverage with p12-4374_62 from S. enterica. Furthermore, to the best of our knowledge, ST162 E. coli is not prevalent among humans [14,114116]. Thus, there is no evidence of ST162 carrying the IncFII/pMLST F40:A−:B− plasmid with blaCMY-2 being transferred between humans and companion dogs.

The other plasmids, 128Sp3 and 128Fc7, were derived from the same dog. They shared 100% sequence homology and belonged to the IncC (formerly IncA/C2) group. Genetic structures surrounding blaCMY-2 consisted of an ISEcp1 element upstream, and blc and IS26 downstream [117,118] (Fig 4G). Fig 4A, 4C and 4E and previous studies [14,107,108] strongly suggest that the blaCMY-2 transposon element is composed of IS-blaCMY-2-blc-sugE. The sugE gene, DR, and IRR’ for ISEcp1 in the 128Sp3 and 128Fc7 plasmids appeared to be disrupted by the insertion of the IS26-mph(A) element, although the duplicated target sequences for IS26 were not identified in this study. Disruption of the sugE gene by the IS26-mph(A) element indicated that blaCMY-2 was present before mph(A) was inserted. The two plasmids also harbored genes for resistance to tetracycline (tetA), aminoglycosides (aph(6)-Id and aph(3′′)-Ib), sulfonamide (sul2), macrolides (mph(A)), and chloramphenicol (floR) (Fig 4H). IncC plasmids are often reported as one of the most common carriers of blaCMY-2 [103,119]. We therefore examined the whole plasmid sequence of our plasmids for homology with those from the NCBI database. The plasmids pAR0429-1 from human-derived E. coli in the USA, p23C50-1 from avian-derived E. coli in Japan in 2011, and p3 from swine-derived E. coli cq9 in China in 2012 were highly homologous (nucleotide identity, 99.94%, 99.61%, and 99.98%, respectively; query coverage, 96%, 100%, and 98%, respectively) with 128Sp3 and 128Fp7 (Fig 4H). All the highly homologous plasmids possessed the intact transposon element ISEcp1-blaCMY-2-blc-sugE, and had conserved resistance genes tetA, aph(6)-Id, aph(3′′)-Ib, sul1, and floR. The most variable position was that contiguous to the blaCMY-2 transposon element. The plasmids pAR-0429-1, p23C50-1, and p3 had a cassette >56 kb that contained genes for conjugative transfer, ferredoxin, DNA methyltransferase, and drug resistance. The role of this cassette remains unclear, and it might restrict plasmid transfer between humans and animals. This type of plasmid possesses multidrug-resistance genes; therefore, continuous surveillance for cross-species transmission is needed.

(v) Virulence genes.

Companion animals are thought to be reservoirs of human ExPEC, because closely related or indistinguishable human ExPEC strains were recovered from humans and their companion animals [71,73,120]. We therefore focused on virulence-associated genes in E. coli isolates for which the whole genome was analyzed (Table 5). Isolation of an E. coli strain from a patient with an extraintestinal infection does not, by itself, confer the designation of ExPEC, because commensal E. coli strains can participate in extraintestinal infections when an aggravating or unfavorable factor is present, such as a foreign body, a compromised host, or a high-density or mixed bacterial species inoculum [63,121]. Therefore, all isolates were screened for human ExPEC status using an established molecular definition consisting of the presence of more than two of five human ExPEC-defining markers, i.e., papAH and/or papC, sfa and/or foc, afa and/or dra, iutA, and kpsMII [122]. We could not determine whether the isolates were canine ExPECs, because virulence-associated genes have different prevalences in humans and animals [123], and there are no molecular definitions for canine ExPEC. The following eight isolates belonged to the human ExPEC group: four ST131 C1-M27 isolates carrying IncF/pMLST F1:A2:B20 with blaCTX-M-27 (026F2, 048F, 056S, and 082S), an ST68 isolate carrying IncI1 with blaCMY-2 (026F3), an ST405 isolate carrying IncFIIK/pMLST K2:A−:B− with blaCTX-M-15 (123S), and two ST998 isolates carrying IncC/pMLST 3 with blaCMY-2 and IncI1/pMLST 16 with blaCTX-M-15 (128S and 128F). Human ExPEC ST131 isolates, in particular, showed extensive spread among human and companion animals [115], and ST131 C1-M27 isolates carrying blaCTX-M-27 in this study were closely related to those derived from humans (Fig 1 and S4 Table), strongly suggesting that certain E. coli isolates that produce ESBLs are transmitted between humans and companion dogs, and cause extraintestinal infections in both host species. The remaining five isolates were not human ExPECs: ST162 carrying blaCMY-2 (019S, 024S, and 066S), ST38 carrying blaCTX-M-14 (049F), and ST10 carrying blaCMY-2 (112S). We consider ST162 carrying the IncFII/pMLST F40:A−:B− plasmid with blaCMY-2 were likely not transmitted between humans and companion dogs as described above under ‘Plasmid harboring blaCMY-2’, because, to the best of our knowledge, there is no report of ST162 E. coli being prevalent in humans [14,114116]. However, there is another possibility; ST162/blaCMY-2 can be transferred between humans and companion dogs, and causes extraintestinal infections in companion dogs but not in humans, because ST162 is not human ExPEC. It is essential to evaluate not only pathogenic but also commensal E. coli to elucidate the precise role of companion dogs in cross-species transmission.

In conclusion, whole-genome analysis revealed that companion dogs in Osaka, Japan possess bla genes that were diverse in genotype, location, and type of carrying plasmid. Certain E. coli isolates producing ESBL/AmpC (i.e., ST131/blaCTX-M-27) appeared capable of cross-species transmission and causing extraintestinal infections in both host species. By contrast, other isolates (e.g., ST162/blaCMY-2) were likely not transmitted between humans and companion dogs. These results suggest that specific E. coli strains can spread between humans and companion dogs. Previous reports on methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus spp. (VRE) reflect this idea; a specific sequence type of MRSA was shared between companion animals and humans with no apparent adaptation to animals [124], and a particular form of a transposon described in only human VRE was found in VRE isolated from a dog [124]. High homologies between the tested dog-derived plasmids and those originating from other hosts, in particular, humans, and their horizontal transfer abilities (Table 6) strongly suggest that plasmids harboring the bla gene (i.e., IncFII/F1:A2:B20/blaCTX-M-27, IncI1/16/blaCTX-M-15, IncI1/blaCMY-2) can spread among humans, companion dogs, and other animals. Further investigation is needed to confirm cross-species transmission of bacteria or plasmids encoding ESBL or AmpC β-lactamase between humans and companion animals.

Supporting information

S1 Fig. Dendrogram of ESBL/AmpC-producing E. coli isolates from companion dogs.

XbaI-digested genomes of 19 ESBL/AmpC-producing E. coli isolates from companion dogs were subjected to PFGE. The dendrogram shows genotypic relatedness. The scale bar represents percent similarity. S, extraintestinal specimen; F, fecal sample.

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

(PDF)

S1 Table. Summary of whole-genome sequences of 14 E. coli isolates from companion dogs.

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

(XLSX)

S2 Table. Pairwise SNP matrix among whole-genome sequenced E. coli isolates from companion dogs.

E. coli isolates 019S, 026F2, 026F3, and 128S were used as reference strains to call SNPs in each ST. Numbers in the distance matrix denote the numbers of SNPs between each sample pair.

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

(XLSX)

S3 Table. Source tracking of whole-genome sequenced E. coli isolates from companion dogs.

The closely-related isolates were screened using the genome-based SNP strategy of the single genome analysis tool on the BacWGST database.

https://doi.org/10.1371/journal.pone.0246482.s004

(XLSX)

S4 Table. Pairwise SNP matrix among ESBL-producing E. coli ST131 isolates.

E. coli JJ1886 was used as a reference strain to call SNPs. Numbers in the distance matrix denote the numbers of SNPs between each sample pair.

https://doi.org/10.1371/journal.pone.0246482.s005

(XLSX)

Acknowledgments

We thank Kaoru Umeda (Osaka Institute of Public Health) for providing bacterial strains. We thank Takashi Hasegawa, Hideo Akiyoshi, Shunsuke Noguchi, Ryoji Kanegi, Marina Funayama, Hirotake Yamazaki, Munetaka Iwata, Tomoyo Nabetani, Shogo Hirata, Keisuke Momoki, Kotaro Shiba, and Yasuma Iimori (Veterinary Medical Center, Osaka Prefecture University) for technical assistance. We thank Awasthi S. Prasad (Osaka Prefecture University) and Yuki Wakabayashi (Osaka Institute of Public Health) for their helpful advice.

References

  1. 1. World Health Organization. Antimicrobial resistance. Global report on surveillance. 2014.
  2. 2. The review on antimicrobial resistance. Tackling drug-resistant infections globally: final report and recommendations. 2016.
  3. 3. Wellington EM, Boxall AB, Cross P, Feil EJ, Gaze WH, Hawkey PM, et al. The role of the natural environment in the emergence of antibiotic resistance in gram-negative bacteria. Lancet Infect Dis. 2013;13(2):155–65. Epub 2013/01/26. pmid:23347633.
  4. 4. Pitout JD, Laupland KB. Extended-spectrum beta-lactamase-producing Enterobacteriaceae: an emerging public-health concern. Lancet Infect Dis. 2008;8(3):159–66. Epub 2008/02/23. pmid:18291338.
  5. 5. Canton R, Novais A, Valverde A, Machado E, Peixe L, Baquero F, et al. Prevalence and spread of extended-spectrum beta-lactamase-producing Enterobacteriaceae in Europe. Clin Microbiol Infect. 2008;14 Suppl 1:144–53. Epub 2007/12/25. pmid:18154538.
  6. 6. Bevan ER, Jones AM, Hawkey PM. Global epidemiology of CTX-M beta-lactamases: temporal and geographical shifts in genotype. J Antimicrob Chemother. 2017;72(8):2145–55. Epub 2017/05/26. pmid:28541467.
  7. 7. Canton R, Coque TM. The CTX-M beta-lactamase pandemic. Curr Opin Microbiol. 2006;9(5):466–75. Epub 2006/09/01. pmid:16942899.
  8. 8. Jacoby GA. AmpC beta-lactamases. Clin Microbiol Rev. 2009;22(1):161–82, Table of Contents. Epub 2009/01/13. pmid:19136439; PubMed Central PMCID: PMC2620637.
  9. 9. Meini S, Tascini C, Cei M, Sozio E, Rossolini GM. AmpC beta-lactamase-producing Enterobacterales: what a clinician should know. Infection. 2019;47(3):363–75. Epub 2019/03/07. pmid:30840201.
  10. 10. Castellanos LR, van der Graaf-van Bloois L, Donado-Godoy P, Mevius DJ, Wagenaar JA, Hordijk J, et al. Phylogenomic investigation of IncI1-Igamma plasmids harboring blaCMY-2 and blaSHV-12 in Salmonella enterica and Escherichia coli in multiple countries. Antimicrob Agents Chemother. 2019;63(7). Epub 2019/05/16. pmid:31085509; PubMed Central PMCID: PMC6591584.
  11. 11. de Been M, Lanza VF, de Toro M, Scharringa J, Dohmen W, Du Y, et al. Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages. PLoS Genet. 2014;10(12):e1004776. Epub 2014/12/19. pmid:25522320; PubMed Central PMCID: PMC4270446.
  12. 12. Economou V, Gousia P. Agriculture and food animals as a source of antimicrobial-resistant bacteria. Infect Drug Resist. 2015;8:49–61. Epub 2015/04/17. pmid:25878509; PubMed Central PMCID: PMC4388096.
  13. 13. Geser N, Stephan R, Hachler H. Occurrence and characteristics of extended-spectrum beta-lactamase (ESBL) producing Enterobacteriaceae in food producing animals, minced meat and raw milk. BMC Vet Res. 2012;8:21. Epub 2012/03/09. pmid:22397509; PubMed Central PMCID: PMC3319423.
  14. 14. Pietsch M, Irrgang A, Roschanski N, Brenner Michael G, Hamprecht A, Rieber H, et al. Whole genome analyses of CMY-2-producing Escherichia coli isolates from humans, animals and food in Germany. BMC Genomics. 2018;19(1):601. Epub 2018/08/11. pmid:30092762; PubMed Central PMCID: PMC6085623.
  15. 15. Bogaerts P, Huang TD, Bouchahrouf W, Bauraing C, Berhin C, El Garch F, et al. Characterization of ESBL- and AmpC-producing Enterobacteriaceae from diseased companion animals in Europe. Microb Drug Resist. 2015;21(6):643–50. Epub 2015/06/23. pmid:26098354.
  16. 16. Ewers C, Bethe A, Semmler T, Guenther S, Wieler LH. Extended-spectrum beta-lactamase-producing and AmpC-producing Escherichia coli from livestock and companion animals, and their putative impact on public health: a global perspective. Clin Microbiol Infect. 2012;18(7):646–55. Epub 2012/04/24. pmid:22519858.
  17. 17. van den Bunt G, Fluit AC, Spaninks MP, Timmerman AJ, Geurts Y, Kant A, et al. Faecal carriage, risk factors, acquisition and persistence of ESBL-producing Enterobacteriaceae in dogs and cats and co-carriage with humans belonging to the same household. J Antimicrob Chemother. 2020;75(2):342–50. Epub 2019/11/12. pmid:31711228; PubMed Central PMCID: PMC6966097.
  18. 18. The AMR One Health Surveillance Committee. Nippon AMR One Health Report (NAOR) 2018. 2018.
  19. 19. The AMR One Health Surveillance Committee. Nippon AMR One Health Report (NAOR) 2019. 2019.
  20. 20. Guardabassi L, Schwarz S, Lloyd DH. Pet animals as reservoirs of antimicrobial-resistant bacteria. J Antimicrob Chemother. 2004;54(2):321–32. Epub 2004/07/16. pmid:15254022.
  21. 21. Lloyd DH. Reservoirs of antimicrobial resistance in pet animals. Clin Infect Dis. 2007;45 Suppl 2:S148–52. Epub 2007/08/19. pmid:17683019.
  22. 22. Albrechtova K, Dolejska M, Cizek A, Tausova D, Klimes J, Bebora L, et al. Dogs of nomadic pastoralists in northern Kenya are reservoirs of plasmid-mediated cephalosporin- and quinolone-resistant Escherichia coli, including pandemic clone B2-O25-ST131. Antimicrob Agents Chemother. 2012;56(7):4013–7. Epub 2012/04/18. pmid:22508313; PubMed Central PMCID: PMC3393422.
  23. 23. Gronthal T, Osterblad M, Eklund M, Jalava J, Nykasenoja S, Pekkanen K, et al. Sharing more than friendship-transmission of NDM-5 ST167 and CTX-M-9 ST69 Escherichia coli between dogs and humans in a family, Finland, 2015. Euro Surveill. 2018;23(27). Epub 2018/07/12. pmid:29991384; PubMed Central PMCID: PMC6152158.
  24. 24. Hong JS, Song W, Park HM, Oh JY, Chae JC, Jeong S, et al. Molecular characterization of fecal extended-spectrum beta-lactamase- and AmpC beta-lactamase-producing Escherichia coli from healthy companion animals and cohabiting humans in South Korea. Front Microbiol. 2020;11:674. Epub 2020/05/01. pmid:32351490; PubMed Central PMCID: PMC7174606.
  25. 25. Toombs-Ruane LJ, Benschop J, French NP, Biggs PJ, Midwinter AC, Marshall JC, et al. Carriage of extended-spectrum beta-lactamase-and AmpC beta-lactamase-producing Escherichia coli from humans and pets in the same households. Appl Environ Microbiol. 2020. Epub 2020/10/11. pmid:33036993.
  26. 26. Clinical and Laboratory Standards Institute. M07-A10. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically; apporoved standard. Tenth edition. 2015.
  27. 27. Clinical and Laboratory Standards Institute. M100S. Performance standards for antimicrobial susceptibility testing.
  28. 28. Arakawa Y, Shibata N, Shibayama K, Kurokawa H, Yagi T, Fujiwara H, et al. Convenient test for screening metallo-beta-lactamase-producing gram-negative bacteria by using thiol compounds. J Clin Microbiol. 2000;38(1):40–3. Epub 2000/01/05. pmid:10618060; PubMed Central PMCID: PMC86013.
  29. 29. Yagi T, Wachino J, Kurokawa H, Suzuki S, Yamane K, Doi Y, et al. Practical methods using boronic acid compounds for identification of class C beta-lactamase-producing Klebsiella pneumoniae and Escherichia coli. J Clin Microbiol. 2005;43(6):2551–8. Epub 2005/06/16. pmid:15956362; PubMed Central PMCID: PMC1151917.
  30. 30. Pomba-Feria C, Canica M. A novel sequence framework blaTEM-1G encoding the parental TEM-1 beta-lactamase. FEMS Microbiol Lett. 2003;220(2):177–80. Epub 2003/04/03. pmid:12670678.
  31. 31. Sutcliffe JG. Nucleotide sequence of the ampicillin resistance gene of Escherichia coli plasmid pBR322. Proc Natl Acad Sci U S A. 1978;75(8):3737–41. Epub 1978/08/01. pmid:358200; PubMed Central PMCID: PMC392861.
  32. 32. Yagi T, Kurokawa H, Shibata N, Shibayama K, Arakawa Y. A preliminary survey of extended-spectrum beta-lactamases (ESBLs) in clinical isolates of Klebsiella pneumoniae and Escherichia coli in Japan. FEMS Microbiol Lett. 2000;184(1):53–6. Epub 2000/02/26. pmid:10689165.
  33. 33. Shibata N, Kurokawa H, Doi Y, Yagi T, Yamane K, Wachino J, et al. PCR classification of CTX-M-type beta-lactamase genes identified in clinically isolated gram-negative bacilli in Japan. Antimicrob Agents Chemother. 2006;50(2):791–5. Epub 2006/01/27. pmid:16436748; PubMed Central PMCID: PMC1366867.
  34. 34. Perez-Perez FJ, Hanson ND. Detection of plasmid-mediated AmpC beta-lactamase genes in clinical isolates by using multiplex PCR. J Clin Microbiol. 2002;40(6):2153–62. Epub 2002/05/31. pmid:12037080; PubMed Central PMCID: PMC130804.
  35. 35. Mena A, Plasencia V, Garcia L, Hidalgo O, Ayestaran JI, Alberti S, et al. Characterization of a large outbreak by CTX-M-1-producing Klebsiella pneumoniae and mechanisms leading to in vivo carbapenem resistance development. J Clin Microbiol. 2006;44(8):2831–7. Epub 2006/08/08. pmid:16891499; PubMed Central PMCID: PMC1594597.
  36. 36. Eckert C, Gautier V, Arlet G. DNA sequence analysis of the genetic environment of various blaCTX-M genes. J Antimicrob Chemother. 2006;57(1):14–23. Epub 2005/11/18. pmid:16291869.
  37. 37. Maeyama Y, Taniguchi Y, Hayashi W, Ohsaki Y, Osaka S, Koide S, et al. Prevalence of ESBL/AmpC genes and specific clones among the third-generation cephalosporin-resistant Enterobacteriaceae from canine and feline clinical specimens in Japan. Vet Microbiol. 2018;216:183–9. Epub 2018/03/10. pmid:29519514.
  38. 38. Canica MM, Lu CY, Krishnamoorthy R, Paul GC. Molecular diversity and evolution of blaTEM genes encoding beta-lactamases resistant to clavulanic acid in clinical E. coli. J Mol Evol. 1997;44(1):57–65. Epub 1997/01/01. pmid:9010136.
  39. 39. Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, et al. Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006;60(5):1136–51. Epub 2006/05/13. pmid:16689791; PubMed Central PMCID: PMC1557465.
  40. 40. Clermont O, Gordon D, Denamur E. Guide to the various phylogenetic classification schemes for Escherichia coli and the correspondence among schemes. Microbiology (Reading). 2015;161(Pt 5):980–8. Epub 2015/02/26. pmid:25714816.
  41. 41. Umeda K, Hase A, Matsuo M, Horimoto T, Ogasawara J. Prevalence and genetic characterization of cephalosporin-resistant Enterobacteriaceae among dogs and cats in an animal shelter. J Med Microbiol. 2019;68(3):339–45. Epub 2019/01/31. pmid:30698519.
  42. 42. Matsumoto Y, Kinjo T, Motooka D, Nabeya D, Jung N, Uechi K, et al. Comprehensive subspecies identification of 175 nontuberculous mycobacteria species based on 7547 genomic profiles. Emerg Microbes Infect. 2019;8(1):1043–53. Epub 2019/07/10. pmid:31287781; PubMed Central PMCID: PMC6691804.
  43. 43. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13(6):e1005595. Epub 2017/06/09. pmid:28594827; PubMed Central PMCID: PMC5481147.
  44. 44. Kolmogorov M, Yuan J, Lin Y, Pevzner PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019;37(5):540–6. Epub 2019/04/03. pmid:30936562.
  45. 45. Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv. 2013;1303. 3997.
  46. 46. Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9(11):e112963. Epub 2014/11/20. pmid:25409509; PubMed Central PMCID: PMC4237348.
  47. 47. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75. Epub 2008/02/12. pmid:18261238; PubMed Central PMCID: PMC2265698.
  48. 48. Tanizawa Y, Fujisawa T, Kaminuma E, Nakamura Y, Arita M. DFAST and DAGA: web-based integrated genome annotation tools and resources. Biosci Microbiota Food Health. 2016;35(4):173–84. Epub 2016/11/22. pmid:27867804; PubMed Central PMCID: PMC5107635.
  49. 49. Feng Y, Zou S, Chen H, Yu Y, Ruan Z. BacWGSTdb 2.0: a one-stop repository for bacterial whole-genome sequence typing and source tracking. Nucleic Acids Res. 2021;49(D1):D644–D50. Epub 2020/10/04. pmid:33010178.
  50. 50. Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, Lund O, Villa L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903. Epub 2014/04/30. pmid:24777092; PubMed Central PMCID: PMC4068535.
  51. 51. 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. Epub 2012/07/12. pmid:22782487; PubMed Central PMCID: PMC3468078.
  52. 52. Joensen KG, Scheutz F, Lund O, Hasman H, Kaas RS, Nielsen EM, et al. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol. 2014;52(5):1501–10. Epub 2014/02/28. pmid:24574290; PubMed Central PMCID: PMC3993690.
  53. 53. Joensen KG, Tetzschner AM, Iguchi A, Aarestrup FM, Scheutz F. Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J Clin Microbiol. 2015;53(8):2410–26. Epub 2015/05/15. pmid:25972421; PubMed Central PMCID: PMC4508402.
  54. 54. Roer L, Tchesnokova V, Allesoe R, Muradova M, Chattopadhyay S, Ahrenfeldt J, et al. Development of a web tool for Escherichia coli subtyping based on fimH alleles. J Clin Microbiol. 2017;55(8):2538–43. Epub 2017/06/09. pmid:28592545; PubMed Central PMCID: PMC5527432.
  55. 55. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. Epub 2011/02/01. pmid:21278367; PubMed Central PMCID: PMC3065679.
  56. 56. Johansson MHK, Bortolaia V, Tansirichaiya S, Aarestrup FM, Roberts AP, Petersen TN. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J Antimicrob Chemother. 2021;76(1):101–9. Epub 2020/10/04. pmid:33009809; PubMed Central PMCID: PMC7729385.
  57. 57. Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16–21. Epub 2016/05/05. pmid:27141966; PubMed Central PMCID: PMC4987931.
  58. 58. Kaas RS, Leekitcharoenphon P, Aarestrup FM, Lund O. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS One. 2014;9(8):e104984. Epub 2014/08/12. pmid:25110940; PubMed Central PMCID: PMC4128722.
  59. 59. Saltykova A, Wuyts V, Mattheus W, Bertrand S, Roosens NHC, Marchal K, et al. Comparison of SNP-based subtyping workflows for bacterial isolates using WGS data, applied to Salmonella enterica serotype Typhimurium and serotype 1,4,[5],12:i. PLoS One. 2018;13(2):e0192504. Epub 2018/02/07. pmid:29408896; PubMed Central PMCID: PMC5800660.
  60. 60. Kudo H, Usui M, Nagafuji W, Oka K, Takahashi M, Yamaguchi H, et al. Inhibition effect of flavophospholipol on conjugative transfer of the extended-spectrum beta-lactamase and vanA genes. J Antibiot (Tokyo). 2019;72(2):79–85. Epub 2018/10/27. pmid:30361635; PubMed Central PMCID: PMC6760635.
  61. 61. Umeda K, Nakamura H, Fukuda A, Yamaguchi T, Matsumoto Y, Motooka D, et al. Molecular characterization of blaKHM-1 encoding plasmid in an Enterobacter hormaechei subsp. hoffmannii isolate from blood culture. PLoS One. 2020;15(1):e0227605. Epub 2020/01/14. pmid:31929575; PubMed Central PMCID: PMC6957179.
  62. 62. Centers for Disease Control and Prevention. Antibiotic resistance threats in the United States. 2019.
  63. 63. Russo TA, Johnson JR. Proposal for a new inclusive designation for extraintestinal pathogenic isolates of Escherichia coli: ExPEC. J Infect Dis. 2000;181(5):1753–4. Epub 2000/05/24. pmid:10823778.
  64. 64. Kawamura K, Nagano N, Suzuki M, Wachino JI, Kimura K, Arakawa Y. ESBL-producing Escherichia coli and its rapid rise among healthy people. Food Saf (Tokyo). 2017;5(4):122–50. Epub 2017/12/29. pmid:32231938; PubMed Central PMCID: PMC6989193.
  65. 65. Johnson TJ, Danzeisen JL, Youmans B, Case K, Llop K, Munoz-Aguayo J, et al. Separate F-type plasmids have shaped the evolution of the H30 subclone of Escherichia coli sequence type 131. mSphere. 2016;1(4). Epub 2016/07/09. pmid:27390780; PubMed Central PMCID: PMC4933990.
  66. 66. Pitout JD, DeVinney R. Escherichia coli ST131: a multidrug-resistant clone primed for global domination. F1000Res. 2017;6. Epub 2017/03/28. pmid:28344773; PubMed Central PMCID: PMC5333602.
  67. 67. Chen Y, Liu Z, Zhang Y, Zhang Z, Lei L, Xia Z. Increasing prevalence of ESBL-producing multidrug resistance Escherichia coli from diseased pets in Beijing, China from 2012 to 2017. Front Microbiol. 2019;10:2852. Epub 2020/01/11. pmid:31921034; PubMed Central PMCID: PMC6915038.
  68. 68. Loncaric I, Misic D, Szostak MP, Kunzel F, Schafer-Somi S, Spergser J. Broad-spectrum cephalosporin-resistant and/or fluoroquinolone-resistant Enterobacterales associated with canine and feline urogenital infections. Antibiotics (Basel). 2020;9(7). Epub 2020/07/11. pmid:32645942; PubMed Central PMCID: PMC7399855.
  69. 69. Underwood AP, Dallman T, Thomson NR, Williams M, Harker K, Perry N, et al. Public health value of next-generation DNA sequencing of enterohemorrhagic Escherichia coli isolates from an outbreak. J Clin Microbiol. 2013;51(1):232–7. Epub 2012/11/09. pmid:23135946; PubMed Central PMCID: PMC3536255.
  70. 70. Eppinger M, Mammel MK, Leclerc JE, Ravel J, Cebula TA. Genomic anatomy of Escherichia coli O157:H7 outbreaks. Proc Natl Acad Sci U S A. 2011;108(50):20142–7. Epub 2011/12/03. pmid:22135463; PubMed Central PMCID: PMC3250189.
  71. 71. Johnson JR, Stell AL, Delavari P. Canine feces as a reservoir of extraintestinal pathogenic Escherichia coli. Infect Immun. 2001;69(3):1306–14. Epub 2001/02/17. pmid:11179292; PubMed Central PMCID: PMC98021.
  72. 72. Johnson JR, Clabots C. Sharing of virulent Escherichia coli clones among household members of a woman with acute cystitis. Clin Infect Dis. 2006;43(10):e101–8. Epub 2006/10/20. pmid:17051483.
  73. 73. Johnson JR, Owens K, Gajewski A, Clabots C. Escherichia coli colonization patterns among human household members and pets, with attention to acute urinary tract infection. J Infect Dis. 2008;197(2):218–24. Epub 2008/01/09. pmid:18179385.
  74. 74. Pightling AW, Pettengill JB, Luo Y, Baugher JD, Rand H, Strain E. Interpreting whole-genome sequence analyses of foodborne bacteria for regulatory applications and outbreak investigations. Front Microbiol. 2018;9:1482. Epub 2018/07/26. pmid:30042741; PubMed Central PMCID: PMC6048267.
  75. 75. Tchesnokova V, Radey M, Chattopadhyay S, Larson L, Weaver JL, Kisiela D, et al. Pandemic fluoroquinolone resistant Escherichia coli clone ST1193 emerged via simultaneous homologous recombinations in 11 gene loci. Proc Natl Acad Sci U S A. 2019;116(29):14740–8. Epub 2019/07/03. pmid:31262826; PubMed Central PMCID: PMC6642405.
  76. 76. Matsumura Y, Johnson JR, Yamamoto M, Nagao M, Tanaka M, Takakura S, et al. CTX-M-27- and CTX-M-14-producing, ciprofloxacin-resistant Escherichia coli of the H30 subclonal group within ST131 drive a Japanese regional ESBL epidemic. J Antimicrob Chemother. 2015;70(6):1639–49. Epub 2015/02/18. pmid:25687644.
  77. 77. Birgy A, Bidet P, Levy C, Sobral E, Cohen R, Bonacorsi S. CTX-M-27-producing Escherichia coli of sequence type 131 and clade C1-M27, France. Emerg Infect Dis. 2017;23(5):885. Epub 2017/04/19. pmid:28418829; PubMed Central PMCID: PMC5403054.
  78. 78. Ghosh H, Doijad S, Falgenhauer L, Fritzenwanker M, Imirzalioglu C, Chakraborty T. blaCTX-M-27-encoding Escherichia coli sequence type 131 lineage C1-M27 clone in clinical isolates, Germany. Emerg Infect Dis. 2017;23(10):1754–6. Epub 2017/09/21. pmid:28930021; PubMed Central PMCID: PMC5621564.
  79. 79. Matsumura Y, Pitout JD, Gomi R, Matsuda T, Noguchi T, Yamamoto M, et al. Global Escherichia coli sequence type 131 clade with blaCTX-M-27 gene. Emerg Infect Dis. 2016;22(11):1900–7. Epub 2016/10/22. pmid:27767006; PubMed Central PMCID: PMC5088012.
  80. 80. Derbyshire KM, Hwang L, Grindley ND. Genetic analysis of the interaction of the insertion sequence IS903 transposase with its terminal inverted repeats. Proc Natl Acad Sci U S A. 1987;84(22):8049–53. Epub 1987/11/01. pmid:2825175; PubMed Central PMCID: PMC299474.
  81. 81. Harmer CJ, Hall RM. An analysis of the IS6/IS26 family of insertion sequences: is it a single family? Microb Genom. 2019;5(9). Epub 2019/09/06. pmid:31486766; PubMed Central PMCID: PMC6807381.
  82. 82. Kim EH, Aoki T. The transposon-like structure of IS26-tetracycline, and kanamycin resistance determinant derived from transferable R plasmid of fish pathogen, Pasteurella piscicida. Microbiol Immunol. 1994;38(1):31–8. Epub 1994/01/01. pmid:8052160.
  83. 83. Jamborova I, Johnston BD, Papousek I, Kachlikova K, Micenkova L, Clabots C, et al. Extensive genetic commonality among wildlife, wastewater, community, and nosocomial isolates of Escherichia coli sequence type 131 (H30R1 and H30Rx subclones) that carry blaCTX-M-27 or blaCTX-M-15. Antimicrob Agents Chemother. 2018;62(10). Epub 2018/08/01. pmid:30061277; PubMed Central PMCID: PMC6153832.
  84. 84. Smet A, Van Nieuwerburgh F, Vandekerckhove TT, Martel A, Deforce D, Butaye P, et al. Complete nucleotide sequence of CTX-M-15-plasmids from clinical Escherichia coli isolates: insertional events of transposons and insertion sequences. PLoS One. 2010;5(6):e11202. Epub 2010/06/30. pmid:20585456; PubMed Central PMCID: PMC2887853.
  85. 85. Dhanji H, Doumith M, Hope R, Livermore DM, Woodford N. ISEcp1-mediated transposition of linked blaCTX-M-3 and blaTEM-1b from the IncI1 plasmid pEK204 found in clinical isolates of Escherichia coli from Belfast, UK. J Antimicrob Chemother. 2011;66(10):2263–5. Epub 2011/07/29. pmid:21795257.
  86. 86. Nicolas E, Lambin M, Dandoy D, Galloy C, Nguyen N, Oger CA, et al. The Tn3-family of replicative transposons. Microbiol Spectr. 2015;3(4). Epub 2015/09/10. pmid:26350313.
  87. 87. Zong Z, Partridge SR, Iredell JR. ISEcp1-mediated transposition and homologous recombination can explain the context of blaCTX-M-62 linked to qnrB2. Antimicrob Agents Chemother. 2010;54(7):3039–42. Epub 2010/04/28. pmid:20421399; PubMed Central PMCID: PMC2897309.
  88. 88. Kim JS, Kim J, Jeon SE, Kim SJ, Kim NO, Hong S, et al. Complete nucleotide sequence of the IncI1 plasmid pSH4469 encoding CTX-M-15 extended-spectrum beta-lactamase in a clinical isolate of Shigella sonnei from an outbreak in the Republic of Korea. Int J Antimicrob Agents. 2014;44(6):533–7. Epub 2014/12/03. pmid:25446906.
  89. 89. Kim JS, Kim JJ, Kim SJ, Jeon SE, Seo KY, Choi JK, et al. Outbreak of ciprofloxacin-resistant Shigella sonnei associated with travel to Vietnam, Republic of Korea. Emerg Infect Dis. 2015;21(7):1247–50. Epub 2015/06/17. pmid:26079171; PubMed Central PMCID: PMC4480405.
  90. 90. Carattoli A, Villa L, Fortini D, Garcia-Fernandez A. Contemporary IncI1 plasmids involved in the transmission and spread of antimicrobial resistance in Enterobacteriaceae. Plasmid. 2018. Epub 2018/12/12. pmid:30529488.
  91. 91. Wong MH, Kan B, Chan EW, Yan M, Chen S. IncI1 plasmids carrying various blaCTX-M genes contribute to ceftriaxone resistance in Salmonella enterica serovar Enteritidis in China. Antimicrob Agents Chemother. 2016;60(2):982–9. Epub 2015/12/09. pmid:26643327; PubMed Central PMCID: PMC4750657.
  92. 92. Bi D, Zheng J, Li JJ, Sheng ZK, Zhu X, Ou HY, et al. In silico typing and comparative genomic analysis of IncFIIK plasmids and insights into the evolution of replicons, plasmid backbones, and resistance determinant profiles. Antimicrob Agents Chemother. 2018;62(10). Epub 2018/07/18. pmid:30012771; PubMed Central PMCID: PMC6153814.
  93. 93. Navon-Venezia S, Kondratyeva K, Carattoli A. Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance. FEMS Microbiol Rev. 2017;41(3):252–75. Epub 2017/05/19. pmid:28521338.
  94. 94. Chavda KD, Chen L, Jacobs MR, Bonomo RA, Kreiswirth BN. Molecular diversity and plasmid analysis of KPC-producing Escherichia coli. Antimicrob Agents Chemother. 2016;60(7):4073–81. Epub 2016/04/27. pmid:27114279; PubMed Central PMCID: PMC4914655.
  95. 95. Ying J, Wu S, Zhang K, Wang Z, Zhu W, Zhu M, et al. Comparative genomics analysis of pKF3-94 in Klebsiella pneumoniae reveals plasmid compatibility and horizontal gene transfer. Front Microbiol. 2015;6:831. Epub 2015/09/09. pmid:26347723; PubMed Central PMCID: PMC4539522.
  96. 96. Yao H, Cheng J, Li A, Yu R, Zhao W, Qin S, et al. Molecular characterization of an IncFIIk plasmid co-harboring blaIMP-26 and tet(A) variant in a clinical Klebsiella pneumoniae isolate. Front Microbiol. 2020;11:1610. Epub 2020/08/15. pmid:32793144; PubMed Central PMCID: PMC7393768.
  97. 97. Peirano G, Bradford PA, Kazmierczak KM, Chen L, Kreiswirth BN, Pitout JD. Importance of clonal complex 258 and IncFK2-like plasmids among a global collection of Klebsiella pneumoniae with blaKPC. Antimicrob Agents Chemother. 2017;61(4). Epub 2017/02/09. pmid:28167556; PubMed Central PMCID: PMC5365689.
  98. 98. Hamamoto K, Ueda S, Toyosato T, Yamamoto Y, Hirai I. High prevalence of chromosomal blaCTX-M-14 in Escherichia coli isolates possessing blaCTX-M-14. Antimicrob Agents Chemother. 2016;60(4):2582–4. Epub 2016/01/27. pmid:26810652; PubMed Central PMCID: PMC4808179.
  99. 99. Harada S, Ishii Y, Saga T, Kouyama Y, Tateda K, Yamaguchi K. Chromosomal integration and location on IncT plasmids of the blaCTX-M-2 gene in Proteus mirabilis clinical isolates. Antimicrob Agents Chemother. 2012;56(2):1093–6. Epub 2011/11/23. pmid:22106217; PubMed Central PMCID: PMC3264238.
  100. 100. Hirai I, Fukui N, Taguchi M, Yamauchi K, Nakamura T, Okano S, et al. Detection of chromosomal blaCTX-M-15 in Escherichia coli O25b-B2-ST131 isolates from the Kinki region of Japan. Int J Antimicrob Agents. 2013;42(6):500–6. Epub 2013/10/05. pmid:24091130.
  101. 101. Hamamoto K, Hirai I. Characterisation of chromosomally-located blaCTX-M and its surrounding sequence in CTX-M-type extended-spectrum beta-lactamase-producing Escherichia coli isolates. J Glob Antimicrob Resist. 2019;17:53–7. Epub 2018/11/18. pmid:30445208.
  102. 102. Borjesson S, Jernberg C, Brolund A, Edquist P, Finn M, Landen A, et al. Characterization of plasmid-mediated AmpC-producing E. coli from Swedish broilers and association with human clinical isolates. Clin Microbiol Infect. 2013;19(7):E309–11. Epub 2013/04/16. pmid:23581796.
  103. 103. Tagg KA, Iredell JR, Partridge SR. Complete sequencing of IncI1 sequence type 2 plasmid pJIE512b indicates mobilization of blaCMY-2 from an IncA/C plasmid. Antimicrob Agents Chemother. 2014;58(8):4949–52. Epub 2014/06/04. pmid:24890591; PubMed Central PMCID: PMC4135994.
  104. 104. Castellanos LR, Donado-Godoy P, Leon M, Clavijo V, Arevalo A, Bernal JF, et al. High heterogeneity of Escherichia coli sequence types harbouring ESBL/AmpC genes on IncI1 plasmids in the Colombian poultry chain. PLoS One. 2017;12(1):e0170777. Epub 2017/01/27. pmid:28125687; PubMed Central PMCID: PMC5268450.
  105. 105. da Silva KC, Cunha MP, Cerdeira L, de Oliveira MG, de Oliveira MC, Gomes CR, et al. High-virulence CMY-2- and CTX-M-2-producing avian pathogenic Escherichia coli strains isolated from commercial turkeys. Diagn Microbiol Infect Dis. 2017;87(1):64–7. Epub 2016/10/25. pmid:27773543.
  106. 106. Hansen KH, Bortolaia V, Nielsen CA, Nielsen JB, Schonning K, Agerso Y, et al. Host-specific patterns of genetic diversity among IncI1-Igamma and IncK plasmids encoding CMY-2 beta-lactamase in Escherichia coli isolates from humans, poultry meat, poultry, and dogs in Denmark. Appl Environ Microbiol. 2016;82(15):4705–14. Epub 2016/05/29. pmid:27235431; PubMed Central PMCID: PMC4984282.
  107. 107. Su LH, Chen HL, Chia JH, Liu SY, Chu C, Wu TL, et al. Distribution of a transposon-like element carrying blaCMY-2 among Salmonella and other Enterobacteriaceae. J Antimicrob Chemother. 2006;57(3):424–9. Epub 2006/01/07. pmid:16396917.
  108. 108. Verdet C, Gautier V, Chachaty E, Ronco E, Hidri N, Decre D, et al. Genetic context of plasmid-carried blaCMY-2-like genes in Enterobacteriaceae. Antimicrob Agents Chemother. 2009;53(9):4002–6. Epub 2009/07/15. pmid:19596889; PubMed Central PMCID: PMC2737857.
  109. 109. Johnson TJ, Shepard SM, Rivet B, Danzeisen JL, Carattoli A. Comparative genomics and phylogeny of the IncI1 plasmids: a common plasmid type among porcine enterotoxigenic Escherichia coli. Plasmid. 2011;66(3):144–51. Epub 2011/08/17. pmid:21843549.
  110. 110. Fricke WF, McDermott PF, Mammel MK, Zhao S, Johnson TJ, Rasko DA, et al. Antimicrobial resistance-conferring plasmids with similarity to virulence plasmids from avian pathogenic Escherichia coli strains in Salmonella enterica serovar Kentucky isolates from poultry. Appl Environ Microbiol. 2009;75(18):5963–71. Epub 2009/08/04. pmid:19648374; PubMed Central PMCID: PMC2747853.
  111. 111. Chiu CH, Su LH, Chu C, Chia JH, Wu TL, Lin TY, et al. Isolation of Salmonella enterica serotype Choleraesuis resistant to ceftriaxone and ciprofloxacin. Lancet. 2004;363(9417):1285–6. Epub 2004/04/20. pmid:15094275.
  112. 112. Lin YT, Pan YJ, Lin TL, Fung CP, Wang JT. Transfer of CMY-2 cephalosporinase from Escherichia coli to virulent Klebsiella pneumoniae causing a recurrent liver abscess. Antimicrob Agents Chemother. 2015;59(8):5000–2. Epub 2015/05/20. pmid:25987637; PubMed Central PMCID: PMC4505199.
  113. 113. Tavakoli N, Comanducci A, Dodd HM, Lett MC, Albiger B, Bennett P. IS1294, a DNA element that transposes by RC transposition. Plasmid. 2000;44(1):66–84. Epub 2000/06/30. pmid:10873528.
  114. 114. Hayashi W, Ohsaki Y, Taniguchi Y, Koide S, Kawamura K, Suzuki M, et al. High prevalence of blaCTX-M-14 among genetically diverse Escherichia coli recovered from retail raw chicken meat portions in Japan. Int J Food Microbiol. 2018;284:98–104. Epub 2018/08/11. pmid:30096596.
  115. 115. Manges AR, Johnson JR. Reservoirs of extraintestinal pathogenic Escherichia coli. Microbiol Spectr. 2015;3(5). Epub 2015/11/07. pmid:26542041.
  116. 116. Riley LW. Pandemic lineages of extraintestinal pathogenic Escherichia coli. Clin Microbiol Infect. 2014;20(5):380–90. Epub 2014/04/29. pmid:24766445.
  117. 117. Iida S, Mollet B, Meyer J, Arber W. Functional characterization of the prokaryotic mobile genetic element IS26. Mol Gen Genet. 1984;198(2):84–9. Epub 1984/01/01. pmid:6097800.
  118. 118. Mollet B, Iida S, Shepherd J, Arber W. Nucleotide sequence of IS26, a new prokaryotic mobile genetic element. Nucleic Acids Res. 1983;11(18):6319–30. Epub 1983/09/24. pmid:6312419; PubMed Central PMCID: PMC326375.
  119. 119. Sidjabat HE, Seah KY, Coleman L, Sartor A, Derrington P, Heney C, et al. Expansive spread of IncI1 plasmids carrying blaCMY-2 amongst Escherichia coli. Int J Antimicrob Agents. 2014;44(3):203–8. Epub 2014/07/24. pmid:25052868.
  120. 120. Harada K, Okada E, Shimizu T, Kataoka Y, Sawada T, Takahashi T. Antimicrobial resistance, virulence profiles, and phylogenetic groups of fecal Escherichia coli isolates: a comparative analysis between dogs and their owners in Japan. Comp Immunol Microbiol Infect Dis. 2012;35(2):139–44. Epub 2012/01/19. pmid:22252047.
  121. 121. Russo TA, Johnson JR. Medical and economic impact of extraintestinal infections due to Escherichia coli: focus on an increasingly important endemic problem. Microbes Infect. 2003;5(5):449–56. Epub 2003/05/10. pmid:12738001.
  122. 122. Johnson JR, Murray AC, Gajewski A, Sullivan M, Snippes P, Kuskowski MA, et al. Isolation and molecular characterization of nalidixic acid-resistant extraintestinal pathogenic Escherichia coli from retail chicken products. Antimicrob Agents Chemother. 2003;47(7):2161–8. Epub 2003/06/25. pmid:12821463; PubMed Central PMCID: PMC161843.
  123. 123. Belanger L, Garenaux A, Harel J, Boulianne M, Nadeau E, Dozois CM. Escherichia coli from animal reservoirs as a potential source of human extraintestinal pathogenic E. coli. FEMS Immunol Med Microbiol. 2011;62(1):1–10. Epub 2011/03/03. pmid:21362060.
  124. 124. Pomba C, Rantala M, Greko C, Baptiste KE, Catry B, van Duijkeren E, et al. Public health risk of antimicrobial resistance transfer from companion animals. J Antimicrob Chemother. 2017;72(4):957–68. Epub 2016/12/22. pmid:27999066.