Skip to main content

Pet cats may shape the antibiotic resistome of their owner’s gut and living environment

Abstract

Background

Companion animals can contribute to the physical and mental health of people and often live in very close association with their owners. However, the antibiotic resistome carried by companion animals and the impact they have on their owners and living environment remain unclear. In this study, we compared the ARG profiles of cats, humans, and their living environments using metagenomic analysis to identify the core ARGs in the cat and human gut and explore the potential impact of cats on ARGs in the human gut through the environment.

Results

Results showed that the abundance of ARGs in the cat gut was significantly higher than that in the human gut (P < 0.0001), with aminoglycoside and tetracycline resistance genes being the dominant ARGs in the cat gut. There was no significant difference in the abundance of total ARGs in the guts of cat owners and non-owners (P > 0.05). However, the abundance of aminoglycoside resistance genes including APH(2'')-IIa and AAC(6')-Im was significantly higher in cat owners than that in non-cat owners (P < 0.001). Also, ARG abundance was positively correlated with the frequency of cat activity in the living environment. Enterobacteriaceae was the dominant ARG host co-occurring in the cat gut, human gut, and living environment.

Conclusions

Our results show that cats may shape the living environment resistome and thus the composition of some ARGs in the human gut, highlighting the importance of companion animal environment health.

Video Abstract

Background

Companion animals can contribute to the physical and mental health of people and often live in very close association with their owners, sometimes even sharing the same bed. However, the potential biosecurity safety risks associated with close contact between companion animals and their owners are often overlooked. Dogs and cats, by far the most popular companion animals, are important sources of zoonotic infections [1]. Studies have shown that pet dogs and cats can carry multiple human-associated pathogens and a variety of multi-drug-resistant bacteria, including methicillin-resistant Staphylococcus aureus and beta-lactam antibiotic-resistant Enterobacteriaceae [2, 3]. These antibiotic-resistant microorganisms may be transmitted from companion animals to humans through direct contact, ectoparasites, and aerosols. In particular, vulnerable populations, including immunocompromised elderly and infants, may be at greater risk. Moreover, as the number of companion animals continues to increase, the use of antibiotics in the clinical consultation of pets has become more common. Antibiotics that are not used appropriately could promote antibiotic resistance [4], emphasizing the fact that the risk of antibiotic resistance during companion animal husbandry cannot be ignored.

The effects of companion animals on the resistome of their owners and living environment are still unclear. Previous studies found that the likelihood of airborne transmission of antibiotic resistance genes (ARGs) was extremely low [5]. However, communication between cats and owners is more likely to be physical than airborne, thus the transmission route of ARGs between cats and owners is more likely to be via the living environment, such as working surfaces or the floor. In addition, there is currently a lot of research focused on ARG abundance changes and a lack of research on mobile-associated ARGs and pathogen-associated ARGs, even though these ARGs pose a greater threat to humans [6, 7].

In this study, we compared the ARG profile of cats, humans, and their living environments using metagenomic analysis to identify the core ARGs in the cat and human gut and explore the potential impact of cats on ARGs in the human gut through the environment.

Materials and methods

Experimental design and sample collection

To investigate the differences in ARG profiles from the intestines of pet cats and humans, 60 fecal samples were collected from 30 pet cats and 30 volunteers from Guangzhou, China, in August 2020 (Fig. 1A). Additional information about the cats and volunteers can be found in Supplementary Table S1 and Table S2. All fecal samples were collected with a sterile fecal collector, delivered to the laboratory within 4 h, and stored in a − 80 ℃ refrigerator for further use.

Fig. 1
figure 1

Comparison of resistome between cat gut and human gut. A Sampling flowchart; the B β-diversity, C abundance, D number, and E α-diversity of ARGs in cat and human gut. The core ARGs in F cat gut and G human gut. H, I ARG abundance in the gut of the cat owner and non-cat owners

Next, the living environment of one of the cat owners was chosen (Sample ID: Human_4). And we collected samples from the cat-keeping environment, including areas such as the living room, dining room, balcony, small bedroom, medium bedroom, large bedroom, and cat litter (located in the toilet). The floor of each sampling point was wiped with a sterile swab in both horizontal and vertical directions. The homeowner cleaned the floor every 3 days, so we collected floor samples before each cleaning, a total of 5 times during April 2022. Ten swabs were collected at each sampling site each time and mixed into one sample, resulting in 35 mixed samples from the living environment. Finally, we downloaded 19 metagenomic sequencing samples from hospital environments from the National Center of Biotechnology Information (NCBI) website to compare the ARGs in the cat-keeping environment and the clinical environment.

Metagenomic sequencing and microbial community analysis

Total DNA was extracted from samples using the QIAamp PowerFecal DNA Kit (Qiagen, Germany) following the manufacturer’s instructions. DNA samples were subjected to metagenomic sequencing at Novogene on the Illumina NovaSeq platform and 6 ~ 8 G of raw data were obtained for each sample. The metagenome raw data from all samples, including those from the hospital environment, were quality-controlled using KneadData (Trimmomatic 0.33) with default parameters (SLIDINGWINDOW:4:20 MINLEN:50) [8]. Subsequently, all human and cat DNA were identified with reference to the human (Homo_sapiens_Bowtie2_v0.1) and cat (GCF_018350175.1) genomes using bowtie2 software (–very-sensitive), and removed [9, 10]. The cleaned data of each sample was used for microbial community analysis based on the standard database of the Kraken2 software [11]. Finally, the clean data were assembled using Megahit v1.2.9 (based on the default kmer) [12]. These contigs were then subjected to open reading frames (ORFs) prediction using Prodigal software (v2.6.3). Thereafter the ORFs were clustered using Cd-hit (v4.8.1) with default parameters (-c 0.95). The ORFs were quantified using Salmon software (v1.4.0) [13,14,15].

Antibiotic resistome analysis

The Comprehensive Antibiotic Resistance Database (CARD v3.2.5, https://card.mcmaster.ca/) was downloaded, collated, and used to build the ARG database using DIAMOND software (v2.0.15.153) [16]. Similarly, mobileOG-db (1.6 v1) was downloaded and used to build the mobile genetic element (MGE) database [17]. The ORFs files were compared with the ARG and MGE databases using DIAMOND software with the following parameters: e value 1e − 7, max-target-seqs 70, id 80. ARG-like-ORFs and MGE-like-ORFs were identified. The number of prokaryotic cells per sample was estimated using ARGs_OAP v3.2, and the abundance of ARGs and MGEs was corrected to copies/cell number [18].

Identification of ARG hosts with MAG binning

The clean metagenome data of all samples were finely assembled using Metaspades (v3.14.0) with default parameters and the contigs were used to assemble metagenome-assemble genomes (MAGs) [19]. The MAGs were generated by Metabat2, Maxbins, and CONCOCT in the MetaWRAP pipeline [20]. Finally, MAGs were output based on > 50% completeness and < 10% contamination. Species classification of MAGs was performed using GTDB-TK v2.1.1 with the GTDB database (R207_v2) [21, 22]. To obtain ARG and MGE hosts, ARGs and MGEs in the MAGs were identified as in “Metagenomic sequencing and microbial community analysis” section. Virulence factors (VFs) in the MAGs were identified by comparing data to the VFDB_setA database (Last update: October 7, 2022), and these MAGs were considered bacterial pathogens in the gut of cats, humans as well as in the living environment.

Identification of ARG public risks

The contigs from the metagenome assembly were compared with ARG, MGE, and VF databases using the BLAST tool to obtain the abundance of ARGs (ARGs-like contigs), mobile-associated ARGs (ARGs-MGEs like contigs) and mobile-associated ARGs in pathogenic bacteria (ARGs-MGEs-VFs like contigs) [7]. MetaCompare pipeline was then used to assess the resistome risk in the gut of cats and humans as well as in the hospital and living environment [6]. And MetaCHIP pipeline was used to assess the horizontal transfer risk of these ARGs [23].

Data analysis and presentation

All data were first collated in WPS office 11.1.0 and significance analysis was performed using SPSS 22.0 (one-way ANOVA). And Tukey’s multiple comparisons test was used. Column graphs were plotted on GraphPad Prism 8.01, while box plots for gene and microbial analysis were drawn on STAMP 2.1.3. Principal coordinates analysis (PCoA) and redundancy analysis (RDA) were performed in R 4.2.2. Phylogenetic trees were drawn using iTOL (https://itol.embl.de/) [24]. Other images were generated in OriginPro 2023. Adobe Illustrator 22.1 was used for the graphic layout.

Results

Comparison of resistome in the gut of cats and humans

Results showed significant differences between the ARG β-diversity in the cat and human gut (P = 0.001) (Fig. 1B). The abundance of ARGs in the cat gut was 1.809 ± 0.070 copies/cell number, which was significantly higher than that in the human gut (1.765 ± 0.185 copies/cell number) (P < 0.0001) (Fig. 1C). The proportion of ARGs in the cat gut was 0.07% of the total number of genes, which was higher than that in the human gut (0.03%) (Fig. 1D). The Richness and Shannon index of ARGs in the cat gut was also higher than that in the human gut (Fig. 1E). The core ARGs in the cat gut were tetracycline, aminoglycoside, and MLS resistance genes, and the core ARGs in the human gut were tetracycline, MLS, and multidrug resistance genes. The abundance of tetracyclines (e.g., tetO and teW) and aminoglycoside resistance genes (e.g., APH(2'')-IIa and ACC(6')-Im) in the cat gut was significantly higher than that in the human gut (P < 0.05) (Fig. 1F, G and Figure S1). In addition, we found that although there was no significant difference between the total ARG abundance in the gut of the cat owner and non-cat owners (Fig. 1H), the abundance of ARGs such as APH(2'')-IIa and ACC(6')-Im were higher in the cat owner than those in non-cat owners (Fig. 1I). The abundance of these specific ARGs in the cat gut was also significantly higher than that in the human gut (P < 0.001).

Risk of resistome transmission from cats to owner and living environment

Results from the MGEs composition analysis in the cat gut showed that the composition of MGEs in the cat gut and human gut differed significantly (P = 0.001) (Fig. 2A). The abundance of MGEs in the human gut was significantly higher than that in the cat gut (P < 0.0001), with 76.61 ± 3.28 copies/cell number and 47.03 ± 1.74 copies/cell number, respectively (Fig. 2B). The mobile-associated ARGs (contigs of ARGs-MGEs) in the human gut accounted for 0.59 ± 0.060% of the total contigs, which was significantly higher than the abundance in the cat gut (0.576 ± 0.068%) (P < 0.05) (Figure S2). The main types of VFs in the cat gut were adherence, immune modulation, and urease, while the main types of VFs in the human gut were adherence, nutritional/metabolic factors, and Immune modulation. The total abundance of VFs in the human gut was 1.319 ± 0.494 copies/cell number, which was significantly higher than that in the cat gut (0.183 ± 0.021 copies/cell number) (P < 0.05) (Fig. 2C, D). However, there was no significant difference in the proportion of pathogen-mobile-associated ARGs (contigs of ARGs-MGEs-VFs) in the cat gut and the human gut. We also found a non-significant difference in resistome risk scores between the human and cat gut (Fig. 2E, F).

Fig. 2
figure 2

Resistome risk in cat, human, and their living environment. The A β-diversity and B abundance of MGEs in cat and human gut. C, D Abundance of VFs. E Proportion of pathogen-mobile-associated ARGs. F Resistome risk score in cat and human gut. G Proportion of pathogen-mobile-associated ARGs in the gut of the cat owner and non-cat owners. H Resistome risk score in cat and human gut in the gut of the cat owner and non-cat owners. I Sampling flowchart of the living environment. The more cats in the figure, the higher the frequency of cat activity. J ARG abundance and K resistome risk score in the living environment

In addition, there was no significant difference in the proportion of pathogen-mobile-associated ARGs (contigs of ARGs-MGEs-VFs) in the guts of the cat owner and the non-cat owner, but the resistome risk score in the cat owner’s gut was significantly higher than that of the non-cat owner’s gut (P < 0.05) (Fig. 2G, H). These suggest that cat-derived ARGs have a risk of affecting the resistome in the human gut. We collected samples from the cat-keeping environment for the study (Fig. 2I). In this study, the large room was the most active area for cat, followed by the dining room, small room, medium room living room, and the balcony. For the living environment analysis, we found the highest abundance of ARGs in the cat litter and the large bedroom with 2.414 ± 0.287 copies/cell number and 2.411 ± 0.330 copies/cell number, respectively, followed by the dining room, small bedroom, medium bedroom and living room (Fig. 2I, J). The lowest abundance of ARGs was found on the balcony with 0.744 ± 0.104 copies/cell number. The abundance of ARGs in the living environment was positively correlated with the frequency of cat activity. This suggests a possible transmission of ARGs from the cat gut to the living environment, which in turn affects the owner. Rather surprisingly, the risk of antibiotic resistome transfer was higher in the living environment than in the hospital (Fig. 2K). reminding us of the importance of focusing on the resistome risk in the keeping-cat living environment.

Host composition characteristics and horizontal transfer risk of ARG

There were significant differences in the bacterial community composition between the cat and human gut, with the highest relative abundance of Actinobacteria (0.45 ± 0.24) and Firmicutes (0.52 ± 0.23) in the cat gut and Firmicutes (0.49 ± 0.20) and Bacteroidetes (0.40 ± 0.23) in the human gut (Fig. 3A). In addition, the abundance of Bifidobacterium in the cat gut was significantly higher than that in the human gut, and the abundance of Bifidobacterium was higher in the gut of the cat owner compared to non-cat owners (Fig. 3B and Supplementary Figure S3). This suggested that cat microbes may influence the microbial composition of the human gut, perhaps including ARG hosts.

Fig. 3
figure 3

Bacterial community compositions in the cat and human gut. Relative abundance of A top 4 phyla and B Bifidobacterium; co-occurrence pattern of ARGs and bacterial communities in C cat gut and D human gut

We mined for ARG hosts using co-occurrence analysis and found that cats and humans share a common host in Enterobacteriaceae, which was significantly and positively correlated with 12 and 46 ARGs (ρ > 0.6, P < 0.01), respectively (Fig. 3C, D). To further identify the major ARG host in the cat and human gut, we then used a metagenome assembly to obtain 828 MAGs, 161 from the cat gut, 313 from the human gut, and 351 from the living environment, composed mainly of Firmicutes, Actinobacteriota, Bacteroidota, and Proteobacteria (Figure S4). Most of the bacteria at different phylum levels identified from cat and human feces were found to carry ARGs, suggesting that ARGs are prevalent in environmental bacteria (Fig. 4 and Figure S5). The Enterobacteriaceae were found to be the major ARG hosts in the cat and human gut, being detected to carry 65 (bin. 62) and 73 (bin. 157) ARGs. More worryingly, the Enterobacteriaceae includes common pathogenic bacteria such as Escherichia, Salmonella, and Shigella, which increase the risk of antibiotic resistance. Further horizontal transfer assessment of these MAGs revealed that bin62 (Enterobacteriaceae), the main ARG host in the cat gut, is at risk of gene transfer to living environmental bin317 (Fig. 5). There was also a risk of horizontal transfer between different MAGs in the living environment, such as between bin197 and bin259. Finally, we identified a risk of horizontal transfer of bacterial genes from the living environment to the human gut, such as between bin197 and bin92 in the environment, and bin138 and bin157 in the human gut. Among them, bin157 was the main ARG host in the human gut.

Fig. 4
figure 4

The number of ARGs carried by MAGs in cat gut

Fig. 5
figure 5

Horizontal gene transfer occurring in MAGs from cat gut, human gut, and environment

Discussion

Companion animals shape the resistome in the human gut

The intestinal tract contains a large number of microorganisms that make up a relatively constant resident bacteria and transient bacteria that come in from the external environment [25]. Resident bacteria are the main components of the intestinal flora and maintain a dynamic balance with the environment and the host by promoting normal physiological functions, such as nutrition, immunity, and digestion in the host [26]. Transient bacteria are harmless or pathogenic foreign bacteria that enter the intestinal tract through contact and aerosols. Normally, resident bacteria adheres, colonizes, and multiplies in the intestine, forming a barrier that inhibits and rejects the transient bacteria through antagonism, thereby protecting intestinal homeostasis and host health. However, through prolonged exposure, microorganisms in the living or working environment also shape the intestinal flora, and changes in the composition of microorganisms, which are important hosts for ARGs, can lead to changes in the resistome [27].

In this study, we found that the composition and diversity of the bacterial community in the cat gut differed significantly from those in the human gut. Also, the α- and β-diversity of ARGs in the cat gut differed significantly from those in the human gut, and the abundance of ARGs in the cat gut was significantly higher than that in the human gut, including all types of ARGs. Bacteria carrying a high risk of mcr-1 resistance gene have also been isolated from cats and dogs [28]. Therefore, the higher risk of antibiotic resistance present in the cat gut cannot be ignored. In addition, this study found that the abundance of APH(2''-IIa and ACC(6')-Im resistance genes in the gut was significantly higher in the cat owner than in non-cat owners. These resistance genes are also core ARGs in the cats tested here. Previous studies have also shown a strong correlation in macrolide resistance genes between pet dogs and their owners [29]. Belas et al. (2020) found that the resistome of companion animals is similar to those of people in close contact with them [30]. Thus, it can be assumed that companion animals shape the gut resistome of their owners to a certain extent.

We also showed that the core ARG hosts in both the cat gut and human gut were Enterobacteriaceae and a risk of horizontal gene transfer between these core ARG hosts and different bacteria in the living environment and human gut. It can be hypothesized that these core ARG hosts in the feline gut can enter the human gut through different pathways as transient bacteria and further alter the resistome of the owner’s gut through horizontal gene transfer. Of course, it is also possible that these ARGs and their hosts may attach to and colonize the human intestine over longer periods. Studies have shown that the gut microbiome and resistome of feeders working on farms for long periods can be reshaped [31, 32]. More worryingly, the core ARG hosts identified in this study, Enterobacteriaceae, also include pathogenic bacteria such as Escherichia, Salmonella, and Shigella, which may pose a more complex risk of antibiotic resistance and pathogenicity.

It is very unfortunate that only 30 volunteers and 30 cats were sampled in the first part of the trial in this study. Our original plan was to collect more than 100 volunteer samples, but COVID-19 pandemic outbreak suddenly became severe and some volunteers had to withdraw from our study. Although the study did not go as smoothly as expected, the results obtained have important reference value.

Pathways by which companion animals influence the human resistome

The main route of transmission of ARGs and their hosts in the environment is through contact, such as air and surfaces [33]. It has been suggested that air is an important pathway for companion animals to influence the human resistome. Because of the distinctive smell in cat-feeding households, it was speculated that antibiotic-resistant microorganisms and ARGs volatilized from cats might be floating in the air of cat-feeding environments. Besides, previous studies have also confirmed that microorganisms and ARGs in the environment can enter humans in the form of aerosols [34, 35]. Surprisingly, our results showed that the microorganisms and ARGs volatilized from the cat gut (feces) into the air were minimal, suggesting that air is not the main route for the transmission of ARGs from cats to owners [5]. Pet cats are curious animals and will cover almost the entire living space in their range of activities. Therefore, we collected floor samples from various areas of the living environment and found that the abundance of ARGs in these areas was positively correlated with the frequency of pet cat activity. Horizontal gene transfer was found to occur between the cat gut, the environment, and the human gut, making it feasible that pet cats can affect the human resistome through their activities in the living environment, and that this effect is greater than the airborne effect. This is especially true for rooms where cats are more active.

Due to homeostatic regulation of the intestine [36], microorganisms that enter the human gut from the cat and the living environment, including ARG hosts, are normally transient for the human gut and rarely attach to and colonize the gut. We found a higher abundance of MGEs in the human gut samples than in the cat and the environmental samples, including phage, integron, and transposon types. These MGEs act as important factors in the horizontal gene transfer of ARGs in different microorganisms [36]. This suggests that there is a greater risk of ARG horizontal gene transfer in the human gut. Also, horizontal gene transfer is an important influencing factor for environmental resistome variation, with MGEs explaining up to 10.3% of the ARG variation in wastewater; 13.9% of the ARG variation in pig intestines, and 13.76% of the ARG variation in dog intestines [29, 37, 38]. Additionally, we then analyzed the contigs of ARGs and MGEs and found that the proportion of the mobile ARGs (MGE-ARG-like contigs) in the cat gut and the human gut was 0.576 ± 0.068‰ and 0.590 ± 0.060‰ of the total contigs, respectively. Horizontal gene transfer was also found between different microorganisms in cat feces, the living environment, and the human gut.

In summary, we mapped the shaping process of the resistome from the cat to the human gut and living environment (Fig. 6). The ARGs and their hosts in cats are brought to the living environment through the activity of the cat, shaping the resistome in the environment, which in turn affects the human gut resistome. There are also a small number of ARG host that enter the human body through the respiratory tract via aerosols [5]. Of course, direct contact with cats may be the main route through which ARG hosts enter the human gut [39]. The ARGs that enter the human gut are then repositioned in the intestine, such as by horizontal gene transfer, to shape the resistome of the owner.

Fig. 6
figure 6

The formation of resistance from the cat to the human gut and living environment

The resistome risk from companion animals is low

ARGs are commonly found in natural environments, even in extreme environments, such as the deep-sea and Arctic permafrost [40, 41]. Therefore, the presence of ARGs or the presence of the horizontal spread of ARGs should not be considered a high resistome risk. There are several methods with which to assess the resistome risk of the environment. The most direct approach is to detect the abundance of ARG score in the environment, both in absolute and normal abundance. It is believed that the higher the abundance of ARGs in the environment, the higher the resistome risk. Resistome risk could also be estimated by combining the analysis of the abundance of total ARGs, mobile ARGs (MGEs-ARGs like contigs), and mobile ARGs associated with pathogenic bacteria (MGEs-VFs-ARGs like contigs) [6, 7]. In this study, we used both methods to assess resistome risk in cat-keeping households and found no significant differences in the abundance of ARGs in cat owners and non-cat owners, as well as no significant differences in the resistome risk scores in the human and cat gut. In addition, we found that the abundance of Bifidobacterium, which is considered to be a beneficial bacterium, was significantly higher in the cat gut than that in the human gut. Also, the abundance of Bifidobacterium in the gut of the cat owner was significantly higher than that of non-cat owners. This evidence suggests that cat ownership can shape the resistome in the owner's gut and living environment, but that the resistome risk is small. Additionally, we could not find any reports of antibiotic-resistant contamination due to pet ownership. Despite this, cleanliness is important when keeping cats. In addition, ARGs can be reduced by adjusting the cat’s diet and the rational use of antibiotics for pets to prevent ARG enrichment.

Conclusion

We found that the ARG abundance in the cat gut was significantly higher than that in the human gut, the resistome risk in the gut of the cat owner was higher than that in non-cat owners, and the ARG abundance in the living environment was positively correlated with the frequency of cat activity. We also showed that the cat gut and the human gut share the same core ARG hosts (Enterobacteriaceae). These results suggest that pet cat may shape the antibiotic resistome in human gut. And the horizontal gene transfer between bacteria in the environment as well as in the human gut can further contribute to resistome shaping. We are reminded of the resistome risk during companion animal keeping and to minimize the transmission pathways of companion animal-derived resistance by managing the living environment and ensuring animal health.

Availability of data and materials

All sequencing data for this study are available from the NCBI database (Accession number: PRJNA944553). The accession number for the hospital environmental metagenome sequencing data used in this study is PRJNA726763.

References

  1. Christou L. The global burden of bacterial and viral zoonotic infections. Clin Microbiol Infect. 2011;17(3):326–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Jung WK, et al. Distribution and antimicrobial resistance profiles of bacterial species in stray cats, hospital admitted cats, and veterinary staff in South Korea. BMC Vet Res. 2020;16(1):109.

  3. Faires MC, et al. Inducible clindamycin-resistance in methicillin-resistant Staphylococcus aureus and methicillin-resistant Staphylococcus pseudintermedius isolates from dogs and cats. Vet Microbiol. 2009;139(3–4):419–20.

    Article  PubMed  CAS  Google Scholar 

  4. Yang Y, et al. Changes in the carbon metabolism of escherichia coli during the evolution of doxycycline resistance. Front Microbiol. 2019;10:2506.

  5. Yang Y, et al. Abundance, diversity and diffusion of antibiotic resistance genes in cat feces and dog feces. Environ Pollut. 2022;292:118364.

    Article  PubMed  CAS  Google Scholar 

  6. Oh M, et al. MetaCompare: a computational pipeline for prioritizing environmental resistome risk. FEMS Microbiol Ecol. 2018;94(7):fiy079.

  7. Zhang A, et al. An omics-based framework for assessing the health risk of antimicrobial resistance genes. Nat Commun. 2021;12(1):4765.

  8. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Langmead B, et al. Scaling read aligners to hundreds of threads on general-purpose processors. Bioinformatics. 2019;35(3):421–32.

    Article  PubMed  CAS  Google Scholar 

  11. Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genom Biol. 2014;15(3):R46.

  12. Li D, et al. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6.

    Article  PubMed  CAS  Google Scholar 

  13. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

    Article  PubMed  CAS  Google Scholar 

  14. Patro R, et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Hyatt D, et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.

    Article  PubMed  CAS  Google Scholar 

  17. Brown CL, et al. mobileOG-db: a manually curated database of protein families mediating the life cycle of bacterial mobile genetic elements. Applied Environ Microbiol. 2022;88(18):e0099122.

  18. Yin X, et al. ARGs-OAP v2.0 with an expanded SARG database and Hidden Markov Models for enhancement characterization and quantification of antibiotic resistance genes in environmental metagenomes. Bioinformatics. 2018;34(13):2263–70.

    Article  PubMed  CAS  Google Scholar 

  19. Nurk S, et al. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome. 2018;6(1):158.

  21. Chaumeil P, et al. GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics. 2022;38(23):5315–16.

  22. Parks DH, et al. GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy. Nucleic Acids Res. 2022;50(D1):D785–94.

    Article  PubMed  CAS  Google Scholar 

  23. Song WZ, et al. MetaCHIP: community-level horizontal gene transfer identification through the combination of best-match and phylogenetic approaches. Microbiome. 2019;7(1):36.

  24. Letunic I, Bork P. Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011;39:W475–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Derrien M, Vlieg J. Fate, activity, and impact of ingested bacteria within the human gut microbiota. Trends Microbiol. 2015;23(6):354–66.

    Article  PubMed  CAS  Google Scholar 

  26. Donaldson GP, Lee SM, Mazmanian SK. Gut biogeography of the bacterial microbiota. Nat Rev Microbiol. 2016;14(1):20–32.

    Article  PubMed  CAS  Google Scholar 

  27. Caballero-Flores G, Pickard JM, Nunez G. Microbiota-mediated colonization resistance: mechanisms and regulation. Nat Rev Microbiol. 2022;21(6):347–60.

  28. Lei L, et al. mcr-1 in Enterobacteriaceae from companion animals, Beijing, China, 2012–2016. Emerg Infect Dis. 2017;23(4):710–1.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhao, R., et al., The co‐occurrence of antibiotic resistance genes between dogs and their owners in families. iMeta, 2022;1(2):e21.

  30. Belas A, et al. Sharing of clinically important antimicrobial resistance genes by companion animals and their human household members. Microb Drug Resist. 2020;26(10):1174–85.

    Article  PubMed  CAS  Google Scholar 

  31. Sun J, et al. Environmental remodeling of human gut microbiota and antibiotic resistome in livestock farms. Nat Commun. 2020;11(1):1427.

  32. Ding D, et al. Effect of cattle farm exposure on oropharyngeal and gut microbial communities and antibiotic resistance genes in workers. Sci Total Environ. 2022;806:150685.

    Article  PubMed  CAS  Google Scholar 

  33. Ginn O, Lowry S, Brown J. A systematic review of enteric pathogens and antibiotic resistance genes in outdoor urban aerosols. Environ Res. 2022;212:113097.

    Article  PubMed  CAS  Google Scholar 

  34. Bai H, et al. Spread of airborne antibiotic resistance from animal farms to the environment: dispersal pattern and exposure risk. Environ Int. 2022;158:106927.

    Article  PubMed  CAS  Google Scholar 

  35. Yang F, et al. Revealing the distribution characteristics of antibiotic resistance genes and bacterial communities in animal-aerosol-human in a chicken farm: From One-Health perspective. Ecotoxicol Environ Saf. 2021;224:112687.

    Article  PubMed  CAS  Google Scholar 

  36. Coyte KZ, Schluter J, Foster KR. The ecology of the microbiome: Networks, competition, and stability. Science. 2015;350(6261):663–6.

    Article  PubMed  CAS  Google Scholar 

  37. Yang Y, et al. Profiles of bacteria/phage-comediated ARGs in pig farm wastewater treatment plants in China: association with mobile genetic elements, bacterial communities and environmental factors. J Hazard Mater. 2021;404:124149.

    Article  PubMed  CAS  Google Scholar 

  38. Yang Y, et al. Short-term cold stress can reduce the abundance of antibiotic resistance genes in the cecum and feces in a pig model. J Hazard Mater. 2021;416:125868.

    Article  PubMed  CAS  Google Scholar 

  39. Wei X, et al. Charting the landscape of the environmental exposome. iMeta. 2022;1(4):e15.

  40. Kim H, et al. Characterization of antimicrobial resistance genes and virulence factor genes in an Arctic permafrost region revealed by metagenomics. Environ Pollut. 2022;294:118634.

    Article  PubMed  CAS  Google Scholar 

  41. Su H, et al. Diversity, abundances and distribution of antibiotic resistance genes and virulence factors in the South China Sea revealed by metagenomic sequencing. Sci Total Environ. 2022;814:152803.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The study was supported by the Project funded by China Postdoctoral Science Foundation (2022T150226) and Guangdong Provincial Science and Technology Plan Project (2020B1212060060).

Author information

Authors and Affiliations

Authors

Contributions

Yiwen Yang: Conceptualization, Investigation, Methodology, Formal analysis, Writing - original draft, Project administration, Funding acquisition. Xinwen Hu: Investigation, Methodology, Writing - review & editing. Shuang Cai: Writing - review & editing, Methodology. Nan Hu: Investigation. Yilin Yuan: Investigation.  Yinbao Wu: Writing - review & editing. Yan Wang: Writing - review & editing. Jiandui Mi: Formal analysis. Xindi Liao: Methodology, Formal analysis, Project administration.

Corresponding authors

Correspondence to Yiwen Yang or Xindi Liao.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was granted by the Ethics Committee of the Second Affiliated Hospital of Guangzhou Medical University (2019-KY-166C). Animal ethical approval was granted by the Ethics Committee of South China Agricultural University.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Information about the cats in this study. Table S2. Information about the volunteers in this study. Figure S1. The core ARGs in (A) cat gut and (B) human gut. Figure S2. Relative abundance of the mobile-associated ARGs. Figure S2. Relative abundance of the mobile-associated ARGs. Figure S4. Relative abundance of the MAGs from cat gut, human gut and living environment. Figure S5. ARGs abundance in MAGs from the human gut.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, Y., Hu, X., Cai, S. et al. Pet cats may shape the antibiotic resistome of their owner’s gut and living environment. Microbiome 11, 235 (2023). https://doi.org/10.1186/s40168-023-01679-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-023-01679-8

Keywords