- 1State Key Laboratory of Crop Biology, College of Horticulture Science and Engineering, Shandong Agricultural University, Tai’an, Shandong, China
- 2Scientific Observing and Experimental Station of Facility Agricultural Engineering (Huang-Huai-Hai Region), Ministry of Agriculture and Rural Affairs, Tai’an, Shandong, China
- 3Shandong Collaborative Innovation Center for Fruit and Vegetable Production with High Quality and Efficiency, Tai’an, Shandong, China
- 4Key Laboratory of Biology and Genetic Improvement of Horticultural Crop (Huang-Huai Region), Ministry of Agriculture and Rural Affairs, Tai’an, Shandong, China
- 5School of Economics, Qingdao University, Qingdao, Shandong, China
- 6College of Plant Protection, Shandong Agricultural University, Tai’an, Shandong, China
Water scarcity is a worldwide problem, and in order to obtain plenty of production, agricultural irrigation water accounts for a large portion. Many studies have shown that the interaction of root microorganisms and soil can promote crop growth. Developing ways to reduce irrigation to maintain soil fertility and ensure crop yield by regulating the root microenvironment is an important research goal. Here, we developed a reasonable irrigation plan for eggplant cultivation in a solar greenhouse. The maximum theoretical amount of water demand during eggplant planting obtained from a previous study was used as the control (CK), and the irrigation in the treatments was reduced by 10, 20 and 30% relative to this amount. The 10% irrigation reduction treatment (T1) significantly improved soil nutrients and increased soil catalase, urease and alkaline phosphatase activities (p < 0.05). Further analysis of rhizosphere microorganisms revealed the highest richness and diversity of the microbial community under the T1 treatment, with Bacilli as the most abundant bacteria and Aspergillaceae as the most abundant fungi and lower relative abundances of Chloroflexi and Acidobacteria (p < 0.05). Changes in microbial community structure under the influence of different irrigation treatments resulted in improvements in rhizosphere N cycling and nutrient catabolism. The plant–microbe interactions led to significant increases in eggplant plant height, root vigour, root surface area, leaf chlorophyll a, leaf net photosynthetic rate, water use efficiency, transpiration rate, and stomatal conductance under the T1 treatment compared to the CK treatment; soluble sugar, soluble protein and free amino acid contents in eggplant fruit increased by 10.8, 12.3 and 6.7%, respectively; and yield increased by 3.9%. Our research proved that the 10% irrigation reduction treatment (T1) could improve microbial community richness and fruit yield, which would improve irrigation efficiency and cost reduction in agriculture.
Introduction
Water scarcity is a major challenge in agriculture. As the population grows and agricultural production expands, the demand for agricultural irrigation water is constantly increasing. However, the anticipation of a climate change and drought has made researchers less optimistic in their predictions for future agricultural water supplies (Ault Toby, 2020; Berdugo et al., 2020). To face the challenge of producing more food in the coming decades, it is necessary to improve the water use efficiency (WUE) in agricultural irrigation (Du et al., 2015).
For crops, especially vegetables, an adequate water supply is essential, but this does not mean that more water is better. Excessive water can cause stress similar to that caused by drought and even plant death (Hirabayashi et al., 2013). Flooding causes hypoxia in plant roots, which can cause irreversible damage (Voesenek and Bailey-Serres, 2015). The amount of irrigation determines the maintenance of nutrients in the soil, which affects the growth and development of plants. Agricultural producers need more water-efficient cropping systems to cope with high water use in irrigated agriculture and high unproductive losses due to runoff and evaporation (Bodner et al., 2015). Among them, water-saving irrigation measures are crucial. Studies have found that a proper reduction in agricultural irrigation is not inevitably detrimental to crops but instead benefits crops in terms of photosynthetic metabolism (Morales et al., 2020). Alternate drip irrigation of 60% field capacity showed no significant difference of tomato yields compared to conventional surface drip irrigation (Wang et al., 2020). The alternate wetting and drying irrigation method, which reduced the total irrigation amount, improved rice grain yield compared with continuous flooding without reducing the quality of grains (Song et al., 2021).
Different irrigation regulation strategies have also been applied in agricultural production, such as regulated deficit irrigation (RDI; Pérez-Álvarez et al., 2021), partial root-zone drying (PRD; Santos et al., 2021) and alternate partial root-zone irrigation (APRI; Kang and Zhang, 2004). Studies on eggplants in greenhouse found that APRI increased the fruit yield and WUE (Du et al., 2005). The same results were found in cucumber research (Fang et al., 2015). However, most studies have focused on the relationship involving yield, quality of the eggplant and irrigation, but few studies have focused on the changes in the soil’s physical and chemical properties caused by different irrigation amounts and their effect on root microorganisms.
The interaction between plant roots and soil microorganisms affects plant growth, including root morphology, the root-to-shoot weight ratio, mineral content and uptake, and the rate of development, thus influencing crop yield and physiological processes (Rovira, 1965). Soil water content has a great impact on microbial communities and affects soil nutrients and pH, which in turn affects the growth and development of plants (Girvan Martina et al., 2003; Frey et al., 2004; Faoro et al., 2010; Rousk et al., 2010). Drought will promote an increase in specific root exudation, which affects soil microorganisms (Naylor and Coleman-Derr, 2018; Williams and de Vries, 2020). Water deficit increased the ratio of Gram-positive to Gram-negative bacteria (Fuchslueger et al., 2014, 2016; Chodak et al., 2015) and increased the proportion of actinobacteria (Hartmann et al., 2017). The reasons for these changes include differences in the substrate preference and metabolic capacity among different microorganisms, the tolerance of microorganisms to water deficit, and the activity level of specific bacteria (Naylor and Coleman-Derr, 2018). Wheat inoculated with inactivation of non-ribosomal peptide and polyketide derived metabolites P. polymyxa bacteria resulted in two times higher plant survival and three times increased biomass under severe drought stress, compared to wild type (Timmusk et al., 2015). Inoculation of Pseudomonas sp. strain GAP-P45 increased the survival, plant biomass, and root adhering soil/root tissue ratio of sunflower seedlings under drought stress (Sandhya et al., 2009). There are also studies suggesting that plant growth-promoting rhizosphere bacteria can produce plant hormones, affect cell division and differentiation, and enhance plant growth and plant pathogen defence (Verbon and Liberman, 2016). However, some studies also imply that the presence of root microorganisms accelerates the spread of soil viruses (Macfarlane, 2003). Biotic stress promotes the defensive response of crops, which affects crop growth and yield (Philippot et al., 2013; Liu et al., 2020; Saijo and Loo, 2020). Therefore, the balance between plant and soil microorganisms needs to be focused on (Bennett and Klironomos, 2019).
Eggplant (Solanum melongena L.) is widely cultivated in facilities worldwide, and its fruit yield depends on the irrigation amount and soil nutrient availability. However, few studies have been conducted on the effect of the irrigation amount on the eggplant rhizosphere microbial community and rhizosphere soil nutrient supply. In this study, we studied the influence of the irrigation amount on the soil microorganisms, rhizosphere environment, yield and quality of eggplants grown in solar greenhouse with the aim of investigating the effects of different irrigation amounts on the following topics: (1) the eggplant yield, soil characteristics, soil microbial community structure and function; (2) the correlation between microorganisms and soil environmental factors; and (3) the effects of microbial community differences on eggplant yield, with the overarching goal of determining the most suitable irrigation amount for eggplant planting in facilities. Our research will benefit agricultural irrigation and water-saving efforts, improve water use efficiency, and have a positive impact on agricultural production.
Materials and methods
Plant material and growth conditions
• Eggplant (Solanum melongena L. cv. ‘Blitha’) was used as the model plant for the assays. The basic physical and chemical properties of the tested soil in the solar greenhouse (N: 36°09′35.67′; E: 117°09′33.81′) were as follows: pH 6.94, EC 1.35 mS cm−1, and field capacity 38.8%. The available nitrogen, phosphorus and potassium contents were 186.3 mg kg−1, 157.0 mg kg−1 and 344.0 mg kg−1, respectively. Different row spacings were used for planting, with a large row spacing of 80 cm, a small row spacing of 60 cm, and a plant spacing of 50 cm. The area of each test plot was 29.4 m2 (L × W = 7 m × 4.2 m), which was replicated three times and arranged randomly. Each treatment was separated by a thick plastic film with a depth of 60 cm to ensure that there was no interference by water exchange among the treated soils. Drip irrigation was adopted, and nutrients and water were supplied by the integrated water and fertilizer solution. The fertilizer supply was N 21.44 kg hm−2 15 days, P2O5 19.29 kg hm−2 15 days, and K2O 56.81 kg hm−2 15 days in the seedling stage and N 41.45 kg hm−2 15 days, P2O5 34.30 kg hm−2 15 days, and K2O 86.83 kg hm−2 15 days in the flowering and fruit setting stage, full fruiting and fruit setting stage and harvesting stage.
• Eggplants in the autumn–winter crop were planted on September 8, 2019, and harvested on January 15, 2020. The eggplant growth period was divided into 4 growth stages: the seedling stage (September 18–October 17, 2019), the flowering and fruit setting stage (October 18–November 16, 2019), the full fruiting stage (November 17–December 16, 2019) and the harvesting stage (December 17 2019–January 15, 2020). In this study, the maximum water demand during eggplant cultivation was calculated by the water balance method (Yuan et al., 2004), which was 200.48 m3 hm−2 15 days−1 at the seedling stage, 334.14 m3 hm−2 15 days−1 at the flowering and fruiting stage, 451.08 m3 hm−2 15 days−1 at the full fruiting stage, 350.84 m3 hm−2 15 days−1 at the harvesting stage. The irrigation amount 100% as the control (CK), and the specific irrigation amounts for the T1 (10% reduction), T2 (20% reduction) and T3 (30% reduction) treatments are shown in Table 1. Drip fertigation was used for irrigation. The 15-day irrigation amount was divided into three times, one time every 5 days. Irrigation water and fertilizers were prepared together every 5 days and stored in buckets, irrigated at once.
Determination of plant height and leaf area index
After planting, 5 eggplant seedlings were marked with the same growth in each treatment. Plant height was measured from 0.5 cm above the graft union to the growth point. The measurement was performed every 15 days.
The effects of different treatments on leaf area were shown by the leaf area index (LAI; Watson, 1947; Song and Huang, 1987). The formula for LAI is as follows:
Lij: Length of the selected leaves (cm).
Bij: Width of the selected leaves (cm).
m: The number of plants.
n: Total leaf number of the i-th (the tested) eggplant.
ρ: The population density of eggplant (plant m−2); the K value is 0.7 (our previous pre-experiments examined K values in leaf area index calculations for different varieties of eggplant by leaf area meter (CI-202, CID Bio-science, US) and leaf length and width calculation method, respectively. The results showed that the coefficient values (K) between the leaf area data measured by the leaf area meter and the calculated by length and width were from 0.66 to 0.71. For convenience of calculation, we used K = 0.7 in this study).
Determination of photosynthetic characteristics
Photosynthetic characteristic data were detected in the full fruiting stage.
The photosynthetic pigments were measured using the following method: 0.2 g of the leaves (without veins) was weighed into a 30 ml test tube, 20 ml of a 1:1 (V:V) mixture of absolute ethanol and acetone was added, and the mixture was sealed and stored in the dark for 48 h until the leaves were white and visible under ultraviolet light. The optical density (OD) values were measured at 663 nm, 646 nm and 470 nm in a spectrophotometer, and the chlorophyll and carotenoid contents were calculated according to the method of Zivcak et al. (2014). This process was repeated 3 times, and the average value was taken.
The photosynthetic parameters were measured using the following method: on a sunny day from 9:00 to 11:30, the third flat functional leaf was selected from the top, and the net photosynthetic rate (Pn), stomatal conductivity (Gs), intercellular CO2 concentration (Ci), and transpiration rate (Tr) were measured with the LI-6400/XT portable photosynthesis instrument produced by Li-cor (United States). A built-in light source with a light intensity of 500 μmol m−2 s−1 was used, the temperature was 27 (±1) °C, and the air CO2 concentration was 420 μmol mol−1. For each treatment, 5 plants were randomly selected, and the average value was calculated.
Determination of root system development index
Root length, root surface area, and root volume were measured by a MICROTEK root analysis system (Shanghai Zhongjing Co. Ltd., China). The data were analysed by LD-WinRHIZO software. Root vigour was measured using the following method: 0.5 g of the root was mixed with 10 ml of a 1:1 mixed 0.4% 2,3,5-triphenyl-2H-tetrazolium chloride (TTC) solution and 0.067 mol L−1 pH 7 phosphate buffer and kept in the dark at 37°C for 4 h, and 2 ml 1 mol L−1 sulfuric acid was added to stop the reaction. The sample was ground with 4 ml ethyl acetate and a small amount of quartz sand. The extract was transferred to a test tube, washed with ethyl acetate and diluted to 10 ml. A blank was used as a control, and the OD485 value was determined (Zhao et al., 2002).
Soil physicochemical characteristics and enzyme activities
Soil samples were collected by removing the plant roots during the harvesting stage. The soil sample was air-dried naturally, ground and passed through a 60-mesh sieve to determine the nutrient content. The available nitrogen content was determined by the alkaline solution diffusion method, the available phosphorus content was determined by the sodium bicarbonate method, and the available potassium content was determined by the ammonium acetate extraction-flame photometer method (Bhandari et al., 2002).
The colorimetric method was used to detect sucrase, urease, and alkaline phosphatase activity (Tabatabai and Bremner, 1969; Wyszkowska and Wyszkowski, 2010). The titration method was used to detect catalase activity (Zaborowska et al., 2019).
For catalase activity, 2 g of air-dried soil was placed in a 100 ml conical flask with 40 ml of distilled water and 5 ml of 0.3% hydrogen peroxide solution and shaken for 20 min, and then 1.5 mol-L−1 H2SO4 solution was added to stabilize the undecomposed hydrogen peroxide. The solution was filtered, and 25 ml of the filtrate was titrated with 0.1 mol ml−1 KMnO4 to the light red endpoint. The hydrogen peroxidase (catalase) activity was expressed as the number of millilitres of 0.1 mol ml−1 KMnO4 in 1 g of soil after 20 min.
For sucrose activity, 5 g of air-dried soil was placed in a 50 ml conical flask with 15 ml of 8% sucrose solution, 5 ml of phosphate buffer (pH 5.5), and 5 ml of toluene, shaken well and put in a constant-temperature chamber at 37°C for 24 h. After filtering, 1 ml of filtrate was mixed with 3 ml of 3,5-dinitrosalicylic acid, placed in a boiling water bath for 5 min, and cooled under running water for 3 min; then, the volume was adjusted to 50 ml, and the absorbance at 508 nm was measured. The 24 h 1 mg of glucose per g soil reflects the sucrase activity.
For urease activity, 5 g of air-dried soil was placed in a 50 ml conical flask with 1 ml of toluene, 10 ml of 10% urea solution and 20 ml of citrate buffer (pH 6.7), shaken well and placed in a constant-temperature chamber at 37°C for 24 h. After filtering, 3 ml of filtrate was diluted to 20 ml with distilled water and then mixed with 4 ml of sodium phenol and 3 ml of sodium hypochlorite solution, shaken well, and incubated for 20 min. The volume was adjusted to 50 ml, and the absorbance at 578 nm was measured. The milligrams of NH3-N in 1 g of soil after 24 h reflected the urease activity.
For alkaline phosphatase activity, 5 g of air-dried soil was placed in a 50 ml conical flask with 5 ml of toluene and 20 ml of 0.5% sodium benzene phosphate solution, shaken well, and then placed in a constant-temperature chamber at 37°C for 2 h. After filtering, 5 ml of filtrate was mixed with 20 ml of distilled water, 0.25 ml of ammonium chloride-ammonium hydroxide buffer (pH 9.8), 0.5 ml of 4-amino antipyrine solution, and 0.5 ml potassium ferricyanide solution, shaken well, and then adjusted to 50 ml. The absorbance was measured at 510 nm within 15 min, and the alkaline phosphatase activity was calculated as milligrams of P2O5 in 100 g of soil after 2 h.
Fruit quality determination
Eggplant fruits that had reached the same level of maturity under different treatments were collected, and the fruit quality indices were measured, including soluble protein, soluble sugar, vitamin C and free amino acid (Gutiérrez-Miceli et al., 2007; Zaller, 2007).
The soluble proteins were detected by the Coomassie Brilliant Blue method. The eggplant flesh was collected, and 1 g was added to 3 ml of phosphate buffer (pH = 7.8, 0.05 mol L−1), ground in an ice bath, and centrifuged at 10000 rpm and 4°C for 20 min. Then, 20 μl supernatant was taken and allowed to stand for 2 min, and its absorbance was measured at 595 nm.
The soluble sugars were detected by the enthrone colorimetric method. Eggplant flesh (0.3 g) was placed in a test tube with 10 ml of distilled water, sealed, and placed in a boiling water bath for 30 min to prepare an extract; the extraction of the flesh was repeated twice. The combined extracts were placed in a 25 ml volumetric flask and filtered, and 0.5 ml of filtrate was mixed with 1.5 ml of distilled water (for the blank without filtrate, 2 ml of distilled water was used), 0.5 ml of ethyl anthraquinone acetate reagent, and 5 ml of concentrated sulfuric acid, shaken thoroughly, sealed, placed in a boiling water bath for 1 min, and allowed to cool naturally; then, the OD value at 630 nm was measured.
Vitamin C (Vc) content was determined by the 2,6-dichlorophenol-indophenol titration method. Eggplant flesh (5.0 g) was placed in a mortar with 5 ml of 2% oxalic acid, ground into a homogenate, and transferred to a 50 ml volumetric flask; the volume was adjusted with 2% oxalic acid, and then it was filtered. The filtrate was pipetted into a 50 ml conical flask and titrated with calibrated 0.1% 2,6-dichlorophenol indophenol solution until it was peachy red and did not fade for 15 s. The amount of dye added was recorded.
The free amino acid amount was detected by the ninhydrin colorimetric method. Eggplant flesh (0.5 g) was placed in a mortar, and 5 ml of 10% acetic acid was added and ground into a homogenate. The residue was rinsed with ammonia-free water, adjusted to 50 ml, filtered and prepared for use. Then, 1 ml of filtrate was placed in a 20 ml stoppered test tube, mixed well with 1 ml of ammonia-free distilled water, 3 ml of ninhydrin reagent, and 100 μl of 0.1% ascorbic acid, sealed, placed in a boiling water bath for 15 min, and cooled under running water until the red colour faded to blue–purple. Then, it was adjusted to 20 ml with 60% ethanol, and the absorbance was detected at 570 nm.
Yield statistics
After planting, each plot was marked with 20 plants for yield statistics. The eggplant plant was subjected to the double-stem pruning method. The individual mature fruit were weighed, and the yield was calculated per plant and converted into yield per hectare.
Soil microorganism DNA extraction and sequencing
Bacterial and fungal DNA was extracted using a PowerSoil DNA Isolation Kit (MoBio, United States). Common primers were used to amplify 338F (5′-ACTCCTAGGGAGCA-3′) and 806R (5′-GGACTCHVGGGTWTTAT-3′) of the V3-V4 region of the bacterial 16S rRNA gene binding adapters and barcode sequences (Quast et al., 2013). Fungal ITS1 primers were used for ITS1 (5′-CTGTCATTAGGGAGAGA-3′), and ITS2 (5′-GCTGCGTTCTTCATCGATGA-3′) was amplified by binding adapter sequences and barcode sequences (Kõljalg et al., 2013). The final products were purified, quantified and homogenized to generate a sequencing library. The constructed libraries for sequencing fragment DNA were quality checked and mass concentration was measured with Qubit fluorometer. Library fragment size was measured with Agilent 2,100., and the qualified libraries were sequenced using Illumina HiSeq 2,500 (Illumina, USA). All raw sequences have been deposited in the NCBI sequence read archive under BioProject accession number PRJNA771427. Each treatment had 6 repeats.
Bioinformatics analysis
Overlapping reads of the samples were spliced using FLASH software (version 1.2.111), and the spliced sequences obtained were the raw label data. The original tags obtained from the splicing were filtered using Trimmomatic software (version 0.33) to obtain high-quality tag data. Uchime software (version 8.1) was used to identify and remove the chimeric sequences to obtain the final data. The tags were clustered using Uparse software (version 7.0.1090) to obtain operational taxonomic units (OTUs) with 97% similarity. The tags were based on bacterial 16S sequences from the SILVA database (Release 138) and fungal ITS sequences from the UNITE database (Release 8.0). The RDP classifier (version 2.24) was used to classify the taxa with a minimum reliability estimate of 80%. The alpha diversity of the microbial communities was analysed by Mothur (version v.1.305) using Chao1 and ACE richness estimators, and the microbial diversity was measured using Shannon–Wiener and Simpson diversity indices. The FAPROTAX database6 linking species taxonomy and functional annotations was constructed and then linked to the OTU taxonomic table and the FAPROTAX database to predict the function of the bacterial communities (Louca et al., 2016). Fungal gene sequence information was obtained from Funguild7 associated with fungal function and used to predict fungal community function (Xie et al., 2020). The relative abundance of the species in different samples was determined using STAMP8 with the G-test and Fisher’s test. A two-sample t test was used to determine significant differences with a p value threshold of 0.05 (Chen et al., 2021).
Data analysis and statistics
The bioinformatics data were analysed on the Majorbio Cloud Platform1 (Shanghai Majorbio Biopharm Technology Co., Ltd.). The data were analysed and plotted (Venn diagram, community bar chart) using R Programming Language (version 3.3.1). Circos was plotted using Circos-0.67-7 software, and significant differences were tested using analysis of variance and the least significant difference test of the agricolae package (version 1.3.3). Linear discriminant analysis (LDA) was performed on six samples based on the taxonomic composition conditioned from the phylum to genus level using LEfSe software to identify communities or species that produced significant differential effects with a log LDA score threshold of 2.0 (Zhou et al., 2019) and to identify the most likely traits to explain the differences between categories (Segata et al., 2011). Principal coordinate analysis (PCoA) based on unweighted UniFrac distance was constructed using R Programming Language (version 3.3.1) to measure the microbial community beta diversity and compare the pairwise differences based on the Wilcoxon rank-sum test. One-way ANOVA was used to test for significant differences between groups, and the stats package of R Programming Language (version 3.3.1) and the SciPy package of Python were used for analysis and plotting. Phylogenetic evolutionary trees were constructed using FastTree software with neighbor-joining at the genus level by Jukes-Cantor distance matrix based on the sequencing, and evolutionary trees were plotted using the R programming language (version 3.3.1). Bacterial metabolic function was predicted, analysed and plotted using the PICRUSt software package. Fungal metabolic function was predicted, analysed and plotted using the FUNGuild package. Prokaryotic (bacterial) ecological functions were predicted using the FAPROTAX database and analysed and plotted using Python. RDA/CCA was based on the experimentally obtained average soil enzyme activity data using R programming language’s vegan package. Correlation heatmaps were analysed and plotted using the R programming language (version 3.3.1; heatmap package). A structural equation model (SEM) of the direct and indirect effects of treatments, physicochemical factors and microorganisms on yield was developed using AMOS software (IBM SPSS Amos 23) by maximum likelihood estimation (MLE). Data from the different treatments were used to construct the SEM. The treatments represent the proportion of organic fertilizer for each treatment: 1, 0.9, 0.8 and 0.7 for the CK, T1, T2 and T3 treatments, respectively. Fungi (Chao1), soil enzymes (catalase, sucrase, urease, and alkaline phosphatase activity) and remaining soil nutrients after harvest (available N, P, and K) represent the RDA1 data for the first axis of the corresponding indicator RDA/CCA. Yields are the sum of the autumn-winter greenhouse eggplant yields.
Statistical analysis was conducted using SPSS 19.0. The figures and table were drawn using GraphPad Prism 7 and Microsoft Excel. All data were analysed by Tukey’s HSD test at significance levels of p < 0.05, 0.01 or 0.001.
Results
The effect of different irrigation amounts on eggplant plant height, leaf area index, root morphological indices and vitality
Eggplant height increased with prolongation of the growth period (Figure 1A). Among all treatments, plants were taller under CK treatment after 60 days, while there was no significant difference between T1 and CK (p > 0.05). The leaf area index of eggplant first increased and then decreased (Figure 1B), reaching its maximum at 75 days after planting. Compared with CK, the leaf area index under all treatments decreased. At 75 days after planting, the leaf area index of T1, T2 and T3 decreased by 5.9, 11.8 and 19.7% compared with that of CK, respectively. There was no significant difference (p > 0.05) between T1 and CK. At 105 days, the leaf area index of T1, T2 and T3 decreased by 3.3, 16.6 and 19.3% compared with that of CK, respectively, and there was no significant difference between T1 and CK (p > 0.05).
Figure 1. The effects of different irrigation amounts on plant height (A), leaf area index (B), root length (C), root area (D), root volume (E), and root vigour (F). n ≥ 3; different letters indicate significant differences (p < 0.05).
The root length of eggplant under T1 and T2 increased significantly, by 16.3 and 31.1%, respectively, compared with that of CK. Under T3, it was significantly reduced by 12.9% compared that of CK (Figure 1C). The root surface area increased significantly by 37.5 and 22.6% under T1 and T2, respectively, compared with that of CK, and the root surface area under T3 was significantly lower than that of CK (Figure 1D). The root volume was the largest under T2, with a significant increase of 13.6% compared with that of CK (p < 0.05). There was no significant difference between T1 and CK (p > 0.05), but the root volume was significantly reduced in T3 significantly, by 15.3% (Figure 1E). The root vigour was highest under T1, which was 18.7% higher than that under CK, and the difference was significant (p < 0.05). The root vigour under T2 was 5.3% lower than that under CK, but the difference was not significant (p > 0.05). The root vigour under T3 was significantly lower than that under CK (Figure 1F).
The effect of different irrigation amounts on photosynthetic parameters in the leaves of eggplants
The chlorophyll a content was higher under T1 and T3, increased significantly by 15.3 and 14.4% compared to CK (Table 2). There was no significant difference in chlorophyll b content between all treatments and CK (p > 0.05). The carotenoid content was highest under T3, 23.3% higher than CK, and the difference was significant, but there was no significant difference between T1, T2 and CK (p > 0.05). The chlorophyll (a + b) content of T1 and T3 was significantly higher (9.1 and 9.9%, respectively) than that of CK.
Table 2. Effects of different irrigation amounts on photosynthetic pigments content in leaves of eggplants.
The net photosynthetic rate (Pn) in leaves of eggplant was significantly higher under T1 and T2, with 33.0 and 26.4% increases in comparison with that under CK, but there was no significant difference between T1 and T2 (p > 0.05; Table 3). There was no significant difference in the transpiration rate (Tr) between T1 and T2 and CK (p > 0.05), but the Tr was significantly lower, by 13.8%, under T3 than under CK. Compared with CK, stomatal conductance (Gs) was significantly highest under T1, with a value 12.1% higher than that under CK. The intercellular carbon dioxide concentration (Gi) under T1 was not significantly different from that under CK, while the values under T2 and T3 were significantly lower than that under CK. All treatments significantly improved the water use efficiency (Pn/Tr), which was 23.0% higher under T1 than under CK, but there was no significant difference among T1, T2 and T3 (p > 0.05).
Table 3. Effects of different irrigation amounts on photosynthetic parameters in leaves of eggplants.
The effect of different irrigation amounts on soil enzyme activity and the contents of available nitrogen, phosphorus, and potassium
As shown in Figure 2, during the flowering and fruit setting stage, soil catalase activity was 5.8% higher under T1 than under CK, and the difference was significant, but there was no significant difference between T1 and T2 (p > 0.05). Soil alkaline phosphatase activity and urease activity were highest under T1, with values that were 59.2 and 21.7% higher than those under CK, respectively, and the differences were significant. Soil sucrase activity was significantly higher under T3, and there was no significant difference between the T1 and T2 treatments and CK (p > 0.05). At the full fruiting stage, compared with CK, soil catalase and alkaline phosphatase activities were 25.1 and 62.0% higher under T1, respectively, and the difference was significant. The urease activity was significantly higher in T1, with a value 22.9% higher than that in CK, and the urease activity in T2 and T3 was also higher than that in CK. The sucrase activity increased with decreasing irrigation amounts and was highest under T3, in which it was 21.5% higher than that under CK, a significant difference. There was no significant difference between T1 and CK (p > 0.05).
Figure 2. The effects of different irrigation amount on catalase activity (CAT; A), sucrase activity (SA; B), urease activity (UA; C), alkaline phosphatase (APA; D), available nitrogen content (E), available phosphorus content (F) and available potassium content (G). n ≥ 3; different letters indicate significant differences (p < 0.05).
The results showed that the contents of available nitrogen (Figure 2E), phosphorus (Figure 2F) and potassium (Figure 2G) were higher under T3 than under CK by 25.9, 51.4 and 62.5%, respectively. T2 followed, but there was no significant difference in available nitrogen content between T2 and T3 or between T1 and CK (p > 0.05). However, the available phosphorus content in T1 was significantly lower than that in CK (p < 0.05). The available potassium content was also reduced by 12.3% in T1 compared to that in CK, and the difference was significant (p < 0.05). Above all, these results implied that the rhizosphere soil nutrient content was reduced with increasing irrigation amounts, which may be due to the high nutrient utilization rate.
Effects of different irrigation amounts on the OTU number and diversity of soil microbial communities
All samples were sequenced, the same number of obtained OTUs were extracted randomly and annotated by species taxonomy, and the corresponding abundance information of each OTU annotation result in each sample was counted and displayed as a rank–abundance curve (Supplementary Figures S1A,D). The pan/core species curve was relatively flat, indicating that the sequencing sample size was sufficient (Supplementary Figures S1B,C,E,F).
Table 4 shows that the bacterial OTU numbers were higher under T3 and lower under T2; the fungal OTU numbers were lower under T3 and higher under T1, but there were no significant differences among the treatments (p > 0.05). A total of 3819 bacterial OTUs were present under all treatments, with approximately 74.3–78.3% of OTUs per treatment, and approximately 1.0–1.5% of the OTUs per treatment were different from the other treatments (Figure 3A). A total of 132 fungal OTUs were present under all treatments, accounting for 52.4–62.6% of OTUs per treatment, and approximately 10–17% of the OTUs per treatment were different from those found in the other treatments (Figure 3B). Thus, compared to bacteria, the fungal OTUs differed more between each treatment.
Table 4. Effects of different irrigation amounts on alpha diversity of bacteria and fungi in OTU level.
Figure 3. The effect of different irrigation amounts on the similarity and overlap of bacterial (A) and fungal (B) OTU compositions in soil; the similarity level was higher than 97% (Venn). Circos reflects the distribution ratio of dominant phyla in different treatments and the distribution ratio of each dominant phylum in different treatments, including bacteria (C) and fungi (D). n = 6.
The results obtained from sequencing were used to analyse all OTUs based on the OTU level, and the same number of obtained OTUs was randomly extracted to analyse the soil microbial community. The flatness of the rarefaction curve indicates that the amount of sequencing data is reasonable (Supplementary Figure S2). The results showed that the alpha diversity of the soil microbial community was affected differently by different irrigation amounts (Table 4), with the fungal ACE index being higher under T1 and significantly lower under CK and T3. The bacterial Chao1 index was higher under T3 and significantly lower under T1, while the fungal Chao1 index was higher under T1 and significantly lower under CK and T3. These results indicated that the fungal community richness was higher under T1 and that the bacterial community richness was higher under T3. The fungal Simpson index was highest under T3 and significantly lower under the other treatments. The fungal Shannon index was higher under T1, but there was no significant difference (p > 0.05) between CK, T1 and T2, and it was significantly lower under T3. This finding indicates that the fungal community diversity was highest in T1. The rest of the indices were not significantly different.
The beta-diversity of soil microbial communities was also affected differently by different irrigation levels. PCoA based on the weighted UniFrac distance algorithm showed that the fungal communities were more dispersed than the bacterial communities under different irrigation treatments, where the fungal communities of T3 were significantly dispersed from those of the other treatments, the fungal communities under T1 and T2 were significantly dispersed from that of CK, and the bacterial communities under T3 were significantly separated from those of CK (Figures 4A,B; Supplementary Figures S3A,B). The results of ANOSIM/Adonis analysis showed that the between-group differences were significantly greater than the within-group differences for all treatments (Supplementary Figures S3C,D).
Figure 4. Principal coordinate analysis (PCoA) of bacterial (A) and fungal (B) OTU levels in soil at different irrigation amounts with an unweighted UniFrac distance algorithm. LEfSe analysis from the phylum to genus level, including bacteria (C) and fungi (D), by the nonparametric factorial Kruskal–Wallis (KW) sum-rank test. n = 6.
Effects of different irrigation amounts on the structure of soil microbial communities
The community distributions of bacteria and fungi are shown in the Circos and community bar diagrams with relative abundance, including differences at the bacterial phylum level (Figure 3C; Supplementary Figure S4A) and at the fungal phylum level (Figure 3D; Supplementary Figure S4B) and genus level (Supplementary Figures S4C,D). Ten bacterial phyla were detected: Proteobacteria, Chloroflexi, Actinobacteriota, Firmicutes, Acidobacteriota, Gemmatimonadota, Bacteroidota, Myxococcota, Patescibacteria and others. Proteobacteria was found to be the most abundant phylum under each treatment except T3, with 21.86, 22.33, 21.61 and 20.19% in CK, T1, T2 and T3, respectively. A higher relative abundance of Chloroflexi was found under T3 (20.54%). Chloroflexi and Acidobacteriota had lower relative abundances under T1, and their relative abundances increased gradually with decreasing irrigation amount. Five fungal phyla were detected, and the relative abundance of Ascomycota was higher under each treatment, with values of 86.91, 79.06, 81.67 and 48.62% under CK, T1, T2 and T3, respectively, with the value under T3 representing a significantly lower value. A total of 17 fungal genera were detected, with a higher relative abundance of Cladosporium under all treatments except T3, and a higher relative abundance of Wallemia was detected in T3.
While the significant differences between species among treatments are shown in Supplementary Figure S5, the differences were significant for Deinococcota at the bacterial phylum level (p < 0.01) and for the rest of the phyla, as shown in the Supplementary Figure S5A. The difference between Ascomycota and Basidiomycota at the level of the fungal phylum was significant (p < 0.001), while the differences in rest of the phyla in Supplementary Figure S5B were significantly different (p < 0.05). Significant differences were found between Norank A4b and Norank SBR1031 (Chloroflexi; p < 0.001), Sphingomonas, Norank Saccharimonadales and Truepera (p < 0.01) at the level of the bacterial genus, and the rest were significantly different (p < 0.05; Supplementary Figure S5C). Significant differences were found for Wallemia (p < 0.001), Fusarium, Phialosimplex and Lecythophora (p < 0.01) at the level of fungal genera, and the rest of the genera in the figure (p < 0.05; Supplementary Figure S5D).
The LEfSe results showed that there were 60 bacterial communities and 35 fungal communities with significant differences (Figures 4C,D). The bacteria family Sphingomonada ceae and the fungus phylum Ascomycota were the most abundant under CK, the bacteria class Bacilli and the fungus family Aspergillaceae were the most abundant under T1, the bacteria genus norank f A4b and fungus class Sordariomycetes the most abundant under T2, and the bacteria phylum Chloroflexi and fungus class Wallemiomycetes were the most abundant under T3 (SM 2 and 3). Phylogenetic tree-and genus-level reads showed the same results (Supplementary Figure S6).
Effects of different irrigation amounts on bacterial functional metabolism and fungal trophic patterns
The bacterial community showed 25 COG functional classifications under all treatments, as predicted by PICRUSt1, with the highest abundance of amino acid transport and metabolism functional classes. However, the differences in each functional class under CK and the three treatments were not significant (Figure 5A). By FAPROTAX annotation, the bacterial community showed 50 functions, among which the relative abundance of chemoheterotrophy and aerobic chemoheterotrophy was higher. The functions with significant differences (p < 0.05) in abundance between treatments are shown in Figure 5. Interestingly, these differences were mainly associated with N cycling and nutrient catabolism. With decreasing irrigation amount, nitrate reduction, nitrogen fixation, manganese oxidation, fumarate respiration, anoxygenic photoautotrophy S oxidation, xylanolysis, and chitinolysis were reduced or significantly reduced, but there was no significant difference between T1 and CK (p > 0.05). With the decrease in irrigation water, the functions of phototrophy, oxygenic photoautotrophy, photoheterotrophy and nitrification increased or significantly increased (Figure 5B).
Figure 5. Clusters of Orthologous Groups (COG) functional classification of PICRUSt in bacteria (A) of different treatments. The functional heatmap predicated by FAPROTAX in bacteria (B) under different treatments. n = 6; different letters indicate significant differences (p < 0.05).
The trophic pattern of the fungal community was dominated by saprotrophs (Figure 6A). The proportion of animal endosymbiont–animal pathogen–endophyte–plant pathogen–undefined saprotrophs decreased, and the proportion of undefined saprotrophs increased with decreasing irrigation amount. The proportion of endophyte–litter saprotroph–soil saprotroph-undefined saprotrophs and animal pathogen–endophyte–lichen parasite–plant pathogen–soil saprotroph–wood saprotrophs was higher in T1 and T2 than in CK and T3, while the proportion of animal pathogens was significantly higher under T1, and the proportion of plant pathogens was significantly higher under T2 than under the other treatments. The FAPROTAX functional prediction heatmap showed that the fungal community functions included cellulolysis, xylanolysis and chemoheterotrophy, and the functional abundance gradually decreased with decreasing irrigation amount, but the difference between treatments was not significant (Figure 6B).
Figure 6. Taxonomic analysis of fungal communities by microecological guild (FUNGuild; A) of different treatments. The functional heatmap predicated by FAPROTAX in fungi (B) under different treatments. n = 6.
Relationship between microbial community structure and soil enzyme activity in different treatments
The soil microbial community structure was closely related to soil enzyme activity. Overall, the correlation between soil microbial community structure and soil enzyme activity on both dimensions (RDA1 and RDA2) was 26.89% for bacterial communities and 58.53% for fungal communities (Figures 7A,B). Table 5 shows that the soil microbial community structure was significantly correlated (p < 0.01) with catalase activity (CAT; bacteria R = 0.6251, fungi R = 0.7587) and sucrase activities (SA; bacteria R = 0.7518, fungi R = 0.9080). There was no significant correlation between soil microbial community structure and urease activity (UA) or alkaline phosphatase activity (APA; p > 0.05; Table 5).
Figure 7. The ranking method developed by redundancy analysis and correspondence analysis (RDA/CCA) in different treatments. The environmental factor adopts the mean value of soil enzyme activity, including bacteria (A) and fungi (B), at the genus level. The relationships and correlations between different phyla in bacteria (C) and fungi (D) with soil enzyme activities. n = 6, different stars indicate significant differences (*p < 0.05; **p < 0.01; ***p < 0.001).
The correlation between the relative abundance of the microbial community at the phylum level and soil enzyme activity under different irrigation treatments was analysed by the Spearman rank correlation coefficient and was shown by a heatmap (Figures 7C,D). Patescibacteria was significantly and positively correlated with SA activity (p < 0.001); Cyanobacteria was significantly and positively correlated with SA activity (p < 0.01); SAR clade Marine Group B, Latescibacteria, Armatimonadota, unclassified k norank d Bacteria, Acidobacteria, and Chloroflexi were significantly and positively correlated with SA activity (p < 0.05); and TX1A-33 was significantly and positively correlated with CAT activity. Margulis bacteria showed a significant negative correlation with CAT and APA activity (p < 0.01), and Chloroflexi showed a significant negative correlation with CAT activity (p < 0.01). TX1A-33, Deinococcota and Firmicutes showed a significant negative correlation with SA activity (p < 0.05), and Bdellovibrionota showed a significant negative correlation with UA activity (p < 0.05). Among the fungal communities, Basidiomycota showed a significant positive correlation with SA and UA activity (p < 0.001). Rozellomycota, Chytridiomycota, and Mortierellomycota showed a significant positive correlation with CAT activity (p < 0.05). Blastocla diomycota showed a significant negative correlation with APA activity (p < 0.01) and with CAT activity or UA activity (p < 0.05). Ascomycota showed a significant negative correlation with SA activity (p < 0.001) and UA activity (p < 0.05).
The effect of different irrigation amounts on the fruit quality and yield of eggplant
As shown in Table 6, the contents of soluble sugars, soluble proteins and free amino acids were all significantly higher under T1 and T2, which were 10.8, 12.3, and 6.7% and 6.6, 6.8, and 3.1% higher than those under CK, respectively. The soluble sugar, soluble protein and free amino acid contents of the fruit under T3 were 4.2, 4.3 and 3.2% higher than those under CK, respectively, values that were not significantly different (p > 0.05). There was no significant difference in the vitamin C (Vc) content among all treatments (p > 0.05). After yield, compared with CK, T1 increased yield by 3.9%, but the difference was not significant (p > 0.05). The yield in T2 and T3 decreased by 14.4 and 16.7%, respectively, compared with that in CK, and the difference between T3 and CK was significant (p < 0.05).
Structural equation modelling
According to the above results, the effect of different irrigation amounts on the fungal community in the soil was greater than that of bacteria; therefore, the structural equation model (SEM) analysis of the direct and indirect effects of different irrigation levels was based on soil fungi and included soil enzyme activity, soil nutrient residuals and autumn–winter crop yield of eggplants. The model fit the data well (χ2 = 1.544, p = 0.213, Goodness of fit index/GFI = 0.930, Normed fit index/NFI = 0.979, Incremental of fit index/IFI = 0.993, Comparative fit index/CFI = 0.992, Root mean square error of approximation/RMSEA = 0.190) and explained 25% of the variation in soil fungi, 93% of the variation in soil nutrients, 60% of the variation in soil enzyme activity and 70% of the variation in yield of eggplants (Figure 8). The results showed that the different irrigation treatments significantly affected the distribution of the fungal community structure (λ = 0.50, p < 0.05) and soil enzymes (λ = 0.50, p < 0.01), but the effect on nutrients was not significant. The fungal community structure had a significant effect on the soil enzymes (λ = 0.39, p < 0.05). More soil nutrients at the harvesting stage indicated lower nutrient utilization, the effect of soil enzymes on nutrients was extremely significant (λ = −1.06, p < 0.001), and the yield was significantly affected by the fungal community structure (λ = 0.46, p < 0.01) and residual soil nutrients (λ = −0.50, p < 0.01).
Figure 8. The values above the structural equation model (SEM) lines indicate path coefficients. The red line represents positive path coefficients, the blue line represents negative path coefficients, and the black line represents nonsignificant path coefficients. The width of the arrow indicates the significance of the standard path coefficient (***p < 0.001; **p < 0.01; *p < 0.05).
Discussion
As climate change intensifies and the population continues to grow rapidly, water scarcity is becoming increasingly pronounced. This study investigated the effects of irrigation water amounts on eggplant growth and development, yield, rhizosphere microbial community structure and function, and rhizosphere soil physicochemical properties. The results showed that, compared with full irrigation (CK), 10% reduced irrigation (T1) treatment improved plant height, fruit yield and quality as well as photosynthetic rate, transpiration rate, stomatal conductance and water use efficiency. The structure and function of the soil microbial community were also significantly altered. Under the 10% reduced irrigation (T1) treatment, the richness and diversity of the microbial community were significantly increased, the activities of soil catalase, urease and alkaline phosphatase increased, and the N cycling rate and nutrient catabolism also increased.
With the development of the agricultural industry, the demand for technological improvements in the field of agricultural research is growing. Much research focuses on the crop itself in efforts to enhance the potential of the crop (Hwang et al., 2017; Baileyserres et al., 2019), but this focus is far from sufficient. Modern agriculture needs more attention, including crop growth environments, crop growth substrates, agricultural facilities, and root microbes (Rouphael et al., 2018; Saiz-Rubio and Rovira Más, 2020). Water saving is one of the important goals of agricultural production (Chai et al., 2015). Therefore, to achieve the goal of saving water, it is not sufficient to develop agriculture only from the perspective of crops (Zhang et al., 2018). The innovation of agricultural facilities and technologies is also very important (Pretty and Bharucha, 2014).
In this research, we used the calculated theoretical greenhouse eggplant water requirements as a control and reduced the irrigation volume by 10, 20, and 30% relative to the control (100%). The results showed that the treatment in which the irrigation amount was reduced by 10% did not significantly decrease the plant height compared with that in CK, and the difference in leaf area was also not significant. However, after the irrigation amount was reduced by 20% or 30%, the plant height and leaf area decreased, which had an obvious adverse effect on eggplant yield. Therefore, an irrigation amount under the 10% reduced irrigation treatment can ensure the water demands of eggplant during different growth stages. The 10% reduced irrigation treatment promoted the development of the root system and increased the root area and root vigour, thereby increasing the ability of eggplant to absorb nutrients. Meanwhile, the soil enzyme activity and soil available nutrient utilization, including available nitrogen, phosphorus, and potassium, were all improved.
Studies have shown that an ebb-and-flow subirrigation system can be used to increase water use efficiency in greenhouse cucumber planting; the enrichment of specific fungi increases under the nutrient solution, but it does not increase the cucumber disease rate (Dong et al., 2020). Specific microorganisms can promote crop resistance to water loss (Rajabi-Khamseh et al., 2020) but may also affect plant root development (Bennett and Klironomos, 2019). The activities of catalase, urease and alkaline phosphatase were all higher under the 10% reduced irrigation treatment in the flowering and fruit setting and full fruiting stages. However, the available nitrogen, phosphorus and potassium contents in the soil were lower under the 10% reduced irrigation treatment than under the control treatment. These results showed that although the microbial community structure under the 10% reduced irrigation treatment was changed, compared to CK, the vitality and fertility of the soil were significantly improved, and the soil nutrients could be utilized by plants to the greatest extent. An improved root system increases the absorption of soil nutrients by crops and promotes growth (Zhang et al., 2021) through changes in the osmosis capacity of root cells or changes in cell wall characteristics (Al-Yasi et al., 2020). This root-soil interaction can significantly increase the yield and quality of crops (Wang et al., 2020). This impact is multifaceted, including root growth (Luo et al., 2016), phytohormone (Puértolas et al., 2020), soil nutrient cycling (Du et al., 2008), root redox balance regulation (Bogale et al., 2016) and other factors.
The impact of irrigation on crops is multidimensional and can directly affect crop water supply and indirectly affect crop yield and quality by improving soil quality and the soil microbial composition (Rousk et al., 2010; Chodak et al., 2015; Santos et al., 2021). In this study, through the sequencing and analysis of soil microorganisms, it was found that the microbial community structure and function in the soil were changed by decreasing irrigation. Some interesting findings were also found, such as the high abundance of Firmicutes, Basidiomycota, and Mortierellomycota and the low abundance of Acidobacteriota under the 10% reduced irrigation treatment. The genus Bacillus of the phylum Firmicutes is considered to be a beneficial microorganism in agroecology, promoting crop growth and enhancing plant immunity (Hashmi et al., 2020). Bacilli (class) are also considered to promote the decomposition of organic fertilizers and enhance the nutrient utilization of rice (Liu et al., 2021). Basidiomycota exhibit the ability to form mycorrhizae in symbiosis with crops and accelerate the decomposition of lignin and cellulose in soil (Vohník et al., 2012; Genre et al., 2020). The abundance of Wallemia was significantly increased under T3, and the genus Wallemia of the phylum Wallemia is a unique drought-and salt-tolerant fungus (Zalar et al., 2005; Kralj Kunčič et al., 2010), which also suggested that the soil moisture deficit of T3 caused greater stress to the crops.
The richness and diversity of the fungal community under the 10% reduced irrigation treatment were higher than those under the other treatments. While the richness and diversity of the fungal community under T3 were also higher, the dispersion between the community structure in T3 and other treatments was larger, showing a notable difference. Under the 10% reduced irrigation treatment or CK treatment, the nutrient metabolism capacity of the microbial community was higher than that under the other treatments, especially the N cycle, and the decomposition capacity of xylan and chitin was improved, which suggests that community catabolism was more active. The 10% reduced irrigation treatment increased soil CAT and SA activities, which may be related to the increase in communities such as Patescibacteria and Basidiomycota. The results of the structural equation model also show that different irrigation treatments can indirectly and significantly influence the yield of eggplants by affecting the microbial community structure, soil enzyme activity and soil nutrients, and the different irrigation treatments explained 70% of the variation in yield.
The regulation of photosynthesis is also an important factor that affects crop growth and yield (Wu et al., 2019). Studies have shown that using an improved irrigation method, alternate drip irrigation, can significantly increase the photosynthetic rate of crops and improve crop quality and yield (Wang et al., 2020). Many studies have shown that partial root zone irrigation can improve the crop transpiration rate (Du et al., 2006), increase the photosynthesis rate (Romero et al., 2012), and lead to the accumulation of more photosynthetic products (Shahnazari et al., 2007). Water deficit affects chlorophyll fluorescence parameters and pigments (Hazrati et al., 2016), photosynthetic electron transport and specific photoprotective responses (Zivcak et al., 2013). However, some soil microorganisms can protect the photosynthetic organs of crops under water deficit (Mathur et al., 2019) and can even increase the stomatal conductance compared with full watering (Augé et al., 2015), thereby improving the ability of plants to resist water loss (Rolli et al., 2015). This ensures that crop fruit quality and yield will not decrease and can even increase under relatively low irrigation. In our research, compared with full irrigation (CK), reducing the irrigation amount by 10% (T1) significantly increased the content of chlorophyll a in eggplant and increased the net photosynthetic rate, transpiration rate, stomatal conductance and water use efficiency.
In conclusion, our research proves that for greenhouse eggplants cultivated in autumn and winter, the optimal irrigation amounts at the seedling stage, flowering and fruit setting stage, full fruiting stage and harvesting stage are 180.43, 300.60, 405.97 and 315.76 m3 hm−2 15 d−1, respectively. The 10% irrigation reduction (T1) treatment enhanced the richness and diversity of eggplant rhizosphere microbial communities and improved N cycling and nutrient catabolism, thereby increasing soil catalase, urease and alkaline phosphatase activities and improving soil nutrients. Plant-microbe interactions in the T1 treatment resulted in a significant increase in eggplant root length, surface area and vigor compared to those in CK. The contents of soluble sugar, soluble protein and free amino acids in fruit and the yield of eggplants were increased In T1, the bacteria Bacilli and the fungus Aspergillaceae were the most abundant, while the relative abundances of Chloroflexi and Acidobacteria were lower than those of the other treatments. However, in T3, the Basidiomycota (genus Wallemia) was the most abundant.
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 in the article/Supplementary material.
Author contributions
FY, MW, and NW designed the study. TJ, XG, and FW carried out the experiments and analysed the data. PJ performed SEM analysis. TJ wrote the manuscript. FY, JL, and NW revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key Research and Development Program (2019YFD1001904), the National Natural Sciences Foundations of China (NO. 32172556), and the Shandong Provincial Key Research and Development Program (2021LZGC017).
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.981288/full#supplementary-material
Footnotes
References
Al-Yasi, H., Attia, H., Alamer, K., Hassan, F., Ali, E., Elshazly, S., et al. (2020). Impact of drought on growth, photosynthesis, osmotic adjustment, and cell wall elasticity in damask rose. Plant Physiol. Biochem. 150, 133–139. doi: 10.1016/j.plaphy.2020.02.038
Augé, R. M., Toler, H. D., and Saxton, A. M. (2015). Arbuscular Mycorrhizal Symbiosis alters stomatal conductance of host plants more under drought than under amply watered conditions: a meta-analysis. Mycorrhiza 25, 13–24. doi: 10.1007/s00572-014-0585-4
Ault Toby, R. (2020). On the essentials of drought in a changing climate. Science 368, 256–260. doi: 10.1126/science.aaz5492
Baileyserres, J., Parker, J. E., Ainsworth, E. A., Oldroyd, G. E. D., and Schroeder, J. I. (2019). Genetic strategies for improving crop yields. Nature 575, 109–118. doi: 10.1038/s41586-019-1679-0
Bennett, J. A., and Klironomos, J. (2019). Mechanisms of plant-soil feedback: interactions among biotic and abiotic drivers. New Phytol. 222, 91–96. doi: 10.1111/nph.15603
Berdugo, M., Delgadobaquerizo, M., Soliveres, S., Hernandezclemente, R., Zhao, Y., Gaitan, J. J., et al. (2020). Global ecosystem thresholds driven by aridity. Science 367, 787–790. doi: 10.1126/science.aay5958
Bhandari, A. L., Ladha, J. K., Pathak, H., Padre, A. T., Dawe, D., and Gupta, R. K. (2002). Yield and soil nutrient changes in A long-term Rice-wheat rotation in India. Soil Sci. Soc. Am. J. 66, 162–170. doi: 10.2136/sssaj2002.1620a
Bodner, G., Nakhforoosh, A., and Kaul, H.-P. (2015). Management of crop water under drought: a review. Agron. Sustain. Dev. 35, 401–442. doi: 10.1007/s13593-015-0283-4
Bogale, A., Nagle, M., Latif, S., Aguila, M., and Müller, J. (2016). Regulated deficit irrigation and partial root-zone drying irrigation impact bioactive compounds and antioxidant activity in two select tomato cultivars. Sci. Hortic. 213, 115–124. doi: 10.1016/j.scienta.2016.10.029
Chai, Q., Gan, Y., Zhao, C., Xu, H.-L., Waskom, R. M., Niu, Y., et al. (2015). Regulated deficit irrigation for crop production under drought stress. A review. Agron. Sustain. Develop. 36:3. doi: 10.1007/s13593-015-0338-6
Chen, M., Zhu, X., Zhao, C., Yu, P., Abulaizi, M., and Jia, H. (2021). Rapid microbial community evolution in initial Carex litter decomposition stages in Bayinbuluk alpine wetland during the freeze-thaw period. Ecol. Indic. 121:107180. doi: 10.1016/j.ecolind.2020.107180
Chodak, M., Gołębiewski, M., Morawska-Płoskonka, J., Kuduk, K., and Niklińska, M. (2015). Soil chemical properties affect the reaction of forest soil bacteria to drought and rewetting stress. Ann. Microbiol. 65, 1627–1637. doi: 10.1007/s13213-014-1002-0
Dong, C. J., Li, Q., Wang, L. L., and Shang, Q. M. (2020). Dynamic changes in bacterial communities in the recirculating nutrient solution of cucumber plug seedlings cultivated in an ebb-and-flow subirrigation system. PLoS One 15:E0232446. doi: 10.1371/journal.pone.0232446
Du, T., Kang, S., Zhang, J., and Davies, W. J. (2015). Deficit irrigation and sustainable water-resource strategies in agriculture for China’s food security. J. Exp. Bot. 66, 2253–2269. doi: 10.1093/jxb/erv034
Du, T., Kang, S., Zhang, J., Li, F., and Hu, X. (2006). Yield and physiological responses of cotton to partial root-zone irrigation in the oasis field of Northwest China. Agric. Water Manag. 84, 41–52. doi: 10.1016/j.agwat.2006.01.010
Du, T., Kang, S., Zhang, J., Li, F., and Yan, B. (2008). Water use efficiency and fruit quality of table grape under alternate partial root-zone drip irrigation. Agric. Water Manag. 95, 659–668. doi: 10.1016/j.agwat.2008.01.017
Du, S. N., Liang, Y., Zhai, S., Xu, F. L., Chen, Z. J., and Chen, Z. H. (2005). Effect of irrigation style on eggplant growth. Chin. Agricult. Sci. Bull. 21, 430–432. doi: 10.3969/j.issn.1000-6850.2005.06.124
Fang, D. P., Zhang, F. C., Li, J., Wang, H. D., Xiang, Y. Z., and Zhang, Y. (2015). Effects of irrigation amount and various fertigation methods on yield and quality of cucumber in greenhouse. J. Appl. Ecol. 26, 1735–1742. doi: 10.13287/j.1001-9332.20150331.016
Faoro, H., Alves, A. C., Souza, E. M., Rigo, L. U., Cruz, L. M., Al-Janabi, S. M., et al. (2010). Influence of soil characteristics on the diversity of bacteria in the Southern Brazilian Atlantic Forest. Appl. Environ. Microbiol. 76, 4744–4749. doi: 10.1128/AEM.03025-09
Frey, S. D., Knorr, M., Parrent, J. L., and Simpson, R. T. (2004). Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests. For. Ecol. Manag. 196, 159–171. doi: 10.1016/j.foreco.2004.03.018
Fuchslueger, L., Bahn, M., Fritz, K., Hasibeder, R., and Richter, A. (2014). Experimental drought reduces the transfer of recently fixed plant carbon to soil microbes and alters the bacterial community composition in a mountain meadow. New Phytol. 201, 916–927. doi: 10.1111/nph.12569
Fuchslueger, L., Bahn, M., Hasibeder, R., Kienzl, S., Fritz, K., Schmitt, M., et al. (2016). Drought history affects grassland plant and microbial carbon turnover during and after a subsequent drought event. J. Ecol. 104, 1453–1465. doi: 10.1111/1365-2745.12593
Genre, A., Lanfranco, L., Perotto, S., and Bonfante, P. (2020). Unique and common traits in mycorrhizal symbioses. Nat. Rev. Microbiol. 18, 649–660. doi: 10.1038/s41579-020-0402-3
Girvan Martina, S., Bullimore, J., Pretty Jules, N., Osborn, A. M., and Ball Andrew, S. (2003). Soil type is the primary determinant of the composition of the total and active bacterial communities in arable soils. Appl. Environ. Microbiol. 69, 1800–1809. doi: 10.1128/AEM.69.3.1800-1809.2003
Gutiérrez-Miceli, F. A., Santiago-Borraz, J., Montes Molina, J. A., Nafate, C. C., Abud-Archila, M., Oliva Llaven, M. A., et al. (2007). Vermicompost as a soil supplement to improve growth, yield and fruit quality of tomato (Lycopersicum Esculentum). Bioresour. Technol. 98, 2781–2786. doi: 10.1016/j.biortech.2006.02.032
Hartmann, M., Brunner, I., Hagedorn, F., Bardgett, R. D., Stierli, B., Herzog, C., et al. (2017). A decade of irrigation transforms the soil microbiome of a semi-arid pine Forest. Mol. Ecol. 26, 1190–1206. doi: 10.1111/mec.13995
Hashmi, I., Bindschedler, S., and Junier, P. (2020). “Chapter 18 – Firmicutes,” in Beneficial Microbes in Agro-Ecology. eds. N. Amaresan, M. Senthil Kumar, K. Annapurna, K. Kumar, and A. Sankaranarayanan (Amsterdam: Academic Press).
Hazrati, S., Tahmasebi-Sarvestani, Z., Modarres-Sanavy, S. A. M., Mokhtassi-Bidgoli, A., and Nicola, S. (2016). Effects of water stress and light intensity on chlorophyll fluorescence parameters and pigments of Aloe vera L. Plant Physiol. Biochem. 106, 141–148. doi: 10.1016/j.plaphy.2016.04.046
Hirabayashi, Y., Mahendran, R., Koirala, S., Konoshima, L., Yamazaki, D., Watanabe, S., et al. (2013). Global flood risk under climate change. Nat. Clim. Chang. 3, 816–821. doi: 10.1038/nclimate1911
Hwang, C., Correll, M. J., Gezan, S. A., Zhang, L., Bhakta, M., Vallejos, C. E., et al. (2017). Next generation crop models: a modular approach to model early vegetative and reproductive development of the common bean (Phaseolus Vulgaris L.). Agric. Syst. 155, 225–239. doi: 10.1016/j.agsy.2016.10.010
Kang, S., and Zhang, J. (2004). Controlled alternate partial root-zone irrigation: its physiological consequences and impact on water use efficiency. J. Exp. Bot. 55, 2437–2446. doi: 10.1093/jxb/erh249
Kõljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F. S., Bahram, M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481
Kralj Kunčič, M., Kogej, T., Drobne, D., and Gunde-Cimerman, N. (2010). Morphological response of the halophilic fungal genus Wallemia to high salinity. Appl. Environ. Microbiol. 76, 329–337. doi: 10.1128/AEM.02318-09
Liu, H., Brettell, L. E., Qiu, Z., and Singh, B. K. (2020). Microbiome-mediated stress resistance in plants. Trends Plant Sci. 25, 733–743. doi: 10.1016/j.tplants.2020.03.014
Liu, J., Shu, A., Song, W., Shi, W., Li, M., Zhang, W., et al. (2021). Long-term organic fertilizer substitution increases rice yield by improving soil properties and regulating soil bacteria. Geoderma 404:115287. doi: 10.1016/j.geoderma.2021.115287
Louca, S., Parfrey Laura, W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507
Luo, H. H., Zhang, Y. L., and Zhang, W. F. (2016). Effects of water stress and rewatering on photosynthesis, root activity, and yield of cotton with drip irrigation under mulch. Photosynthetica 54, 65–73. doi: 10.1007/s11099-015-0165-7
Macfarlane, S. A. (2003). Molecular determinants of the transmission of plant viruses by nematodes. Mol. Plant Pathol. 4, 211–215. doi: 10.1046/j.1364-3703.2003.00164.x
Mathur, S., Tomar, R. S., and Jajoo, A. (2019). Arbuscular mycorrhizal fungi (Amf) protects photosynthetic apparatus of wheat under drought stress. Photosynth. Res. 139, 227–238. doi: 10.1007/s11120-018-0538-4
Morales, F., Anc, N. M., Fakhet, D., Gonz Lez-Torralba, J. G., Mez, A. L., Seminario, A., et al. (2020). Photosynthetic metabolism under stressful growth conditions as a bases for crop breeding and yield improvement. Plants 9:88. doi: 10.3390/plants9010088
Naylor, D., and Coleman-Derr, D. (2018). Drought stress and root-associated bacterial communities. Frontiers in. Plant Sci. 8:2223. doi: 10.3389/fpls.2017.02223
Pérez-Álvarez, E. P., Intrigliolo, D. S., Almajano, M. P., Rubio-Bret, N. P., and Garde-Cerd, N. T. (2021). Effects of water deficit irrigation on phenolic composition and antioxidant activity of monastrell grapes under semiarid conditions. Antioxidants 10:1301. doi: 10.3390/antiox10081301
Philippot, L., Raaijmakers, J. M., Lemanceau, P., and Van Der Putten, W. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Pretty, J., and Bharucha, Z. P. (2014). Sustainable intensification in agricultural systems. Ann. Bot. 114, 1571–1596. doi: 10.1093/aob/mcu205
Puértolas, J., Pardos, M., De Ollas, C., Albacete, A., and Dodd, I. C. (2020). Soil moisture heterogeneity regulates water use in Populus Nigra L. by altering root and xylem sap Phytohormone concentrations. Tree Physiol. 40, 762–773. doi: 10.1093/treephys/tpaa037
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
Rajabi-Khamseh, S., Danesh-Shahraki, A., and Rafieiolhossaini, M. (2020). Stress tolerance in flax plants inoculated with bacillus and Azotobacter species under deficit irrigation. Physiol. Plant 170, 269–279.
Rolli, E., Marasco, R., Vigani, G., Ettoumi, B., Mapelli, F., Deangelis, M. L., et al. (2015). Improved plant resistance to drought is promoted by the root-associated microbiome as a water stress-dependent trait. Environ. Microbiol. 17, 316–331. doi: 10.1111/1462-2920.12439
Romero, P., Dodd, I. C., and Martinez-Cutillas, A. (2012). Contrasting physiological effects of partial root zone drying in field-grown grapevine (Vitis vinifera L. Cv. Monastrell) according to Total soil water availability. J. Exp. Bot. 63, 4071–4083. doi: 10.1093/jxb/ers088
Rouphael, Y., Kyriacou, M. C., Petropoulos, S. A., De Pascale, S., and Colla, G. (2018). Improving vegetable quality in controlled environments. Sci. Hortic. 234, 275–289. doi: 10.1016/j.scienta.2018.02.033
Rousk, J., B Th, E., Brookes, P. C., Lauber, C. L., Lozupone, C., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a Ph gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Rovira, A. D. (1965). Interactions between plant roots and soil microorganisms. Annu. Rev. Microbiol. 19, 241–266. doi: 10.1146/annurev.mi.19.100165.001325
Saijo, Y., and Loo, E. P.-I. (2020). Plant immunity in signal integration between biotic and abiotic stress responses. New Phytol. 225, 87–104. doi: 10.1111/nph.15989
Saiz-Rubio, V., and Rovira Más, F. (2020). From smart farming towards agriculture 5.0: A review on crop data management. Agronomy 10:207. doi: 10.3390/agronomy10020207
Sandhya, V., Shaik, Z. A., Grover, M., Reddy, G., and Venkateswarlu, B. (2009). Alleviation of drought stress effects in sunflower seedlings by the exopolysaccharides producing pseudomonas putida strain gap-P45. Biol. Fertil. Soils 46, 17–26. doi: 10.1007/s00374-009-0401-z
Santos, D. L., Coelho, E. F., De Oliveira, R. A., Filgueiras, R., Da Silva Alves, M., Rodrigues, W. P., et al. (2021). Impact of soil water regimes and partial root-zone drying in field-grown papaya in semi-arid conditions. Sci. Rep. 11:10638.
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Shahnazari, A., Liu, F., Andersen, M. N., Jacobsen, S.-E., and Jensen, C. R. (2007). Effects of partial root-zone drying on yield, tuber size and water use efficiency in potato under field conditions. Field Crop Res. 100, 117–124. doi: 10.1016/j.fcr.2006.05.010
Song, T., Das, D., Hu, Q., Yang, F., and Zhang, J. (2021). Alternate wetting and drying irrigation and phosphorus rates affect grain yield and quality and heavy metal accumulation in Rice. Sci. Total Environ. 752:141862. doi: 10.1016/j.scitotenv.2020.141862
Song, Y., and Huang, C. (1987). Study on measuring method of eggplant leaf area. Northern Horticulture. 4, 23–28.
Tabatabai, M. A., and Bremner, J. M. (1969). Use of P-Nitrophenyl phosphate for assay of soil phosphatase activity. Soil Biol. Biochem. 1, 301–307. doi: 10.1016/0038-0717(69)90012-1
Timmusk, S., Kim, S.-B., Nevo, E., Abd El Daim, I., Ek, B., Bergquist, J., et al. (2015). Sfp-type PPtase inactivation promotes bacterial biofilm formation and ability to enhance wheat drought tolerance. Front. Microbiol. 6:387. doi: 10.3389/fmicb.2015.00387
Verbon, E. H., and Liberman, L. M. (2016). Beneficial microbes affect endogenous mechanisms controlling root development. Trends Plant Sci. 21, 218–229. doi: 10.1016/j.tplants.2016.01.013
Voesenek, L. A. C. J., and Bailey-Serres, J. (2015). Flood adaptive traits and processes: an overview. New Phytol. 206, 57–73. doi: 10.1111/nph.13209
Vohník, M., Sadowsky, J. J., Kohout, P., Lhot Kov, Z., Nestby, R., and Kolařík, M. (2012). Novel root-fungus symbiosis in ericaceae: sheathed ericoid mycorrhiza formed by a hitherto undescribed basidiomycete with affinities to trechisporales. PLoS One 7:E39524. doi: 10.1371/journal.pone.0039524
Wang, J., Li, Y., and Niu, W. (2020). Deficit alternate drip irrigation increased root-soil-plant interaction, tomato yield, and quality. Int. J. Environ. Res. Public Health 17:781.
Watson, D. J. (1947). Comparative physiological studies on the growth of field crops: I. Variation in net assimilation rate and leaf area between species and varieties, and within and between years. Ann. Bot. 11, 41–76.
Williams, A., and De Vries, F. T. (2020). Plant root exudation under drought: implications for ecosystem functioning. New Phytol. 225, 1899–1905. doi: 10.1111/nph.16223
Wu, A., Hammer, G. L., Doherty, A., Von Caemmerer, S., and Farquhar, G. D. (2019). Quantifying impacts of enhancing photosynthesis on crop yield. Nature Plants 5, 380–388. doi: 10.1038/s41477-019-0398-8
Wyszkowska, J., and Wyszkowski, M. (2010). Activity of soil dehydrogenases, urease, and acid and alkaline phosphatases in soil polluted with petroleum. J. Toxic. Environ. Health A 73, 1202–1210. doi: 10.1080/15287394.2010.492004
Xie, F., Ma, A., Zhou, H., Liang, Y., Yin, J., Ma, K., et al. (2020). Revealing fungal communities in alpine wetlands through species diversity, functional diversity and ecological network diversity. Microorganisms 8:632. doi: 10.3390/microorganisms8050632
Yuan, B. Z., Sun, J., and Nishiyama, S. (2004). Effect of drip irrigation on strawberry growth and yield inside a plastic greenhouse. Biosyst. Eng. 87, 237–245. doi: 10.1016/j.biosystemseng.2003.10.014
Zaborowska, M., Wyszkowska, J., and Kucharski, J. (2019). Soil enzyme response to Bisphenol F contamination in the soil bioaugmented using bacterial and mould fungal consortium. Environ. Monit. Assess. 192:20. doi: 10.1007/s10661-019-7999-6
Zalar, P., Sybren De Hoog, G., Schroers, H.-J., Frank, J. M., and Gunde-Cimerman, N. (2005). Taxonomy and phylogeny of the xerophilic genus Wallemia (Wallemiomycetes and Wallemiales, Cl. Et Ord. Nov.). Antonie Van Leeuwenhoek 87, 311–328. doi: 10.1007/s10482-004-6783-x
Zaller, J. G. (2007). Vermicompost As A substitute for peat in potting media: effects on germination, biomass allocation, yields and fruit quality of three tomato varieties. Sci. Hortic. 112, 191–199. doi: 10.1016/j.scienta.2006.12.023
Zhang, H., Li, Y., and Zhu, J.-K. (2018). Developing naturally stress-resistant crops for a sustainable agriculture. Nature Plants 4, 989–996. doi: 10.1038/s41477-018-0309-4
Zhang, Y., Zhao, Y., Sun, L., Han, P., Bai, X., Lin, R., et al. (2021). The N uptake-associated physiological processes at late growth stage in wheat (Triticum aestivum) under N deprivation combined with deficit irrigation condition. Plant Physiol. Biochem. 164, 160–172. doi: 10.1016/j.plaphy.2021.04.023
Zhao, S., Shi, G., and Dong, X. (2002). Plant physiology experiment guide. Beijing: China Agriculture Press.
Zhou, H., Zhang, D., Jiang, Z., Sun, P., Xiao, H., Yuxin, W., et al. (2019). Changes in the soil microbial communities of alpine steppe at Qinghai-Tibetan plateau under different degradation levels. Sci. Total Environ. 651, 2281–2291. doi: 10.1016/j.scitotenv.2018.09.336
Zivcak, M., Brestic, M., Balatova, Z., Drevenakova, P., Olsovska, K., Kalaji, H. M., et al. (2013). Photosynthetic electron transport and specific photoprotective responses in wheat leaves under drought stress. Photosynth. Res. 117, 529–546. doi: 10.1007/s11120-013-9885-3
Zivcak, M., Brestic, M., Kalaji, H. M., and Govindjee, (2014). Photosynthetic responses of sun-and shade-grown barley leaves to high light: is the lower Psii connectivity in shade leaves associated with protection against excess of light? Photosynth. Res. 119, 339–354. doi: 10.1007/s11120-014-9969-8
Keywords: irrigation amount, microbial community, soil physicochemical factors, yield, eggplant
Citation: Ji T, Guo X, Wu F, Wei M, Li J, Ji P, Wang N and Yang F (2022) Proper irrigation amount for eggplant cultivation in a solar greenhouse improved plant growth, fruit quality and yield by influencing the soil microbial community and rhizosphere environment. Front. Microbiol. 13:981288. doi: 10.3389/fmicb.2022.981288
Edited by:
Marco Nuti, Sant’Anna School of Advanced Studies, ItalyReviewed by:
Roberto Mancinelli, University of Tuscia, ItalySaurav Das, University of Nebraska-Lincoln, United States
Copyright © 2022 Ji, Guo, Wu, Wei, Li, Ji, Wang and Yang. 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: Ningxin Wang, bnh3YW5nQHNkYXUuZWR1LmNu; Fengjuan Yang, Zmp5YW5nQHNkYXUuZWR1LmNu
†These authors have contributed equally to this work