<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1238708</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Metagenomics reveals the abundance and accumulation trend of antibiotic resistance gene profile under long-term no tillage in a rainfed agroecosystem</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Weiyan</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
<xref rid="aff3" ref-type="aff"><sup>3</sup></xref>
<xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1418299/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shen</surname>
<given-names>Pengfei</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Lu</surname>
<given-names>Zhiqiang</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mo</surname>
<given-names>Fei</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liao</surname>
<given-names>Yuncheng</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wen</surname>
<given-names>Xiaoxia</given-names>
</name>
<xref rid="aff1" ref-type="aff"><sup>1</sup></xref>
<xref rid="aff2" ref-type="aff"><sup>2</sup></xref>
<xref rid="c002" ref-type="corresp"><sup>&#x002A;</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>College of Agronomy, Northwest A&#x0026;F University</institution>, <addr-line>Yangling, Shaanxi</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Key Laboratory of Low-carbon Green Agriculture in Northwestern China, Ministry of Agriculture and Rural Affairs</institution>, <addr-line>Yangling, Shaanxi</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Key Laboratory of Crop Physi-ecology and Tillage Science in Northwestern Loess Plateau, Ministry of Agriculture, Northwest A&#x0026;F University</institution>, <addr-line>Yangling, Shaanxi</addr-line>, <country>China</country></aff>
<author-notes>
<fn id="fn0010" fn-type="edited-by"><p>Edited by: Zhigang Qiu, Tianjin Institute of Environmental and Operational Medicine, China</p></fn>
<fn id="fn0011" fn-type="edited-by"><p>Reviewed by: Zheng-Rong Kan, Nanjing Agricultural University, China; Wensheng Fang, Chinese Academy of Agricultural Sciences (CAAS), China; Yilai Lou, Chinese Academy of Agricultural Sciences, China</p></fn>
<corresp id="c001">&#x002A;Correspondence: Weiyan Wang, <email>wweiyan@nwsuaf.edu.cn</email></corresp>
<corresp id="c002">Xiaoxia Wen, <email>wenxiaoxia6811@163.com</email>; <email>wenxx66@nwsuaf.edu.cn</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>20</day>
<month>07</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1238708</elocation-id>
<history>
<date date-type="received">
<day>12</day>
<month>06</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>07</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Wang, Shen, Lu, Mo, Liao and Wen.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Wang, Shen, Lu, Mo, Liao and Wen</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Widespread soil resistance can seriously endanger sustainable food production and soil health. Conservation tillage is a promising practice for improving soil structure and health. However, the impact of long-term no-tillage on the presence of antibiotic resistance genes in agricultural soils remains unexplored. Based on the long-term (&#x003E;11&#x2009;yr) tillage experimental fields that include both conservation tillage practices [no tillage (ZT)] and conventional tillage practices [plough tillage (PT)], we investigated the accumulation trend of antibiotic resistance genes (ARGs) in farmland soils under long-term no-tillage conditions. We aimed to provide a scientific basis for formulating agricultural production strategies to promote ecological environment safety and human health. In comparison to PT, ZT led to a considerable reduction in the relative abundance of both antibiotic resistance genes and antibiotic target gene families in the soil. Furthermore, the abundance of all ARGs were considerably lower in the ZT soil. The classification of drug resistance showed that ZT substantially decreased the relative abundance of Ethambutol (59.97%), &#x03B2;-lactams (44.87%), Fosfomycin (35.82%), Sulfonamides (34.64%), Polymyxins (33.67%), MLSB (32.78%), Chloramphenicol (28.57%), Multi-drug resistance (26.22%), Efflux pump (23.46%), Aminoglycosides (16.79%), Trimethoprim (13.21%), Isoniazid (11.34%), Fluoroquinolone (6.21%) resistance genes, compared to PT soil. In addition, the abundance of the bacterial phyla Proteobacteria, Actinobacteria, Acidobacteria, and Gemmatimonadetes decreased considerably. The Mantel test indicated that long-term ZT practices substantially increased the abundance of beneficial microbial flora and inhibited the enrichment of ARGs in soil by improving soil microbial diversity, metabolic activity, increasing SOC, TN, and available Zn, and decreasing pH. Overall, long-term no-tillage practices inhibit the accumulation of antibiotic resistance genes in farmland soil, which is a promising agricultural management measure to reduce the accumulation risk of soil ARGs.</p>
</abstract>
<kwd-group>
<kwd>conservation tillage</kwd>
<kwd>antibiotic resistance genes</kwd>
<kwd>soil drug resistance</kwd>
<kwd>microbial metabolic activity</kwd>
<kwd>metagenomic sequencing</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="4"/>
<equation-count count="0"/>
<ref-count count="59"/>
<page-count count="12"/>
<word-count count="8204"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Antimicrobials, Resistance and Chemotherapy</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="sec1" sec-type="intro">
<label>1.</label>
<title>Introduction</title>
<p>Antibiotic resistance genes (ARGs), a new type of environmental pollutant, are widely enriched in farmland soils, promoting the spread of drug resistance and seriously affecting soil safety (<xref ref-type="bibr" rid="ref17">Huijbers et al., 2015</xref>; <xref ref-type="bibr" rid="ref35">Rosenkrantz et al., 2019</xref>). The soil&#x2013;plant system is a significant source of exposure to environmental resistance genes, which can spread from the environment to humans (<xref ref-type="bibr" rid="ref9001">Wang et al., 2022a</xref>,<xref ref-type="bibr" rid="ref47">b</xref>). The food chain is the most direct and main pathway for resistant bacteria and genes to enter the human body in soil&#x2013;plant systems (<xref ref-type="bibr" rid="ref3">Blau et al., 2018</xref>). With the rapid development of the global economy, resistance genes in agricultural products can spread rapidly worldwide through food processing, preservation, and global transportation, triggering public safe risks and affecting global health (<xref ref-type="bibr" rid="ref28">Nowakowska et al., 2019</xref>; <xref ref-type="bibr" rid="ref58">Zhu et al., 2019</xref>). Therefore, to remove and control ARGs in farmland soils, it is particularly important to clarify the influence of soil management practices on ARGs and their relationship with microbial communities.</p>
<p>Cropping system soil, one of the habitats with the highest microbial diversity on Earth, is a huge reservoir of ARGs (<xref ref-type="bibr" rid="ref12">Forsberg et al., 2012</xref>; <xref ref-type="bibr" rid="ref37">Sanderson et al., 2016</xref>; <xref ref-type="bibr" rid="ref43">Van Goethem et al., 2018</xref>). Internal resistance and exogenous input are the main sources of ARGs in agricultural soil (<xref ref-type="bibr" rid="ref34">Riber et al., 2014</xref>; <xref ref-type="bibr" rid="ref27">N&#x00F5;lvak et al., 2016</xref>). Internal resistance refers to the inherent genome of microorganisms, that can be inherited by subsequent generations. Exogenous input involves the horizontal transfer and enrichment of ARGs directly into the soil through the addition of external sources, such as manure, sludge, and reclaimed water. The composition and abundance of ARGs are closely linked to the microbial community and can be affected by changes in microbial community diversity and composition (<xref ref-type="bibr" rid="ref11">Forsberg et al., 2014</xref>; <xref ref-type="bibr" rid="ref24">Liu et al., 2019</xref>; <xref ref-type="bibr" rid="ref54">Zhang et al., 2021</xref>). Studies have shown that <italic>Actinobacteria</italic> produce the most natural antibiotics, and <italic>Firmicutes</italic> and <italic>Actinobacteria</italic> have been identified as the main groups carrying ARGs (<xref ref-type="bibr" rid="ref57">Zhou et al., 2015</xref>; <xref ref-type="bibr" rid="ref52">Xun et al., 2016</xref>). <xref ref-type="bibr" rid="ref6">Chen et al. (2016)</xref> found that the abundance of tetracyclines ARGs increased as the abundance of <italic>Actinomycete</italic> increased and the abundance of <italic>Chloroflexi</italic> decreased. A study have also shown that the increased abundance of Proteobacteria was accompanied by increasing &#x03B2;-lactam-resistance genes (<xref ref-type="bibr" rid="ref46">Wang et al., 2016</xref>). In addition, <xref ref-type="bibr" rid="ref50">Xiang et al. (2018)</xref> and <xref ref-type="bibr" rid="ref42">Urra et al. (2019)</xref> showed that the decrease of soil microbial &#x03B1;-diversity led to an increase in MLSB, aminoglycosides, &#x03B2;-lactams, amide alcohols, sulfonamides ARGs and transposase gene <italic>tnpA</italic>. In summary, previous studies mainly focuses on manure application, sewage-irrigated, and typically contaminated soils (<xref ref-type="bibr" rid="ref1">Achiba et al., 2010</xref>; <xref ref-type="bibr" rid="ref51">Xie et al., 2018</xref>), however, there are no reports on the influence mechanism of tillage management measures on agricultural soil ARGs.</p>
<p>The introduction and enrichment of ARGs in farmland soils are commonly attributed to agricultural practices (<xref ref-type="bibr" rid="ref49">Wu et al., 2022</xref>). For example, manure application can considerably enrich resistant bacteria and ARGs in agricultural soils (<xref ref-type="bibr" rid="ref26">Munk et al., 2018</xref>; <xref ref-type="bibr" rid="ref51">Xie et al., 2018</xref>). Fertilizer application can significantly increase the abundance of important resistance genes, including <italic>intI1</italic>, <italic>tetW</italic>, <italic>tetO</italic>, <italic>sul1</italic>, and <italic>sul2</italic> (<xref ref-type="bibr" rid="ref51">Xie et al., 2018</xref>). <xref ref-type="bibr" rid="ref56">Zhao et al. (2019)</xref> found that the abundance of tetracycline- and ampicillin-resistant bacteria in greenhouse soils was substantially higher than that in field soils. Besides the change of soil physical and chemical properties (i.e., pH, soil organic carbon (SOC), total nitrogen (TN), NH<sub>4</sub><sup>+</sup>, NO<sub>3</sub><sup>&#x2212;</sup>, soil texture, etc) may create different niches for bacterial populations and alter the microbial communities, thereby affecting the composition and abundance of ARGs (<xref ref-type="bibr" rid="ref9">Dumbrell et al., 2010</xref>). In addition, <xref ref-type="bibr" rid="ref55">Zhao et al. (2017)</xref> found that the concentration of some heavy metals (such as Cu, Zn, Mn, Fe) in soil may also drive the distribution and spread of ARGs in the environment. Conservation tillage, that is, reduced or no tillage with straw return, can decrease soil disturbance and improve its physical and chemical properties, thereby leading to changes in microbial habitats and communities, consequently affecting likely the abundance of ARGs in agricultural soil (<xref ref-type="bibr" rid="ref45">Wang et al., 2017</xref>, <xref ref-type="bibr" rid="ref44">2021</xref>). This practice has been widely used in dry farmlands in northern China (<xref ref-type="bibr" rid="ref48">Wang et al., 2020</xref>). However, the impact of long-term no-tillage on ARGs in agricultural soil remains unknown, hampering a complete understanding of ARGs in this context.</p>
<p>Therefore, based on a long-term (&#x003E;&#x2009;11&#x2009;yr) conservation tillage experimental field, we investigated the effect on the abundance of ARGs and microbial community composition using shotgun metagenomic sequencing under two different tillage systems, no tillage and plough tillage, in a winter wheat&#x2013;summer maize rotation system. We hypothesized that long-term no-till application could reduce the abundance and accumulation of ARGs by improving soil properties and microbial community structure. Therefore, three main objectives were established: (1) to analyze the effects of different tillage treatments on the relative abundance of ARGs; (2) to study the contribution of the microbial community to ARGs under long-term no tillage vs. plough tillage; and (3) to determine the correlations between soil variables and changes in ARGs under different tillage treatments.</p>
</sec>
<sec id="sec2" sec-type="materials|methods">
<label>2.</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1.</label>
<title>Experimental site and design</title>
<p>This study was conducted at the conservation tillage long-term experimental field (34&#x00B0;16&#x2032;48&#x2033; N and 108&#x00B0;4&#x2032;12&#x2033; E) in the north campus of Northwest A&#x0026;F University, Yangling, Shaanxi, China. The site is located at the southern part of the Loess Plateau and has a temperate monsoon climate, with high summer temperatures and cold winters, annual average temperature of 12.9&#x00B0;C, and the frost-free period of 221&#x2009;d. The annual rainfall is approximately 620&#x2009;mm, and 60&#x2013;70% of the rainfall occurs between July and October. The soil texture was classified as silty clay (<xref ref-type="bibr" rid="ref38">Soil Survey Staff, 2010</xref>). The initial soil properties at a depth of 0&#x2013;20&#x2009;cm was reported in a previous study (<xref ref-type="bibr" rid="ref44">Wang et al., 2021</xref>). This is the main winter wheat-summer maize double-cropping area in China.</p>
<p>The long-term tillage experiment in this study began in 2009, with a total of two conservation tillage treatments [i.e., no tillage with straw return (ZT) and conventional tillage plough with straw return (PT)]. In this study, the ZT treatment, the most representative conservation tillage practice, was selected to explore the effect of long-term conservation tillage on ARGs and the driving mechanism that ensures summer maize growth. A randomized block design was adopted with three replicates for each treatment. The plot area was 125&#x2009;m<sup>2</sup> (25&#x2009;&#x00D7;&#x2009;5&#x2009;m). Tillage operations were performed before sowing winter wheat and summer maize. The specific tillage steps were as follows: in ZT, to minimize soil disturbance and create a good seedling bed to ensure seedling emergence, a rotary machine was used for shallow rotation with a depth of 5&#x2009;cm, and then the maize precision seeder was used for sowing; in PT, the plot was ploughed to a depth of 20&#x2013;30&#x2009;cm by a mouldboard plough machine, the seedbed was arranged by a rotary tiller, and the maize precision seeder was used for sowing.</p>
<p>Summer maize (cultivar Shandan 609) was sown on 17 June 2018 and harvested on 5 October 2018. The row spacing and plant spacing of summer maize plant were 70&#x2009;cm and 23&#x2009;cm, respectively. Wheat straw was crushed in pieces less than 5&#x2009;cm in length and evenly distributed on the ground after harvest and then tilled and sown. The fertilizers were applied as basal fertilizers; which involved spreading them on the soil surface by hand before tillage. The amount of urea (N, 46%) and diammonium phosphate (N, 18%; P<sub>2</sub>O<sub>5</sub>, 46%) applied was 375&#x2009;kg hm<sup>&#x2212;1</sup>.</p>
</sec>
<sec id="sec4">
<label>2.2.</label>
<title>Soil sampling and physicochemical properties analysis</title>
<p>During the 2018 summer maize growing season, soil samples were collected on the first day after tillage, seedling stage, and silking stage under two tillage treatments, with three replicates. Using the S-type sampling method, five soil cores (0&#x2013;20&#x2009;cm depth) were collected from each plot to form a sample, and large rocks and roots were screened using a 2&#x2009;mm mesh. Representative samples from each plot (0&#x2013;20&#x2009;cm depth) were divided into two parts: one was stored at &#x2212;80&#x00B0;C for metagenomic sequencing; and the other was stored at &#x2212;20&#x00B0;C for determining soil properties and available metal ion concentrations. In addition, the soil bulk density of 0&#x2013;10&#x2009;cm and 10&#x2013;20&#x2009;cm was measured by ring knife method.</p>
<p>The soil particle size distribution was analyzed using a Malvern laser particle size analyzer (Malvern Ms 2000, Malvern Instruments Ltd., United Kingdom). The moisture content of the soil samples was determined using the gravimetric method, and soil temperature was determined using a digital thermometer sensor Prober (EC-TP101, Henglu Biological Technology Co., Ltd., China). Soil pH at a water-soil ratio of 1:5 was measured using an electronic pH detector (CyberScan pH 510, Thermo Fisher Scientific, United States). Soil inorganic nitrogen was extracted by 1&#x2009;mol&#x2009;L<sup>&#x2212;1</sup> KCl solution, and then the concentration of NO<sub>3</sub><sup>&#x2212;</sup> and NH<sub>4</sub><sup>+</sup> was determined using an AA3 HR Autoanalyzer (Seal Analytical GmbH, Germany). Soil total nitrogen (TN) content was measured using the Kjeldahl digestion method (<xref ref-type="bibr" rid="ref39">Sparks, 1996</xref>). Soil organic carbon (SOC) content was determined using the potassium dichromate-concentrated sulfuric acid plus heat capacity method (<xref ref-type="bibr" rid="ref44">Wang et al., 2021</xref>). Microbial biomass carbon (MBC) was calculated from the difference of SOC between chloroform fumigated samples and control samples (<xref ref-type="bibr" rid="ref33">Ren et al., 2021</xref>). Microbial respiration (MR: CO<sub>2</sub> flux) was measured by alkali absorption. The metabolic quotient (qCO<sub>2</sub>)&#x2009;=&#x2009;MR/MBC, and microbial quotient (Cmic:Corg)&#x2009;=&#x2009;MBC/SOC (<xref ref-type="bibr" rid="ref33">Ren et al., 2021</xref>). The concentrations of available iron (Fe), manganese (Mn), copper (Cu) and zinc (Zn) were measured by DTPA leaching atomic absorption spectrophotometry method (<xref ref-type="bibr" rid="ref1">Achiba et al., 2010</xref>).</p>
</sec>
<sec id="sec5">
<label>2.3.</label>
<title>DNA extraction</title>
<p>Genomic DNA of soil microorganisms was extracted the FastDNA<sup>&#x00AE;</sup> SPIN Kit for Soil (MP Biomedicals Co., Ltd., Santa Ana, CA, United States) as per the manufacturer&#x2019;s instructions. Each soil sample was extracted four times and mixed. DNA concentration and purification were determined using a TBS-380 microfluorometer. DNA integrity was detected using 1% agarose gel electrophoresis.</p>
</sec>
<sec id="sec6">
<label>2.4.</label>
<title>Shotgun metagenomic sequencing, assembly, and ARGs annotation</title>
<p>Six soil samples were collected at the seedling stage (two tillage practices&#x2009;&#x00D7;&#x2009;with three replicates) and selected for shotgun metagenomic sequencing. M220 Focused-ultrasonicator (Covaris M220, Gene Company Limited, Hong Kong, China) was used to break the microbial genomic DNA into 400&#x2009;bp fragments. The PE library was constructed using a NEXTFLEX Rapid DNA-Seq Library Building Kit (BioScientific, United States). After bridge PCR amplification, libraries were sequenced using an Illumina HiSeq 3000/4000 platform (Illumina, United States) at Majorbio Bio-Pharm Technology Co., Ltd. Approximately 120&#x2009;GB original sequences were obtained using shotgun metagenome sequencing with 292,349,780 paired reads. Raw reads were uploaded to the NCBI for Biotechnology Information Sequence Reading Archive database (accession no. PRJNA641749).</p>
<p>The adaptors were stripped using SeqPrep.<xref rid="fn0001" ref-type="fn"><sup>1</sup></xref> Raw sequencing reads were subjected to a series of quality controls using fastp<xref rid="fn0002" ref-type="fn"><sup>2</sup></xref>, in which reads with a length&#x2009;&#x003C;&#x2009;50&#x2009;bp, quality value &#x003C;&#x2009;20, or having N bases were filtered and removed. IDBA-UD/Megahit<xref rid="fn0003" ref-type="fn"><sup>3</sup></xref> and Newbler<xref rid="fn0004" ref-type="fn"><sup>4</sup></xref> were used to maximize the use of high-quality clean reads for multiple hybrid splicing and assembly (<xref ref-type="bibr" rid="ref30">Peng et al., 2012</xref>; <xref ref-type="bibr" rid="ref19">Li D. H. et al., 2015</xref>), contigs&#x2009;&#x003E;&#x2009;500&#x2009;bp in length were retained for gene prediction. MetaGene<xref rid="fn0005" ref-type="fn"><sup>5</sup></xref> was used to predict the ORF of the retained contigs. Genes with nucleic acid lengths &#x003E;100&#x2009;bp were selected and translated into amino acid sequences to obtain gene prediction tables for each sample. Subsequently, using SOAPaligner software<xref rid="fn0006" ref-type="fn"><sup>6</sup></xref>, high-quality reads and non-redundant gene catalog of each sample were compared (95% identity) to obtain non-redundant gene catalog for subsequent analysis (<xref ref-type="bibr" rid="ref13">Fu et al., 2012</xref>).</p>
<p>The comprehensive antibiotic resistance database<xref rid="fn0007" ref-type="fn"><sup>7</sup></xref> contains reference genes related to antibiotic resistance in various organisms, genomes, and plasmids that can be used to guide research on drug and antibiotic resistance mechanisms in environmental flora (<xref ref-type="bibr" rid="ref25">Mart&#x00ED;nez et al., 2007</xref>; <xref ref-type="bibr" rid="ref18">Jia et al., 2016</xref>). To obtain antibiotic resistance genes annotation tables, DIAMOND software<xref rid="fn0008" ref-type="fn"><sup>8</sup></xref> was used to compare non-redundant gene catalog with CARD database (<italic>e</italic>-value &#x2264;1e-5) (<xref ref-type="bibr" rid="ref4">Buchfink et al., 2021</xref>). In this step, we performed functional screening at the Class level of the CARD database, in which we constructed gene subsets of the antibiotic biosynthesis (AB), antibiotic resistance (AR), antibiotic sensitive (AS), and antibiotic target (AT) families, respectively, for subsequent analysis. Then, to assess the abundance of these genes, high-quality sequences from each sample were mapped onto the predicted gene sequences using Salmon (<xref ref-type="bibr" rid="ref29">Patro et al., 2017</xref>) and the transcripts per kilobase per million mapped reads (TPM) were used to normalize the abundance values in the metagenomes. The data were analyzed on the online tool of Majorbio Cloud Platform<xref rid="fn0009" ref-type="fn"><sup>9</sup></xref> (<xref ref-type="bibr" rid="ref32">Ren et al., 2022</xref>).</p>
</sec>
<sec id="sec7">
<label>2.5.</label>
<title>Statistical analysis</title>
<p>Microsoft Office Professional Plus 2021 software (Microsoft Co., Redmond, WA, United States) was used for data collection and calculations. One-way ANOVA was performed using the SPSS software (version 19.0; SPSS Inc., Armonk, NY, United States). The significant differences between the treatments were tested by Tukey&#x2019;s HSD test (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). Nonmetric multi-dimensional scaling (NMDS) and redundancy analysis (RDA) were performed using &#x201C;vegan&#x201D; packages in the R 3.6.1 for Windows (<xref ref-type="bibr" rid="ref31">R Core Team, 2022</xref>). Comparative analysis of the two groups were conducted to analyze ARGs with significant difference between ZT and PT using the &#x201C;stats&#x201D; package in the R 3.6.1 for Windows. Furthermore, based on Pearson correlation analysis and Mantel test, the relationship between antibiotic resistance genes carrying community and microbial metabolic coefficient and soil properties was analyzed using &#x201C;cor&#x201D; function and &#x201C;vegan&#x201D; package in the R 3.6.1 for Windows.</p>
</sec>
</sec>
<sec id="sec8" sec-type="results">
<label>3.</label>
<title>Results</title>
<sec id="sec9">
<label>3.1.</label>
<title>Long-term tillage practice contributes to differences in soil characteristics</title>
<p>Long-term ZT significantly altered soil physical properties, chemical properties, and available soil concentrations of Fe, Mn, Cu, and Zn (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). Compared with the PT treatment, the ZT treatment greatly increased the fraction of sand and silt, soil moisture, soil bulk in the 0&#x2013;10&#x2009;cm layer, SOC, TN, NO<sub>3</sub><sup>&#x2212;</sup>, and NH<sub>4</sub><sup>+</sup>. Conversely, the clay fraction, pH, and soil bulk (10&#x2013;20&#x2009;cm layer) were significantly lower in ZT soil than that in PT soil (<xref rid="tab1" ref-type="table">Table 1</xref>). In addition, the available concentrations of Fe, Mn, Cu, and Zn were enriched in the ZT soil for up to a decade and were considerably higher than those in the PT soil (<xref rid="tab2" ref-type="table">Table 2</xref>).</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption><p>The physical and chemical properties of topsoil, as affected by the tillage practices.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="middle" rowspan="2">Tillage</th>
<th align="center" valign="middle" colspan="11">Physical and chemical properties of topsoil</th>
</tr>
<tr>
<th align="center" valign="middle">Sand (%)</th>
<th align="center" valign="middle">Slit (%)</th>
<th align="center" valign="middle">Clay (%)</th>
<th align="center" valign="middle">pH</th>
<th align="center" valign="middle">TN (g&#x2009;kg<sup>&#x2212;1</sup>)</th>
<th align="center" valign="middle">SOC (g&#x2009;kg<sup>&#x2212;1</sup>)</th>
<th align="center" valign="middle">Moisture (%)</th>
<th align="center" valign="middle">NO<sub>3</sub><sup>&#x2212;</sup> (mg&#x2009;kg<sup>&#x2212;1</sup>)</th>
<th align="center" valign="middle">NH<sub>4</sub><sup>+</sup> (mg&#x2009;kg<sup>&#x2212;1</sup>)</th>
<th align="center" valign="middle">B10 (g&#x2009;cm<sup>&#x2212;3</sup>)</th>
<th align="center" valign="middle">B20 (g&#x2009;cm<sup>&#x2212;3</sup>)</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top" colspan="12">The first day after tillage</td>
</tr>
<tr>
<td align="left" valign="top">Zero-till</td>
<td align="center" valign="middle">0.97&#x2009;&#x00B1;&#x2009;0.02a</td>
<td align="center" valign="middle">64.72&#x2009;&#x00B1;&#x2009;0.10a</td>
<td align="center" valign="middle">34.31&#x2009;&#x00B1;&#x2009;0.12b</td>
<td align="center" valign="middle">7.60&#x2009;&#x00B1;&#x2009;0.04 b</td>
<td align="center" valign="middle">1.35&#x2009;&#x00B1;&#x2009;0.05a</td>
<td align="center" valign="middle">13.41&#x2009;&#x00B1;&#x2009;0.96a</td>
<td align="center" valign="middle">20.76&#x2009;&#x00B1;&#x2009;0.52a</td>
<td align="center" valign="middle">2.93&#x2009;&#x00B1;&#x2009;0.61a</td>
<td align="center" valign="middle">26.19&#x2009;&#x00B1;&#x2009;1.94a</td>
<td align="center" valign="middle">1.30&#x2009;&#x00B1;&#x2009;0.02a</td>
<td align="center" valign="middle">1.30&#x2009;&#x00B1;&#x2009;0.02b</td>
</tr>
<tr>
<td align="left" valign="top">Plow&#x2013;till</td>
<td align="center" valign="middle">0.88&#x2009;&#x00B1;&#x2009;0.01b</td>
<td align="center" valign="middle">63.45&#x2009;&#x00B1;&#x2009;0.056b</td>
<td align="center" valign="middle">35.68&#x2009;&#x00B1;&#x2009;0.07a</td>
<td align="center" valign="middle">7.89&#x2009;&#x00B1;&#x2009;0.03 a</td>
<td align="center" valign="middle">1.17&#x2009;&#x00B1;&#x2009;0.03b</td>
<td align="center" valign="middle">10.44&#x2009;&#x00B1;&#x2009;0.46b</td>
<td align="center" valign="middle">18.39&#x2009;&#x00B1;&#x2009;1.21b</td>
<td align="center" valign="middle">0.3&#x2009;&#x00B1;&#x2009;0.32b</td>
<td align="center" valign="middle">9.95&#x2009;&#x00B1;&#x2009;2.65b</td>
<td align="center" valign="middle">1.26&#x2009;&#x00B1;&#x2009;0.02b</td>
<td align="center" valign="middle">1.53&#x2009;&#x00B1;&#x2009;0.083a</td>
</tr>
<tr>
<td align="left" valign="top" colspan="12">The seedling stage (Metagenomic sequencing period)</td>
</tr>
<tr>
<td align="left" valign="middle">Zero-till</td>
<td align="center" valign="middle">1.19&#x2009;&#x00B1;&#x2009;0.01a</td>
<td align="center" valign="middle">63.70&#x2009;&#x00B1;&#x2009;0.19b</td>
<td align="center" valign="middle">35.20&#x2009;&#x00B1;&#x2009;0.41b</td>
<td align="center" valign="middle">7.50&#x2009;&#x00B1;&#x2009;0.09 b</td>
<td align="center" valign="middle">1.37&#x2009;&#x00B1;&#x2009;0.04a</td>
<td align="center" valign="middle">11.99&#x2009;&#x00B1;&#x2009;1.89a</td>
<td align="center" valign="middle">17.31&#x2009;&#x00B1;&#x2009;0.64a</td>
<td align="center" valign="middle">8.83&#x2009;&#x00B1;&#x2009;0.82a</td>
<td align="center" valign="middle">4.21&#x2009;&#x00B1;&#x2009;0.56a</td>
<td align="center" valign="middle">1.33&#x2009;&#x00B1;&#x2009;0.03a</td>
<td align="center" valign="middle">1.39&#x2009;&#x00B1;&#x2009;0.02b</td>
</tr>
<tr>
<td align="left" valign="middle">Plow&#x2013;till</td>
<td align="center" valign="middle">0.90&#x2009;&#x00B1;&#x2009;0.04b</td>
<td align="center" valign="middle">62.45&#x2009;&#x00B1;&#x2009;0.24b</td>
<td align="center" valign="middle">36.55&#x2009;&#x00B1;&#x2009;0.24a</td>
<td align="center" valign="middle">7.78&#x2009;&#x00B1;&#x2009;0.07 a</td>
<td align="center" valign="middle">1.18&#x2009;&#x00B1;&#x2009;0.06b</td>
<td align="center" valign="middle">9.63&#x2009;&#x00B1;&#x2009;0.53b</td>
<td align="center" valign="middle">15.37&#x2009;&#x00B1;&#x2009;0.62b</td>
<td align="center" valign="middle">7.33&#x2009;&#x00B1;&#x2009;0.28a</td>
<td align="center" valign="middle">2.03&#x2009;&#x00B1;&#x2009;0.52b</td>
<td align="center" valign="middle">1.29&#x2009;&#x00B1;&#x2009;0.02b</td>
<td align="center" valign="middle">1.54&#x2009;&#x00B1;&#x2009;0.11a</td>
</tr>
<tr>
<td align="left" valign="top" colspan="12">The silking stage</td>
</tr>
<tr>
<td align="left" valign="middle">Zero-till</td>
<td align="center" valign="middle">1.23&#x2009;&#x00B1;&#x2009;0.01a</td>
<td align="center" valign="middle">63.61&#x2009;&#x00B1;&#x2009;0.38a</td>
<td align="center" valign="middle">35.07&#x2009;&#x00B1;&#x2009;0.13b</td>
<td align="center" valign="middle">7.65&#x2009;&#x00B1;&#x2009;0.03 b</td>
<td align="center" valign="middle">1.34&#x2009;&#x00B1;&#x2009;0.12a</td>
<td align="center" valign="middle">13.26&#x2009;&#x00B1;&#x2009;2.50a</td>
<td align="center" valign="middle">14.16&#x2009;&#x00B1;&#x2009;0.68a</td>
<td align="center" valign="middle">0.32&#x2009;&#x00B1;&#x2009;0.04a</td>
<td align="center" valign="middle">2.25&#x2009;&#x00B1;&#x2009;0.30a</td>
<td align="center" valign="middle">1.38&#x2009;&#x00B1;&#x2009;0.01a</td>
<td align="center" valign="middle">1.39&#x2009;&#x00B1;&#x2009;0.05b</td>
</tr>
<tr>
<td align="left" valign="middle">Plow&#x2013;till</td>
<td align="center" valign="middle">0.79&#x2009;&#x00B1;&#x2009;0.03b</td>
<td align="center" valign="middle">62.55&#x2009;&#x00B1;&#x2009;0.22b</td>
<td align="center" valign="middle">36.71&#x2009;&#x00B1;&#x2009;0.23a</td>
<td align="center" valign="middle">7.92&#x2009;&#x00B1;&#x2009;0.05 a</td>
<td align="center" valign="middle">1.15&#x2009;&#x00B1;&#x2009;0.01b</td>
<td align="center" valign="middle">9.66&#x2009;&#x00B1;&#x2009;2.64b</td>
<td align="center" valign="middle">12.03&#x2009;&#x00B1;&#x2009;0.82b</td>
<td align="center" valign="middle">0.18&#x2009;&#x00B1;&#x2009;0.04b</td>
<td align="center" valign="middle">2.03&#x2009;&#x00B1;&#x2009;0.27a</td>
<td align="center" valign="middle">1.31&#x2009;&#x00B1;&#x2009;0.02b</td>
<td align="center" valign="middle">1.55&#x2009;&#x00B1;&#x2009;0.13a</td>
</tr>
<tr>
<td align="left" valign="top" colspan="12">Analysis of variance</td>
</tr>
<tr>
<td align="left" valign="top">Tillage (T)</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
</tr>
<tr>
<td align="left" valign="top">Time (Ti)</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;</td>
<td align="center" valign="top">ns</td>
</tr>
<tr>
<td align="left" valign="top">T&#x2009;&#x00D7;&#x2009;Ti</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">&#x002A;&#x002A;</td>
<td align="center" valign="top">ns</td>
<td align="center" valign="top">ns</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Results are reported as means&#x2009;&#x00B1;&#x2009;SD (<italic>n</italic>&#x2009;=&#x2009;3). Different letters within a row indicate significant differences between the treatments by Tukey&#x2019;s HSD test (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). &#x002A; and &#x002A;&#x002A; Indicates significant difference at <italic>p</italic>&#x2009;&#x003C;&#x2009;0.05, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.01, respectively; ns, indicates no significant difference. B10, soil 0&#x2009;~&#x2009;10&#x2009;cm bulk density; B10, soil 10&#x2009;~&#x2009;20&#x2009;cm bulk density. In addition to soil bulk density, the sampling depth of other soil properties was 0&#x2009;~&#x2009;20&#x2009;cm.</p>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption><p>The available concentrations of Fe, Mn, Cu, and Zn, as affected by the tillage practices.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Tillage</th>
<th align="center" valign="top" colspan="4">Soil microelement concentration (mg&#x2009;kg<sup>&#x2212;1</sup>)</th>
</tr>
<tr>
<th align="center" valign="top">Fe</th>
<th align="center" valign="top">Mn</th>
<th align="center" valign="top">Cu</th>
<th align="center" valign="top">Zn</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Zero-till</td>
<td align="center" valign="middle">5.79&#x2009;&#x00B1;&#x2009;0.11 a</td>
<td align="center" valign="middle">20.53&#x2009;&#x00B1;&#x2009;0.98 a</td>
<td align="center" valign="middle">1.03&#x2009;&#x00B1;&#x2009;0.04 a</td>
<td align="center" valign="middle">0.95&#x2009;&#x00B1;&#x2009;0.04 a</td>
</tr>
<tr>
<td align="left" valign="middle">Plow&#x2013;till</td>
<td align="center" valign="middle">5.39&#x2009;&#x00B1;&#x2009;0.15 b</td>
<td align="center" valign="middle">17.95&#x2009;&#x00B1;&#x2009;0.86 b</td>
<td align="center" valign="middle">0.94&#x2009;&#x00B1;&#x2009;0.05 b</td>
<td align="center" valign="middle">0.81&#x2009;&#x00B1;&#x2009;0.1 b</td>
</tr>
<tr>
<td align="left" valign="middle" colspan="5">ANOVA results</td>
</tr>
<tr>
<td align="left" valign="middle">Tillage (T)</td>
<td align="center" valign="middle">&#x002A;&#x002A;&#x002A;</td>
<td align="center" valign="middle">&#x002A;&#x002A;</td>
<td align="center" valign="middle">&#x002A;&#x002A;</td>
<td align="center" valign="middle">&#x002A;&#x002A;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Results are reported as means&#x2009;&#x00B1;&#x2009;SD (<italic>n</italic>&#x2009;=&#x2009;3). Different letters within a row indicate significant differences between the treatments by Tukey&#x2019;s HSD test (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). &#x002A;&#x002A; and &#x002A;&#x002A;&#x002A; Indicates significant difference at <italic>p</italic> &#x003C; 0.01, <italic>p</italic> &#x003C; 0.001, respectively; ns, indicates no significant difference.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec10">
<label>3.2.</label>
<title>Trends of antibiotic resistance gene families under long-term no tillage vs. plough tillage</title>
<p>Compared with PT, the application of long-term ZT significantly decreased the abundance of ARs and ATs (<xref rid="fig1" ref-type="fig">Figure 1</xref>, <italic>p</italic>&#x2009;&#x003C;&#x2009;0.05), along with ABSs and ASs. The NMDS further illustrated significant differences in ABSs (ANOSIM: <italic>r</italic>&#x2009;=&#x2009;0.56, <italic>p</italic>&#x2009;=&#x2009;0.031, <xref rid="fig2" ref-type="fig">Figure 2A</xref>), ARs (ANOSIM: <italic>r</italic>&#x2009;=&#x2009;0.89, <italic>p</italic>&#x2009;=&#x2009;0.002, <xref rid="fig2" ref-type="fig">Figure 2B</xref>), ASs (ANOSIM: <italic>r</italic>&#x2009;=&#x2009;0.24, <italic>p</italic>&#x2009;=&#x2009;0.049, <xref rid="fig2" ref-type="fig">Figure 2C</xref>), and ATs (ANOSIM: <italic>r</italic>&#x2009;=&#x2009;0.68, <italic>p</italic>&#x2009;=&#x2009;0.018, <xref rid="fig2" ref-type="fig">Figure 2D</xref>) families between the ZT and PT treatments.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption><p>The abundance (TPM: transcripts per kilobase per million mapped reads) of microbial antibiotic biosynthesis (ABSs, <bold>A</bold>), antibiotic resistance (ARs, <bold>B</bold>), antibiotic sensitive (ASs, <bold>C</bold>), and antibiotic target (ATs, <bold>D</bold>) following tillage practices. ZT, no tillage; PT, plough tillage.</p></caption>
<graphic xlink:href="fmicb-14-1238708-g001.tif"/>
</fig>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption><p>Nonmetric multidimensional scaling of the antibiotic biosynthesis (ABSs, <bold>A</bold>), antibiotic resistance (ARs, <bold>B</bold>), antibiotic sensitive (ASs, <bold>C</bold>), and antibiotic target (ATs, <bold>D</bold>) following tillage practices. ZT, no tillage; PT, plough tillage.</p></caption>
<graphic xlink:href="fmicb-14-1238708-g002.tif"/>
</fig>
<p>A total of 703,162 ARGs were identified from the 6,395,732 non-redundant gene catalogs of the whole metagenome (1.00%, average value), which were matched 13,104 ARO accession and assigned to 434 ARO names in the CARD database. Among these, 93.44, 0.05, and 0.56% were assigned to bacteria, fungi, and archaea, respectively, whereas the remaining 5.95% were unassigned. Moreover, most ABSs, ARs, ASs, and ATs ARGs were largely explained by the bacterial community, which accounted for 95.45% of the ABSs reads (Actinobacteria: 32.25%, Proteobacteria: 26.84%, Chloroflexi: 13.19%, and Gemmatimonadetes: 8.94%, etc), 93.43% of the ARs reads (Proteobacteria: 41.45%, Actinobacteria: 12.81%, and Acidobacteria: 12.34%, etc), 96.24% of the ASs reads (Proteobacteria: 41.49%, Actinobacteria: 15.90%, and Acidobacteria: 12.47%, and Gemmatimonadetes: 8.93%, etc), and 88.65% of the ATs reads (Proteobacteria: 36.16%, Acidobacteria: 14.72%, and Actinobacteria: 14.09%, etc) (<xref rid="fig1" ref-type="fig">Figure 1</xref>).</p>
<p>Among the functional groups of ARGs, metagenomic analysis revealed significant differences in the abundance of ARGs between ZT and PT soils (<xref rid="fig3" ref-type="fig">Figure 3</xref>). These trends in antibiotic biosynthesis, resistance, sensitivity, and target families were largely assigned to the dominant bacterial phyla, including Proteobacteria, Actinobacteria, Acidobacteria, and Gemmatimonadetes (<xref rid="fig4" ref-type="fig">Figure 4</xref>). In particular, Proteobacteria, Actinobacteria, Acdobacteria exhibited a higher number of TMP, which was higher in the PT than in the ZT soil, whereas Gemmatimonadetes and Chloroflexi were significantly lower in the PT than in the ZT soil.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption><p>Analysis for significant differences between the groups based on the abundance of antibiotic biosynthesis <bold>(A)</bold>, antibiotic resistance <bold>(B)</bold>, antibiotic sensitive <bold>(C)</bold>, and antibiotic target <bold>(D)</bold>. Gene abundances were based on TPM: transcripts per kilobase per million mapped reads of metagenomic reads. The ordinate represents the names of the genes that encode different enzymes, the abscisic represents the percentage value of a certain gene abundance of the sample, and different colors represent different groups. The right reflected by the abundance response ratios of different gene groups. Error bars represent 90% confidence intervals. &#x002A;Significant at <italic>p</italic>&#x2009;&#x003C;&#x2009;0.05; &#x002A;&#x002A;significant at <italic>p</italic>&#x2009;&#x003C;&#x2009;0.01.</p></caption>
<graphic xlink:href="fmicb-14-1238708-g003.tif"/>
</fig>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption><p>Contribution of microbial (bacterial and fungal) phyla to microbial ARGs genes for the Antibiotic Biosynthesis, Antibiotic Resistance, Antibiotic Sensitive, and Antibiotic Target classes. <bold>(A)</bold> Contribution of microbial (bacterial and fungal) phyla to Antibiotic Biosynthesis. <bold>(B)</bold> Contribution of microbial (bacterial and fungal) phyla to Antibiotic Resistance. <bold>(C)</bold> Contribution of microbial (bacterial and fungal) phyla to Antibiotic Sensitive. <bold>(D)</bold> Contribution of microbial (bacterial and fungal) phyla to Antibiotic Target.</p></caption>
<graphic xlink:href="fmicb-14-1238708-g004.tif"/>
</fig>
</sec>
<sec id="sec11">
<label>3.3.</label>
<title>Significant differences of antibiotic resistance gene profile under long-term no tillage vs. plough tillage</title>
<p>For the antibiotic resistance gene family, significantly different gene abundances, such as <italic>oqxB</italic>, <italic>kdpE</italic>, <italic>mexK</italic>, <italic>dfrE</italic>, <italic>tlyA</italic>, etc., were significantly inhibited in ZT soil (<xref rid="fig3" ref-type="fig">Figure 3A</xref>). In addition, antibiotic-sensitive and antibiotic target families were enriched in the PT treatment (<xref rid="fig3" ref-type="fig">Figures 3C</xref>,<xref rid="fig3" ref-type="fig">D</xref>). There were no significant differences in the abundance of antibiotic biosynthesis genes between the two tillage treatments (<xref rid="fig3" ref-type="fig">Figure 3B</xref>). These results indicate that long-term ZT can effectively reduce the relative abundance of ARGs.</p>
<p>Furthermore, after analyzing the antibiotics encoded by the significant difference antibiotic resistance genes family, we found that the abundance of <italic>fluoroquinolone</italic>, <italic>efflux pump</italic>, <italic>aminoglycosides</italic>, <italic>isoniazid</italic>, <italic>trimethoprim</italic>, <italic>multi-drug resistance</italic>, <italic>polymyxins</italic>, <italic>fosfomycin</italic>, <italic>MLSB</italic>, <italic>&#x03B2;-lactams</italic>, <italic>chloramphenicol</italic>, <italic>sulfonamides</italic>, <italic>ethambutol</italic> in the two soils was significantly different (<xref rid="tab3" ref-type="table">Table 3</xref>). Compared with PT soil, the abundance of these 13 antibiotic genes was significantly decreased by 6.21, 23.46, 16.79, 11.34, 13.21, 26.22, 33.67, 35.82, 32.78, 44.87, 28.57, 34.64, and 59.97% in ZT soil, respectively, indicating that long-term ZT application has a positive effect on reducing the cumulative risk of antibiotics in agricultural soil.</p>
<table-wrap position="float" id="tab3">
<label>Table 3</label>
<caption><p>The abundance of drug resistance, as affected by the tillage practices.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="middle">Drug resistance type</th>
<th align="center" valign="middle">Zero-till</th>
<th align="center" valign="middle">Plow&#x2013;till</th>
<th align="center" valign="middle"><italic>F</italic></th>
<th align="center" valign="middle"><italic>P</italic></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Fluoroquinolone</td>
<td align="center" valign="middle">6935.33&#x2009;&#x00B1;&#x2009;168.40 b</td>
<td align="center" valign="middle">7394.67&#x2009;&#x00B1;&#x2009;273.18 a</td>
<td align="center" valign="middle">8.441</td>
<td align="center" valign="middle">0.044</td>
</tr>
<tr>
<td align="left" valign="middle">Efflux pump</td>
<td align="center" valign="middle">2259.33&#x2009;&#x00B1;&#x2009;76.38 b</td>
<td align="center" valign="middle">2952.00&#x2009;&#x00B1;&#x2009;203.45 a</td>
<td align="center" valign="middle">30.497</td>
<td align="center" valign="middle">0.005</td>
</tr>
<tr>
<td align="left" valign="middle">Aminoglycosides</td>
<td align="center" valign="middle">1189.33&#x2009;&#x00B1;&#x2009;76.27 b</td>
<td align="center" valign="middle">1429.33&#x2009;&#x00B1;&#x2009;51.63 a</td>
<td align="center" valign="middle">20.371</td>
<td align="center" valign="middle">0.011</td>
</tr>
<tr>
<td align="left" valign="middle">Isoniazid</td>
<td align="center" valign="middle">1022.00&#x2009;&#x00B1;&#x2009;80.07 b</td>
<td align="center" valign="middle">1152.67&#x2009;&#x00B1;&#x2009;46.88 a</td>
<td align="center" valign="middle">18.363</td>
<td align="center" valign="middle">0.013</td>
</tr>
<tr>
<td align="left" valign="middle">Trimethoprim</td>
<td align="center" valign="middle">910.67&#x2009;&#x00B1;&#x2009;15.01 b</td>
<td align="center" valign="middle">1049.33&#x2009;&#x00B1;&#x2009;91.22 a</td>
<td align="center" valign="middle">32.652</td>
<td align="center" valign="middle">0.005</td>
</tr>
<tr>
<td align="left" valign="middle">Multi-drug resistance</td>
<td align="center" valign="middle">435.33&#x2009;&#x00B1;&#x2009;22.30 b</td>
<td align="center" valign="middle">590.00&#x2009;&#x00B1;&#x2009;57.58 a</td>
<td align="center" valign="middle">18.820</td>
<td align="center" valign="middle">0.012</td>
</tr>
<tr>
<td align="left" valign="middle">Polymyxins</td>
<td align="center" valign="middle">174.67&#x2009;&#x00B1;&#x2009;25.79 b</td>
<td align="center" valign="middle">263.33&#x2009;&#x00B1;&#x2009;19.73 a</td>
<td align="center" valign="middle">22.363</td>
<td align="center" valign="middle">0.009</td>
</tr>
<tr>
<td align="left" valign="middle">Fosfomycin</td>
<td align="center" valign="middle">57.33&#x2009;&#x00B1;&#x2009;18.15 b</td>
<td align="center" valign="middle">89.33&#x2009;&#x00B1;&#x2009;9.45 a</td>
<td align="center" valign="middle">31.135</td>
<td align="center" valign="middle">0.005</td>
</tr>
<tr>
<td align="left" valign="middle">MLSB</td>
<td align="center" valign="middle">54.67&#x2009;&#x00B1;&#x2009;15.01 b</td>
<td align="center" valign="middle">81.33&#x2009;&#x00B1;&#x2009;9.02 a</td>
<td align="center" valign="middle">20.000</td>
<td align="center" valign="middle">0.011</td>
</tr>
<tr>
<td align="left" valign="middle">&#x03B2;-lactams</td>
<td align="center" valign="middle">28.67&#x2009;&#x00B1;&#x2009;2.31 b</td>
<td align="center" valign="middle">52.00&#x2009;&#x00B1;&#x2009;8.00 a</td>
<td align="center" valign="middle">23.558</td>
<td align="center" valign="middle">0.008</td>
</tr>
<tr>
<td align="left" valign="middle">Chloramphenicol</td>
<td align="center" valign="middle">30.00&#x2009;&#x00B1;&#x2009;5.29 b</td>
<td align="center" valign="middle">42.00&#x2009;&#x00B1;&#x2009;4.00 a</td>
<td align="center" valign="middle">9.818</td>
<td align="center" valign="middle">0.035</td>
</tr>
<tr>
<td align="left" valign="middle">Sulfonamides</td>
<td align="center" valign="middle">2.00&#x2009;&#x00B1;&#x2009;0.00 b</td>
<td align="center" valign="middle">3.06&#x2009;&#x00B1;&#x2009;3.06 a</td>
<td align="center" valign="middle">9.143</td>
<td align="center" valign="middle">0.039</td>
</tr>
<tr>
<td align="left" valign="middle">Ethambutol</td>
<td align="center" valign="middle">2.67&#x2009;&#x00B1;&#x2009;1.15 b</td>
<td align="center" valign="middle">6.67&#x2009;&#x00B1;&#x2009;1.15 a</td>
<td align="center" valign="middle">18.000</td>
<td align="center" valign="middle">0.013</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Results are reported as means&#x2009;&#x00B1;&#x2009;SD (<italic>n</italic> =&#x2009;3). The unit of drug resistance abundance is transcripts per kilobase per million mapped reads (TPM). Different letters within a row indicate significant differences between the treatments by Tukey&#x2019;s HSD test (<italic>p</italic> &#x003C;&#x2009;0.05).</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec12">
<label>3.4.</label>
<title>Changes in soil microbial metabolic coefficients following long-term no tillage vs. plough tillage</title>
<p>We found that the long-term application of no-tillage practices can significantly change soil microbial diversity and metabolic coefficients (<xref rid="tab4" ref-type="table">Table 4</xref>). Compared with PT, the soil organic content, microbial diversity, and microbial biomass carbon significantly increased after continuous application (11&#x2009;yr) of ZT (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). Moreover, microbial respiration and qCO<sub>2</sub> decreased significantly with increasing age under no-tillage application (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). Whereas, the Cmic:Corg ratio increased significantly with ZT application (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05).</p>
<table-wrap position="float" id="tab4">
<label>Table 4</label>
<caption><p>Changes in soil microbial metabolic coefficients following long-term no tillage.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Parameters</th>
<th align="center" valign="top">Zero-till</th>
<th align="center" valign="top">Plow&#x2013;till</th>
<th align="center" valign="top"><italic>F</italic></th>
<th align="center" valign="top"><italic>P</italic></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">SOC (g/kg)</td>
<td align="center" valign="middle">13.98&#x2009;&#x00B1;&#x2009;0.56 a</td>
<td align="center" valign="middle">10.83&#x2009;&#x00B1;&#x2009;0.37 b</td>
<td align="center" valign="middle">59.339</td>
<td align="center" valign="middle">&#x003C;0.001</td>
</tr>
<tr>
<td align="left" valign="middle">MBC (mg/kg)</td>
<td align="center" valign="middle">513.18&#x2009;&#x00B1;&#x2009;22.5 a</td>
<td align="center" valign="middle">353.03&#x2009;&#x00B1;&#x2009;8.51 b</td>
<td align="center" valign="middle">177.291</td>
<td align="center" valign="middle">&#x003C;0.001</td>
</tr>
<tr>
<td align="left" valign="middle">Shannon index</td>
<td align="center" valign="middle">6.85&#x2009;&#x00B1;&#x2009;0.03 a</td>
<td align="center" valign="middle">6.74&#x2009;&#x00B1;&#x2009;0.03 b</td>
<td align="center" valign="middle">15.876</td>
<td align="center" valign="middle">0.011</td>
</tr>
<tr>
<td align="left" valign="middle">MR (mg/kg)</td>
<td align="center" valign="middle">113.14&#x2009;&#x00B1;&#x2009;1.51 b</td>
<td align="center" valign="middle">149.45&#x2009;&#x00B1;&#x2009;3.37 a</td>
<td align="center" valign="middle">86.450</td>
<td align="center" valign="middle">&#x003C;0.001</td>
</tr>
<tr>
<td align="left" valign="middle">qCO<sub>2</sub></td>
<td align="center" valign="middle">0.22&#x2009;&#x00B1;&#x2009;0.01 b</td>
<td align="center" valign="middle">0.42&#x2009;&#x00B1;&#x2009;0.02 a</td>
<td align="center" valign="middle">458.68</td>
<td align="center" valign="middle">&#x003C;0.001</td>
</tr>
<tr>
<td align="left" valign="middle">C<sub>mic</sub>: C<sub>org</sub></td>
<td align="center" valign="middle">36.71&#x2009;&#x00B1;&#x2009;1.28 a</td>
<td align="center" valign="middle">32.61&#x2009;&#x00B1;&#x2009;0.93 b</td>
<td align="center" valign="middle">26.859</td>
<td align="center" valign="middle">0.002</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Values are represented as the mean value followed by standard error in parentheses (<italic>n</italic>&#x2009;=&#x2009;3). Different lowercase letters indicate a significant difference (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05) between different tillage practices, based on a one-way ANOVA followed by Tukey&#x2019;s HSD test (<italic>p</italic>&#x2009;&#x003C;&#x2009;0.05). SOC: soil organic carbon, MBC: microbial biomass carbon; Shannon index: Microbial diversity; MR: microbial respiration; qCO<sub>2</sub> stands for the metabolic quotient, which was estimated by MR: MBC; C<sub>mic</sub>: C<sub>org</sub> stands for the microbial quotient, which was estimated by MBC: SOC.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec13">
<label>3.5.</label>
<title>Specific ARGs family responses for microbial metabolic activity and soil properties</title>
<p>The Mantel test indicated that soil texture (clay, silt, and sand), soil moisture, pH, TN, MBC, SOC, and microbial metabolic activity (diversity, MR, qCO<sub>2</sub>, and Cmic:Corg ratio) correlated significantly with various genes of the antibiotic biosynthesis family (<xref rid="fig5" ref-type="fig">Figure 5</xref>), in which the effect values of MBC and SOC were the largest. The significantly different genes of the antibiotic resistance family exhibited a strong association with soil texture (i.e., clay, silt, and sand), soil temperature, soil moisture, pH, TN, available Zn, MBC, SOC, and microbial metabolic activity (i.e., diversity, MR, qCO<sub>2</sub>, and Cmic:Corg ratio), in which pH, diversity, MR, and qCO<sub>2</sub> were the main regulatory factors. Soil temperature, Fe, and MR had a significant impact on the differentially expressed genes of the antibiotic-sensitive family, whereas soil temperature and MR were significantly correlated with differentially expressed genes of the antibiotic target family. However, there was no significant correlation between the antibiotic target gene family and the soil properties.</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption><p>Relationship among the abundance of the ARGs genes involved in antibiotic biosynthesis (ABSs), antibiotic resistance (ARs), antibiotic sensitive (ASs), and antibiotic target (ATs), microbial metabolic activity and soil properties following tillage practices.</p></caption>
<graphic xlink:href="fmicb-14-1238708-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="sec14" sec-type="discussions">
<label>4.</label>
<title>Discussion</title>
<p>Agricultural management practices, such as fertilization types, agricultural planting patterns, and planting crops, are assumed to introduce or enrich a large number of ARGs and exert selection pressure on indigenous soil microbial communities (<xref ref-type="bibr" rid="ref11">Forsberg et al., 2014</xref>; <xref ref-type="bibr" rid="ref45">Wang et al., 2017</xref>; <xref ref-type="bibr" rid="ref21">Li et al., 2022</xref>). However, there is a dearth of knowledge regarding how long-term conservation tillage influences the patterns of a wide spectrum of ARGs in dry agricultural soils. In this study, non-metric multidimensional scaling (NMDS) revealed that the soil ARGs clustered into two groups based on the two tillage treatments (<xref rid="fig2" ref-type="fig">Figure 2</xref>). These findings indicate that the community structure and abundance of ARGs in the soil were influenced by long-term ZT and PT applications. The total abundance of microbial antibiotic resistance, sensitive, and target genes was observed to be lower after long-term ZT application than after PT (<xref rid="fig1" ref-type="fig">Figure 1</xref>). This suggests that the abundance of soil ARGs may decrease with reduced soil disturbance intensity.</p>
<p>Among the general mechanisms of antibiotic resistance, multidrug efflux pumps excluding drugs serve as critical one in multidrug resistance (<xref ref-type="bibr" rid="ref16">Hernando-Amado et al., 2016</xref>). Our findings indicate that long-term ZT significantly reduced the abundance of soil efflux pump conferring ARGs, with a significant difference of 23.46% compared to PT soil. Simultaneously, a large number and wide distribution of broad-spectrum efflux pumps, such as <italic>oqxB</italic>, <italic>mexK</italic>, <italic>efpA</italic>, <italic>adeH</italic>, <italic>acrD</italic>, <italic>lrfA vgaALC</italic>, <italic>amrB</italic>, <italic>mexT</italic>, <italic>acrE</italic>, <italic>lmrB</italic>, and <italic>mdtL</italic> in no-tillage soils than in plough tillage soils (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>). These results suggest that long-term ZT can reduce the accumulation of antibiotic-resistant microorganisms by weakening the generalized role of broad-spectrum efflux pumps in protecting cells from multiple antibiotics encountered in the agricultural environment (<xref ref-type="bibr" rid="ref54">Zhang et al., 2021</xref>). In addition, compared to PT soil, the abundance of <italic>fluoroquinolone, aminoglycosides, isoniazid, trimethoprim, multi-drug resistance, polymyxins, fosfomycin, MLSB, &#x03B2;-lactams, chloramphenicol, sulfonamides, and ethambutol</italic> resistance genes under long ZT soil were also significantly decreased. This may be because long-term ZT creates a soil environment suitable for microbial proliferation, which significantly increases the abundance of beneficial microorganisms, thereby reducing the abundance of antibiotic-resistant microbial hosts (<xref ref-type="bibr" rid="ref10">Dunivin and Shade, 2018</xref>). These results indicate that long-term ZT application for 11 consecutive years can inhibit the extensive enrichment of ARG abundance in agricultural soil and reduce the spread of resistance (<xref ref-type="bibr" rid="ref14">Gullberg et al., 2014</xref>; <xref ref-type="bibr" rid="ref2">Blanco et al., 2016</xref>). This can greatly contribute to the control of soil resistance gene pollution and achieve sustainable agricultural production.</p>
<p>We used mantel test to identify the key factors affecting the accumulation of ARGs, including physical and chemical soil properties, microbial diversity, and microbial metabolic activity (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S2</xref>). The conservation tillage has important effects on the accumulation of ARGs in soil by regulating soil microhabitats. Specifically, the abundance of antibiotic biosynthesis, resistance, sensitivity, and target genes was significantly regulated by soil pH (<xref rid="fig5" ref-type="fig">Figure 5</xref>). This is consistent with previous studies showing that pH is the most important factor affecting the microbial community and accumulation of antibiotic resistance genes (<xref ref-type="bibr" rid="ref36">Rousk et al., 2010</xref>; <xref ref-type="bibr" rid="ref20">Li et al., 2020</xref>). Changes in soil pH led to significant changes in the soil physicochemical properties and microbial communities, resulting in changes in the accumulation of ARGs in the soil (<xref ref-type="bibr" rid="ref23">Liu et al., 2017</xref>). In addition, mental analysis showed that the SOC content had a significant effect on AGR accumulation (<xref rid="fig5" ref-type="fig">Figure 5</xref>). The ARGs in the no-tillage soil were significantly lower than those in the plough tillage soil, mainly because long-term no-tillage formed a carbon-rich microenvironment that was conducive to the balanced proliferation of various microbial populations, thereby inhibiting the niche of the ARGs host population and reducing the accumulation of ARGs in the soil (<xref ref-type="bibr" rid="ref11">Forsberg et al., 2014</xref>). Besides, our result also found that the concentration of soil available Zn had a significant effect on the Ars family genes, and the concentration of available Fe had a significant effect on the ASs family genes (<xref rid="fig5" ref-type="fig">Figure 5</xref>). This was mainly due to long-term no-tillage significantly changed the concentration of available Zn and Fe, which had an effect on soil bacteria and bacteriophages, and thus inhibited the accumulation of ARGs in no-tillage soil (<xref ref-type="bibr" rid="ref8">Dickinson et al., 2019</xref>).</p>
<p>In addition to the influence of soil properties on ARGs, the diversity and structural composition of the microbial community also significantly affect ARGs (<xref ref-type="bibr" rid="ref10">Dunivin and Shade, 2018</xref>). Our previous study found that long-term ZT significantly increased soil bacterial and fungal alpha diversity (<xref ref-type="bibr" rid="ref44">Wang et al., 2021</xref>), which may be another key biological factor in reducing ARGs accumulation in the soil. A higher soil microbial diversity can significantly inhibit the spread of resistance genes (<xref ref-type="bibr" rid="ref7">Chen et al., 2019</xref>). Previous studies have confirmed that the decrease of &#x03B1;-diversity of soil microbial communities can lead to the increase the abundance of <italic>MLSB</italic>, <italic>aminoglycosides</italic>, <italic>&#x03B2;-lactams</italic>, <italic>amide alcohols</italic>, <italic>sulfonamides</italic> and <italic>tnpA</italic> resistance genes (<xref ref-type="bibr" rid="ref15">Han et al., 2018</xref>; <xref ref-type="bibr" rid="ref42">Urra et al., 2019</xref>). These results are consistent with those of this study (<xref rid="tab3" ref-type="table">Table 3</xref>), which may be because long-term ZT increased the content of organic matter and mineral nutrients in soil, thus increasing the niche breadth of microorganisms in the ZT soil (<xref ref-type="bibr" rid="ref11">Forsberg et al., 2014</xref>), alleviating the selective pressure on microorganisms (<xref ref-type="bibr" rid="ref5">Chase and Myers, 2011</xref>), and compressing the living space of antibiotic-carrying host microorganisms (<xref ref-type="bibr" rid="ref10">Dunivin and Shade, 2018</xref>). Furthermore, <xref ref-type="bibr" rid="ref11">Forsberg et al. (2014)</xref> and <xref ref-type="bibr" rid="ref22">Li B. et al., (2015)</xref> found robust correlations between bacterial phyla and ARGs, suggesting that the microbial structure composition could be one of the factors affecting the soil resistome. Our previous studies have also shown that long-term ZT application has significant effects on soil bacterial and fungal community structures (<xref ref-type="bibr" rid="ref45">Wang et al., 2017</xref>, <xref ref-type="bibr" rid="ref44">2021</xref>), suggesting that ZT may lead to changes in the abundance and composition of host populations carrying antibiotic resistance genes. For examples, compared with PT, the relative abundance of <italic>Firmicutes</italic> in ZT soils increased significantly, whereas the relative abundance of <italic>Alphaproteobacteria</italic>, <italic>Betaproteobacteria</italic> and <italic>&#x03B3;-proteobacteria</italic> decreased significantly (<xref rid="fig4" ref-type="fig">Figure 4</xref>). The average number of ARGs in the genome of <italic>Proteobacteria</italic> is the largest, and their proliferation can promote the growth and creation of drug-resistant bacteria in the soil (<xref ref-type="bibr" rid="ref57">Zhou et al., 2015</xref>; <xref ref-type="bibr" rid="ref52">Xun et al., 2016</xref>). <xref ref-type="bibr" rid="ref41">Udikovic-Kolic et al. (2014)</xref> found that the alteration of the soil resistome can mainly be attributed to an increase in certain bacterial populations, such as <italic>Pseudomonas</italic> spp., <italic>Janthinobacterium</italic> sp. and <italic>Psychrobacter pulmonis</italic>, which harbor ARGs in farmland soil.</p>
<p>Microbial metabolic activity also had a significant effect on the accumulation of antibiotic resistance genes (<xref rid="fig5" ref-type="fig">Figure 5</xref>). The results showed that the long-term ZT soil microbial metabolic activity was significantly higher than that of the PT soil, which may be another key factor in reducing the accumulation of ARGs in long-term ZT soils. This is mainly because the higher the metabolic activity of microorganisms, the higher their versatility and health, thereby occupying a niche and crowding out the living space of antibiotic-carrying host microorganisms (<xref ref-type="bibr" rid="ref53">Yang et al., 2022</xref>). <xref ref-type="bibr" rid="ref40">Su et al. (2015)</xref> reported similar findings, demonstrating that soil with stronger microbial metabolic activity and a higher relative balance of community structure composition can inhibit the abundance and activity of pathogenic microorganisms and antibiotic gene hosts more effectively.</p>
</sec>
<sec id="sec15" sec-type="conclusions">
<label>5.</label>
<title>Conclusion</title>
<p>The long-term application of ZT significantly reduced the accumulation of fluoroquinolone, aminoglycosides, isoniazid, trimethoprim, multi-drug resistance, polymyxins, fosfomycin, MLSB, &#x03B2;-lactams, chloramphenicol, sulfonamides, and ethambutol resistance genes under dryland conditions, compared to PT soil. Such antibiotic resistance gene reduction is closely associated with soil properties, microbial community composition, and metabolic activity. Our findings suggest that long-term ZT practices can reduce the accumulation of antibiotic resistance genes through the following mechanisms: decreased pH and increased SOC content are key physicochemical factors that inhibit ARG accumulation in no-tillage soil, while decreased abundance of Proteobacteria, Actinobacteria, and Acidobacteria communities, and increased abundance of Gemmatimonadetes and Chloroflexi communities, as well as microbial metabolic activity, are key biological factors that inhibit the reproduction of antibiotic-resistant host microbes. This implies that long-term no-tillage practices inhibit the accumulation of antibiotic resistance genes in farmland soil, which is a promising agricultural management measure to reduce the accumulation risk of soil ARGs. Therefore, combining conservation tillage with other high-yield and efficient agricultural practices will greatly improve the sustainability of agricultural production and make a significant contribution to human health in future agricultural production. Future research should identify the risks of active antibiotic resistance gene host microbial populations in the agricultural environment at the transcriptome level, based on RNA sequencing, to ensure that ZT is used sustainably.</p>
</sec>
<sec id="sec16" sec-type="data-availability">
<title>Data availability statement</title>
<p>Raw reads were uploaded to the NCBI Sequence Read Archive database (accession no. PRJNA641749).</p>
</sec>
<sec id="sec17">
<title>Author contributions</title>
<p>WW: conceptualization, investigation, data curation, writing&#x2014;original draft, review, editing, and funding acquisition. PS and ZL: partial data investigation. FM: conceptualization, writing&#x2014;reviewing and editing. YL: conceptualization and funding acquisition. XW: conceptualization, funding acquisition, supervision, writing&#x2014;reviewing and editing. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec sec-type="funding-information" id="sec18">
<title>Funding</title>
<p>This study was financially supported by the Science and Technology Plan Project of Shaanxi Province (Key Industry Innovation Chain Project) (grant number 2023-ZDLNY-03), Starting Research Fund from the Northwest A&#x0026;F University (grant number 2452021112), and National Natural Science Foundation of China (grant number 31971860).</p>
</sec>
<sec sec-type="COI-statement" id="sec19">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="sec143" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="sec20" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2023.1238708/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2023.1238708/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Achiba</surname> <given-names>W. B.</given-names></name> <name><surname>Lakhdar</surname> <given-names>A.</given-names></name> <name><surname>Gabteni</surname> <given-names>N.</given-names></name> <name><surname>Laing</surname> <given-names>G. D.</given-names></name> <name><surname>Verloo</surname> <given-names>M.</given-names></name> <name><surname>Boeckx</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>Accumulation and fractionation of trace metals in a Tunisian calcareous soil amended with farmyard manure and municipal solid waste compost</article-title>. <source>J. Hazard. Mater.</source> <volume>176</volume>, <fpage>99</fpage>&#x2013;<lpage>108</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jhazmat.2009.11.004</pub-id>, PMID: <pub-id pub-id-type="pmid">19945790</pub-id></citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blanco</surname> <given-names>P.</given-names></name> <name><surname>Hernando-Amado</surname> <given-names>S.</given-names></name> <name><surname>Reales-Calderon</surname> <given-names>J. A.</given-names></name> <name><surname>Corona</surname> <given-names>F.</given-names></name> <name><surname>Lira</surname> <given-names>F.</given-names></name> <name><surname>Alcalde-Rico</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Bacterial multidrug efflux pumps: much more than antibiotic resistance determinants</article-title>. <source>Microorganisms</source> <volume>4</volume>, <fpage>2076</fpage>&#x2013;<lpage>2607</lpage>. doi: <pub-id pub-id-type="doi">10.3390/microorganisms4010014</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blau</surname> <given-names>K.</given-names></name> <name><surname>Bettermann</surname> <given-names>A.</given-names></name> <name><surname>Jechalke</surname> <given-names>S.</given-names></name> <name><surname>Fornefeld</surname> <given-names>E.</given-names></name> <name><surname>Vanrobaeys</surname> <given-names>Y.</given-names></name> <name><surname>Stalder</surname> <given-names>T.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>The transferable resistome of produce</article-title>. <source>MBio</source> <volume>9</volume>, <fpage>2150</fpage>&#x2013;<lpage>7511</lpage>. doi: <pub-id pub-id-type="doi">10.1128/mBio.01300-18</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Buchfink</surname> <given-names>B.</given-names></name> <name><surname>Reuter</surname> <given-names>K.</given-names></name> <name><surname>Drost</surname> <given-names>H. G.</given-names></name></person-group> (<year>2021</year>). <article-title>Sensitive protein alignments at tree-of-life scale using DIAMOND</article-title>. <source>Nat. Methods</source> <volume>18</volume>, <fpage>366</fpage>&#x2013;<lpage>368</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41592-021-01101-x</pub-id>, PMID: <pub-id pub-id-type="pmid">33828273</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chase</surname> <given-names>J. M.</given-names></name> <name><surname>Myers</surname> <given-names>J. A.</given-names></name></person-group> (<year>2011</year>). <article-title>Disentangling the importance of ecological niches from stochastic processes across scales</article-title>. <source>Philos. T. R. Soc. B.</source> <volume>366</volume>, <fpage>2351</fpage>&#x2013;<lpage>2363</lpage>. doi: <pub-id pub-id-type="doi">10.1098/rstb.2011.0063</pub-id>, PMID: <pub-id pub-id-type="pmid">21768151</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>Q.</given-names></name> <name><surname>An</surname> <given-names>X.</given-names></name> <name><surname>Li</surname> <given-names>H.</given-names></name> <name><surname>Su</surname> <given-names>J.</given-names></name> <name><surname>Ma</surname> <given-names>Y.</given-names></name> <name><surname>Zhu</surname> <given-names>Y. G.</given-names></name></person-group> (<year>2016</year>). <article-title>Long-term field application of sewage sludge increases the abundance of antibiotic resistance genes in soil</article-title>. <source>Environ. Int.</source> <volume>92-93</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envint.2016.03.026</pub-id>, PMID: <pub-id pub-id-type="pmid">27043971</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>Q. L.</given-names></name> <name><surname>An</surname> <given-names>X. L.</given-names></name> <name><surname>Zheng</surname> <given-names>B. X.</given-names></name> <name><surname>Gillings</surname> <given-names>M.</given-names></name> <name><surname>Pe&#x00F1;uelas</surname> <given-names>J.</given-names></name> <name><surname>Cui</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Loss of soil microbial diversity exacerbates spread of antibiotic resistance</article-title>. <source>Soil Ecol. Lett.</source> <volume>1</volume>, <fpage>3</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s42832-019-0011-0</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dickinson</surname> <given-names>A. W.</given-names></name> <name><surname>Power</surname> <given-names>A.</given-names></name> <name><surname>Hansen</surname> <given-names>M. G.</given-names></name> <name><surname>Brandt</surname> <given-names>K. K.</given-names></name> <name><surname>Piliposian</surname> <given-names>G.</given-names></name> <name><surname>Appleby</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Heavy metal pollution and co-selection for antibiotic resistance: a microbial palaeontology approach</article-title>. <source>Environ. Int.</source> <volume>132</volume>:<fpage>105117</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envint.2019.105117</pub-id>, PMID: <pub-id pub-id-type="pmid">31473413</pub-id></citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dumbrell</surname> <given-names>A. J.</given-names></name> <name><surname>Nelson</surname> <given-names>M.</given-names></name> <name><surname>Helgason</surname> <given-names>T.</given-names></name> <name><surname>Dytham</surname> <given-names>C.</given-names></name> <name><surname>Fitter</surname> <given-names>A. H.</given-names></name></person-group> (<year>2010</year>). <article-title>Relative roles of niche and neutral processes in structuring a soil microbial community</article-title>. <source>ISME J.</source> <volume>4</volume>, <fpage>337</fpage>&#x2013;<lpage>345</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2009.122</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dunivin</surname> <given-names>T. K.</given-names></name> <name><surname>Shade</surname> <given-names>A.</given-names></name></person-group> (<year>2018</year>). <article-title>Community structure explains antibiotic resistance gene dynamics over a temperature gradient in soil</article-title>. <source>FEMS Microbiol. Ecol.</source> <volume>94</volume>:<fpage>fiy016</fpage>. doi: <pub-id pub-id-type="doi">10.1093/femsec/fiy016</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Forsberg</surname> <given-names>K. J.</given-names></name> <name><surname>Patel</surname> <given-names>S.</given-names></name> <name><surname>Gibson</surname> <given-names>M. K.</given-names></name> <name><surname>Lauber</surname> <given-names>C. L.</given-names></name> <name><surname>Knight</surname> <given-names>R.</given-names></name> <name><surname>Fierer</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Bacterial phylogeny structures soil resistomes across habitats</article-title>. <source>Nature</source> <volume>509</volume>, <fpage>612</fpage>&#x2013;<lpage>616</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nature13377</pub-id>, PMID: <pub-id pub-id-type="pmid">24847883</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Forsberg</surname> <given-names>K. J.</given-names></name> <name><surname>Reyes</surname> <given-names>A.</given-names></name> <name><surname>Bin</surname> <given-names>W.</given-names></name> <name><surname>Selleck</surname> <given-names>E. M.</given-names></name> <name><surname>Sommer</surname> <given-names>M. O. A.</given-names></name> <name><surname>Dantas</surname> <given-names>G.</given-names></name></person-group> (<year>2012</year>). <article-title>The shared antibiotic resistome of soil bacteria and human pathogens</article-title>. <source>Science</source> <volume>337</volume>, <fpage>1095</fpage>&#x2013;<lpage>9203</lpage>. doi: <pub-id pub-id-type="doi">10.1126/science.1220761</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fu</surname> <given-names>L. M.</given-names></name> <name><surname>Niu</surname> <given-names>B. F.</given-names></name> <name><surname>Zhu</surname> <given-names>Z. W.</given-names></name> <name><surname>Wu</surname> <given-names>S. T.</given-names></name> <name><surname>Li</surname> <given-names>W. Z.</given-names></name></person-group> (<year>2012</year>). <article-title>CD-HIT: accelerated for clustering the next-generation sequencing data</article-title>. <source>Bioinformatics</source> <volume>28</volume>, <fpage>3150</fpage>&#x2013;<lpage>3152</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/bts565</pub-id>, PMID: <pub-id pub-id-type="pmid">23060610</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gullberg</surname> <given-names>E.</given-names></name> <name><surname>Albrecht Lisa</surname> <given-names>M.</given-names></name> <name><surname>Karlsson</surname> <given-names>C.</given-names></name> <name><surname>Sandegren</surname> <given-names>L.</given-names></name> <name><surname>Andersson Dan</surname> <given-names>I.</given-names></name></person-group> (<year>2014</year>). <article-title>Selection of a multidrug resistance plasmid by sublethal levels of antibiotics and heavy metals</article-title>. <source>mBio</source> <volume>5</volume>:<fpage>e01918-14</fpage>. doi: <pub-id pub-id-type="doi">10.1128/mBio.01918-14</pub-id>, PMID: <pub-id pub-id-type="pmid">25293762</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Han</surname> <given-names>X. M.</given-names></name> <name><surname>Hu</surname> <given-names>H. W.</given-names></name> <name><surname>Chen</surname> <given-names>Q. L.</given-names></name> <name><surname>Yang</surname> <given-names>L. Y.</given-names></name> <name><surname>Li</surname> <given-names>H. L.</given-names></name> <name><surname>Zhu</surname> <given-names>Y. G.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Antibiotic resistance genes and associated bacterial communities in agricultural soils amended with different sources of animal manures</article-title>. <source>Soil Biol. Biochem.</source> <volume>126</volume>, <fpage>91</fpage>&#x2013;<lpage>102</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2018.08.018</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hernando-Amado</surname> <given-names>S.</given-names></name> <name><surname>Blanco</surname> <given-names>P.</given-names></name> <name><surname>Alcalde-Rico</surname> <given-names>M.</given-names></name> <name><surname>Corona</surname> <given-names>F.</given-names></name> <name><surname>Reales-Calderon</surname> <given-names>J. A.</given-names></name> <name><surname>Sanchez</surname> <given-names>M. B.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Multidrug efflux pumps as main players in intrinsic and acquired resistance to antimicrobials</article-title>. <source>Drug Resist. Update</source> <volume>28</volume>, <fpage>13</fpage>&#x2013;<lpage>27</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.drup.2016.06.007</pub-id>, PMID: <pub-id pub-id-type="pmid">27620952</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huijbers</surname> <given-names>P. M. C.</given-names></name> <name><surname>Blaak</surname> <given-names>H.</given-names></name> <name><surname>de Jong</surname> <given-names>M. C. M.</given-names></name> <name><surname>Graat</surname> <given-names>E. A. M.</given-names></name> <name><surname>Vandenbroucke-Grauls</surname> <given-names>C. M. J. E.</given-names></name> <name><surname>de Roda Husman</surname> <given-names>A. M.</given-names></name></person-group> (<year>2015</year>). <article-title>Role of the environment in the transmission of antimicrobial resistance to humans: a review</article-title>. <source>Environ. Sci. Technol.</source> <volume>49</volume>, <fpage>11993</fpage>&#x2013;<lpage>12004</lpage>. doi: <pub-id pub-id-type="doi">10.1021/acs.est.5b02566</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jia</surname> <given-names>B. F.</given-names></name> <name><surname>Raphenya</surname> <given-names>A. R.</given-names></name> <name><surname>Alcock</surname> <given-names>B.</given-names></name> <name><surname>Waglechner</surname> <given-names>N.</given-names></name> <name><surname>Guo</surname> <given-names>P. Y.</given-names></name> <name><surname>Tsang</surname> <given-names>K. K.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database</article-title>. <source>Nucleic Acids Res.</source> <volume>45</volume>, <fpage>D566</fpage>&#x2013;<lpage>D573</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkw1004</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>D. H.</given-names></name> <name><surname>Liu</surname> <given-names>C. M.</given-names></name> <name><surname>Luo</surname> <given-names>R. B.</given-names></name> <name><surname>Sadakane</surname> <given-names>K.</given-names></name> <name><surname>Lam</surname> <given-names>T. W.</given-names></name></person-group> (<year>2015</year>). <article-title>MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph</article-title>. <source>Bioinformatics</source> <volume>31</volume>, <fpage>1674</fpage>&#x2013;<lpage>1676</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btv033</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>M. M.</given-names></name> <name><surname>Ray</surname> <given-names>P.</given-names></name> <name><surname>Knowlton</surname> <given-names>K. F.</given-names></name> <name><surname>Pruden</surname> <given-names>A.</given-names></name> <name><surname>Xia</surname> <given-names>K.</given-names></name> <name><surname>Teets</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Fate of pirlimycin and antibiotic resistance genes in dairy manure slurries in response to temperature and pH adjustment</article-title>. <source>Sci. Total Environ.</source> <volume>710</volume>:<fpage>136310</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2019.136310</pub-id>, PMID: <pub-id pub-id-type="pmid">32050366</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>H.</given-names></name> <name><surname>Wang</surname> <given-names>X.</given-names></name> <name><surname>Tan</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>Q.</given-names></name> <name><surname>Zhang</surname> <given-names>C.</given-names></name> <name><surname>Wei</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Coconut shell and its biochar as fertilizer amendment applied with organic fertilizer: efficacy and course of actions on eliminating antibiotic resistance genes in agricultural soil</article-title>. <source>J. Hazard. Mater.</source> <volume>437</volume>:<fpage>129322</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jhazmat.2022.129322</pub-id>, PMID: <pub-id pub-id-type="pmid">35728320</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>B.</given-names></name> <name><surname>Yang</surname> <given-names>Y.</given-names></name> <name><surname>Ma</surname> <given-names>L.</given-names></name> <name><surname>Ju</surname> <given-names>F.</given-names></name> <name><surname>Guo</surname> <given-names>F.</given-names></name> <name><surname>Tiedje</surname> <given-names>J. M.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes</article-title>. <source>ISME J.</source> <volume>9</volume>, <fpage>2490</fpage>&#x2013;<lpage>2502</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2015.59</pub-id>, PMID: <pub-id pub-id-type="pmid">25918831</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>P.</given-names></name> <name><surname>Jia</surname> <given-names>S.</given-names></name> <name><surname>He</surname> <given-names>X.</given-names></name> <name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Ye</surname> <given-names>L.</given-names></name></person-group> (<year>2017</year>). <article-title>Different impacts of manure and chemical fertilizers on bacterial community structure and antibiotic resistance genes in arable soils</article-title>. <source>Chemosphere</source> <volume>188</volume>, <fpage>455</fpage>&#x2013;<lpage>464</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.chemosphere.2017.08.162</pub-id>, PMID: <pub-id pub-id-type="pmid">28898777</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Sun</surname> <given-names>H.</given-names></name> <name><surname>Zhang</surname> <given-names>M.</given-names></name> <name><surname>Liu</surname> <given-names>Y.</given-names></name></person-group> (<year>2019</year>). <article-title>Dynamics of microbial community and tetracycline resistance genes in biological nutrient removal process</article-title>. <source>J. Environ. Manag.</source> <volume>238</volume>, <fpage>84</fpage>&#x2013;<lpage>91</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jenvman.2019.02.123</pub-id>, PMID: <pub-id pub-id-type="pmid">30849601</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mart&#x00ED;nez</surname> <given-names>J. L.</given-names></name> <name><surname>Baquero</surname> <given-names>F.</given-names></name> <name><surname>Andersson</surname> <given-names>D. I.</given-names></name></person-group> (<year>2007</year>). <article-title>Predicting antibiotic resistance</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>5</volume>, <fpage>958</fpage>&#x2013;<lpage>965</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nrmicro1796</pub-id></citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Munk</surname> <given-names>P.</given-names></name> <name><surname>Knudsen</surname> <given-names>B. E.</given-names></name> <name><surname>Lukjancenko</surname> <given-names>O.</given-names></name> <name><surname>Duarte</surname> <given-names>A. S. R.</given-names></name> <name><surname>Van Gompel</surname> <given-names>L.</given-names></name> <name><surname>Luiken</surname> <given-names>R. E. C.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries</article-title>. <source>Nat. Microbiol.</source> <volume>3</volume>, <fpage>898</fpage>&#x2013;<lpage>908</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41564-018-0192-9</pub-id></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>N&#x00F5;lvak</surname> <given-names>H.</given-names></name> <name><surname>Truu</surname> <given-names>M.</given-names></name> <name><surname>Kanger</surname> <given-names>K.</given-names></name> <name><surname>Tampere</surname> <given-names>M.</given-names></name> <name><surname>Espenberg</surname> <given-names>M.</given-names></name> <name><surname>Loit</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Inorganic and organic fertilizers impact the abundance and proportion of antibiotic resistance and integron-integrase genes in agricultural grassland soil</article-title>. <source>Sci. Total Environ.</source> <volume>562</volume>, <fpage>678</fpage>&#x2013;<lpage>689</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2016.04.035</pub-id>, PMID: <pub-id pub-id-type="pmid">27115621</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nowakowska</surname> <given-names>M.</given-names></name> <name><surname>van Staa</surname> <given-names>T.</given-names></name> <name><surname>M&#x00F6;lter</surname> <given-names>A.</given-names></name> <name><surname>Ashcroft</surname> <given-names>D. M.</given-names></name> <name><surname>Tsang</surname> <given-names>J. Y.</given-names></name> <name><surname>White</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Antibiotic choice in UK general practice: rates and drivers of potentially inappropriate antibiotic prescribing</article-title>. <source>J. Antimicrob. Chemother</source> <volume>74</volume>, <fpage>3371</fpage>&#x2013;<lpage>3378</lpage>. doi: <pub-id pub-id-type="doi">10.1093/jac/dkz345</pub-id>, PMID: <pub-id pub-id-type="pmid">31430365</pub-id></citation></ref>
<ref id="ref29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Patro</surname> <given-names>R.</given-names></name> <name><surname>Duggal</surname> <given-names>G.</given-names></name> <name><surname>Love</surname> <given-names>M. I.</given-names></name> <name><surname>Irizarry</surname> <given-names>R. A.</given-names></name> <name><surname>Kingsford</surname> <given-names>C.</given-names></name></person-group> (<year>2017</year>). <article-title>Salmon provides fast and bias-aware quantification of transcript expression</article-title>. <source>Nat. Methods</source> <volume>14</volume>, <fpage>417</fpage>&#x2013;<lpage>419</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.4197</pub-id>, PMID: <pub-id pub-id-type="pmid">28263959</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peng</surname> <given-names>Y.</given-names></name> <name><surname>Leung</surname> <given-names>H. C. M.</given-names></name> <name><surname>Yiu</surname> <given-names>S. M.</given-names></name> <name><surname>Chin</surname> <given-names>F. Y. L.</given-names></name></person-group> (<year>2012</year>). <article-title>IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth</article-title>. <source>Bioinformatics</source> <volume>28</volume>, <fpage>1420</fpage>&#x2013;<lpage>1428</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/bts174</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="book"><person-group person-group-type="author"><collab id="coll1">R Core Team</collab></person-group> (<year>2022</year>). <source>R: A language and environment for statistical computing</source>. <publisher-name>R Foundation for Statistical Computing</publisher-name>, <publisher-loc>Vienna, Austria</publisher-loc>.</citation></ref>
<ref id="ref32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ren</surname> <given-names>Y.</given-names></name> <name><surname>Yu</surname> <given-names>G.</given-names></name> <name><surname>Shi</surname> <given-names>C. P.</given-names></name> <name><surname>Liu</surname> <given-names>L. M.</given-names></name> <name><surname>Guo</surname> <given-names>Q.</given-names></name> <name><surname>Han</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Majorbio cloud: a one-stop, comprehensive bioinformatic platform for multi-omics analyses</article-title>. <source>iMeta</source> <volume>1</volume>:<fpage>e12</fpage>. doi: <pub-id pub-id-type="doi">10.1002/imt2.12</pub-id></citation></ref>
<ref id="ref33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ren</surname> <given-names>C. J.</given-names></name> <name><surname>Zhang</surname> <given-names>X. Y.</given-names></name> <name><surname>Zhang</surname> <given-names>S. H.</given-names></name> <name><surname>Wang</surname> <given-names>J. Y.</given-names></name> <name><surname>Xu</surname> <given-names>M. P.</given-names></name> <name><surname>Guo</surname> <given-names>Y. X.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Altered microbial CAZyme families indicated dead biomass decomposition following afforestation</article-title>. <source>Soil Biol. Biochem.</source> <volume>160</volume>:<fpage>108362</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2021.108362</pub-id></citation></ref>
<ref id="ref34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Riber</surname> <given-names>L.</given-names></name> <name><surname>Poulsen</surname> <given-names>P. H. B.</given-names></name> <name><surname>Al-Soud</surname> <given-names>W. A.</given-names></name> <name><surname>Skov Hansen</surname> <given-names>L. B.</given-names></name> <name><surname>Bergmark</surname> <given-names>L.</given-names></name> <name><surname>Brejnrod</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Exploring the immediate and long-term impact on bacterial communities in soil amended with animal and urban organic waste fertilizers using pyrosequencing and screening for horizontal transfer of antibiotic resistance</article-title>. <source>FEMS Microbiol. Ecol.</source> <volume>90</volume>, <fpage>206</fpage>&#x2013;<lpage>224</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1574-6941.12403</pub-id>, PMID: <pub-id pub-id-type="pmid">25087596</pub-id></citation></ref>
<ref id="ref35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rosenkrantz</surname> <given-names>L.</given-names></name> <name><surname>Amram</surname> <given-names>O.</given-names></name> <name><surname>Caudell</surname> <given-names>M. A.</given-names></name> <name><surname>Schuurman</surname> <given-names>N.</given-names></name> <name><surname>Call</surname> <given-names>D. R.</given-names></name></person-group> (<year>2019</year>). <article-title>Spatial relationships between small-holder farms coupled with livestock management practices are correlated with the distribution of antibiotic resistant bacteria in northern Tanzania</article-title>. <source>One Health</source> <volume>8</volume>:<fpage>100097</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.onehlt.2019.100097</pub-id>, PMID: <pub-id pub-id-type="pmid">31249856</pub-id></citation></ref>
<ref id="ref36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rousk</surname> <given-names>J.</given-names></name> <name><surname>Baath</surname> <given-names>E.</given-names></name> <name><surname>Brookes</surname> <given-names>P. C.</given-names></name> <name><surname>Lauber</surname> <given-names>C. L.</given-names></name> <name><surname>Lozupone</surname> <given-names>C.</given-names></name> <name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>Soil bacterial and fungal communities across a pH gradient in an arable soil</article-title>. <source>ISME J.</source> <volume>4</volume>, <fpage>1340</fpage>&#x2013;<lpage>1351</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2010.58</pub-id>, PMID: <pub-id pub-id-type="pmid">20445636</pub-id></citation></ref>
<ref id="ref37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sanderson</surname> <given-names>H.</given-names></name> <name><surname>Fricker</surname> <given-names>C.</given-names></name> <name><surname>Brown</surname> <given-names>R. S.</given-names></name> <name><surname>Majury</surname> <given-names>A.</given-names></name> <name><surname>Liss</surname> <given-names>S. N.</given-names></name></person-group> (<year>2016</year>). <article-title>Antibiotic resistance genes as an emerging environmental contaminant</article-title>. <source>Environ. Rev.</source> <volume>24</volume>, <fpage>205</fpage>&#x2013;<lpage>218</lpage>. doi: <pub-id pub-id-type="doi">10.1139/er-2015-0069</pub-id></citation></ref>
<ref id="ref38"><citation citation-type="book"><person-group person-group-type="author"><collab id="coll2">Soil Survey Staff</collab></person-group> (<year>2010</year>). <source>Keys to soil taxonomy</source>. <edition>11</edition>. <publisher-name>USDA-NRCS</publisher-name>: <publisher-loc>Washington D.C.</publisher-loc></citation></ref>
<ref id="ref39"><citation citation-type="book"><person-group person-group-type="editor"><name><surname>Sparks</surname> <given-names>D. L.</given-names></name></person-group> (Ed.) (<year>1996</year>). <source>Methods of soil analysis Part 3: chemical methods</source>. <publisher-name>Soil Science Society of America, American Society of Agronomy</publisher-name>: <publisher-loc>Madison</publisher-loc>.</citation></ref>
<ref id="ref40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Su</surname> <given-names>J. Q.</given-names></name> <name><surname>Wei</surname> <given-names>B.</given-names></name> <name><surname>Ou-Yang</surname> <given-names>W. Y.</given-names></name> <name><surname>Huang</surname> <given-names>F. Y.</given-names></name> <name><surname>Zhao</surname> <given-names>Y.</given-names></name> <name><surname>Xu</surname> <given-names>H. J.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Antibiotic resistome and its association with bacterial communities during sewage sludge composting</article-title>. <source>Environ. Sci. Technol.</source> <volume>49</volume>, <fpage>7356</fpage>&#x2013;<lpage>7363</lpage>. doi: <pub-id pub-id-type="doi">10.1021/acs.est.5b01012</pub-id>, PMID: <pub-id pub-id-type="pmid">26018772</pub-id></citation></ref>
<ref id="ref41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Udikovic-Kolic</surname> <given-names>N.</given-names></name> <name><surname>Wichmann</surname> <given-names>F.</given-names></name> <name><surname>Broderick</surname> <given-names>N. A.</given-names></name> <name><surname>Handelsman</surname> <given-names>J.</given-names></name></person-group> (<year>2014</year>). <article-title>Bloom of resident antibiotic-resistant bacteria in soil following manure fertilization</article-title>. <source>PNAS</source> <volume>111</volume>, <fpage>15202</fpage>&#x2013;<lpage>15207</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.1409836111</pub-id>, PMID: <pub-id pub-id-type="pmid">25288759</pub-id></citation></ref>
<ref id="ref42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Urra</surname> <given-names>J.</given-names></name> <name><surname>Alkorta</surname> <given-names>I.</given-names></name> <name><surname>Lanz&#x00E9;n</surname> <given-names>A.</given-names></name> <name><surname>Mijangos</surname> <given-names>I.</given-names></name> <name><surname>Garbisu</surname> <given-names>C.</given-names></name></person-group> (<year>2019</year>). <article-title>The application of fresh and composted horse and chicken manure affects soil quality, microbial composition and antibiotic resistance</article-title>. <source>Appl. Soil Ecol.</source> <volume>135</volume>, <fpage>73</fpage>&#x2013;<lpage>84</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.apsoil.2018.11.005</pub-id></citation></ref>
<ref id="ref43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Van Goethem</surname> <given-names>M. W.</given-names></name> <name><surname>Pierneef</surname> <given-names>R.</given-names></name> <name><surname>Bezuidt</surname> <given-names>O. K. I.</given-names></name> <name><surname>Van De Peer</surname> <given-names>Y.</given-names></name> <name><surname>Cowan</surname> <given-names>D. A.</given-names></name> <name><surname>Makhalanyane</surname> <given-names>T. P.</given-names></name></person-group> (<year>2018</year>). <article-title>A reservoir of &#x2018;historical&#x2019; antibiotic resistance genes in remote pristine Antarctic soils</article-title>. <source>Microbiome</source> <volume>6</volume>:<fpage>40</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s40168-018-0424-5</pub-id>, PMID: <pub-id pub-id-type="pmid">29471872</pub-id></citation></ref>
<ref id="ref44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>W. Y.</given-names></name> <name><surname>Hou</surname> <given-names>Y. T.</given-names></name> <name><surname>Pan</surname> <given-names>W. H.</given-names></name> <name><surname>Vinay</surname> <given-names>N.</given-names></name> <name><surname>Mo</surname> <given-names>F.</given-names></name> <name><surname>Liao</surname> <given-names>Y. C.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Continuous application of conservation tillage affects in situ N<sub>2</sub>O emissions and nitrogen cycling gene abundances following nitrogen fertilization</article-title>. <source>Soil Biol. Biochem.</source> <volume>157</volume>:<fpage>108239</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2021.108239</pub-id></citation></ref>
<ref id="ref45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Z. T.</given-names></name> <name><surname>Li</surname> <given-names>T.</given-names></name> <name><surname>Wen</surname> <given-names>X. X.</given-names></name> <name><surname>Liu</surname> <given-names>Y.</given-names></name> <name><surname>Han</surname> <given-names>J.</given-names></name> <name><surname>Liao</surname> <given-names>Y. C.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Fungal communities in rhizosphere soil under conservation tillage shift in response to plant growth</article-title>. <source>Front. Microbiol.</source> <volume>8</volume>:<fpage>1301</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2017.01301</pub-id>, PMID: <pub-id pub-id-type="pmid">28744278</pub-id></citation></ref>
<ref id="ref46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>F.</given-names></name> <name><surname>Stedtfeld</surname> <given-names>R. D.</given-names></name> <name><surname>Kim</surname> <given-names>O. S.</given-names></name> <name><surname>Chai</surname> <given-names>B.</given-names></name> <name><surname>Yang</surname> <given-names>L.</given-names></name> <name><surname>Stedtfeld</surname> <given-names>T. M.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Influence of soil characteristics and proximity to antarctic research stations on abundance of antibiotic resistance genes in soils</article-title>. <source>Environ. Sci. Technol.</source> <volume>50</volume>, <fpage>12621</fpage>&#x2013;<lpage>12629</lpage>. doi: <pub-id pub-id-type="doi">10.1021/acs.est.6b02863</pub-id>, PMID: <pub-id pub-id-type="pmid">27797533</pub-id></citation></ref>
<ref id="ref47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>F.</given-names></name> <name><surname>Sun</surname> <given-names>R.</given-names></name> <name><surname>Hu</surname> <given-names>H.</given-names></name> <name><surname>Duan</surname> <given-names>G.</given-names></name> <name><surname>Meng</surname> <given-names>L.</given-names></name> <name><surname>Qiao</surname> <given-names>M.</given-names></name></person-group> (<year>2022b</year>). <article-title>The overlap of soil and vegetable microbes drives the transfer of antibiotic resistance genes from manure-amended soil to vegetables</article-title>. <source>Sci. Total Environ.</source> <volume>828</volume>:<fpage>154463</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2022.154463</pub-id>, PMID: <pub-id pub-id-type="pmid">35276164</pub-id></citation></ref>
<ref id="ref9001"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>T. T.</given-names></name> <name><surname>Sun</surname> <given-names>S. L.</given-names></name> <name><surname>Xu</surname> <given-names>Y. X.</given-names></name> <name><surname>Waigi</surname> <given-names>M. G.</given-names></name> <name><surname>Odinga</surname> <given-names>E. S.</given-names></name> <name><surname>Vasilyeva</surname> <given-names>G. K.</given-names></name></person-group> (<year>2022a</year>). <article-title>Nitrogen regulates the distribution of antibiotic resistance genes in the soil-vegetable system</article-title>. <source>Front. Microbiol.</source> <volume>13</volume>:<fpage>848750</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2022.848750</pub-id></citation></ref>
<ref id="ref48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>W. Y.</given-names></name> <name><surname>Yuan</surname> <given-names>J. C.</given-names></name> <name><surname>Gao</surname> <given-names>S. M.</given-names></name> <name><surname>Li</surname> <given-names>T.</given-names></name> <name><surname>Li</surname> <given-names>Y. J.</given-names></name> <name><surname>Vinay</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Conservation tillage enhances crop productivity and decreases soil nitrogen losses in a rainfed agroecosystem of the Loess Plateau, China</article-title>. <source>J. Clean. Prod.</source> <volume>274</volume>:<fpage>122854</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jclepro.2020.122854</pub-id></citation></ref>
<ref id="ref49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>J. Y.</given-names></name> <name><surname>Li</surname> <given-names>Z. T.</given-names></name> <name><surname>Guo</surname> <given-names>S. M.</given-names></name> <name><surname>Li</surname> <given-names>K. J.</given-names></name> <name><surname>Xu</surname> <given-names>P. S.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Antibiotics and antibiotic resistance genes in agricultural soils: a systematic analysis</article-title>. <source>Crit. Rev. Env. Sci. Technol.</source> <volume>53</volume>, <fpage>847</fpage>&#x2013;<lpage>864</lpage>. doi: <pub-id pub-id-type="doi">10.1080/10643389.2022.2094693</pub-id></citation></ref>
<ref id="ref50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiang</surname> <given-names>Q.</given-names></name> <name><surname>Chen</surname> <given-names>Q. L.</given-names></name> <name><surname>Zhu</surname> <given-names>D.</given-names></name> <name><surname>An</surname> <given-names>X. L.</given-names></name> <name><surname>Yang</surname> <given-names>X. R.</given-names></name> <name><surname>Su</surname> <given-names>J. Q.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Spatial and temporal distribution of antibiotic resistomes in a peri-urban area is associated significantly with anthropogenic activities</article-title>. <source>Environ. Pollut.</source> <volume>235</volume>, <fpage>525</fpage>&#x2013;<lpage>533</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envpol.2017.12.119</pub-id>, PMID: <pub-id pub-id-type="pmid">29324382</pub-id></citation></ref>
<ref id="ref51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xie</surname> <given-names>W. Y.</given-names></name> <name><surname>Yuan</surname> <given-names>S. T.</given-names></name> <name><surname>Xu</surname> <given-names>M. G.</given-names></name> <name><surname>Yang</surname> <given-names>X. P.</given-names></name> <name><surname>Shen</surname> <given-names>Q. R.</given-names></name> <name><surname>Zhang</surname> <given-names>W. W.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Long-term effects of manure and chemical fertilizers on soil antibiotic resistome</article-title>. <source>Soil Biol. Biochem.</source> <volume>122</volume>, <fpage>111</fpage>&#x2013;<lpage>119</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2018.04.009</pub-id></citation></ref>
<ref id="ref52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xun</surname> <given-names>W.</given-names></name> <name><surname>Zhao</surname> <given-names>J.</given-names></name> <name><surname>Xue</surname> <given-names>C.</given-names></name> <name><surname>Zhang</surname> <given-names>G.</given-names></name> <name><surname>Ran</surname> <given-names>W.</given-names></name> <name><surname>Wang</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Significant alteration of soil bacterial communities and organic carbon decomposition by different long-term fertilization management conditions of extremely low-productivity arable soil in South China</article-title>. <source>Environ. Microbiol.</source> <volume>18</volume>, <fpage>1907</fpage>&#x2013;<lpage>1917</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1462-2920.13098</pub-id>, PMID: <pub-id pub-id-type="pmid">26486414</pub-id></citation></ref>
<ref id="ref53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>L. Y.</given-names></name> <name><surname>Zhou</surname> <given-names>S. Y. D.</given-names></name> <name><surname>Lin</surname> <given-names>C. S.</given-names></name> <name><surname>Huang</surname> <given-names>X. R.</given-names></name> <name><surname>Neilson</surname> <given-names>R.</given-names></name> <name><surname>Yang</surname> <given-names>X. R.</given-names></name></person-group> (<year>2022</year>). <article-title>Effects of biofertilizer on soil microbial diversity and antibiotic resistance genes</article-title>. <source>Sci. Total Environ.</source> <volume>820</volume>:<fpage>153170</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2022.153170</pub-id>, PMID: <pub-id pub-id-type="pmid">35051473</pub-id></citation></ref>
<ref id="ref54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>N.</given-names></name> <name><surname>Juneau</surname> <given-names>P.</given-names></name> <name><surname>Huang</surname> <given-names>R.</given-names></name> <name><surname>He</surname> <given-names>Z.</given-names></name> <name><surname>Sun</surname> <given-names>B.</given-names></name> <name><surname>Zhou</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Coexistence between antibiotic resistance genes and metal resistance genes in manure-fertilized soils</article-title>. <source>Geoderma</source> <volume>382</volume>:<fpage>114760</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.geoderma.2020.114760</pub-id></citation></ref>
<ref id="ref55"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>Z. L.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Han</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>J. W.</given-names></name> <name><surname>Liu</surname> <given-names>G. F.</given-names></name> <name><surname>Lu</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Nutrients, heavy metals and microbial communities co-driven distribution of antibiotic resistance genes in adjacent environment of mariculture</article-title>. <source>Environ. Pollut.</source> <volume>220</volume>, <fpage>909</fpage>&#x2013;<lpage>918</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envpol.2016.10.075</pub-id>, PMID: <pub-id pub-id-type="pmid">27814984</pub-id></citation></ref>
<ref id="ref56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>X.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Zhu</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name></person-group> (<year>2019</year>). <article-title>Field-based evidence for enrichment of antibiotic resistance genes and mobile genetic elements in manure-amended vegetable soils</article-title>. <source>Sci. Total Environ.</source> <volume>654</volume>, <fpage>906</fpage>&#x2013;<lpage>913</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2018.10.446</pub-id>, PMID: <pub-id pub-id-type="pmid">30453260</pub-id></citation></ref>
<ref id="ref57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>J.</given-names></name> <name><surname>Guan</surname> <given-names>D.</given-names></name> <name><surname>Zhou</surname> <given-names>B.</given-names></name> <name><surname>Zhao</surname> <given-names>B.</given-names></name> <name><surname>Ma</surname> <given-names>M.</given-names></name> <name><surname>Qin</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Influence of 34-years of fertilization on bacterial communities in an intensively cultivated black soil in Northeast China</article-title>. <source>Soil Biol. Biochem.</source> <volume>90</volume>, <fpage>42</fpage>&#x2013;<lpage>51</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2015.07.005</pub-id></citation></ref>
<ref id="ref58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>Y. G.</given-names></name> <name><surname>Zhao</surname> <given-names>Y.</given-names></name> <name><surname>Zhu</surname> <given-names>D.</given-names></name> <name><surname>Gillings</surname> <given-names>M.</given-names></name> <name><surname>Penuelas</surname> <given-names>J.</given-names></name> <name><surname>Ok</surname> <given-names>Y. S.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Soil biota, antimicrobial resistance and planetary health</article-title>. <source>Environ. Int.</source> <volume>131</volume>:<fpage>105059</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.envint.2019.105059</pub-id></citation></ref>
</ref-list>
<glossary>
<def-list>
<title>Abbreviations</title>
<def-item><term>ZT, No tillage.</term></def-item>
<def-item><term>PT, Plough tillage</term></def-item>
<def-item><term>ARG, Antibiotic resistance genes</term></def-item>
<def-item><term>SOC, Soil organic carbon</term></def-item>
<def-item><term>TN, Soil total nitrogen</term></def-item>
<def-item><term>Zn, zinc.</term></def-item>
</def-list>
</glossary>
<fn-group>
<fn id="fn0001"><p><sup>1</sup><ext-link xlink:href="https://github.com/jstjohn/SeqPrep" ext-link-type="uri">https://github.com/jstjohn/SeqPrep</ext-link></p></fn>
<fn id="fn0002"><p><sup>2</sup><ext-link xlink:href="https://github.com/OpenGene/fastp" ext-link-type="uri">https://github.com/OpenGene/fastp</ext-link></p></fn>
<fn id="fn0003"><p><sup>3</sup><ext-link xlink:href="https://github.com/voutcn/megahit" ext-link-type="uri">https://github.com/voutcn/megahit</ext-link></p></fn>
<fn id="fn0004"><p><sup>4</sup><ext-link xlink:href="https://ngs.csr.uky.edu/Newbler" ext-link-type="uri">https://ngs.csr.uky.edu/Newbler</ext-link></p></fn>
<fn id="fn0005"><p><sup>5</sup><ext-link xlink:href="http://metagene.cb.k.u-tokyo.ac.jp" ext-link-type="uri">http://metagene.cb.k.u-tokyo.ac.jp</ext-link></p></fn>
<fn id="fn0006"><p><sup>6</sup><ext-link xlink:href="http://soap.genomics.org.cn" ext-link-type="uri">http://soap.genomics.org.cn</ext-link></p></fn>
<fn id="fn0007"><p><sup>7</sup><ext-link xlink:href="http://arpcard.Mcmaster.ca" ext-link-type="uri">http://arpcard.Mcmaster.ca</ext-link></p></fn>
<fn id="fn0008"><p><sup>8</sup><ext-link xlink:href="https://github.com/bbuchfink/diamond" ext-link-type="uri">https://github.com/bbuchfink/diamond</ext-link></p></fn>
<fn id="fn0009"><p><sup>9</sup><ext-link xlink:href="https://cloud.majorbio.com/page/tools/" ext-link-type="uri">https://cloud.majorbio.com/page/tools/</ext-link></p></fn>
</fn-group>
</back>
</article>
