- 1Key Laboratory of Ecology and Environment in Minority Areas, Minzu University of China, National Ethnic Affairs Commission, Beijing, China
- 2College of Life and Environmental Sciences, Minzu University of China, Beijing, China
The response of soil microorganisms to altered nitrogen (N) and rainfall patterns plays an important role in understanding ecosystem carbon and nitrogen cycling processes under global change. Previous studies have separately focused on the effects of N addition and rainfall on soil microbial diversity and community composition. However, the combined and interactive impact of N addition and rainfall on soil microbial diversity and function mediated by plant and soil processes have been poorly investigated for grassland ecosystems. Here, we conducted a field experiment with simulated N addition (N addition: 10 g N m–2yr–1) and altered rainfall pattern [control, rainfall reduction (compared to control –50%); rainfall addition (compared to control + 50%)] to study their interactive effects on soil microbial diversity and function in a temperate steppe of Inner Mongolia. Our results showed that N addition and rainfall addition significantly increased soil bacterial diversity, and the bacterial diversity was positively correlated with soil microbial biomass nitrogen, inorganic nitrogen, and Stipa krylovii root exudate C:N ratio, Allium polyrhizum root exudate C and N, and A. polyrhizum root exudate C:N ratio. N addition and rainfall reduction significantly reduced fungal diversity, which correlated closely with soil microbial biomass carbon and the C:N ratio of A. polyrhizum root exudates. Bacteria were mainly eutrophic r-strategists, and the responses of bacterial function guilds to the interaction between N addition and rainfall pattern were not significant. However, the arbuscular mycorrhizal fungi (AMF), in the functional classification of fungi, were significantly reduced under the condition of N addition and rainfall reduction, and the absolute abundance of the phylum Glomeromycota increased under rainfall addition, suggesting that AMFs are sensitive to altered N and rainfall patterns over short timescales (1 year). Collectively, our results have important implications for understanding the plant–soil–microbe system of grasslands under climate change.
Introduction
Due to the unreasonable use of resources and the burning of fossil fuels, atmospheric nitrogen (N) deposition has increased by 2–7 fold compared to pre-industrial revolution levels (Galloway et al., 2004). N deposition increased soil nitrogen content, which can affect the growth of plants, quality of grassland forage, and nutrient cycling within ecosystems (Waldrop et al., 2004; Turner and Henry, 2010). Meanwhile, due to natural factors and human activities, global rainfall patterns have also changed at different latitudes (Allen et al., 2018). Previous studies have shown water availability is an important factor limiting the productivity of arid and semi-arid grasslands (Füzy et al., 2008); furthermore, the grassland of Inner Mongolia is an area that is more sensitive to rainfall. Therefore, a better understanding of the impact of climate change on the grassland ecosystem could be gained from studying the interaction between N addition and rainfall.
As decomposers, the diversity, composition, and function of soil microorganisms correlate closely with the biogeochemical cycles and homeostasis of ecosystems (Bardgett and Putten, 2014; Fierer, 2017; Tian et al., 2021). N deposition enhanced the photosynthetic capacity of plants (Wang et al., 2012), leading to an increase in plant biomass and the quantity of plant litter, and then acting as a sufficient nutrient source for soil microorganisms and affecting microbial diversity, community structure, and function. Previous studies have shown that N deposition increased the competition among plant species and decreased plant species diversity, thus affecting the food source and living microenvironments of soil microorganisms and resulting in the variation of soil microbial diversity (Fierer and Jackson, 2006; Fierer et al., 2007; Rousk et al., 2010). The variation of rainfall is an important factor affecting global climate change (Liu et al., 2013; Marvel and Bonfils, 2013). Rainfall addition increases soil moisture and affects the relative abundance of soil fungi and Gram-negative bacteria, resulting in altered bacteria: fungi ratios, which further impact the bacteria and fungi that participate in soil carbon and nitrogen cycling (Williams and Rice, 2007). Also, Jie et al. (2011) have documented that changes in rainfall patterns influence soil N availability. In brief, N and water are two important factors that limit soil microbial activity in the grasslands of northern China (Zhang et al., 2015). A previous study showed that precipitation was an important environmental factor affecting soil microbial carbon utilization profiles based on BIOLOG carbon substrate utilization under N addition (Sun et al., 2015). However, the functional classification of microorganisms involved in carbon cycling is still insufficient. Plants secrete a variety of root exudates during the process of growth, and these exudates act as the interaction medium between plants and microorganisms. The composition of root exudates is known to affect the soil microbial diversity and function (Baudoin et al., 2003). However, the mechanism by which the C:N ratio of root exudates affects soil microbe is still unclear. Understanding the C:N ratio of root exudates lays a foundation for explaining how plant–soil–microbial interactions play significant roles in ecosystem processes, such as carbon and nitrogen nutrient cycling.
Liang et al. (2020) have reported long-term fertilization significantly changed the composition of soil microbial community, which in turn affects the functional characteristics of the community (Fukami, 2015). It was reported that Proteobacteria, Acidobacteria, and Actinobacteria dominated the bacterial communities in the typical steppe of northern China, and that the fungal community was mainly composed of Ascomycota, Basidiomycota, and Glomeromycota (Yu et al., 2021a). Actinobacteria, Proteobacteria, and Bacteroidetes are classified as r-strategists, mainly decomposing labile C, while most Acidobacteria are classified as K-strategists and decompose recalcitrant C (Sun et al., 2021). However, it is not clear whether N addition, rainfall, and their interaction are dominated by r-strategist bacteria and accelerated soil carbon and nitrogen cycling. At present, the effects of fertilization on fungal functional groups mainly focused on arbuscular mycorrhizal fungi (AMF). AMF play important roles in plant growth, by providing multiple nutrients and assisting the host plants to resist pathogen invasion. Chen et al. (2017) have found that fertilization has decreased AMF α-diversity and altered community composition and phylogenetic patterns, suggesting that environmental factors are important driving forces of AMF community assembly. Nevertheless, the effect of interaction between N addition and rainfall on fungal functional classification is still unknown.
In order to understand how microbial diversity and function are affected by altering N and rainfall patterns, we conducted a short-term (1 year) field N addition and rainfall interaction experiment in a temperate steppe of Inner Mongolia. The main objectives of this study were (1) to assess the effects of N addition, rainfall, and their interaction on soil microbial diversity and community composition, and to confirm whether the community composition mainly belongs to r-strategists, and (2) to examine the functions of soil microorganisms in response to N addition, rainfall, and their interaction. Since bacteria are more sensitive to changes in soil moisture than fungi (Gordon et al., 2008), we hypothesized that N and rainfall interaction would (i) increase soil bacterial diversity, while not affecting fungal diversity; (ii) may increase the soil inorganic nitrogen content and soil water content and create sufficient food sources that, in turn, will increase bacterial diversity and affect bacterial community composition, mainly belonging to r-strategies; and (iii) the increased soil nutrient input will accelerate soil bacterial cycling of soil carbon and nitrogen, and enhance the functional classes of microorganisms involved in soil carbon and nitrogen cycling. Thus, the interaction between N and rainfall may increase the relative abundance of the AMF.
Materials and Methods
Site Description
We conducted a field experiment in a temperate steppe at the Hulun Lake Reserve in Inner Mongolia (113.21°E, 48.75°N), China. This site is located in a temperate zone with a mean annual temperature of –0.4°C and mean annual precipitation of 283 mm, respectively. The mean temperature was 1.8°C, and the precipitation during the growing season was 169 mm in 2018. The soil is classified as Calcic Luvisols (FAO, 1957). The vegetation is dominated by perennial species, primarily Stipa krylovii and Allium polyrhizum.
Experimental Design
The experimental design has been described in detail by Li C. et al. (2021). In brief, the experiment used a factorial, nested design with rainfall treatment [control, reduction (–50%), and addition (+ 50%)] and nitrogen (N) treatments (control and addition). Each treatment had five replicate blocks, for a total of thirty 2 m × 2 m plots that were set up on 1 June 2018. The rainfall treatments were covered by transparent rainfall shelters supported on an iron frame roof, the rainfall reduction was collected through a PVC pipe into a barrel, and rainfall addition was applied by respraying the collected rainwater from the bucket onto the treatment area. The N addition rate (in the form of NH4NO3) was 10 g N m–2 yr–1, which was applied in June, July, and August.
Field Sampling
Soil samples were collected on 10 August 2018. Five soil cores (10 cm in diameter) were randomly collected from 0 to 10 cm topsoil and then mixed to create a composite sample for each plot. The composite soil samples were sieved through a 2.0-mm sieve to remove plant debris and roots. The sieved soil samples were then separated into two parts, with one stored at 4°C for testing of physical and chemical properties and one at –80°C for high-throughput sequencing of microorganisms. The aboveground biomass was measured by cutting living plants at ground level within 0.5 m × 0.5 m quadrat. The belowground biomass was estimated by using the soil cores method (Vogt and Persson, 1991). The root exudation of S. krylovii and A. polyrhizum was collected by using a modified method in August 2018 (Phillips et al., 2008). The methods of collecting root exudates and measuring the carbon and nitrogen of root exudates were described in detail by Li C. et al. (2021).
Lab Analysis
Soil samples were oven-dried for 48 h at 105°C to obtain soil moisture content. Ten grams of air-dried soil was used to calculate the soil ammonia nitrogen and nitrate nitrogen with 50 ml of 2 M KCl, shaken for 30 min, and the filtered solution was measured by using a continuous flow analyzer (Futura, Alliance Instruments, Frépillon, France). Soil microbial biomass carbon and nitrogen contents (SMBC and SMBN, respectively) were measured using the fumigation extraction method (Vance et al., 1987). Twenty-five grams of fresh soil was extracted with 100 ml of 0.5 M K2SO4, and the extracted solutions were analyzed by a TOC analyzer (Analytik, Jena, Germany). Aboveground and belowground biomass samples were oven-dried at 65°C to a constant weight.
DNA Extraction and Polymerase Chain Reaction Amplification and Sequencing
Soil DNA was extracted from 0.5 g of fresh soil using the PowerSoil expand DNA Isolation Kit (MO BIO) according to the manufacturer’s instructions. Nanodrop Spectrophotometer (Nano-100, Aosheng Instrument Co Ltd.) assays qualified genomic DNA concentration and polymerase chain reaction (PCR) amplification products. The primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) were used to amplify the V4 region of bacteria 16S rRNA gene (Caporaso et al., 2012). In order to distinguish the different samples, a barcode sequence was added to the 5′ – end of the reverse primer. The PCR was performed in a 50-μl reaction mixture containing 20 μl of nuclease-free water, 25 μl of 2 × Primer Taq (Takara, RR902A), 1 μl of forward primer (10 mM), 1 μl of reverse primer (10 mM), and 3 μl of g-DNA. PCR cycling conditions were set for 95°C for 5 min, followed by 40 cycles of denaturation at 94°C for 30 s, annealing at 53°C for 30 s, extension at 72°C for 40 s, with a final extension of 72°C for 8 min. For fungi, the primers ITS3 (5′-GCATCGATGAAGAACGCAGC-3′) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) (Ihrmark et al., 2012) were used to amplify the ITS2 region (Gardes and Bruns, 2010). The PCR reaction mixture was the same as that used for the bacterial 16S rRNA gene. PCR cycling conditions were set at 95°C for 5 min, followed by 40 cycles of denaturation at 94°C for 30 s, annealing at 58°C for 30 s, extension at 72°C for 40 s, with a final extension of 72°C for 8 min. The PCR products were purified using the Omega kit following the manufacturer’s protocol. The purified PCR products were paired-end sequenced using Illumina HiSeq 2500 platform. We also quantified the absolute abundance of 16S rRNA and ITS genes. The reaction mixture for qPCR included 10 μl of 2 × qPCR mix, 2 μl of g-DNA, 2 μl of 2.5 μM (515F/806R and ITS3/ITS4) primers, and 6 μl of ddH2O. PCR cycling conditions were set at 95°C for 5 min, followed by 31 cycles of denaturation at 95°C for 30 s, annealing at 53°C for 15 s (16s rRNA) and 58°C for 15 s (ITS RNA), and extension at 72°C for 30 s. The melt curve was set for temperature from 60 to 95°C, and the curve was checked for every increase of 0.3°C.
Sequence Data Processing
Raw sequencing data were processed using an in-house Galaxy Pipeline1 (Feng et al., 2017). We detected barcodes used Trimmomatic-0.33 to discard the low-quality sequences and obtain high-quality sequences. We use Flash to combine the forward sequences and reverse sequences. Reads were assigned to the same OTU by using the UPARSE method according to the 97% identity threshold. All the analysis is based on the generated OTU table. 16S and ITS OTU representative sequences were taxonomically assigned by using RDP Classifier against the 16S ref. Silva database and ITS Refs database (Zhou et al., 2011). On average, we obtained high-quality sequences with 48,783 16S rRNA sequences and 42,487 ITS sequences per sample across all soil samples. We then resampled 33,371 16S rRNA sequences and 27,731 ITS sequences per sample for further microbial analyses.
Data Analysis
Based on the resampled OTU table, we calculated the bacterial and fungal α-diversity (Shannon index and Inv-Simpson Index). We measured the influences of N and rainfall patterns on the α-diversity index of bacteria and fungi by using one-way ANOVA analysis. The relative abundance of species was estimated through the data analysis platform (see text footnote 1) (Louca et al., 2016). Pearson correlation analysis was used to evaluate the relationship between soil physicochemical properties, plant root exudation, and bacteria α-diversity by using ggcorplot and ggthemes packages in RStudio (3.6.2). LefSe analysis was used to test the bacterial and fungal biomarker groups for N addition and altered rainfall patterns. LDA (linear discriminant analysis) measured the relative abundance of bacterial and fungal species (Yu et al., 2021a). We used the FAPROTAX database and FUNGuild pipeline to predict the bacterial and fungal functions by data process platform (see text footnote 1) (Yu et al., 2021b). Mantel and partial Mantel tests were used to examine the relationships between the Glomeromycota community composition and environmental variables.
Results
Plant Biomass and Soil Chemo-Physical Properties
The results of aboveground biomass, belowground biomass, and root C and N exudation rates and C:N ratios of root exudates, soil moisture, soil nitrate, and ammonium were previously published in Li C. et al. (2021). The SMBC significantly increased under nitrogen (N) addition and N addition and rainfall addition interaction (RA + N), and significantly decreased under rainfall reduction (RR) (Supplementary Table 1). The SMBN significantly increased under N addition, N addition and rainfall reduction interaction (RR + N), and RA + N (Supplementary Table 1). The SMBC:SMBN ratio significantly decreased in response to RR (Supplementary Table 1).
Factors Influencing Bacterial and Fungal Diversity
By themselves, N and rainfall did not significantly affect bacterial and fungal Shannon diversity and Inv-Simpson indices (Figure 1). Bacterial Shannon diversity and Inv-Simpson diversity indices significantly increased under the condition of RA + N (Figures 1A,B), while the fungal Shannon diversity index significantly decreased under rainfall reduction and N addition (RR + N) (Figures 1C,D). Pearson’s correlation analysis indicated that SMBN, inorganic nitrogen (NO3– and NH4+), S. krylovii root exudate C:N ratio, and A. polyrhizum root exudate C, N, and C:N ratio all significantly affected bacterial Shannon diversity index (Figure 2A). SMBC and A. polyrhizum root exudate C:N ratio significantly affected the fungal Shannon diversity index (Figure 2B).
Figure 1. The effects of altered nitrogen addition and rainfall, and their interaction, on bacterial Shannon index (A) and Inv-Simpson index (B), fungal Shannon index (C) and Inv-Simpson index (D) in August 2018. Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition. Values are mean ± SE (n = 5). Different lowercase letters indicate significant differences (p < 0.05) according to Duncan’s post hoc test.
Figure 2. Pearson correlations among bacterial Shannon index (A), fungi Shannon index (B), aboveground biomass, belowground biomass, soil moisture, soil microbial biomass carbon (SMBC), soil microbial biomass nitrogen (SMBN), SMBC:SMBN, NO3–, NH4+, root C exudates of Stipa krylovii, root N exudates of S. krylovii, root C: N exudates of S. krylovii, root C exudates of Allium polyrhizum, root N exudates of A. polyrhizum, and root C: N exudates of A. polyrhizum. The number represents the correlation coefficient. Coefficients that are not correlated have been hidden. Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition.
Bacterial and Fungi Community Composition and Biomarkers
At the phylum level, the bacterial community was dominated by Actinobacteria (mean⋅20%), Acidobacteria (mean 16%), and Proteobacteria (mean 13%) (Figure 3A). The fungal community was mainly constituted of Ascomycota (mean 78%), Basidiomycota (mean 10%), and Glomeromycota (mean 0.70%) (Figure 3B). LefSe showed that the bacterial biomarkers of RR mainly belonged to the phylum Proteobacteria (genus Ensifer, family Rhizobiaceae). The biomarkers of RA included the phyla Gemmatimonadetes (class Longimicrobia), Firmicutes (class Clostridia and order Bacillales), Bacteroidetes (order Cytophagales), and class Deltaproteobacteria. The biomarkers of N addition mainly belonged to the phylum Proteobacteria (order Neisseriales). The biomarker bacteria of RR + N mainly belonged to the phylum Proteobacteria, while the biomarkers of RA + N belonged to the phylum Proteobacteria (genus Inquilinus, family Rhodospirillaceae, and order Rhodospirillales) (LDA = 2, Figure 4A). The biomarker fungi of RR mainly belonged to the phyla Ascomycota (class Dothideomycetes) and Basidiomycota (genus Clavulinopsis, family Clavariaceae, and order Sebacinales). The fungal biomarkers of N treatment mainly belonged to the phylum Basidiomycota (genus Bovista and family Lycoperdaceae). The indicator fungi of RR + N treatment mainly belonged to the phylum Ascomycota (genus Nigrospora, family Trichosphaeriaceae, and order Trichosphaeriales). The indicator fungi of RA + N treatment belonged to the phylum Ascomycota (class Pezizomycetes, order Botryosphaeriales, genus Aureobasidium, family Aureobasidiaceae, and order Dothideales) (LDA = 4, Figure 4B).
Figure 3. The effects of altered nitrogen addition and rainfall, and their interaction, on the relative abundance of dominant members of the bacterial (A) and fungal (B) communities in August 2018. Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition.
Figure 4. Cladogram indicates the phylogenetic distribution of bacterial lineages (P = 0.05, LDA = 2.0) (A) and fungal lineages (P = 0.05, LDA = 4.0) (B) under the treatment of altered nitrogen addition and rainfall patterns. Each circle’s diameter is proportional to the given taxon’s relative abundance (n = 5). Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition.
Functional Prediction of Soil Bacteria and Fungi
There was no significant response by soil bacterial function prediction to N addition and rainfall patterns (Figure 5A). For fungal function prediction, only arbuscular mycorrhizal significantly decreased under the RR + N treatment (Figure 5B). As the arbuscular mycorrhizal (AM) fungi are dominated by the members of the phylum Glomeromycota, within belowground plant–microbe systems (Johnson et al., 2013), by calculating the absolute abundance of this phylum, it was shown that RA treatment significantly increased the absolute abundance of Glomeromycota as compared to RR + N treatment (Figure 6). Mantel test showed that the root C exudation rates of A. polyrhizum significantly affected the abundance of Glomeromycota (Table 1).
Figure 5. The effects of altered nitrogen addition and rainfall, and their interaction, on the relative abundance of bacterial (A) and fungal (B) functional groups based on the FAPROTAX tool and FUNGuild pipeline, respectively. Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition. Values are mean ± SE (n = 5). Different lowercase letters indicate the significant differences (p < 0.05) according to Duncan’s post hoc test. Different lowercase letters indicate significant differences (p < 0.05) according to Duncan’s post hoc test.
Figure 6. The effects of nitrogen addition and rainfall and their interaction on the absolute abundance of Glomeromycota. Different treatments: Con: Control; RR: Rainfall reduction; RA: Rainfall addition; N: Nitrogen addition; RR + N: Rainfall reduction and Nitrogen addition; RA + N: Rainfall addition and Nitrogen addition. Values are mean ± SE (n = 5). Different lowercase letters indicate significant differences (p < 0.05) according to Duncan’s post hoc test.
Discussion
Effects of Nitrogen Addition and Rainfall Interaction on Soil Bacterial and Fungal Diversity
Our results showed that nitrogen (N) addition did not affect the diversity of bacteria and fungi, which was consistent with the non-linear response of bacterial diversity to N addition found by Liu et al. (2021a), and supports our first hypothesis. A previous study showed that the diversity of bacteria decreased with N addition when the level of N input was between 16 and 32 g N m–2 yr–1 (Liu et al., 2021a), while Zhang et al. (2018) found that fungal diversity did not significantly change under N addition alone (< 100 kg N⋅ha–1⋅yr–1). In our study, N addition was 10 g N m–2 yr–1, which was below the threshold level required to affect bacteria and fungi, and therefore the α-diversities of bacteria and fungi were not significantly increased by N addition. However, the treatment of RA + N significantly increased the bacterial Shannon diversity index. The interaction between N and rainfall may provide better conditions for bacterial survival, which may be the main reason for the increased bacterial α-diversity under RA + N (Shade et al., 2012). When water is sufficient, N addition has a positive effect on bacterial diversity (Zhang et al., 2015). Intriguingly, our results revealed that N addition alone did not affect fungal diversity and that RR + N decreased fungal diversity, therefore ignoring the factor that rainfall may influence judgment of the effect of N addition on soil fungal diversity. Consistent with our research, Wang H. et al. (2020). showed that rainfall played a key role in influencing the composition of the fungal community under N addition. Therefore, rainfall can alleviate the effects of N deposition on fungal diversity. Previous research found that N addition increased fungal α-diversity, and that the interaction between N and irrigation did not affect fungal diversity (Zhang et al., 2018). In total, soil fungal diversity might be overestimated if variations in moisture under N deposition are ignored (Zhang et al., 2018). Also, previous research, using PLFA technology, provided evidence that N deposition decreased the relative abundance and biomass of bacteria and fungi (Ma et al., 2016; Tian et al., 2017), and that under a high level of N addition (200 kg N ha–1 yr–1) in the absence of water leaching, the growth of fungi was significantly inhibited (Sun et al., 2015). Therefore, under the interaction of RR + N, the α-diversity of fungi was significantly reduced. Moreover, some studies have indicated that the fungal community was more sensitive than bacteria under N deposition (Herold et al., 2012; Wang J. et al., 2020; Liu et al., 2021b). The long-term response of bacterial and fungal diversity to N addition and rainfall interaction requires further experimental verification.
In grassland ecosystems, soil microbial diversity was influenced by biological and abiotic factors, including soil nutrients, aboveground biomass, and litter quality (Contosta et al., 2015). Using meta-analysis, Wang et al. (2018) showed that decreased soil microbial diversity was associated with the decreased microbial biomass under N deposition. Our results showed that there was a significant correlation between the bacterial diversity and SMBN, ammonia nitrogen, and nitrate nitrogen. A possible explanation for the correlation may be that N addition increases ammonia nitrogen and nitrate nitrogen, stimulating the soil bacteria to uptake inorganic nitrogen and increase bacterial biomass and diversity. Due to the carbon–nitrogen coupling mechanism of bacteria, N input accelerates bacterial carbon uptake and increases bacterial activity, which stimulates the ability of bacteria to degrade litter (Schleuss et al., 2019). This may be a possible explanation for the increase in bacterial diversity under RR + N and RA + N conditions. Our results indicated that the main factors affecting fungal diversity were SMBC and the A. polyrhizum root exudate C:N ratio. In general, fungal communities are less affected by the external environment in the short term when they are mycelium rich. However, our study results indicated that soil abiotic and biological factors had significant effects on bacterial and fungal diversity under N, rainfall, and their interaction. It also highlights the impact of climate change on plant–soil–microbe interactions.
Effects of Nitrogen Addition and Rainfall Interaction on Biomarker Organisms and the Survival Strategies of Soil Bacteria and Fungi
Actinobacteria and Acidobacteria dominated all of the bacterial communities in our study. Based on Pianka’s r/K-selection theory, both Actinobacteria (20%) and Proteobacteria (13%) grow faster and are eutrophic r-strategists (Fierer et al., 2007; Schleuss et al., 2019; Sun et al., 2021), while Acidobacteria (16%) are oligotrophic K-strategists (Sun et al., 2021). The majority of soil bacteria are r-strategists, with a few being K-strategists (Sun et al., 2021). Through Lefse analysis across treatments, we found nitrogen (N) addition mainly increased the relative abundance of the Proteobacteria, which are r-strategists, and therefore, N addition increased the bacterial r-strategist community. A potential explanation for this might be that N addition increased inorganic nitrogen content and stimulated microbial input of nitrogen and labile C (Fontaine et al., 2003; Fierer et al., 2012), resulting in faster microbial growth rates. This is in agreement with Sun et al. (2015) and is consistent with our second hypothesis that copiotrophic groups of bacteria will be mainly observed in a nutrient-rich environment under N addition treatment. Rainfall addition mainly increased the relative abundance of the phyla Gemmatimonadetes, Firmicutes, Bacteroidetes. Gemmatimonadetes, Firmicutes, and Bacteroidetes, all of which are copiotrophic r-strategists (Li H. et al., 2021; Sun et al., 2021). The effect of N and rainfall interaction also increased the relative abundance of the eutrophic bacterial community, which is in agreement with previous research that N and rainfall interaction promoted plant growth, plant litter quantity, dissolved carbon and nitrogen, and root exudate components, therefore accelerated bacterial growth, especially that of copiotrophic r-strategists (Meier and Bowman, 2008; Millard and Singh, 2010). As for fungal community structure, N addition mainly increased the relative abundance of the phylum Basidiomycota, while the effects of N and rainfall interaction increased the relative abundance of Ascomycota and Basidiomycota. These two phyla mainly represent copiotrophic r-strategists and oligotrophic k-strategists, respectively (Yao et al., 2017). The fungal communities in the control soil were also dominated by Ascomycota and Basidiomycota, and N addition and rainfall had little effect on fungi over the short term of our study. In fact, due to the rich mycelium of fungal communities, they are not very sensitive to external environmental stimuli (Berlemont, 2017; Schleuss et al., 2019).
Effects of Nitrogen and Rainfall Interaction on the Functions of Soil Bacteria and Fungi
In this study, the interaction between N addition and rainfall had no significant effect on the bacteria participating in soil carbon and nitrogen cycling. This was inconsistent with our third hypothesis. Our results suggested that although the bacterial community saw an increased copiotrophic r-strategy groups, there was no obvious effect on the degradation of soil organic matter. The reason for this inconsistency may be that the short timescale (only 1 year) of the treatment did not significantly affect soil physical and chemical properties in our study, and bacteria had functional redundancy that they did not affect the soil carbon and nitrogen cycle. Moreover, soil bacteria can flexibly adapt to environmental changes, leading to bacterial functions not being significantly affected (Shade et al., 2012). Previous studies found that long-term fertilization significantly decreased microbial functional guilds involved in carbon fixation, degradation, N fixation, and mineralization. Because changes in soil functional community structure lead to a decrease in stochastic processes in the community assembly, microbial community structure tends to develop in a deterministic manner (Liang et al., 2020). Based on the functions of fungi predicted by using the FUNGuild database, we found that the AMF fungi community functional classification was significantly reduced under RR + N. AMF play an important role in ecosystem functions and regulate ecosystem stability by providing plants with N and P nutrients through symbiotic association with plant roots (Cotton et al., 2015; Yang et al., 2018). Previous studies found that N enrichment negatively affected AMF spore abundance, species richness, and biodiversity (Egerton-Warburton and Allen, 2000; Liu et al., 2013; Chen et al., 2017). A possible explanation for the significant decrease of AMF under RR + N treatment observed in the current study may be that N addition increased the inorganic nitrogen content (ammonia and nitrate nitrogen) of the soil. Plants prefer to utilize inorganic nitrogen in the soil for their own growth, and fertilization would reduce the demand for inorganic nitrogen produced by microbial decomposition and mineralization. Moreover, the AMF are a significant fungal community that participates in soil nitrogen mineralization (Cavagnaro et al., 2012; Bender et al., 2015). The reduction of soil water content also reduces the leaching of soil nutrient elements. Plants would require fewer nutrient elements from AMF community decomposition. Correspondently, the relative abundance of AMF was significantly reduced under RR + N conditions in our study.
The AMF of belowground ecosystem communities is mainly dominated by the phylum Glomeromycota (Johnson et al., 2013). We calculated the absolute abundance of the Glomeromycota through the absolute quantification of the ITS gene, and our results found that the absolute abundance of Glomeromycota did not significantly change between control and RR + N. However, it still exhibited a decreasing tendency, which was somewhat consistent with our fungal function prediction results. Our results further indicated that rainfall addition had positive effects on the phylum Glomeromycota. Rainfall addition increased the leaching of soil nutrients, especially NO3–. This would greatly increase plant demand for nitrogen, phosphorus, and other nutrients, and may cause AMF to symbiotically interact with plant roots to provide more of these nutrients to the host plants. Therefore, the absolute abundance of AMF significantly increased under rainfall addition. Mantel test showed that the rates of root exudate C from A. polyrhizum significantly affected the relative abundance of Glomeromycota. Because the Glomeromycota may have a symbiotic relationship with A. polyrhizum roots, the relative abundance of these fungi may be affected to a greater degree by plant root exudates. In addition, regarding the functional microbial groups involved in the soil carbon and nitrogen cycles, the response of AMF under the interaction of N and rainfall is more sensitive than that of other groups, suggesting that the interaction between these two factors has a significant impact on AMF over short timescales.
Conclusion
In this study, we studied how the diversity and function of soil bacteria and fungi responded to the interaction between N addition and rainfall patterns in a temperate steppe of Inner Mongolia. The interaction between N addition and altered rainfall patterns had obvious effects on the diversity of these microbial communities. Additionally, the C:N ratios of root exudates from A. polyrhizum affected not only the diversity of soil bacteria, but also the diversity of fungi. The functions of bacteria did not significantly change, while the AMF in fungal function classification was obviously reduced under RR + N, suggesting that plants will first use the inorganic nitrogen in the soil for their growth, reducing the demand for inorganic nitrogen produced by AMF decomposition and mineralization. The absolute abundance of the phylum Glomeromycota significantly increased under the treatment of rainfall addition, potentially illustrating that AMF were more sensitive to changes in nitrogen and rainfall patterns. Our findings suggest that this may provide novel and important insights for understanding and predicting microbial-mediated ecological processes under global climate change. Lastly, our experiment examined the short-term response of bacterial and fungal diversity and functions to changes in soil nitrogen and moisture. However, there is still a need for further research about the long-term response mechanism of bacterial and fungal functions and plant–soil–microbe interactions.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The name of the repository and accession number(s) can be found below: National Center for Biotechnology Information (NCBI) BioProject, https://www.ncbi.nlm.nih.gov/bioproject/, PRJNA826987 and PRJNA826724.
Author Contributions
CX designed this study. YY and LL conducted the study. YY, CX, JZ, SW, and YZ drafted the manuscript. All authors contributed to revisions and have read and approved the manuscript.
Funding
This study was financially supported by the National Natural Science Foundation of China (31770501) and Key Laboratory of Ecology and Environment in Minority Areas (Minzu University of China), and National Ethnic Affairs Commission (KLEEMA202103).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to Lang Zheng for their help in samplings during the field experiment and to Hulun Lake Reserve Grassland Ecology Research Station of Minzu University of China for logistics and permission to access the study sites.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.906818/full#supplementary-material
Footnotes
References
Allen, M., Babiker, M., Chen, Y., and Zickfeld, K. (2018). “Summary for Policymakers. In: Global warming of 1.5°C,” in An IPCC Special Report on the impacts of global warming 1.5°C, Vol. 1, (IPCC).
Bardgett, R. D., and Putten, W. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi: 10.1038/nature13855
Baudoin, E., Benizri, E., and Guckert, A. (2003). Impact of artificial root exudates on the bacterial community structure in bulk soil and maize rhizosphere. Soil Biol. Biochem. 35, 1183–1192. doi: 10.1016/S0038-0717(03)00179-2
Bender, S. F., Conen, F., and Van, D. (2015). Mycorrhizal effects on nutrient cycling, nutrient leaching and N2O production in experimental grassland. Soil Biol. Biochem. 80, 283–292. doi: 10.1016/j.soilbio.2014.10.016
Berlemont, R. (2017). Distribution and diversity of enzymes for polysaccharide degradation in fungi. Sci. Rep-UK. 7:222. doi: 10.1038/s41598-017-00258-w
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
Cavagnaro, T. R., Barrios-Masias, F. H., and Jackson, L. E. (2012). Arbuscular mycorrhizas and their role in plant growth, nitrogen interception and soil gas efflux in an organic production system. Plant Soil 353, 181–194.
Chen, Y. L., Xu, Z. W., Xu, T. L., Veresoglou, S. D., Yang, G. W., and Chen, B. D. (2017). Nitrogen deposition and precipitation induced phylogenetic clustering of arbuscular mycorrhizal fungal communities. Soil Biol. Biochem. 115, 233–242. doi: 10.1016/j.soilbio.2017.08.024
Contosta, A. R., Frey, S. D., and Cooper, A. B. (2015). Soil microbial communities vary as much over time as with chronic warming and nitrogen additions. Soil Biol. Biochem. 88, 19–24. doi: 10.1016/j.soilbio.2015.04.013
Cotton, T. E. A., Fitter, A. H., Miller, R. M., Dumbrell, A. J., and Helgason, T. (2015). Fungi in the future: interannual variation and effects of atmospheric change on arbuscular mycorrhizal fungal communities. New Phytol. 205, 1598–1607. doi: 10.1111/nph.13224
Egerton-Warburton, L. M., and Allen, E. B. (2000). Shifts in arbuscular mycorrhizal communities along an anthropogenic nitrogen deposition gradient. Ecol. Appl. 10, 484–496. doi: 10.1890/1051-07612000010[0484:siamca]2.0.co;2
Feng, K., Zhang, Z., Cai, W., Liu, W., Xu, M., Yin, H., et al. (2017). Biodiversity and species competition regulate the resilience of microbial biofilm community. Mol. Ecol. 26, 6710–6082. doi: 10.1111/mec.14356
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological calssification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. P. Natl. A Sci. 103, 626–631. doi: 10.1073/pnas.0507535103
Fierer, N., Leff, J. W., Adams, B. J., Nielsen, U. N., Bates, S. T., Lauber, C. L., et al. (2012). Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. P. Natl. A. Sci. 109:21390. doi: 10.1073/pnas.1215210110
Fontaine, S., Mariotti, A., and Abbadie, L. (2003). The priming effect of organic matter: a question of microbial competition? Soil Biol. Biochem. 35, 837–843. doi: 10.1016/S0038-0717(03)00123-8
Fukami, T. (2015). Historical contingency in community assembly: integrating niches, species pools, and priority effects. Annu. Rev. Ecol. Evol. S. 46, 1–23. doi: 10.1146/annurev-ecolsys-110411-160340
Füzy, A., Biró, B., Tóth, T., Hildebrandt, U., and Bothe, H. (2008). Drought, but not salinity, determines the apparent effectiveness of halophytes colonized by arbuscular mycorrhizal fungi. J. Plant Physiol. 165, 1181–1192. doi: 10.1016/j.jplph.2007.08.010
Galloway, J., Dentener, F., Capone, D., Boyer, E., Howarth, R., Seitzinger, S., et al. (2004). Nitrogen cycles: past, present, and future. Biogeochemistry 70, 153–226.
Gardes, M., and Bruns, T. D. (2010). ITS primers with enhanced specificity for basidiomycetes–application to the identification of mycorrhizae and rusts. Mol. Ecol. 2, 113–118. doi: 10.1111/j.1365-294X.1993.tb00005.x
Gordon, H., Haygarth, P. M., and Bardgett, R. D. (2008). Drying and rewetting effects on soil microbial community composition and nutrient leaching. Soil Biol. Biochem. 40, 302–311. doi: 10.1016/j.soilbio.2007.08.008
Herold, M. B., Baggs, E. M., and Daniell, T. J. (2012). Fungal and bacterial denitrification are differently affected by long-term pH amendment and cultivation of arable soil. Soil Biol. Biochem. 54, 25–35. doi: 10.1016/j.soilbio.2012.04.031
Ihrmark, K., Bödeker, I. T. M., Cruz-Martinez, K., Friberg, H., Kubartova, A., Schenck, J., et al. (2012). New primers to amplify the fungal ITS2 region – evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 82, 666–677. doi: 10.1111/j.1574-6941.2012.01437.x
Jie, B., Zhang, N., Yu, L., Yang, H., and Ma, K. (2011). Interactive effects of water and nitrogen addition on soil microbial communities in a semiarid steppe. J. Plant Ecol. 5:3. doi: 10.1093/jpe/rtr046
Johnson, N. C., Angelard, C., Sanders, I. R., and Kiers, E. T. (2013). Predicting community and ecosystem outcomes of mycorrhizal responses to global change. Ecol. Lett. 16, 140–153. doi: 10.1111/ele.12085
Li, C., Liu, L., Zheng, L., Yu, Y., Mushinski, R. M., Zhou, Y., et al. (2021). Greater soil water and nitrogen availability increase C: N ratios of root exudates in a temperate steppe. Soil Biol. Biochem. 161:108384. doi: 10.1016/j.soilbio.2021.108384
Li, H., Yang, S., Semenov, M. V., Yao, F., Ye, J., Bu, R., et al. (2021). Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community. Glob. Change Biol. 27, 2763–2779. doi: 10.1111/gcb.15593
Liang, Y., Ning, D., Lu, Z., Zhang, N., Hale, L., Wu, L., et al. (2020). Century long fertilization reduces stochasticity controlling grassland microbial community succession. Soil Biol. Biochem. 151:108023. doi: 10.1016/j.soilbio.2020.108023
Liu, J., Wang, B., Cane, M. A., Yim, S. Y., and Lee, J. Y. (2013). Divergent global precipitation changes induced by natural versus anthropogenic forcing. Nature 493, 656–659. doi: 10.1038/nature11784
Liu, W., Jiang, L., Yang, S., Wang, Z., Tian, R., Peng, Z., et al. (2021a). Critical transition of soil bacterial diversity and composition triggered by nitrogen enrichment. Ecology 101:e03053. doi: 10.1002/ecy.3053
Liu, W., Liu, L., Yang, X., Deng, M., Wang, Z., Wang, P., et al. (2021b). Long-term nitrogen input alters plant and soil bacterial, but not fungal beta diversity in a semiarid grassland. Glob. Change Biol. 27, 3939–3950. doi: 10.1111/gcb.15681
Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353:1272. doi: 10.1126/science.aaf4507
Ma, H., Bai, G., Sun, Y., Kostenko, O., Zhu, X., Lin, S., et al. (2016). Opposing effects of nitrogen and water addition on soil bacterial and fungal communities in the Inner Mongolia steppe: a field experiment. Appl. Soil Ecol. 108, 128–135. doi: 10.1016/j.apsoil.2013.03.006
Marvel, K., and Bonfils, C. (2013). Identifying external influences on global precipitation. P. Natl. A. Sci. 110, 19301–19306. doi: 10.1073/pnas.1314382110
Meier, C. L., and Bowman, W. (2008). Links between plant litter chemistry, species diversity, and below-ground ecosystem function. P. Natl. A. Sci. 105, 19780–19785. doi: 10.1073/pnas.0805600105
Millard, P., and Singh, B. K. (2010). Does grassland vegetation drive soil microbial diversity? Nutr. Cycl. Agroecosys. 88, 147–158. doi: 10.1007/s10705-009-9314-3
Phillips, R. P., Erlitz, Y., Bier, R., and Bernhardt, E. S. (2008). New approach for capturing soluble root exudates in forest soils. Funct. Ecol. 22, 990–999. doi: 10.1111/j.1365-2435.2008.01495.x
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
Schleuss, P. M., Widdig, M., Heintz-Buschart, A., Guhr, A., Martin, S., Kirkman, K., et al. (2019). Stoichiometric controls of soil carbon and nitrogen cycling after long-term nitrogen and phosphorus addition in a mesic grassland in South Africa. Soil Biol. Biochem. 135, 294–303. doi: 10.1016/j.soilbio.2019.05.018
Shade, A., Peter, H., Allison, S. D., Baho, D. L., Berga, M., Bürgmann, H., et al. (2012). Fundamentals of Microbial Community Resistance and Resilience. Front. Microbiol. 3:417. doi: 10.3389/fmicb.2012.00417
Sun, L., Qi, Y., Dong, Y., He, Y., Peng, Q., Liu, X., et al. (2015). Interactions of water and nitrogen addition on soil microbial community composition and functional diversity depending on the inter-annual precipitation in a Chinese steppe. J. Integr. Agr. 14, 788–799. doi: 10.1016/S2095-3119(14)60773-5
Sun, Y., Wang, C., Yang, J., Liao, J., Chen, H., and Ruan, H. (2021). Elevated CO2 shifts soil microbial communities from K- to r-strategists. Glob. Ecol. Biogeogr. 30, 961–972. doi: 10.1111/geb.13281
Tian, B., Zhu, M., Pei, Y., Ran, G., Shi, Y., and Ding, J. (2021). Climate warming alters the soil microbial association network and role of keystone taxa in determining wheat quality in the field. Agr. Ecosyst. Environ. 326:107817. doi: 10.1016/j.agee.2021.107817
Tian, D., Jiang, L., Ma, S., Fang, W., Schmid, B., Xu, L., et al. (2017). Effects of nitrogen deposition on soil microbial communities in temperate and subtropical forests in China. Sci. Total Environ. 607- 608, 1367–1375. doi: 10.1016/j.scitotenv.2017.06.057
Turner, M. M., and Henry, H. A. L. (2010). Interactive effects of warming and increased nitrogen deposition on 15N tracer retention in a temperate old field: seasonal trends. Glob. Change Biol. 15, 2885–2893. doi: 10.1111/j.1365-2486.2009.01881.x
Vance, E. D., Brookes, P. C., and Jenkinson, D. S. (1987). An extraction method for measuring soil microbial biomass C. Soil Biol. Biochem. 19, 703–707. doi: 10.1016/0038-0717(87)90052-6
Vogt, K. A., and Persson, H. (1991). Measuring growth and development of roots. Techniques and approaches in forest tree ecophysiology. Boca Raton: CRC Press, 477–501.
Waldrop, M. P., Zak, D. R., Sinsabaugh, R. L., Gallo, M., and Lauber, C. (2004). Nitrogen deposition modifies soil carbon storage through changes in microbial enzymatic activity. Ecol. Appl. 14, 1172–1177. doi: 10.1890/03-5120
Wang, C., Liu, D., and Bai, E. (2018). Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biol. Biochem. 120, 126–133. doi: 10.1016/j.soilbio.2018.02.003
Wang, H., Na, T., Jin, K., Ji, B., Schellenberg, M., Wei, Z., et al. (2020). Interactive effects of nitrogen fertilizer and altered precipitation on fungal communities in arid grasslands of northern china. J. Soil Sediment. 20, 1–13.
Wang, J., Shi, X., Zheng, C., Suter, H., and Huang, Z. (2020). Different responses of soil bacterial and fungal communities to nitrogen deposition in a subtropical forest. Sci. Total Environ. 755:142449. doi: 10.1016/j.scitotenv.2020.142449
Wang, M., Shi, S., Lin, F., Hao, Z., Jiang, P., and Dai, G. (2012). Effects of soil water and nitrogen on growth and photosynthetic response of Manchurian ash (Fraxinus mandshurica) seedlings in Northeastern China. PLoS One. 7:e30754. doi: 10.1371/journal.pone.0030754
Williams, M. A., and Rice, C. W. (2007). Seven years of enhanced water availability influences the physiological, structural, and functional attributes of a soil microbial community. Appl. Soil Ecol. 35, 535–545. doi: 10.1016/j.apsoil.2006.09.014
Yang, G., Wagg, C., Veresoglou, S. D., Hempel, S., and Rillig, M. C. (2018). How soil biota drive ecosystem stability. Trends Plant Sci. 23, 1057–1067. doi: 10.1016/j.tplants.2018.09.007
Yao, F., Shan, Y., Wang, Z., Wang, X., Ye, J., Wang, X., et al. (2017). Microbial Taxa Distribution Is Associated with Ecological Trophic Cascades along an Elevation Gradient. Front. Microbiol. 8:2071. doi: 10.3389/fmicb.2017.02071
Yu, Y., Zang, L., Zhou, Y., Sang, W., Zhao, J., Liu, L., et al. (2021a). Changes in soil microbial community structure and function following degradation in a temperate grassland. J. Plant Ecol. 14, 384–397. doi: 10.1093/jpe/rtaa102
Yu, Y., Liu, L., Wang, J., Zhang, Y., and Xiao, C. (2021b). Effects of warming on the bacterial community and its function in a temperate steppe. Sci. Total Environ. 792:148409. doi: 10.1016/j.scitotenv.2021.148409
Zhang, H., Wang, L., Liu, H., Zhao, J., Li, G., Wang, H., et al. (2018). Nitrogen deposition combined with elevated precipitation is conducive to maintaining the stability of the soil fungal diversity on the Stipa baicalensis steppe. Soil Biol. Biochem. 117, 135–138. doi: 10.1016/j.soilbio.2017.11.004
Zhang, N., Wan, S., Guo, J., Han, G., Gutknecht, J., Schmid, B., et al. (2015). Precipitation modifies the effects of warming and nitrogen addition on soil microbial communities in northern Chinese grasslands. Soil Biol. Biochem. 89, 12–23. doi: 10.1016/j.soilbio.2015.06.022
Keywords: nitrogen addition, rainfall addition, rainfall reduction, bacterial diversity, fungal diversity, functional prediction, temperate steppe
Citation: Yu Y, Liu L, Zhao J, Wang S, Zhou Y and Xiao C (2022) The Diversity and Function of Soil Bacteria and Fungi Under Altered Nitrogen and Rainfall Patterns in a Temperate Steppe. Front. Microbiol. 13:906818. doi: 10.3389/fmicb.2022.906818
Received: 29 March 2022; Accepted: 29 April 2022;
Published: 28 June 2022.
Edited by:
José A. Siles, Center for Edaphology and Applied Biology of Segura (CSIC), SpainReviewed by:
Longyu Hou, Institute of Botany (CAS), ChinaLinna Ma, Institute of Botany (CAS), China
Copyright © 2022 Yu, Liu, Zhao, Wang, Zhou and Xiao. 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: Chunwang Xiao, cwxiao@muc.edu.cn; Yijun Zhou, zhouyijun@muc.edu.cn