Next Article in Journal
Comparative Performance of Line Probe Assay and GeneXpert in the Detection of Rifampicin Monoresistance in a TB-Endemic African Country
Next Article in Special Issue
Concordance between Antimicrobial Resistance Phenotype and Genotype of Staphylococcus pseudintermedius from Healthy Dogs
Previous Article in Journal
Antifungal Activity of a Library of Aminothioxanthones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Canine Saliva as a Possible Source of Antimicrobial Resistance Genes

1
Centre for Bioinformatics, University of Veterinary Medicine, 1078 Budapest, Hungary
2
Department of Operative Tecniques and Surgical Research, Faculty of Medicine, University of Debrecen, 4032 Debrecen, Hungary
3
Department of Thoracic Surgery, Borsod-Abaúj-Zemplén County Hospital and University Teaching Hospital, 3526 Miskolc, Hungary
4
Department of Small Animal Clinical Sciences, Western College of Veterinary Medicine, Saskatoon, SK S7N 5B4, Canada
5
Department of Surgery, Paracelsus Medical University, 90419 Nuremberg, Germany
6
Department of Surgery, Transplantation and Gastroenterology, Semmelweis University, 1085 Budapest, Hungary
7
Interdisciplinary Gastroenterology Working Group, Semmelweis University, 1085 Budapest, Hungary
8
Department and Clinic of Surgery and Ophthalmology, University of Veterinary Medicine, 1078 Budapest, Hungary
9
Department of Agricultural Biotechnology, Tekirdag Namik Kemal University, Tekirdag 59030, Turkey
10
Department of Microbiology and Infectious Diseases, University of Veterinary Medicine, 1143 Budapest, Hungary
11
Department of Phyisics of Complex Systems, Eötvös Loránd University, 1117 Budapest, Hungary
*
Author to whom correspondence should be addressed.
Antibiotics 2022, 11(11), 1490; https://doi.org/10.3390/antibiotics11111490
Submission received: 30 September 2022 / Revised: 23 October 2022 / Accepted: 24 October 2022 / Published: 27 October 2022

Abstract

:
While the One Health issues of intensive animal farming are commonly discussed, keeping companion animals is less associated with the interspecies headway of antimicrobial resistance. With the constant advance in veterinary standards, antibiotics are regularly applied in companion animal medicine. Due to the close coexistence of dogs and humans, dog bites and other casual encounters with dog saliva (e.g., licking the owner) are common. According to our metagenome study, based on 26 new generation sequencing canine saliva datasets from 2020 and 2021 reposited in NCBI SRA by The 10,000 Dog Genome Consortium and the Broad Institute within Darwin’s Ark project, canine saliva is rich in bacteria with predictably transferable antimicrobial resistance genes (ARGs). In the genome of potentially pathogenic Bacteroides, Capnocytophaga, Corynebacterium, Fusobacterium, Pasteurella, Porphyromonas, Staphylococcus and Streptococcus species, which are some of the most relevant bacteria in dog bite infections, ARGs against aminoglycosides, carbapenems, cephalosporins, glycylcyclines, lincosamides, macrolides, oxazolidinone, penams, phenicols, pleuromutilins, streptogramins, sulfonamides and tetracyclines could be identified. Several ARGs, including ones against amoxicillin–clavulanate, the most commonly applied antimicrobial agent for dog bites, were predicted to be potentially transferable based on their association with mobile genetic elements (e.g., plasmids, prophages and integrated mobile genetic elements). According to our findings, canine saliva may be a source of transfer for ARG-rich bacteria that can either colonize the human body or transport ARGs to the host bacteriota, and thus can be considered as a risk in the spread of antimicrobial resistance.

1. Introduction

Antimicrobial resistance (AMR) is a threat of utmost significance that constantly raises medical challenges all around the globe. It is a fact that the development of AMR and the increase in its prevalence are due to the use or misuse of antimicrobials. Moreover, according to the concept of One Health, the effects of antimicrobial use in human, animal and environmental sectors are interconnected, and thus interdependent. In the case of AMR, the relatedness of human and animal dimensions described by the One Health approach can be the best referred to by considering that most antimicrobial classes are co-used in both sectors. Human health antimicrobial use has been overshadowed for years by farm animal mass medication, although this tendency has recently changed in some parts of the world [1]. While the appearance and advance of AMR and, as an underlying cause, the enrichment and transmission of antimicrobial resistance genes (ARGs) in antibiotic-dense environments such as intensive animal production farms, is a well-examined phenomenon [2,3], the spread of AMR may also derive from other animal-borne routes.
Over the past decades, the number of companion animals has been tendentially and steadily rising [4]. Between 2000 and 2017, the number of dogs in the United States escalated from 68 million to 89.7 million [5]. In total, 67.9% of all households in the U.S. were associated with the ownership of various pet species and 48% of all with dogs in 2016 [4]. In the years 2019–2020, 50% of the U.S. population owned a dog [6]. The coronavirus disease (COVID-19) pandemic outbreak has resulted in elevated companion animal acquisition rates, albeit often followed by retention or replacement [7,8]. In addition to the popularity of keeping small animals, the quality of human–pet bonds has also changed. According to the survey of the American Veterinary Medical Association, 70% of pet owners consider their pets as family members, 17% as companions and 3% as property [9]. The role of pets can principally be defined as social companionship. Nowadays, having physical proximity is very common for pet–owner co-existences; pets often sleep together with their owners and lick their face or wounds [10]. Unfortunately and unsurprisingly, with such high dog numbers, the occurrence of dog bites is also common. Between 2001 and 2003, approximately 4.5 million dog bites were registered yearly in the United States, 19% of which necessitated medical intervention [11]. In the years 2005–2013, an average of 337,103 dog bite injuries were treated at U.S. emergency departments [12], although dog bites in general are under-reported [13]. In English hospitals, the number of dog bite admissions rose from 6.34 (95%CI 6.12–6.56) in 1998 to 14.99 (95%CI 14.67–15.31) admissions per 100,000 inhabitants in 2018 [14]. In parallel, a study on the database of The International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10) reveals that dog attack fatalities (excluding deaths caused by postattack infectious complications) have been constantly rising between 1995 and 2016. Incidence rates reached 0.009 per 100 000 inhabitants in Europe, 0.011 in the U.S., 0.007 in Canada and 0.004 in Australia [15]. Interestingly, 3 of 5 bites are executed by family dogs, which is more common than attacks by strays [16].
The modern mindset of providing regular veterinary healthcare services to our pets and keeping them in our closest surroundings may contribute to the interspecies transmission of AMR. Several studies have already turned their attention to the role of companion animals in the headway of AMR [17,18,19,20,21,22]. Nevertheless, the significance of the direct pet-borne AMR spread route has been given less attention when compared with the rather indirect, mostly food-transmitted farm-animal-associated route. The routes of ARG transmission can be assessed by analyzing the genetic surroundings of ARGs. Genetic elements that facilitate horizontal gene transfer (HGT) including plasmids, bacteriophages and integrative mobile genetic elements (iMGEs) may contribute to three different HGT processes, namely conjugation, transduction or transformation. By conjugation, genes are delivered to a recipient cell on plasmids if bacterial cells are physically binded, while by transduction, direct cell-to-cell contact is disclaimed due to the presence of bacteriophages as means of gene conduit. Transformation negates the need for particular delivery processes. In this case, bacteria take up genetic fragments from their environment [23]. After dog bites or close encounters with saliva from dogs that undergo veterinary treatments, and thus may carry bacteria with a possibly enriched ARG content, resistant bacteria may be introduced to the human body, and later the HGT of antimicrobial resistance determinants may be exchanged with the host bacteriota. Our study aimed to reveal the ARG content of canine saliva samples, attach the ARGs with the bacterial species that they derive from and report the ARGs’ spreading capabilities to weigh the above-mentioned phenomenon. For this purpose, freely accessible next-generation sequencing (NGS) shotgun datasets were downloaded and bioinformatically analyzed using an advanced metagenomic pipeline.

2. Materials and Methods

Deep-sequenced canine saliva datasets were searched in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) repository. In May 2021, two shotgun metagenomic BioProjects (PRJNA648123 [24]—The 10,000 Dog Genome Consortium and PRJNA683923 [25]—Broad Institute, Darwin’s Ark project) with more than 100,000,000 paired-end reads per sample were identified (Table 1). Both projects collected and sequenced dog saliva samples to investigate polymorphisms in the dog genome from which the samples were derived. The median read count (interquartile range, IQR) of the samples was 177.7 × 10 6 ( 26.6 × 10 6 ) and 417.7 × 10 6 ( 90.1 × 10 6 ) in datasets PRJNA648123 and PRJNA683923, respectively.

Bioinformatic Analysis

Quality-based filtering and trimming of the raw short reads was performed by TrimGalore (v.0.6.6, https://github.com/FelixKrueger/TrimGalore, accessed on 24 September 2022), setting 20 as a quality threshold. Only reads longer than 50 bp were retained and taxonomically classified using Kraken2 (v2.1.1) [26] and a database created (24 March 2021) from the NCBI RefSeq complete archaeal, bacterial, viral and plant genomes. For this taxon assignment, the–confidence 0.5 parameter was used to obtain more precise species-level hits. The taxon classification data were managed in R [27] using functions of the packages phyloseq (v1.36.0) [28] and microbiome (v1.14.0) [29]. Reads classified as origin of bacteria were assembled to contigs by MEGAHIT (v1.2.9) [30] using default settings. The contigs were also classified taxonomically by Kraken2 with the same database as above. All possible open reading frames (ORFs) were gathered by Prodigal (v2.6.3) [31] from the contigs. The protein-translated ORFs were aligned to the ARG sequences of the Comprehensive Antibiotic Resistance Database (CARD, v.3.1.3) [32,33] by Resistance Gene Identifier (RGI, v5.2.0) with Diamond [34]. The ORFs classified as perfect or strict were further filtered with 90% identity and 90% coverage. All nudged hits were excluded. The iMGE content of the ARG-harboring contigs was analyzed by MobileElementFinder (v1.0.3) and its database (v1.0.2) [35]. Following the distance concept of Johansson et al. [35] for each bacterial species, only those with a distance threshold defined within iMGEs and ARGs were considered associated. In the MobileElementFinder database (v1.0.2) for Bacteroides, the longest composite transposon (cTn) was the Tn6186. In the case of this genus, its length (8505 bp) was taken as the cut-off value. For the genera Enterococcus and Klebsiella, Tn6246 (5147 bp) and Tn125 (10,098 bp) provided the thresholds, respectively. In the case of Escherichia coli, this limit was the length of the Tn1681 transposon, namely 24,488 bp, while for Pseudomonas aeruginosa Tn6060 (25,440 bp). As the database neither contains species-level, nor genus-level cTn data for the rest of the species, a general cut-off value was chosen for the contigs of these species. This value was declared as the median of the longest cTns per species in the database (10,098 bp). The plasmid origin probability of the contigs was estimated by PlasFlow (v.1.1) [36]. The prophage content of the assembled contigs was prediced by VirSorter2 (v2.2.3) [37]. The findings were filtered for dsDNAphages and ssDNAs. All data management procedures, analyses and plottings were performed in R environment (v4.1.0) [27].

3. Results

After presenting the bacteriome and the identified AGRs (resistome), predictions regarding the mobility potential of ARGs were also resumed based on genetic characteristics that may play a significant role in HGT.

3.1. Bacteriome

By taxon classification, the number of reads aligning to bacterial genomes differed in the various samples. In the saliva, median bacterial read count of the samples was 4.3 × 10 6 (IQR: 3.4 × 10 6 ). A total of 16 major bacterial genera were detected within the saliva samples, out of which several aerobic and anaerobic genera often become isolated from infected dog bite wounds. The relative abundances of genera that achieved more than 1% of the bacterial hits in any of the saliva samples is shown in Figure 1. In the saliva samples, the dominant genera (with mean prevalence) in descending order were Porphyromonas (49%), Prevotella (15%), Pasteurella (12%), Neisseria (10%), Capnocytophaga (9%), Conchiformibius (7%), Frederiksenia (7%), Cutibacterium (6%), Actinomyces (5%), Campylobacter (4%), Desulfomicrobium (4%), Bacteroides (3%), Fusobacterium (3%), Mycoplasmopsis (3%), Treponema (3%) and Streptococcus (2%). In the sample No. 20, no reads were classified to bacteria.

3.2. Resistome

Applying the above-mentioned filtering conditions, we identified 69 ARGs that are presented together with their prevalence rates within the samples and the drug classes that these ARGs affect adversely in Table 2. These ARGs per sample with their coverage and the sequence identity rate are shown in Figure 2. As a result of the taxon classification on the contigs harboring the ARGs, it was possible to predict the bacterial species of putative origin for all but ten genes (Table 3). The dominant mechanisms of identified ARGs were antibiotic inactivation (47.69%), antibiotic target protection (23.41%), antibiotic target alteration (15.90%), antibiotic efflux (7.80%) and antibiotic target replacement (5.20%).

3.3. Mobilome

Many of the identified ARGs are harbored by iMGEs, prophages or plasmids. The frequencies of iMGEs, prophages and plasmids associated with ARGs by bacteria of origin are summarized in Figure 3. Some genes could have been attached to two of the above-mentioned mobility groups in the genome of one species, including the iMGE and prophage co-appearance of aminoglycoside resistance encoding aad(6) and aph(3’)-IIIa in E. faecium; the iMGE and plasmid co-appearance of aph(3’)-Ia in Corynebacterium sp. 1959 and K. quasipneumoniae; aph(3’)-Ia and aph(6)-Id in Variovorax sp. SRS16; aph(3”)-Ib in Variovorax sp. PAMC28562; tetracycline resistance encoding tetM in E. faecalis; and prophage and plasmid co-appearance of macrolide, lincosamide and streptogramin resistance encoding ermB in Enterococcus sp. FDAARGOS_375. The blaOXA-2, blaOXA-347 and blaTEM-116 genes associated with amoxicillin–clavulanate resistance all appeared in plasmids in various species; moreover, blaOXA-2 was associated with both an iMGE and a plasmid in the genome of P. aeruginosa.

4. Discussion

During the bacteriome, resistome and mobilome analysis of the canine saliva samples, a large set of results was obtained that can be examined from a One Health point of view, merging the small animal veterinary sector with the perspective of the human healthcare system.
ARGs were identified in all but one sample (No. 20). No reads of bacterial origin were found in this sample. We speculate that this may be since only those reads generated from sequencing were uploaded to the SRA mapped to the dog genome.
Many of the detected aerobic and anaerobic bacterial genera often become isolated from infected dog bite wounds. Dog bite infections are normally polymicrobial, and the bite wound bacteriota consist of bacteria from the animals’ oral cavity, the recipients’ skin and the environment. The most common pathogens in dog bites are Pasteurella spp. (P. canis and P. multocida), Staphylococcus spp., Streptococcus spp. and Capnocytophaga spp., Porphyromonas spp., Bacteroides spp., Fusobacterium spp. and Corynebacterium spp. [38], which all appeared in the analyzed saliva samples. Some other bacterial groups of a relatively higher clinical significance that were detected in the saliva samples, including Enterococcus spp., Moraxella spp., Neisseria spp., Prevotella spp. and Pseudomonas spp., are also often isolated from dog bite wound infections. The vast majority of other genera isolated in the samples have been mentioned to appear in dog saliva in previous publications with variable abundance rates [13,39]. Even though membersof Clostridium spp. were detected in the samples, genome fragments of C. tetani, the bacterium responsible for tetanus, were not identified. The number of detected ARGs was relatively high in the salivary bacteriome. Examining eight genera (Pasteurella spp., Staphylococcus spp., Streptococcus spp., Capnocytophaga spp., Porphyromonas spp., Bacteroides spp., Fusobacterium spp. and Corynebacterium spp.) that were indicated to be the most relevant ones in dog bite infections by other authors [38,39], we could identify genes that confer resistance against aminoglycosides, carbapenems, cephalosporins, glycylcyclines, lincosamides, macrolides, oxazolidinone, penams, phenicols, pleuromutilins, streptogramins, sulfonamides and tetracyclines, while other antimicrobial groups including fluoroquinolones appeared in E. coli, one of the six leading pathogens responsible for the deaths associated with resistance in 2019 worldwide [40].
Such a great number and broad spectrum of ARGs and potentially affected antimicrobial groups associated with the canine saliva samples may be related to the use of antibiotics at small-animal veterinary practices. Antibiotic consumption rates in the companion animal sector are rather difficult to evaluate. However, some systems exist for the surveillance of magnitude of companion animal antibiotic consumption, such as the European Surveillance of Veterinary Antimicrobial Consumption (ESVAC) [41], VetCompass [42] or the Small Animal Veterinary Surveillance Network (SAVSNET) [43], these rates are still less well-documented. Moreover, in many countries, antimicrobial use is often just estimated by rough sales data [21]. Nevertheless, according to the two UK-based surveillance systems (VetCompass from Royal Veterinary College, and SAVSNET from Liverpool University) and one EU report (ESVAC), antibiotics are rather frequently prescribed at small-animal clinics. A study states 1 in 4 UK dogs (25.2%, 95% CI: 25.1–25.3%) were treated with antibiotics in a two-year period [44]. Even though the vast majority of veterinarians are aware of the fact that improper AMU contributes to selection for AMR, and that it is a significant problem according to nationwide surveys [45,46], there are many factors that influence the antibiotic prescription preferences of veterinarians in addition to the perspectives of antimicrobial stewardship [47,48,49,50,51,52].
Broad-spectrum amoxicillin–clavulanate is the flagship of antimicrobial agents applied in dogs in many countries, while first-generation cephalosporins are also routinely used [21,53,54]. Lincosamides (clindamycin), macrolides, tetracyclines (doxycycline), nitroimidazoles and trimethoprim/sulphonamides have also been reported to be frequently used in small-animal practices [21]. Third- and fourth-generation cephalosoprines, fluoroquinolones and polymixins that belong to cathegory B, ‘last resort’, or highest-priority Critically Important Antibiotics (HPCIAs) according to the European Medicines Agency [41] should be avoided unless sensitivity testing is conducted and no other antibiotics would be effective. Nevertheless, HPCIAs have been estimated to be prescribed in around 5–6% of total antimicrobial agent usage events. Of the HPCIA category, fluoroquinolones are the most common in dogs, constituting ∼4 to 5% of total antibiotic prescriptions [55].
Neverheless, AMR determinants against the above-mentioned antimicrobial compound groups have been detected in and associated with many bacterial species of the examined canine saliva samples. Some of these ARG-associated bacteria can also exert pathogenic effects and are often isolated by dog bite infections. ARGs against cephalosporins were identified in many, often clinically significant bacteria, including but not limited to Bacteroides spp., Capnocytophaga spp., E. coli, F. ulcerans, Porphyromonas spp. and P. aeruginosa. Likewise, ARGs against lincosamides appeared in Bacteroides spp., C. stomatis, Enterococcus spp., P. cangingivalis, P. intermedia and Streptococcus spp.; ARGs against macrolides in Bacteroides spp., C. stomatis, Enterococcus spp., Porphyromonas spp., P. intermedia and Streptococcus spp.; ARGs against tetracyclines in Bacteroides spp., C. stomatis, Enterococcus spp., P. crevioricanis, Prevotella spp., P. putida and Streptococcus spp.; ARGs against sulfonamides in Corynebacterium spp., M. bovis, P. multocida and P. aeruginosa; and ARGs against fluoroquinolones in E. coli (non-exhaustive lists of detected ARG-bacteria associations in the samples, with an emphasis on clinically relevant bacteria). Bacterial associations and the significance of AMR determinants affecting amoxicillin–clavulanate, the most commonly prescribed antibiotic in veterinary medicine, are discussed later on. Even though nitroimidazoles (e.g., metronidazole) are described to be often used in small-animal veterinary practices, no ARGs could be detected in the canine saliva samples against this antibiotic group.
In the current literature, human infections associated to dog bites are better and more frequently documented than the transmission route of licking. Three to thirty percent of dog bites lead to infection [13]. The management of animal bites rests on two pillars: local wound care and adequately applied systemic treatment. Essentials of local therapy include inspection, debridement of the wound accompanied by the removal of possible foreign bodies, e.g., teeth, and irrigation with saline solution. As for the systematic therapy, tetanus boosters (if none given in the past year) and rabies prophylaxis should always be considered. In our study, genome fragments of C. tetani, the causative agent of tetanus, were not detected in any of the examined saliva samples. No consensus has yet been found in the use of antibiotics for animal bite wound care. Prophylactic antibiotics should be considered unless the wound is very superficial and clean. Explicit indications for antibiotic prophylaxis or therapy include presentation at least 8 h after the bite, clear signs of superinfection, moderate or severe wounds with crush injuries or devitalized tissues requiring surgery, deep puncture wounds (exceeding the layer of epidermis), wounds close to joints, diabetes mellitus, asplenic or immunocompromised state, alcohol abuse, or involvement of the genital area, face or hand [56,57,58,59,60]. In the absence of the above reasons, antibiotic therapy may not be necessary. Interestingly, injuries are normally located on the head, neck and face in children and on the hand or upper extremity in adults due to height ratios with the attacking dog [13,61]. An adequately chosen antibiotic agent is expected to be effective against anaerobe bacteria (Bacteroides spp., Fusobacterium spp., Porphyromonas spp., Prevotella spp. etc.), in addition to the Staphylococcus, Streptococcus and Pasterurella species. Prophylactic treatment is normally 3 to 5 days long, while medication for 10 days or longer is recommended if the wound is infected. The first-line choice for oral therapy is amoxicilin–clavulanate, accompanied with a first dose of intravenous antibiotic (e.g., ampicillin/sulbactam, ticarcillin–clavulanate, piperacillin–tazobactam, or a carbapenem) in high-risk patients. Amoxicillin–clavunalate is often combined with metronidazole or clindamycin and is also sometimes replaced with cephalosporins, e.g., cefuroxime, cefotaxime, ceftriaxone or amoxicillin, fluoroquinolones, sulfamethoxazole and trimethoprim, and alternatively, although less effective, azithromycin or doxycycline in this combination [58,60]. Due to high resistance rates, flucloxacillin, erythromycin and cephalosporins are often ineffective in Pasteurella infections, and thus should rather be avoided [57]. In our case, no genes conferring resistance to these agent groups could be identified in Pasturella spp.
Data on the outcome of antibiotic prophylaxis in animal bite management by humans is limited and rather controversial and conflicting. While a meta-analysis of eight randomized trials indicated a benefit of antibiotic prophylaxis [62], some studies concluded that antibiotic prophylaxis does not result in a statistically significant difference in the frequency of wound infections among treated and untreated patient groups, except for wounds to the hand [63]. Based on other publications, antibiotic prophylaxis should be recommended for high-risk patient groups only [64,65].
Based on antibiotic prescription data from human and veterinary medical practices described above, amoxicillin–clavulanate and cephalosporins are the most commonly used agents in the treatment of animal patients and dog bite infections, while lincosamides (mostly clindamycin), sulfonamides (mostly potentiated sulfonamides) and fluoroquinolones also appear in both sectors [21,53,54,55,57,58,61].
Amoxicillin–clavulanate, the most commonly used antibiotic in small-animal medicine, and the first choice for canine bite wounds, is a member of broad-spectrum penicillins that have been a frequently consumed key antibiotic group in the high-income super-region between 2000 and 2018 according to a global study [66]. All in all, six ARG types were detected in the dog saliva samples that may confer resistance against amoxicillin–clavulanate, which were either the members of the blaTEM or OXA family [67,68]. blaTEM-116 was identified in E. coli, while various members of the OXA family appeared in in many genera, including A. baumannii, Bacteroides spp., Capnocytophaga spp., F. ulcerans and Pseudomonas spp., which can have a high clinical relevance in dog bite infections. Moreover, plasmid-borne blaOXA-2, blaOXA-347 and blaTEM-116 genes in different bacteria, and OXA-2 associated with both an iMGE and a plasmid from P. aeruginosa, all confer resistance against amoxicillin–clavulanate and have a higher potential to spread from bacteria to bacteria. The accumulation of various mobility factors around the genes may increase the chance of the horizontal transfer of the given ARG. The canine saliva-borne transmission of bacteria harboring mobile ARGs may hamper antimicrobial use in human clinical settings and can also contribute to the spread of AMR among the bacteria derived from pets to the bacteriota appearing in humans.
Cephalosporins, which are also commonly used both in companion animals medicine [21,54] and human medicine, including in protocols for dog bite infections [58,61], have been associated with all in all twenty-four bacterial species harboring related ARGs in the canine saliva samples. These bacteria include pathogens such as A. baumanni, Bacteroides spp., Capnocytophaga spp., F. ulcerans and P. aeruginosa that contain blaOXA-2, blaOXA-347, blaOXA-347, blaOXA-85 and blaOXA-2, respectively. Importantly, P. aeruginosa-associated blaOXA-2, which was associated with both an iMGE and a plasmid, has a high potential for HGT and can be considered as a major public health concern.
Lincosamides, including clindamycin, are also significant in both veterinary medicine [21,54] and the treatment of dog bite cases in humans [58,61]. Genes affecting lincosamides appeared in seventeen bacterial species, many of which, e.g., B. fragilis (ErmF), C. stomatis (ErmF), Porphyromonas spp. (ErmF) and Streptococcus spp. (ErmB, lnuB and lnuC), are potentially pathogenic. While the above-mentioned species, which are often associated with dog bite infections, contained no genetic element around these ARGs that could facilitate their transfer, an Enterococcus species from the canine saliva samples is linked to carrying ermB with prophage and plasmid co-appearance. Thus, the possibility of the transfer of this gene to other bacteria with higher clinical significance is given in the case of lincosamides as well.
Sulfonamide resistance genes appeared in relatively fewer, namely eight, bacterial species, including P. multocida and P. aeruginosa harboring sul2 and sul1, respectively. In addition to the low ARG counts, none of the sulfonamide resistance genes appeared to be mobile, while, in contrast, HGT has been found to be highly characteristic for these two genes in other publications [69].
In the case of fluoroquinolones, only E. coli harbored related ARGs, namely the non-mobile acrA, gadW and gadX. Interestingly, while both gadW and gadX are AraC family regulators that promote mdtEF expression to confer multidrug resistance, when they co-occur, gadW inhibits gadX-dependent activation by repressing gadX (CARD).
Furthermore, fosfomycin and tygecyclin, which are often used as last-resort antimicrobial agents [70,71] and are involved in the list of Critically Important Antibiotics for Human Medicine by WHO [72], also appeared among the affected antibiotic groups, namely due to the presence of fosA2 and tetX4, tetX5, respectively. Nonetheless, while the members of the tet(X) family are often plasmid-associated [73,74,75,76], their plasmid relatedness was not predicted within this study.
While the above-mentioned findings may raise awareness of the potential public health risk associated with canine saliva, this material has been used to promote rapid healing and to reduce bacterial contamination in the past according to the reports of ethnoveterinary and ethnomedicinal practices [77,78]. Antimicrobial and anti-imflammatory activity of canine saliva induced by thiocyanate, lysozyme and, indirectly, nitrate, among others [79,80], can even appear at low concentrations [81]. However, according to our findings, canine saliva can also be associated with public health risks, since salival bacteria may contaminate the surroundings of people and may also colonize human skin and mucous membranes. Thus, ARG-rich bacteria present in and around humans do not even necessarily need to transfer their ARGs to potentially cause severe harm to various groups of people with weaknesses of the immune system, e.g., extremities in age or diseased state.
As a common trend among many nations, veterinary use of antibiotics is gradually declining [49,55,81,82,83]. In human medicine, antibiotic sales elevated by 65% in low- and middle-income countries and decreased slightly by 4% in high-income countries between 2000 and 2015, which adds up to a rise in global antibiotic consumption rates [66,84].As a presumable conclusion, several genes conferring resistance against clinically important antimicrobial groups are present in the salivary bacteriome of dogs that may drift to the genome of bacteria in humans. Encounters with dog saliva and dog bites may serve as an interspecies platform for the migration of bacteria and ARGs. Transmitted bacteria may cause clinical symptoms, and ARGs that they harbor may confer resistance against antimicrobial agents of a clinical relevance.

Author Contributions

N.S. takes responsibility for the integrity of the data and the accuracy of the data analysis. A.G.T., L.M. and N.S. conceived the concept of the study. E.G.K. and N.S. participated in the sample selection. A.G.T. and N.S. participated in the bio-informatic analysis. A.G.T., B.R., I.T. and N.S. participated in the drafting of the manuscript. A.D., A.G.T., Á.V.P., B.R., I.T., L.M., N.S., S.K. and T.N. carried out the critical revision of the manuscript for important intellectual content. All authors read and approved the final manuscript.

Funding

The research was supported by the European Union’s Horizon 2020 research and innovation program under Grant Agreement No. 874735 (VEO).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The short-read data of the sample data are publicly available and can be accessed through the PRJNA648123 [24] and PRJNA683923 [25] from the NCBI Sequence Read Archive (SRA).

Acknowledgments

The authors would like to thank the providers of BioProject PRJNA648123 [24] and, furthermore, Kathleen Morrill; Elinor Karlsson; Kerstin Lindblad-Toh, (Karlsson Laboratory, Broad Institute) and the Darwin’s Ark project (darwinsark.org, accessed on 24 September 2022) for providing the datasets and metadata for BioProject PRJNA683923 [25].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Centre for Disease Prevention and Control (ECDC); European Food Safety Authority (EFSA); European Medicines Agency (EMA). Third joint inter-agency report on integrated analysis of consumption of antimicrobial agents and occurrence of antimicrobial resistance in bacteria from humans and food-producing animals in the EU/EEA. EFSA J. 2021, 19, e06712. [Google Scholar] [CrossRef]
  2. Manyi-Loh, C.; Mamphweli, S.; Meyer, E.; Okoh, A. Antibiotic use in agriculture and its consequential resistance in environmental sources: Potential public health implications. Molecules 2018, 23, 795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tóth, A.G.; Csabai, I.; Krikó, E.; Tőzsér, D.; Maróti, G.; Patai, Á.V.; Makrai, L.; Szita, G.; Solymosi, N. Antimicrobial resistance genes in raw milk for human consumption. Sci. Rep. 2020, 10, 7464. [Google Scholar] [CrossRef] [PubMed]
  4. Rowan, A.N. Companion Animal Statistics in the USA. 2018. Available online: https://animalstudiesrepository.org/cgi/viewcontent.cgi?article=1002&context=demscapop (accessed on 6 August 2022).
  5. Bedford, E. Number of Dogs in the United States from 2000 to 2017. 2019. Available online: https://www.statista.com/statistics/198100/dogs-in-the-united-states-since-2000 (accessed on 6 March 2022).
  6. American Pet Products Association. 2021–2022 APPA National Pet Owners Survey. 2022. Available online: https://www.americanpetproducts.org/puBs_survey.asp (accessed on 6 March 2022).
  7. Ho, J.; Hussain, S.; Sparagano, O. Did the COVID-19 pandemic spark a public interest in pet adoption? Front. Vet. Sci. 2021, 8, 647308. [Google Scholar] [CrossRef] [PubMed]
  8. Hoffman, C.L.; Thibault, M.; Hong, J. Characterizing pet acquisition and retention during the COVID-19 pandemic. Front. Vet. Sci. 2021, 8, 781403. [Google Scholar] [CrossRef] [PubMed]
  9. Jensen, H. AVMA Pet Ownership and Demographics Sourcebook, 2017–2018 ed.; American Veterinary Medical Association: Schaumburg, IL, USA, 2017. [Google Scholar]
  10. Overgaauw, P.A.; Vinke, C.M.; van Hagen, M.A.; Lipman, L.J. A one health perspective on the human–companion animal relationship with emphasis on zoonotic aspects. Int. J. Environ. Res. Public Health 2020, 17, 3789. [Google Scholar] [CrossRef] [PubMed]
  11. Gilchrist, J.; Sacks, J.J.; White, D.; Kresnow, M.j. Dog bites: Still a problem? Inj. Prev. 2008, 14, 296–301. [Google Scholar] [CrossRef]
  12. Loder, R.T. The demographics of dog bites in the United States. Heliyon 2019, 5, e01360. [Google Scholar] [CrossRef] [Green Version]
  13. Dhillon, J.; Hoopes, J.; Epp, T. Scoping decades of dog evidence: A scoping review of dog bite-related sequelae. Can. J. Public Health 2019, 110, 364–375. [Google Scholar] [CrossRef]
  14. Tulloch, J.S.; Owczarczak-Garstecka, S.C.; Fleming, K.M.; Vivancos, R.; Westgarth, C. English hospital episode data analysis (1998–2018) reveal that the rise in dog bite hospital admissions is driven by adult cases. Sci. Rep. 2021, 11, 1767. [Google Scholar] [CrossRef]
  15. Sarenbo, S.; Svensson, P.A. Bitten or struck by dog: A rising number of fatalities in Europe, 1995–2016. Forensic Sci. Int. 2021, 318, 110592. [Google Scholar] [CrossRef] [PubMed]
  16. Overall, K.L.; Love, M. Dog bites to humans—Demography, epidemiology, injury, and risk. J. Am. Vet. Med. Assoc. 2001, 218, 1923–1934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Guardabassi, L.; Schwarz, S.; Lloyd, D.H. Pet animals as reservoirs of antimicrobial-resistant bacteria. J. Antimicrob. Chemother. 2004, 54, 321–332. [Google Scholar] [CrossRef] [PubMed]
  18. Lloyd, D.H. Reservoirs of antimicrobial resistance in pet animals. Clin. Infect. Dis. 2007, 45, S148–S152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Pedersen, K.; Pedersen, K.; Jensen, H.; Finster, K.; Jensen, V.F.; Heuer, O.E. Occurrence of antimicrobial resistance in bacteria from diagnostic samples from dogs. J. Antimicrob. Chemother. 2007, 60, 775–781. [Google Scholar] [CrossRef] [Green Version]
  20. Damborg, P.; Sørensen, A.H.; Guardabassi, L. Monitoring of antimicrobial resistance in healthy dogs: First report of canine ampicillin-resistant Enterococcus Faecium Clonal Complex 17. Vet. Microbiol. 2008, 132, 190–196. [Google Scholar] [CrossRef]
  21. Pomba, C.; Rantala, M.; Greko, C.; Baptiste, K.E.; Catry, B.; Van Duijkeren, E.; Mateus, A.; Moreno, M.A.; Pyörälä, S.; Ružauskas, M.; et al. Public health risk of antimicrobial resistance transfer from companion animals. J. Antimicrob. Chemother. 2017, 72, 957–968. [Google Scholar]
  22. Li, Y.; Fernández, R.; Durán, I.; Molina-López, R.A.; Darwich, L. Antimicrobial resistance in bacteria isolated from cats and dogs from the Iberian Peninsula. Front. Microbiol. 2021, 11, 621597. [Google Scholar] [CrossRef] [PubMed]
  23. Goh, S. Clostridium Difficile: Methods and Protocols; Chapter Phage Transduction; Springer: New York, NY, USA, 2016; pp. 177–185. [Google Scholar]
  24. Ostrander, E.A.; Wang, G.D.; Larson, G.; Vonholdt, B.M.; Davis, B.W.; Jagannathan, V.; Hitte, C.; Wayne, R.K.; Zhang, Y.P. Dog10K: An international sequencing effort to advance studies of canine domestication, phenotypes and health. Natl. Sci. Rev. 2019, 6, 810–824. [Google Scholar] [CrossRef]
  25. Morrill, K.; Hekman, J.; Li, X.; McClure, J.; Logan, B.; Goodman, L.; Gao, M.; Dong, Y.; Alonso, M.; Carmichael, E.; et al. Ancestry-inclusive dog genomics challenges popular breed stereotypes. Science 2022, 376, eabk0639. [Google Scholar] [CrossRef]
  26. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [PubMed]
  27. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  28. McMurdie, P.J.; Holmes, S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Lahti, L.; Shetty, S. Microbiome R Package, 2012–2019. Available online: https://www.bioconductor.org/packages/release/bioc/html/microbiome.html (accessed on 29 September 2022).
  30. Li, D.; Liu, C.M.; Luo, R.; Sadakane, K.; Lam, T.W. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015, 31, 1674–1676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Hyatt, D.; Chen, G.L.; LoCascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010, 11, 119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. McArthur, A.G.; Waglechner, N.; Nizam, F.; Yan, A.; Azad, M.A.; Baylay, A.J.; Bhullar, K.; Canova, M.J.; De Pascale, G.; Ejim, L.; et al. The comprehensive antibiotic resistance database. Antimicrob. Agents Chemother. 2013, 57, 3348–3357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef] [PubMed]
  34. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods. 2015, 12, 59–60. [Google Scholar] [CrossRef]
  35. Johansson, M.H.; Bortolaia, V.; Tansirichaiya, S.; Aarestrup, F.M.; Roberts, A.P.; Petersen, T.N. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J. Antimicrob. Chemother. 2021, 76, 101–109. [Google Scholar] [CrossRef]
  36. Krawczyk, P.S.; Lipinski, L.; Dziembowski, A. PlasFlow: Predicting plasmid sequences in metagenomic data using genome signatures. Nucleic Acids Res. 2018, 46, e35. [Google Scholar] [CrossRef] [Green Version]
  37. Guo, J.; Bolduc, B.; Zayed, A.A.; Varsani, A.; Dominguez-Huerta, G.; Delmont, T.O.; Pratama, A.A.; Gazitúa, M.C.; Vik, D.; Sullivan, M.B.; et al. VirSorter2: A multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses. Microbiome 2021, 9, 37. [Google Scholar] [CrossRef]
  38. Thomas, N.; Brook, I. Animal bite-associated infections: Microbiology and treatment. Expert Rev. Anti. Infect. Ther. 2011, 9, 215–226. [Google Scholar] [CrossRef] [PubMed]
  39. Abrahamian, F.M.; Goldstein, E.J. Microbiology of animal bite wound infections. Clin. Microbiol. Rev. 2011, 24, 231–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Murray, C.J.; Ikuta, K.S.; Sharara, F.; Swetschinski, L.; Aguilar, G.R.; Gray, A.; Han, C.; Bisignano, C.; Rao, P.; Wool, E.; et al. Global burden of bacterial antimicrobial resistance in 2019: A systematic analysis. Lancet 2022, 399, 629–655. [Google Scholar] [CrossRef]
  41. European Medicines Agency. European Surveillance of Veterinary Antimicrobial Consumption. Available online: https://www.ema.europa.eu/en/veterinary-regulatory/overview/antimicrobial-resistance/european-surveillance-veterinary-antimicrobial-consumption-esvac (accessed on 6 March 2022).
  42. McGreevy, P.; Thomson, P.; Dhand, N.K.; Raubenheimer, D.; Masters, S.; Mansfield, C.S.; Baldwin, T.; Soares Magalhaes, R.J.; Rand, J.; Hill, P.; et al. VetCompass Australia: A national big data collection system for veterinary science. Animals 2017, 7, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Radford, A.; Noble, P.; Coyne, K.; Gaskell, R.; Jones, P.; Bryan, J.; Setzkorn, C.; Tierney, Á.; Dawson, S. Antibacterial prescribing patterns in small animal veterinary practice identified via SAVSNET: The small animal veterinary surveillance network. Vet. Rec. 2011, 169, 310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Buckland, E.L.; O’Neill, D.; Summers, J.; Mateus, A.; Church, D.; Redmond, L.; Brodbelt, D. Characterisation of antimicrobial usage in cats and dogs attending UK primary care companion animal veterinary practices. Vet. Rec. 2016, 179, 489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Odoi, A.; Samuels, R.; Carter, C.N.; Smith, J. Antibiotic prescription practices and opinions regarding antimicrobial resistance among veterinarians in Kentucky, USA. PLoS ONE 2021, 16, e0249653. [Google Scholar] [CrossRef]
  46. Sobierajski, T.; Mazińska, B.; Chajęcka-Wierzchowska, W.; Śmiałek, M.; Hryniewicz, W. Antimicrobial and Antibiotic Resistance from the Perspective of Polish Veterinary Students: An Inter-University Study. Antibiotics 2022, 11, 115. [Google Scholar] [CrossRef]
  47. Hughes, L.A.; Williams, N.; Clegg, P.; Callaby, R.; Nuttall, T.; Coyne, K.; Pinchbeck, G.; Dawson, S. Cross-sectional survey of antimicrobial prescribing patterns in UK small animal veterinary practice. Prev. Vet. Med. 2012, 104, 309–316. [Google Scholar] [CrossRef]
  48. Jacob, M.E.; Hoppin, J.A.; Steers, N.; Davis, J.L.; Davidson, G.; Hansen, B.; Lunn, K.F.; Murphy, K.M.; Papich, M.G. Opinions of clinical veterinarians at a US veterinary teaching hospital regarding antimicrobial use and antimicrobial-resistant infections. J. Am. Vet. Med. Assoc. 2015, 247, 938–944. [Google Scholar] [CrossRef]
  49. Hardefeldt, L.Y.; Selinger, J.; Stevenson, M.A.; Gilkerson, J.R.; Crabb, H.; Billman-Jacobe, H.; Thursky, K.; Bailey, K.E.; Awad, M.; Browning, G.F. Population wide assessment of antimicrobial use in dogs and cats using a novel data source—A cohort study using pet insurance data. Vet. Microbiol. 2018, 225, 34–39. [Google Scholar] [CrossRef] [PubMed]
  50. Hopman, N.E.; Mughini-Gras, L.; Speksnijder, D.C.; Wagenaar, J.A.; van Geijlswijk, I.M.; Broens, E.M. Attitudes and perceptions of Dutch companion animal veterinarians towards antimicrobial use and antimicrobial resistance. Prev. Vet. Med. 2019, 170, 104717. [Google Scholar] [CrossRef] [PubMed]
  51. Norris, J.M.; Zhuo, A.; Govendir, M.; Rowbotham, S.J.; Labbate, M.; Degeling, C.; Gilbert, G.L.; Dominey-Howes, D.; Ward, M.P. Factors influencing the behaviour and perceptions of Australian veterinarians towards antibiotic use and antimicrobial resistance. PLoS ONE 2019, 14, e0223534. [Google Scholar] [CrossRef]
  52. Samuels, R.; Qekwana, D.N.; Oguttu, J.W.; Odoi, A. Antibiotic prescription practices and attitudes towards the use of antimicrobials among veterinarians in the City of Tshwane, South Africa. PeerJ 2021, 9, e10144. [Google Scholar] [CrossRef] [PubMed]
  53. Kvaale, M.; Grave, K.; Kristoffersen, A.; Norström, M. The prescription rate of antibacterial agents in dogs in Norway–geographical patterns and trends during the period 2004–2008. J. Vet. Pharmacol. Ther. 2013, 36, 285–291. [Google Scholar] [CrossRef] [PubMed]
  54. Tompson, A.C.; Mateus, A.L.; Brodbelt, D.C.; Chandler, C.I. Understanding antibiotic use in companion animals: A literature review identifying avenues for future efforts. Front. Vet. Sci. 2021, 8, 719547. [Google Scholar] [CrossRef] [PubMed]
  55. Singleton, D.; Sánchez-Vizcaíno, F.; Dawson, S.; Jones, P.; Noble, P.; Pinchbeck, G.; Williams, N.; Radford, A. Patterns of antimicrobial agent prescription in a sentinel population of canine and feline veterinary practices in the United Kingdom. Vet. J. 2017, 224, 18–24. [Google Scholar] [CrossRef] [Green Version]
  56. Ellis, R.; Ellis, C. Dog and cat bites. Am. Fam. Physician 2014, 90, 239–243. [Google Scholar]
  57. Malahias, M.; Jordan, D.; Hughes, O.; Khan, W.S.; Hindocha, S. Bite injuries to the hand: Microbiology, virology and management. Open Orthop. J. 2014, 8, 157–161. [Google Scholar] [CrossRef] [Green Version]
  58. Aziz, H.; Rhee, P.; Pandit, V.; Tang, A.; Gries, L.; Joseph, B. The current concepts in management of animal (dog, cat, snake, scorpion) and human bite wounds. J. Trauma Acute Care Surg. 2015, 78, 641–648. [Google Scholar] [CrossRef]
  59. Jakeman, M.; Oxley, J.A.; Owczarczak-Garstecka, S.C.; Westgarth, C. Pet dog bites in children: Management and prevention. BMJ Paediatr. Open 2020, 4, e000726. [Google Scholar] [CrossRef] [PubMed]
  60. Zangari, A.; Cerigioni, E.; Nino, F.; Guidi, R.; Gulia, C.; Piergentili, R.; Ilari, M.; Mazzoni, N.; Cobellis, G. Dog bite injuries in a tertiary care children’s hospital: A seven-year review. Pediatr. Int. 2021, 63, 575–580. [Google Scholar] [CrossRef] [PubMed]
  61. Monroy, A.; Behar, P.; Nagy, M.; Poje, C.; Pizzuto, M.; Brodsky, L. Head and neck dog bites in children. Otolaryngol. Head Neck Surg. 2009, 140, 354–357. [Google Scholar] [CrossRef]
  62. Cummings, P. Antibiotics to prevent infection in pateints with dog bite wounds: A meta-analysis of randomized trials. Ann. Emerg. Med. 1994, 23, 535–540. [Google Scholar] [CrossRef]
  63. Medeiros, I.M.; Saconato, H. Antibiotic prophylaxis for mammalian bites. Cochrane Database Syst. Rev. 2001. [Google Scholar] [CrossRef] [PubMed]
  64. Stevens, D.L.; Bisno, A.L.; Chambers, H.F.; Everett, E.D.; Dellinger, P.; Goldstein, E.J.C.; Gorbach, S.L.; Hirschmann, J.V.; Kaplan, E.L.; Montoya, J.G.; et al. Practice Guidelines for the Diagnosis and Management of Skin and Soft-Tissue Infections. Clin. Infect. Dis. 2005, 41, 1373–1406. [Google Scholar] [CrossRef] [PubMed]
  65. Morgan, M.; Palmer, J. Dog bites. BMJ 2007, 334, 413–417. [Google Scholar] [CrossRef] [PubMed]
  66. Browne, A.J.; Chipeta, M.G.; Haines-Woodhouse, G.; Kumaran, E.P.; Hamadani, B.H.K.; Zaraa, S.; Henry, N.J.; Deshpande, A.; Reiner Jr, R.C.; Day, N.P.; et al. Global antibiotic consumption and usage in humans, 2000–18: A spatial modelling study. Lancet Planet. Health 2021, 5, e893–e904. [Google Scholar] [CrossRef]
  67. Di Conza, J.A.; Badaracco, A.; Ayala, J.; Rodriguez, C.; Famiglietti, Á.; Gutkind, G.O. β-lactamases produced by amoxicillin-clavulanate-resistant enterobacteria isolated in Buenos Aires, Argentina: A new blaTEM gene. Rev. Argent Microbiol. 2014, 46, 210–217. [Google Scholar] [CrossRef] [Green Version]
  68. Davies, T.J.; Stoesser, N.; Sheppard, A.E.; Abuoun, M.; Fowler, P.; Swann, J.; Quan, T.P.; Griffiths, D.; Vaughan, A.; Morgan, M.; et al. Reconciling the potentially irreconcilable? Genotypic and phenotypic amoxicillin-clavulanate resistance in Escherichia Coli. Antimicrob. Agents Chemother. 2020, 64, e02026-19. [Google Scholar] [CrossRef] [Green Version]
  69. Jiang, H.; Cheng, H.; Liang, Y.; Yu, S.; Yu, T.; Fang, J.; Zhu, C. Diverse mobile genetic elements and conjugal transferability of sulfonamide resistance genes (sul1, sul2, and sul3) in Escherichia coli isolates from Penaeus vannamei and pork from large markets in Zhejiang, China. Front. Microbiol. 2019, 10, 1787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Knafl, D.; Winhofer, Y.; Lötsch, F.; Weisshaar, S.; Steininger, C.; Burgmann, H.; Thalhammer, F. Tigecycline as last resort in severe refractory Clostridium difficile infection: A case report. J. Hosp. Infect. 2016, 92, 296–298. [Google Scholar] [CrossRef] [PubMed]
  71. Wong, M.; Wong, D.; Malhotra, S. Intravenous fosfomycin as salvage therapy for osteomyelitis caused by multidrug-resistant Pseudomonas aeruginosa. Am. J. Health-Syst. Pharm. 2021, 78, 2209–2215. [Google Scholar] [CrossRef] [PubMed]
  72. World Health Organization. Critically Important Antimicrobials for Human Medicine; World Health Organization: Geneva, Switzerland, 2019.
  73. Chen, C.; Cui, C.Y.; Wu, X.T.; Fang, L.X.; He, Q.; He, B.; Long, T.F.; Liao, X.P.; Chen, L.; Liu, Y.H.; et al. Spread of tet (X5) and tet (X6) genes in multidrug-resistant Acinetobacter baumannii strains of animal origin. Vet. Microbiol. 2021, 253, 108954. [Google Scholar] [CrossRef] [PubMed]
  74. Fu, Y.; Chen, Y.; Liu, D.; Yang, D.; Liu, Z.; Wang, Y.; Wang, J.; Wang, X.; Xu, X.; Li, X.; et al. Abundance of tigecycline resistance genes and association with antibiotic residues in Chinese livestock farms. J. Hazard. Mater. 2021, 409, 124921. [Google Scholar] [CrossRef]
  75. He, T.; Wang, R.; Liu, D.; Walsh, T.R.; Zhang, R.; Lv, Y.; Ke, Y.; Ji, Q.; Wei, R.; Liu, Z.; et al. Emergence of plasmid-mediated high-level tigecycline resistance genes in animals and humans. Nat. Microbiol. 2019, 4, 1450–1456. [Google Scholar] [CrossRef]
  76. Wang, L.; Liu, D.; Lv, Y.; Cui, L.; Li, Y.; Li, T.; Song, H.; Hao, Y.; Shen, J.; Wang, Y.; et al. Novel plasmid-mediated tet (X5) gene conferring resistance to tigecycline, eravacycline, and omadacycline in a clinical Acinetobacter baumannii isolate. Antimicrob. Agents Chemother. 2019, 64, e01326-19. [Google Scholar] [CrossRef]
  77. Verrier, L. Dog licks man. Lancet 1970, 295, 615. [Google Scholar] [CrossRef]
  78. Vallejo, J.R.; Santos-Fita, D.; González, J.A. The therapeutic use of the dog in Spain: A review from a historical and cross-cultural perspective of a change in the human-dog relationship. J. Ethnobiol. Ethnomed. 2017, 13, 47. [Google Scholar] [CrossRef] [Green Version]
  79. Benjamin, N.; Pattullo, S.; Weller, R.; Smith, L.; Ormerod, A. Wound licking and nitric oxide. Lancet 1997, 349, 1776. [Google Scholar] [CrossRef]
  80. Dhasarathan, P.; Arunkumar, R.; Blessy Jesubell, R.; Sowmya, P. Analysis of compounds and screening for anti-microbial, anti-inflammatory activity from saliva of Canis lupus Fam. Int. J. Ethnomed. Pharm. Res. 2013, 1, 1–6. [Google Scholar]
  81. Akpomie, O.; Ukoha, P.; Nwafor, O.; Umukoro, G. Saliva of different dog breeds as antimicrobial agents against microorganisms isolated from wound infections. Anim. Sci. J. 2011, 2, 18–22. [Google Scholar]
  82. Bager, F.; Bortolaia, V.; Ellis-Iversen, J.; Hendriksen, R.; Borck Høg, B.; Jensen, L.; Korsgaard, H.; Pedersen, K.; Dalby, T.; Træholt Franck, K.; et al. DANMAP 2016—Use of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Bacteria from Food Animals, Food and Humans in Denmark; Statens Serum Institut: Copenhagen, Denmark; National Veterinary Institute: Uppsala, Sweden; Technical University of Denmark National Food Institute: Kongens Lyngby, Denmark; Technical University of Denmark: Kongens Lyngby, Denmark, 2017.
  83. Hopman, N.E.; Portengen, L.; Heederik, D.J.; Wagenaar, J.A.; Van Geijlswijk, I.M.; Broens, E.M. Time trends, seasonal differences and determinants of systemic antimicrobial use in companion animal clinics (2012–2015). Vet. Microbiol. 2019, 235, 289–294. [Google Scholar] [CrossRef] [PubMed]
  84. Sutherland, M.E. Antibiotic use across the globe. Nat. Hum. Behav. 2018, 2, 373. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Saliva core bacteriome. The relative abundances of genera that achieved more than 1% of the bacterial hits in any of the samples. In the sample No. 20, no reads were classified to bacteria.
Figure 1. Saliva core bacteriome. The relative abundances of genera that achieved more than 1% of the bacterial hits in any of the samples. In the sample No. 20, no reads were classified to bacteria.
Antibiotics 11 01490 g001
Figure 2. Identifed ARGs by samples. For each sample–ARG combination, only the best finding is plotted. The size and the color of the dots correspond to the coverage and the sequence identity of hits on reference genes, respectively. In sample No. 20, there was no identifiable ARG. The gene names that are too long are abbreviated (acrA: E. coli acrA; emrE: E. coli emrE).
Figure 2. Identifed ARGs by samples. For each sample–ARG combination, only the best finding is plotted. The size and the color of the dots correspond to the coverage and the sequence identity of hits on reference genes, respectively. In sample No. 20, there was no identifiable ARG. The gene names that are too long are abbreviated (acrA: E. coli acrA; emrE: E. coli emrE).
Antibiotics 11 01490 g002
Figure 3. Mobile ARG frequency by bacteria of origin. The size of the dots indicates the occurrence frequency of the given gene flanked by iMGE, positioned in plasmid or prophage.
Figure 3. Mobile ARG frequency by bacteria of origin. The size of the dots indicates the occurrence frequency of the given gene flanked by iMGE, positioned in plasmid or prophage.
Antibiotics 11 01490 g003
Table 1. The list of analyzed samples was obtained from the National Center for Biotechnology Information Sequence Read Archive. Column Run contains the NCBI SRA run identifiers. Bacterial read count represents the number of reads that were classified taxonomically to any bacteria.
Table 1. The list of analyzed samples was obtained from the National Center for Biotechnology Information Sequence Read Archive. Column Run contains the NCBI SRA run identifiers. Bacterial read count represents the number of reads that were classified taxonomically to any bacteria.
IDBioProjectRunBacterial Read Count
1PRJNA648123SRR123300292,900,387
2 SRR1233004116,153,172
3 SRR1233004213,072,781
4 SRR1233004313,774,332
5 SRR123300446,123,646
6 SRR1233004516,707,766
7 SRR1233009818,826,266
8 SRR1233010427,598,592
9 SRR123302209,938,948
10 SRR1233026017,642,933
11 SRR1233029817,277,697
12 SRR1233035613,988,719
13 SRR1233036417,378,513
14 SRR1233037712,155,726
15 SRR1233037834,183,357
16 SRR1233038222,353,314
17 SRR1233038322,886,951
18 SRR1233038418,328,656
19 SRR123303856,631,504
20PRJNA683923SRR133405340
21 SRR133405356,752,169
22 SRR133405378,245,374
23 SRR1334053841,212,470
24 SRR1334053913,028,655
25 SRR133405406,964,460
26 SRR133405416,279,921
Table 2. Identified ARGs and the drug classes affected by them. The frequency columns show how many samples the genes occurred in.
Table 2. Identified ARGs and the drug classes affected by them. The frequency columns show how many samples the genes occurred in.
ARG(s)FrequencyDrug Class
n%
aac(6’)-Im27.7aminoglycoside
aad(6)27.7aminoglycoside
aadA213.8aminoglycoside
aadA313.8aminoglycoside
aadA513.8aminoglycoside
aadA1513.8aminoglycoside
aadS1246.2aminoglycoside
acrA13.8cephalosporin, fluoroquinolone, glycylcycline, penam, phenicol, rifamycin, tetracycline, triclosan
ant(2”)-Ia13.8aminoglycoside
ant(3”)-IIa13.8aminoglycoside
ant(6)-Ib13.8aminoglycoside
aph(2”)-IIa311.5aminoglycoside
aph(3”)-Ib1038.5aminoglycoside
aph(3’)-Ia519.2aminoglycoside
aph(3’)-IIa13.8aminoglycoside
aph(3’)-IIIa415.4aminoglycoside
aph(6)-Id1038.5aminoglycoside
bacA13.8peptide
blaACT-1213.8carbapenem, cephalosporin, cephamycin, penam
blaOXA-21246.2carbapenem, cephalosporin, penam
blaOXA-8513.8carbapenem, cephalosporin, penam
blaOXA-11913.8carbapenem, cephalosporin, penam
blaOXA-3471661.5carbapenem, cephalosporin, penam
blaROB-12180.8cephalosporin, penam
blaROB-913.8cephalosporin, penam
blaROB-10311.5cephalosporin, penam
blaTEM-11627.7cephalosporin, monobactam, penam, penem
catIII27.7phenicol
cfxA22076.9cephamycin
cmlA913.8phenicol
dfrA14311.5diaminopyrimidine
emrE13.8macrolide
emrK13.8tetracycline
ereA13.8macrolide
ermB934.6lincosamide, macrolide, streptogramin
ermF1869.2lincosamide, macrolide, streptogramin
ermG27.7lincosamide, macrolide, streptogramin
ermX13.8lincosamide, macrolide, streptogramin
fosA213.8fosfomycin
gadW13.8fluoroquinolone, macrolide, penam
gadX13.8fluoroquinolone, macrolide, penam
lnuB27.7lincosamide
lnuC27.7lincosamide
lsaE27.7lincosamide, macrolide, oxazolidinone, phenicol, pleuromutilin, streptogramin, tetracycline
mdtN13.8acridine dye, disinfecting agents and intercalating dyes, nucleoside
mef(En2)1246.2macrolide
mel519.2lincosamide, macrolide, oxazolidinone, phenicol, pleuromutilin, streptogramin, tetracycline
oqxA13.8diaminopyrimidine, fluoroquinolone, glycylcycline, nitrofuran, tetracycline
pgpB2388.5peptide
qacL27.7disinfecting agents and intercalating dyes
SAT-413.8nucleoside
sul1726.9sulfonamide
sul2830.8sulfonamide
tet321973.1tetracycline
tet4413.8tetracycline
tetH311.5tetracycline
tetM519.2tetracycline
tetO1869.2tetracycline
tetQ2076.9tetracycline
tetS13.8tetracycline
tetW519.2tetracycline
tetWNW519.2tetracycline
tetX1038.5glycylcycline, tetracycline
tetX113.8tetracycline
tetX413.8glycylcycline, tetracycline
tetX51038.5tetracycline
tetY13.8tetracycline
tetZ13.8tetracycline
ugd13.8peptide
Table 3. Identified ARGs and the predicted bacterial species of origin. For ten genes (aadA5, aadA15, ant(2”)-Ia, bacA, blaACT-12, cmlA9, fosA2, oqxA, tetX1 and tetY), no species-level prediction was obtained as to which bacterium the contig carrying the gene might have originated from.
Table 3. Identified ARGs and the predicted bacterial species of origin. For ten genes (aadA5, aadA15, ant(2”)-Ia, bacA, blaACT-12, cmlA9, fosA2, oqxA, tetX1 and tetY), no species-level prediction was obtained as to which bacterium the contig carrying the gene might have originated from.
ARG(s)Bacteria of Origin
aac(6’)-ImClostridioides difficile
aad(6)Enterococcus faecium, Staphylococcus aureus
aadA2Acinetobacter baumannii
aadA3Neisseria animaloris
aadSBacteroides fragilis, Capnocytophaga sp. H2931, Capnocytophaga sp. H4358, Capnocytophaga stomatis, Chryseobacterium indologenes, Riemerella anatipestifer
acrAE. coli
ant(3”)-IIaAeromonas hydrophila
ant(6)-IbAmedibacterium intestinale
aph(2”)-IIaC. difficile
aph(3”)-IbCorynebacterium sp. 1959, Haemophilus parahaemolyticus, Klebsiella michiganensis, Moraxella bovis, Variovorax sp. SRS16
aph(3’)-IaCorynebacterium sp. 1959, E. coli, Klebsiella quasipneumoniae, Variovorax sp. PAMC28562
aph(3’)-IIaE. coli
aph(3’)-IIIaE. faecium, S. aureus, Streptococcus agalactiae
aph(6)-IdCorynebacterium sp. 1959, K. michiganensis, Neisseria shayeganii, Providencia rettgeri, Variovorax sp. SRS16
blaOXA-2A. baumannii, A. hydrophila, P. aeruginosa
blaOXA-85Fusobacterium ulcerans
blaOXA-119Geobacter sulfurreducens
blaOXA-347Alistipes shahii, B. fragilis, Bacteroides heparinolyticus, Capnocytophaga sp. H2931, Capnocytophaga sp. H4358, C. stomatis, Chryseobacterium sp. POL2, Elizabethkingia anophelis, Empedobacter brevis, Myroides odoratimimus, R. anatipestifer
blaROB-1Actinobacillus pleuropneumoniae, Conchiformibius steedae, Glaesserella parasuis, Haemophilus haemolyticus
blaROB-9G. parasuis
blaROB-10Bibersteinia trehalosi
blaTEM-116E. coli
catIIIK. michiganensis
cfxA2Capnocytophaga cynodegmi, Parabacteroides distasonis, Porphyromonas cangingivalis, Porphyromonas crevioricanis, Porphyromonas gingivalis, Tannerella forsythia
dfrA14K. michiganensis, Ochrobactrum anthropi
emrEE. coli
emrKE. coli
ereAGeobacter daltonii
ermBEnterococcus gilvus, Enterococcus sp. FDAARGOS_375, S. agalactiae, Streptococcus suis
ermFA. shahii, B. fragilis, C. stomatis, C. indologenes, P. distasonis, P. cangingivalis, Prevotella intermedia, R. anatipestifer
ermGC. difficile
ermXTrueperella pyogenes
gadWE. coli
gadXE. coli
lnuBS. suis
lnuCStreptococcus equi, Streptococcus gwangjuense
lsaES. suis
mdtNE. coli
mef(En2)B. fragilis, P. cangingivalis, P. gingivalis, P. intermedia
melStreptococcus pluranimalium
pgpBP. gingivalis
qacLP. aeruginosa
SAT-4S. aureus
sul1A. hydrophila, P. aeruginosa
sul2Corynebacterium sp. 1959, H. parahaemolyticus, K. michiganensis, M. bovis, Pasteurella multocida, P. rettgeri
tet32Blautia hansenii, Bulleidia sp. zg-1006, C. difficile, Clostridium cellulovorans, Eubacterium maltosivorans, Eubacterium sp. NSJ-61, Faecalibacterium prausnitzii, Lachnoanaerobaculum umeaense, Peptoclostridium acidaminophilum, Roseburia intestinalis, Streptococcus anginosus, Streptococcus constellatus, S. equi
tet44A. intestinale
tetHProteus vulgaris, Pseudomonas putida
tetMC. difficile, Enterococcus faecalis, Mogibacterium pumilum, Streptococcus sp. FDAARGOS_521
tetOC. difficile, Enterococcus hirae, Murdochiella vaginalis, Streptococcus acidominimus, S. anginosus, S. constellatus, S. equi, S. suis
tetQAlistipes indistinctus, Bacteroides dorei, B. heparinolyticus, Bacteroides ovatus, Bacteroides sp. HF-5287, Phocaeicola coprophilus, P. crevioricanis, Prevotella fusca, P. intermedia
tetSStreptococcus parauberis
tetWEnterocloster bolteae, F. prausnitzii, Megasphaera stantonii, S. suis
tetWNWFilifactor alocis, M. pumilum
tetXB. fragilis, P. distasonis, P. intermedia, R. anatipestifer
tetX4R. anatipestifer
tetX5B. fragilis, C. stomatis, R. anatipestifer
tetZRothia nasimurium
ugdE. coli
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tóth, A.G.; Tóth, I.; Rózsa, B.; Dubecz, A.; Patai, Á.V.; Németh, T.; Kaplan, S.; Kovács, E.G.; Makrai, L.; Solymosi, N. Canine Saliva as a Possible Source of Antimicrobial Resistance Genes. Antibiotics 2022, 11, 1490. https://doi.org/10.3390/antibiotics11111490

AMA Style

Tóth AG, Tóth I, Rózsa B, Dubecz A, Patai ÁV, Németh T, Kaplan S, Kovács EG, Makrai L, Solymosi N. Canine Saliva as a Possible Source of Antimicrobial Resistance Genes. Antibiotics. 2022; 11(11):1490. https://doi.org/10.3390/antibiotics11111490

Chicago/Turabian Style

Tóth, Adrienn Gréta, Imre Tóth, Bernadett Rózsa, Attila Dubecz, Árpád V. Patai, Tibor Németh, Selçuk Kaplan, Eszter Gabriella Kovács, László Makrai, and Norbert Solymosi. 2022. "Canine Saliva as a Possible Source of Antimicrobial Resistance Genes" Antibiotics 11, no. 11: 1490. https://doi.org/10.3390/antibiotics11111490

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