Rong Hu1
Nafang Zhou
Jun Wang- 1Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, College of Urban and Environmental Science, Northwest University, Xi’an, China
- 2Shaanxi Key Laboratory for Carbon Neutral Technology, Northwest University, Xi’an, China
Crop rotation enhances agricultural productivity and soil fertility but may also contribute to the accumulation of antibiotic resistance genes (ARGs). However, the changes in soil ARGs and their associated resistance mechanisms under different crop rotation regimes are not well understood. In this study, we employed metagenomics to comprehensively investigate soil ARGs under different crop rotation regimes and complexity. Our findings revealed that soil properties varied significantly with crop rotation regime and complexity. Specifically, soil pH and the total carbon/nitrogen ratio (C/N) were the highest in bare land (BL) and gradually decreased in the order non-rotation, simple rotation, and complex rotation systems. The composition of soil ARGs exhibited significant differentiation by crop rotation complexity. Furthermore, differential gene analysis identified four specific types of ARGs—glycopeptide, multidrug, fluoroquinolone, and macrolide-lincosamide-streptogramin B (MLSB)—and two resistance mechanisms—cellular protection and efflux pump. Notably, soil microbial biomass carbon, soil microbial biomass nitrogen, and soil organic carbon are significantly correlated with ARGs in complex crop rotation systems, whereas soil pH and C/N ratio show significant associations in BL. The C/N ratio was identified as the most relevant determinant for glycopeptide, multidrug, fluoroquinolone, and MLSB resistance genes. Overall, these findings elucidate key factors associated with ARGs under long-term crop rotation, thereby providing valuable insights into the influence of crop rotation regimes on soil ARGs and enhancing soil fertility by improving soil properties.
1 Introduction
ARGs have become ubiquitous in various environmental compartments including airborne particulate matter, wastewater treatment sludge (Xie et al., 2016), and soil, posing significant threats to public health through their widespread dissemination (Wang et al., 2023). In recent years, a diverse array of ARGs has been detected in agricultural soils worldwide (Ondon et al., 2021; Shen et al., 2023). Notably, as a critical reservoir for ARGs (Khalid et al., 2023), soil acquires resistance genes through dual pathways (Su et al., 2014; Wang et al., 2023): intrinsic microbial carriage (Han et al., 2023) and exogenous inputs from agricultural activities (Wang et al., 2018, 2023; Gao et al., 2023; Han et al., 2023). Fertilization and irrigation have been demonstrated to influence the abundance and dissemination of ARGs by altering microbial communities, nutrient dynamics, and cross-ecosystem transmission pathways (Xie et al., 2016; Wang et al., 2023). However, the long-term effects of crop rotation—a crucial anthropogenic intervention—on ARGs remain insufficiently understood.
Research has demonstrated that crop rotation regimes profoundly influence the soil ARGs (Zhang et al., 2021). The abundance of sulfonamide genes was found to be higher in maize fields compared to wheat fields under wheat–maize rotation systems (Xu et al., 2021). In monoculture cropping systems, the abundance of resistance genes also varied depending on the cropping system (Huang et al., 2021). Crop diversification typically enhances microbial abundance (Liu et al., 2024), providing potential microbial carriers for soil ARGs. Alterations in the composition and function of microbial communities—a crucial component of soil ecosystems—determine the distribution, dissemination, and evolution of soil ARGs (Zhang et al., 2023). Crop rotation systems alter the diversity and distribution of soil ARGs by cultivating a variety of plant species (Venter et al., 2016; Huang et al., 2021). Most existing studies focus on short-term effects in simple rotation systems. However, the cumulative effects in complex (multi-species, long-term) systems remain unclear, limiting our ability to predict the ecological risks of ARGs under realistic agricultural practices. Although complex and diverse crop rotation systems substantially enhance agricultural productivity and soil fertility (Bowles et al., 2020, 2022), they also pose the potential risk of accelerating the accumulation of ARGs. A comprehensive investigation into the impacts of long-term, complex rotation systems on ARGs is crucial for optimizing agricultural benefits while mitigating the associated ecological and public health risks.
Environmental factors are closely associated with ARGs. In surface soils with long-term pig manure application, heavy metals were found to be significantly positively correlated with the abundance of ARGs (Guo et al., 2018). Microbial carbon and nitrogen (MBC and MBN) have been identified as important predictors of changes in ARG patterns (Li et al., 2022). In greenhouse soils, total nitrogen and organic matter emerge as the two major soil properties affecting the abundance of ARGs (Wang et al., 2020). Moreover, the C/N ratio has been proven to regulate ARG removal during the composting process (Zhu et al., 2021). These studies highlight the influence of environmental factors on ARGs in different soils. However, the relationship between environmental factors and ARGs under crop rotation systems requires further investigation. Notably, most previous studies have relied on quantitative PCR and amplicon sequencing to determine the abundance and distribution of soil ARGs (Li et al., 2020). In contrast, metagenomic sequencing provides a more comprehensive approach by capturing the genomic diversity of environmental microbial communities, enabling the identification of novel ARGs and tracing their transmission pathways (Marić and Šikić, 2019). However, studies combining metagenomics with long-term field experiments to investigate ARGs in arid agricultural ecosystems have not yet been reported.
The Loess Plateau of China is a typical arid and semi-arid region. Its unique soil and water conditions may create distinctive distribution patterns of ARGs. Crop rotation and diversification are important agricultural practices in this region. However, existing studies in this region have primarily focused on the effects of crop rotation on soil fertility (Liu Y. et al., 2023), with limited attention to the impact of crop rotation on ARGs as emerging contaminants. Therefore, this study aims to: (1) evaluate the effects of different crop rotation systems on soil properties and microbial biomass; (2) investigate the abundance, diversity, and resistance mechanisms of ARGs under long-term complex crop rotation regimes using metagenomic analysis; and (3) identify the factors associated with ARGs under crop rotation regimes and complexity.
2 Materials and methods
2.1 Experimental site description and soil sampling
The crop rotation experiment was initiated in September 1984 in a wheat field, with the site located at the Loess Plateau Agricultural Ecology Experimental Station (35°12′N, 107°44′E) of the Chinese Academy of Sciences (Changwu County, Shaanxi Province, China) (Cai and Hao, 2015; Fu et al., 2019; Liu Y. et al., 2023). This site is situated in the East Asian continental monsoon climate zone, with an average annual temperature of 9.1°C and 580 mm of precipitation. The soil is classified as Heilutu silt loam (Calcarid Regosol according to the Food and Agriculture Organization of the United Nations classification system), with 45 g kg−1 sand, 656 g kg−1 silt and 309 g kg−1 clay at depth of 0–20 cm at the beginning of the experiment (Fu et al., 2019).
In this study, we set up five different crop rotation regimes: bare land (BL), winter wheat monoculture (W), pea–winter wheat–winter wheat–millet (PWWM), corn–winter wheat–winter wheat–millet (CWWM), and alfalfa (4 yr) –potato (1 yr) –winter wheat (3 yr) (A4PoW3). We calculated the rotation complexity index (RCI) for these regimes, which was defined as the square root of the number of cash and cover crop species (crop species richness) in a rotation multiplied by the length of the rotation (crop cover duration). The five rotation regimes were classified into four crop rotation complexity levels based on RCI: bare land (BL), non-rotation (W, RCI = 1), simple rotation (PWWM and CWWM, RCI = 3), and complex rotation (A4PoW3, RCI = 5) (Bowles et al., 2020, 2022). The experiment was arranged in a completely randomized block design with three replications over five crop rotations with a rotation length of 1–10 years (Fu et al., 2019; Liu Y. et al., 2023).
Samples were collected in September 2021 from a wheat field and individual soil samples were obtained by mixing six randomly collected soil cores with a depth of 0–20 cm from the middle row of each plot. After homogenization, soil samples were separated into two subsamples using a 2 mm sieve. One subsample was allocated for the analysis of soil properties and microbial biomass, while the other subsample were preserved at −80°C until microbial DNA could be extracted.
2.2 Analyses of soil properties and microbial biomass
Soil samples were air-dried for additional analysis of other soil properties after the soil water content (SWC) was determined. The soil pH was measured in a soil suspension (soil:water, 1:2.5, w/v) made in distilled water. The soil samples were ground to < 0.50 mm, and the soil was pretreated with 6 mol L– 1 HCl to remove inorganic carbon. The combustion method was used to determine the amount of soil organic carbon (SOC). The C/N ratio was consistent with previous studies (Liu et al., 2020). Measurement of potential carbon mineralization (PCM) was based on CO2–C accumulation after capillary rewetting and a 24-h incubation period (Haney et al., 2004). Particulate organic nitrogen (PON) was analyzed using the sodium hexametaphosphate method (Liu Y. et al., 2023; Liu Z. S. et al., 2023), and the modified incubation method was used to determine potential nitrogen mineralization (PNM) concentration (Haney et al., 2004). The ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3−-N) concentrations in the extracts were measured by a modified Griess-Ilosvay method using an automated analyzer (Wang et al., 2021). Soil microbial biomass carbon (MBC) and nitrogen (MBN) were measured by the chloroform-fumigation extraction method (Wu et al., 1990).
2.3 DNA extraction and bioinformatic analyses of metagenomic sequencing
The total DNA was extracted from soil samples (0.5 g each) using the Fast DNA®SPIN Kit (MP Biochemicals, Solon, United States) following the manufacturer’s procedures. Concentration and purity of extracted DNA was determined with TBS-380 and NanoDrop2000, respectively. DNA extract quality was checked on 1% agarose gel. Sequencing was performed by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China), generating 150 bp paired-end reads at enhanced depth. Reads that aligned with human or plant genomes were excluded, and remaining sequences were trimmed using Sickle and subjected to quality filtering to ensure data integrity. Raw metagenomic reads were converted to FASTQ format, followed by multiple filtering steps: adapter removal, read trimming, and elimination of low-quality reads. Clean reads were generated using fastp on the Majorbio Cloud Platform,1 and high-quality reads were assembled using MEGAHIT,2 which made use of succinct de Bruijn graphs, retaining contigs ≥ 300 bp for downstream analysis. Open reading frames (ORFs) in each assembled contig were predicted using Prodigal. ORFs ≥ 100 bp were retrieved and translated into amino acid sequences. A non-redundant gene catalog was constructed using CD-HIT3 with a sequence identity of 90 and 90% coverage. High-quality reads were aligned to the non-redundant gene catalogs to calculate gene abundance with 95% identity using SOAPaligner,4 and gene abundance in each sample was evaluated. We utilized per-million-base transcript (TPM) mapping readings to standardize abundance values in metagenomes. The non-redundant gene set was aligned against the Comprehensive Antibiotic Resistance Database (CARD v3.0.9) using Diamond5 to annotate ARGs with an e-value cutoff of 1e−5. Details of the 212 identified ARGs, including gene names, types, and resistance mechanisms, are listed in Supplementary Table S1.
2.4 Statistical analyses
Statistical analyses were conducted using the R software platform (v4.3.26). Soil properties, ARGs and their resistance mechanisms were subjected to multiple comparisons using analysis of variance (ANOVA), followed by multiple comparison of means through Tukey’s HSD test, implemented using the “multcomp” package (Hothorn et al., 2008). Principal component analysis (PCA) was performed to assess soil property distributions across a matrix of proportional parameters, visualized using the “FactoMineR” package (Lê et al., 2008). A circos plot, constructed using the “circlize” package (Gu et al., 2014), depicted the abundance of soil ARGs and their resistance mechanisms. The composition of soil ARGs was represented through principal coordinate analysis (PCoA) using the “ggplot2” package (Wickham, 2009) and multifactor ANOVA (Adonis) based on Bray-Curtis distance was performed by the “vegan” package (Dixon, 2003). The volcano plot was based on the “ropls” package (Thev, 2015) for pls-da analysis and visualized by “ggplot2” package. The Venn diagram was created with the “Venn diagram” package (Wen et al., 2023) and the “ggplot2” package on the basis of significantly differential genes screened by the volcano plot.
Redundancy analysis (RDA) with 999 permutations was conducted to examine ARG-soil property relationships, implemented in vegan. Mantel tests, also via “vegan,” assessed correlations between specific ARG types and resistance mechanisms (identified from volcano and Venn analyses) and soil properties. Random forest (RF) regression quantified the influence of soil properties on ARG abundance and resistance mechanisms. Pearson correlation coefficients between specific ARG types and soil properties were calculated and visualized using the “corrplot” package (Wei and Simko, 2024). The “Hmisc” and “igraph” packages were used to identify the potential associations (spearman’s correlation |r| > 0.7, p < 0.05) between the subtypes of soil ARGs and properties and the networks were illustrated with the collaborative platform Gephi.
3 Results
3.1 Variations in soil properties and microbial biomass under distinct crop rotation regimes
The ANOVA results indicated that most soil properties showed significant differences (p < 0.05) across different crop rotation regimes and complexity levels (Supplementary Table S2). Specifically, soil pH, SOC, PON, MBN, C/N ratio, and PCM/PNM differed significantly (p < 0.05) between the BL regimes and other crop rotation regimes. Similarly, significant differences (p < 0.05) in soil pH, MBN, and C/N ratio were observed between the W regimes and other crop rotation regimes (PWWM, CWWM, and A4PoW3). In comparison, the values of NO₃−-N, PCM/PNM, and MBN were the highest under the CWWM regime, with the PCM/PNM value almost double that under other treatments and exhibiting a significant difference (p < 0.05) compared with that under the A4PoW3 regime. Conversely, NH₄+-N and MBC showed the highest values under PWWM, with MBC being approximately twice the level observed in BL, demonstrating a significant difference (p < 0.05) relative to BL. Nevertheless, NO₃−-N and NH₄+-N did not show significant variations among the different rotation regimes. Furthermore, soil pH and C/N ratio were highest in BL, followed by those under non-rotation and simple rotation regimes, and lowest in the complex rotation regime.
Additionally, principal component analysis (PCA) revealed that the first two principal components accounted for 70.8% of the percent variability (PC1 = 49.9% and PC2 = 20.9%) in soil properties, with the measured soil properties indicators clustered in distinct groups (Figure 1A). The contribution of each variable to the principal components varied, with their interrelationships—reflected by the angles between vectors—supported by subsequent correlation analysis. Soil pH, C/N, PON, MBN, and SOC were the dominant contributors to Dim1, while NO₃−-N, NH₄+-N, and MBC exhibited stronger associations with Dim2. Notably, the C/N ratio was negatively correlated with most soil properties, except for soil pH and soil water content (SWC). SOC was significantly positively correlated with PON and MBN (Figure 1B). Overall, soil properties displayed clear differentiation across crop rotation regimes and complexity, with their inter-dependencies varying accordingly.
Figure 1. Principal component analysis (PCA) of soil properties under different crop rotation systems. (A) Ordinal distribution pattern. The same color indicates the same group. The horizontal axis represents the first principal component coordinate and its contribution rate, while the vertical axis represents the second principal component coordinate and its contribution rate. Ellipses denote 95% confidence intervals, with distinct ellipses corresponding to different treatment groups. (B) Contribution of each variable to the principal components. The depth of the arrow line color represents the degree of contribution. Darker hues indicate higher contributions. The angle of the arrow line represents the correlation. Acute angles (<90°) indicate strong positive correlations. Obtuse angles (>90°) reflect strong negative correlations. Right angles (≈90°) suggest no significant correlation. BL, bare ground; W, continuous winter wheat; PWWM, pea-winter wheat-winter wheat-wheat; CWWM, corn-winter wheat-winter wheat-wheat; A4PoW3, alfalfa-potato-winter wheat. pH, soil pH; SWC, soil water content; NH4+-N, ammonium nitrogen; NO3−-N, nitrate nitrogen; and SOC, total organic carbon; MBC, soil microbial biomass carbon; MBN, soil microbial biomass nitrogen; C/N, total carbon to nitrogen ratio; PCM/PNM, potential carbon mineralization/potential nitrogen mineralization ratio; PON, particulate organic nitrogen.
3.2 Distribution and composition of soil microbial antibiotic resistance genes and their resistance mechanisms
The diversity of soil ARG subtypes was assessed, revealing a total of 212 distinct ARG subtypes across five crop rotation regimes (Supplementary Table S3). These ARG subtypes were classified into 12 resistance categories, including aminocoumarin, aminoglycoside, beta-lactam, elfamycin, fluoroquinolone, glycopeptide, MLSB, multidrug, rifampin, sulfonamide, tetracycline, and others (Supplementary Table S4). The abundance of soil ARGs exhibited significant variations (p < 0.05) among different crop rotation regimes (Figure 2A). Notably, fluoroquinolone and multidrug resistance genes were the most abundant ARGs across all regimes examined (Supplementary Figure S1). To investigate the effect of different crop rotation regimes and complexity on the distribution of soil ARGs, PCoA was conducted using the abundance data of soil ARGs (Figure 3). The PCoA revealed distinct separation among the A4PoW3, PWWM, and BL regimes along the axes, indicating that crop rotation complexity significantly shapes the composition of soil ARGs. Furthermore, ADONIS analysis (R2 = 0.61, p < 0.05) confirmed the statistical significance of this separation. These findings demonstrate that crop rotation regimes significantly affected the distribution and composition of soil ARGs.
Figure 2. Abundance and composition of soil ARGs and their resistance mechanisms. (A) Circos plots showing the distribution of 12 resistance gene types under different crop rotation systems. The thickness of each ribbon represents the abundance of crop rotation regimes assigned to different resistance gene types. (B) Circos plots showing the distribution of 7 resistance mechanisms under different crop rotation systems. The thickness of each ribbon represents the abundance of crop rotation regimes assigned to different 7 resistance mechanisms. A4PoW3, alfalfa-potato-winter wheat. BL, bare land; CWWM, corn-winter wheat-winter wheat-millet; PWWM, pea-winter wheat-winter wheat-millet; W, continuous winter wheat.
Figure 3. Principal Coordinate Analysis (PCoA) of ARGs under different crop rotation systems. PCOA1 and PCOA2 are two principal coordinate components. Each point in the figure represents a sample, and the color of the point represents the group of the sample. The R2 (ADONIS) quantifies the explained variance (higher values indicate greater explanatory power). p < 0.05 indicates a significant difference. A4PoW3, alfalfa-potato-winter wheat. BL, bare land; CWWM, corn-winter wheat-winter wheat-millet; PWWM, pea-winter wheat-winter wheat-millet; W, continuous winter wheat.
Compared to BL, the number of significantly differential genes (p < 0.05) varied across treatments: W had 11 genes, PWWM had 18 genes, CWWM had 15 genes, and A4PoW3 had 30 genes (Supplementary Figure S2). This indicated that various treatments had distinct impacts on ARGs. In the A4PoW3 regime, differentially abundant genes were most numerous, with a higher count of specific ARG subtypes compared to other regimes (Supplementary Table S5), with the largest number of ARGs in the glycopeptide and multidrug categories, followed by those in the fluoroquinolone category. Some differential genes, such as oleR and oleB (both MLSB), were consistently observed across distinct contrast groups, including W vs. BL, PWWM vs. BL, CWWM vs. BL, and A4PoW3 vs. BL (Supplementary Figure S3).
In this study, resistance mechanisms under different crop rotation regimes were analyzed, including antibiotic target replacement, antibiotic deactivation, antibiotic target modification, antibiotic target protection, cellular protection, and efflux pump. The results showed that cellular protection and efflux pump were the two predominant resistance mechanisms across all crop rotation regimes (Figure 2B) and the abundance of soil ARGs associated with these resistance mechanisms showed significant differences among the different treatment groups (p < 0.05; Supplementary Table S6). In A4PoW3, the abundance of genes related to cellular protection and efflux pump was significantly higher than that in other treatments, with multidrug resistance genes particularly enriched in the efflux pump mechanism (Supplementary Table S1).
3.3 Correlations between soil properties, microbial biomass, and antibiotic resistance genes and their resistance mechanisms
The RDA was employed to assess the influence of soil properties and microbial biomass on ARG abundance (Supplementary Figure S4). The results indicated that C/N ratio, pH, MBC, and MBN can significantly explain the variation in the abundance of resistance genes, and their effects on the abundance of soil ARGs were statistically significant (Supplementary Table S7). MBC, MBN, and SOC were significantly correlated with soil ARGs under A4PoW3, and the correlation between soil C/N ratio and pH and soil ARGs was particularly significant under BL. Our findings reveal that differential crop rotation regimes and complexity potentially modify edaphic factors, which are significantly associated with ARG abundance and dispersal characteristics in soils.
A Mantel test (Supplementary Tables S8, S9) was conducted to examine the correlation between specific ARG types (glycopeptide, multidrug, fluoroquinolones, and MLSB) and resistance mechanisms (efflux pump and cellular protection) with various soil properties and microbial biomass—soil pH, SWC, NH4+-N, NO3−-N, SOC, MBC, MBN, C/N ratio, and PCM/PNM—with analyses performed across different crop rotation regimes. The results revealed soil pH, SOC, C/N ratio, and MBN were significantly correlated with specific ARG types and resistance mechanisms (p < 0.05; Figure 4). Notably, of these soil properties, soil C/N consistently explained the highest proportion of variance in specific ARG types and resistance mechanisms, which was further confirmed by the RF analysis (Tables 1, 2). In crop rotation regimes, soil properties and microbial biomass were identified in ARG interaction networks (Supplementary Figure S5). C/N ratio recorded the highest node connectivity in the ARG network (node degree = 88; Supplementary Table S10), indicating the crucial role of the C/N ratio in maintaining the stability of the resistance gene network. Overall, the complexity of different crop rotation regimes, particularly their impact on soil C/N ratio, may be associated with the abundance and distribution of ARG types and resistance mechanisms.
Figure 4. Relationship between soil properties and ARGs and their resistance mechanisms. (A) Correlation analysis between soil properties and ARGs based on Mantel test. Mantel r and p values (*p < 0.05; **p < 0.01), as well as Pearson r values, are indicated by the color and width of the connection lines as specified in the legend. (B) Correlation analysis between soil properties and resistance mechanisms based on Mantel test. Mantel r and p values (*p < 0.05; **p < 0.01), as well as Pearson r values, are indicated by the color and width of the connection lines as specified in the legend. pH, soil pH; SWC, soil water content; NH4+-N, ammonium nitrogen; NO3−-N, nitrate nitrogen; and SOC, total organic carbon; MBC, soil microbial biomass carbon; MBN, soil microbial biomass nitrogen; C/N, total carbon to nitrogen ratio; PCM/PNM, potential carbon mineralization/potential nitrogen mineralization ratio; PON, particulate organic nitrogen.
Table 1. Random forest mean predictor importance of the soil properties for the abundance of the representative types of antibiotic resistance genes.
Table 2. Random forest mean predictor importance of the soil properties for the abundance of the representative types of antibiotic resistance mechanisms.
4 Discussion
4.1 Effects of crop rotation regime and complexity on soil properties and microbial biomass
Crop rotation plays a critical role in enhancing soil fertility (Giacometti et al., 2021). Crop rotation effectively improves the chemical and biological characteristics of soil through mechanisms such as plant root growth, rhizosphere microbial activity, and the accumulation of organic matter (McDaniel and Grandy, 2016; Giacometti et al., 2021). Crop rotation regime and complexity significantly enhanced soil microbial activity and nitrogen mineralization. We found that PCM/PNM (9.95 ± 0.68) and MBN (63.83 ± 1.19) in the CWWM regime were significantly higher values than in other treatment regimes and MBC value was the highest in the PWWM group. This phenomenon may be attributed to maize, a carbon-rich and nitrogen-demanding crop. When integrated into diversified rotations, maize enhances microbial activity and accelerates organic matter decomposition, thereby promoting nitrogen mineralization (McDaniel and Grandy, 2016). Additionally, the input of high C/N ratio residues (such as corn) drives microorganisms to prioritize the mineralization of nitrogen-containing compounds to meet their nitrogen demands and maintain a stable C/N ratio in crop rotation systems. This process suppresses the secretion of carbon-degrading enzymes, delaying carbon mineralization. Consequently, the difference in carbon and nitrogen mineralization rates becomes more pronounced, leading to a faster accumulation of PCM compared to PNM (Islam et al., 2022; Jiang et al., 2025). The introduction of leguminous plants enhances soil microbial activity by providing additional nitrogen sources and organic matter, thereby increasing the soil MBC content. This discovery provides direct evidence that legume crop rotation improves soil fertility. Although the contents of NH4+-N and NO3−-N varied across different cropping systems, these differences were not significant (p > 0.05). The increase in soil temperature and water content enhanced the mineralization of soil organic nitrogen (Fu et al., 2019). Combined with the lack of nitrogen uptake by crops, this leads to an increase in NH4+-N levels (Fu et al., 2019). Additionally, the application of legume residues and nitrogen fertilizer can also influence the NH4+-N content (Hu et al., 2023). Soil NO3−-N is lost through runoff and leaching, and high nitrogen fertilizer input affects NO3−-N levels as well (Fu et al., 2019; Hu et al., 2023). Overall, the variations in NH4+-N and NO3−-N were regulated by multiple factors, and the effects of different cropping systems were offset, resulting in no significant differences. Furthermore, soil pH were the highest in BL and the lowest in A4PoW3. Compared to the BL treatment, the A4PoW3 treatment introduced more diverse crops, and the root exudates of these crops mostly possessed acidic properties (Lv et al., 2023). The roots may have lowered the soil pH by increasing the release of acidic compounds, leading to differences in pH across the different crop rotations (Wen et al., 2017). Furthermore, changes in soil pH altered nutrient availability, ultimately affecting soil fertility and its capacity to support crop growth. In this study, long-term diversified crop rotation reduced the C/N ratio. Previous studies have shown that introducing legumes into long-term diversified planting may effect nitrogen balance, but it may not have a significant impact on carbon (Islam et al., 2022).
Additionally, different crop rotation regimes affect the soil environmental condition through the input of post-harvest crop residues and various root exudates secreted during the different growth stages of plants (Crotty et al., 2016), leading to the significant differences in individual physicochemical properties observed in the PCA. Moreover, strong interdependencies observed between the C/N ratio, other soil properties, and microbial biomass. C/N ratio influences the decomposition rate of organic matter, microbial community activity, and nutrient release and fixation, thereby affecting soil structural stability, nutrient availability, and pH buffering capacity (Wang et al., 2024).
4.2 Effects of crop rotation regimes on the dynamics of antibiotic resistance genes and their resistance mechanisms
Agricultural management practices are closely associated with the contamination of soil with microbial ARGs (Pruden et al., 2013). In this investigation, fluoroquinolone and multidrug resistance genes emerged as the most abundant ARG types. The widespread use of fluoroquinolones and multidrug antibiotics in agriculture (Sun et al., 2017; Huang et al., 2021; Wu et al., 2023; Zhang, 2024) has led to the adaptation of microbial communities to antibiotic selection pressure (Forsberg et al., 2014; Jechalke et al., 2014). This adaptation process leads to the selective enrichment and persistence of corresponding resistance genes in soil. Additionally, under diversified crop rotation systems, there are complex interactions between root exudates and microorganisms in the rhizosphere soil (Mishra et al., 2022). These interactions may promote the faster proliferation and transmission of some ARGs in the soil (Haichar et al., 2008). Fertilization reshapes the structure of rhizosphere microbial communities, driving the adaptive evolution of environmental microorganisms toward a higher resistance gene profile (Guan et al., 2025). These changes promote the expression and transmission of ARGs (Gao et al., 2023; Wang et al., 2023), which may lead to increased resistant microorganisms and ARG accumulation. The PCoA results showed a significant separation of ARGs among crop rotation regimes with varying complexity. Crop rotation indirectly shapes ARG diversity and abundance by changing the structure of soil microbial communities (Zhang et al., 2023) and affecting stability of their ecological networks (Fan et al., 2022).
Furthermore, different crop rotation regimes significantly affected the abundance of soil ARGs (p < 0.05). The enrichment of agricultural soil ARGs is influenced by the cultivation of different crop types (Cerqueira et al., 2019; Huang et al., 2021). Crop root exudates play an important role in this process by regulating the structure and function of rhizosphere microbial communities, indirectly promoting the growth of ARGs in the soil–plant system (Gao et al., 2023). In A4PoW3, the number of differential genes was the highest among all treatments. This likely resulted from the use of a more complex crop combination and extended cultivation duration, which resulted in richer and more complex root exudates and crop residue inputs (Bowles et al., 2022). These inputs may provide diverse carbon and nitrogen sources for the soil microbial communities, thereby enhancing microbial activity and metabolic process diversity. Such stimulation of microbial functions could significantly influence distribution patterns of ARGs. This phenomenon highlights the vital role of crop diversity in regulating the distribution of soil ARGs (Cerqueira et al., 2019). Additionally, differences in oleR and oleB across regimes indicates the transfer of relevant potential genes in the microbial community and exposure to selection pressure under different crop rotation regimes. This finding further emphasizes the dynamic role of microbial communities in the spread of soil ARGs (Zhang et al., 2022).
Our findings indicate that efflux pumps and cellular protection were significant resistance mechanisms under crop rotation. The resistance mechanisms efflux pump and cellular protection impart resistance to various antibiotics. Efflux pumps expel various types of antibiotic molecules from the cells (Zack et al., 2024), and the cellular protection mechanism protects numerous target proteins and nucleic acids from antibiotics (Darby et al., 2023). Efflux pumps and cellular protection were identified as the main resistance mechanisms (Alshehri et al., 2023) in the root nodule soils associated with the naturally growing Abutilon fruticosum. Root secretions of different crops affect the soil microbial communities under different crop rotation regimes (Li et al., 2025). In microbial communities, different resistance mechanisms are selected for through multiple interactions (Bottery et al., 2021), which increases the soil microbial antibiotic resistance (Xiao et al., 2023).
4.3 Contribution of soil properties and microbial biomass to soil microbial antibiotic resistance genes
Both RDA and the Mantel test results indicated that soil properties and microbial biomass were significantly correlated with the abundance and distribution of soil ARGs (p < 0.05). Soil properties and microbial biomass may impact the abundance and distribution of soil ARGs by affecting the transmission of soil ARGs regulated by the composition and species diversity of soil microbial communities (Wu et al., 2023). In this study, MBC, MBN, and SOC showed a more significant correlation with the distribution of soil ARGs under the A4PoW3 treatment, suggesting that the relationship between microbial communities and ARGs is more pronounced under a complex crop rotation regime. Soil pH, SOC, MBN, C/N were found to be significantly correlated with ARGs types as well as their resistance mechanisms (p < 0.05). Soil pH emerged as a critical driver of microbial community evolution. In general, bacteria were most abundant in neutral environments, and soil pH influences microbial selection by altering nutrient availability or affecting physiological activity. Certain pH values might apply direct pressure on bacterial cells and lead to the selection of certain bacterial populations, indirectly affecting the production of ARGs (Han et al., 2022). SOC serves as a carbon and nutrient source for bacteria and influence the abundance of soil ARGs by promoting the growth of ARG-harboring bacteria (Guo et al., 2020; Shen et al., 2023). Through further network analysis and prediction using random forests, identified that the C/N ratio might be the most relevant factor. The C/N ratio directly affects microbial nutritional status and metabolic efficiency (Gaudel et al., 2024). According to the stoichiometric decomposition theory (Chen et al., 2014), microorganisms tend to consume resources that meet their stoichiometric requirements to maintain their ideal elemental stoichiometric ratios. The C/N ratio may influence the activity and function of microbial communities by regulating microbial energy metabolism. The C/N ratio not only regulates microbial metabolic activities but may also further alter the efficiency of horizontal gene transfer by affecting resource utilization efficiency and microbial community structure. It enhances the ecological adaptability of certain microorganisms with competitive advantages in resource acquisition (Ding and Sun, 2025), promoting the gene exchange efficiency among microorganisms carrying specific resistance genes (Shen et al., 2023). In summary, the C/N ratio significantly affects microbial community structure and its ecological functions by regulating microbial energy metabolism, horizontal gene transfer efficiency, and carbon/nutrient competition. These mechanisms collectively explain why the C/N ratio becomes a key factor related to the distribution of ARGs.
In this study, the correlation between specific ARG types and soil properties as well as microbial biomass was analyzed using the Mantel test. Soil pH, SOC, MBN, C/N ratio were found to be significantly correlated with ARGs types as well as their resistance mechanisms (p < 0.05). These four factors may exhibited selective effects on specific ARG types. For example, they influences ARGs by promoting the growth of bacteria carrying specific ARGs (Zhang et al., 2024), shaping the resistome (Liu Z. S. et al., 2023), and enhancing the competitive advantage of specific ARG-carrying bacteria (Zheng et al., 2022). Future research could further explore the mechanisms of changes in microbial metabolic pathways under different C/N ratios, and how these changes further regulate soil functions and the spread of resistance genes by affecting microbial community dynamics and gene exchange efficiency.
5 Conclusion
In this study, we examined the effects of crop rotation regimes and complexity on soil properties, microbial biomass, and the distribution, composition, and resistance mechanisms of soil ARGs using metagenomics and multivariate statistical analyses. Our findings demonstrated that crop rotation regimes and complexity significantly influenced most soil properties and microbial biomass, resulting in distinct separation and significant differences among regimes. Similarly, the abundance and distribution of soil ARGs exhibited significant variation across crop regimes, reflecting the impact of rotation complexity. Furthermore, analysis of differently abundant genes further enabled the identification of specific ARG types and their associated resistance mechanisms. The results suggest that soil properties are correlated with the selection of ARGs and their resistance mechanisms, with the C/N ratio being the most significant correlated factor. Collectively, these findings offer a comprehensive insight into the changes in soil ARGs under long-term crop rotation regimes, highlighting the interplay between long-term crop rotation systems, soil environmental conditions, and microbial resistance.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/sra/PRJNA985043.
Author contributions
RH: Investigation, Writing – original draft. YL: Supervision, Writing – review & editing. MW: Methodology, Writing – review & editing. NZ: Methodology, Writing – review & editing. JW: Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Grant Nos. 42277322, 42407403), the China Postdoctoral Science Foundation (2023M732844), and the Shaanxi Province Science Foundation for Youth (2024JC-YBQN-0275).
Conflict of interest
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1603518/full#supplementary-material
Footnotes
1. ^https://github.com/OpenGene/fastp
2. ^https://github.com/voutcn/megahit
3. ^http://www.bioinformatics.org/cd-hit/
4. ^http://soap.genomics.org.cn/, version 2.21.
5. ^http://www.diamondsearch.org/index.php, version 0.8.35.
References
Alshehri, W. A., Abulfaraj, A. A., Alqahtani, M. D., Alomran, M. M., Alotaibi, N. M., Alwutayd, K., et al. (2023). Abundant resistome determinants in rhizosphere soil of the wild plant Abutilon fruticosum. AMB Express 13:92. doi: 10.1186/s13568-023-01597-w
Bottery, M. J., Pitchford, J. W., and Friman, V. P. (2021). Ecology and evolution of antimicrobial resistance in bacterial communities. ISME J. 15, 939–948. doi: 10.1038/s41396-020-00832-7
Bowles, T. M., Jilling, A., Morán-Rivera, K., Schnecker, J., and Grandy, A. S. (2022). Crop rotational complexity affects plant-soil nitrogen cycling during water deficit. Soil Biol. Biochem. 166:108552. doi: 10.1016/j.soilbio.2022.108552
Bowles, T. M., Mooshammer, M., Socolar, Y., Calderón, F., Cavigelli, M. A., Culman, S. W., et al. (2020). Long-term evidence shows that crop-rotation diversification increases agricultural resilience to adverse growing conditions in North America. One Earth 2, 284–293. doi: 10.1016/j.oneear.2020.02.007
Cai, Y., and Hao, M. D. (2015). Effects of rotation model and period on wheat yield, nutrient uptake and soil fertility in the Loess Plateau. J. Plant Nutr. Fertil. 21, 864–872. doi: 10.11674/zwyf.2015.0405
Cerqueira, F., Matamoros, V., Bayona, J. M., Berendonk, T. U., Elsinga, G., Hornstra, L. M., et al. (2019). Antibiotic resistance gene distribution in agricultural fields and crops: a soil-to-food analysis. Environ. Res. 177:108608. doi: 10.1016/j.envres.2019.108608
Chen, R. R., Senbayram, M., Blagodatsky, S., Myachina, O., Dittert, K., Lin, X. G., et al. (2014). Soil C and N availability determine the priming effect: microbial N mining and stoichiometric decomposition theories. Glob. Change Biol. 20, 2356–2367. doi: 10.1111/gcb.12475
Crotty, F. V., Fychan, R., Sanderson, R., Rhymes, J. R., Bourdin, F., Scullion, J., et al. (2016). Understanding the legacy effect of previous forage crop and tillage management on soil biology, after conversion to an arable crop rotation. Soil Biol. Biochem. 103, 241–252. doi: 10.1016/j.soilbio.2016.08.018
Darby, E. M., Trampari, E., Siasat, P., Gaya, M. S., Alav, I., Webber, M. A., et al. (2023). Molecular mechanisms of antibiotic resistance revisited. Nat. Rev. Microbiol. 21, 280–295. doi: 10.1038/s41579-022-00820-y
Ding, C. L., and Sun, J. (2025). The potential contribution of microbial communities to carbon fixation and nitrogen cycle in the eastern Indian Ocean. Mar. Environ. Res. 207:107056. doi: 10.1016/j.marenvres.2025.107056
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Fan, Y. Q., Wang, Y. N., Huo, R. X., Yao, T., Yang, Z. P., Qiao, Y. J., et al. (2022). Research status of interaction between soil microorganisms and nematode communities under rotation mode. Chin. Agric. Sci. Bull. 38, 108–113. doi: 10.11924/j.issn.1000-6850.casb2021-0859
Forsberg, K. J., Patel, S., Gibson, M. K., Lauber, C. L., Knight, R., Fierer, N., et al. (2014). Bacterial phylogeny structures soil resistomes across habitats. Nature 509, 612–616. doi: 10.1038/nature13377
Fu, X., Wang, J., Sainju, U. M., and Liu, W. Z. (2019). Soil nitrogen fractions under long-term crop rotations in the Loess Plateau of China. Soil Tillage Res. 186, 42–51. doi: 10.1016/j.still.2018.10.004
Gao, Y. J., Luo, W., Zhang, H. F., Chen, Y. Y., Li, Z. B., Wei, G. H., et al. (2023). Enrichment of antibiotic resistance genes in roots is related to specific bacterial hosts and soil properties in two soil–plant systems. Sci. Total Environ. 886:163933. doi: 10.1016/j.scitotenv.2023.163933
Gaudel, G., Xing, L., Shrestha, S., Poudel, M., Sherpa, P., Raseduzzaman, M., et al. (2024). Microbial mechanisms regulate soil organic carbon mineralization under carbon with varying levels of nitrogen addition in the above-treeline ecosystem. Sci. Total Environ. 917:170497. doi: 10.1016/j.scitotenv.2024.170497
Giacometti, C., Mazzon, M., Cavani, L., Triberti, L., Baldoni, G., Ciavatta, C., et al. (2021). Rotation and fertilization effects on soil quality and yields in a long term field experiment. Agronomy 11:636. doi: 10.3390/agronomy11040636
Gu, Z. G., Gu, L., Eils, R., Schlesner, M., and Brors, B. (2014). Circlize implements and enhances circular visualization in R. Bioinformatics 30, 2811–2812. doi: 10.1093/bioinformatics/btu393
Guan, X. J., Li, Y. H., Yang, Y. Y., Liu, Z. H., Shi, R. G., and Xu, Y. (2025). Root exudates regulate soil antibiotic resistance genes via rhizosphere microbes under long-term fertilization. Environ. Int. 195:109180. doi: 10.1016/j.envint.2024.109180
Guo, T., Lou, C. L., Zhai, W. W., Tang, X. J., Hashmi, M. Z., Murtaza, R., et al. (2018). Increased occurrence of heavy metals, antibiotics and resistance genes in surface soil after long-term application of manure. Sci. Total Environ. 635, 995–1003. doi: 10.1016/j.scitotenv.2018.04.194
Guo, X. P., Zhao, S., Chen, Y. R., Yang, J., Hou, L. J., Liu, M., et al. (2020). Antibiotic resistance genes in sediments of the Yangtze estuary: from 2007 to 2019. Sci. Total Environ. 744:140713. doi: 10.1016/j.scitotenv.2020.140713
Haichar, F. Z., Marol, C., Berge, O., Rangel-Castro, J. I., Prosser, J. I., Balesdent, J., et al. (2008). Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2, 1221–1230. doi: 10.1038/ismej.2008.80
Han, B. H., Ma, L., Yu, Q. L., Yang, J. W., Su, W. H., Hilal, M. G., et al. (2022). The source, fate and prospect of antibiotic resistance genes in soil: a review. Front. Microbiol. 13:976657. doi: 10.3389/fmicb.2022.976657
Han, B. J., Yang, F. X., Shen, S. Z., Mu, M. R., and Zhang, K. Q. (2023). Effects of soil habitat changes on antibiotic resistance genes and related microbiomes in paddy fields. Sci. Total Environ. 895:165109. doi: 10.1016/j.scitotenv.2023.165109
Haney, R. L., Franzluebbers, A. J., Porter, E. B., Hons, F. M., and Zuberer, D. A. (2004). Soil carbon and nitrogen mineralization. Soil Sci. Soc. Am. J. 68, 489–492. doi: 10.2136/sssaj2004.4890
Hothorn, T., Bretz, F., and Westfall, P. (2008). Simultaneous inference in general parametric models. Biom. J. 50, 346–363. doi: 10.1002/bimj.200810425
Hu, K. B., Zhao, P., Wu, K. X., Yang, H. L., Yang, Q. X., Fan, M. P., et al. (2023). Reduced and deep application of controlled-release urea maintained yield and improved nitrogen-use efficiency. Field Crop Res. 295:108876. doi: 10.1016/j.fcr.2023.108876
Huang, C. H., Zhou, S. Y. D., Li, H., Su, J. Q., and Wang, J. N. (2021). Profiling of antibiotic resistance genes in different croplands. Environ. Sci. 42, 2975–2980. doi: 10.13227/j.hjkx.202009091
Islam, K. R., Dick, W. A., Watts, D. B., Gonzalez, J. M., Fausey, N. R., Flanagan, D. C., et al. (2022). Gypsum, crop rotation, and cover crop impacts on soil organic carbon and biological dynamics in rainfed transitional no-till corn-soybean systems. PLoS One 17:e0275198. doi: 10.1371/journal.pone.0275198
Jechalke, S., Heuer, H., Siemens, J., Amelung, W., and Smalla, K. (2014). Fate and effects of veterinary antibiotics in soil. Trends Microbiol. 22, 536–545. doi: 10.1016/j.tim.2014.05.005
Jiang, Y. L., Xiao, L., Liu, J. H., Chen, Y., Deng, X. P., Duan, P. P., et al. (2025). Responses of soil nitrogen and carbon mineralization rates to fertilization and crop rotation. J. Soils Sediments 24, 1289–1301. doi: 10.1007/s11368-023-03694-6
Khalid, M., Liu, X. X., Zheng, B. X., Su, L. T., Kotze, D. J., Setälä, H., et al. (2023). Distinct climatic regions drive antibiotic resistance genes dynamics across public parks and pristine soil ecosystems. J. Clean. Prod. 409:137275. doi: 10.1016/j.jclepro.2023.137275
Lê, S., Josse, J., and Husson, F. (2008). Factominer: an R package for multivariate analysis. J. Stat. Soft. 25, 1–18. doi: 10.18637/jss.v025.i01
Li, L. G., Huang, Q., Yin, X. L., and Zhang, T. (2020). Source tracking of antibiotic resistance genes in the environment—challenges, progress, and prospects. Water Res. 185:116127. doi: 10.1016/j.watres.2020.116127
Li, Z. X., Sun, A. X., Liu, X., Chen, Q. L., Bi, L., Ren, P. X., et al. (2022). Climate warming increases the proportions of specific antibiotic resistance genes in natural soil ecosystems. J. Hazard. Mater. 430:128442. doi: 10.1016/j.jhazmat.2022.128442
Li, H. Y., Yang, Q. F., Liu, W., Li, R. L., Zhang, D. D., Zhang, G. L., et al. (2025). Plant secretions and volatiles contribute to the evolution of bacterial antibiotic resistance in soil-crop system. J. Environ. Sci. 152, 516–526. doi: 10.1016/j.jes.2024.05.039
Liu, B., Ahnemann, H., Arlotti, D., Huyghebaert, B., Cuperus, F., and Tebbe, C. C. (2024). Impact of diversified cropping systems and fertilization strategies on soil microbial abundance and functional potentials for nitrogen cycling. Sci. Total Environ. 932:172954. doi: 10.1016/j.scitotenv.2024.172954
Liu, Y., Liu, R., Ghimire, R., Zhang, N. N., Zhou, S., Zhao, F. Z., et al. (2023). Linking soil phosphorus fractions to associated microbial functional profiles under crop rotation on the Loess Plateau of China. Soil Tillage Res. 233:105809. doi: 10.1016/j.still.2023.105809
Liu, Y., Zhang, L. P., Lu, J. K., Chen, W. M., Wei, G. H., and Lin, Y. B. (2020). Topography affects the soil conditions and bacterial communities along a restoration gradient on loess-plateau. Appl. Soil Ecol. 150:103471. doi: 10.1016/j.apsoil.2019.103471
Liu, Z. S., Zhao, Y. X., Zhang, B. F., Wang, J. Q., Zhu, L. Z., and Hu, B. L. (2023). Deterministic effect of pH on shaping soil resistome revealed by metagenomic analysis. Environ. Sci. Technol. 57, 985–996. doi: 10.1021/acs.est.2c06684
Lv, M., Huang, M., Zhao, K. N., Gu, X. X., Li, S. Q., Wang, J. T., et al. (2023). Effects of partial substitution of rrganic fertilizer for synthetic N fertilizer on yield and N use efficiencies in a semiarid winter wheat–summer maize rotation. Agronomy 13:2281. doi: 10.3390/agronomy13092281
Marić, J., and Šikić, M. (2019). “Approaches to metagenomic classification and assembly” in 2019 42nd International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO), 348–356.
McDaniel, M. D., and Grandy, A. S. (2016). Soil microbial biomass and function are altered by 12 years of crop rotation. Soil 2, 583–599. doi: 10.5194/soil-2-583-2016
Mishra, A. K., Sudalaimuthuasari, N., Hazzouri, K. M., Saeed, E. E., Shah, I., and Amiri, K. M. A. (2022). Tapping into plant–microbiome interactions through the lens of multi-omics techniques. Cells 11:3254. doi: 10.3390/cells11203254
Ondon, B. S., Li, S. G., Zhou, Q. X., and Li, F. X. (2021). Sources of antibiotic resistant Bacteria (ARB) and antibiotic resistance genes (ARGs) in the soil: a review of the spreading mechanism and human health risks. Rev. Environ. Contam. Toxicol. 256, 121–153. doi: 10.1007/398_2020_60
Pruden, A., Larsson, D. J., Amézquita, A., Collignon, P., Brandt, K. K., Graham, D. W., et al. (2013). Management options for reducing the release of antibiotics and antibiotic resistance genes to the environment. Environ. Health Perspect. 121, 878–885. doi: 10.1289/ehp.1206446
Shen, Y. P., Liu, Y. B., Du, Y. T., Wang, X., Guan, J. N., Jia, X. H., et al. (2023). Transfer of antibiotic resistance genes from soil to wheat: role of host bacteria, impact on seed-derived bacteria, and affecting factors. Sci. Total Environ. 905:167279. doi: 10.1016/j.scitotenv.2023.167279
Su, J. Q., Wei, B., Xu, C. Y., Qiao, M., and Zhu, Y. G. (2014). Functional metagenomic characterization of antibiotic resistance genes in agricultural soils from China. Environ. Int. 65, 9–15. doi: 10.1016/j.envint.2013.12.010
Sun, J. T., Zeng, Q. T., Tsang, D. C. W., Zhu, L. Z., and Li, X. D. (2017). Antibiotics in the agricultural soils from the Yangtze River Delta, China. Chemosphere 189, 301–308. doi: 10.1016/j.chemosphere.2017.09.040
Thev, E. A. (2015). Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J. Proteome Res. 14, 3322–3335. doi: 10.1021/acs.jproteome.5b00354
Venter, Z. S., Jacobs, K., and Hawkins, H. J. (2016). The impact of crop rotation on soil microbial diversity: a meta-analysis. Pedobiologia 59, 215–223. doi: 10.1016/j.pedobi.2016.04.001
Wang, J., Fu, X., Ghimire, R., Sainju, U. M., Jia, Y., and Zhao, F. Z. (2021). Responses of soil bacterial community and enzyme activity to organic matter components under long-term fertilization on the Loess Plateau of China. Appl. Soil Ecol. 166:103992. doi: 10.1016/j.apsoil.2021.103992
Wang, L. J., Wang, J. H., Wang, J., Zhu, L. S., Conkle, J. L., and Yang, R. (2020). Soil types influence the characteristic of antibiotic resistance genes in greenhouse soil with long-term manure application. J. Hazard. Mater. 392:122334. doi: 10.1016/j.jhazmat.2020.122334
Wang, F., Xu, M., Stedtfeld, R. D., Sheng, H., Fan, J. B., Liu, M., et al. (2018). Long-term effect of different fertilization and cropping systems on the soil antibiotic resistome. Environ. Sci. Technol. 52, 13037–13046. doi: 10.1021/acs.est.8b04330
Wang, Z. R., Zhang, N., Li, C. M., and Shao, L. (2023). Diversity of antibiotic resistance genes in soils with four different fertilization treatments. Front. Microbiol. 14:1291599. doi: 10.3389/fmicb.2023.1291599
Wang, J. R., Zhou, J. J., and Zhu, G. F. (2024). A review of soil carbon, nitrogen, phosphorus and stoichiometry studies in China. Chin. J. Ecol. 43, 2493–2501. doi: 10.13292/j.1000-4890.202408.011
Wei, T. Y., and Simko, V. (2024). Corrplot: visualization of a correlation matrix. doi: 10.32614/CRAN.package.corrplot
Wen, Z. H., Li, H. G., Shen, J. B., and Rengel, Z. (2017). Maize responds to low shoot P concentration by altering root morphology rather than increasing root exudation. Plant Soil 416, 377–389. doi: 10.1007/s11104-017-3214-0
Wen, T., Niu, G. Q., Chen, T., Shen, Q. R., Yuan, J., and Liu, Y. X. (2023). The best practice for microbiome analysis using R. Protein Cell 14, 713–725. doi: 10.1093/procel/pwad024
Wu, J., Joergensen, R. G., Pommerening, B., Chaussod, R., and Brookes, P. C. (1990). Measurement of soil microbial biomass C by fumigation-extraction-an automated procedure. Soil Biol. Biochem. 22, 1167–1169. doi: 10.1016/0038-0717(90)90046-3
Wu, J., Wang, J. Y., Li, Z. T., Guo, S. M., Li, K. J., Xu, P. S., et al. (2023). Antibiotics and antibiotic resistance genes in agricultural soils: a systematic analysis. Crit. Rev. Environ. Sci. Technol. 53, 847–864. doi: 10.1080/10643389.2022.2094693
Xiao, R. H., Huang, D. L., Du, L., Song, B., Yin, L. S., Chen, Y., et al. (2023). Antibiotic resistance in soil-plant systems: a review of the source, dissemination, influence factors, and potential exposure risks. Sci. Total Environ. 869:161855. doi: 10.1016/j.scitotenv.2023.161855
Xie, W. Y., McGrath, S. P., Su, J. Q., Hirsch, P. R., Clark, I. M., Shen, Q., et al. (2016). Long-term impact of field applications of sewage sludge on soil antibiotic resistome. Environ. Sci. Technol. 50, 12602–12611. doi: 10.1021/acs.est.6b02138
Xu, F., Liu, Y. L., Du, W. C., Li, C. L., Xu, M. L., Xie, T. C., et al. (2021). Response of soil bacterial communities, antibiotic residuals, and crop yields to organic fertilizer substitution in North China under wheat–maize rotation. Sci. Total Environ. 785:147248. doi: 10.1016/j.scitotenv.2021.147248
Zack, K. M., Sorenson, T., and Joshi, S. G. (2024). Types and mechanisms of efflux pump systems and the potential of efflux pump inhibitors in the restoration of antimicrobial susceptibility, with a special reference to Acinetobacter baumannii. Pathogens 13:197. doi: 10.3390/pathogens13030197
Zhang, L. W. (2024). Temporal variation and risk assessment of typical fluoroquinolone antibiotic residues in a river-reservoir system during the reservoir discharge and filling periods. Water Air Soil Pollut. 235:573. doi: 10.1007/s11270-024-07406-9
Zhang, D. D., Li, H. Y., Yang, Q. F., and Xu, Y. (2024). Microbial-mediated conversion of soil organic carbon co-regulates the evolution of antibiotic resistance. J. Hazard. Mater. 471:134404. doi: 10.1016/j.jhazmat.2024.134404
Zhang, H. F., Luo, G. W., Wang, Y. Z., Fei, J. C., Rong, X. M., Peng, J. W., et al. (2023). Crop rotation-driven change in physicochemical properties regulates microbial diversity, dominant components, and community complexity in paddy soils. Agric. Ecosyst. Environ. 343:108278. doi: 10.1016/j.agee.2022.108278
Zhang, W. G., Wen, T., Liu, L. Z., Li, J. Y., Gao, Y., Zhu, D., et al. (2021). Agricultural land-use change and rotation system exert considerable influences on the soil antibiotic resistome in Lake Tai Basin. Sci. Total Environ. 771:144848. doi: 10.1016/j.scitotenv.2020.144848
Zhang, Y., Zhou, J., Wu, J., Hua, Q. W., and Bao, C. X. (2022). Distribution and transfer of antibiotic resistance genes in different soil–plant systems. Environ. Sci. Pollut. Res. 29, 59159–59172. doi: 10.1007/s11356-021-17465-8
Zheng, D. S., Yin, G. Y., Liu, M., Hou, L. J., Yang, Y., Boeckel, T. P. V., et al. (2022). Global biogeography and projection of soil antibiotic resistance genes. Sci. Adv. 8:eabq8015. doi: 10.1126/sciadv.abq8015
Keywords: antibiotic resistance genes, crop rotation complexity, crop rotation regimes, Loess Plateau, resistance mechanisms, winter wheat field
Citation: Hu R, Liu Y, Wen M, Zhou N and Wang J (2025) Crop rotation complexity affects soil properties shaping antibiotic resistance gene types and resistance mechanisms. Front. Microbiol. 16:1603518. doi: 10.3389/fmicb.2025.1603518
Edited by:
Jian Yang, China University of Geosciences Wuhan, ChinaReviewed by:
Liang Shen, Anhui Normal University, ChinaJialing Li, Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), China
Lianhai Ren, Beijing Technology and Business University, China
Copyright © 2025 Hu, Liu, Wen, Zhou and Wang. 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.
*Correspondence: Yang Liu, MjAyMjkwMzhAbnd1LmVkdS5jbg==; Jun Wang, d2FuZ2pAbnd1LmVkdS5jbg==