Grazing disturbance increased the mobility, pathogenicity and host microbial species of antibiotic resistance genes, and multidrug resistance genes posed the highest risk in the habitats of wild animals

There are frequent exchanges of antibiotic-resistant bacteria and their antibiotic resistance genes (ARGs) between the external environment or livestock and wild animals. Grazing disturbance (GD) is a major factor that causes dramatic changes in China’s nature reserves. Studying the risk of ARGs in China’s nature reserves under GD has great significance for assessing the health of the habitats where wild animals live. In our study, the potential ecological risks of ARGs in soil microorganisms of wild animal habitats under GD were analyzed. Our results showed that the diversities of the ARGs in GD were higher than those in check control (CK) that not be disturbed by grazing, and 46 types of ARGs were only checked in GD. The ARGs were only carried by plasmids, and their relative abundances and the numbers of ARGs carried by plasmids were increased by GD, six types of transferred open reading frames (ORFs) carrying ARGs were detected only in GD. GD increased the relative abundances of drug-resistant pathogens, and the pathogen, Acinetobacter baumannii ACICU, that was only found in GD. Our results indicated that GD increased the risk of ARGs to the habitat of wild animals by reducing the total soil microbial species, increasing the numbers, species and mobility of ARGs in soil microorganisms and the species and abundance of drug-resistant pathogens in soil, as well as the ARG carrying capacity of microorganisms. Multidrug resistance genes (MRGs) posed the highest risk in the habitats of wild animals, and GD increased its risk because the largest proportion of ARGs carried by microorganisms were MRGs, which were also the most abundant ARGs carried by plasmids, and the highest proportion of ARGs carried by Proteobacteria (the dominant host bacterium of ARGs) were MRGs. Thus, GD may cause microorganisms in the habitats of wild animals become resistant to many antibiotics. MRGs in soil microorganisms could be used as an indicator for predicting the risk of GD in the habitats of wild animals.


Introduction
In China, a great deal of money and resources have been spent to build nature reserves containing wild animals, and they have achieved great success . However, most of the nature reserve areas are mountainous, their economic development level is limited, the cultural level of the indigenous residents is low, and animal husbandry is the economic source for the people to survive in these nature reserves . Human disturbances are the main factors affecting the survival of wild animals in nature reserves because the fundamental contradiction between humans and wild animals in using nature reserves cannot be reconciled. Results showed that the number of livestock in 2015 increased by 43.37 percent compared to 2003 , and there were 46% and 56% of the roaming areas of horses and cattle are distributed in the habitat of wild animals (Lia et al., 2017). Among the various human disturbances, grazing disturbances (GD) have the highest occurrence rate in the habitats of wild animals in Southwest China (State Forestry Administration, 2021). Previous studies have shown that the feces produced by grazing animals could be used as an index to evaluate the nutrition of soil in the wild animal habitat because grazing could affect the characteristics of the physical, chemical and microbial communities in the soil through excrement discharged into the soil, which affects the microecological environment of the soil in nature reserves (Orellana et al., 2019). However, due to changes in the human living environment, the extensive use of antibiotics in animal husbandry has led to an increase in the abundances of antibiotic resistance genes (ARGs) in microorganisms in livestock manure and has polluted the surrounding environment (Jadeja & Worrich, 2022).
ARGs can be accumulated in the food chain and then affects the composition, structure and function of the intestinal microorganisms in the animals. Because ARGs can be carried by mobile genetic elements (MGEs) of microorganisms and will transfer along with the food chain by plant microbiome, the plant microbiome represents a major pathway by which animals are exposed to microbes and genes consumed with food, and affect the treatment of animal diseases (Cerqueira et al., 2020;Mei et al., 2021). To further demonstrate the transfer of ARGs in the food chain, the Folsomia candida and Hypoaspis aculeifer were used as a model to study the transfer of ARGs in food chain, which showed that ARGs were transferred from the microbial community of soil to Hypoaspis aculeifer. Eventually, the ARGs accumulate in the intestinal microorganisms of Hypoaspis aculeifer and their characteristics of intestinal microbes were changed . Livestock, wild animals, humans, and their food have extensive and rapid exchanges of ARGs, and they are interconnected sources of ARGs (Brealey et al., 2021). The nature reserves of wild animals are special ecosystems that are relatively less disturbed by humans than human ecosystems, and the ARGs in the intestinal microorganisms of wild animals could be used as an ecological indicator to indicate the ARG pollution in the habitats of wild animals (Brealey et al., 2021). That is, when the habitats of wild animals are polluted by ARGs from grazing, the risk posed by ARGs to wild animals may increase. Our previous study also showed that the total abundance of ARGs in the intestinal microorganisms of wild animals in GD was significantly higher than that in the control check (CK), which was not disturbed by grazing (Xia et al., 2022). Therefore, it is necessary to detect the characteristics of ARGs in the habitats of wild animals under GD to evaluate the potential ecological risks of ARGs to wild animals. Previous results have already directly addressed the ARGs in intestinal microorganisms of wild animals to reflect the risk of contamination from antibiotics and their ARGs to wild animals under human interference (Huang et al., 2022). Unfortunately, there are few studies that use ARGs from soil microorganisms in wild animal habitats to assess the risks of ARGs to their survival under GD.
The accessibility, mobility, pathogenicity and clinical availability are the factors affecting ARGs' ecological risk, which indicated that whether the ARGs can be accessed, the ARGs can be moved, or whether the host bacteria with ARGs is pathogenic bacteria and the danger of pathogenic bacteria must be considered Zhang et al., 2022). Thus, it is difficult to assess the risks of ARGs in the habitats of wild animals under GD. Large amounts of microbial information can be reflected by metagenomic analysis (Guo et al., 2017). Recently, researchers have used metagenome sequencing technology not only to obtain the species and abundances of ARGs and calculate the ability of ARGs to be transferred by MGEs but also to analyze pathogenic genes in soil microorganisms (Ju et al., 2019). If a VF and a ARG occurrence in one open reading frame (ORF) of the contigs, the host bacteria of the ORF is defined as potential drug-resistant pathogens, and it represent the highest risk of ARGs . The species, abundances and the ability of ARGs together with pathogenic genes could be used to evaluate the ecological risk of ARGs (Liu et al., 2021). Thus, the ecological risks of ARGs could be comprehensively evaluated by metagenome sequencing technology. Additively, microorganisms are the hosts of ARGs, and it is necessary to consider the ability of microorganisms to transfer and store ARGs when their risks in ecosystems containing wild animals are evaluated under GD.
In conclusion, a hypothese be put forward that ARGs in soil microorganisms under GD is one of the potential factors that may threat the survival of wild animals. In this study, we investigated the ARGs of soil in GD and CK using 16S rRNA sequencing technology and metagenomic sequencing, with the aim of 1) analyze the characteristics and horizontal transfer potential of ARGs; 2) study the pathogenicity of their host bacteria in the soil microorganisms of the wild animal habitats under GD and CK; Frontiers in Environmental Science frontiersin.org 3) research the ability of soil microorganisms to transfer and store ARGs. Our results could reflect the risk of ARGs in the soil microorganisms of the habitats under GD, which could provide a theoretical basis for supervising this risk in the habitats of wild animals and to formulate emergency measures for health monitoring of wild animals.
2 Materials and methods

Sample collection
Baihe National Nature Reserve is a habitat for endangered wild animals such as giant pandas (Ailuropoda melanoleuca) and snubnosed golden monkeys (Rhinopithecus roxellana) in China, but it is not included in China's Giant Panda National Park, so there are some problems such as lax management and control at present (Liu et al., 2021). The original residents of the nature reserve depended on animal husbandry for their livelihoods, and 19.69% of the total income of the original residents came from animal husbandry, but the feeding and management methods of animal husbandry are extensive (Yuan, 2018). In the nature reserve, the grazing traces even overlap with the home ranges of the Sichuan snub-nosed monkeys. Therefore, our study was conducted in the Baihe National Nature Reserve.
In this study, soil samples were taken from the habitats of Xiapingdi monkey groups with GD and Tongjiashan monkey groups without GD (CK), and the soil samples were selected according to Figure 1. The specific sampling method was as follows: First, survey lines were set up in the area under GD and CK, and the starting and end points of the sample lines referred to Frontiers in Environmental Science frontiersin.org the results of Ran et al. (2003). After that, five large quadrats of 20 m × 20 m were set up in each sample line, and three soil samples in each large quadrat were collected randomly; then, the litter layer on the soil surface was removed. After that, a sieve with a pore size of 2 mm was used to remove roots and stones. Lastly, the treated soils were stored at −80°C for DNA extraction.

DNA extraction and sequencing
Frozen soil (2 g) was used to extract DNA using the FastDNA ®  (Xia et al., 2022). Shotgun metagenomic sequencing of DNA extracted from soil samples was performed simultaneously . In brief, approximately 1 μg of DNA was used to construct a library with a 300 bp insert size, followed by sequencing using the Zhu et al., 2021 platform (Majorbio Company in Shanghai) with the PE150 (paired end sequencing 150 × 2) strategy. The quality control method used on the raw data is described in detail in . There were 384443 reads of 10 soil samples. The read data were deposited into the NCBI Sequence Read Archive (SRA) database under accession number PRJNA906922. Following quality control, the clean reads were assembled using IDBAUD , and the assembled contigs longer than 500 bp were reserved. We obtained 7635863 contigs with an average length of 828.42 bp each from 10 metagenomic samples. The (ORFs) of the contigs were predicted using Prodigal v2.6.3 with a meta model. The coverage of each contig was calculated by mapping clean reads to the contigs using bbmap (https://sourceforge. net/projects/bbmap/) with the default parameters (Ju et al., 2019).

Identification of ARG-like ORFs and drug-resistant pathogens
The ORFs were identified against the CARD database (https:// card.mcmaster.ca/) using DIAMOND with an E-value ≤10 -10 (Buchfink et al., 2015), and the results that identified ≥80% were identified as ARG-like ORFs (Ma et al., 2015).
To compare coverage between different samples, the coverage of ARG-like ORFs was normalized using the data size of each sample (copies/Gb). The calculation formula was as follows: where n is the total number of ARG-like ORFs in one sample, N is the number of clean reads mapped to ARG-like ORFs, 150 is the length of clean reads, L is the length of the target ARG-like ORFs, and G is the data size (Gb) of clean reads per sample (Ma, et al., 2015;Xiong, et al., 2018). The ORFs were detected using VFanalyzer on VFDB in setB (http://www.mgc.ac.cn/VFs/) with an E-value ≤10 -10 to identify the virulence factors (VFs), and the results that identified ≥80% were identified as ARG-like ORFs. If ARG-like ORFs also had VFs, ORFs were drug resistance genes that could reflect the risk of drugresistant pathogens (Liu et al., 2018).

Risk of ARGs in transferring and spreading 2.4.1 Ability of MGEs to transfer ARGs
To estimate the potential horizontal gene transfer of ARGs in soil samples, the ORFs were searched against the MGE database (http:// aclame.ulb.ac.be/) using BLASTP with an e-value ≤10 -5 , and the results with a shared identity ≥80% were identified as MGEs. The ORF annotation results were used to identify the MGEs on plasmids, chromosomes and viruses. To compare the MGE abundance among different samples, the coverage of the MGE-like ORFs was also normalized by the data size of each sample (copies/Gb) according to Eq. 1. If the ARGlike ORFs also contained MGEs, the ARGs could be moved by the MGEs.

Ability of microorganisms to transfer ARGs
Network analysis was used to explore the potential correlations between bacteria and ARGs . Based on the hypothesis that a significant correlation between sOTUs and ARG subtypes was observed, this observation indicated that the bacterial community structure determined the distribution of the resistome and that the horizontal gene transfer of ARGs was not sufficient to obscure their association with bacterial genomes. The abundance and species of the host bacteria carrying ARGs could also be used to indicate the ability of microorganisms to transfer the ARGs. The higher the abundance and species of host bacteria carrying ARGs was, the stronger the ability of microorganisms to transfer ARGs.
The ORF sequences that carried ARGs were searched against the NCBI RefSeq database using BLASTP with an e-value ≤10 -5 . The results were parsed using MEGAN (version 5) (Huson et al., 2016), and the ORFs were annotated as the taxon if more than 50% of the ORFs were annotated as the same taxon; the taxons are the host bacteria of ARGs (Ishii et al., 2013).

Statistical analysis
The diversities and relative abundances of the soil microorganisms and their ARGs as well as the abundances of MGEs were their mean values ± standard error in GD or CK. The differences in the diversities and relative abundances of the soil microorganisms and their ARGs as well as the abundances of MGEs between GD and CK were analyzed by a t-test with a significance level of p < 0.05. The average abundances of Frontiers in Environmental Science frontiersin.org OTUs were higher than 0.1% in 10 soil samples, and the relative abundances of all types of ARGs were used to calculate the correlation coefficients (ρ) according to Spearman between soil microorganisms and ARGs using the vegan package in R4.2.1. If ρ ≥ 0.8 and p value ≤0.01, the relationship of OTUs and ARGs was used to perform a network analysis using the igraph package in R4.2.1 and then visualized in Gephi (version 0.9.5).
3 Results and discussion

Effects of GD on the characteristics of ARGs in soil microorganisms
Microbial diversity has a positive effect on the cycling of nutrients, such as the cycling of carbon and nitrogen (Zhang H et al., 2020). Simpson's index gives more weight to the more abundant species in a sample. It takes into account the number of species present, as well as the abundance of each species, while the ACE diversities shows more species present in a sample (Schmidt et al., 2019). Supplementary Figure S1 shows that the Simpson's diversities of soil microorganisms in GD were significantly higher than that of the CK (p = 0.032), while the ACE diversities in GD were significantly lower than that of the CK (p = 0.012), which indicated that GD reduced the total soil microbial species, but increased the dominance of dominant microbial species. GD reduced the total soil microbial species, which could be explained by the fact that GD reduces the diversity of plants and the content of humic substances in wildlife habitats , which results in the reduction of mineralization rates of C and N in soil and the decrease of available nutrients for microbial reproduction by GD (Wu et al., 2022). However, the relative abundances of dominant microorganisms were significantly increased by GD, they were Chloroflexi and Verrucomicrobia (p = 0.001 and p = 0.049), while

FIGURE 2
Effects of GD on characteristics of ARGs in soil microorganisms. Note: (A-C) represent the types of ARGs, the relative abundance of ARGs classes, and Simpson indexes of ARGs in GD and CK, respectively. (D) shows the relative abundances of ARGs with significant differences between GD and CK. * means that the Simpson indexes of ARGs or the relative abundance of the ARGs are significant difference between CK and GD (p < 0.05), "GD" means the ARGs only be checked in GD.

Frontiers in Environmental Science
frontiersin.org GD had no significant effect on the abundances of other microorganisms (p > 0.05), thus the soil Simpson's diversities were increased by GD, even resulted in a rise in the function of dominant microbial species in the soil of wild animal habitats.
Chloroflexi and Verrucomicrobia are associated with the carbon and nitrogen contents of the soil. Chloroflexi is a photosynthetic bacterium that can synthesize carbon and nitrogen (Singleton et al., 2020), while the abundance of Verrucomicrobia is also significantly positively correlated with the total carbon and nitrogen contents in soils (Freitas et al., 2012). Thus, our results indicated that the cycling nutrients in the soils of the wildlife habitat were reduced by microorganisms under GD (Supplementary Figure S1C). A total of 650 ORFs containing ARGs were detected in the soil microorganisms, and 291 types of ARGs were classified into 18 kinds of resistance classes in GD and CK (Figure 2A), which indicated that there was a potential risk from ARGs in the habitats of wild animals. Because ARGs in the surroundings of the habitat can affect the health of wild animal (Xu et al., 2020) and the accessibility of ARGs to animals is one of the factors that affect the risk from ARG . Further analysis showed there were no significant difference in the total relative abundances of ARGs between the GD and CK (p > 0.05, Supplementary Figure S2A). However, the Simpson indexes of ARGs were increased significantly by GD (p < 0.05, Supplementary Figure S2B). Simpson's diversity is a measure of diversity which takes into account both species and the relative abundance of the different species (Schmidt et al., 2019), thus combined with our analysis, we found two reasons for the increased the diversity of ARGs: First, there were 287 and 245 types of ARGs in GD and CK. Among them, 46 types of ARGs and triclosan and bicyclomycin were checked in GD but not in CK (Figures 2A, D), which showed that GD could increase types of ARGs. Second, the relative abundances of MLS and mupirocin in GD were higher than those in CK (p < 0.05, Figure 2B). The diversity of ARGs is higher, the risk from ARGs is greater . Thus, our result indicated that the risk of the ARGs was increased by increasing types and the abundances of ARGs under GD. ARGs in the surroundings of the habitat can affect the health of wild animal (Xu et al., 2020), as well as humans or other animals could promote the transfer of ARGs, and increases the relative abundance of ARGs Lu, et al., 2022). Thus, we concluded that the wild animals may have a risk of ARGs in acessibility of ARGs in the habitat, especially, the ARGs of MLS, mupirocin, triclosan and bicyclomycin that only be checked in GD but not in CK. MLS, mupirocin, triclosan and bicyclomycin, can be used to treat infections caused by a variety of sensitive bacteria (Patteson, et al., 2017). We speculated that if the wild animals who lives in the habitats under GD are infected by pathogenic bacteria, the risk of treatment to wild animals will increase, for the pathogenic bacteria have increased their resistance to the four classes ARGs of MLS, mupirocin, triclosan and bicyclomycin under GD (Figure 2).

Effects of GD on the horizontal transfer potential of ARGs in soil microorganisms
MGEs play an important role in the horizontal transfer of ARGs, the higher the abundance of MGEs is in microorganisms, the stronger the ability of ARGs to be carried by MGEs, and the

FIGURE 3
Effects of GD on the characteristics of MGEs in soil microorganisms. Note: (A, B) show the Simpson's indexes and relative abundances of MGEs of soil microorganisms, respectively. The same lowercase indicates there are no significant difference between CK and GD (p ≥ 0.05).** means that the relative abundance of the MGEs are significant difference between CK and GD (p < 0.01).

Frontiers in Environmental Science
frontiersin.org greater the risk of horizontal transfer of ARGs . Our results showed that there were 12160 MGE-like ORFs in both CK and GD, the relative abundances of MGEs in GD were higher than those in CK ( Figure 3B; p = 0.000162), which indicated that GD could increase the horizontal transfer ability of MGEs to ARGs. Gut microbiome could acquire exogenous ARGs via horizontal gene transfer (Mei et al., 2021). Table 1 also shows that there are ARG-like ORFs carried by plasmids only be checked in GD, namely, tetM, MexD, ceoB, MuxC, sul1 and catII, which indicated that if the carried ARGs are transferred to gut microbiome of wild animals under GD, the wild animals they will develop resistance to some antibiotics. The relative abundances of MGEs carried by plasmid, prophage and virus were as follows: plasmid>prophage>virus, and the relative abundances of plasmids in the GD were significantly higher than those in the CK ( Figure 3B; p = 0.000160), while there were no significant differences in the relative abundances of MGEs carried by prophages and virus between the CK and GD ( Figure 3B; p > 0.05). If both ARGs and MGEs were present in one ORF, the ARGs could be moved by MGEs (Ju et al., 2019). Table 1 shows that 55 ARG-like ORFs could be carried by plasmids, and the total number of ORFslike plasmids in GD was significantly higher than that in CK (Supplementary Figure S3; p = 0.034). Combined with Figure 3B; Table 1, we found that plasmids had the strongest ability to carry ARGs and GD increased their ability to carry ARGs in wildlife habitats. In nature, ARGs carried by plasmids are the main method by which ARGs are transferred horizontally, studies have shown that the ARGs in the microorganisms of water, activated sludge and soil are mainly carried by plasmids (Che et al., 2019;Fan et al., 2019;Li et al., 2020;, which supports our results.

Effects of GD on drug-resistant pathogens
If both an ARG and a virulence factor (VF) are present in one ORF, the bacteria in the ORF are called drug-resistant pathogens that reflected the highest level of risk of ARGs . There were 26 ARG-like ORFs that were detected in three types of pathogenic bacteria: Klebsiella pneumoniae subsp. pneumoniae NTUH-K2044, Acinetobacter baumannii ACICU and Neisseria meningitidis MC58. The total number of the 26 types of ORFs in GD was significantly higher than that in CK (p = 0.032, Supplementary Figure S4). Among them, six types of ORFs were detected only in GD but not in CK (Supplementary Table S1), they were found to be two types of pathogenic bacteria: Klebsiella pneumoniae NTUH-K2044 and Acinetobacter baumannii ACICU. Thus, not only the total abundance but also the types of drugresistant pathogens were increased by GD. Previous studies have shown that the wide application of animal manure could introduce a large number of ARGs and drug-resistant pathogen into the soil surroundings (Duan et al., 2017;Liu et al., 2017;, because antibiotic stress in animal manure under GD can not only increase the abundance and species of ARGs but also increase the abundance of VFs due to the adaptive evolution of pathogenic bacteria (Durrant et al., 2020). So, we speculate that once drug-

Effects of GD on the host bacteria of ARGs
Although these findings have demonstrated that the risk of ARGs was increased by GD in the habitats of wildlife, as the host, the microbial community determines the pattern and environmental behavior of ARGs, and the transfer of ARGs is also related to microorganisms (Bing et al., 2015). The co-occurrence network correlation between ARGs and bacterial taxa is regarded as an effective way to track potential ARG hosts in various environments (Hu et al., 2017). To investigate the association between ARGs and the bacterial community, a correlation analysis was performed between the relative abundance of ARGs and bacteria. There were 271 points and 520 edges in the CK, of which 56 edges were significantly negative and 462 edges were significantly positive (|r|>0.8, p < 0.01), while there were 305 points and 1,150 edges in the GD, of which 60 edges were significantly negative and 1,090 edges were significantly positive (|r|>0.8, p < 0.01; Figures 4A, B). If there is a significant correlation between the abundance of ARGs and the abundance of microorganisms, the microorganisms can be considered potential hosts of ARGs (Zhu et al., 2018), and the network relationship between ARGs and microorganisms can explain the transfer potential of Host bacteria of ARGs and relationships between ARGs and soil microorganisms. Note: (A, B) were the network of spearman correlation relationships between ARGs calasses and OTUs in soil microorganisms in CK and GD, while (C, D) were the host bacteria of ARGs in CK and GD, respectively. The red lines indicated that there were significant positive correlations between the relative abundances of ARGs and the relative abundances of dominant OTUs (p < 0.01; r > 0.8). The green line indicated that the relative abundances of ARGs were significantly negatively correlated with the relative abundances of dominant OTUs (p < 0.01; r < −0.8).

Frontiers in Environmental Science
frontiersin.org microorganisms to ARGs in the environment to some extent . Thus, our results indicated that the bacterial community structure determined the distribution of the resistome and that the horizontal transfer of ARGs was not sufficient to obscure their association with bacterial genomes. The species and numbers of the host bacteria of ARGs were increased by GD, which indicated that the carrying capacity of ARGs in soil microorganisms to move was increased by GD. Although the network relationship results indicated a significant correlation between ARGs and bacteria in the habitats of wild animals, they could not reflect the ability of microorganisms to carry ARGs. If a microorganism and ARG are on one same ORF, the microorganism can carry the ARGs to transfer (Zeng et al., 2019). Figure 4 shows that there were 232 types of ARG-like ORFs that were detected in seven types of bacteria, namely, Bacteroidetes, Firmicutes, Actinobacteria, Proteobacteria, Spirochaetes, Tenericutes and unclassified bacteria, they carried 5.70%, 17.34%, 8.08%, 54.63%, 23.75%, 23.75% and 13.78% of the ARG-like ORFs in the CK, respectively. While there were 242 types of ARG-like ORFs that were detected in eight types of bacteria, namely, Actinobacteria, Bacteroidetes, Deferribacteres, Firmicutes, Proteobacteria, Spirochaetes, Tenericutes and unclassified bacteria, they carried 6.02%, 4.47%, 1.94%, 14.56%, 63.69%, 1.94%, 1.94% and 10.87% of the ARG-like ORFs in the GD, respectively. Among the ARG-like ORFs, 18 types and 16 antibiotic classes could be carried by bacteria in GD and CK, respectively, and triclosan and bicyclomycin (79 types of ARG-like ORFs were not detected in the CK) were not carried in CK. Results have showed that the bacterial hosts of ARGs were preferably distributed in Proteobacteria (Zeng et al., 2019), which was consistent with our research results that Proteobacteria was the main host of ARGs in the habitats of wild are no significant difference between CK and GD (p ≥ 0.05), while the same lowercase indicates there are significant difference between CK and GD (p < 0.05). "**" and "**" means there are significant difference between CK group and GD on 0.01 level and 0.05 level.

Frontiers in Environmental Science
frontiersin.org animals. Proteobacteria is one of the largest groups of bacteria in nature, and it includes many pathogenic bacteria (Rizzatti et al., 2017). In our results, the highest proportion (98.73%) of the host of the new ARGs in GD was Proteobacteria, and GD increased the ability of Proteobacteria to carry ARGs ( Figures 4A, B). We concluded that the ecological risk of ARGs carried by Proteobacteria in the habitats of wild animals is the largest, which was increased by GD, as supported by the results of Supplementary Table S2 because Supplementary  Table S2 shows that all three types of drug-resistant pathogens were Proteobacteria, including the two types of new drug-resistant pathogens in GD.
3.5 Effects of GD on the risk of multiresistance genes (MRGs) Figure 5 showed that the Simpson's diversity of the MRGs was higher than that of the other ARG classes, and GD increased the Simpson's diversity of the MRGs, as well as the MRGs accounted for the largest proportion in the ARG-like ORFs that were only checked in the GD ( Figures 5A, B), which indicated that the risk of MRGs was the greatest, and GD increased this risk. MRGs was the most abundant and types and human disturbance increases the risk of MRGs in soil samples (Mei et al., 2021), which supported our result. Our result showed that the MRGs-like ORFs were the largest proportion that could be transferred by plasmids and they were increased by GD ( Figure 5C; p = 0.007), the largest proportions of ARGs that the host bacteria carried were MRGs, GD increased this proportions ( Figures 4C, D), as well as the MRGs had the maximum distribution degree in the network, GD increased their degree distribution, (Supplementary Table S2). Our results indicated that GD increased mobility of MRG and the storage capacity of MRGs in multidrug-resistant bacteria, which could be explained by this reason that multidrug-resistant bacteria are mainly generated by transfer of MRGs in MGEs and the accumulating MRGs in microorganisms (Kim et al., 2021). The risk of MRGs favors the development of multidrug-resistant bacteria through horizontal gene transfer (Kang et al., 2022), thus we concluded that GD increased the risk of MRGs by increasing the mobility of MRGs and the storage capacity of microorganisms to MRGs.
MRGs present in microorganisms, which could reflect the risk level of ARGs in a region (Zeng et al., 2019). Today, infection caused by multidrug-resistant bacteria is a major threat to ecosystem health (Fiona and Brion, 2013). The bad thing is that, GD increased the abundance of multidrug-resistant pathogens ( Figure 5D; p = 0.033), and the multidrug-resistant pathogens-like ORFs of Klebsiella pneumoniae NTUH-K2044 only check in GD. As early as 2008, result had shown that Acinetobacter baumannii ACICU could cause a variety of bacterial infectious diseases, and it had resistance to some (Iacono et al., 2008). The propensity of this bacterium to rapidly acquire antibiotic resistance could lead to the emergence and spread of Acinetobacter baumannii strains with MRGs, and could develop resistance to latest antibiotics (Law and Tan, 2022), which supported our findings that the ARGs in Acinetobacter baumannii ACICU are MRGs, it could develop resistance to a variety of antibiotics (Supplementary Table S1). More seriously, there was one MGE in the Acinetobacter baumannii ACICU under GD, it was pGMI1000MP, which indicated that the Acinetobacter baumannii ACICU could be moved by plasmid. Acinetobacter baumannii is a pathogen responsible for several serious infections, including pneu-monia, sepsis, and meningitis (Law and Tan, 2022). Thus, we concluded that GD increased the risk of MRGs, multidrugresistant pathogens may appear in the habitats of wild animals after GD, there is also a risk of horizontal transfer of drug-resistant pathogens to wild animals, which has a profound impact on survival of wild animals.

Conclusion
The occurrence, characteristics and transfer risk of soil microorganism ARGs were studied under GD, which demonstrated that ARGs could be essential indicators for the risk of GD in the habitats of wild animals since GD had great influences on the characteristics, mobility, pathogenicity and microbial storage capacity of ARGs in the habitats of wild animals. GD increased the diversity of ARGs by increasing new types of ARGs, GD increased the carrying capacity of MGEs to ARGs by increasing the species and abundance of the MGEs, and the ARGs in the habitats of wild animals could only be transferred horizontally through plasmids. GD increased the risk of host bacteria pathogens by increasing the variety and abundance of drugresistant pathogens, and Acinetobacter baumannii ACICU was the new drug-resistant pathogen after GD. GD increased the storage capacity of microorganisms, especially the dominant microorganism storage of ARGs, rather than increasing the microbial diversity in the soil. All three types of drug-resistant pathogens, Klebsiella pneumoniae subsp. pneumoniae NTUH-K2044, Acinetobacter baumannii ACICU and Neisseria meningitidis MC58, were in the Proteobacteria, and the ecological risk of ARGs carried by Proteobacteria in the habitats of wild animals was the largest, which was increased by GD. The risk of MRGs was the highest, and GD increased this risk. Moreover, multidrugresistant pathogens may have a risk of transfer to wild animals. We worry that ARGs may spread widely in wildlife or be transferred to the gut microbiome of wildlife, which will use different resistance mechanisms to develop resistance to different antibiotics and increase the risk of disease treatment in wildlife under GD. We suggest that to protect the health of wild animals and their habitats, it is necessary to control the scale of grazing and monitor drug-resistant pathogens from GD.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions
Conceptualization, DL; Formal analysis, SZ; Funding acquisition, SZ; Investigation, TL and CH; Supervision, DL; Writing-original draft, SZ and TL; Writing-review and editing, DL. All authors have read and agreed to the published version of the manuscript.

Funding
This work was supported by National Natural Science Foundation of China (32200403 and 32270548), and National