Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 06 October 2022
Sec. Terrestrial Microbiology
This article is part of the Research Topic Soil Microbiome Community and Functional Succession Mechanism Driven by Different Factors in Agricultural Ecology View all 43 articles

Co-analysis of cucumber rhizosphere metabolites and microbial PLFAs under excessive fertilization in solar greenhouse

Shuang Wang,,,Shuang Wang1,2,3,4Ting Bian,,,Ting Bian1,2,3,4Tong Wu,,,Tong Wu1,2,3,4Yidi Zhang,,,Yidi Zhang1,2,3,4Muhammad Awais,,,Muhammad Awais1,2,3,4Hongdan Fu,,,
Hongdan Fu1,2,3,4*Zhouping Sun,,,
Zhouping Sun1,2,3,4*
  • 1College of Horticulture, Shenyang Agricultural University, Shenyang, China
  • 2Key Laboratory of Protected Horticulture of the Education Ministry, Liaoning, China
  • 3National and Local Joint Engineering Research Center of Northern Horticultural Facilities Design and Application Technology, Shenyang, China
  • 4Collaborative Innovation Center of Protected Vegetable Surround Bohai Gulf Region, Shenyang, China

Fertilizer application is the most common measure in agricultural production, which can promote the productivity of crops such as cucumbers, but the problem of excessive fertilization occurs frequently in solar greenhouses. However, the effects of fertilization levels on cucumber rhizosphere soil microbes and metabolites and their relationships are still unclear. In order to determine how fertilization levels affect the rhizosphere microenvironment, we set up four treatments in the solar greenhouse: no-fertilization (N0P0K0), normal fertilization (N1P1K1), slight excessive fertilization (N2P2K2), and extreme excessive fertilization (N3P3K3). The results showed that fertilization treatments significantly increased cucumber yield compared to no-fertilization, but, the yield of N3P3K3 was significantly lower than that of N1P1K1 and N2P2K2. Fertilization levels had significant effects on rhizosphere microorganisms, and pH, NH4+-N and AP were the main environmental factors that affected the changes in microbial communities. The total PLFAs, the percentages of fungi and arbuscular mycorrhizal fungi (AMF) were significantly reduced and bacteria percentage was significantly increased in N3P3K3 compared to other fertilization treatments. Differential metabolites under different fertilization levels were mainly organic acids, esters and sugars. Soil phenols with autotoxic effect under fertilization treatments were higher than that of N0P0K0. In addition, compared with soil organic acids and alkanes of N0P0K0, N2P2K2 was significantly increased, and N3P3K3 was not significantly different. This suggested that cucumber could maintain microbial communities by secreting beneficial metabolites under slight excessive fertilization (N2P2K2). But under extremely excessive fertilization (N3P3K3), the self-regulating ability of cucumber plants and rhizosphere soil was insufficient to cope with high salt stress. Furthermore, co-occurrence network showed that 16:1ω5c (AMF) was positively correlated with 2-palmitoylglycerol, hentriacontane, 11-octadecenoic acid, decane,4-methyl- and D-trehalose, and negatively correlated with 9-octadecenoic acid at different fertilization levels. This indicated that the beneficial microorganisms in the cucumber rhizosphere soil promoted with beneficial metabolites and antagonized with harmful metabolites. But with the deepening of overfertilization, the content of beneficial microorganisms and metabolites decreased. The study provided new insights into the interaction of plant rhizosphere soil metabolites and soil microbiomes under the different fertilization levels.

Introduction

With the constantly increasing of greenhouse area in recent years, greenhouse vegetable cultivation has gradually become an important component and future trend in vegetable production (Zhang et al., 2017). However, some problems have also arisen accordingly. One of them is the overuse of fertilizer by farmers to increase agricultural output (Zhang et al., 2012). China accounts for about 9% of global cropland, but consumes more than 30% of the world’s fertilizers (Chen et al., 2019). Excessive fertilization can lead to soil salinization, acidification and groundwater pollution (Tang et al., 2018; Li et al., 2019; Bindraban et al., 2020; Chynoweth et al., 2021). Muhammad et al. (2008) and Yuan et al. (2007) explained that the accumulation of salinity in soil not only has a direct negative impact on crops. It also has a significant inhibitory effect on microbial respiration, resulting in impaired microbial metabolic function and indirectly affecting crop growth. Inappropriate soil inputs had been reported to disturb the balance between beneficial microbes and metabolites and pathogenic microbes and harmful metabolites (Shen et al., 2010; Chen et al., 2019). It was reported that soil microbial communities and soil metabolites composition were closely related to plant–soil health (Reinhold-Hurek et al., 2015; Dong et al., 2018). Therefore, research on how chemical fertilizer application affects soil microorganisms and metabolites in the rhizosphere is crucial for the long-term sustainability of solar greenhouses.

Soil microbes are key drivers of biogeochemical cycles, energy turnover, and plant growth (Mohanram and Kumar, 2019). Microbial diversity is one of the most important soil health indicators (Wang et al., 2020), it varies with plant, season, soil type, and fertilizer management mode (Garbeva et al., 2004; Zhang et al., 2012; Li et al., 2021). The key responses of soil microbiota to fertilizers can be revealed by analyzing microorganisms (Islam et al., 2011). The 10 year field trial results clearly show that long-term nitrogen fertilization affects soil pH, thereby significantly altering most dominant soil bacterial species (Ren et al., 2020). As well as excessive phosphorus fertilizers reduced microbial diversity and significantly altered microbial community structures and compositions, phosphorus levels also interfered with the complexity of soil bacterial and fungal symbiosis networks (Cheng et al., 2022). Moreover, soil metabolomics provides a window into microbial behavior, as the metabolites secreted by microbes and the residual decomposition products of microbes are the main components of the soil metabolome (Jones et al., 2014). Microbial metabolism includes many simultaneous anabolic and catabolic reactions that inform biochemical activities within the entire soil microbiota (Muhammad et al., 2008; Poeplau et al., 2016). Previous studies had shown that fertilization affects the diversity of soil microbial functions (Shen et al., 2010; Hao et al., 2020), and certain microbial metabolic processes, such as carbohydrate and amino acid metabolism (Li et al., 2021).

In addition to the metabolites secreted by microorganisms and the decomposition of microbial residues, root exudates of plants are also the main components of rhizosphere soil metabolites (Badri and Vivanco, 2009; Song et al., 2020; Smercina et al., 2021). Plant roots secrete a large number of low molecular weight and high molecular weight compounds into the surrounding soil, including sugars, organic acids, amino acids and other secondary metabolites. These rhizosphere sediments account for 10% of the plant net photosynthetically fixed carbon (Mohanram and Kumar, 2019). They are important carbon substrates for soil microorganisms near the roots (Wallace, 1994; Yuan et al., 2016), and as signaling molecules to influence the complex interactions between plants and rhizosphere microbes (Badri and Vivanco, 2009; Mohanram and Kumar, 2019; Phour et al., 2020). Plants can secrete root exudates and change root morphology to influence and recruit the specific microbial communities (Reinhold-Hurek et al., 2015; Dong et al., 2018). It was recently found that rhizosphere fungal community assemblage and species coexistence differ from bulk soils, caused by changes in root exudates under different nitrogen inputs (Wang et al., 2022). Therefore, understanding the effects of fertilizer usage on rhizosphere metabolism is helpful to further understand the effects of fertilizer on soil microbial community and its interactions with plants.

In order to clarify the effect of excessive NPK fertilization on soil microorganisms and metabolites, and to provide a basis for reasonable fertilization of cucumbers in solar greenhouse. In this study, the main cultivated crop in solar greenhouse, cucumber was used as the experimental material, and four fertilization levels were set to conduct pot experiments. The aim was to (1) elucidate the effect of excessive chemical fertilizer application on cucumber rhizosphere soil microbial PLFAs and metabolites in solar greenhouse, and to (2) reveal the interaction between microbial communities and metabolites under different fertilization levels.

Materials and methods

Soil preparation and sampling

This experiment was carried out in a solar greenhouse (41°310 N–123°240E) of Shenyang Agricultural University, China. The cucumber variety was “Jinyou 30,” and the soil used in the experiment was garden soil that has never been planted with cucurbit vegetables. The soil chemical properties are as follow: pH = 7.21, soil organic matter (SOM) = 17.8 g kg−1, available N (AN) = 105 mg kg−1, available P (AP) = 112.8 mg kg−1, and available K (AK) = 121.5 mg kg−1, were determined by the Bao’s Method (Bao, 2000). According to the national survey and test results of vegetable fertilization (Li et al., 2019). Four fertilization levels were set in the pot experiment (each pot contained 14 kg air-dried soil), including no-fertilization (N0P0K0), normal fertilization (N1P1K1), slight excessive fertilization (N2P2K2), and extreme excessive fertilization (N3P3K3). The four treatments were implemented in a randomized design with five replications. The specific fertilization amounts were shown in Table 1. Except for the fertilization levels, other cultivation and management measures were the same.

TABLE 1
www.frontiersin.org

Table 1. The amount of fertilizer used in the experiment treatments.

The first crop was planted in September 2017, the second in March 2018, and the third in September 2018. In this experiment, samples were collected from the five replicates for each treatment after the third crop. Leaf, stem, and root of cucumber were separated and dried in an oven at 65°C for 48 h until reaching a constant weight for measuring dry weights. The fourth leaf of cucumber was collected for determining photosynthetic pigments according to Lichtenrhaler’s method (Lichtenthaler and Welbum, 1983). Yield was the total fruit fresh weight per plant during cucumber growth. In this study, we took out the entire root system from the pot, after removing the bulk soil around the root system, we use a brush to gently collect the soil attached to the cucumber root. After picking out the plant residues, the soil was sieved through a 2 mm mesh and divided into three subsamples. The first subsample was freeze-dried and stored at −80°C for microbial PLFAs and metabolites measured, the second subsample was placed at 4°C to measure soil enzymes activities, and the third subsamples was air-dried to determine the soil chemical properties.

Soil chemical properties and enzymes activities determination

Soil pH (soil:H2O 1:2.5) was measured with pH meter (sartorius, Germany). SOM was determined by potassium dichromate volumetric method. AP was measured by molybdenum blue spectrophotometry method. AK was determined by atomic adsorption spectrophotometer. The total N (TN) was measured by Kjeldahl method. After the soil digestion with HNO3, HCl and HF, total phosphorus (TP) was determined with an ultraviolet spectrophotometer and total potassium (TK) was determined with atomic absorption flame spectrophotometer (Bao, 2000). NH4+-N and NO3-N were extracted with 1 mol·L−1 KCl and measured using a continuous colorimetric flow system (Skalar SAN++ system, Netherlands) (Zhou et al., 2010). Urease was measured by indophenol blue colorimetric method. Polyphenol oxidase was determined by iodine titration and neutral phosphatase activities were determined by phenyl disodium phosphate colorimetry (Guan, 1986; Zhang et al., 2011).

PLFA analysis

The microbial PLFAs were measured according to previous methods (Bossio, 1998; Briar et al., 2011). Briefly, lipids were extracted from 8 g of freeze-dried soil using a mixed extract (phosphate buffer:chloroform:methanol mixture, 0.8:1:2, v/v/v), and separated PLFAs from neutral lipids and glycolipids fatty acids using a SPE silica column (Supelco Inc, Bellefonte, PA). Using 19:0 methyl ester as internal standard. PLFAs were measured using gas chromatograph (Agilent 6850 Series, USA) equipped with the Sherlock microbial identification system (MIDI Inc., Newark, USA). Specific PLFA fingers could represent specific microbial communities, and the classification results of PLFAs were shown in Supplementary Table S1 (McKinley et al., 2005; Aciego Pietri and Brookes, 2009; Bach et al., 2010).

Soil metabolite profiling analysis

The soil metabolite extraction method was modified according to Song et al. (2020). Briefly, 1 g freeze-dried soil was placed in 5 ml Eppendorf tube, added 1 ml methanol:H2O = 3:1 (v/v), 1 ml ethyl acetate, and 10 μl 0.5 mg ml−1 adonitol dissolved in deionized water (internal standard) sequentially. Samples were homogenized at 45 Hz for 4 min and sonicated in ice water for 5 min. Then centrifuged at 12000r for 15 min at 4°C. All the supernatants were dried completely in vacuo at 28°C, and then dissolved in 60 μl of methoxylamine hydrochloride (15 mg ml−1 in pyridine) followed by 30 min of incubation at 80°C. Then, 80 μl of BSTFA (with 1% TMCS) was added to each sample and incubated for 90 min at 70°C. Finally, The derivatized sample was analyzed using GC–MS (7890A-5795C, Agilent, Palo Alto, Calif., USA) system equipped with Agilent HP-5MS column (30 m × 250 μm × 0.25 μm). The injection volume was 1 μl, the helium gas flow rate was 1 ml min−1. The injector temperature and ion source temperature were at 250 and 230°C, respectively. The oven program was as follows: the initial temperature was kept at 50°C for 2 min, raised to 120°C at a rate of 5°C min−1, then hold for 2 min at 120°C, increased to 180°C at a rate of 6°C min−1 and hold for 2 min, and finally increased to 250°C at a rate of 4°C min−1 and hold for 12 min.

The quality control (QC) sample was prepared by pooling 10 μl aliquots of all samples to check the stability of the instrument throughout the analysis. Compound identification was performed by mass comparison with NIST 08 (http://www.nist.gov) mass spectrometry library. The detailed classification information of metabolites was shown in Supplementary Table S2.

Statistical analysis

Significant differences in plant growth parameters, soil chemical properties, soil enzymes activities, soil microbial PLFAs and relative abundance of classified soil metabolites among fertilization treatments were examined using one-way analysis of variance (ANOVA) and mean values were compared with Duncan’s multiple range test (p < 0.05) using the SPSS version 19.0 software package (SPSS, Inc., Chicago, IL, United States). The histograms were drawn using origin 2019b. The partial least-square discriminant analysis (PLS-DA) based on metabolomics data were performed using SMICA-P software (version 14.0). The differential metabolites were selected by a variable importance (VIP) > 1 and p < 0.05 (Jia et al., 2019). Using vegan (Oksanen et al., 2013) and ggplot2 (Wickham, 2011) packages to draw redundancy analysis (RDA) and heatmap based on microbial PLFAs data in R v.4.1.0. The correlation and visualization between microbial PLFAs and differential metabolites were performed using R and Cytoscape 3.7.2 software.

Results

Plant growth

Different fertilization levels had an obvious effects on cucumber plant morphology (Figure 1A). Compared with no-fertilization, fertilization led to a significant increase in cucumber yield, and the yield of N1P1K1 and N2P2K2 were significantly higher than that of N3P3K3 (Figure 1B). Fertilization significantly increased the chlorophyll contents in leaves, and the chlorophyll a under N2P2K2 treatment and the chlorophyll b under the N2P2K2 and N3P3K3 treatments were significantly higher than those under the other treatments (Figure 1C). In addition, the plants dry weights under fertilization treatments were significantly higher than those under no-fertilization treatment. However, from N1P1K1 to N3P3K3 treatments, leaf and root dry weights decreased significantly overall (Figure 1D).

FIGURE 1
www.frontiersin.org

Figure 1. The effects of different fertilization levels on (A) plant images, (B) cucumber yields, (C) leaf photosynthetic pigments, and (D) plant dry weights. Values are means ± SEM (n = 5). Different letters above the bars indicate statistically significant differences at p < 0.05.

Soil chemical properties and enzymes activities

Soil pH decreased significantly, and contents of NH4+-N, NO3-N, AP, AK, TN, TP and TK significantly increased with the increase of chemical fertilization. In addition, the SOM contents under N3P3K3 treatment were significantly higher than those under the N0P0K0 and N1P1K1 treatments (Table 2).

TABLE 2
www.frontiersin.org

Table 2. The effects of different fertilization levels on soil chemical properties.

The effect of different fertilization levels on soil enzyme activities were shown in Figure 2. The activity of polyphenol oxidase decreased first and then increased with the fertilization rates. Among the four treatments, the activity of polyphenol oxidase in N1P1K1 was the lowest and that in N2P2K2 was the highest (Figure 2A). With the increase of fertilization, the activities of urease and neutral phosphatase both increased first and then decreased significantly, and the highest points of enzyme activities were N2P2K2 and N1P1K1, respectively (Figures 2B,C).

FIGURE 2
www.frontiersin.org

Figure 2. The effects of different fertilization levels on soil enzymes activities. Values are means ± SEM (n = 5). Different letters above the bars indicate statistically significant differences at p < 0.05.

Soil microbial community structure

In this study, total PLFAs first increased and then decreased with fertilization rate, and among the four fertilization treatments, total PLFA was significantly highest under N1P1K1 and lowest under N3P3K3 (Table 3). The percentage of bacteria was significantly higher under N3P3K3 treatment than other treatments. The percentages of fungi and arbuscular mycorrhizal fungi (AMF) were significantly lower under N3P3K3 treatment than other treatments. In addition, the percentages of actinomycetes under N1P1K1 and N2P2K2 treatments were significantly higher than those under N0P0K0 and N3P3K3 treatments (Table 3). A total of 22 PLFAs were detected in rhizosphere soils under different fertilization gradients (Figure 3A). Amongst them, the relative abundance Gram-negative bacteria (16:0, 18:1ω7c, 18:1ω5c and 16:1ω7c), Gram-positive bacteria (i15:0, i17:2ω9c, cy19:0, a15:0 and i16: 0), actinomycetes (10Me16:0), and fungi (18:1ω9c, 16:1ω5c) were all more than 3%, and they accounted for 82.67–83.17% of the total PLFAs. Notably, the relative abundances of 16:1ω5c, 18:1ω7c and 18:1ω9c decreased with increasing fertilizer application (Figure 3A).

TABLE 3
www.frontiersin.org

Table 3. The effects of different fertilization levels on soil total PLFA and percentage of microbial PLFAs.

FIGURE 3
www.frontiersin.org

Figure 3. The effects of different fertilization levels on soil microbial communities. (A) Heatmap of relative abundance of microbial PLFAs, reddish color indicates increased metabolites content while bluish color indicates decreased metabolites content. (B) Redundancy analysis (RDA) of the rhizosphere microbial PLFAs with soil chemical properties, Significant variables via forward selection are labeled with asterisk (**, and * represent p < 0.01, and 0.05, respectively).

RDA showed that no-fertilization and fertilization treatments were clearly distinguished along the first axis with described 55.75% of variations in the microbial PLFAs. Furthermore, the microbial PLFAs of N1P1K1 and N2P2K2 treatments were clustered together, and along the second axis which described 14.61% of variations significantly distinguished with the N3P3K3 treatment (Figure 3B). With the increment of fertilization, the pH (p = 0.001) decreased significantly, the NH4+-N (p = 0.0464) and AP (p = 0.016) increased significantly, which had the obviously effects on the changes of microbial communities (Figure 3B).

Soil metabolites

A total of 69 metabolites were detected and identified in all treatments, the specific information of soil metabolites listed in Supplementary Table S2. The PLS-DA score plot showed the soil metabolites under four fertilization levels were clearly distinguished, the first and second axis, respectively, described 17.8% and 14.1% of variations in the metabolites (Figure 4), Metabolites were divided into 10 categories according to their molecular structures (Figure 5A). Soil sugars and esters contents were significantly lowest under N3P3K3 treatment than other treatments. The content of organic acids of N2P2K2 treatment was significantly higher than those of N0P0K0, N1P1K1 and N3P3K3 treatments. The content of alkanes under N2P2K2 treatment was significantly higher than those under N0P0K0 and N1P1K1 treatments. In addition, phenols and nitriles contents were significantly higher under fertilization treatments than N0P0K0 treatment, and others content were significantly higher under N3P3K3 treatment compared to N0P0K0 treatment (Figure 5A).

FIGURE 4
www.frontiersin.org

Figure 4. The partial least squares-discriminant analysis (PLS-DA) of soil metabolites under different fertilization levels.

FIGURE 5
www.frontiersin.org

Figure 5. (A) Relative abundance of classified metabolites under different fertilization levels. Values are means ± SEM (n = 5). Different letters above the bars indicate statistically significant differences at p < 0.05. (B) Heatmap of differential metabolites under different fertilization levels, reddish color indicates increased metabolites content while bluish color indicates decreased metabolites content, * represent differential metabolites for three fertilization treatments compared to no-fertilization (VIP > 1 and p < 0.05).

The relative abundance of differential metabolites among four fertilization treatments were shown in heatmap (Figure 5B). In our study, the main differential metabolites were organic acids, esters, sugars, etc. Specifically, the relative abundances of formamide, N,N-diethyl-, octanoic acid, decane,4-methyl- and dibutyl phthalate were gradually increased, while the relative abundances of D-trehalose, 2-palmitoylglycerol and 2-methyloctacosane were gradually down-regulated with fertilization addition rates (Figure 5B).

Correlations between the soil metabolism and microbial PLFAs

To further elucidate the relationship between microbial PLFAs and metabolites under different fertilization levels, a correlation network was constructed (Figure 6). There were more negative correlations than positive correlations between soil microbial PLFAs and differential metabolites (22 positive correlations and 27 negative correlations). The most relevant metabolites with microbial PLFAs were oleamide, phosphoric acid and decane,3,8-dimethyl. The microbial PLFAs closely related to metabolites were 18:1ω9c, 18:1ω5c, 16:1ω5c and i17:2ω9c. Specifically, 18:1ω9c was positively correlated with decane,3,8-dimethyl-, 11-octadecenoic acid, benzaldehyde,3,5-dimethyl- and dibutyl phthalate, and negatively correlated with heptadecanenitrile, O-toluic acid, octanoic acid, decane,4-methyl-, acetamide,N,N-diethyl-, oleamide and hexanoic acid. 18:1ω5c was positively correlated with heptadecanenitrile, O-toluic acid and octanoic acid, and negatively correlated with oleamide, hexadecane,7,9-dimethyl- and decane,3,8-dimethyl-. i17:2ω9c was positively correlated with 11-octadecenoic acid and heneicosane, and negatively correlated with formamide,N,N-diethyl-, dibutyl phthalate, oleamide and decane,3,8-dimethyl-. And 16:1ω5c was positively correlated with 2-palmitoylglycerol, hentriacontane, 11-octadecenoic acid, decane,4-methyl- and D-trehalose, and negatively correlated with 9-octadecenoic acid.

FIGURE 6
www.frontiersin.org

Figure 6. Co-occurrence network of the microbial PLFAs and differential metabolites under different fertilization levels. A connection indicates a strong correlation (spearman correlation analysis, p < 0.05), the stronger the spearman correlation, the thicker the line. The size of each node is proportional to the number of connections, the blue and red lines indicate negative and positive relationships, respectively.

Discussion

Appropriate nutrient inputs in crop production systems can help to increase crop yields. It has been reported that at least 30 to 50% of crop yield can be attributed to commercial fertilizer nutrient inputs (Stewart et al., 2005). In this study, the yields of N1P1K1 and N2P2K2 were significantly increased by 5.69- and 5.40-fold compared to N0P0K0, respectively (Figure 1). Chemical fertilizers are frequently applied in much bigger amounts than crops are actually consumed to improve crop yields. This resulted in excess of N, P and K in the soil, which could cause soil obstacles and affect crop productivity (Dong and Lin, 2020). The results of this study are consistent with this, the cucumber yield of N3P3K3 was much lower than those of N1P1K1 and N2P2K2 (Figure 1).

Excessive fertilization led to soil nutrient accumulation and soil acidification (Table 2). In addition, polyphenol oxidase activity of N2P2K2 was significantly higher than other fertilization treatments (Figure 2A). Polyphenol oxidase is common redox enzyme related to decomposing toxic substances in soil, and plays a role in soil environmental remediation (Fioretto et al., 2000; Borowik et al., 2014; Grosso et al., 2016). It was reported that in addition to adversely affecting the soil, excessive fertilization could stimulate some self-healing abilities of soil, but this ability decreased with the severity of excessive fertilization. Urease and phosphatase were common indicators that were used to evaluate the ability of microorganisms to utilize nitrogen and phosphorus (Adetunji et al., 2017). The urease and neutral phosphatase activities under fertilization treatments were significantly higher than under N0P0K0 treatment (Figures 2B,C), this might be caused by the low soil nutrient content of N0P0K0 treatment. However, the activities of urease and neutral phosphatase of N3P3K3 were significantly lower than those of N2P2K2 and N1P1K1 treatments, indicating that the soil nutrient utilization ability would decrease under excessive fertilization.

Fertilizer addition had a significant effect on rhizosphere microbial PLFAs (Figure 3A). Total PLFAs of N3P3K3 were significantly lower than other treatments (Table 3). This might be related to soil acidification and nitrate accumulation caused by excessive fertilization (Table 2; Figure 3B). It was consistent with previous studies that soil acidification and high salinity stress could negatively impact the composition and function of soil microbial communities, leading to the attenuation of microbial biomass (He et al., 2007; Kamble et al., 2014; Wei et al., 2017; Becker et al., 2020). The majority of bacteria were better suited to nutrient-rich and highly available carbon sources, while fungus appeared to be more capable of utilizing stubborn carbon sources (He et al., 2007; Grosso et al., 2016). Accordingly, as fertilizer input increased, the percentage of bacteria increased, while the percentage of fungi decreased (Table 3). In addition, the percentage of AMF was significantly lower under N3P3K3 treatment than other treatments (Table 3). AMF is a ubiquitous soil microorganism that can form a mutualistic symbiotic relationship with plant roots, transferring immobile phosphorus resources from great distances to the vicinity of roots (Lin et al., 2012). However, in the case of nutrient sufficiency or even overuse, plants do not need to rely on AMF symbiosis to obtain nutrients, so the relative abundance of AMF decreases with the amount of fertilization (Table 3). In the study, the relative abundances of 16:1ω5c, 18:1ω7c, 18:1ω9c decreased with increasing fertilization rates, suggesting that they could serve as main microbial markers sensitive to excessive fertilization (Figure 3A). Similarly, 18:1ω7c (gram-negative bacteria) has been reported to respond negatively to nitrogen addition (Liu et al., 2015).

Soil metabolites are substances produced by plant photosynthesis, secreted by roots into the rhizosphere soil, and transformed by soil microbial communities. Metabolites play important roles in regulating both plant growth and soil microbial communities (Son et al., 2016; Lu et al., 2022). Previous studies have shown that different fertilization levels have significant effects on soil rhizosphere metabolite profiles (Figure 5). Fertilization affects the root morphology and root exudates of grassland plants to improve plant fitness (Lemanski and Scheu, 2014). Phenols are generally considered to be highly active autotoxins. It has been reported that the accumulation of phenols in soil can inhibit the growth and development of plants, which are among the key reasons for continuous cropping obstacles (Deng et al., 2017; Wu et al., 2017). In the study, soil phenol contents were higher under fertilization treatments than that under no-fertilization treatment (Figure 5A), elucidating that fertilizer inputs might accelerate the progression of continuous cropping obstacles. In addition, compared with soil organic acids and alkanes of N0P0K0, N2P2K2 was significantly increased, and N3P3K3 was not significantly different (Figure 5A). Previous report suggested that the secretion of organic acids might be a strategy for plants to recruit beneficial microorganisms in the root zone in response to high N input (Chen et al., 2019). Cucumbers tend to increase disease resistance by secreting more organic acids to aggregate beneficial microorganisms (Wen et al., 2020). Alkanes could fight microbial pathogens and inhibit the reproduction of pathogenic bacteria (Lu et al., 2022). Plants adapt to biotic and abiotic stresses by altering their root-secreted chemicals to aggregate health-promoting microbiomes. This so-called “cry-for-help” hypothesis provides an explanation for the feedback response of rhizosphere soils to stress (Rolfe et al., 2019). Therefore, it was speculated that organic acids and alkanes might play a key role in recruiting rhizosphere growth-promoting bacteria (PGPR) and affecting the microbial communities to alleviate the stress of slight excessive fertilization (N2P2K2). However, in extreme excessive fertilization (N3P3K3), the self-regulation of plants and soils might no longer be able to cope with high salinity stress. Because of high fertilizer input, the underdeveloped root system of crops communicates less efficiently with soil microbial communities. And in the case of poor soil quality, beneficial metabolites in the rhizosphere soil are more likely to be hijacked by parasitic microorganisms and arthropods (Rolfe et al., 2019). In addition, it was worth noting that the content of dibutyl phthalate (DBP) was up-regulated and the content of D-trehalose was down-regulated with fertilizer addition levels (Figure 5B). DBP might threaten the stability of soil microbial communities and functions, sustainable development of agriculture and human health (Xu et al., 2008). Extensive application of chemical fertilizers and pesticides in agriculture might lead to the accumulation of DBP in soil (He et al., 2015). Trehalose as a protective agent could positively affect plant growth and overcome the adverse effects of salt stress (Jain and Roy, 2009; Sadak et al., 2019). However, with the increase of fertilizer input, trehalose content also decreased, which might be related to the decrease in the abundance of nitrogen-fixing microorganisms capable of synthesize trehalose with higher fertilizer application (Garg and Chandel, 2011).

The correlation between rhizosphere microorganisms and metabolites is very important for the homeostasis of soil environment. In the study, it was found that there were more negative correlations than positive correlations between microorganisms and differential metabolites under different fertilization levels (Figure 6). Because rhizosphere metabolites are the most easily consumed and utilized carbon sources by microorganisms (Yuan et al., 2016). In this study, 18:1ω9c, 18:1ω5c, 16:1ω5c and i17:2ω9c in rhizosphere microbial PLFAs were more closely related to metabolites under different fertilization levels (Figure 6). Notably, 16:1ω5c (AMF) is a beneficial microorganism that forms a symbiotic relationship with plant roots (Lin et al., 2012). Its relative abundance decreased significantly with increasing fertilizer inputs (Figure 3A). In this study, AMF was positively correlated with esters (2-palmitoylglycerol), alkanes (hentriacontane, decane,4-methyl-), organic acids (11-octadecenoic acid) and sugars (D-trehalose), and negatively correlated with organic acids (9-octadecenoic acid) (Figure 6). It has been reported that plant esters are most likely transferred from host plants to AMF in the form of 2-palmitoylglycerol (Kameoka et al., 2019), so 2-palmitoylglycerol was significantly positively correlated with AMF. Alkanes was known to be strong resistant to plant pathogens (Lu et al., 2022). 11-Octadecenoic acid has broad-spectrum antibacterial activity, showing dose-dependent antioxidant activity (Alqahtani et al., 2019). And sugars were not only carbon source that utilized by microorganisms, but also modulates the chemotaxis of rhizosphere bacteria (Eilers et al., 2010; Jin et al., 2019). In addition, 9-octadecenoic acid, which was significantly negatively correlated with AMF, had been shown to be autotoxin, affecting germination and seedling growth in wheat and Sinapis arvensis (Tahir et al., 2018). This indicated that with the increase of fertilization amount, the contents of beneficial microorganisms and metabolites in the rhizosphere soil decreased. Excessive fertilization had an adverse effect on the rhizosphere soil environment.

Conclusion

Excessive fertilization led to enrichment of soil nutrients and soil pH decreased. Compared with other treatments, total PLFAs content, arbuscular mycorrhizal fungi (AMF) content and fungi percentage were significantly decreased, while percentage of bacteria was significantly increased under N3P3K3. Compared with soil organic acids and alkanes of N0P0K0, N2P2K2 was significantly increased, and N3P3K3 was not significantly different. In addition, the yield of N1P1K1 and N2P2K2 were significantly higher than that of N3P3K3. This suggested that the fertilization of N1P1K1 was more reasonable, as N2P2K2 did not significantly improve yield. The non-decreased yield of N2P2K2 might be due to the fact that cucumber maintained normal growth and yield by secreting beneficial metabolites (organic acids and alkanes) to alleviate the adverse effects on microorganisms of slight excessive fertilization. The significant decrease in N3P3K3 yield was closely related to the decrease of beneficial microorganisms and metabolites in the rhizosphere soil with the amount of fertilization. These results not only clarified the effect of excessive NPK fertilization on soil microorganisms and metabolites, but also provided a theoretical basis for guiding the rational use of chemical fertilizers and promoting the sustainable development of cucumber production in solar greenhouse.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

SW and TB: data curation. ZS and HF: funding acquisition. SW, TB, TW, and YZ: experiments. SW: writing-original draft. ZS, HF, and MA: writing-review and editing. All authors contributed to the article and approved the submitted version.

Funding

This study was financially supported by the National Key R&D Program of China (2020YFD1000301),China Agriculture Research System (CARS-23-D01), the National Natural Science Foundation of China (31902093) and Shenyang Science and technology project (21109308).

Acknowledgments

The authors would like to thank Muhammad Sohail Khan for his contributions in revising the manuscript.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1004836/full#supplementary-material

References

Aciego Pietri, J. C., and Brookes, P. C. (2009). Substrate inputs and pH as factors controlling microbial biomass, activity and community structure in an arable soil. Soil Biol. Biochem. 41, 1396–1405. doi: 10.1016/j.soilbio.2009.03.017

CrossRef Full Text | Google Scholar

Adetunji, A. T., Lewu, F. B., Mulidzi, R., and Ncube, B. (2017). The biological activities of beta-glucosidase, phosphatase and urease as soil quality indicators: a review. J. Soil Sci. Plant Nutr. 17, 794–807. doi: 10.4067/S0718-95162017000300018

CrossRef Full Text | Google Scholar

Alqahtani, F. Y., Aleanizy, F. S., Mahmoud, A. Z., Farshori, N. N., Alfaraj, R., Al-Sheddi, E. S., et al. (2019). Chemical composition and antimicrobial, antioxidant, and anti-inflammatory activities of Lepidium sativum seed oil. Saudi J. Biol. Sci. 26, 1089–1092. doi: 10.1016/j.sjbs.2018.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Bach, E. M., Baer, S. G., Meyer, C. K., and Six, J. (2010). Soil texture affects soil microbial and structural recovery during grassland restoration. Soil Biol. Biochem. 42, 2182–2191. doi: 10.1016/j.soilbio.2010.08.014

CrossRef Full Text | Google Scholar

Badri, D. V., and Vivanco, J. M. (2009). Regulation and function of root exudates. Plant Cell Environ. 32, 666–681. doi: 10.1111/j.1365-3040.2008.01926.x

CrossRef Full Text | Google Scholar

Bao, S. D. (2000). Soil Agro-Chemical Analyze. China Agricultural Press: Beijing

Google Scholar

Becker, F., Mac Laren, C., Brink, C. J., Jacobs, K., Roux, M. R., and Swanepoel, P. A. (2020). High nitrogen rates do not increase canola yield and may affect soil bacterial functioning. Agron. J. 112, 523–536. doi: 10.2134/agronj2019.05.0359

CrossRef Full Text | Google Scholar

Bindraban, P. S., Dimkpa, C. O., and Pandey, R. (2020). Exploring phosphorus fertilizers and fertilization strategies for improved human and environmental health. Biol. Fertil. Soils 56, 299–317. doi: 10.1007/s00374-019-01430-2

CrossRef Full Text | Google Scholar

Borowik, A., Wyszkowska, J., Kucharski, M., and Kucharski, J. (2014). Resistance of dehydrogenases, catalase, urease and plants to soil contamination with zinc. J. Elem. 19, 929–946. doi: 10.5601/jelem.2013.18.4.566

CrossRef Full Text | Google Scholar

Bossio, S. (1998). Impacts of carbon and flooding on soil microbial communities: phospholipid fatty acid profiles and substrate utilization patterns. Microb. Ecol. 35, 265–278. doi: 10.1007/s002489900082

PubMed Abstract | CrossRef Full Text | Google Scholar

Briar, S. S., Fonte, S. J., Park, I., Six, J., Scow, K., and Ferris, H. (2011). The distribution of nematodes and soil microbial communities across soil aggregate fractions and farm management systems. Soil Biol. Biochem. 43, 905–914. doi: 10.1016/j.soilbio.2010.12.017

CrossRef Full Text | Google Scholar

Chen, S., Waghmode, T. R., Sun, R., Kuramae, E. E., Hu, C., and Liu, B. (2019). Root-associated microbiomes of wheat under the combined effect of plant development and nitrogen fertilization. Microbiome 7:136. doi: 10.1186/s40168-019-0750-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, H., Yuan, M., Tang, L., Shen, Y., Yu, Q., and Li, S. (2022). Integrated microbiology and metabolomics analysis reveal responses of soil microorganisms and metabolic functions to phosphorus fertilizer on semiarid farm. Sci. Total Environ. 817:152878. doi: 10.1016/j.scitotenv.2021.152878

PubMed Abstract | CrossRef Full Text | Google Scholar

Chynoweth, R., Jimenez, D., Liberti, D., Bellon-Dona, D., Carralero-González, A., Crespo-Sempere, A., et al. (2021). First report of cucurbit chlorotic yellows virus infecting cucumber plants in Spain. Plant Dis. 105:2258. doi: 10.1094/PDIS-12-20-2553-PDN

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, J. J., Zhang, Y. L., Hu, J. W., Jiao, J. G., Hu, F., Li, H. X., et al. (2017). Autotoxicity of phthalate esters in tobacco root exudates: effects on seed germination and seedling growth. Pedosphere 27, 1073–1082. doi: 10.1016/S1002-0160(17)60374-6

CrossRef Full Text | Google Scholar

Dong, L., Xu, J., Li, Y., Fang, H., Niu, W., Li, X., et al. (2018). Manipulation of microbial community in the rhizosphere alleviates the replanting issues in Panax ginseng. Soil Biol. Biochem. 125, 64–74. doi: 10.1016/j.soilbio.2018.06.028

CrossRef Full Text | Google Scholar

Dong, N.-Q., and Lin, H.-X. (2020). Higher yield with less nitrogen fertilizer. Nat. Plants 6, 1078–1079. doi: 10.1038/s41477-020-00763-3

CrossRef Full Text | Google Scholar

Eilers, K. G., Lauber, C. L., Knight, R., and Fierer, N. (2010). Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biol. Biochem. 42, 896–903. doi: 10.1016/j.soilbio.2010.02.003

CrossRef Full Text | Google Scholar

Fioretto, A., Papa, S., Curcio, E., Sorrentino, G., and Fuggi, A. (2000). Enzyme dynamics on decomposing leaf litter of Cistus incanus and Myrtus communis in a Mediterranean ecosystem. Soil Biol. Biochem. 32, 1847–1855. doi: 10.1016/S0038-0717(00)00158-9

CrossRef Full Text | Google Scholar

Garbeva, P., Veen, J. A., and Elsas, J. D. (2004). Microbial diversity in soil: selection of microbial populations by plant and soil type and implications for disease suppressiveness. Annu. Rev. Phytopathol. 42, 243–270. doi: 10.1146/annurev.phyto.42.012604.135455

PubMed Abstract | CrossRef Full Text | Google Scholar

Garg, N., and Chandel, S. (2011). The effects of salinity on nitrogen fixation and Trehalose metabolism in Mycorrhizal Cajanus cajan (L.) Millsp plants. J. Plant Growth Regul. 30, 490–503. doi: 10.1007/s00344-011-9211-2

CrossRef Full Text | Google Scholar

Grosso, F., Bååth, E., and De Nicola, F. (2016). Bacterial and fungal growth on different plant litter in Mediterranean soils: effects of C/N ratio and soil pH. Appl. Soil Ecol. 108, 1–7. doi: 10.1016/j.apsoil.2016.07.020

CrossRef Full Text | Google Scholar

Guan, S. Y. (1986). Study Way of Soil Enzymes. Agriculture Press: Beijing

Google Scholar

Hao, Z., Chen, B., and Li, X. (2020). Relationship between soil chemical properties and microbial metabolic patterns in intensive greenhouse tomato production systems. Arch. Agron. Soil Sci. 66, 1334–1343. doi: 10.1080/03650340.2019.1666369

CrossRef Full Text | Google Scholar

He, J. Z., Shen, J.-P., Zhang, L.-M., Zhu, Y.-G., Zheng, Y.-M., Xu, M.-G., et al. (2007). Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea of a Chinese upland red soil under long-term fertilization practices. Environ. Microbiol. 9, 2364–2374. doi: 10.1111/j.1462-2920.2007.01358.x

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Gielen, G., Bolan, N. S., Zhang, X., Qin, H., Huang, H., et al. (2015). Contamination and remediation of phthalic acid esters in agricultural soils in China: a review. Agron. Sustain. Dev. 35, 519–534. doi: 10.1007/s13593-014-0270-1

CrossRef Full Text | Google Scholar

Islam, M. R., Chauhan, P. S., Kim, Y., Kim, M., and Sa, T. (2011). Community level functional diversity and enzyme activities in paddy soils under different long-term fertilizer management practices. Biol. Fertil. Soils 47, 599–604. doi: 10.1007/s00374-010-0524-2

CrossRef Full Text | Google Scholar

Jain, N. K., and Roy, I. (2009). Effect of trehalose on protein structure. Protein Sci. 18:36. doi: 10.1002/pro.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, G., Sha, K., Feng, X., and Liu, H. (2019). Post-thawing metabolite profile and amino acid oxidation of thawed pork tenderloin by HVEF-A short communication. Food Chem. 291, 16–21. doi: 10.1016/j.foodchem.2019.03.154

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Y., Zhu, H., Luo, S., Yang, W., Zhang, L., Li, S., et al. (2019). Role of maize root exudates in promotion of colonization of bacillus velezensis strain S3-1 in Rhizosphere soil and root tissue. Curr. Microbiol. 76, 855–862. doi: 10.1007/s00284-019-01699-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, O. A. H., Sdepanian, S., Lofts, S., Svendsen, C., Spurgeon, D. J., Maguire, M. L., et al. (2014). Metabolomic analysis of soil communities can be used for pollution assessment. Environ. Toxicol. Chem. 33, 61–64. doi: 10.1002/etc.2418

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamble, P. N., Gaikwad, V. B., Kuchekar, S. R., and Bååth, E. (2014). Microbial growth, biomass, community structure and nutrient limitation in high pH and salinity soils from Pravaranagar (India). Eur. J. Soil Biol. 65, 87–95. doi: 10.1016/j.ejsobi.2014.10.005

CrossRef Full Text | Google Scholar

Kameoka, H., Tsutsui, I., Saito, K., Kikuchi, Y., Handa, Y., Ezawa, T., et al. (2019). Stimulation of asymbiotic sporulation in arbuscular mycorrhizal fungi by fatty acids. Nat. Microbiol. 4, 1654–1660. doi: 10.1038/s41564-019-0485-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemanski, K., and Scheu, S. (2014). Fertilizer addition lessens the flux of microbial carbon to higher trophic levels in soil food webs of grassland. Oecologia 176, 487–496. doi: 10.1007/s00442-014-3037-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B. B., Roley, S. S., Duncan, D. S., Guo, J. R., Quensen, J. F., Yu, H. Q., et al. (2021). Long-term excess nitrogen fertilizer increases sensitivity of soil microbial community to seasonal change revealed by ecological network and metagenome analyses. Soil Biol. Biochem. 160:108349. doi: 10.1016/j.soilbio.2021.108349

CrossRef Full Text | Google Scholar

Li, J. G., Wan, X., Liu, X. X., Chen, Y., Slaughter, L. C., Weindorf, D. C., et al. (2019). Changes in soil physical and chemical characteristics in intensively cultivated greenhouse vegetable fields in North China. Soil Tillage Res. 195:104366. doi: 10.1016/j.still.2019.104366

CrossRef Full Text | Google Scholar

Lichtenthaler, H. K., and Welburn, A. R. (1983). Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents. Biochem. Soc. Trans. 11, 591–592. doi: 10.1042/bst0110591

CrossRef Full Text | Google Scholar

Lin, X., Feng, Y., Zhang, H., Chen, R., Wang, J., Zhang, J., et al. (2012). Long-term balanced fertilization decreases Arbuscular Mycorrhizal fungal diversity in an arable soil in North China revealed by 454 pyrosequencing. Environ. Sci. Technol. 46, 5764–5771. doi: 10.1021/es3001695

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Gundersen, P., Zhang, W., Zhang, T., Chen, H., and Mo, J. (2015). Effects of nitrogen and phosphorus additions on soil microbial biomass and community structure in two reforested tropical forests. Sci. Rep. 5, 1–10. doi: 10.1038/srep14378

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, S., He, Y. H., Chen, Y. Q., Chen, L. J., Wang, Z. Y., Yuan, J., et al. (2022). Co-analysis of rhizosphere metabolomics and bacterial community structures to unfold soil ecosystem health in Camellia oleifera land under long-term cultivation. Appl. Soil Ecol. 171:104336. doi: 10.1016/j.apsoil.2021.104336

CrossRef Full Text | Google Scholar

McKinley, V. L., Peacock, A. D., and White, D. C. (2005). Microbial community PLFA and PHB responses to ecosystem restoration in tallgrass prairie soils. Soil Biol. Biochem. 37, 1946–1958. doi: 10.1016/j.soilbio.2005.02.033

CrossRef Full Text | Google Scholar

Mohanram, S., and Kumar, P. (2019). Rhizosphere microbiome: revisiting the synergy of plant-microbe interactions. Ann. Microbiol. 69, 307–320. doi: 10.1007/s13213-019-01448-9

CrossRef Full Text | Google Scholar

Muhammad, S., Müller, T., and Joergensen, R. (2008). Relationships between soil biological and other soil properties in saline and alkaline arable soils from the Pakistani Punjab. J. Arid Environ. 72, 448–457. doi: 10.1016/j.jaridenv.2007.06.016

CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., Ohara, R., et al. (2013). Package ‘vegan’. Community Ecology Package, Version 2, 1–295

Google Scholar

Phour, M., Sehrawat, A., Sindhu, S. S., and Glick, B. R. (2020). Interkingdom signaling in plant-rhizomicrobiome interactions for sustainable agriculture. Microbiol. Res. 241:126589. doi: 10.1016/j.micres.2020.126589

PubMed Abstract | CrossRef Full Text | Google Scholar

Poeplau, C., Herrmann, A. M., and Kätterer, T. (2016). Opposing effects of nitrogen and phosphorus on soil microbial metabolism and the implications for soil carbon storage. Soil Biol. Biochem. 100, 83–91. doi: 10.1016/j.soilbio.2016.05.021

CrossRef Full Text | Google Scholar

Reinhold-Hurek, B., Bunger, W., Burbano, C. S., Sabale, M., and Hurek, T. (2015). Roots shaping their microbiome: global hotspots for microbial activity. Global hotspots for microbial activity. Annu. Rev. Phytopathol. 53, 403–424. doi: 10.1146/annurev-phyto-082712-102342

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, N., Wang, Y., Ye, Y. L., Zhao, Y. A., Huang, Y. F., Fu, W., et al. (2020). Effects of continuous nitrogen fertilizer application on the diversity and composition of Rhizosphere soil bacteria. Front. Microbiol. 11:1948. doi: 10.3389/fmicb.2020.01948

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolfe, S. A., Griffiths, J., and Ton, J. (2019). Crying out for help with root exudates: adaptive mechanisms by which stressed plants assemble health-promoting soil microbiomes. Curr. Opin. Microbiol. 49, 73–82. doi: 10.1016/j.mib.2019.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadak, M. S., El-Bassiouny, H. M. S., and Dawood, M. G. (2019). Role of trehalose on antioxidant defense system and some osmolytes of quinoa plants under water deficit. Bull. Natl. Res. Cent. 43, 1–11. doi: 10.1186/s42269-018-0039-9

CrossRef Full Text | Google Scholar

Shen, W. S., Lin, X. G., Shi, W. M., Min, J., Gao, N., Zhang, H. Y., et al. (2010). Higher rates of nitrogen fertilization decrease soil enzyme activities, microbial functional diversity and nitrification capacity in a Chinese polytunnel greenhouse vegetable land. Plant Soil 337, 137–150. doi: 10.1007/s11104-010-0511-2

CrossRef Full Text | Google Scholar

Smercina, D. N., Bowsher, A. W., Evans, S. E., Friesen, M. L., Eder, E. K., Hoyt, D. W., et al. (2021). Switchgrass Rhizosphere metabolite chemistry driven by nitrogen availability. Phytobiomes J. 5, 88–96. doi: 10.1094/Pbiomes-09-19-0055-Fi

CrossRef Full Text | Google Scholar

Son, S. Y., Kim, N. K., Lee, S., Singh, D., Kim, G. R., Lee, J. S., et al. (2016). Metabolite fingerprinting, pathway analyses, and bioactivity correlations for plant species belonging to the Cornaceae, Fabaceae, and Rosaceae families. Plant Cell Rep. 35, 1917–1931. doi: 10.1007/s00299-016-2006-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Y., Li, X. N., Yao, S., Yang, X. L., and Jiang, X. (2020). Correlations between soil metabolomics and bacterial community structures in the pepper rhizosphere under plastic greenhouse cultivation. Sci. Total Environ. 728:138439. doi: 10.1016/j.scitotenv.2020.138439

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, W., Dibb, D., Johnston, A., and Smyth, T. (2005). The contribution of commercial fertilizer nutrients to food production. Agron. J. 97, 1–6. doi: 10.2134/agronj2005.0001

CrossRef Full Text | Google Scholar

Tahir, N., Qader, K., Azeez, H., and Rashid, J. (2018). Inhibitory allelopathic effects of Moringa oleifera Lamk plant extracts on wheat and Sinapis arvensis L. Allelopath. J. 44, 53–66. doi: 10.26651/allelo.j./2018-44-1-1152

CrossRef Full Text | Google Scholar

Tang, J., Wang, J. J., Li, Z. Y., Wang, S. N., and Qu, Y. K. (2018). Effects of irrigation regime and nitrogen fertilizer management on CH4, N2O and CO2 emissions from saline-alkaline Paddy fields in Northeast China. Sustainability 10, 475. doi: 10.3390/su10020475

CrossRef Full Text | Google Scholar

Wallace, A. (1994). Soil acidification from use of too much fertilizer. Commun. Soil Sci. Plant Anal. 25, 87–92.

Google Scholar

Wang, Q. Q., Liu, L. L., Li, Y., Qin, S., Wang, C. J., Cai, A. D., et al. (2020). Long-term fertilization leads to specific PLFA finger-prints in Chinese Hapludults soil. J. Integr. Agric. 19, 1354–1362. doi: 10.1016/S2095-3119(19)62866-2

CrossRef Full Text | Google Scholar

Wang, J., Liao, L., Wang, G., Liu, H., Wu, Y., Liu, G., et al. (2022). N-induced root exudates mediate the rhizosphere fungal assembly and affect species coexistence. Sci. Total Environ. 804:150148. doi: 10.1016/j.scitotenv.2021.150148

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, M., Hu, G. Q., Wang, H., Bai, E., Lou, Y. H., Zhang, A. J., et al. (2017). 35 years of manure and chemical fertilizer application alters soil microbial community composition in a Fluvo-aquic soil in northern China. Eur. J. Soil Biol. 82, 27–34. doi: 10.1016/j.ejsobi.2017.08.002

CrossRef Full Text | Google Scholar

Wen, T., Yuan, J., He, X. M., Lin, Y., Huang, Q. W., and Shen, Q. R. (2020). Enrichment of beneficial cucumber rhizosphere microbes mediated by organic acid secretion. Hortic. Res. 7:154. doi: 10.1038/s41438-020-00380-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H. (2011). ggplot 2. Wiley Interdiscip. Rev. Comput. Stat. 3, 180–185. doi: 10.1002/wics.147

CrossRef Full Text | Google Scholar

Wu, H., Xu, J., Wang, J., Qin, X., Wu, L., Li, Z., et al. (2017). Insights into the mechanism of proliferation on the special microbes mediated by phenolic acids in the radix pseudostellariae Rhizosphere under continuous monoculture regimes. Front. Plant Sci. 8:659. doi: 10.3389/fpls.2017.00659

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, G., Li, F., and Wang, Q. (2008). Occurrence and degradation characteristics of dibutyl phthalate (DBP) and di-(2-ethylhexyl) phthalate (DEHP) in typical agricultural soils of China. Sci. Total Environ. 393, 333–340. doi: 10.1016/j.scitotenv.2008.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, B.-C., Li, Z.-Z., Liu, H., Gao, M., and Zhang, Y.-Y. (2007). Microbial biomass and activity in salt affected soils under arid conditions. Appl. Soil Ecol. 35, 319–328. doi: 10.1016/j.apsoil.2006.07.004

CrossRef Full Text | Google Scholar

Yuan, H., Zhu, Z., Liu, S., Ge, T., Jing, H., Li, B., et al. (2016). Microbial utilization of rice root exudates: 13C labeling and PLFA composition. Biol. Fertil. Soils 52, 615–627. doi: 10.1007/s00374-016-1101-0

CrossRef Full Text | Google Scholar

Zhang, C., Liu, G., Xue, S., and Song, Z. (2011). Rhizosphere soil microbial activity under different vegetation types on the loess plateau, China. Geoderma 161, 115–125. doi: 10.1016/j.geoderma.2010.12.003

CrossRef Full Text | Google Scholar

Zhang, Q. C., Shamsi, I. H., Xu, D. T., Wang, G. H., Lin, X. Y., Jilani, G., et al. (2012). Chemical fertilizer and organic manure inputs in soil exhibit a vice versa pattern of microbial community structure. Appl. Soil Ecol. 57, 1–8. doi: 10.1016/j.apsoil.2012.02.012

CrossRef Full Text | Google Scholar

Zhang, Z., Gao, Z., Li, S., and Zhao, Y. (2017). The method of distribution network fault location based on improved Dempster-Shafer theory of evidence. In 2017 IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC), 1–5.

Google Scholar

Zhou, J.-B., Chen, Z.-J., Liu, X.-J., Zhai, B.-N., and Powlson, D. S. (2010). Nitrate accumulation in soil profiles under seasonally open ‘sunlight greenhouses’ in Northwest China and potential for leaching loss during summer fallow. Soil Use Manag. 26, 332–339. doi: 10.1111/J.1475-2743.2010.00284.X

CrossRef Full Text | Google Scholar

Keywords: excessive chemical fertilizer, rhizosphere, microbials PLFAs, metabolites, cucumber, solar greenhouse

Citation: Wang S, Bian T, Wu T, Zhang Y, Awais M, Fu H and Sun Z (2022) Co-analysis of cucumber rhizosphere metabolites and microbial PLFAs under excessive fertilization in solar greenhouse. Front. Microbiol. 13:1004836. doi: 10.3389/fmicb.2022.1004836

Received: 27 July 2022; Accepted: 20 September 2022;
Published: 06 October 2022.

Edited by:

Qin Gu, Nanjing Agricultural University, China

Reviewed by:

Fangwei Yu, Jiangsu Academy of Agricultural Sciences (JAAS), China
Yaqing Chang, Dalian Ocean University, China

Copyright © 2022 Wang, Bian, Wu, Zhang, Awais, Fu and Sun. 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: Hongdan Fu, fuhongdan@syau.edu.cn; Zhouping Sun, sunzp@syau.edu.cn

Disclaimer: 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.