- 1College of Life Sciences, Sichuan Normal University, Chengdu, China
- 2Leshan Haitang Experimental Middle School, Leshan, China
Drought restricts the growth of alpine grassland vegetation. This study aimed to explore a new technical system to improve the drought resistance of forage grass. Qinghai cold-land Poa pratensis seedlings were used in the drought stress experiment. A combination of abscisic acid (ABA) and polyacrylamide (PAM) were used to affect the growth, leaf physiology, soil enzyme activity, and rhizosphere microbial diversity of P. pratensis. The fresh leaf weight and root surface area were significantly increased after ABA-PAM combined treatment, while root length was significantly reduced. Besides, the leaf catalase (CAT) and superoxide dismutase (SOD) enzyme activity, proline and chlorophyll content, increased after the treatment, while malondialdehyde (MDA) content decreased. The treatment also increased sucrase, urease, and alkaline protease activities in rhizosphere soil, while decreasing acid phosphatase and neutral phosphatase enzyme activities. ABA-PAM combined treatment enhanced the rhizosphere microbial community and forage drought resistance by altering the abundance of various dominant microorganisms in the rhizosphere soil. The relative abundances of Actinobacteria, Chloroflexi, and Acidobacteria decreased, while Proteobacteria, Firmicutes, and Ascomycota increased. Unlike the relative abundance of Gibberella that decreased significantly, Komagataeibacter, Lactobacillus, Pichia, and Dekkera were significantly increased. Single-factor collinearity network analysis revealed a close relationship between the different rhizosphere microbial communities of forage grass, after ABA-PAM treatment. This study implies that ABA-PAM combined treatment can improve the drought resistance of forages. Therefore, it provides a theoretical and practical basis for restoring drought-induced grassland degradation.
Introduction
Grassland is one of the most widely distributed terrestrial ecosystems. It regulates ecosystem balance, biodiversity, regional economy, climate change, and animal husbandry development (Du et al., 2004; Curlock and Hall, 2010; Zheng et al., 2018). The total grassland area in China is about 3.9 × 108 hm2, accounting for 13% of the global grassland area and about 41% of the land area in China (Kang et al., 2007). However, interference by natural factors and human activities in the recent decades has severely destroyed grasslands (Gang et al., 2014; Yang et al., 2016). Grassland degradation often leads to soil desertification, with large pores in the sandy soil particles, loose soil texture, high permeability, poor water retention, and low water content. These culminate into grassland desertification, ecological deterioration, and even sandstorms. The northwestern part of Sichuan Province is a high-latitude, high-altitude alpine grassland ecosystem. Besides, it is an important part of the plateau ecosystem in the entire Qinghai-Tibet area and the Three River Source region (Wen et al., 2013). Drought and cold stress have affected the vegetation growth in this area for a long time due to the high altitude and scarce precipitation, thus decreasing biodiversity and biomass (Gang et al., 2014; Wang et al., 2016; Yang et al., 2016). These have consequently resulted in a low robustness and resilience. As a result, nearly 38.8% of the grassland has been degraded (Xia et al., 2014; Wang et al., 2016).
Soil provides nutrients, water, and root-fixing substrate for plant growth. Several studies have shown that plant-soil-microbe interaction is crucial in maintaining the balance of the terrestrial ecosystems. Soil microorganisms, including bacteria, fungi, actinomycetes, algae, and protozoa, are essential in the soil ecosystem. Active microorganisms in the rhizosphere micro-ecological environment provide good nutritional conditions for plants. They also directly or indirectly affect soil nutrient cycling, energy flow, etc. Therefore, microorganisms can be used as a key parameter for evaluating soil environmental quality (Koide et al., 2011; Pratscher et al., 2011; Eisenhauer et al., 2012; Petersen et al., 2012; Purahong and Krüger, 2012). However, rhizosphere microorganisms are sensitive to the subtle changes in the soil micro-ecological environment, including drought stress (Berendsen et al., 2012; Mendes et al., 2013). Drought stress disturbs the normal physiological and biochemical reactions of plants, such as the water content, net photosynthetic rate, transpiration rate and stomatal conductance of leaves (Lafitte et al., 2007; Chaves et al., 2009; Li et al., 2015; Yuan et al., 2016). Drought stress also directly causes the lysis and death of some microbial cells (Turner et al., 2003) and changes the physical structure, chemical properties, nutrient status, element stoichiometric ratio, and enzyme activity of soil (Cregger et al., 2014; Wang et al., 2018). Moreover, drought stress indirectly changes the use of carbon sources by microorganisms by affecting the photosynthesis, respiration and other metabolic processes of plants (Lemanceau et al., 2006).
Polyacrylamide (PAM) is a high molecular polymer polymerized by Acrylamide (AMD) and its derivatives. It can be used as a water-retaining agent to reduce soil evaporation. This enhances soil water holding capacity (Wen et al., 2013), thus increasing the microbial biomass carbon, nitrogen mass fraction, and crop yield (Li et al., 2017). Chen found that PAM as a water-retaining agent, increases the soil wettability of artificial soil (Chen Z. et al., 2018). Researchers have also found that soaking millet seeds with an appropriate concentration of PAM effectively reduces the drought stress damage to millet seed germination and seedling growth, thus enhancing plant resistance to drought (Ke et al., 2015). Young also found that PAM increases the water storage capacity during drought period, significantly enhancing the drought tolerance of high roof greening plants (Young et al., 2017). Abscisic acid (ABA) is an endogenous hormone that responds to biotic and abiotic stresses in plants, thus can significantly improve plant resistance to drought. ABA regulates plant growth, inhibits seed germination, and increases plant resistance to stress (Apel and Hirt, 2004; De Smet et al., 2006; Wasilewska et al., 2008; Ding and De Smet, 2013). ABA can induce the formation of protective enzymes in the biomembrane system, reduce the degree of membrane lipid peroxidation, and protect the integrity of the membrane structure and photosynthetic characteristics. These are initiated by increasing the content of proline and soluble sugars and controlling the closing of stomata to reduce leaf water evaporation (Liu et al., 2013; Li et al., 2014; Minardi et al., 2014; Kaur et al., 2015; Han et al., 2016; Wu and Liang, 2017; Wang et al., 2020). Although studies have reported the effects of ABA on crop resistance (Lu et al., 2003; Souza et al., 2013; Chae et al., 2015; Estrada-Melo et al., 2015) and the restoration of degraded grassland using water-retaining agents (Yuan, 2010; Chu et al., 2012), none has reported the effect of abscisic acid-polyacrylamide (ABA-PAM) on the restoration of degraded grassland vegetation and soil microbial diversity.
This study systematically investigated the effects of ABA-PAM combined treatment on plant growth and the rhizosphere microbial diversity under drought stress conditions. The study aimed to improve the adaptability of forages to arid environments. Therefore, this study may provide a theoretical and practical basis for restoring degraded grassland and managing undegraded grassland via ABA-PAM combined treatment.
Materials and methods
Experimental materials
This study used a local grass species, P. pratensis, a native perennial herb mainly distributed in Qinghai Tibet Plateau. The severe desertified soil was collected from heavy-desertification grasslands in Xiaman Town, Zoige County, Sichuan Province, at an altitude of 3,500 m. The soil had a 90.7% sediment content. The area is cold and windy in winter and spring; and rainy and humid in summer and autumn. It has an average annual temperature and rainfall of 1.4°C and 643.3 mm, respectively.
Forage planting and treatment
A 40-mesh sieve was used to remove impurities from the severe desertified soil. Each square flowerpot (13 × 13 × 15 cm) was filled with 2.5 kg of desertified soil. The pots were divided into four groups (six pots per group). Only groups 1 (ABA-PAM) and 2 (PAM) were treated with PAM (8 gm⋅–2 PAM mixed with 0.5 kg upper soil). The healthy and plump P. pratensis seeds were evenly mixed into the upper soil of the flowerpot (0.3 g per flowerpot), watered (100 mL), and cultivated in a light incubator at 25 °C at a 16/8 h photoperiod. The height of the forage grass plants was measured every seven days. Each pot was watered (100 mL) every five days after sowing. Drought stress was applied when the forage had grown to a height of 5 cm (seven days after sowing) and watered (100 mL) every seven days. ABA (2 mg⋅L–1) was sprayed (thrice) on the leaves at 450 L⋅hm–2 after 12 days of stress treatment for groups 1 and 3 (ABA), once every 12 days, for a total of 3 sprays. Groups 2 and 4 (CK) received equal volumes of distilled water. Soil water content were measured with WET-2 Soil Water Sensor (Delta-T Devices Ltd., United Kingdom) at different time point.
Determination of leaf and soil enzyme activity
After the treatment period (63 days), the whole plant was uprooted, washed, and the fresh weight, root weight, root length, and plant height were measured. A root scanner (WinRHIZO STD4800 LA2400, Regent, Canada) was used to determine the root surface area. Forage leaves and soil samples were collected to determine enzyme activity. Leaf catalase, peroxidase, superoxide dismutase, malondialdehyde, H2O2, proline, and chlorophyll contents were determined using the corresponding kits (A007-1-1, A084-3-1, A001-1-2, A003-1-2, A064-1-1, A107-1-1, and A147-1-1, Nanjing Jiancheng Bio-Engineering Institute Co., Ltd., China). The activity of soil alkaline protease (ALPT), urease (UE), neutral phosphatase (NP), acid phosphatase (ACP), polyphenol oxidase (PPO), and sucrase (SC) were determined using soil enzyme activity kits (BC0885, BC0125, BC0465, BC0145, BC0115, and BC0245, Beijing Solarbio Science & Technology Co., Ltd., China).
Microbial diversity analysis
Forage plants were planted in 48 pots and divided into eight groups. As mentioned above, the groups were treated with drought stress, PAM, and ABA, and sampled at different time points. Correspondingly, eight sample groups were collected, including severe desertified soil grew grass without PAM and ABA treatment (CK); severe desertified soil grew grass with drought (D10, D40); severe desertified soil grew grass with PAM treatment, without (P0) or with drought (PD10, PD40); severe desertified soil grew grass, with drought and ABA treatment (AD40); severe desertified soil grew grass, with drought, PAM, and ABA treatment (PAD40) (Table 1). The whole plants were uprooted and rhizosphere soil samples were collected after the treatment period (63 days). The soil was divided into three parts. One part was stored in a sterile EP tube (three replicates for each sample), another part was stored at –80°C after quick freezing in liquid nitrogen, and the last part was sent to Shanghai Majorbio Bio-pham Technology Co., Ltd. Powersoil DNA Isolation Kit (MO-BIO, United States) was used to extract the total soil DNA. NanoDrop 2000 (NanoDrop, United States) and 1% agarose gel electrophoresis were used to determine the purity, concentration, and integrity of DNA samples. The DNA samples were then amplified with FastPfu polymerase (TransGen, China). The primers 338F (5′-ACTCCTACGGGAGGCAGCAG -3′) and 806R (5′-GGACTACHVGGGTWTCTAAT -3′) were used to amplify the bacterial 16S rRNA gene V3-V4 region. ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′) were used to amplify the fungal ITS1-ITS2 region. Each sample was conducted in replicate with 20 μL PCR reaction system, including 4 μL 5 × FastPfu Buffer, 2 μL 2.5 mM dNTPs, 0.8 μL Forward Primer (5 μM), 0.8 μL Reverse Primer (5 μM), 0.4 μL FastPfu Polymerase, 0.2 μL BSA, 10 ng Template DNA, added ddH2O to 20 μL. The PCR program for 16S rRNA gene included initial denaturation at 95°C for 3 min, followed by 30 cycles of 95°C for 30 s, 56°C for 30 s and 72°C for 45 s, and a final extension at 72°C for 10 min. The PCR condition for ITS was 95°C for 3 min, followed by 30 cycles of 95°C for 30 s, 55°C for 30 s and 72°C for 45 s, and a final extension at 72°C for 10 min. The amplicons were then purified with AxyPrep DNA Gel Extraction Kit (Axygen, United States) and quantified with Quantus™ Fluorometer (Promega, United States). Sequencing libraries were prepared with NEXTflex® Rapid DNA-Seq Kit (Bioo Scientific, United States). Illumina MiSeq 2500 platform (Illumina, United States) was used for pair-end sequencing of the amplified library. The sequencing data were uploaded to the Majorbio Cloud platform1 for analysis. The sequencing reads were deposited into the NCBI Sequence Read Archive (SRA) database under the accession number PRJNA805202. Raw fastq files were quality-filtered by Fastp (v 0.19.2) (Chen S. et al., 2018) and assemblies by FLASH (v1.2.113) (Magoč and Salzberg, 2011). Operational taxonomic units (OTUs) were clustered with a 97% similarity cut-off using UPARSE (v7.0.10904) (Edgar, 2013) using a novel ‘greedy’ algorithm that simultaneously performs chimera filtering. Then the OTUs were classified using RDP classifier (v2.115) (Wang et al., 2007) against the Silva 16S rRNA database (v1386) and UNITE ITS database (v8.07) using a confidence threshold of 70% to obtain the species classification information.
Data processing and analysis
Excel 2010 was used to pre-process and analyze the test data. SPSS 19.0 was used to analyze variance and test of significance of difference (LSD method). Origin 2018 software was used for mapping based on the values calculated by SPSS. The community column chart was drawn using the R language (version 3.3.1) according to the community composition data at different taxonomic levels. The Kruskal Wallis h test and one-way ANOVA were used to evaluate the significance of the observed differences and draw the relevant histogram based on the community abundance data in the sample. Network software was used for single factor correlation network analysis.
Results
Effect of abscisic acid-polyacrylamide combined treatment on plant growth
In preliminary experiments, we firstly added 60 and 90% sand to the culture substrate, and then added 0, 4, 8, and 16 g⋅m–2 PAM. The water contents were analyzed. When added 8 g⋅m–2 PAM into the soil and the soil was treated with drought stress for 80 days, the water content was 3.18 ± 0.06% and 10.99 ± 0.97%, respectively, corresponding to 60 and 90% sediment content. When added 0, 4, and 16 g⋅m–2 PAM to the soil with 90% sediment content, the water contents were 1.48 ± 0.22%, 2.16 ± 0.24%, and 1.88 ± 0.23%, while that for 60% sediment content were 8.89 ± 0.55%, 8.96 ± 0.61%, and 9.69 ± 1.08%, respectively. These results suggested that 8 g⋅m–2 was the optimal PAM ratio. So, we selected 8 g⋅m–2 PAM in this following study.
The results showed that the growth situation increased after ABA-PAM treatment under drought stress compared with the control, ABA, and PAM groups (Figure 1A). On the 28th day, the average plant heights in the ABA-PAM, ABA, PAM, and control groups were 5.99 ± 0.08, 5.85 ± 0.10, 6.09 ± 0.04 cm, and 5.83 ± 0.43 cm, respectively (Figure 1B). Moreover, the fresh leaf weight and root surface area were significantly higher in the ABA-PAM treatment group (0.38 ± 0.01 g and 0.39 ± 0.02 cm2, respectively) than in the control group (0.10 ± 0.04 and 0.30 ± 0.09, respectively). However, root length was significantly lower in the ABA-PAM group (5.11 ± 0.21 cm) than in the control group (7.34 ± 0.09 cm) (Figure 1C).
Figure 1. Effect of ABA-PAM combined treatment on the growth of grasses. Drought-stressed forage grasses were treated with ABA and/or PAM (A), and plant height were measured at different time points (B). After the treatment period (63 days), the whole plants were uprooted and washed. The fresh weight, root weight, and root length were measured (C). Leaf catalase (CAT), peroxidase (POD), superoxide dismutase (SOD), malondialdehyde (MDA), H2O2, proline, and chlorophyll contents were determined (D,E). The units for MDA, H2O2, proline, and chlorophyll were mmol⋅g–1 protein, mmol⋅g–1 protein, μg⋅g–1, and mg⋅g–1, respectively.
Water contents were measured when the forage grass plants were treated with drought stress. Results showed that the soil water contents were increased after PAM treatment (Supplementary Table 1). Soil water contents for PAM group were higher than the control group at all time points. However, when the grass plants were sprayed with ABA, the soil water content were lower than that in the control group, which could be because the ABA treatment group had higher transpiration rate. ABA-PAM combined treatment not only promoted the growth of grass plants, but also increased the soil water contents than the control, although the improvement rate were lower than the PAM treatment group. Besides, the soil water contents of drought-stressed grass plants were lower than 8% at most time points, and lower than 5% on every seventh day, which means moderate and severe drought, respectively.
Effect of abscisic acid-polyacrylamide combined treatment on leaf physiology
Antioxidative enzymes, such as catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD), are well known to be involved in protecting plants against oxidative stress. CAT, POD, and SOD contents in the leaves of the drought-stressed P. pratensis were 14.73 ± 0.42, 58.31 ± 2.94, and 145.66 ± 0.85 U⋅mg–1 protein, respectively, in the ABA-PAM group, and 9.51 ± 1.19, 41.19 ± 0.74, and 101.06 ± 2.11 U⋅mg–1 protein, respectively, in the control group (Figure 1D). The contents of toxic substances (malondialdehyde and H2O2) were significantly lower in the ABA-PAM group (1.33 ± 0.12 and 19.26 ± 2.06 mmol⋅g–1 protein, respectively) than in control (2.35 ± 0.07 and 31.94 ± 2.99 mmol⋅g–1 protein, respectively) and the PAM groups (1.78 ± 0.08, 26.51 ± 0.40 mmol⋅g–1 protein, respectively). The MDA levels were significantly lower in the ABA-PAM treatment group (1.87 ± 0.001 mmol⋅g–1 protein) than in the ABA group. Besides, the content of the two toxic substances were also lower in ABA groups than in control (Figure 1E). However, proline and chlorophyll content were elevated in the ABA-PAM (301.95 ± 3.30 μg⋅g–1, 6.65 ± 0.26 mg⋅g–1) and ABA (293.43 ± 3.33 μg⋅g–1, 6.49 ± 0.35 mg⋅g–1) group than in the PAM (273.08 ± 5.78 μg⋅g–1, 5.47 ± 0.27 mg⋅g–1) and control (257.40 ± 6.60 μg⋅g–1, 4.52 ± 0.23 mg⋅g–1) group (Figure 1E). The above results indicate that ABA-PAM combined treatment enhances antioxidant activity, reduces the content of toxic substances, degradation of chlorophyll, and the damage of drought stress, thus providing adequate nutrition guarantee for the growth and development of plants.
Effects of abscisic acid-polyacrylamide combined treatment on soil enzyme activities
Abscisic acid (ABA) and PAM significantly increased soil invertase (SC) activity. However, SC activity was highest in the ABA-PAM treatment group (123.22 ± 1.05 U⋅g–1, Figure 2A). The acid phosphatase (ACP) activity was inhibited after PAM treatment (4.36 ± 0.64 nmoldg–1) and ABA-PAM treatment (4.73 ± 0.30 nmoldg–1). In contrast, ABA alone had no significant effect on ACP activity (7.17 ± 1.35 nmoldg–1). Neutral phosphatase (NP) activity was inhibited in ABA and ABA-PAM treatment groups (1.43 ± 0.06 and 1.07 ± 0.06 nmoldg–1, respectively) (Figure 2B). However, PAM treatment did not significantly affect NP activity (1.94 ± 0.09 nmoldg–1). Moreover, ABA and PAM had no significant effect on polyphenol oxidase (PPO) activity (Figure 2A). Urease (UE) activity was highest in the ABA-PAM treatment group, reaching 0.038 ± 0.001 U⋅g–1 (Figure 2B). PAM enhanced alkaline protease (ALPT) activity (0.85 ± 0.05 U⋅g–1). However, ABA-PAM significantly enhanced ALPT activity (1.34 ± 0.10 U⋅g–1) (Figure 2B).
Figure 2. Effect of ABA-PAM combined treatment on soil invertase (SC), acid phosphatase (ACP), neutral phosphatase (NP), polyphenol oxidase (PPO), urease (UE) and protease (ALPT). After the trea tment period (63 days), soil samples were collected for the analyses of soil enzyme activities. (A) Activities of SC, PPO, and ACP; (B) Activities of ALPT, NP, and UE.
Effects of abscisic acid-polyacrylamide combined treatment on microbial diversity under drought stress
16S rDNA/ITS sequencing analysis
Plant rhizosphere soil microbiome is essential for plant growth and development, nutrient acquisition, and adaptation to stress (Philippot et al., 2013; Trivedi et al., 2020). Soil samples were collected from the alpine degraded grassland in Xiaman Town, Zoige County, Sichuan Province. P. pratensis plants were grown and treated with drought stress, followed by the sequencing of the 16S rDNA/ITS (CK, P0, D10, PD10, D40, AD40, PD40, PAD40, see Table 1 for sample treatment methods, see part 2.2 for soil treat and P. pratensis planting methods) to analyze the effect of ABA and PAM on the soil microbial community structure of forage rhizosphere. A total of 1,527,099 16S rDNA reads and 1,620,295 ITS DNA reads were obtained (average lengths, 415 and 221 bp, respectively).
A total of 4,812 bacterial OTUs and 2,278 fungal OTUs were obtained using 97% sequence similarity as the classification threshold for OTUs clustering (Supplementary Table 2). The bacterial OTUs were distributed in 41 phyla, 106 classes, 256 orders, 453 families, 925 genera, and 1,802 species, while fungal OTUs were distributed in eight phyla, 31 classes, 92 orders, 205 families, 446 genera, and 750 species. Alpha diversity analysis found that PAM and ABA alone had no significant effect on the bacterial Chao1 index and Shannon index under drought stress compared with the control group (Figures 3A,B and Supplementary Table 3). In contrast, ABA-PAM treatment decreased the bacterial Chao1 index (3152.69), Shannon index (5.3562), and bacterial diversity. Moreover, PAM, ABA, and ABA-PAM had no significant effect on the Chao1 and fungal Shannon indices compared with the control group (Figures 3C,D and Supplementary Table 3).
Figure 3. Effects of PAM and ABA on soil microbial alpha diversity. (A) Chao index of bacterial communities; (B) Shannon index of bacterial communities; (C) Chao index of fungal communities; (D) Shannon index of fungal communities. CK: severe desertified soil grew grass without PAM and ABA treatment; D10, D40: severe desertified soil grew grass with 10 or 40 days drought; P0: severe desertified soil grew grass with PAM but without drought; PD10, PD40: severe desertified soil grew grass with 10 or 40 days drought; AD40: severe desertified soil grew grass with 40 days drought and ABA treatment; PAD40: severe desertified soil grew grass, with 40 days drought, PAM, and ABA treatment. *p < 0.05 and **p < 0.01.
Effect of polyacrylamide on microbial community composition
Analysis of six samples without ABA (CK, P0, D10, PD10, D40, and PD40) found that the relative abundance of nine phyla bacteria was more than 1%. Actinobacteria had the highest abundance, followed by Proteobacteria, Chloroflex, Acidobacteria, Firmicutes, Gemmatimonadetes, Bacteroidetes, Rokuhacteria, and Planctomycetes. The number of sequences of the five phyla with the highest abundance accounted for more than 89.4% of the total sequences and was the dominant flora in the bacterial community. The relative abundance of Actinobacteria, Chloroflexi, Acidobacteria, Gemmatimonadetes, Rokuhacteria, and Planctomycetes in soil was decreased by 2.1%, 1.25%, 5.99%, 0.16%, 0.17%, and 0.51%, respectively, in the absence of drought stress. PAM increased the relative abundance of other dominant bacteria with an abundance > 1% (P0 vs. CK). Moreover, the relative abundances of Actinobacteria, Chloroflexi, Rokuhacteria, and Planctomycetes in PAM-treated soil were decreased by 2.63%, 1.94%, 0.28%, and 0.43%, respectively, after 10 days of drought stress (PD10 vs. D10). The relative abundances of Actinobacteria, Chloroflexi, Acidobacteria, Gemmatimonadetes, Rokuhacteria and Planctomycetes in PAM treatment group was decreased by 9.80%, 3.89%, 3.61%, 1.16%, 0.42% and 0.14%, respectively, after 40 days of drought treatment (PD40 vs. D40) (Supplementary Figure 1A and Supplementary Table 4). In conclusion, PAM increased the relative abundance of dominant bacteria Proteobacteria, Firmicutes, and Bacteroidetes, while it decreased Actinobacteria, Chloroflex, and Acidobacteria.
Fungal community structure analysis (Supplementary Figure 1C and Supplementary Table 4) found that the relative abundance of three phyla in six samples was greater than 1% (Ascomycota, Zygomycota, and Basidiomycota). Ascomycota accounted for more than 82.4% of the total sequences and was the dominant flora of the fungal community. PAM reduced the relative abundance of Ascomycota and Basidiomycota by 1.29% and 3.77%, respectively, in the absence of drought stress (P0 vs. CK). Moreover, the relative abundance of Basidiomycota in the PAM treatment group was decreased by 0.02% after 10 days of drought stress treatment (PD10 vs. D10). PAM reduced the relative abundance of Ascomycota and Basidiomycota by 5.64% and 0.71%, respectively, after 40 days of drought treatment (PD40 vs. D40). In conclusion, although PAM decreased the relative abundance of Ascomycota and Basidiomycota, Ascomycota was the dominant bacteria in rhizosphere soil.
Further analysis showed that Arthrobacter (2.34%) and RB41 (2.69%) were the dominant genera in the control group (CK) in the absence of drought stress. Komagataeibacter (4.59%), Lactobacillus (3.25%), Arthrobacter (3.01%) were the dominant genera in the PAM (P0)- treated soil (Supplementary Figure 1B and Supplementary Table 4). However, PAM did not affect the dominant genera in rhizosphere soil (Komagataeibacter and Lactobacillus) after 10 days of drought stress treatment, but it altered the relative abundance. The relative abundances of Komagataeibacter and Lactobacillus in group D10 were 4.84% and 3.57%, respectively, and 5.20% and 4.20%, respectively, in group PD10. Arthrobacter (3.74%) and RB41 (2.32%) were the dominant genera in the D40 group after 40 days of drought stress treatment, while Komagataeibacter (8.80%) and Lactobacillus (6.63%) were the dominant genera in the PD40 group. The top 10 genera with significant differences in average relative abundance were screened (Supplementary Figure 2A). The results showed that PAM significantly increased the relative abundance of Komagataeibacter and Lactobacillus, making them the dominant bacteria, and thus improving the synthesis efficiency of cellulose and the conversion efficiency of nitrogen (Chen et al., 2016).
Legthophora (19.97%), Gibberella (18.16%), Pichia (8.90%), and Dekkera (7.21%) were the dominant fungi in the non-PAM-treated soil in the absence of drought stress treatment. Besides, Gibberella (20.29%), Mortierella (5.58%), and Fusarlum (4.84%) were the dominant genera in PAM-treated soil. Gibberella (27.33%), Fusarlum (16.89%), and Mortierella (6.52%) were the dominant genera in the rhizosphere soil of the D10 group after 10 days of drought stress treatment. Gibberella (30.08%), Pichia (11.59%), Dekkera (7.89%), and Mortierella (6.76%) were the dominant genera in the rhizosphere soil of the PD10 group. Furthermore, Gibberella (19.15%), Pichia (8.60%), Mortierella (7.12%), and Dekkera (6.61%) were the dominant genera in the rhizosphere soil of the D40 group after 40 days of drought stress treatment. Gibberella (32.31%), Mortierella (11.88%), and Fusarium (7.75%) were the dominant genera in the rhizosphere soil of the PD40 group. These results suggest that PAM increased the relative abundance of Gibberella and Mortierella. The top 10 groups with the highest relative abundance and significant differences were also screened at the genus level. Fusarium had the highest relative abundance. Moreover, the relative abundance of Fusarium and Myrothecium was increased significantly in group D10 (Supplementary Figure 2B).
Effect of abscisic acid on microbial community composition
Analysis of D40, AD40, PD40, and PAD40 showed that the relative abundance of dominant bacterial groups was decreased by 2.98% (Actinobacteria), 2.15% (Chloroflexi), and 2.55% (Acidobacteria) in non-PAM-treated soil (AD40) when compared with D40. After ABA-PAM treatment, the relative abundance of dominant bacterial groups was significantly reduced in PAD40 compared with D40, AD40, and PD40. In contrast, ABA-PAM increased the relative abundance of dominant populations Proteobacteria and Firmicutes (Figure 4A and Supplementary Table 5).
Figure 4. Relative abundance of the rhizosphere microbial communities of ABA treated samples at phylum and genus level. (A) Bacteria phylum level; (B) bacteria genus level; (C) fungi phylum level; (D) fungi genus level. D40: severe desertified soil grew grass with 40 days drought; PD40: severe desertified soil grew grass with 40 days drought; AD40: severe desertified soil grew grass with 40 days drought and ABA treatment; PAD40: severe desertified soil grew grass, with 40 days drought, PAM, and ABA treatment.
The relative abundance of Komagataeibacter and Lactobacillus in rhizosphere soil was significantly increased after ABA (AD40) and ABA-PAM (PAD40) treatments at the genus level compared with the control group D40, making them the dominant genera (Figure 4B). The top 10 genera with significant differences in average relative abundance were also screened (Figure 5). The relative abundance of Komagataeibacter and Lactobacillus was higher in ABA-PAM (PAD40) and ABA alone (AD40) groups than in the control group (PD40). Besides, the relative abundance of Komagataeibacter and Lactobacillus was higher in the ABA-PAM group than in the ABA alone group (Figure 5A). The effect of ABA on fungi is shown in Figure 4C. ABA (AD40) significantly decreased the relative abundance of Ascomycota by 7.15%, while it increased the relative abundance of other phyla compared with D40, AD40, and PD40.
Figure 5. Community differences of bacteria (A) and fungi (B) at the genera level caused by ABA (95% confidence interval). D40: severe desertified soil grew grass with 40 days drought; PD40: severe desertified soil grew grass with 40 days drought; AD40: severe desertified soil grew grass with 40 days drought and ABA treatment; PAD40: severe desertified soil grew grass, with 40 days drought, PAM, and ABA treatment.
Abscisic acid (ABA) significantly increased the relative abundance of Gibberella (32.45%), Mortierella (8.88%), and Fusarium (5.99%) compared with the control group (PD40), thus making them the dominant bacteria (Figure 4D). Pichia (22.8%), Dekkera (16.61%), and Gibberella (15.60%) were the dominant genera after ABA-PAM (PAD40) treatment. ABA-PAM significantly increased the relative abundance of Pichia and Dekkera, while it decreased that of Gibberella. Among the 10 genera with the highest average relative abundance and significant difference, ABA-PAM decreased the relative abundance of Fusarium (Figure 5B).
Single-factor collinear network analysis
A species correlation network map was used to reveal the relationship between species and reflect the overall trend. The 50 bacterial genera with the highest total abundance were analyzed using a single-factor correlation network. Actinobacteria and Proteobacteria accounted for 35.4 and 25%, respectively, of all nodes without PAM addition. However, PAM altered the proportions to 37.1 and 34.3%, respectively (Supplementary Figure 3). Actinobacteria and Proteobacteria accounted for 43.5 and 26%, respectively, without ABA addition, and the whole network diagram was divided into two modules. There were no corresponding values of “Diameter” and “Average shortest path length,” and the network had no connectivity. However, ABA treatment changed the proportions to 38.8 and 28.6%, respectively (Figure 6). There was a connection between the two separated modules, the network is more stable, indicating that ABA enhanced the connectivity between bacterial genera, thus improving the drought stress resistance. Actinobacteria formed the largest proportion of all nodes under the control group and ABA-PAM treatments, accounting for 42% of all network nodes, and was the dominant phylum in the treatments (Supplementary Figure 4). The number of nodes with Node Degree of 15–30 in the control group was 20, and that with Node Degree of 16–31 in the ABA-PAM group was 32, which indicated that the connectivity among bacterial genera in ABA-PAM treatment was higher, the network was more complex and stable, and the interaction among microorganisms was stronger.
Figure 6. Effect of ABA-PAM on the collinear network of soil microorganism. (A,B) Bacteria; (C,D) fungi; (A,C) ABA-PAM absent; (B,D) ABA-PAM added.
A single-factor correlation network analysis was performed on 50 fungal genera with the highest total abundance. Ascomycota accounted for the largest proportion of all nodes in control and ABA-PAM groups, accounting for 88% of all network nodes, and was the dominant microflora in the two treatments. However, ABA-PAM treatment had a more complex and stronger inter-species interaction correlation networks. The number of nodes with Node Degree 16–24 in the control group was 10, while that with Node Degree of 16–24 in the ABA-PAM treatment group was 23. These results indicate that ABA-PAM treatment enhances linkage of fungi in the microflora; making the network more complex and stable, thus enhancing the resistance from the interaction between microorganisms to drought stress.
Discussion
The northwest plateau of Sichuan is located in the eastern part of the Qinghai-Tibet Plateau. It is an essential water conservation area in the upper reaches of the Yangtze and Yellow Rivers, with a special ecological niche and important ecological functions. The Zoige grassland is located in the northwest plateau of Sichuan Province and has scarce rainfall. Human activities and global warming have caused drought in the grassland ecosystem, rapidly increasing desertification, thus affecting the ecological environment and economic development of the area (Meng et al., 2013). Therefore, improving soil water holding and the resistance of surface plants can effectively alleviate the degradation of high-cold grasslands.
A water retention agent is a high molecular polymer made of strong absorbent resin with ultra-high water absorption ability. It is insoluble in water, can quickly absorb a large amount of water, thus enhancing soil water holding capacity. Besides, it slowly releases its absorbed water when water shortage occurs in its surroundings, thus improving soil physical and chemical structure and increasing water utilization rate (Schröfl et al., 2012). In recent years, polyacrylamide has been a popular soil water retention agent due to its stable performance, long water retention cycle, natural decomposition without residues, and low cost. As one of the six phytohormones, ABA regulates plant seed germination, growth and development, and fruit ripening. Besides, it promotes resistance to stress caused by abiotic factors (Apel and Hirt, 2004; Wasilewska et al., 2008; Umezawa et al., 2010). Plants can enhance resistance under adversity stresses by increasing the activity of antioxidant enzymes, increasing endogenous ABA synthesis, and reducing MDA levels. Such strategies protect against plant damage (Atici et al., 2005). ABA can significantly improve the expression of various transcription factors and genes associated with signaling pathways, such as salicylic acid, jasmonic acid, and ethylene, significantly increasing the content of endogenous ABA to improve plant stress resistance (Xiong and Zhu, 2003; Souza et al., 2013; Wang et al., 2013). In addition, ABA enhances plant resistance by closing plant pores to reduce vaporization and maintain moisture in plants (Jiang and Zhang, 2002; Castro Neto, 2003; Peeva and Cornic, 2009). ABA also regulates the expression of many genes associated with plant dehydration tolerance, such as late-embryonic redundant proteins with high hydrophilicity, and can protect biomolecules and maintain moisture in dehydrated cells (Jiang and Zhang, 2002).
A strong root system is essential for plant growth and development. The growth of the root system is closely related to soil moisture conditions. ABA-PAM combined treatment increased the soil moisture content under drought stress. As a result, plants did not seek water through root elongation, explaining the shorter root length and increased root surface area and water absorption in the ABA-PAM group than in the control group. ABA also increased the metabolic rate and height of plants compared with the control group. Exogenous ABA enhances ABA synthesis in the plant and the redistribution within the plant under drought stress to promote the closure of the stoma, thereby reducing the loss of water (Wang et al., 2011) and significantly increasing the fresh weight of the leaves.
Adversity stress induces several reactive oxygen species (ROS) production, causing a series of harmful biochemical reactions. Plants form an antioxidant protection system to remove ROS in the long-term evolution to prevent ROS damage in the body, thus improving plant resistance (Anjum et al., 2017). Studies have shown that ABA enhances the activity of phenylalanine ammonialyase (PAL), superoxide dismutase (SOD), polyphenol oxidase (PPO), and peroxidase (POD) in tomatoes with Alternaria solani (Song et al., 2011). CAT, POD, and SOD activities in the leaves were significantly increased after ABA-PAM treatment under drought stress. Therefore, ABA-PAM treatment improves ROS removal in plants, protects enzymes involved in anabolic correlation in the body, and improves drought stress resistance. Malondialdehyde (MDA) is a membrane fat peroxidation product that can indirectly reflect the degree of damage to a plant cell membrane. Increased MDA content indicates an increased degree of membrane damage. ABA-PAM treatment significantly reduced MDA and H2O2 content compared with the control group, indicating that ABA-PAM treatment reduces the damage caused by toxic substances to plant cells. Plants quickly synthesize various osmotic regulatory substances, including proline and betaine, to resist osmotic stress when affected by abiotic factors, such as drought and high salt (Yoshiba et al., 1997). Proline is the most widely distributed osmotic regulatory substance (Zhang et al., 1997; Knipp and Honermeier, 2006; Singh and Singh, 2006; Vendruscolo et al., 2007; Dobra et al., 2010). Several studies have shown that proline removes ROS either by itself or by stimulating CAT, PPO, SOD, and POD activities (Han, 2006). Proline works in synergy with antioxidant systems, such as CAT, POD, SOD, to regulate ROS balance in cells (Rodriguez and Redman, 2005). Proline accumulation is positively correlated with plant tolerance to abiotic stress (Knipp and Honermeier, 2006; Slama et al., 2006; Trovato et al., 2008). In this study, ABA-PAM treatment promoted the accumulation of grass proline and enhanced plant resistance to drought.
Soil enzymes are mainly involved in the circulation of matter and energy and are susceptible to environmental factors. Soil sucrase can break down sucrose into glucose and fructose for plants and microorganisms and is an important indicator of soil fertility (Pancholy and Rice, 1973). ABA-PAM combined treatment increased SC activity, enhancing soil fertility. Phosphatase can hydrolyze organophosphorus in soil into inorganic phosphorus and can be used as an indicator of soil phosphorus conversion and plant phosphorus absorption (Wang, 2008; Nannipieri et al., 2011). PAM and ABA-PAM decreased ACP activity. However, ABA had no significant effect on ACP activity. ABA and ABA-PAM decreased NP activity, while PAM had no significant effect on NP activity. Therefore, PAM inhibits ACP activity, while ABA inhibits NP activity. Phosphorus can improve the root ratio surface area of the plant, reduce the respiratory rate of the root system, and increase the potential of root water under drought stress. This enhances the absorption and utilization of water and nutrients of the root system, improves the moisture condition in the plant, and increases the resistance to plant stresses (Zhang and Li, 1996). The combined action of ABA-PAM enhances soil moisture content and plant drought resistance. Plants do not need extra phosphorus to resist adversity, resulting in reduced phosphatase activity and the breakdown of organophosphorus in the soil. The activity of soil UE, as a key enzyme in urea hydrolysis, reflects the inorganic nitrogen supply capacity of the soil. Soil proteases are involved in protein hydrolysis in soil, providing nitrogen sources for plant growth (Jian et al., 2016). ABA-PAM significantly increased UE and ALPT activities. ABA and PAM had different effects on the activities of various enzymes in this study. Besides, ABA-PAM treatment enhances soil enzyme activity to a certain extent, improves nutrient supply and conversion levels in root systems under drought stress, and promotes forage grass growth.
The high Chao1 index and Shannon index indicate more species in the sample and higher species diversity of the sample, respectively (Grice et al., 2009). PAM and ABA had no significant effect on the Chao1 and Shannon indexes of bacteria. However, ABA-PAM treatment decreased bacteria’s Chao1 and Shannon indexes. ABA-PAM treatment reduced the richness of bacteria and the diversity of species in the root-soil. Moreover, ABA-PAM, PAM, and ABA had no significant effect on Chao1 and Shannon indexes of fungi, fungal diversity, and richness. Therefore, these results show that ABA and PAM changed the physical and chemical properties of the soil, promoting plant growth, and thus affecting the production of root secretions, thus decreasing bacterial abundance and diversity. The spearman correlation coefficients were also calculated between the physiological characteristics and the rhizosphere dominant microorganisms of sample D40, AD40, PD40, and PAD40, results showed that most of the abundance changes of rhizosphere dominant microorganisms were correlated to the changes of physiological characteristics (Supplementary Figure 5).
Increased Proteobacteria can effectively fix nitrogen sources. Some members of Actinobacteria regulate symbiotic nitrogen fixation and phosphorus solution (He et al., 2009; Tu et al., 2015; Cui et al., 2018; Zhao et al., 2020), Chloroflexi can produce energy through photosynthesis but cannot regulate nitrogen fixation. Moreover, Chloroflexi has a potential phosphorus biolysis function (Kragelund et al., 2007; Yu et al., 2016). Bacteroidetes have a strong metabolic capacity for complex organic compounds, proteins, and lipids (Hill et al., 2007). ABA and PAM made the Actinobacteria, Proteobacteria, Chloroflexi, Acidobacteria, Firmicutes dominant in the root-soil. This improved the soil environment, maintained the balance of root nutrient absorption and microenvironment, and enhanced the drought stress resistance of the plants. The increase in the relative abundance of Firmicutes and Proteobacteria may have inhibited the growth of Actinobacteria, Chloroflexi, and Acidobacteria, decreasing their relative abundance. The dominance of Komagataeibacter and Fusarium gradually declined as drought stress time increased, while Komagataeibacter abundance significantly increased in PAM and ABA groups. These results indicate that drought causes the decline of Komagataeibacter and Fusarium, while PAM and ABA increase the relative abundance of Komagataeibacter and Fusarium in the soil.
Rhizosphere soil fungi were classified into eight phyla, 31 classes, 92 orders, and 205 families. Ascomycota, Zygomycota, and Basidiomycota were the dominant groups in soil fungi. Basidiomycetes and ascomycetes play different roles in the decomposition of soil organic matter. For instance, basidiomycetes degrade the refractory lignocellulose in plant litter, a dominant compound in the soil with high lignin content. However, ascomycetes mainly decomposes the rotten complex organic matter in the soil (Lienhard et al., 2014; De Araujo et al., 2017; Yao et al., 2017). Therefore, drought stress can increase the lignin content of plants, change the quality of litter input into the soil, increase the difficulty of degradation, and increase the relative abundance of ascomycetes (García-Palacios et al., 2016; Risch et al., 2019), explaining why ascomycetes had the highest relative abundance in all fungal communities. Although the addition of ABA and PAM reduced the drought stress effect, ascomycetes was still the dominant class in the fungal community. Moreover, the addition of ABA, PAM or ABA-PAM can increase the relative abundance of fungi in rhizosphere soil at the phylum and genus levels. ABA and PAM treatment mainly altered the structure instead of rhizosphere soil fungal community composition.
Wang simulated soil drought stress through pot water control experiment and found that nitrogen-fixing bacteria can decrease under prolonged drought stress, indicating that drought stress can reduce the content of nitrogen-fixing bacteria in the soil (Wang et al., 2005). ABA and PAM treatment increased the relative abundance of nitrogen-fixing bacteria Proteus in the soil. Proteus can effectively fix nitrogen and slow down the impact of drought stress on plants. Some members of actinomycetes have symbiotic nitrogen fixation and phosphorus removal function (He et al., 2009; Tu et al., 2015; Cui et al., 2018; Zhao et al., 2020). Actinomycetes have high tolerance to drought (Zheng, 2021). This study revealed that Curvularia viridis can produce energy through photosynthesis, but cannot fix nitrogen. Moreover, it has good biological phosphorus hydrolysis ability (Kragelund et al., 2007; Yu et al., 2016). Bacteroidetes have a strong metabolic capacity for complex organics, proteins and lipids (Hill et al., 2007). ABA and PAM combined treatment increased the dominance of Actinomycetes, Proteus, Curvularia, Acidobacteria, and Firmicutes in rhizosphere soil. This enhanced soil phosphorus solubilization and nitrogen-fixation, improved soil environment, and maintained the balance of root nutrient absorption and microenvironment, thus enhancing plant resistance to drought stress.
Conclusion
When the Poa pratensis seedlings were treated with drought and ABA-PAM simultaneously, the fresh leaf weight and root surface area were increased, while root length was reduced. The leaf enzyme activity of CAT and SOD, and content of proline and chlorophyll increased, MDA content decreased. The sucrase, urease, and alkaline protease activities in rhizosphere soil increased, while acid phosphatase and neutral phosphatase enzyme activities decreased. ABA-PAM combined treatment improved the rhizosphere microbial community and forage drought resistance by altering the abundance of various dominant microorganisms in the rhizosphere soil, resulted in a more complex and stable network between different rhizosphere microbial communities of forage grass. All these results imply that ABA-PAM combined treatment can improve the drought resistance of forages.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA805202.
Author contributions
XTao and BY conceived this study, designed the experimental plan, analyzed data, and drafted and revised the manuscript. XTang, XF, and YS participated in sample preparing, treating, data analyses, and manuscript draft and revise. HS participated in data analyses and manuscript revise. JZ, XH, and XW participated in data analyses. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (2017YFC0504903), the Project of Science and Technology in Chengdu (2016-HM01-00046-SF), and the Scientific Research Foundation of Sichuan Normal University (024-341906001).
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/fpls.2022.973665/full#supplementary-material
Supplementary Figure 1 | Relative abundance of the rhizosphere microbial communities of PAM treated samples at phylum and genus level. (A) Bacteria phylum level; (B) bacteria genus level; (C) fungi phylum level; (D) fungi genus level. CK: severe desertified soil grew grass without PAM and ABA treatment; D10, D40: severe desertified soil grew grass with 10 or 40 days drought; P0: severe desertified soil grew grass with PAM but without drought; PD10, PD40: severe desertified soil grew grass with 10 or 40 days drought.
Supplementary Figure 2 | Community differences of bacteria (A) and fungi (B) at the genera level caused by PAM. CK: severe desertified soil grew grass without PAM and ABA treatment; D10, D40: severe desertified soil grew grass with 10 or 40 days drought; P0: severe desertified soil grew grass with PAM but without drought; PD10, PD40: severe desertified soil grew grass with 10 or 40 days drought.
Supplementary Figure 3 | Effect of PAM on the collinear network of soil microorganism. (A,B) Bacteria; (C,D) fungi; (A,C) PAM absent; (B,D) PAM added.
Supplementary Figure 4 | Effect of ABA on the collinear network of soil microorganism. (A,B) Bacteria; (C,D) fungi; (A,C) ABA absent; (B,D) ABA added.
Supplementary Figure 5 | Correlation analysis between physiological characteristics and rhizosphere microbes. Spearman correlation coefficients were calculated between the physiological characteristics and rhizosphere dominant microorganisms for sample D40, AD40, PD40, and PAD40.
Footnotes
- ^ https://cloud.majorbio.com/
- ^ https://github.com/OpenGene/fastp
- ^ http://www.cbcb.umd.edu/software/flash
- ^ http://www.drive5.com/uparse/
- ^ https://sourceforge.net/projects/rdp-classifier/
- ^ https://www.arb-silva.de/
- ^ https://unite.ut.ee/
References
Anjum, S. A., Ashraf, U., Tanveer, M., Khan, I., Hussain, S., Shahzad, B., et al. (2017). Drought induced changes in growth, osmolyte accumulation and antioxidant metabolism of three maize hybrids. Front. Plant Sci. 8:69. doi: 10.3389/fpls.2017.00069
Apel, K., and Hirt, H. (2004). Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Rev. Plant Biol. 55, 373–399. doi: 10.1146/annurev.arplant.55.031903.141701
Atici, Ö., Ağar, G., and Battal, P. (2005). Changes in phytohormone contents in chickpea seeds germinating under lead or zinc stress. Biol. Plant. 49, 215–222. doi: 10.1007/s10535-005-5222-9
Berendsen, R. L., Pieterse, C. M., and Bakker, P. A. H. M. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. doi: 10.1016/j.tplants.2012.04.001
Castro Neto, M. T. (2003). Effect of water deficit on the transpiration and stomatal resistance of mango tree. Rev. Bras. Frutic. 25, 93–95. doi: 10.1590/S0100-29452003000100027
Chae, W. L., Woonhee, B., Jangho, J., Jung-Hyun, K., and Sung, C. L. (2015). Function of ABA in stomatal defense against biotic and drought stresses. Int. J. Mol. Sci. 16, 15251–15270. doi: 10.3390/ijms160715251
Chaves, M. M., Flexas, J., and Pinheiro, C. (2009). Photosynthesis under drought and salt stress: Regulation mechanisms from whole plant to cell. Ann. Bot. 103, 551–560. doi: 10.1093/aob/mcn125
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Chen, H. M., Li, C. F., Bi, J. C., Zhang, C., Hu, Q. S., Wang, M. Y., et al. (2016). Study on optimal carbon sources of Komagataeibacter spp. for bacterial cellulose synthesis. Food Mach. 32, 37–41. doi: 10.13652/j.issn.1003-5788.2016.04.009
Chen, Z., Wang, R., Han, P., Sun, H., Sun, H., Li, C., et al. (2018). Soil water repellency of the artificial soil and natural soil in rocky slopes as affected by the drought stress and polyacrylamide. Sci. Total Environ. 619, 401–409. doi: 10.1016/j.scitotenv.2017.11.146
Chu, L. M., Yan, G. E., Hui, J. Y. J. S., and China, W. C. (2012). Effect of polyacrylamide on plant recovery of degraded grassland. Soil Water Conserv. China 2012:05. doi: 10.14123/j.cnki.swcc.2012.05.001
Cregger, M. A., McDowell, N. G., Pangle, R. E., Pockman, W. T., and Classen, A. T. (2014). The impact of precipitation change on nitrogen cycling in a semi-arid ecosystem. Funct. Ecol. 28, 1534–1544. doi: 10.1111/1365-2435.12282
Cui, J. T., Li, Y. A., Wang, C. Y., Kim, K. S., Wang, T. Y., and Liu, S. X. (2018). Characteristics of the rhizosphere bacterial community across different cultivation years in saline–alkaline paddy soils of Songnen Plain of China. Can. J. Microbiol. 64, 925–936. doi: 10.1139/cjm-2017-0752
Curlock, J. M. O., and Hall, D. O. (2010). The global carbon sink: A grassland perspective. Glob. Change Biol. 4, 229–233. doi: 10.1046/j.1365-2486.1998.00151.x
De Araujo, A. S. F., Bezerra, W. M., Dos Santos, V. M., Nunes, L. A. P. L., De Lyra, M. D. C. P., Figueiredo, M. D. B., et al. (2017). Fungal diversity in soils across a gradient of preserved Brazilian Cerrado. J. Microbiol. 55, 273–279. doi: 10.1007/s12275-017-6350-6
De Smet, I., Zhang, H., Inze, D., and Beeckman, T. (2006). A novel role for abscisic acid emerges from underground. Trends Plant Sci. 11, 434–439. doi: 10.1016/j.tplants.2006.07.003
Ding, Z., and De Smet, I. (2013). Localised ABA signalling mediates root growth plasticity. Trends Plant Sci. 18, 533–535. doi: 10.1016/j.tplants.2013.08.009
Dobra, J., Motyka, V., Dobrev, P., Malbeck, J., Prasil, I. T., Haisel, D., et al. (2010). Comparison of hormonal responses to heat, drought and combined stress in tobacco plants with elevated proline content. Plant Physiol. 167, 1360–1370. doi: 10.1016/j.jplph.2010.05.013
Du, M. Y., Kawashima, S., Yonemura, S., Zhang, X. Z., and Chen, S. B. (2004). Mutual influence between human activities and climate change in the Tibetan Plateau during recent years. Glob. Planet Change 41, 241–249. doi: 10.1016/j.gloplacha.2004.01.010
Edgar, R. C. (2013). UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Eisenhauer, N., Scheu, S., and Jousset, A. (2012). Bacterial diversity stabilizes community productivity. PLoS One 7:e34517. doi: 10.1371/journal.pone.0034517
Estrada-Melo, A. C., Ma, C., Reid, M. S., and Jiang, C. Z. (2015). Overexpression of an ABA biosynthesis gene using a stress-inducible promoter enhances drought resistance in petunia. Hortic. Res. 2:15013. doi: 10.1038/hortres.2015.13
Gang, C. C., Chen, Y. Z., Wang, Z. Q., Sun, Z. G., Li, J. L., Qi, J. G., et al. (2014). Quantitative assessment of the contributions of climate change and human activities on global grassland degradation. Carrpath. J. Earth Environ. 72, 4273–4282. doi: 10.1007/s12665-014-3322-6
García-Palacios, P., Prieto, I., Ourcival, J.-M., and Hättenschwiler, S. (2016). Disentangling the litter quality and soil microbial contribution to leaf and fine root litter decomposition responses to reduced rainfall. Ecosystems 19, 490–503. doi: 10.1007/s10021-015-9946-x
Grice, E. A., Kong, H. H., Conlan, S., Deming, C. B., Davis, J., Young, A. C., et al. (2009). Topographical and temporal diversity of the human skin microbiome. Science 324, 1190–1192. doi: 10.1126/science.1171700
Han, X. L. (2006). In vitro selection of L-hydroxyproline resistant variants of Coronilla varia land studies on their salt tolerance. Xi’an: Northwestern University. doi: 10.7666/d.y894256
Han, Y. S., Wang, S. J., Zhao, N., Deng, S. R., Zhao, C. J., Li, N. F., et al. (2016). Exogenous abscisic acid alleviates cadmium toxicity by restricting Cd2+ influx in Populus euphratica cells. Plant Growth Regul. 35, 827–837. doi: 10.1007/s00344-016-9585-2
He, C. N., Gao, W. W., Yang, J. X., Bi, W., Zhang, X. S., and Zhao, Y. J. (2009). Identification of autotoxic compounds from fibrous roots of Panax quinquefolium L. Plant Soil 318, 63–72. doi: 10.1007/s11104-008-9817-8
Hill, V. R., Kahler, A. M., Jothikumar, N., Johnson, T. B., Hahn, D., and Cromeans, T. L. (2007). Multistate evaluation of an ultrafiltration-based procedure for simultaneous recovery of enteric microbes in 100-liter tap water samples. Appl. Environ. Microbiol. 73, 4218–4225. doi: 10.1128/AEM.02713-06
Jian, S. Y., Li, J. W., Chen, J., Wang, G. S., Mayes, M. A., Dzantor, K. E., et al. (2016). Soil extracellular enzyme activities, soil carbon and nitrogen storage under nitrogen fertilization: A meta-analysis. Soil Biol. Biochem. 101, 32–43. doi: 10.1016/j.soilbio.2016.07.003
Jiang, M., and Zhang, J. H. (2002). Water stress-induced abscisic acid accumulation triggers the increased generation of reactive oxygen species and up-regulates the activities of antioxidant enzymes in maize leaves. J. Exp. Bot. 53, 2401–2410. doi: 10.1093/jxb/erf090
Kang, L., Han, X. G., Zhang, Z. B., and Sun, O. J. (2007). Grassland ecosystems in China: Review of current knowledge and research advancement. Philos. Trans. R. Soc. 362, 997–1008. doi: 10.1098/rstb.2007.2029
Kaur, S., Gupta, A. K., and Zhawar, V. K. (2015). ABA-dependent sucrose regulation of antioxidant metabolism in wheat cultivars varying in ABA-sensitivity. Biologia 70, 165–173. doi: 10.1515/biolog-2015-0029
Ke, Z. J., Yin, M. Q., Wen, Y. Y., Huang, M. J., Huang, X. F., Guo, P. Y., et al. (2015). Effects of polyacrylamide seed soaking on seed germination and drought resistance of millet (Setaria italica) seedlings under drought stress. J. Nucl. Agr. Sci. 29:563. doi: 10.11869/j.issn.100-8551.2015.03.0563
Knipp, G., and Honermeier, B. (2006). Effect of water stress on proline accumulation of genetically modified potatoes (Solanum tuberosum L.) generating fructans. Plant Physiol. 163, 392–397. doi: 10.1016/j.jplph.2005.03.01
Koide, R. T., Fernandez, C., and Petprakob, K. (2011). General principles in the community ecology of ectomycorrhizal fungi. Ann. For. Sci. 68, 45–55. doi: 10.1007/s13595-010-0006-6
Kragelund, C., Caterina, L., Borger, A., Thelen, K., Eikelboom, D., Tandoi, V., et al. (2007). Identity, abundance and ecophysiology of filamentous Chloroflexi species present in activated sludge treatment plants. FEMS Microbiol. Ecol. 59, 671–682. doi: 10.1111/j.1574-6941.2006.00251.x
Lafitte, H., Yongsheng, G., Yan, S., and Li, Z. K. (2007). Whole plant responses, key processes, and adaptation to drought stress: The case of rice. J. Exp. Bot. 58, 169–175. doi: 10.1093/jxb/erl101
Lemanceau, P., Offre, P., Mougel, C., Gamalero, B. G., Dessaux, Y., Meönne-Loccoz, Y., et al. (2006). Plant-microbe interactions and soil quality. Microb. Ecol. Rhizosphere CABI. doi: 10.1079/9780851990989.0228
Li, Q., Ba, T., Li, Y. L., Xie, L., Liu, J. H., Yu, Z., et al. (2017). Effects of application methoths of super absorbent polymers on soil moisture content and microbial biomass and potato yield. Acta Agr. Boreali-Occidentalis Sinica. 26, 1453–1460. doi: 10.7606/j.issn.1004-1389.2017.10.006
Li, S., Zhang, L., and Yao, Y. Q. (2015). Effects of different water stress on active oxygen, stoma and photosynthesis characteristics of wheat. J. Hebei Univ. 35, 487–493. doi: 10.3969/j.issn.1000-1565.2015.05.008
Li, S. W., Leng, Y., Feng, L., and Zeng, X.-Y. (2014). Involvement of abscisic acid in regulating antioxidative defense systems and IAA-oxidase activity and improving adventitious rooting in mung bean [Vigna radiata (L.) Wilczek] seedlings under cadmium stress. Environ. Sci. Pollut. Res. Int. 21, 525–537. doi: 10.1007/s11356-013-1942-0
Lienhard, P., Terrat, S., Prevost-Boure, N. C., Nowak, V., Regnier, T., Sayphoummie, S., et al. (2014). Pyrosequencing evidences the impact of cropping on soil bacterial and fungal diversity in Laos tropical grassland. Agron. Sustain. Dev. 34, 525–533. doi: 10.1007/s13593-013-0162-9
Liu, Y. J., Ji, X. Y., Zheng, L., Nie, X. G., and Wang, Y. C. (2013). Microarray analysis of transcriptional responses to abscisic acid and salt stress in Arabidopsis thaliana. Int. J. Mol. Sci. 14, 9979–9998. doi: 10.3390/ilms14059979
Lu, S., Guo, Z., and Peng, X. X. (2003). Effects of ABA and S-3307 on drought resistance and antioxidative enzyme activity of turfgrass. J. Hortic. Sci. Biotechnol. 78, 663–666. doi: 10.1080/14620316.2003.11511681
Magoč, T., and Salzberg, S. L. (2011). FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663.
Meng, J. W., Zuo, L., Cai, Y. J., Deng, D. Z., and Yan, W. X. (2013). An analysis of desertification status and its preventive strategies in Ruoergai County. J. Sichuan For. Sci. Technol. 34, 42–46. doi: 10.16779/j.cnki.1003-5508.2013.04.012
Minardi, B. D., Voytena, A. P. L., Santos, M., and Randi, A. M. (2014). Water stress and abscisic acid treatments induce the CAM pathway in the epiphytic fern Vittaria lineata (L.) Smith. Photosynthetica 52, 404–412. doi: 10.1007/s11099-014-0047-4
Nannipieri, P., Giagnoni, L., Landi, L., and Renella, G. (2011). “Role of phosphatase enzymes in soil,” in Phosphorus in action, eds E. Bünemann, A. Oberson, and E. Frossard (Berlin: Springer), 215–243. doi: 10.1007/978-3-642-15271-9_9
Pancholy, S. K., and Rice, E. L. (1973). Soil enzymes in relation to old field succession: Amylase, cellulase, invertase, dehydrogenase, and urease. Soil. Sci. Soc. Am. J. 37, 47–50.
Peeva, V., and Cornic, G. (2009). Leaf photosynthesis of Haberlea rhodopensis before and during drought. Environ. Exp. Bot. 65, 310–318. doi: 10.1016/j.envexpbot.2008.09.009
Petersen, D. G., Blazewicz, S. J., Firestone, M., Herman, D. J., Turetsky, M., and Waldrop, M. (2012). Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska. Environ. Microbiol. 14, 993–1008. doi: 10.1111/j.1462-2920.2011.02679.x
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
Pratscher, J., Dumont, M. G., and Conrad, R. (2011). Ammonia oxidation coupled to CO2 fixation by archaea and bacteria in an agricultural soil. Proc. Natl. Acad. Sci. U.S.A. 108, 4170–4175. doi: 10.1073/pnas.1010981108
Purahong, W., and Krüger, D. (2012). A better understanding of functional roles of fungi in the decomposition process: Using precursor rRNA containing ITS regions as a marker for the active fungal community. Ann. For. Sci. 69, 659–662. doi: 10.1007/s13595-012-0210-7
Risch, A. C., Zimmermann, S., Ochoa-Hueso, R., Schütz, M., Frey, B., Firn, J. L., et al. (2019). Soil net nitrogen mineralisation across global grasslands. Nat. Commun. 10:4981.
Rodriguez, R., and Redman, R. (2005). Balancing the generation and elimination of reactive oxygen species. Proc. Natl. Acad. Sci. U.S.A. 102, 3175–3176. doi: 10.1073/pnas.0500367102
Schröfl, C., Mechtcherine, V., and Gorges, M. (2012). Relation between the molecular structure and the efficiency of superabsorbent polymers (SAP) as concrete admixture to mitigate autogenous shrinkage. Cem. Concr. Res. 42, 865–873. doi: 10.1016/j.cemconres.2012.03.011
Singh, B., and Singh, G. (2006). Effects of controlled irrigation on water potential, nitrogen uptake and biomass production in Dalbergia sissoo seedlings. Environ. Exp. Bot. 55, 209–219. doi: 10.1016/j.envexpbot.2004.11.001
Slama, I., Messedi, D., Ghnaya, T., Savoure, A., and Abdelly, C. (2006). Effects of water deficit on growth and proline metabolism in Sesuvium portulacastrum. Environ. Exp. Bot. 56, 231–238. doi: 10.1016/j.envexpbot.2005.02.007
Song, W. W., Ma, X. R., Tan, H., and Zhou, J. Y. (2011). Abscisic acid enhances resistance to Alternaria solani in tomato seedlings. Plant Physiol. Biochem. 49, 693–700. doi: 10.1016/j.plaphy.2011.03.018
Souza, T., Magalhães, P. C., Castro, E., Albuquerque, P., and Marabesi, M. A. (2013). The influence of ABA on water relation, photosynthesis parameters, and chlorophyll fluorescence under drought conditions in two maize hybrids with contrasting drought resistance. Acta Physiol. Plant. 35, 515–527. doi: 10.1007/s11738-012-1093-9
Trivedi, P., Leach, J. E., Tringe, S. G., Sa, T. M., and Singh, B. K. (2020). Plant-microbiome interactions: From community assembly to plant health. Nat. Rev. Microbiol. 18, 607–621. doi: 10.1038/s41579-020-0412-1
Trovato, M., Mattioli, R., and Costantino, P. (2008). Multiple roles of proline in plant stress tolerance and development. Rend. Lincei 19, 325–346. doi: 10.1007/s12210-008-0022-8
Tu, Y., Yang, W. Y., Liu, W. G., Yong, T. W., Jiang, L. Q., and Wang, X. C. (2015). Effects of relay strip intercropping years between flue-cured tobacco and soybean on rhizospheric microbes quantities. Agron. Sin. 41, 733–742. doi: 10.3724/SP.J.1006.2015.00733
Turner, B. L., Driessen, J. P., Haygarth, P. M., and Mckelvie, I. D. (2003). Potential contribution of lysed bacterial cells to phosphorus solubilisation in two rewetted Australian pasture soils. Soil Biol. Biochem. 35, 187–189. doi: 10.1016/S0038-0717(02)00244-4
Umezawa, T., Nakashima, K., Miyakawa, T., Kuromori, T., Tanokura, M., Shinozaki, K., et al. (2010). Molecular basis of the core regulatory network in ABA responses: Sensing, signaling and transport. Plant Cell Physiol. 51, 1821–1839. doi: 10.1093/pcp/pcq156
Vendruscolo, E. C. G., Schuster, I., Pileggi, M., Scapim, C. A., Correa Molinari, H. B., and Marur, C. J. (2007). Stress-induced synthesis of proline confers tolerance to water deficit in transgenic wheat. J. Plant Physiol. 164, 1367–1376. doi: 10.1016/j.jplph.2007.05.001
Wang, D., Rui, Y., Ding, K., Cui, X., Hao, Y., Tang, L., et al. (2018). Precipitation drives the biogeographic distribution of soil fungal community in Inner Mongolian temperate grasslands. J. Soils Sediments 18, 222–228. doi: 10.1007/s11368-017-1727-z
Wang, H., Blakeslee, J. J., Jones, M. L., Chapin, L. J., and Dami, I. E. (2020). Exogenous abscisic acid enhances physiological, metabolic, and transcriptional cold acclimation responses in greenhouse-grown grapevines. Plant Sci. 293:110437. doi: 10.1016/j.plantsci.2020.110437
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, R. S., Pandey, S., Li, S., Gookin, T. E., Zhao, Z. X., and Albert, R. (2011). Common and unique elements of the ABA-regulated transcriptome of Arabidopsis guard cells. BMC Genomics 12:216. doi: 10.1186/1471-2164-12-216
Wang, S. J. (2008). Impact of long-term chemical fertilizer application on nitrogen, phosphorus and enzyme activities of black soil. Harbin: Northeast Agricultural University. doi: 10.7666/d.y1403911
Wang, Y., Tao, X., Tang, X. M., Xiao, L., Sun, J. L., Yan, X. F., et al. (2013). Comparative transcriptome analysis of tomato (Solanum lycopersicum) in response to exogenous abscisic acid. BMC Genomics 14:841. doi: 10.1186/1471-2164-14-841
Wang, Y. X., Reng, L. H., Weng, B. Q., Ying, Z. Y., and Huang, Y. B. (2005). Ecological effects of soil drought stress on soil microorganisms of Chamaecrista rotundifolia. J. Xiamen Univ. Nat. Sci. Edn. 44:3. doi: 10.3321/j.issn:0438-0479.2005.z1.015
Wang, Z. Q., Zhang, Y. Z., Yang, Y., Zhou, W., Gang, C. C., Zhang, Y., et al. (2016). Quantitative assess the driving forces on the grassland degradation in the Qinghai-Tibet Plateau. Ecol. Inform. 33, 32–44. doi: 10.1016/j.ecoinf.2016.03.006
Wasilewska, A., Vlad, F., Sirichandra, C., Redko, Y., Jammes, F., Valon, C., et al. (2008). An update on abscisic acid signaling in plants and more. Mol. Plant 1, 198–217. doi: 10.1093/mp/ssm022
Wen, L., Dong, S. K., Li, Y. Y., Li, X. Y., Shi, J. J., Wang, Y. L., et al. (2013). Effect of degradation intensity on grassland ecosystem services in the alpine region of Qinghai-Tibetan Plateau, China. PLoS One 8:e58432. doi: 10.1371/journal.pone.0058432
Wu, X., and Liang, C. J. (2017). Enhancing tolerance of rice (Oryza sativa) to simulated acid rain by exogenous abscisic acid. Environ. Sci. Pollut. Res. 24, 4860–4870. doi: 10.1007/s11356-016-8219-3
Xia, J. Z., Liu, S. G., Liang, S. L., Chen, Y., Xu, W. F., and Yuan, W. P. (2014). Spatio-temporal patterns and climate variables controlling of biomass carbon stock of global grassland ecosystems from 1982 to 2006. Remote Sen. 6, 1783–1802. doi: 10.3390/rs6031783
Xiong, L. M., and Zhu, J. K. (2003). Regulation of abscisic acid biosynthesis. Plant Physiol. 133, 29–36. doi: 10.1104/pp.103.025395
Yang, Y., Wang, Z. Q., Li, J. L., Gang, C. C., Zhang, Y. Z., Zhang, Y., et al. (2016). Comparative assessment of grassland degradation dynamics in response to climate variation and human activities in China, Mongolia, Pakistan and Uzbekistan from 2000 to 2013. J. Arid. Environ. 135, 164–172. doi: 10.1016/j.jaridenv.2016.09.004
Yao, M., Rui, J., Niu, H., Hedìnec, P., Li, J., He, Z., et al. (2017). The differentiation of soil bacterial communities along a precipitation and temperature gradient in the eastern Inner Mongolia steppe. Catena 152, 47–56. doi: 10.1016/j.catena.2017.01.007
Yuan, R., Li, P., Hu, X. X., Zong, Y. Z., Sun, M., Dong, Q., et al. (2016). Effects of drought stress on physiological characteristics and yield of wheat. Shanxi Agric. Sci. 44, 1446–1449. doi: 10.3969/j.issn.1002-2481.2016.10.06
Yoshiba, Y., Kiyosue, T., Nakashima, K., Yamaguchi-Shinozaki, K., Shinozaki, K. J. P., and Physiology, C. (1997). Regulation of levels of proline as an osmolyte in plants under water stress. Plant Cell Physiol. 38, 1095–1102. doi: 10.1093/oxfordjournals.pcp.a029093
Young, T. M., Cameron, D. D., and Phoenix, G. K. (2017). Increasing green roof plant drought tolerance through substrate modification and the use of water retention gels. Urban Water J. 14, 551–560. doi: 10.1080/1573062X.2015.1036761
Yu, H., Si, P., Shao, W., Qiao, X., Yang, X., Gao, D., et al. (2016). Response of enzyme activities and microbial communities to soil amendment with sugar alcohols. MicrobiologyOpen 5, 604–615. doi: 10.1002/mbo3.355
Yuan, H. Y. (2010). Effects of water retaining agent on soil water content and several warm-season turf grasses seed germination. J. Anhui Agric. Sci. 38:19979. doi: 10.13989/j.cnki.0517-6611.2010.35.123
Zhang, C. S., Lu, Q., and Verma, D. P. S. (1997). Characterization of Δ1-pyrroline-5-carboxylate synthetase gene promoter in transgenic Arabidopsis thaliana subjected to water stress. Plant Sci. 129, 81–89. doi: 10.1016/S0168-9452(97)00174-x
Zhang, S. Q., and Li, Y. Y. (1996). Study on effects of fertilizing on crop yield and its mechanism to raise water use efficiency. Res. Soil Water Conserv. 3, 185–191.
Zhao, J., Xie, H. J., and Zhang, J. (2020). Microbial diversityand physicochemical properties of rhizosphere microenvironment in saline-alkali soils of the Yellow River Delta. Environ. Sci. 41, 1449–1455. doi: 10.13227/j.hjkx.201908044
Zheng, W. L. (2021). Comparative analysis of soil microbial community structure in typical degraded grassland areas in Inner Mongolia under different forbidden grazing time and precipitation. Hohhot: Inner Mongolia Agricultural University. doi: 10.27229/d.cnki.gnmnu.2021.000712
Keywords: soil bacterial community, abscisic acid, polyacrylamide, drought stress, forage grass
Citation: Tang X, Fei X, Sun Y, Shao H, Zhu J, He X, Wang X, Yong B and Tao X (2022) Abscisic acid-polyacrylamide (ABA-PAM) treatment enhances forage grass growth and soil microbial diversity under drought stress. Front. Plant Sci. 13:973665. doi: 10.3389/fpls.2022.973665
Received: 20 June 2022; Accepted: 15 August 2022;
Published: 02 September 2022.
Edited by:
Pawan Shukla, Seribiotech Research Laboratory (SBRL), IndiaReviewed by:
Subha Narayan Das, Indira Gandhi National Tribal University, IndiaChao He, Chinese Academy of Medical Sciences & Peking Union Medical College, China
Copyright © 2022 Tang, Fei, Sun, Shao, Zhu, He, Wang, Yong and Tao. 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: Xiang Tao, taoxiang@sicnu.edu.cn; Bin Yong, binyong1225@163.com
†These authors have contributed equally to this work