Jiawei Ren
Hongbin Lu
Shaoyong Lu1*- 1State Key Laboratory of Environmental Criteria and Risk Assessment, National Engineering Laboratory for Lake Pollution Control and Ecological Restoration, Chinese Research Academy of Environmental Sciences, Beijing, China
- 2Lake Forest Academy, Lake Forest, IL, United States
- 3State Key Laboratory of Coal Conversion, Institute of Coal Chemistry, Chinese Academy of Sciences, Taiyuan, China
Antibiotics are an important pharmaceutical class excessively used by humans. Its presence in the soil can impact plant growth and induce antibiotic resistance. This research studies the effect of sulfamethoxazole (SMX) on plant growth, rhizosphere bacteria composition, and resistance genes. Two sets of vegetables (basil, cilantro, and spinach) were treated separately with water and SMX solution. The plant growth data and soil samples were collected and analyzed. The results revealed that SMX increased spinach leaf length (34.0%) while having no significant impacts on basil and cilantro. On the other hand, SMX improved the bacterial diversity in all samples. The shifts in the abundance of plant growth-promoting bacteria could indirectly affect vegetable stem and leaf length. SMX also significantly increased the abundance of resistance genes Sul1 and Sul2. A further study into the correlation between bacteria highlights the importance of Shingomonas and Alfipia for inhibiting the spread of key resistance gene hosts, namely, Pseudomonas, Stenotrophomonas, and Agrobacterium. This research provides insight into SMX’s impact on vegetable growth and microbial diversity. It also points out important microbial interactions that could potentially be utilized to mitigate ARG proliferation.
1 Introduction
Antibiotics have been widely used for medication and animal husbandry (Patangia et al., 2022). They are the most important pharmaceutical class today and their consumption has increased by 46% between 2000 and 2018 (Browne et al., 2021; Uddin et al., 2021). However, antibiotics are not completely metabolized by humans and animals. As much as 30%–90% of antibiotics are excreted through urine and feces (Zhang et al., 2015). They end up in the environment through runoffs, manure treatments, and effluents of wastewater plants (Zhao et al., 2020).
Sulfamethoxazole (SMX) belongs to a group of antibiotics called sulfonamide (SA) (Coleman, 2022). Because of its affordability and versatility, SA has been commonly applied and is the second most concentrated antibiotic in the soil (approximately 25 ppm) (Robles-Jimenez et al., 2021). SA can affect plant growth and ultimately threaten public health through the food chain (Baran et al., 2011). SMX in particular has caused antibiotic resistance in the past decade worldwide (Du et al., 2019; Wang et al., 2019; Chia-Wei et al., 2022). Antibiotic-resistant genes (ARGs) are the main reasons that bacteria develop antibiotic resistance (Zhao et al., 2020). Antibiotic presence in the soil induces ARGs by exerting selective pressures. Bacteria develop resistance using mechanisms that tolerate, extrude, interfere, inactivate, or destruct the antibiotics. Those modifications are facilitated by mutational resistance and passed between bacterial species through horizontal gene transfer (Skandalis et al., 2021). ARB can be transmitted from farms to the human body via the food chain and environment, which causes prolonged illness, increased hospitalization, and treatment failure (World Health Organization). Even if ARB is killed or damaged, ARGs can be released into the water, soil, or air and transferred into other bacteria (Zhang et al., 2020; Zhuang et al., 2021). The World Health Organization considers antibiotic resistance to be one of the most urgent public health concerns (World Health Organization, 2023). According to a United Nations report, in 2019 an estimated 1.27 million global deaths were directly caused by drug-resistant disease (United Nations Environmental Programme, 2023).
Because of the limited access to clean water, wastewater is being increasingly used and recycled in crop irrigation worldwide (Slobodiuk et al., 2021). However, due to insufficient treatment, this practice exposes crops to antibiotics in the wastewater and affects plant growth (Kampouris et al., 2021). Pan and Chu’s study in Hong Kong found that the presence of sulfamethazine (157 mg/L) can significantly inhibit lettuce root and stem elongation (Pan and Chu, 2016). Several studies have shown that high antibiotic concentrations can lead to delayed germination, reduced root and stem length, decreased biomass, and even death among crops. Most of the research used starchy crops (maze, wheat, pea, etc.) or fruity crops such as cucumbers.
Not many studies have been done with vegetables such as spinach, cilantro, and basil (Migliore et al., 1995; Liu et al., 2009; Yang et al., 2010; Michelini et al., 2012). Therefore, investigating the effect of common antibiotics on vegetables could be valuable to modern agricultural development. Considering the connection between soil antibiotics, plants, and public health, it is crucial to understand how antibiotics impact vegetable growth, soil bacteria composition, and ARGs. To simulate the effects of SMX in wastewater irrigation, in this research three vegetables (spinach, cilantro, and basil) were treated with SMX solution. The plant growth analysis revealed the vegetable’s physiological response to SMX. The metagenomic analysis explored changes under SMX stress on rhizosphere bacteria and ARGs. These results provide a reference for SMX’s impact on farmland through wastewater irrigation and also offer suggestions to mitigate ARGs.
2 Materials and methods
2.1 Site setup
Two sets of ten pots were installed. One set was treated with sulfamethoxazole solution (75 mg/L) and the other with water. Within each set, three pots were assigned to each vegetable (spinach as Spinacia oleracea, cilantro as Coriandrum sativum, and basil as Ocimum basilicum), and one pot was designated for soil as control. All pots (22 cm in diameter, 21 cm in height) were filled with a layer of 18 cm organic soil as the substrate for plant growth. SMX solution (75 mg/L) was prepared using 98% SMX from Shanghai Aladdin Biochemical Technology Co., Ltd. The experiment was conducted under a temperature range of 27°C–33°C. Two sets of pots were treated with respective solutions (avg. 60mL/day) to keep them moist in an environment with a relative humidity of 57%.
The spinach seeds were directly dispersed into the pots. But basil and cilantro seeds had to break dormancy for all plants to sprout at a relatively similar time. Each cilantro seed was a pod of two or three individual seeds. The pods were broken gently using the handle of a metal spoon and soaked in 38°C water for 24 h. The basil seeds were placed on a sheet of moist paper towels for 24 h. After the pre-germination process, 12 seeds were dispersed evenly in each pot, different by plant type. The seeds were covered with an additional 3 mm layer of soil. Shortened names are assigned to each of the eight experimental groups (Table 1).
TABLE 1. Group name with corresponding treatment and vegetable.
2.2 Growth data and soil sample collection
After 19 days, four plants of medium height in each group were selected to track height and leaf growth. Heights were measured from the surface of the soil to the highest point of the plant. An original soil sample was collected. Three weeks after sowing, soil samples were collected from each group by either gathering soil of the plant rhizosphere or if no plant was present they were directly taken from the pot. All samples were kept at −80°C before being processed and analyzed at Wekemo Tech Co., Ltd. Shenzhen, China.
2.3 DNA extraction
DNA was extracted from the soil samples using the CTAB method. First, the soil samples were transferred into 2.0-mL Eppendorf tubes containing 1,000 µL of CTAB extraction buffer. The tubes were then bathed at 65°C and gently mixed until adequate homogenization. After centrifugation, the supernatants were mixed with an equal volume of phenol chloroform isoamyl alcohol (25:24:1) which were centrifuged at 12,000 rpm for 10 min. Chloroform isoamyl alcohol (24:1) was added to supernatants, followed by mixing and 10 min centrifugation at 12,000 rpm. Next, the supernatants were transferred into 1.5-mL tubes with isopropyl alcohol, mixed by inversion, and incubated at −20°C. After the precipitation and solution separated, the tubes were centrifuged at 12,000 rpm for 10 min and the precipitations were washed with 1 mL 75% ethanol twice. After drying the precipitation, ddH2O and 1ul of RNase were added. The mixtures were incubated at 37°C for 15 min.
2.4 Metagenomic analysis
2.4.1 DNA library constructing, sequencing, and preprocessing
NEBNext® UltraTM DNA Library Prep Kit for Illumina (NEB, USA) was used for DNA extraction to construct sequencing libraries, following the manufacturer’s protocols. The qualified library was sequenced on the Illumina NovaSeq PE150 platform at Wekemo Tech Co., Ltd. Shenzhen, China.
Raw sequenced libraries were processed at Wekmo Tech Co. using Kneaddata software. First, sequencing adaptors, amino acid bases with quality scores less than 20, and sequences shorter than 50 bp were trimmed based on the Trimmomatic algorithm (Bolger et al., 2014). To reduce host pollution, genes of soil bacteria origin were deleted by aligning the sequenced data with bacterial genes on Bowtie2 (Langmead and Salzberg, 2012). Finally, clean data were obtained.
2.4.2 Taxonomy annotation
To understand the taxonomic distribution of the DNA library, kraken2 and a self-built database (Wekemo Tech Co.) were used to classify and annotate genes in the sample. Kraken2 that stores data of k-mer and its least common ancestor (LCA) matches the gene sequence with the most specific node of the taxonomic tree (Wood and Salzberg, 2014). After identifying the taxonomic distribution, Bracken was used to predict the actual relative abundance of species in the samples. By redistributing reads to species level through probability calculation, this method generated a species estimation of the dataset (Lu et al., 2017).
2.4.3 Functional gene analysis
The clean sequences were compared with the protein database UniRed90 using HUMAnN3 software (Franzosa et al., 2018). After filtering out failed reads, the relative abundance of the functional genes was calculated using the TPM metric, which counts all matching reads in a fragment as one fragment count. This method calculated the proportion over an arbitrary total count, which is useful for comparing different samples (Puente-Sánchez et al., 2020; Zhao et al., 2021). Annotated genes were then sent to the ARDB database for ARG identification and to the KEGG database for host organism annotation (Liu and Pop, 2009; Kanehisa et al., 2017).
2.4.4 Statistical analysis
Basic data analysis was performed on Microsoft Excel 2016 to obtain averages and standard deviations. One way ANOVA test carried out on Wekemo (Wekemo Tech Co. was used to analyze the statistical significance of changes caused by sulfamethoxazole. Any differences with a p-value lower than 0.05 were considered significant.
3 Results and discussions
3.1 Antibiotic impact on plant growth
Vegetable growth under SMX stress was accessed through two parameters visualized in Figure 1: stem height (graph a-c) and leaf length (graph d-f). At the start of the measuring period, all plants treated with SMX solution had similar lengths compared to those under normal conditions. A linear regression analysis showed that SMX solution improved basil growth on average by 0.05 cm/day and its final medium height (10.5 cm) surpassed the height under control conditions (9.3 cm) by 12.9%. However, SMX only changed the growth of cilantro and spinach by about 0.01 cm/day. A further significance analysis corresponded to the linear regression, revealing that 75 mg/L SMX solution did not pose significant impacts on basil, cilantro, and spinach growth (p-value >0.05 for all three plants).
FIGURE 1. (A–C): Stem heights of basil, cilantro and spinach. (D–F): Leaf lengths of basil, cilantro, and spinach at the start, middle, and end of the measuring period.
Compared to the effect on stem growth, SMX had a clearer impact on leaf growth in all vegetables. Graph d-f in Figure 1 shows the leaf lengths at the start, middle, and end of the measuring period. Opposite to SMX’s facilitation in basil stem development, antibiotics led to smaller basil leaves throughout the measuring period. Its final medium length (3.92 cm) under SMX treatment was 13.8% shorter than the average leaf length (4.55 cm) of the control group. SMX also adversely affected the development of cilantro leaves with a final inhibition rate of 5.9%. Spinach leaf lengths were similar between the antibiotic group and the control group at the start and middle of the growth period. However, by the end of the period, Spinach leaves in the antibiotics group (12.13 cm) significantly outperformed the control group (9.05 cm) by 34.0% (p-value = 0.042 < 0.05).
Stem and leaf growth are indicators to assess the phytotoxicity of environmental pollutants (Mei et al., 2021). The data showed that SMX insignificantly promoted the basil stem growth. This result corresponded to Zhang et al.‘s study on pakchoi (B. chinensis) where SMX improved pakchoi height (Zhang et al., 2017). Interestingly, our findings showed the opposite of the general understanding that SMX adversely affects plant growth (Cheong et al., 2020). For instance, high SMX concentration in the soil could lead to a 44.7% decrease in tomato plant height and a 38.2% decrease in cucumber plant height. Similarly, Lv et al. found that exposure to 100 mg/L SMX solution resulted in a 28.13% decrease in strawberry (Fragaria ananassa) stem weight (Lv et al., 2021). This dichotomous effect of SMX between our results and the general knowledge could be explained by the plant’s hormetic response to antibiotics. Researchers found that while high concentrations of antibiotics tend to inhibit plant growth, low concentrations actually facilitate stem development (Migliore et al., 2000). This phenomenon was supported by recent studies. The low dosages of SMX led to longer stems in both pakchoi and strawberry, while higher concentrations reduced the stem length (Zhang et al., 2017; Lv et al., 2021). Therefore, it is possible that under SMX stress basil exhibited an insignificant height increase, while cilantro and spinach were not affected at all.
SMX impacted the leaf differently compared to its effect on plant height. SMX’s inhibition effect on basil and cilantro leaves could be due to the antibiotic’s higher concentration in leaves than stems. Although the sulfamethoxazole uptake mechanism is not fully understood, the majority of studies have pointed to passive diffusion as the main pathway (Franklin et al., 2016). Thus, the transpiration stream plays a crucial role in SMX’s translocation. Since stems are the main antibiotic channels rather than reservoirs, stems do not retain as many antibiotics as leaves (Liu et al., 2013). Several studies have shown that leaves accumulate more antibiotics than stems (Herklotz et al., 2010; Liu et al., 2013; Zhao et al., 2019). Therefore, basil and cilantro leaves could accumulate higher antibiotic concentrations than stems, which inhibits leaf development. On the other hand, spinach tends to accumulate fewer antibiotics in its tissues than other plants (Kodešová et al., 2019). Thus, SMX could exist in spinach leaves at a relatively lower concentration just enough to stimulate leaf growth.
The results showed that plants had different antibiotic tolerance. While SMX led to longer basil stems and larger spinach leaves, it suppressed cilantro’s growth. However, SMX’s impact on plant growth was overall insignificant, which suggested that vegetable stem and length lengths were not reliable indicators of SMX presence in the soil. Despite these insignificant impacts, SMX could still noticeably affect rhizobacteria diversity. Therefore, a further investigation into the microbial community under SMX stress is needed.
3.2 Analysis of microbial diversity
Raw sequence libraries were constructed based on the soil sample from eight pots. The number of raw reads ranged from 22884329 to 33657674 (avg. 26714901). Low-quality or short sequences were filtered out with the Trimmomatic algorithm, followed by the removal of host gene. As a result, 21521875 to 32085255 clean reads (avg. 25269515) were obtained for each group. A comparison with kraken2 and a self-built database (Wekemo Tech Co.) revealed the taxonomic distribution of each sample.
Alpha diversity indices, calculating each sample based on 5329 identified operational taxonomic units (OTUs), are shown in Table 2. Both the abundance-based coverage estimator (ace) and chao1 assessed the species richness. Compared to ace, chao1 gave more weight to the species of low abundance (Bo-Ra et al., 2017). Higher ace and chao1 values were observed for soil, basil, and cilantro treated with SMX. This result suggested that the presence of SMX in the sample increased the species richness except for spinach. Shannon and Simpson’s indices reflected bacterial diversity through richness and evenness. While Shannon emphasized more on bacterial species richness, Simpson focused on evenness (Suturina et al., 2022). As both Shannon and Simpson indices were generally higher for SMX groups than control, SMX in this experiment increased bacterial richness and evenness in all samples.
TABLE 2. Abundance indices of microbial community composition in each group.
The microbial composition at phylum, order, family, and genus levels are presented in Figure 2. SMX changed microbial composition in pots with the vegetables. At the phylum level, Pseudomonadota was identified as the dominant phylum for all samples, followed by Actinomycetota. Pseudomonadota, better known as proteobacteria, is the most prevalent bacterial phylum in the environment (Leão et al., 2023; Selmani et al., 2023). This phylum includes many rhizobacteria associated with plant growth promoting activities, including nitrogen fixation, phosphorus uptake, and antifungal effect (Ramakrishna et al., 2019; Gómez-Godínez et al., 2023). At order rank, Xanthomonadales, Burkholderiales, and Pseudomonadales are the top three most dominant orders all of which belong to the phylum Psedomonadota. SMX increased the richness of Xanthomonadales in BAA (16.4%) and CIA (16.8%), compared to Xanthomonadales proportion in BAT (9.2%) and CIT (10.3%). At the family level, Xanthomonadacae, under order Xanthomonadales, was the most dominant family. It existed in a larger abundance in BAA and CIA than BAT and CIT. Although SPA had a similar soil profile to other plants, SPT altered the microbial composition significantly. Xanthomonadacae took up 42.1% in SPT, compared to the average abundance of Xanthomonadacae in the rest of the samples (13.8%). Spinach also facilitated the proliferation of Rhizobiaceae (15.9%), a family under Hyphomicrobiales, which is a symbiotic nitrogen-fixing bacteria (Ormeño-Orrillo et al., 2015). Gómez et al. found that Rhizobium could colonize spinach roots and promote root, shoot, and leaf growth (Jimenez Gómez et al., 2016).
FIGURE 2. Relative bacterial abundance of eight groups at phylum, order, family, and genus level. (A) Relative bacterial abundances of the eight groups at the phylum level. (B) Relative bacterial abundances of the eight groups at the order level. (C) Relative bacterial abundances of the eight groups at the family level. (D) Relative bacterial abundances of the eight groups at the genus level.
At the genus level, the effect of antibiotics was registered across all plant types. Lysobacter under the family Xanthomonadacae was the dominant genus. Its abundance in basil, cilantro, and spinach with SMX (14.9%, 14.8%, 9.0% respectively) was consistently higher than in samples treated with water (7.2%, 8.8%, 7.7%). Lysobacter protects plants through the release of antimicrobial compounds and facilitates the degradation of macromolecules with extracellular enzymes (Gómez Expósito et al., 2015; Brescia et al., 2020). Its high concentration in sulfamethoxazole-contaminated soil could be explained by plants’ need for antibiotic degradation. Pseudomonas and Hydrogenophaga were the second and third most abundant genus. While Pseudomonas promotes plant growth and phytopathogen control, Hydrogenophaga modifies pollutants through enzymic activities and transforms them into less toxic substances (Sah et al., 2021; Majeed et al., 2022). SMX did not have noticeable impacts on the abundance of Pseudomonas and Hydrogenophaga. Three genera emerged from SPT, including Pseudoxanthomonas, Sphingomonas, and Neorhizobium. Peseudoxanthomonas and Neorhizobium facilitate plant growth through molecule degradation and nitrogen-fixing respectively. Sphingomonas increase plant’s stress tolerance and production of phytohormone (Sah et al., 2021; Lombardino et al., 2022; Marzba et al., 2022).
From the heat map in Figure 3, SPA had high abundances of Sphigobium, Hydrogenphaga, Acidvorax, Stutzerimonas, and Bordetella, which could induce spinach’s larger leaf length. On the other hand, SMX led to decreases in the abundance of Streptomyces and Sulfuricaulis in the basil and cilantro rhizosphere. Streptomyces aids the plant with nitrogen, phosphorus, and potassium uptake while serving as a biocontrol agent against plant pathogens (Amaresan et al., 2018). Sulfuricaulis plays the role of a sulfur-oxidizer. The decrease in Streptomyces and Sulfuricaulis under SMX stress might eventually lead to smaller basil and cilantro leaves.
FIGURE 3. Heatmap comparing bacterial relative abundance across eight groups at the genus level.
Overall, SMX impacted the bacterial composition and increased the diversity in all samples of this experiment. Past studies have shown that antibiotics altered the bacteria diversity through the elimination of sensitive bacteria and the outgrowth of resistant bacteria (Cycoń et al., 2019). Corresponding to the study by Reichel et al., Pseudomonas were insignificantly impacted by SMX, possibly because of the prevalence of resistance genes in Pseudomonas (Reichel et al., 2014). SMX’s impact on bacterial composition could contribute to the changes in plant growth. For instance, a decrease in the abundance of Streptomyces and Sulfuricaulis could have led to smaller basil and cilantro leaves, while the increase in Sphigobium, Hydrogenphaga, Acidvorax, Stutzerimonas, Bordetella might facilitate spinach leaf growth. Although the analysis of microbial diversity demonstrated SMX’s impact on soil bacteria, an examination of ARGs was still needed to study the correlation between antibiotics and ARB.
3.3 Functional gene host identification and interbacterial relationships
The cleaned sequence library was compared with the protein database UniRef 90. ARGs in each sample were identified through the ARBD database. According to Figure 4, MacB was the most abundant gene in every rhizosphere sample (average 19.4%). This gene encodes for an ATP-binding cassette (ABC) transporter that exports antibiotics specifically macrolide through mechanotransmission (Greene et al., 2018). Among the top 10 genes, 7 genes encoded for antibiotics efflux pump (MacB, BcrA, MexF, CeoB, MexW, RosB, AcrB), and 3 genes encoded for antibiotics target alternation (VanRA, BacA, Pbp1a) (Charlebois et al., 2012; Stogios and Savchenko, 2020; Fauzia et al., 2023). The presence of SMX spiked the amount of Sul1 on average by 108% and increased the amount of Sul2 by 33.7% (p-value = 0.006). Both Sul1 and Sul2 protect the bacteria through antibiotic target alternation by encoding for dihydropteroate synthase (DHPS) with low affinity to sulfamethoxazole. This discovery supported the fact that antibiotics could establish selective pressure that led to the proliferation of specific ARGs (Skandalis et al., 2021).
FIGURE 4. Abundance of ARGs in each sample.
KEGG orthology database identified host bacteria and annotated functional genes similar to known resistance genes (Martínez et al., 2015). Figure 5 reveals the concentration of functional genes among bacterial genera. Due to its greatest abundance, Lysobacter hosted most functional genes. Among the top 200 species, genus Pseudomonas held 13,778 reads on average, while Hydrogenophaga contained 84,529 reads. Pseudomonas’s abundance was only 3.6% higher than Hydrogenophaga’s, but Pseudomonas held 62.4% more functional genes. This finding revealed that Pseudomonas tend to host more functional genes than Hydrogenphaga. Correspondent to our discovery, Ashy et al. also concluded that Pseudomonas was a key host for viral genes and antibiotic resistance genes in the rhizosphere of moringa (Moringa oleifera) (Ashy et al., 2023). In the soil sample of spinach treated with water, Stenotrophomonas and Agrobacterium were not abundant genera, but hosted relatively large amounts of functional genes. The presence of those functional gene hosts, namely, Pseudomonas, Stenotrophomonas, and Agrobacterium, could facilitate the proliferation of ARGs. Figure 6 visualizes the interactions among bacterial genera. Shingomonas and Alfipia are two common bacteria that did not carry many functional genes but hindered other key hosts. Shingomonas competed with Pseudomonas, while Alfipia limited the growth of Stenotrophomonas and Agrobacterium. Therefore, it would be potentially beneficial to monitor and facilitate the growth of Shingomonas and Alfipia to outcompete functional gene hosts.
FIGURE 5. Heatmap of key functional gene hosts in each sample.
FIGURE 6. Correlation between bacteria genera. The triangle symbolizes phylum Pseudomonadota, circle symbolizes the phylum Actinomycetota. The green line for positive and the red line for negative correlation.
4 Conclusion
This study analyzed sulfamethoxazole’s impact on vegetable growth, soil bacteria composition, and ARG profile. In the experiment, SMX promoted the development of spinach leaves and basil stems, but had no significant impact on cilantro growth. A further study into the rhizosphere bacterial composition showed an increased microbial diversity under SMX stress. Most of the soil bacteria identified promoted plant growth through nutrient uptake and antimicrobial defense. SMX eliminated sensitive bacteria and increased the proportions of resistant species. These shifts in the bacterial community could contribute to the observed changes in plant growth. SMX also led to the proliferation of the Sul1 and Sul2 genes. Pseudomonas, Stenotrophomonas, and Agrobacterium were identified as key virulent gene hosts. By analyzing bacterial interactions, this study recommended promoting the growth of Shingomonas and Alfipia in the soil to outcompete the resistance gene hosts, which could reduce the amount of ARGs. However, future investigation is required to verify and consolidate the method.
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 a the following link: https://doi.org/10.6084/m9.figshare.25115306.v1.
Author contributions
JR: Formal Analysis, Investigation, Software, Writing–original draft. HL: Data curation, Funding acquisition, Methodology, Resources, Supervision, Writing–review and editing. SL: Funding acquisition, Writing–review and editing. ZH: Formal analysis, Writing–review and editing
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the Open Research Fund of State Environmental Protection Key Laboratory for Lake Pollution Control (2022HPYB-09), ICC-CAS(SCJC-WRW-2023-18), Fundamental Research Program of Shanxi Province (202203021222404).
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.
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.
References
Amaresan, N., Kumar, K., Naik, J. H., Bapatla, K. G., and Mishra, R. K. (2018). “Chapter 8 - Streptomyces in plant growth promotion: mechanisms and role,” in New and future developments in microbial Biotechnology and bioengineering. Editors B. P. Singh, V. K. Gupta, and A. K. Passari (Elsevier), 125–135.
Ashy, R. A., Jalal, R. S., Sonbol, H. S., Alqahtani, M. D., Sefrji, F. O., Alshareef, S. A., et al. (2023). Functional annotation of rhizospheric phageome of the wild plant species Moringa oleifera. Front. Microbiol. 14, 1166148. doi:10.3389/fmicb.2023.1166148
Baran, W., Adamek, E., Ziemiańska, J., and Sobczak, A. (2011). Effects of the presence of sulfonamides in the environment and their influence on human health. J. Hazard Mater 196, 1–15. doi:10.1016/j.jhazmat.2011.08.082
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30 (15), 2114–2120. doi:10.1093/bioinformatics/btu170
Bo-Ra, K., Jiwon, S., Robin, B. G., Lee, J. H., Kim, D. W., Seol, K. H., et al. (2017). Deciphering diversity indices for a better understanding of microbial communities. J. Microbiol. Biotechnol. 27 (12), 2089–2093. doi:10.4014/jmb.1709.09027
Brescia, F., Pertot, I., and Puopolo, G. (2020). “Chapter 16 - lysobacter,” in Beneficial microbes in agro-ecology. Editors N. Amaresan, M. Senthil Kumar, K. Annapurna, K. Kumar, and A. Sankaranarayanan (Academic Press), 313–338.
Browne, A. J., Chipeta, M. G., Haines-Woodhouse, G., Kumaran, E. P. A., Hamadani, B. H. K., Zaraa, S., et al. (2021). Global antibiotic consumption and usage in humans, 2000-18: a spatial modelling study. Lancet Planet Health 5 (12), e893–e904. doi:10.1016/s2542-5196(21)00280-1
Charlebois, A., Jalbert, L. A., Harel, J., Masson, L., and Archambault, M. (2012). Characterization of genes encoding for acquired bacitracin resistance in Clostridium perfringens. PLoS One 7 (9), e44449. doi:10.1371/journal.pone.0044449
Cheong, M. S., Seo, K. H., Chohra, H., Yoon, Y. E., Choe, H., Kantharaj, V., et al. (2020). Influence of sulfonamide contamination derived from veterinary antibiotics on plant growth and development. Antibiot. (Basel) 9 (8), 456. doi:10.3390/antibiotics9080456
Chia-Wei, L., Cheng, J.-F., Tung, K.-C., Hong, Y. K., Lin, J. H., Lin, Y. H., et al. (2022). Evolution of trimethoprim/sulfamethoxazole resistance in Shewanella algae from the perspective of comparative genomics and global phylogenic analysis. J. Microbiol. Immunol. Infect. 55 (6), 1195–1202. doi:10.1016/j.jmii.2021.09.014
Coleman, T. R. K. M. (2022). Trimethoprim sulfamethoxazole. St. Petersburg, FL, United States: StatPearls Publishing.
Cycoń, M., Mrozik, A., and Piotrowska-Seget, Z. (2019). Antibiotics in the soil environment—degradation and their impact on microbial activity and diversity. Front. Microbiol. 10, 338. doi:10.3389/fmicb.2019.00338
Du, B., Yang, Q., Wang, R., and Xin, Y. (2019). Evolution of antibiotic resistance and the relationship between the antibiotic resistance genes and microbial compositions under long-term exposure to tetracycline and sulfamethoxazole. Int. J. Environ. Res. Public Health 16 (23), 4681. doi:10.3390/ijerph16234681
Fauzia, K. A., Aftab, H., Tshibangu-Kabamba, E., Alfaray, R. I., Saruuljavkhlan, B., Cimuanga-Mukanya, A., et al. (2023). Mutations related to antibiotics resistance in Helicobacter pylori clinical isolates from Bangladesh. Antibiot. (Basel) 12 (2), 279. doi:10.3390/antibiotics12020279
Franklin, A. M., Williams, C. F., Andrews, D. M., Woodward, E. E., and Watson, J. E. (2016). Uptake of three antibiotics and an antiepileptic drug by wheat crops spray irrigated with wastewater treatment plant effluent. J. Environ. Qual. 45 (2), 546–554. doi:10.2134/jeq2015.05.0257
Franzosa, E. A., McIver, L. J., Rahnavard, G., Thompson, L. R., Schirmer, M., Weingart, G., et al. (2018). Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 15 (11), 962–968. doi:10.1038/s41592-018-0176-y
Gómez Expósito, R., Postma, J., Raaijmakers, J. M., and De Bruijn, I. (2015). Diversity and activity of lysobacter species from disease suppressive soils. Front. Microbiol. 6, 1243. doi:10.3389/fmicb.2015.01243
Gómez-Godínez, L. J., Aguirre-Noyola, J. L., Martínez-Romero, E., Arteaga-Garibay, R. I., Ireta-Moreno, J., and Ruvalcaba-Gómez, J. M. (2023). A look at plant-growth-promoting bacteria. Plants 12 (8), 1668. doi:10.3390/plants12081668
Greene, N. P., Kaplan, E., Crow, A., and Koronakis, V. (2018). Antibiotic resistance mediated by the MacB ABC transporter family: a structural and functional perspective. Front. Microbiol. 9, 950. doi:10.3389/fmicb.2018.00950
Herklotz, P. A., Gurung, P., Vanden Heuvel, B., and Kinney, C. A. (2010). Uptake of human pharmaceuticals by plants grown under hydroponic conditions. Chemosphere 78 (11), 1416–1421. doi:10.1016/j.chemosphere.2009.12.048
Jimenez Gómez, A., Menendez, E., Flores-Félix, J. D., Garcia-Fraile, P., Mateos, P., and Rivas, R. (2016). “Effective colonization of spinach root surface by Rhizobium,” in Biological nitrogen fixation and beneficial plant-microbe interaction. Editors E. James, and F. González-Andrés (Springer Cham), 109–122.
Kampouris, I. D., Klümper, U., Agrawal, S., Orschler, L., Cacace, D., Kunze, S., et al. (2021). Treated wastewater irrigation promotes the spread of antibiotic resistance into subsoil pore-water. Environ. Int. 146, 106190. doi:10.1016/j.envint.2020.106190
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 45 (D1), 353–361. doi:10.1093/nar/gkw1092
Kodešová, R., Klement, A., Golovko, O., Fér, M., Nikodem, A., Kočárek, M., et al. (2019). Root uptake of atenolol, sulfamethoxazole and carbamazepine, and their transformation in three soils and four plants. Environ. Sci. Pollut. Res. Int. 26 (10), 9876–9891. doi:10.1007/s11356-019-04333-9
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9 (4), 357–359. doi:10.1038/nmeth.1923
Leão, I., de Carvalho, T. B., Henriques, V., Ferreira, C., Sampaio-Maia, B., and Manaia, C. M. (2023). Pseudomonadota in the oral cavity: a glimpse into the environment-human nexus. Appl. Microbiol. Biotechnol. 107 (2-3), 517–534. doi:10.1007/s00253-022-12333-y
Liu, B., and Pop, M. (2009). ARDB--Antibiotic resistance genes database. Nucleic Acids Res. 37 (D), 443–447. doi:10.1093/nar/gkn656
Liu, F., Ying, G. G., Tao, R., Zhao, J. L., Yang, J. F., and Zhao, L. F. (2009). Effects of six selected antibiotics on plant growth and soil microbial and enzymatic activities. Environ. Pollut. 157 (5), 1636–1642. doi:10.1016/j.envpol.2008.12.021
Liu, L., Liu, Y.-h., Liu, C.-x., Wang, Z., Dong, J., Zhu, G. f., et al. (2013). Potential effect and accumulation of veterinary antibiotics in Phragmites australis under hydroponic conditions. Ecol. Eng. 53, 138–143. doi:10.1016/j.ecoleng.2012.12.033
Lombardino, J., Bijlani, S., Singh Nitin, K., Wood, J. M., Barker, R., Gilroy, S., et al. (2022). Genomic characterization of potential plant growth-promoting features of Sphingomonas strains isolated from the international space station. Microbiol. Spectr. 10 (1), 0199421–201921. doi:10.1128/spectrum.01994-21
Lu, J., Breitwieser, F., Thielen, P., and Salzberg, S. L. (2017). Bracken: estimating species abundance in metagenomics data. PeerJ Comput. Sci. 3, e104. doi:10.7717/peerj-cs.104
Lv, Y., Li, Y., Liu, X., and Xu, K. (2021). Effect of soil sulfamethoxazole on strawberry (Fragaria ananassa): growth, health risks and silicon mitigation. Environ. Pollut. 286, 117321. doi:10.1016/j.envpol.2021.117321
Majeed, A., Muhammad, Z., Ullah, R., Ullah, K., Ali, H., and Inayat, N. (2022). “Chapter 18 - plant growth-promoting rhizobacteria as bioremediators of polluted agricultural soils: challenges and prospects,” in Hazardous and trace materials in soil and plants. Editors M. Naeem, T. Aftab, A. Ali Ansari, S. S. Gill, and A. Macovei (Academic Press), 265–275.
Martínez, J. L., Coque, T. M., and Baquero, F. (2015). What is a resistance gene? Ranking risk in resistomes. Nat. Rev. Microbiol. 13 (2), 116–123. doi:10.1038/nrmicro3399
Marzban, M., Pourbabaee, A., Amoozegar, M., Naghavi, M. R., and Abbasi, A. (2022). Biodiversity and phylogenetic study of root nodule bacteria isolated from some wildly grown legume species growing in rangeland areas. Biharean Biol. 16, 27–37.
Mei, L., Zhu, Y., Zhang, X., Zhou, X., Zhong, Z., Li, H., et al. (2021). Mercury-induced phytotoxicity and responses in upland cotton (Gossypium hirsutum L.) seedlings. Plants (Basel) 10 (8), 1494. doi:10.3390/plants10081494
Michelini, L., Reichel, R., Werner, W., Ghisi, R., and Thiele-Bruhn, S. (2012). Sulfadiazine uptake and effects on Salix fragilis L. And Zea mays L. Plants. Water, air. Soil Pollut. 223 (8), 5243–5257. doi:10.1007/s11270-012-1275-5
Migliore, L., Brambilla, G., and Cozzolino, S. (1995). Effect on plants of sulphadimethoxine used in intensive farming (Panicum miliaceum, Pisum sativum and Zea mays). Agric. Ecosyst. Environ. 52 (2), 103–110. doi:10.1016/0167-8809(94)00549-t
Migliore, L., Cozzolino, S., and Fiori, M. (2000). Phytotoxicity to and uptake of flumequine used in intensive aquaculture on the aquatic weed, Lythrum salicaria L. Chemosphere 40 (7), 741–750. doi:10.1016/s0045-6535(99)00448-8
Ormeño-Orrillo, E., Servín-Garcidueñas, L. E., Rogel, M. A., González, V., Peralta, H., Mora, J., et al. (2015). Taxonomy of rhizobia and agrobacteria from the Rhizobiaceae family in light of genomics. Syst. Appl. Microbiol. 38 (4), 287–291. doi:10.1016/j.syapm.2014.12.002
Pan, M., and Chu, L. M. (2016). Phytotoxicity of veterinary antibiotics to seed germination and root elongation of crops. Ecotoxicol. Environ. Saf. 126, 228–237. doi:10.1016/j.ecoenv.2015.12.027
Patangia, D. V., Anthony Ryan, C., Dempsey, E., Paul Ross, R., and Stanton, C. (2022). Impact of antibiotics on the human microbiome and consequences for host health. Microbiologyopen 11 (1), e1260. doi:10.1002/mbo3.1260
Puente-Sánchez, F., García-García, N., and Tamames, J. (2020). SQMtools: automated processing and visual analysis of 'omics data with R and anvi'o. BMC Bioinforma. 21 (1), 358. doi:10.1186/s12859-020-03703-2
Ramakrishna, W., Yadav, R., and Li, K. (2019). Plant growth promoting bacteria in agriculture: two sides of a coin. Appl. Soil Ecol. 138, 10–18. doi:10.1016/j.apsoil.2019.02.019
Reichel, R., Radl, V., Rosendahl, I., Albert, A., Amelung, W., Schloter, M., et al. (2014). Soil microbial community responses to antibiotic-contaminated manure under different soil moisture regimes. Appl. Microbiol. Biotechnol. 98 (14), 6487–6495. doi:10.1007/s00253-014-5717-4
Robles-Jimenez, L. E., Aranda-Aguirre, E., Castelan-Ortega, O. A., Shettino-Bermudez, B. S., Ortiz-Salinas, R., Miranda, M., et al. (2021). Worldwide traceability of antibiotic residues from livestock in wastewater and soil: a systematic review. Anim. (Basel) 12 (1), 60. doi:10.3390/ani12010060
Sah, S., Krishnani, S., and Singh, R. (2021). Pseudomonas mediated nutritional and growth promotional activities for sustainable food security. Curr. Res. Microb. Sci. 2, 100084. doi:10.1016/j.crmicr.2021.100084
Selmani, Z., Attard, E., Lauga, B., Barakat, M., Ortet, P., Tulumello, J., et al. (2023). Culturing the desert microbiota. Front. Microbiol. 14, 1098150. doi:10.3389/fmicb.2023.1098150
Skandalis, N., Maeusli, M., Papafotis, D., Miller, S., Lee, B., Theologidis, I., et al. (2021). Environmental spread of antibiotic resistance. Antibiot. (Basel) 10 (6), 640. doi:10.3390/antibiotics10060640
Slobodiuk, S., Niven, C., Arthur, G., Thakur, S., and Ercumen, A. (2021). Does irrigation with treated and untreated wastewater increase antimicrobial resistance in soil and water: a systematic review. Int. J. Environ. Res. Public Health 18 (21), 11046. doi:10.3390/ijerph182111046
Stogios, P. J., and Savchenko, A. (2020). Molecular mechanisms of vancomycin resistance. Protein Sci. 29 (3), 654–669. doi:10.1002/pro.3819
Suturina, L., Belkova, N., Igumnov, I., Lazareva, L., Danusevich, I., Nadeliaeva, I., et al. (2022). Polycystic ovary syndrome and gut microbiota: phenotype matters. Life (Basel) 13 (1), 7. doi:10.3390/life13010007
Uddin, T. M., Chakraborty, A. J., Khusro, A., Zidan, B. R. M., Mitra, S., Emran, T. B., et al. (2021). Antibiotic resistance in microbes: history, mechanisms, therapeutic strategies and future prospects. J. Infect. Public Health 14 (12), 1750–1766. doi:10.1016/j.jiph.2021.10.020
United Nations Environmental Programme (2023). Bracing for Superbugs: strengthening environmental action in the One Health response to antimicrobial resistance. Nairobi, Kenya: United Nation Environmental Programme.
Wang, H.-T., Chi, Q.-Q., Zhu, D., Li, G., Ding, J., An, X. L., et al. (2019). Arsenic and sulfamethoxazole increase the incidence of antibiotic resistance genes in the gut of earthworm. Environ. Sci. Technol. 53 (17), 10445–10453. doi:10.1021/acs.est.9b02277
Wood, D. E., and Salzberg, S. L. (2014). Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 15 (3), R46. doi:10.1186/gb-2014-15-3-r46
World Health Organization (2023). Antimicrobial resistance: key facts. Available at: https://www.who.int/news-room/fact-sheets/detail/antimicrobial-resistance.
World Health Organization (2017). Antimicrobial resistance: the food chain. Available at: https://www.who.int/news-room/questions-and-answers/item/antimicrobial-resistance-in-the-food-chain.
Yang, Q., Zhang, J., Zhang, W., Wang, Z., Xie, Y., and Zhang, H. (2010). Influence of tetracycline exposure on the growth of wheat seedlings and the rhizosphere microbial community structure in hydroponic culture. J. Environ. Sci. Health B 45 (3), 190–197. doi:10.1080/03601231003613492
Zhang, H., Du, M., Jiang, H., Zhang, D., Lin, L., Ye, H., et al. (2015). Occurrence, seasonal variation and removal efficiency of antibiotics and their metabolites in wastewater treatment plants, Jiulongjiang River Basin, South China. Environ. Sci. Process Impacts 17 (1), 225–234. doi:10.1039/c4em00457d
Zhang, H., Li, X., Yang, Q., Sun, L., Yang, X., Zhou, M., et al. (2017). Plant growth, antibiotic uptake, and prevalence of antibiotic resistance in an endophytic system of pakchoi under antibiotic exposure. Int. J. Environ. Res. Public Health 14 (11), 1336. doi:10.3390/ijerph14111336
Zhang, S., Lu, Y. X., Zhang, J. J., Liu, S., Song, H. L., and Yang, X. L. (2020). Constructed wetland revealed efficient sulfamethoxazole removal but enhanced the spread of antibiotic resistance genes. Molecules 25 (4), 834. doi:10.3390/molecules25040834
Zhao, F., Chen, L., Yen, H., Li, G., Sun, L., and Yang, L. (2020). An innovative modeling approach of linking land use patterns with soil antibiotic contamination in peri-urban areas. Environ. Int. 134, 105327. doi:10.1016/j.envint.2019.105327
Zhao, F., Yang, L., Chen, L., Li, S., and Sun, L. (2019). Bioaccumulation of antibiotics in crops under long-term manure application: occurrence, biomass response and human exposure. Chemosphere 219, 882–895. doi:10.1016/j.chemosphere.2018.12.076
Zhao, Y., Li, M. C., Konaté, M. M., Chen, L., Das, B., Karlovich, C., et al. (2021). TPM, fpkm, or normalized counts? A comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository. J. Transl. Med. 19 (1), 269. doi:10.1186/s12967-021-02936-w
Keywords: sulfamethoxazole, basil, cilantro, spinach, antibiotic resistance
Citation: Ren J, Lu H, Lu S and Huang Z (2024) Impacts of sulfamethoxazole stress on vegetable growth and rhizosphere bacteria and the corresponding mitigation mechanism. Front. Bioeng. Biotechnol. 12:1303670. doi: 10.3389/fbioe.2024.1303670
Received: 28 September 2023; Accepted: 02 January 2024;
Published: 08 February 2024.
Edited by:
Yong Jiang, Fujian Agriculture and Forestry University, ChinaReviewed by:
Ziyang Zhao, Xi’an University of Posts and Telecommunications, ChinaQiuheng Zhu, Nanjing Hydraulic Research Institute Nanjing, China
Jiali Chang, Leshan Normal University, China
Copyright © 2024 Ren, Lu, Lu and Huang. 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: Hongbin Lu, luhongbin@sxicc.ac.cn; Shaoyong Lu, sylu@craes.org.cn