With global nitrogen (N) deposition continuously increasing, few reports exist describing how soil bacteria respond at the molecular level to long-term addition of excess N in variously aged forests. To reveal the responses of soil bacteria to the elevated N, an experiment was conducted with a chronic nutrient addition within differently aged stands (46-, 78-, and 200-years-old) in the northern subtropical China since 2011, including three treatments, CK (no N nor phosphorus (P) additions), N treatment (N, 100 kg N ha-1), and N with P (N+P, 100 kg N ha-1 + 50 kg P ha-1) to examine potential P limitation under N deposition. Metagenomic sequencing was used to examine the snapshot responses of soil bacterial communities. Soil moisture and texture, ammonium, nitrate, SOC (soil organic carbon), TN (soil total nitrogen), TP (total phosphorus), DOC (dissolved organic carbon), DON (dissolved organic nitrogen) were measured to explain the influence mechanism of forest age and fertilization on changes of microbial community. Following N addition, soil bacterial community diversity and most dominant phyla increased, but they showed a decrease with increasing stand age. The effects of fertilization on the same taxa were variable across forest ages. Soil bacterial community responded differently in 7-year fertilization, with distinct shift in 46-year-old forest and adaptability to long-term N addition in the 200-year-old forest. Soil texture and moisture, DOC, DON, pH, SOC/TN and TP were significantly correlated with bacterial community across stand ages, while N fertilization affected the bacterial community mostly via inducing soil moisture, NO3--N, DOC and pH in the 46-year-old forest, whose effects decreased with increasing stand age. Our results suggest that due to the variation of soil physicochemical properties among forest ages, soil bacterial communities are more stable and resilient to N deposition with increasing the age of stands. Soil bacterial communities might not encounter P limitation following the long-term addition of N in the subtropics.
Nitrogen (N) is one of the key elements involved in ecosystem functioning. Currently, due to the use and emission of active N, China has become the third largest N deposition area after Europe and North America. Total nitrogen oxide (NOx) emissions increased from 8.4 Tg yr-1 in 1990 to 11.3 Tg yr-1 in 2000, whereas total NH3 emissions increased from 10.8 to 13.6 Tg during the same period (Lue & Tian 2007, Liu et al. 2013). Long-term inputs of inorganic N via atmospheric deposition may have a series of negative impacts on forest ecosystems, such as soil acidification (Tian et al. 2017, Li et al. 2019, Wu et al. 2019), conspicuous changes in soil nutrient (Li et al. 2019), and loss of plant community diversity (Zhang et al. 2017).
Forest soils maintain the most diverse microbial habitats on earth, in which bacteria are the most abundant group of microorganisms (Hardoim et al. 2015). These bacteria contribute to a range of essential soil processes involved in the cycling of soil elements, and in mineral weathering that leads to the release of inorganic nutrients (Uroz et al. 2011). The composition and functional group definition of bacterial phyla, genera and species have been reported in previous studies (Weand et al. 2010, Marketa et al. 2011, Llado et al. 2017). The addition of N is considered to increase the availability of plant-derived C and/or N resources to meet the demands of microbial growth in the soils of most terrestrial ecosystems (Weand et al. 2010, Marketa et al. 2011). Additionally, changes in soil physicochemical properties resulted directly or indirectly from the addition of N could affect soil microbial communities (Nie et al. 2018, Li et al. 2019).
Due to bacteria exhibiting an extremely broad variety of metabolic types, the mechanisms underlying the microbial feedback of soils in response to global climate change remain elusive (Frey et al. 2014, Llado et al. 2017). Nitrogen enters ecosystems largely through fixation, which is dominated by bacteria, accounting for more than 95% of the N input in undisturbed forests. The diversity of N-fixing bacteria was found to change rapidly with exogenous N inputs (Reed et al. 2011). Ramirez et al. (2012) found that the abundance of Acidobacteria and Verrucomicrobia decreased, while Actinobacteria and Firmicutes increased, following the addition of N to a variety of soils. Apart from the modified composition of soil microbial communities, the biomass of soil bacteria reduced by up to 50% (Frey et al. 2014).
Several studies have revealed that bacterial diversity declines with the rate of N addition (Kaspari et al. 2017, Wang et al. 2018b), while a reduction in oligotrophic bacteria was often concomitantly found with decreased enzyme activity related to cellulose breakdown. This further supports that the N addition negatively affects microbial SOM (soil organic matter) decomposition. However, when the N addition rate was lower than 100 kg N ha-1 year-1 (about ten times higher than the normal global rate), a positive response on microbial growth was found (Zhou et al. 2017). Distinct influences on the community composition of soil bacteria suggest that the effects of N deposition on soil microbes are likely ecosystem-specific (Llado et al. 2017, Zhang et al. 2018b).
Significant changes in the composition of microbial communities have been reported during the development of forest stands (Jangid et al. 2013). Among Eucalyptus plantations, the bacterial community diversity differed significantly with different ages, other than different vegetation species (Qu et al. 2020). Such shifts in bacterial community structures are largely driven by changes in resource availability and alteration of soil chemistry in response to forest succession (Li et al. 2019, Zhao et al. 2019). Due to that the responses of bacteria to future high N deposition rates might be primarily contingent on bacterial composition, background acidification status or nutrient availability (Mao et al. 2017, Zhao et al. 2019), the responses of soil bacterial communities to N addition in differently aged forests may be significantly variable, and the initially younger sites with less-fertility may show greater compositional changes (Gravuer et al. 2020). However, few studies have been carried out to verify this hypothesis or specify which environmental factors exert the strongest effects on bacterial communities at different ages.
Although N and phosphorus (P) supplies are relatively closely balanced in most natural ecosystems (Elser et al. 2007), single element enrichment can lead to a limitation by the alternative nutrient. Soil microorganisms may also be limited by the availability of N or P (Vitousek et al. 2010, Wardle 2010). As the result of high N deposition, P limitation may occur in soil bacterial communities. The addition of P could partially mitigate the negative effects of N addition on Shannon’s indices (Wang et al. 2018b), as well as the biomass of soil bacteria (Cui et al. 2017).
Determining how relationships between site factors and bacterial communities affect the equilibrium between soil organic matter decomposition and C sequestration is vital for the prediction of responses of forest ecosystems to N deposition (Singh et al. 2010). Our study aimed to reveal the effects of N addition and N+P interactions on the richness and community composition of soil bacteria, as well as the underlying mechanisms in differently aged subtropical forests. We initially hypothesized that: (i) the heterogeneity of microbial community was distinct under different forest ages, meanwhile the effect of N-induced was partly offset by the change of environment with forest age, and the negative response of bacterial community in the older forest with stronger acidity and higher nutrient content to N deposition would be weaker; (ii) the addition of P to N fertilization would cause a relatively reduced decline in soil bacterial abundance, and lead to a change in soil bacterial community composition, by alleviating nutrient limitation for bacterial communities. Furthermore, based on the major physicochemical factors affecting soil bacterial community, we explored whether the mediators of microbial community variation were the same in different stands following a continuous seven-year nutrient additions.
Materials and methodsStudy area and sites
The study area encompasses subtropical evergreen broadleaved forests of different ages in Dashan, Shitai County, Anhui Province in China (DS, 200-year-old - 30° 01′ 47″ N, 117° 21′ 23″ E), Zhawan Forest Farm, Qimen County (QM, 78-year-old - 29° 37′ 05″ N, 117° 32′ 30″ E), and Jitan, Shitai County (JT, 46-year-old - 30° 16′ 28″ N, 117° 26′ 23″ E) were selected (Fig. S1 in Supplementary material). The forest ages at the beginning of the experiment were determined based on the forest sub-compartment data of Forest Management Inventory (personal communication). The distance among sites is about 40 km, with an average elevation of 200-500 m a.s.l. This region belongs to the northern part of subtropical monsoon climate, and has distinct seasons (long winters and summers and a short spring and autumn). The annual average temperature is 16.0-16.6 °C, with the annual maximum of 38.8 °C and the minimum of -8.9 °C. The average annual precipitation ranges from 1521-1790 mm, of which 71% is primarily distributed during April-September. The experimental forests are secondary evergreen broad-leaved forest.
The dominant tree species were Castanopsis eyrei, C. sclerophylla, and Cyclobalanopsis glauca. The understorey vegetation was dominated by Camellia fraterna, Symplocos sumuntia, Lindera aggregate, and Woodwardia japonica. The detailed description of the experimental forests are given in Dai et al. (2018) and Xu et al. (2016). The soils of different plots are reddish yellow soils developed from the parent material of granite weathering. Moreover, to ensure the effects of the inherent soil properties (i.e., soil texture) and forest age, we conducted ribbon analysis using the procedure described by Thien (1979) and sampled soil and conducted the fractionation analysis as described by Hunt & Gilkes (1992). Both tests confirmed that soil texture of the all stands are sandy loam and silty loam (Tab. S1 in Supplementary material).
Experimental design
Based on the N deposition rate in the study area (bulk deposition 13 kg N ha-1 year-1), the treatments were designed as: CK (no N and P addition), N treatment (100 kg N ha-1 yr-1), and N+P treatment (100 kg N ha-1 yr-1 + 50 kg P ha-1 yr-1), with three replicate plots for each treatment. The experiment began in 2011 at the abovementioned three sites, with a completely randomized design.
The plot size was 20 × 20 m, with a 10-m wide buffer zone between adjacent plots. The N was applied in the form of NH4NO3, and P was in the form of Ca(H2PO4)2. The nutrients were dissolved in 20 L of stream water and evenly sprayed onto the corresponding plots near the soil surface using a backpack sprayer every two months from June 2011. The controls received the same amount of stream water (20 L) in each treatment event.
Sampling and measurements
In late July 2018, fresh soil samples from the 0-10 cm soil layers were extracted from each plot according to the five-point method, sieved (2 mm), thoroughly homogenized, and divided into two subsamples: one for soil analysis and another being stored at -80 °C for DNA extraction and subsequent molecular analysis. Soil bulk density (BD) was determined using a core sampler. The soil water content, ammonium (NH4+-N), nitrate (NO3--N), dissolved organic carbon (DOC), and dissolved organic nitrogen (DON) were measured with fresh soil.
Soil moisture was quantified by oven-drying at 105 °C to a constant weight. NH4+-N and NO3--N were measured by a FIAstar™ 5000 flow injection analyzer (FOSS Analytical, Hilleroed, Denmark). DOC and DON were extracted with a 0.5 M K2SO4 solution, and the extract was measured using a Multi 3100® C/N analyzer (Analytik Jena AG, Jena, Germany). The soil samples were air-dried for two weeks to a constant weight, and ground and passed though using a 0.154-mm sieve for measuring the soil organic carbon (SOC), total nitrogen (TN), pH, and total phosphorus (TP). SOC and TN were determined using an elemental analyzer (EA3000®, EuroVector, Pavia, Italy). After mixing water and soil to a ratio of 2.5:1, the pH of the soil solution was measured using a pH meter. Following wet digestion, the TP was quantified using a FIAstar™ 5000 flow injection analyzer. All measurements were repeated in triplicate.
Soil DNA extraction
The total microbial group DNA from different samples was extracted using an EZNA® Stool DNA Kit (D5625-02, Omega Bio-Tek, Inc., Norcross, GA, USA) according to the manufacturer’s instructions. The total DNA was eluted in 50 μL of elution buffer and stored at -80 °C until measured by PCR from LC-Bio Technology (Hangzhou, China), with the isolation being confirmed by 1.2% agarose gel electrophoresis. After amplified by PCR, followed by purification, quantification, and library quality control. High throughput sequencing was performed with HiSeq4000® (Illumina Inc., San Diego, CA, USA) once the library was qualified. The sequencing mode was PE150.
Sequence analysesPreprocessing of sequence data
Raw data was preprocessed by the software Cutadapt ver. 1.9 (Martin 2011) for removing contained joints and low-quality data to obtain clean data for subsequent analysis. Subsequently, the sequence reads were scanned by the window method with 6 bp defaults. When the average quality value in the window was less than 20, the reads were cut off from window start to the 3′ end. Sequences of less than 60 bp in length, and sequences with an N content > 5% following truncation were abandoned. Finally, the host sequence contamination was removed (Kanehisa & Goto 2000).
Genome assembly
We initially assembled the genome of each sample using the software metaSPAdes ver. 3.10.0 (Bankevich et al. 2012). After a single sample was assembled, the unused reads of each sample were merged into a MIX sample for mixed assembly. The assembled result was a FASTA format file, in which a sequence (contig) was reserved for subsequent clustering analysis with contigs longer than 500 bp. Once the genome assembly was complete, we proceed to assemble the resulting evaluation. According to the assembly result, we employed the QUAST program matched with SPAdes to evaluate the result of the assembly (Chen & Pachter 2005, Gurevich et al. 2013).
Gene prediction and abundance analysis
We used MetaGeneMark® for Linux 64bit (Gene Probe Inc., Atlanta, GA, USA) to predict the CDS (Coding Region) of each sample and mixed assembly contigs (≥ 500 bp), and filtered the sequence with CDS lengths of less than 100 nt using the software CD-HIT ver. 4.6.1 (Li & Godzik 2006). Next, according to the CDS prediction results, non-redundant genes were obtained using CD-HIT software. Simultaneously, 95% identity and 90% coverage were used for clustering and the selection of the longest sequence as representative sequence. The software Bowtie2 ver. 2-2.2.0 (Langmead & Salzberg 2012 - http://bowtie-bio.sourceforge.net/bowtie2/index.shtml) was used to align the clean data of each sample to gene sequences, and the number of reads of each gene of every sample was calculated. Unigenes for subsequent analysis were obtained once the number of reads ≤ 2 genes were filtered out and compared between the samples.
According to the number of reads and the length of each gene, the abundance information (TPM value) of each gene in each sample is calculated. The formula was as follows (eqn. 1):
where G represents the abundance of a unigene, K represents the K-th unigene, r is the number of reads for comparison of the above genes, and L is the length of the gene.
Species annotation
We compared the detected unigenes with the NR_mate Library to obtain species annotation information using the software Bowtie2. A rarefaction curve was developed to compare the species richness of the samples with different quantities of sequencing data, and to show the rationality of the amount of sequencing data. When the extraction sequence was 468.935, the rarefaction curves of different soil samples exhibited a tendency to flatten, which indicated that these sample sizes were sufficient to cover all bacteria, and the sequencing data covered most of the bacterial information in the sample (Fig. S2 in Supplementary material).
Data analysis
The alpha diversity indices included richness, Shannon-Weiner index (H′ - eqn. 2), and Simpson’s index (D - eqn. 3), and Chao1 index (eqn. 4), which were calculated at the species level. Richness refers to the sum of bacterial species in the soil. The diversity indices were calculated as follows:
where Ni is the number of individuals of the i-th species, N is the number of individuals of all species in the community, S is the overall number of species, and Pi is the proportion of individuals belonging to species i in all individuals; Sobs is the number of individuals of all species, F1 and F2 are the singleton and doubleton counts, respectively.
We tested the effects of the CK, N, and N+P treatments and forest age, and their interactions on the abundance of soil bacteria, alpha diversity (richness, Chao1, Shannon’s index, and Simpson’s index), and the abundance of groups with the relative abundances greater than 1% in the total bacterial community by performing two-way analysis of variance (ANOVA). We assessed the assumptions of normality for all models using the Shapiro test, histograms, and QQ normal plots of model residuals. We employed Levene’s test for the assumption of homogeneity. These analyses confirmed that both assumptions were met for all models.
To examine the effects of treatment and age on bacterial composition, we conducted a permutational multivariate analysis of variance (perMANOVA) implemented by the “adonis” function of the “vegan” package in R (Oksanen et al. 2017). In perMANOVA, we used a Bray-Curtis dissimilarity matrix to summarize species composition and employed 999 permutations to determine statistical significance. We visualized the compositional data using non-metric multidimensional scaling with a Bray-Curtis dissimilarity measure.
Further, we examined how the bacterial composition was linked to soil properties by using the function “envfit” of the “vegan” package (Oksanen et al. 2017). Soil parameters, which have significant effects (P<0.05) on bacteria as determined by the “envfit” function, were further used to examine the associations with soil bacterial abundance, alpha diversity and dominant phylum by Pearson’s correlation analysis. All analyses were performed in R ver. 3.6.0 (R Development Core Team 2019).
ResultsEffects of N and P addition on the abundance of bacterial communities
The effect of fertilization on the abundance of soil bacteria was strongly dependent on forest age (Tab. 1, Fig. 1). For CK, the bacterial abundance was lower in the 78-year-old stand than in the 46- and 200-year-old stands (Fig. 1). Compared with CK, both N and N+P treatments decreased bacterial abundance in the 46-year-old stand, with a greater reduction under the N+P treatment (Fig. 1). In the 78-year-old stand, the N+P treatment increased the bacterial abundance, while no effect appeared in the N treatment. However, in the 200-year-old stand, neither N nor N+P treatment apparently affected bacterial abundance.
Effects of N and P addition on alpha diversity of bacteria
The species richness and Chao1 in CK increased with stand age, while Shannon’s and Simpson’s indices were lower in the 46-year-old stands compared with the older stands (Tab. 1, Fig. 2). Both N and N+P treatments increased species richness and Chao1 in the 46- and 78-year-old stands, with the greatest increases observed in the 46-year-old stand; however, increases in 200-year-old stands were far less pronounced. Similarly, both N and N+P treatments increased Shannon’s and Simpson’s indices in the 46-year-old stands but had little effect in the older stands (Tab. 1, Fig. 2).
Effects of N and P additions on the composition of bacterial communities
In different stands (Fig. 3a, Fig. 3b, Fig. 3c), the bacterial composition was dominated (>1%) by the phylum Proteobacteria (52-77%), followed by Acidobacteria (11-29%). The most frequent genera (>1%) for all stands were Bradyrhizobim (9-15%), Burkholderia (5-8%, except 6-34% in CK of the 46-year-old). The most prevalent bacterial species were Candidatus_Koribacter_versatilis (2-5%), and Candidatus_Solibacter_usitatus (2-4%).
The overall dominant bacterial composition, at the phylum, genus, and species levels, varied strongly with forest age and treatment (Tab. 2; Fig. S3, Fig. S4, Fig. S5 in Supplementary material). The abundance of Proteobacteria in 46-year-old stands was 58.4% and 46.5% higher than in 78- and 200-year-old stands. The phylum Verrucomicrobia, Cyanobacteria, Chloroflexi, Gemmatimonadetes, and genera Candidatus_Koribacter, Candidatus_Solibacter, Edaphobacter, Chthoniobacter in 78-year-old were higher than in the other two stands.
In CK, the bacterial communities were more concentrated by the most abundant plyla and genera (e.g., Proteobacteria, Bradyrhizobyum, Burkholderia) in the 46-year-old stand than in the older stands. The responses of the same dominant bacterial communities to fertilization was contingent on forest age (Tab. S2, Tab. S3 in Supplementary material). It was observed that fertilization decreased Proteobacteria, and increased Actinomycetes, Verrucomicrobia, Firmicutes, Cyanobacteria, and Chloroflexi in the 46-year-old stands, while decreased Verrucomicrobia and increased Planctomycetes in the 78-year-old stands, and did not change at the phylum level in the 200-year-old stands.
At the genus level, fertilization increased Candidatus_Koribacter, Candidatus_Solibacter in 46- and 200-year-old stands, while both decreased in the 78-year-old. Moreover, fertilization increased Acidobacterium, Edaphobacter, Granulicella and decreased Burkholderia in the 46-year-old stands, reduced Bradyrhizobim, Burkholderia, Edaphobacter, and Chthoniobacter in 78-year-old stands, and decreased Acidobacterium and Granulicella in 200-year-old stands.
At the species level, fertilization increased Acidobacteriaceae_bacterium_URHE0068, Rhodospirillales_bacterium_URHD0088, and decreased Bradyrhizobium_sp_Tv2a2 in the 46- and 200-year-old stands. Moreover, Candidatus_Koribacter_versatilis and Candidatus_Solibacter_usitatus were reduced in 78-year-old stands but increased in 200-year-old stands by fertilization.
In CK, bacterial communities (phylum, genus, species level) formed independent groups by forest age, and the 46-year-old were grouped on the left and older sites on the right of axis 1 in the ordination space (Fig. 4a, Fig. 4b, Fig. 4c), which indicated that the younger stands possessed different overall species composition from the older stands. Moreover, three distinct fertilization type groups were prominently separated in the 46- and 78-year-old stands. For the same forest age, the degree of cluster separation formed by the treatments gradually decreased, which suggested the fertilization-induced effect was reduced with forest age.
Effects of soil physicochemical and vegetation factors on bacterial communities
The SOC/TN, pH, BD and soil sand content decreased, while the soil moisture, DOC, DON and soil silt content increased over time with more pronounced changes between stand ages (Fig. S6 in Supplementary material). Our result also indicated that changes of soil properties by fertilization varied between stand ages (Fig. S6 in Supplementary material). Further analysis revealed that in terms of bacterial species distribution, the interpretation rates of different physical and chemical factors were as follows: soil moisture (R2 = 0.70, P < 0.01) > DON (R2 = 0.61, P < 0.01) > DOC (R2 = 0.57, P < 0.01) > sand, silt (R2 = 0.58, P < 0.01) > pH, BD (R2 = 0.54, P < 0.001 - Tab. 3, Fig. 3d, Fig. 3e, Fig. 3f). Based on the different background, the effect of soil physiochemical properties on bacterial community changed among forest ages under fertilization. Correlation analysis (Tab. S2, Tab. S3 in Supplementary material) showed that soil moisture, NO3--N, DOC and pH were the main influencing factors in the 46-year-old, but these factors had little effect on 78- and 200-year-old stands among treatments.
Discussion
Molecular analyses provided insights into the fertilization effects on soil bacteria in forest ecosystems. In this study, the soil bacterial diversity in the CK showed an overall upward trend with forest age, which might be because older forests, due to the higher vegetation diversity, contain more plentiful available niches to promote the coexistence of more species of soil bacteria (Li et al. 2019, Qu et al. 2020). Many previous studies reported that the bacterial community diversity clearly declines with increasing N addition (Zeng et al. 2016, Kaspari et al. 2017, Wang et al. 2018b). In this study, such decline was only reflected by the Simpson’s and Shannon’s indices of the 200-year-old within the alpha index we considered, but this negative effect was relatively small (P > 0.05). Our result confirmed previous reports on the occurrence of different responses of soil bacterial abundance to fertilization between forest ages, with a significant decrease in young stands whilst the older stands had more stable bacterial communities (Frey et al. 2014, Maaroufi et al. 2015, Wu et al. 2019).
The majority of the unigenes detected fell into four bacterial phylum categories. Proteobacteria constantly maintained the absolute dominance across different stand ages (Zhang et al. 2018a), followed by Acidobacteria and Actinobacteria, all of which are considered to be universally present in acidic soils (Wang et al. 2018c). In our study, a decrease of Proteobacteria with increasing stand age was found, which was in contrast with previous studies reporting an increasing abundance of Proteobacteria along the succession (Jangid et al. 2013). Previous studies revealed that soil bacteria of the same trophic type did not show uniformity subsequent to the nutrient addition for different vegetation types (Yuan & Chen 2016, Nie et al. 2018, Li et al. 2019), which was further confirmed by our result. Although, a nutritional theory was proposed that the enhanced N availability increased the abundance of copiotrophic bacterial taxa, including Proteobacteria and Bacteroidetes, but lowered the proportion of oligotrophic Acidobacteria (Nie et al. 2018), which is further confirmed in the existing study of evergreen broad-leaved forests (Yang et al. 2020). However, the changes in bacterial species following fertilization in our study did not completely support the nutritional theory among forest ages, as we observed a decrease in the abundance of Proteobacteria (copiotrophic) in 46-year-old stands, and no significant change in older stands following N addition.
It was reported that fertilization results in a decline in diversity by stimulating the expansion of nitrophilous species and the competitive exclusion of others (Bobbink et al. 2010, Nie et al. 2018). Compared with major phylum, the N-transformation-related phyla (e.g., Proteobacteria, Planctomycetes, Cyanobacteria) had no advantages as the result of N addition. The point was supported by several studies that the availability of NH4+-N rather than NO3--N was primarily responsible for the modification of bacterial diversity (Zeng et al. 2016, Nie et al. 2018). Although NH4+-N is the preferred nutrient for most bacteria (Merrick & Edwards 1995), the change of NO3--N was larger than that of NH4+-N after fertilization in this study, which may also be the reason for the different results. Chen et al. (2016) observed that soil total N and NO3--N were the main factors that influenced most variations of bacterial community composition in the agricultural ecosystems.
In this study, the effect of fertilization on soil bacterial community decreased with forest age. Bacterial community structure was altered by N treatment in 46-year-old stands, and the effects of N and N + P treatments differed significantly from CK. A weak effect of nutrient additions on soil bacterial diversity and community composition was found in 200-year-old stands in this study, which is similar to the results between primary and secondary tropical rainforests (Li et al. 2019). Soil texture, nutrient availability and acidification may account for the variable response of soil bacteria to N addition (Mao et al. 2017). The differences in soil environmental regime between the different successional stages rather than chronic N fertilization might be the major control in soil bacterial composition in forest ecosystems. At the same time, the response of the dominant bacteria phyla to long-term N addition was distinct from that in short-term between forest ages (Cui et al. 2017, Han et al. 2018, Teng et al. 2018). At our sites, the abundance of Proteobacteria was inhibited after 4-year N and N+P additions in both 46- and 200-year-old stands, while only occurred in the 46-year-old stands after 7-year N additions. This might reflect the acclimation to long-term high N deposition (Han et al. 2018) and the elevated P availability with P addition, thus alleviating the disadvantages of chronic N enrichment (Wang et al. 2018b).
Soil texture, moisture, NO3--N, DOC, and pH are important factors for predicting the distribution of soil bacterial composition, along with fertilization among forest ages (Wang et al. 2018a, Li et al. 2019, Zhao et al. 2019). Our data revealed that fertilization-influenced soil bacteria mainly via change of soil moisture, NO3--N, DOC and pH in the 46-year-old stands, which had less effects in the older stands. Results from Wu et al. (2019) confirmed a strongly negative correlation between soil bacterial alpha diversity, main phyla and pH. Furthermore, soil microbial communities were likely to be altered with the variations in soil moisture (Wu et al. 2019). The significant increase of soil moisture and decrease of soil pH with forest age might alleviate the effects to soil bacteria by fertilization (Wang et al. 2018b). In the 46-year-old, the increase of bacterial diversity induced by soil moisture may be due to the decrease of the abundance of Proteobacteria and the increase of the abundance of other dominant flora, while the influence of pH is converse, but this regulation was not found in the older forest.
Background differences in soil C and N status between forest types have been thought to be possible causes for dissimilarity in soil bacterial community composition in tropical rainforests (Li et al. 2019). Fertilization might indirectly affect soil microbes by modifying the soil C availability (Wang et al. 2018a). When SOC/TN was >10, TN could be limited to soil microbes (Wang et al. 2018a). In our sites, the SOC/TN was ranged from 12.32 to 15.38, so that TN might be a limiting nutrient to soil bacteria. Results from Wang et al. (2018a) showed that TP significantly affected soil bacterial communities, while our study showed that N+P addition did not cause pronounced changes in soil bacterial communities compared with N addition (Wang et al. 2018c). The responses of soil microbes to P addition might vary with availability of other soil nutrients, as Esberg et al. (2010) showed that microbial biomass and activity were more likely to be limited by P when C availability is high. Our results revealed that compared with N addition, the excessive P addition (N+P) exhibited negative effects for Proteobacteria in 46-year-old stands, as well as Cyanobacteria and Gemmatimonadetes, Verrucomicrobia, and Acidobacteria in 78-year-old stands. These results suggested that soil bacterial communities might not encounter P limitation following the long-term addition of P in N-rich forest ecosystems (Wang et al. 2018b, Wang et al. 2018c).
Conclusions
The abundance of bacteria showed a concave distribution with forest age, while bacterial diversity showed an overall uptrend with ages in a chronosequence of subtropical forests. Soil bacterial community became more stable and resilient to nutrient additions with increasing of the stand age, with a significant decrease occurring only in the youngest stand within the age-sequence from 46- to 200-year-old. Moreover, the same soil bacterial group did not respond to nutrient addition consistently at different stand ages. Nitrogen transformation-related phylum had no advantages in N enrichment compared to others. Further, a prolonged N deposition can likely induce soil bacteria responses that are different from the short-term deposition. In this study, the driving factors of soil bacterial communities were forest age, followed by nutrient addition. The core bacterial communities were mainly controlled by soil texture, moisture, DOC, DON, pH, SOC/TN and TP among different stand ages.
Acknowledgements
This study was financially supported by the National Natural Science Foundation of China (NSFC, No. 31770672 and 31370626) and the National Basic Research Program of China (2010CB950602). We gratefully acknowledge Professor. Han Y. H. Chen for his help in writing the manuscript. Thank the Forestry Bureau of Shitai County, and the Zhawan Nature Reserve, Anhui Province, for their providing fertilization experiment forests. We are thankful to Mr. Sun MY, Li PX, Wang F and Ma RH for sample collection and/or laboratory analysis. We also thank anonymous reviewers for comments which substantially improved this manuscript.
ReferencesBankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Son P, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PASpades: a new genome assembly algorithm and its applications to single-cell sequencing. Journal of Computational Biology 19: 455-477.2012Bobbink R, Hicks K, Galloway J, Spranger T, Alkemade R, Ashmore M, Bustamante M, Cinderby S, Davidson E, Dentener F, Emmett B, Erisman JW, Fenn M, Gilliam F, Nordin A, Pardo L, De Vries WGlobal assessment of nitrogen deposition effects on terrestrial plant diversity: a synthesis. Ecological Applications 20: 30-59.2010Chen C, Zhang JN, Lu M, Qin C, Chen YH, Yang L, Huang QW, Wang JC, Shen ZG, Shen QRMicrobial communities of an arable soil treated for 8 years with organic and inorganic fertilizers. Biology and Fertility of Soils 52: 455-467.2016Chen K, Pachter LBioinformatics for whole-genome shotgun sequencing of microbial communities. PLoS Computational Biology 1: 106-112.2005Cui J, Wang J, Xu J, Xu C, Xu XChanges in soil bacterial communities in an evergreen broad-leaved forest in east China following 4 years of nitrogen addition. Journal of Soils and Sediments 17: 2156-2164.2017Dai YZ, Wang HR, Li AQ, Zhu YY, Wang XP, Fang Z, Xu XNDynamical pattern of soil microbial biomass and its response to nitrogen and phosphorus additions in a subtropical evergreen broad-leaved forest. Ecology and Environmental Sciences 27: 1395-1404. [in Chinese with English Abstract]2018Elser JJ, Bracken MES, Cleland EE, Gruner DS, Harpole WS, Hillebrand H, Ngai JT, Seabloom EW, Shurin JB, Smith JEGlobal analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems. Ecological Letters 10: 1135-1142.2007Esberg C, Toit B, Olsson R, Ilstedt U, Giesler RMicrobial responses to P addition in six south African forest soils. Plant and Soil 329: 209-225.2010Frey SD, Ollinger S, Nadelhoffer K, Bowden R, Brzostek E, Burton A, Caldwell BA, Crow S, Goodale CL, Grandy AS, Finzi AC, Kramer MG, Lajtha K, LeMoine J, Martin M, McDowell WH, Minocha R, Sadowsky JJ, Templer PH, Wickings KChronic nitrogen additions suppress decomposition and sequester soil carbon in temperate forests. Biogeochemistry 121: 305-316.2014Gravuer K, Eskelinen A, Winbourne JB, Harrison SPVulnerability and resistance in the spatial heterogeneity of soil microbial communities under resource additions. Proceedings of the National Academy of Sciences USA 117: 7263-7270.2020Gurevich A, Saveliev V, Vyahhi N, Tesler GQuast: quality assessment tool for genome assemblies. Bioinformatics 29: 1072-1075.2013Han XG, Shen WJ, Zhang JB, Muller CMicrobial adaptation to long-term N supply prevents large responses in N dynamics and N losses of a subtropical forest. Science of the Total Environment 626 (971): 1175-1187.2018Hardoim PR, van Overbeek LS, Berg G, Pirttilä AM, Compant S, Campisano A, Döring M, Sessitsch AThe hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiology and Molecular Biology Reviews 79 (3): 293-320.2015Hunt N, Gilkes RFarm monitoring handbook - a practical down-to-earth manual for farmers and other land users. University of Western Australia, Nedlands, WA, and Land Management Society Como, WA, Australia, vol. 496, pp. 198-205.1992Jangid K, Whitman WB, Condron LM, Turner BL, Williams MASoil bacterial community succession during long-term ecosystem development. Molecular Ecology 22: 3415-3424.2013Kanehisa M, Goto SKEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Research 28: 27-30.2000Kaspari M, Bujan J, Weiser MD, Ning D, Michaletz ST, He Z, Enquist BJ, Waide RB, Zhou J, Turner BL, Wright SJBiogeochemistry drives diversity in the prokaryotes, fungi, and invertebrates of a Panama forest. Ecology 98: 2019-2028.2017Langmead B, Salzberg SLFast gapped-read alignment with Bowtie 2. Nature Methods 9: 357-359.2012Li W, Godzik ACd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658-1659.2006Li P, Shen CC, Jiang L, Feng Z, Fang JDifference in soil bacterial community composition depends on forest type rather than nitrogen and phosphorus additions in tropical montane rainforests. Biology and Fertility of Soils 55: 313-323.2019Liu X, Zhang Y, Han W, Tang A, Shen J, Cui Z, Vitousek P, Erisman JW, Goulding K, Christie P, Fangmeier A, Zhang FEnhanced nitrogen deposition over China. Nature 494: 459-462.2013Llado S, Lopez-Mondejar R, Baldrian PForest soil bacteria: diversity, involvement in ecosystem processes, and response to global change. Microbiology and Molecular Biology Reviews 81: e00063-16.2017Lue C, Tian HSpatial and temporal patterns of nitrogen deposition in China: synthesis of observational data. Journal of Geophysical Research 112: D22S05.2007Maaroufi NI, Nordin A, Hasselquist NJ, Bach LH, Palmqvist K, Gundale MJAnthropogenic nitrogen deposition enhances carbon sequestration in boreal soils. Global Change Biology 21: 3169-3180.2015Mao Q, Lu X, Zhou K, Chen H, Zhu X, Mori T, Mo JEffects of long-term nitrogen and phosphorus additions on soil acidification in an N-rich tropical forest. Geoderma 285: 57-63.2017Marketa SM, Marek O, Ladislav C, Zdenek K, Jana O, Evelyn H, Jan K, Franz HJA, Hadacek EMicrobial communities show parallels at sites with distinct litter and soil characteristics. Applied and Environmental Microbiology 77: 7560-7567.2011Martin MCutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet Journal 17 (1): 10-12.2011Merrick MJ, Edwards RANitrogen control in bacteria. Microbiological Reviews 59 (4): 604-622.1995Nie YX, Wang MC, Zhang W, Ni Z, Hashidoko Y, Shen WAmmonium nitrogen content is a dominant predictor of bacterial community composition in an acidic forest soil with exogenous nitrogen enrichment. Science of the Total Environment 624: 407-415.2018Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Henry M, Stevens H, Szoecs E, Wagner HVegan: community ecology package. Web site.2017Qu ZL, Liu B, Ma Y, Sun HDifferences in bacterial community structure and potential functions among Eucalyptus plantations with different ages and species of trees. Applied Soil Ecology 149: 103515.2020Ramirez KS, Craine JM, Fierer NConsistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Global Change Biology 18: 1918-1927.2012Reed SC, Cleveland CC, Townsend ARFunctional ecology of free-living nitrogen fixation: a contemporary perspective. Annual Review of Ecology Evolution and Systematics 42: 489-512.2011R Development Core TeamR: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.2019Singh BK, Bardgett RD, Smith P, Reay DSMicroorganisms and climate change: terrestrial feedbacks and mitigation options. Nature Reviews Microbiology 8: 779-790.2010Teng Z, Cui J, Wang JJ, Fu XL, Xu XXEffect of exogenous nitrogen and phosphorus inputs on the microbe-soil interaction in the secondary Castanopsis sclerophylla forest in east China. iForest 11: 794-801.2018Thien SJA flow diagram for teaching texture by feel analysis. Journal of Agronomy Education 8: 54-55.1979Tian D, Jiang L, Ma SH, Fang WJ, Schmid B, Xu LC, Zhu JX, Li P, Losapio G, Jing X, Zheng C, Shen H, Xu X, Zhu B, Fang JYEffects of nitrogen deposition on soil microbial communities in temperate and subtropical forests in China. Science of the Total Environment 607: 1367-1375.2017Uroz S, Oger P, Lepleux C, Collignon C, Frey-Klett P, Turpault MPBacterial weathering and its contribution to nutrient cycling in temperate forest ecosystems. Reserch in Microbiology 162: 820-831.2011Vitousek PM, Porder S, Houlton BZ, Chadwick OATerrestrial phosphorus limitation: mechanisms, implications, and nitrogen-phosphorus interactions. Ecological Applications 20: 5-15.2010Wang C, Liu D, Bai EDecreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biology and Biochemistry 120: 126-133.2018aWang H, Liu S, Zhang X, Mao Q, Li X, You Y, Wang J, Zheng M, Zhang W, Lu X, Mo JNitrogen addition reduces soil bacterial richness, while phosphorus addition alters community composition in an old-growth N-rich tropical forest in southern China. Soil Biology and Biochemistry 127: 22-30.2018bWang QF, Cong W, Yu WW, Turak A, Chen D, Ying H, Ao J, Yong J, Huang ZEffects of nitrogen and phosphorus inputs on soil bacterial abundance, diversity, and community composition in Chinese fir plantations. Frontiers in Microbiology 9: 1543.2018cWardle DAA comparative assessment of factors which influence microbial biomass carbon and nitrogen levels in soil. Biological Reviews 67: 321-358.2010Weand MP, Arthur MA, Lovett GM, Mcculley RL, Weathers KCEffects of tree species and N additions on forest floor microbial communities and extracellular enzyme activities. Soil Biology and Biochemistry 42: 2161-2173.2010Wu J, Liu W, Zhang W, Shao Y, Duan H, Chen B, Wei X, Fan HLong-term nitrogen addition changes soil microbial community and litter decomposition rate in a subtropical forest. Applied Soil Ecology 142: 43-51.2019Xu J, Li M, Liu L, Tian D, Xu XNEffects of nitrogen deposition on litterfall production and its nutrients in subtropical evergreen broad-leaved forests. Bulletin of Soil and Water Conservation 36: 113-118. [in Chinese with English Abstract]2016Yang L, Wang N, Chen Y, Yang W, Tian DS, Zhang CY, Zhao XH, Wang JS, Niu SLCarbon management practices regulate soil bacterial communities in response to nitrogen addition in a pine forest. Plant and Soil 452: 137-151.2020Yuan ZY, Chen HYNegative effects of fertilization on plant nutrient resorption. Ecology 926: 373-380.2016Zeng J, Liu XJ, Song L, Lin XG, Zhang HY, Shen CC, Chu HYNitrogen fertilization directly affects soil bacterial diversity and indirectly affects bacterial community composition. Soil Biology and Biochemistry 92: 41-49.2016Zhang HF, Liu HM, Zhao JN, Wang LL, Li G, Huang FC, Wang H, Lai X, Li J, Yang DElevated precipitation modifies the relationship between plant diversity and soil bacterial diversity under nitrogen deposition in Stipa baicalensis steppe. Applied Soil Ecology 119: 345-353.2017Zhang KR, Cheng XL, Shu X, Liu Y, Zhang QFLinking soil bacterial and fungal communities to vegetation succession following agricultural abandonment. Plant and Soil 431: 19-36.2018aZhang TA, Chen HYH, Ruan HGlobal negative effects of nitrogen deposition on soil microbes. The ISME Journal 12: 1817-1825.2018bZhao FZ, Bai L, Wang JY, Deng J, Ren CJ, Han XH, Yang GH, Wang JChange in soil bacterial community during secondary succession depend on plant and soil characteristics. Catena 173: 246-252.2019Zhou Z, Wang C, Zheng M, Jiang L, Luo YPatterns and mechanisms of responses by soil microbial communities to nitrogen addition. Soil Biology and Biochemistry 115: 433-441.2017
Interactive effects of treatment and stand age on the abundance of soil bacteria (in terms of number of unigenes detected). Treatments are control (CK), nitrogen addition (N), and N and phosphorus addition (N+P). Values are bootstrapped means ± 95% confidence intervals (CIs) for individual samples and their sum (n=9). Differences between stand ages or treatments are significant (α = 0.05) when their CIs do not overlap.
Interactive effects of treatment and stand age on soil bacterial alpha diversity. Treatments are control (CK), nitrogen addition (N), and N and phosphorus addition (N+P). Values are bootstrapped means ± 95% confidence intervals (CIs) for individual samples and their sum (n=9). Differences between stand ages or treatments are significant (α = 0.05) when their CIs do not overlap.
Relative abundances of the soil bacterial groups (a: phylum level; b: genus level; c: species level) in each plot (n=9). Only groups with relative abundance > 1% are shown, while those with < 1% relative abundance are grouped into “others”.
Non-metric multidimensional scaling ordination of the composition of all bacterial phyla, genus and species of soil microbial communities. Close sites in the ordination space have similar bacterial assemblages, whereas those farther apart are less similar. Contours indicate time since reclamation. The length of vectors indicate the strength of associations; red vectors are significant (P<0.05).
Interactive effects of treatment and stand age on soil bacterial abundance, alpha diversity, and composition (n=9). P-values from ANOVA are displayed.
Characteristic
Age
Treatment
Age × Treatment
Abundance
<0.001
0.076
<0.001
Richness
<0.001
<0.001
0.001
Chao1
<0.001
<0.001
0.111
Shannon’s index
0.148
<0.001
<0.001
Simpson’s index
0.110
<0.001
<0.001
Interactive effects of treatment and stand age on soil bacterial phylum, genus, and species (n=9). P-values from ANOVA are displayed.
Taxa
Age
Treatment
Age × Treatment
Phylum
<0.001
0.009
<0.001
Genus
<0.001
<0.001
<0.001
Species
<0.001
<0.001
0.111
Results of soil properties linked to the soil bacterial composition (at phylum, genus, and species level) using the function “envfit” in the vegan package (n=9). R2 and Pr(>r) represent the coefficient of determination and probability of significance, respectively.
Variable
Phylum
Genus
Species
R2
Pr (>r)
R2
Pr (>r)
R2
Pr (>r)
Moisture content
0.75
0.001
0.64
0.001
0.7
0.001
BD
0.59
0.001
0.51
0.002
0.54
0.001
Sand
0.42
0.004
0.54
0.001
0.58
0.001
Silt
0.43
0.002
0.5
0.001
0.58
0.001
Clay
0.04
0.633
0.01
0.98
0.01
0.97
NH4+-N
0.02
0.835
0.06
0.442
0.3
0.692
NO3--N
0.12
0.216
0.2
0.068
0.23
0.047
DOC
0.65
0.001
0.6
0.001
0.57
0.001
DON
0.6
0.001
0.58
0.001
0.61
0.001
TP
0.24
0.033
0.43
0.002
0.47
0.002
TN
0.07
0.432
0.17
0.99
0.13
0.192
SOC
0.15
0.168
0.24
0.046
0.19
0.083
SOC/TN
0.37
0.007
0.3
0.015
0.33
0.006
pH
0.63
0.001
0.51
0.001
0.54
0.001
Fig. S1 - The map of plot distribution.
Fig. S2 - Rarefaction curves of bacterial communities.
Fig. S3 - Interactive effects of treatment and stand age on the unigenes abundance of soil bacterial dominate phylum.
Fig. S4 - Interactive effects of treatment and stand age on the unigenes abundance of soil bacterial dominate genus.
Fig. S5 - Interactive effects of treatment and stand age on the unigenes abundance of soil bacterial dominate species.
Fig. S6 - Interactive effects of treatment and stand age on soil parameters.
Tab. S1 - Soil texture analysed by the fractionation test described by Hunt & Gilkes (1992).
Tab. S2 - Pearson correlation coefficients between soil properties, bacterial abundance and its alpha diversity indices.
Tab. S3 - Pearson correlation analysis between soil bacterial dominate phylum (the relative abundance of Unigenes >0.01) and soil factors.