Jun Liu
Wen-Yu Qi1
Chao Song
Shu-Guang Wang- 1Shandong Key Laboratory of Water Pollution Control and Resource Reuse, School of Environmental Science and Engineering, Shandong University, Qingdao, China
- 2College of Agriculture and Forestry Science, Linyi University, Linyi, China
- 3State Key Laboratory of Microbial Technology, Shandong University, Qingdao, China
Selenium (Se) is an essential trace element in the human body. Se-enriched agricultural products, obtained by applying Se fertilizer, are important sources of Se supplement. However, Se fertilizer may cause a series of environmental problems. This study investigated the transformation of exogenous selenium nanoparticles (SeNPs) and selenite (SeO32–) in soil and explored their effects on soil microbial community and typical microorganisms. SeNPs exhibited a slow-release effect in soil, which promoted the growth of soil microorganisms and enriched soil probiotics. SeO32– was converted to a stable and low toxic state in soil, increasing persistent free radicals and decreasing microbial abundance and diversity. The influences of SeNPs and SeO32– on two typical soil microorganisms (Bacillus sp. and Escherichia coli) were also evaluated, and SeNPs were more difficult to enter into microorganisms directly, with lower toxicity and higher safety. These results indicated that SeNPs were a more environment-friendly Se additive for agriculture applications. This work provides useful information for better understanding the environmental fate and behavior of Se fertilizer in the soil.
Introduction
Selenium (Se) is a micronutrient essential for human health (Rayman, 2000). The deficiency of Se can lead to cardiovascular disease, oxidative stress, decreased thyroid and immune functions, and an increased risk of various cancers (Zhou et al., 2018; Kumar and Prasad, 2021). Dietary uptake of Se from agricultural products is a major source of Se for humans (Evans et al., 2016). However, natural Se is insufficient in most areas of China, and billions of people are estimated to suffer Se deficiency (Terry et al., 2000; Tan J. et al., 2002). Therefore, Se-rich agricultural products are usually considered a feasible and effective way to improve the Se content in human bodies (Bhatia et al., 2013).
Se fertilizers, including inorganic [e.g., selenate, selenite (SeO32–)] and organic (e.g., Se amino acid) varieties, have been widely used in soils to produce Se-rich crops (Alfthan et al., 2015). Generally, Se fertilizers effectively increase the Se content in plants (Kumar and Prasad, 2021). In addition, Se fertilizers could also promote the growth and development of plants, especially under stressful conditions such as high salinity and high temperature (Ashraf et al., 2018; Djanaguiraman et al., 2018). The early form of Se fertilizers was mainly SeO32–. However, SeO32– is highly toxic in nature due to high bioavailability and high mobility (Rosenfeld et al., 2017). As new Se supplements, selenium nanoparticles (SeNPs), elemental form of Se at the nanoscale, have been demonstrated to be advantageous to the environment for their low toxicity and high efficiency (Skalickova et al., 2017). Owing to the application of exogenous Se fertilizers, the amount of Se in plow layer soil is increasing (Tan L.C. et al., 2016).
Soil microorganisms, especially rhizosphere microorganisms, are active in nutrient cycling and are key to the functioning of agricultural systems (Panke-Buisse et al., 2015). They play important roles in the degradation of organic matter, nitrogen fixation, and phosphate solubilization, all of which provide nutrients for plants (Martin et al., 2016; Valyi et al., 2016; Rosenblueth et al., 2018). In addition, soil microbes could facilitate plants to resist environmental stresses, such as drought (Guerena et al., 2019), extreme temperature (Liu et al., 2019), soil salinity (Ben Laouane et al., 2019), acidity (Tullio et al., 2019), alkalinity (Abd-Alla et al., 2014), and heavy metals (Benjelloun et al., 2019). Some microorganisms can also protect crops against diseases via producing antifungal metabolites (Weller et al., 2002). Moreover, soil microorganisms play an important role in the transformation and cycle of Se in nature. For example, some soil microorganisms (e.g., Bacillus oryziterrae, Streptomyces sp.) have been reported to metabolize and transform SeO32– to SeNPs (Bao et al., 2016; Tan Y. et al., 2016). Therefore, a stable and healthy soil microbial community is crucial for both plant growth and Se cycle. However, the effect of exogenous Se fertilizer on local soil microorganisms remains unknown. Therefore, it is necessary to explore the effect of Se fertilizer on soil microorganisms.
In this study, we investigated the transformation of exogenous SeNPs and SeO32– in soil and their effects on local soil microbial community and typical microorganisms. The transformation of Se in soil was determined by a continuous extraction method, and the changes in the abundance and diversity of soil microorganisms were analyzed using quantitative PCR (qPCR) and high-throughput sequencing. The core microbial population in the presence of exogenous Se was analyzed using linear discriminant analysis (LDA) effect size (LEfSe) and correlation network analysis. To further explore the influencing mechanism of SeNPs and SeO32– on the soil microbial community, persistent free radicals were detected using an electron paramagnetic resonance (EPR). In addition, the responses of two typical soil microorganisms (Bacillus sp. and Escherichia coli) to exogenous Se were investigated to further explore the biological effects of SeNPs and SeO32–. This work completed a limited study on the transformation and effects of Se in soil, but the approach could be helpful for usage of SeNPs as a Se fertilizer for Se-rich agricultural products.
Materials and Methods
Reagents
Sodium SeO32–, sodium hydroxide (NaOH), dipotassium phosphate (K2HPO4), monopotassium phosphate (KH2PO4), potassium persulfate (K2S2O8), sodium sulfite (Na2SO3), n-caprylic alcohol, acetic acid, and hydrochloric acid were obtained from Sinopharm (Shanghai, China). SeNPs were synthesized and purified using a recently described method (Shakibaie et al., 2010). In brief, the SeNPs were biologically reduced from SeO32– by microorganisms. The reduction product was washed and disrupted in liquid nitrogen. The suspension was then liquid–liquid extracted by n-caprylic alcohol. The purified SeNPs were collected from water phase. More details of synthesis purification and SEM observation were shown in Supplementary Text 1.
Experimental Design
The soil used in this study was collected from a farmland in Jimo, Shandong province, and the physical and chemical properties of soil were shown in Supplementary Table 1. Five experiments with different Se treatment groups were established: (1) Control group, in which ddH2O was applied; (2) Low SeO32– group, in which 1 mmol/kg soil SeO32– solution was applied; (3) High SeO32– group, in which 10 mmol/kg soil SeO32– solution was applied; (4) Low SeNPs group, in which 1 mmol/kg soil SeNP suspension was applied; (5) High SeNPs group, in which 10 mmol/kg soil SeNP suspension was applied. The experiments were conducted in a greenhouse at 22°C, and the photoperiod was 16 h light and 8 h dark. Three soil samples were collected before all treatments were applied and designated as samples for Day 0. Soil samples of each treatment were collected in triplicate on Days 7 and 30 for further analysis.
Determination of Available Forms of Selenium
The contents and existing forms of Se in soil samples were determined using a modified continuous extraction method (Martens and Suarez, 1997; Wright et al., 2003; Kulp and Pratt, 2004; Keskinen et al., 2009). Briefly, 1,000-g soil sample was placed in a 50-ml centrifuge tube and extracted using different solutions with a 1:10 solid/liquid ratio for each step. The water-soluble Se was extracted using H2O, and the ligand-exchangeable Se was extracted using 0.1 mol/L K2HPO4-KH2PO4 (pH 7.0). The alkali-soluble Se was extracted by 0.1 mol/L NaOH, and the elemental Se was extracted with the combination of 1 mol/L Na2SO3 (pH was adjusted to 7.0 by hydrochloric acid) and ultrasonication. The acid-soluble Se was extracted using acetic acid. The Se concentration in the extracts was detected with an atomic fluorescence photometer (AFS-933, Titan Instrument, Beijing, China). More details were described in Supplementary Text 2.
DNA Extraction
Here, 0.25-g fresh soil sample was weighed to extract the total soil DNA. The extraction was performed using a TIANamp Soil DNA Kit (TIANGEN Biotech, Beijing, China). The DNA extracts were quantified and quality checked using a NanoDrop One UV–VIS Spectrophotometer (Thermo Fisher Scientific, United States). The DNA extracts were stored at −80°C for further analysis.
Quantification of Total Bacteria in Soil and High-Throughput Sequencing
The abundance of total bacteria was assessed by quantification of 16S rRNA gene copies through qPCR. A universal primer set K90q/K94q (Supplementary Table 2) targeting bacterial 16S rRNA was employed. The protocol of qPCR followed a previously described method (Liu et al., 2020).
The V4–V5 region of the bacteria 16S rRNA gene was amplified with the primer set 515F/907R (details shown in Supplementary Table 2). In this step, an 8-bp-long barcode unique to each sample was appended. The PCR product was purified and quantified for the library generation. Then, the amplicon library was sequenced using an Illumina MiSeq platform at Shanghai Biozeron Biological Technology Co., Ltd. (Shanghai, China). Raw sequencing data were uploaded in NCBI Sequence Read Archive (SRA) database under accession number PRJNA664643.
Bioinformatics and Statistical Analysis
High-throughput sequences were processed using the QIIME2 platform (version 2019.10) (Bolyen et al., 2019). Briefly, raw sequencing data were demultiplexed, filtered, and denoised. The denoised sequences were then clustered into operational taxonomic units (OTUs) with a 97% similarity. The representative sequences of each OTU were aligned against the SILVA database (SSU 135) with a confidence threshold of 70% to identify the taxonomy (Quast et al., 2013).
Alpha diversity indexes including Chao1, Shannon, and Simpson indices were calculated following the Vegan package in R. Principal coordinates analysis (PCoA) was calculated following the ape and GUniFrac package in R. LEfSe and correlation network analysis was conducted through the MicrobiomeAnalyst platform (Chong et al., 2020). LEfSe analysis was calculated using a significance level of p < 0.05 and an LDA score of 2. Correlation network analysis was performed at the genus level using the SparCC method with a significance level of p < 0.05 and a correlation threshold of 0.8 (Friedman and Alm, 2012; Jiang et al., 2017). According to the species classification data, the effects of Se on the function of soil microorganisms were predicted using the Tax4Fun2 algorithm.
Measurements of Soil-Persistent Free Radicals
Fresh soil samples of 0.25 g were collected for EPR analysis. The samples were loaded into a quartz EPR tube, placed into the EPR resonator, and analyzed with a Bruker EMXnano spectrometer (Karlsruhe, Germany). The parameters were set as follows: center field, 3,420 G; sweep width, 80 G; sweep time, 40.06 s; receiver gain, 40 dB; modulation amplitude, 2.000 G; microwave frequency, 9.8 GHz; and microwave attenuation, 20.00 dB.
Effect of Selenium Fertilizers on Typical Soil Microorganisms
Two typical soil microorganisms, namely, Gram-negative bacteria Bacillus sp. and Gram-positive bacteria E. coli, were chosen to investigate responses of different soil bacteria to exogenous Se. Briefly, the two bacteria were inoculated in a Luria-Bertani (LB) medium at a ratio of 1:100. SeNPs or SeO32– at a concentration of 0, 0.1, 0.5, 1, 5, and 10 mmol/L were added, respectively, and the mixed media were incubated at 37°C, 180 rpm for 24 h. OD600 was measured to characterize the biomass. Reactive oxygen species level was determined using Reactive Oxygen Species Assay Kit (Solarbio, Shanghai, China) following the manufacturer’s instructions. Glutathione peroxidase (GSH-Px) activity was determined using a Micro Glutathione Peroxidase Assay Kit according to the manufacturer’s protocol (Solarbio, Shanghai, China). Protein concentrations were determined following the bicinchoninic acid assay method (Smith et al., 1985).
Results and Discussion
Transformation of Selenium Nanoparticles and Selenite in Soil
The transformation of SeNPs and SeO32– in soil was analyzed with a continuous extraction method. As shown in Figure 1, the element Se decreased significantly along time in the High SeNP group from 248.58 mg/kg on Day 7 to 229.93 mg/kg on Day 30 (p < 0.001). In contrast, there was no significant difference in Se content in the Low SeNP group (29.98 mg/kg on Day 7, 28.99 mg/kg on Day 30). No soluble Se was detected in treatments with SeNPs. In both High SeO32– group and Low SeO32– group, the dominant form of Se was water-soluble SeO32–, and its concentration decreased significantly along time (p < 0.001). Alkali-soluble Se (-II) was the second predominant Se form, which was reported to be mainly presented as organic-associated Se in soil (Keskinen et al., 2009). In addition, a relatively high elemental Se was observed in the High SeO32– group. Studies have shown that SeO32– is affected by redox potential in soil and can be reduced gradually by soil microorganisms, where the toxicity of Se decreases with decreasing valence (Piacenza et al., 2017; Wang et al., 2017; Zhai et al., 2021).
Figure 1. Variation of soil-available selenium (Se) fractions. (A) High selenium nanoparticles (SeNPs): the samples collected from the group treated with 10 mmol/kg soil SeNPs. (B) High selenite (SeO32–): the samples collected from the group treated with 10 mmol/kg soil SeO32–. (C) Low SeNPs: the samples collected from the group treated with 1 mmol/kg soil SeNPs. (D) Low SeO32–: the samples collected from the group treated with 1 mmol/kg soil SeO32–. Significant differences were analyzed by ANOVA; ∗∗∗p < 0.001.
Effects of Selenium Nanoparticles and Selenite on the Abundance and Diversity of Soil Microorganisms
To explore the effect of SeNPs and SeO32– on the soil microbial community, soil samples were collected on Days 7 and 30. The abundance of soil bacteria was quantified by qPCR. The abundance of bacteria in the Control group remained relatively stable at 2.78–2.37 ∗ 1011 copies/g soil (Figure 2A). In the Low SeNP group, the abundance of soil bacteria increased by 146% (Day 7) and 171% (Day 30) compared with that of the Control group. A relatively smaller increase in the number of bacteria was observed in the High SeNP group, which was 129 and 136% on Days 7 and 30, respectively. However, the abundance of soil bacteria significantly decreased in treatments with SeO32–, which was decreased by 45% (Day 7) and 55% (Day 30) in the Low SeO32– group. In the High SeO32– group, the higher concentration of SeO32– exhibited a stronger adverse effect on soil bacteria, where the abundance of bacteria was decreased by 70% on Day 30 compared with that of the Control group.
Figure 2. Abundance and diversity of the microbiota in each experimental group. (A) Total soil bacteria abundance. (B) Shannon Index. (C) Principal coordinates analysis (PCoA) plots of Bray–Curtis distances at the genus level. (D) PCoA plots of weighted UniFrac distances at the genus level. Control: the samples collected from the control group without selenium (Se) treatment. High selenite (SeO32–): the samples collected from the group treated with 10 mmol/kg soil SeO32–. Low SeO32–: the samples collected from the group treated with 1 mmol/kg soil SeO32–. High selenium nanoparticles (SeNPs): the samples collected from the group treated with 10 mmol/kg soil SeNPs. Low SeNPs: the samples collected from the group treated with 1 mmol/kg soil SeNPs. ∗ in (C,D) shows the data point of untreated soil. Significant differences were analyzed by ANOVA; ∗∗p < 0.01. Different letters indicate a significant difference (p < 0.05).
In this study, SeNPs promoted the development of soil microorganisms. It has been reported that SeNPs can play an antioxidant role in organisms by enhancing the activities of key enzymes in the antioxidant defense system (Zhai et al., 2017), activating antioxidant selenase (Guo et al., 2020), inhibiting protein glycosylation (Yu et al., 2015), or directly scavenging free radicals in vivo, therein promoting microbial growth. SeNPs are more stable in the environment and exert lower toxicity to microorganisms as a biological reduction product. However, the transformation of SeO32– in and out of microorganisms usually produces reactive oxygen species (Kessi and Hanselmann, 2004) that cause the toxicity of SeO32–.
The effects of SeNPs and SeO32– on soil microbial community structure were further analyzed by 16S rRNA amplicon sequencing. A total of 1,488,614 sequences were obtained and assigned to 4,069 OTUs according to 97% similarity. No significant difference in α diversity was observed between the SeNP groups and the Control group (Figure 2B and Supplementary Table 3). According to PCoA of β diversity distance between samples (Figures 2C,D), the distance between the High/Low SeNP groups and the Control group was smaller than that between the SeO32– groups and the Control group. These results indicated that SeNPs exhibited no significant disturbance to the diversity of the soil microbial community.
However, there was a significant difference in the α diversity indexes between the High/Low SeO32– groups and the other groups (p < 0.01). The observed OTU number, Chao1 index, Simpson index, and Shannon index were all obviously reduced (Figure 2B and Supplementary Table 3). The High/Low SeO32– groups were separated from the other groups (Figures 2C,D). These results indicated that SeO32– affected the species and abundance of soil microorganisms. Further analysis of β diversity revealed a significant separation between the SeO32– groups and the Control group in the PCoA results based on the Bray–Curtis distance, which correlates with species abundance (Bray and Curtis, 1957; Figure 2C). However, PCoA results based on unweighted UniFrac distance, which is related to the developmental composition of species (Lozupone et al., 2011), showed no significant grouping differences, and the degree of interpretation among the principal components was low (Figure 2D). The different results of the two diversity distances implied that the SeO32– caused an abundance decrease but not extinction of soil microorganisms. It suggests that the effect of Se on soil microflora may not aim at specific species but be broad spectrum to all the soil microorganisms.
Species Evolution and Soil Function Prediction in the Presence of Selenium Nanoparticles and Selenite
As shown in Figure 3A, Proteobacteria dominated in soil samples from all treatment groups, with a relative abundance exceeding 50%. Significant differences were observed between the experimental groups at the genus level (Figure 3B). To quantitatively analyze the species changes, LEfSe method was used to analyze the species abundance at the genus level on Days 7 and 30 (Figures 3C,D), aiming to screen the species with significant differences in each experimental group. Under the given conditions, 205 and 267 significantly different species were identified, respectively. The relative abundance of several genera, e.g., Acidibacter, Acidipila, and Bradyrhizobium, was higher in the Control group, indicating that their relative abundance decreased under the stress of SeNPs and SeO32–. The relative abundance of Rhodanobacter, Bacillus, Alicyclobacillus, and Granulicella increased in the SeNP groups, while the relative abundance of Chujaibacter, Candidimonas, Jatrophihabitans, and Nevskia increased in the SeO32– groups. The increased relative abundance may be due to the promotion effect of SeNPs or Se resistance of microorganisms.
Figure 3. Taxonomy composition of the microbiota in each experimental group. (A) The relative abundance of major bacteria phyla. (B) The relative abundance of major bacteria genera. (C,D) The top 15 significantly changed genera identified by linear discriminant analysis effect size (LEfSe) on Days 7 (C) and 30 (D). Control: the samples collected from the control group without selenium (Se) treatment. High selenite (SeO32–): the samples collected from the group treated with 10 mmol/kg soil SeO32–. Low SeO32–: the samples collected from the group treated with 1 mmol/kg soil SeO32–. High selenium nanoparticles (SeNPs): the samples collected from the group treated with 10 mmol/kg soil SeNPs. Low SeNPs: the samples collected from the group treated with 1 mmol/kg soil SeNPs.
The influence of SeNPs and SeO32– on the soil microbial community was studied with a co-occurrence network method (Figure 4). Based on SparCC method, the correlation of soil microbial abundance was calculated, and then the graph inference of the network and the estimation of several topological characteristics were carried out. In the SeNP groups, the association network has 89 nodes and 308 connections, which can be classified into 16 constituent units (Figures 4A,C). In the SeO32– groups, there were 119 nodes, 442 connections, and 25 basic units (Figures 4B,D). To summarize, the impact of SeNPs on soil microbial network is less than that of SeO32–. Based on intermediate centrality, the key microorganisms in each network were identified. In SeNP groups, Tuberibacillus and Telmatospirillum were the microorganisms located in the core of the network, while Bryobacter and Mizugakibacter were the microorganisms located in the core of the network in SeO32– groups.
Figure 4. Co-occurrence network analysis showing the biological interactions at the family level under the treatment of selenite (SeO32–) (A) and selenium nanoparticles (SeNPs) (B). The size of the nodes is relative abundance. The scatter plot shows criteria of selecting for the keystone taxa in SeO32– (C) and SeNP (D) networks.
The enriched microorganisms (e.g., Bacillus, Alicyclobacillus, Granulicella) in SeNP groups have been previously reported as soil probiotics (Zelenikhin et al., 2013; Long et al., 2018). In addition, Tuberibacillus and Telmatospirillum at the core of the microbial network are involved in the decomposition of soil organic matter and utilization of carbon, nitrogen, phosphorus, and other elements (Hatayama et al., 2006; Xue et al., 2017; Ding et al., 2019), indicating that the addition of SeNPs induces active probiotics related to basic soil functions and improves the soil microbial community. Chujaibacter, Candidimonas, Jatrophihabitans, Nevskia, Bryobacter, and Mizugakiibacter, which were enriched in the SeO32– groups or lived in the core of the network, have the ability to treat heavy metals and resist UV, antibiotics, and other biological or abiotic stresses (Vaz-Moreira et al., 2011; Větrovský and Baldrian, 2015; Prakash et al., 2020; Sánchez-Soto et al., 2020). These results imply that SeO32– appeared to be toxic to the soil microbial community.
The effects of exogenous Se on soil functions were predicted based on the Tax4Fun2 algorithm (Supplementary Figure 2). There was no significant difference between treatments, indicating that exogenous Se exhibited little impact on soil function. Therefore, the soil microbial community can effectively resist the external interface and maintain the ecological functions of the soil, reflecting the robustness of the soil microbial community function.
Persistent Free Radicals Release From Soil Induced by Selenium Nanoparticles and Selenite
In this study, exogenous SeNPs and SeO32– could change soil microbial community obviously. To explore the influencing mechanism of SeNPs and SeO32– on the soil microbial community, persistent free radicals in the soil were determined, as the migration and transformation of Se may induce the release of free radicals (Kessi and Hanselmann, 2004). By analyzing the G-factor of free radical signal, it was found that the persistent free radicals in soil were mainly oxygen-center free radicals and semi-quinone free radicals (Khachatryan and Dellinger, 2011; Figure 5 and Supplementary Table 4). There was no significant difference in free radical signals between the SeNP groups and the Control group. However, in the SeO32– groups, the persistent free radical signal in the soil was significantly enhanced on Day 7, but there was no significant difference on Day 30. Nevertheless, the bacterial community showed a more significant difference on Day 30, indicating that the corresponding response of the soil microbial community to exogenous Se compounds lagged behind the migration and transformation of exogenous Se in soil. The migration and transformation of exogenous Se in the environment, especially SeO32–, is mainly due to chemical reactions, but the changes of soil microbial community need a longer period compared with chemical reactions.
Figure 5. Electron paramagnetic resonance (EPR) spectra of each experimental group on Days 7 and 30.
Effects of Selenium Nanoparticles and Selenite on Typical Soil Microorganisms
In this study, the application of SeNPs and SeO32– exhibited different effects on the soil microbial community. SeNPs promoted the propagation of soil probiotics, while SeO32– induced the enrichment of microorganisms with a higher tolerance to environmental stress. In addition, we found that the relative abundance of Gram-negative bacteria increased more than that of Gram-positive bacteria after SeNP and SeO32– treatments. Therefore, we selected two representative soil microorganisms (Bacillus sp. and E. coli) to investigate their response to exogenous Se. As shown in Figure 6A, the microbial biomass increased with the increase of SeNP concentration, while SeO32– at different concentrations had little influence on the microbial biomass. Regardless of the Se concentration, both tested bacteria showed higher growth in treatments with SeNPs compared with that of SeO32– groups. In addition, experimental groups with E. coli always showed higher biomass than that with Bacillus sp. Similar to the results of biomass, the total protein concentration in SeNP treatment groups was higher than that of SeO32– treatment group, and the total protein content of E. coli was higher than that of Bacillus sp. (Figure 6B). The above results indicated that SeNPs promoted the growth of microorganisms at a certain concentration compared with SeO32–, and the promotion was more significant in E. coli. We further determined the intracellular reactive oxygen species of different treatment groups and found that the intracellular reactive oxygen species of Bacillus sp. were lower than that of E. coli (Figure 6D). It implies that E. coli had a higher ability to resist oxidative stress, which is consistent with its higher biomass.
Figure 6. Effects of selenium nanoparticles (SeNPs) and selenite (SeO32–) on typical soil microorganisms. (A) OD600. (B) Total protein concentration. (C) Glutathione peroxidase (GSH-Px) activity. (D) Reactive oxygen species.
GSH-Px is a kind of antioxidant enzyme with Se atom as its active center, and its activity can characterize the bioavailability of Se (Yu and Wei, 2013). The GSH-Px activity in the SeNP groups was significantly lower than that in the SeO32– treatment groups (Figure 6C). Although GSH-Px activity was significantly different (p < 0.05), there was no significant correlation between enzyme activity and Se concentration (Figure 6C and Supplementary Figure 3). The GSH-Px activity of Bacillus sp. was higher than that of E. coli, partly resulting in the lower amount of intracellular reactive oxygen species in Bacillus sp.
The results indicated that SeNPs were less bioavailable than SeO32–. It has been reported that SeO32– can be easily transported across the membrane via binding to the sulfhydryl sites on the cell membrane or by phosphate transfer protein (Yu et al., 2018; Zhu et al., 2020). In contrast, SeNPs are in a reduced state. The Se atom is coated with a layer of extracellular polymer (Avendaño et al., 2016), making it difficult to enter the microorganism directly. It may partly explain the lower toxicity of SeNPs.
Conclusion
In this study, the migration and transformation of two Se fertilizers (i.e., SeNPs and SeO32–) in soil were determined. Their effects on the soil microbial community were assessed, and the possible mechanism was analyzed. It was found that SeNPs were more stable in the soil, while SeO32– would be gradually transformed into a more stable and lower toxicity state in the soil. SeNPs promoted the propagation of soil probiotics, while SeO32– reduced the abundance and diversity of microorganisms possibly due to the toxicity of the persistent free radicals released during the transformation of SeO32–. Experiments with two typical soil microorganisms showed that SeNPs showed lower toxicity and higher safety because SeNPs enter microorganisms slowly. These results indicated that SeNPs were a more environment-friendly Se additive for agricultural application. This work provides useful information for better understanding of the environmental fate and behavior of Se fertilizer in the soil.
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 below: https://www.ncbi.nlm.nih.gov/bioproject/ PRJNA664643.
Author Contributions
JL and QL conceptualized and designed the study. JL and W-YQ carried out the experimental work and described in the article. JL analyzed, interpreted, and visualized the data. JL, CS, and HC drafted the manuscript. S-GW acquired the funding. All authors contributed to revise the manuscript and approved the submitted version, with proofreading performed by JL.
Funding
This work was supported by the National Natural Science Foundation of China (U20A20146) and the National Key R&D Program of China (2016YFC0502802).
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.
Acknowledgments
We thank Haiyan Yu, Xiaomin Zhao, and Sen Wang from the State Key Laboratory of Microbial Technology, Shandong University, for the assistance in microimaging of SEM analysis. We extend their appreciation to Wenshan Qi from Biozeron Co., Ltd., for his technical assistance of high-throughput sequencing in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.746046/full#supplementary-material
References
Abd-Alla, M. H., El-Enany, A.-W. E., Nafady, N. A., Khalaf, D. M., and Morsy, F. M. (2014). Synergistic interaction of Rhizobium leguminosarum bv. viciae and arbuscular mycorrhizal fungi as a plant growth promoting biofertilizers for faba bean (Vicia faba L.) in alkaline soil. Microbiol. Res. 169, 49–58. doi: 10.1016/j.micres.2013.07.007
Alfthan, G., Eurola, M., Ekholm, P., Venäläinen, E.-R., Root, T., Korkalainen, K., et al. (2015). Effects of nationwide addition of selenium to fertilizers on foods, and animal and human health in Finland: from deficiency to optimal selenium status of the population. J. Trace Elem. Med. Biol. 31, 142–147. doi: 10.1016/j.jtemb.2014.04.009
Ashraf, M. A., Akbar, A., Parveen, A., Rasheed, R., Hussain, I., and Iqbal, M. (2018). Phenological application of selenium differentially improves growth, oxidative defense and ion homeostasis in maize under salinity stress. Plant Physiol. Biochem. 123, 268–280. doi: 10.1016/j.plaphy.2017.12.023
Avendaño, R., Chaves, N., Fuentes, P., Sánchez, E., Jiménez, J. I., and Chavarría, M. (2016). Production of selenium nanoparticles in Pseudomonas putida KT2440. Sci. Rep. 6:37155. doi: 10.1038/srep37155
Bao, P., Xiao, K.-Q., Wang, H.-J., Xu, H., Xu, P.-P., Jia, Y., et al. (2016). Characterization and potential applications of a selenium nanoparticle producing and nitrate reducing bacterium Bacillus oryziterrae sp. nov. Sci. Rep. 6:34054. doi: 10.1038/srep34054
Ben Laouane, R., Meddich, A., Bechtaoui, N., Oufdou, K., and Wahbi, S. (2019). Effects of arbuscular mycorrhizal fungi and rhizobia symbiosis on the tolerance of Medicago Sativa to salt stress. Gesunde Pflanz. 71, 135–146. doi: 10.1007/s10343-019-00461-x
Benjelloun, I., Alami, I. T., Douira, A., and Udupa, S. M. (2019). Phenotypic and genotypic diversity among symbiotic and non-symbiotic bacteria present in chickpea nodules in Morocco. Front. Microbiol. 10:1885. doi: 10.3389/fmicb.2019.01885
Bhatia, P., Aureli, F., D’Amato, M., Prakash, R., Cameotra, S. S., Nagaraja, T. P., et al. (2013). Selenium bioaccessibility and speciation in biofortified Pleurotus mushrooms grown on selenium-rich agricultural residues. Food Chem. 140, 225–230. doi: 10.1016/j.foodchem.2013.02.054
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of southern wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268
Chong, J., Liu, P., Zhou, G., and Xia, J. (2020). Using microbiomeanalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat. Protoc. 15, 799–821. doi: 10.1038/s41596-019-0264-1
Ding, Y., Xiong, J., Zhou, B., Wei, J., Qian, A., Zhang, H., et al. (2019). Odor removal by and microbial community in the enhanced landfill cover materials containing biochar-added sludge compost under different operating parameters. Waste Manag. 87, 679–690. doi: 10.1016/j.wasman.2019.03.009
Djanaguiraman, M., Belliraj, N., Bossmann, S. H., and Prasad, P. V. V. (2018). High-temperature stress alleviation by selenium nanoparticle treatment in grain sorghum. ACS Omega 3, 2479–2491. doi: 10.1021/acsomega.7b01934
Evans, R. D., Grochowina, N. M., Basu, N., O’Connor, E. M., Hickie, B. E., Rouvinen-Watt, K., et al. (2016). Uptake of selenium and mercury by captive mink: results of a controlled feeding experiment. Chemosphere 144, 1582–1588. doi: 10.1016/j.chemosphere.2015.09.096
Friedman, J., and Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687. doi: 10.1371/journal.pcbi.1002687
Guerena, D. T., Lehmann, J., Thies, J. E., Vanek, S., Karanja, N., and Neufeldt, H. (2019). Nodulation of beans with inoculant carriers from pyrolyzed and non-pyrolyzed sugarcane bagasse in response to different pre-planting water availability. Appl. Soil Ecol. 143, 126–133. doi: 10.1016/j.apsoil.2019.06.010
Guo, L., Xiao, J., Liu, H., and Liu, H. (2020). Selenium nanoparticles alleviate hyperlipidemia and vascular injury in ApoE-deficient mice by regulating cholesterol metabolism and reducing oxidative stress. Met. Integr. Biometal Sci. 12, 204–217. doi: 10.1039/c9mt00215d
Hatayama, K., Shoun, H., Ueda, Y., and Nakamura, A. (2006). Tuberibacillus calidus gen. nov., sp nov., isolated from a compost pile and reclassification of Bacillus naganoensis Tomimura et al. 1990 as Pullulanibacillus naganoensis gen. nov., comb. nov and Bacillus laevolacticus Andersch et al. 1994 as Sporolactobacillus laevolacticus comb. nov. Int. J. Syst. Evol. Microbiol. 56, 2545–2551. doi: 10.1099/ijs.0.64303-0
Jiang, Y., Liu, M., Zhang, J., Chen, Y., Chen, X., Chen, L., et al. (2017). Nematode grazing promotes bacterial community dynamics in soil at the aggregate level. ISME J. 11, 2705–2717. doi: 10.1038/ismej.2017.120
Keskinen, R., Ekholm, P., Yli-Halla, M., and Hartikainen, H. (2009). Efficiency of different methods in extracting selenium from agricultural soils of Finland. Geoderma 153, 87–93. doi: 10.1016/j.geoderma.2009.07.014
Kessi, J., and Hanselmann, K. W. (2004). Similarities between the abiotic reduction of selenite with glutathione and the dissimilatory reaction mediated by Rhodospirillum rubrum and Escherichia coli. J. Biol. Chem. 279, 50662–50669. doi: 10.1074/jbc.M405887200
Khachatryan, L., and Dellinger, B. (2011). Environmentally persistent free radicals (EPFRs)-2. Are free hydroxyl radicals generated in aqueous solutions? Environ. Sci. Technol. 45, 9232–9239. doi: 10.1021/es201702q
Kulp, T. R., and Pratt, L. M. (2004). Speciation and weathering of selenium in upper cretaceous chalk and shale from South Dakota and Wyoming, USA. Geochim. Cosmochim. Acta 68, 3687–3701. doi: 10.1016/j.gca.2004.03.008
Kumar, A., and Prasad, K. S. (2021). Role of nano-selenium in health and environment. J. Biotechnol. 325, 152–163. doi: 10.1016/j.jbiotec.2020.11.004
Liu, J., Qi, W., Li, Q., Wang, S.-G., Song, C., and Yuan, X. (2020). Exogenous phosphorus-solubilizing bacteria changed the rhizosphere microbial community indirectly. 3 Biotech 10:164. doi: 10.1007/s13205-020-2099-4
Liu, Y.-S., Geng, J.-C., Sha, X.-Y., Zhao, Y.-X., Hu, T.-M., and Yang, P.-Z. (2019). Effect of rhizobium symbiosis on low-temperature tolerance and antioxidant response in alfalfa (Medicago sativa L.). Front. Plant Sci. 10:538. doi: 10.3389/fpls.2019.00538
Long, X.-E., Yao, H., Huang, Y., Wei, W., and Zhu, Y.-G. (2018). Phosphate levels influence the utilisation of rice rhizodeposition carbon and the phosphate-solubilising microbial community in a paddy soil. Soil Biol. Biochem. 118, 103–114. doi: 10.1016/j.soilbio.2017.12.014
Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J., and Knight, R. (2011). UniFrac: an effective distance metric for microbial community comparison. ISME J. 5, 169–172. doi: 10.1038/ismej.2010.133
Martens, D. A., and Suarez, D. L. (1997). Selenium speciation of soil/sediment determined with sequential extractions and hydride generation atomic absorption spectrophotometry. Environ. Sci. Technol. 31, 133–139. doi: 10.1021/es960214+
Martin, F., Kohler, A., Murat, C., Veneault-Fourrey, C., and Hibbett, D. S. (2016). Unearthing the roots of ectomycorrhizal symbioses. Nat. Rev. Microbiol. 14, 760–773. doi: 10.1038/nrmicro.2016.149
Panke-Buisse, K., Poole, A. C., Goodrich, J. K., Ley, R. E., and Kao-Kniffin, J. (2015). Selection on soil microbiomes reveals reproducible impacts on plant function. ISME J. 9, 980–989. doi: 10.1038/ismej.2014.196
Piacenza, E., Bulgarini, A., Lampis, S., Vallini, G., and Turner, R. J. (2017). “Biogenic SeNPs from Bacillus mycoides SelTE01 and Stenotrophomonas maltophilia SelTE02: Characterization with reference to their associated organic coating,” in Proceedings of the AIP Conference on NANOINNOVATION, (University of Calgary, Canada), 020005. doi: 10.1063/1.4997134
Prakash, O., Green, S. J., Singh, P., Jasrotia, P., and Kostka, J. E. (2020). Stress-related ecophysiology of members of the genus Rhodanobacter isolated from a mixed waste contaminated subsurface. Front. Environ. Sci. Eng. 15:23. doi: 10.1007/s11783-020-1315-0
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic. Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Rayman, M. P. (2000). The importance of selenium to human health. Lancet 356, 233–241. doi: 10.1016/S0140-6736(00)02490-9
Rosenblueth, M., Ormeno-Orrillo, E., Lopez-Lopez, A., Rogel, M. A., Reyes-Hernandez, B. J., Martinez-Romero, J. C., et al. (2018). Nitrogen fixation in cereals. Front. Microbiol. 9:1794. doi: 10.3389/fmicb.2018.01794
Rosenfeld, C. E., Kenyon, J. A., James, B. R., and Santelli, C. M. (2017). Selenium (IV,VI) reduction and tolerance by fungi in an oxic environment. Geobiology 15, 441–452. doi: 10.1111/gbi.12224
Sánchez-Soto, D., Agüero-Chapin, G., Armijos-Jaramillo, V., Perez-Castillo, Y., Tejera, E., Antunes, A., et al. (2020). ShadowCaster: compositional methods under the shadow of phylogenetic models to detect horizontal gene transfers in prokaryotes. Genes 11:756. doi: 10.3390/genes11070756
Shakibaie, M., Khorramizadeh, M. R., Faramarzi, M. A., Sabzevari, O., and Shahverdi, A. R. (2010). Biosynthesis and recovery of selenium nanoparticles and the effects on matrix metalloproteinase-2 expression. Biotechnol. Appl. Biochem. 56, 7–15. doi: 10.1042/BA20100042
Skalickova, S., Milosavljevic, V., Cihalova, K., Horky, P., Richtera, L., and Adam, V. (2017). Selenium nanoparticles as a nutritional supplement. Nutrition 33, 83–90. doi: 10.1016/j.nut.2016.05.001
Smith, P. K., Krohn, R. I., Hermanson, G. T., Mallia, A. K., Gartner, F. H., Provenzano, M. D., et al. (1985). Measurement of protein using bicinchoninic acid. Anal. Biochem. 150, 76–85. doi: 10.1016/0003-2697(85)90442-7
Tan, J., Zhu, W., Wang, W., Li, R., Hou, S., Wang, D., et al. (2002). Selenium in soil and endemic diseases in China. Sci. Total Environ. 284, 227–235. doi: 10.1016/S0048-9697(01)00889-0
Tan, L. C., Nancharaiah, Y. V., van Hullebusch, E. D., and Lens, P. N. L. (2016). Selenium: environmental significance, pollution, and biological treatment technologies. Biotechnol. Adv. 34, 886–907. doi: 10.1016/j.biotechadv.2016.05.005
Tan, Y., Yao, R., Wang, R., Wang, D., Wang, G., and Zheng, S. (2016). Reduction of selenite to Se (0) nanoparticles by filamentous bacterium Streptomyces sp. ES2-5 isolated from a selenium mining soil. Microb. Cell Fact. 15:157. doi: 10.1186/s12934-016-0554-z
Terry, N., Zayed, A. M., de Souza, M. P., and Tarun, A. S. (2000). Selenium in higher plants. Annu. Rev. Plant Physiol. Plant Mol. Biol. 51, 401–432. doi: 10.1146/annurev.arplant.51.1.401
Tullio, L. D., Gomes, D. F., Silva, L. P., Hungria, M., and da Silva Batista, J. S. (2019). Proteomic analysis of Rhizobium freirei PRF 81(T) reveals the key role of central metabolic pathways in acid tolerance. Appl. Soil Ecol. 135, 98–103. doi: 10.1016/j.apsoil.2018.11.014
Valyi, K., Mardhiah, U., Rillig, M. C., and Hempel, S. (2016). Community assembly and coexistence in communities of arbuscular mycorrhizal fungi. ISME J. 10, 2341–2351. doi: 10.1038/ismej.2016.46
Vaz-Moreira, I., Figueira, V., Lopes, A. R., De Brandt, E., Vandamme, P., Nunes, O. C., et al. (2011). Candidimonas nitroreducens gen. nov., sp. nov. and Candidimonas humi sp. nov., isolated from sewage sludge compost. Int. J. Syst. Evol. Microbiol. 61, 2238–2246. doi: 10.1099/ijs.0.021188-0
Větrovský, T., and Baldrian, P. (2015). An in-depth analysis of actinobacterial communities shows their high diversity in grassland soils along a gradient of mixed heavy metal contamination. Biol. Fertil. Soils 51, 827–837. doi: 10.1007/s00374-015-1029-9
Wang, X., Zhang, D., Pan, X., Lee, D.-J., Al-Misned, F. A., Mortuza, M. G., et al. (2017). Aerobic and anaerobic biosynthesis of nano-selenium for remediation of mercury contaminated soil. Chemosphere 170, 266–273. doi: 10.1016/j.chemosphere.2016.12.020
Weller, D. M., Raaijmakers, J. M., Gardener, B. B. M., and Thomashow, L. S. (2002). Microbial populations responsible for specific soil suppressiveness to plant pathogens. Annu. Rev. Phytopathol 40, 309–348. doi: 10.1146/annurev.phyto.40.030402.110010
Wright, M. T., Parker, D. R., and Amrhein, C. (2003). Critical evaluation of the ability of sequential extraction procedures to quantify discrete forms of selenium in sediments and soils. Environ. Sci. Technol. 37, 4709–4716. doi: 10.1021/es0342650
Xue, L., Ren, H., Li, S., Leng, X., and Yao, X. (2017). Soil bacterial community structure and co-occurrence pattern during vegetation restoration in karst rocky desertification area. Front. Microbiol. 8:2377. doi: 10.3389/fmicb.2017.02377
Yu, Q., Boyanov, M. I., Liu, J., Kemner, K. M., and Fein, J. B. (2018). Adsorption of selenite onto Bacillus subtilis: the overlooked role of cell envelope sulfhydryl sites in the microbial conversion of Se (IV). Environ. Sci. Technol. 52, 10400–10407. doi: 10.1021/acs.est.8b02280
Yu, S., Zhang, W., Liu, W., Zhu, W., Guo, R., Wang, Y., et al. (2015). The inhibitory effect of selenium nanoparticles on protein glycation in vitro. Nanotechnology 26:145703. doi: 10.1088/0957-4484/26/14/145703
Yu, Y., and Wei, J. (2013). Glutathione peroxidase and its synthesis mechanism. Acta Biophys. Sin. 29:724. doi: 10.3724/SP.J.1260.2013.30119
Zelenikhin, P. V., Makeeva, A. V., Lozhkin, A. P., Rodionov, A. A., Nguen, N., and Ilinskaya, O. N. (2013). Effect of Bacillus pumilus ribonuclease on the paramagnetic centers of microbial cells. Microbiology 82, 772–777. doi: 10.1134/S0026261714010172
Zhai, H., Kleawsampanjai, P., Wang, M., Qi, M., Liu, Y., Liu, N., et al. (2021). Effects of soil moisture on aging of exogenous selenate in three different soils and mechanisms. Geoderma 390:114966. doi: 10.1016/j.geoderma.2021.114966
Zhai, X., Zhang, C., Zhao, G., Stoll, S., Ren, F., and Leng, X. (2017). Antioxidant capacities of the selenium nanoparticles stabilized by chitosan. J. Nanobiotechnol. 15:4. doi: 10.1186/s12951-016-0243-4
Zhou, H., Wang, T., Li, Q., and Li, D. (2018). Prevention of keshan disease by selenium supplementation: a systematic review and meta-analysis. Biol. Trace Elem. Res. 186, 98–105. doi: 10.1007/s12011-018-1302-5
Keywords: selenium nanoparticles, selenite, selenium fertilizer, soil microorganism, environmental persistent free radicals
Citation: Liu J, Qi W-Y, Chen H, Song C, Li Q and Wang S-G (2021) Selenium Nanoparticles as an Innovative Selenium Fertilizer Exert Less Disturbance to Soil Microorganisms. Front. Microbiol. 12:746046. doi: 10.3389/fmicb.2021.746046
Received: 23 July 2021; Accepted: 20 August 2021;
Published: 13 September 2021.
Edited by:
Fayuan Wang, Qingdao University of Science and Technology, ChinaReviewed by:
Ting Hu, China Agricultural University, ChinaKe-Qing Xiao, University of Leeds, United Kingdom
Li Yanhui, Qingdao University, China
Copyright © 2021 Liu, Qi, Chen, Song, Li 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: Qiang Li, liqiang67@lyu.edu.en; Shu-Guang Wang, wsg@sdu.edu.cn