Nitrogen deposition and management practices increase soil microbial biomass carbon but decrease diversity in Moso bamboo plantations.

Because microbial communities play a key role in carbon (C) and nitrogen (N) cycling, changes in the soil microbial community may directly affect ecosystem functioning. However, the effects of N deposition and management practices on soil microbes are still poorly understood. We studied the effects of these two factors on soil microbial biomass carbon (MBC) and community composition in Moso bamboo plantations using high-throughput sequencing of the 16S rRNA gene. Plantations under conventional (CM) or intensive management (IM) were subjected to one of four N treatments for 30 months. IM and N addition, both separately and in combination, significantly increased soil MBC while decreasing bacterial diversity. However, increases in soil MBC were inhibited when N addition exceeded 60 kg N∙ha(-1)∙yr(-1). IM increased the relative abundances of Actinobacteria and Crenarchaeota but decreased that of Acidobacteria. N addition increased the relative abundances of Acidobacteria, Crenarchaeota, and Actinobacteria but decreased that of Proteobacteria. Soil bacterial diversity was significantly related to soil pH, C/N ratio, and nitrogen and available phosphorus content. Management practices exerted a greater influence over regulation of the soil MBC and microbial diversity compared to that of N deposition in Moso bamboo plantations.

As the most abundant organisms on earth, microbes play a key role in natural ecosystems, including the biogeochemical cycles of n class="Chemical">carbon (C) aclass="Chemical">nd class="Chemical">n class="Chemical">nitrogen (N) and the biodegradation or stabilization of environmental contaminants123. The soil microbial biomass has an important role in nutrient cycling and is therefore essential for plant growth4. Changes in the soil microbial community may directly affect soil ecosystem functioning, particularly C and N cycling5. In recent years, high-throughput sequencing technologies have been widely used to analyze the species composition and functional diversity of microbial populations under different fertilization regimes678. These technologies provide a more detailed description of microbial communities than traditional methods, such as cloning libraries, denaturing gradient gel electrophoresis (DGGE), and phospholipid fatty acid (PLFA) analysis and has proven to be a very powerful technique in microbial ecology research9. Because soil microbes are sensitive to environmental change1011, a comprehensive analysis of the microbial community structure and microbial responses to environmental change using high-throughput sequencing technologies could improve our understanding of biogeochemical cycles in natural or managed ecosystems12.

n class="Species">Moso bamboo (class="Chemical">n class="Species">Phyllostachys pubescens Mazel ex H. de Lehaie), distributed widely in subtropical China13, grows rapidly and is the most important source of non-wood forest products in China1415. Under typical conventional management (CM), trunks and shoots of Moso bamboo are harvested regularly without any other management interventions15. In recent decades, intensive management (IM) practices, such as removing understory vegetation, plowing, and fertilization, have been widely adopted in order to maximize economic returns1315. It has been reported that these IM practices alter the soil microbial biomass and community diversity161718. For example, Liu et al.19 determined via polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) that fertilizer application significantly increased rhizosphere microbial diversity in a Moso bamboo forest. In contrast, Sun et al.20 observed that long-term conventional management of bamboo plantations did not significantly alter soil microbial biomass and diversity. However, few studies have used high-throughput sequencing technologies to assess the effects of management practices on soil microbial biomass and community diversity in Moso bamboo plantations.

Over the last three decades, the average bulk deposition of N has rapidly increased across China21, particularly in subtropical China2223. n class="Species">Moso bamboo placlass="Chemical">ntatioclass="Chemical">ns are primarily distributed across the regioclass="Chemical">n with the greatest N depositioclass="Chemical">n iclass="Chemical">n Chiclass="Chemical">na, accordiclass="Chemical">ng to both curreclass="Chemical">nt estimates aclass="Chemical">nd future predictioclass="Chemical">ns24. Research has showclass="Chemical">n that N additioclass="Chemical">n may reduce soil microbial biomass2526 aclass="Chemical">nd commuclass="Chemical">nity diversity27. class="Chemical">n class="Species">Compton et al.28 and Frey et al.29 found that chronic N supplementation decreased soil microbial biomass at the Harvard Forest. However, Paul et al.30 found that the microbial biomass was greater under high N input (300 kg N∙ha−1∙yr−1) than under low N input (100 kg N∙ha−1∙yr−1). Johnson et al.31 also observed that N addition (80 kg N∙ha−1∙yr−1) increased the soil microbial biomass over seven years in a simulated N deposition study. Boxman et al.32 found that N deposition did not affect soil microbial biomass. The reasons for these differences in results remain unclear and require further study. The lack of understanding of the effects of N deposition on soil microbial quantity, community composition, and microbial diversity in Moso bamboo plantations with different management practices limits our ability to predict the complex response of plantation ecosystems to global environmental changes.

In the present study, we employed high-throughput sequencing to analyze the effects of n class="Chemical">nitrogen depositioclass="Chemical">n aclass="Chemical">nd maclass="Chemical">nagemeclass="Chemical">nt practices oclass="Chemical">n soil microbial commuclass="Chemical">nity aclass="Chemical">nd diversity iclass="Chemical">n class="Chemical">n class="Species">Moso bamboo plantations. The objectives of this study were to test the following hypotheses: (1) IM practices increase soil microbial biomass carbon (MBC) and microbial diversity; (2) N deposition increases soil MBC and diversity; and (3) the combination of IM and N deposition exerts a greater effect on soil MBC and community composition than each practice independently.

Results

Soil microbial biomass carbon

Through assessment of soil samples from n class="Disease">CM aclass="Chemical">nd IM plots iclass="Chemical">n the class="Chemical">n class="Species">Moso bamboo plantation, we found that soil MBC content was significantly higher under IM than under CM when no N was added (Fig. 1). Soil MBC content increased significantly with the addition of N under both CM and IM, but decreased significantly when the N addition rate exceeded 30 kg N∙ha−1∙yr−1 (N30 treatment) under IM and 60 kg N∙ha−1∙yr−1 (N60 treatment) under CM (Fig. 1).

Figure 1

Soil microbial biomass carbon (MBC) contents under conventional management (CM) or intensive management (IM) and three N addition treatments.N30: 30 kg N∙ha−1∙yr−1; N60: 60 kg N∙ha−1∙yr−1; N90: 90 kg N∙ha−1∙yr−1.

Lowercase letters indicate significant differences between N addition treatments under IM(P < 0.05). Capital letters indicate significant differences between N addition treatments under CM (P < 0.05).Asterisks indicate significant differences between CM and IM at identical N addition rates (*P < 0.05, **P < 0.01). Values represent means of three replicates, and error bars indicate standard errors.

In the IM plots, the maximum soil n class="Disease">MBC coclass="Chemical">nteclass="Chemical">nt (2237.3 mg∙kg−1) was observed iclass="Chemical">n the N30 treatmeclass="Chemical">nt aclass="Chemical">nd was 105% greater thaclass="Chemical">n that observed iclass="Chemical">n the coclass="Chemical">ntrol. Iclass="Chemical">n the class="Chemical">n class="Disease">CM plots, the maximum soil MBC content (1583.0 mg∙kg−1) was observed in the N60 treatment and was 79.6% greater than that in the control. While the soil MBC contents of the N addition treatments differed significantly from each other, all were significantly greater than that of the control treatment.

Soil microbial community composition

More than 35,000 valid reads were obtained for each replicate via a sequence optimization process and quality filtering. The median sequence length of each read was 253 bp. A total of 68,531 operational taxonomic units (OTUs) representing 41 phyla and 266 genera were detected using 97% identity as the cutoff. The dominant phyla were Proteobacteria (34.9%), Acidobacteria (27.7%), Verrucomicrobia (8.5%), and Actinobacteria (7%) (Fig. 2). In addition, 22 phyla were detected at relatively low abundances and accounted for less than 1% of the observed phyla. The dominant genera were DA101 (4.3%), Rhodoplanes (1.8%), Candidatus Koribacter (1.7%), Candidatus Solibacter (1.4%), and Rhodanobacter (1.4%).

Figure 2

Taxonomic distribution of soil microbial communities in the Moso bamboo plantation under all treatments.

Numbers indicate the ten most prevalent phyla: (1) Proteobacteria, (2) Acidobacteria, (3) Verrucomicrobia, (4) Actinobacteria, (5) Chloroflexi, (6) Planctomycetes, (7) Gemmatimonadetes, (8) AD3, (9) Crenarchaeota, and (10) Others, including sequences that could not be classified into any known group.

Although fewer OTUs were detected in the n class="Disease">CM plots (54,439) thaclass="Chemical">n iclass="Chemical">n the IM plots (58,187), a greater class="Chemical">number of phyla were observed uclass="Chemical">nder the class="Chemical">n class="Disease">CM treatments (39) than under the IM treatments (36). Under both CM and IM treatments, Proteobacteria, Acidobacteria, Verrucomicrobia, Actinobacteria, Chloroflexi, Planctomycetes, Gemmatimonadetes, AD3, Crenarchaeota, and Bacteroidetes contributed to a large proportion of phyla (Fig. 3). The relative abundances of Actinobacteria and Proteobacteria were higher under IM than under CM, whereas Acidobacteria, Crenarchaeota, and Verrucomicrobia exhibited the opposite trend.

Figure 3

Comparison of the soil bacterial communities at the phylum level in the Moso bamboo plantation under all treatments.

The relative abundances of the dominant bacterial groups in the soil differed depending on the nitrogen deposition and management practices. Relative abundances are based on the proportional frequencies of DNA sequences that could be classified. CN0: conventional management with no added N; CN30: conventional management with 30 kg N∙ha−1∙yr−1; CN60: conventional management with 60 kg N∙ha−1∙yr−1; CN90: conventional management with 90 kg N∙ha−1∙yr−1; IN0: intensive management with no added N; IN30: intensive management with 30 kg N∙ha−1∙yr−1; IN60: intensive management with 60 kg N∙ha−1∙yr−1; IN90: intensive management with 90 kg N∙ha−1∙yr−1.

No significant differences were observed in the number of OTUs among the N addition treatments. Under n class="Disease">CM, N additioclass="Chemical">n decreased the relative abuclass="Chemical">ndaclass="Chemical">nces of Proteobacteria aclass="Chemical">nd Bacteroidetes but iclass="Chemical">ncreased the relative abuclass="Chemical">ndaclass="Chemical">nce of Acidobacteria. Uclass="Chemical">nder IM, N additioclass="Chemical">n iclass="Chemical">ncreased the relative abuclass="Chemical">ndaclass="Chemical">nces of Creclass="Chemical">narchaeota aclass="Chemical">nd Acticlass="Chemical">nobacteria but decreased the relative abuclass="Chemical">ndaclass="Chemical">nce of Proteobacteria, with the exceptioclass="Chemical">n of the N90 (90 kg N∙ha−1∙yr−1) treatmeclass="Chemical">nt.

Soil microbial community and diversity

The Chao1 index, which reflects the species richness of a community, was significantly greater under n class="Disease">CM thaclass="Chemical">n uclass="Chemical">nder IM treatmeclass="Chemical">nt, regardless of N treatmeclass="Chemical">nt. Iclass="Chemical">ndex values were sigclass="Chemical">nificaclass="Chemical">ntly lower iclass="Chemical">n treatmeclass="Chemical">nts with added N compared to that of the coclass="Chemical">ntrol uclass="Chemical">nder both class="Chemical">n class="Disease">CM and IM (Fig. 4).

Figure 4

Chao 1 index of soil microbial diversity under conventional management (CM) or intensive management (IM) and three N addition treatments.

N30: 30 kg N∙ha−1∙yr−1; N60: 60 kg N∙ha−1∙yr−1; N90: 90 kg N∙ha−1∙yr−1. Lowercase letters indicate significant differences between N addition treatments under IM(P < 0.05). Capital letters indicate significant differences between N addition treatments under CM(P < 0.05).Asterisks indicate significant differences between CM and IM at the same N addition rate (*P < 0.05, **P < 0.01). Values represent means of three replicates, and error bars indicate standard errors.

Cluster analysis of the soil microbial communities divided the bacterial communities into two groups (Fig. 5). All samples from the n class="Disease">CM treatmeclass="Chemical">nts clustered iclass="Chemical">n group 1, aclass="Chemical">nd all samples from the IM treatmeclass="Chemical">nts clustered iclass="Chemical">n group II, regardless of N additioclass="Chemical">n. This groupiclass="Chemical">ng was further coclass="Chemical">nfirmed by priclass="Chemical">ncipal coordiclass="Chemical">nate aclass="Chemical">nalysis (PCoA) (Fig. 6). Differeclass="Chemical">nces iclass="Chemical">n microbial commuclass="Chemical">nity structure were primarily due to a combiclass="Chemical">natioclass="Chemical">n of class="Chemical">n class="Chemical">nitrogen deposition and management practices (57.73%), with management practices alone accounting for 36.26% of the variation and nitrogen addition accounting for 21.47%.

Figure 5

UniFrac UPGMA (unweighted pair group method with arithmetic mean) cluster analysis of microbial communities in soil samples from Moso bamboo plantation under different nitrogen deposition and management practice treatments.

Figure was constructed based on Illumina sequencing data. CN0: conventional management with no added N; CN30: conventional management with 30 kg N∙ha−1∙yr−1; CN60: conventional management with 60 kg N∙ha−1∙yr−1; CN90: conventional management with 90 kg N∙ha−1∙yr−1; IN0: intensive management with no added N; IN30: intensive management with 30 kg N∙ha−1∙yr−1; IN60: intensive management with 60 kg N∙ha−1∙yr−1; IN90: intensive management with 90 kg N∙ha−1∙yr−1.

Figure 6

Principal coordinate analysis (PCoA) plot based on 16S rRNA gene sequencing of 24 samples.

Scatter plot shows principal coordinate 1 (PC1) versus principal coordinate 2 (PC2). Percentages shown are percentages of variation explained by the components. CN0: conventional management with no added N; CN30: conventional management with 30 kg N∙ha−1∙yr−1; CN60: conventional management with 60 kg N∙ha−1∙yr−1; CN90: conventional management with 90 kg N∙ha−1∙yr−1; IN0: intensive management with no added N; IN30: intensive management with 30 kg N∙ha−1∙yr−1; IN60: intensive management with 60 kg N∙ha−1∙yr−1; IN90: intensive management with 90 kg N∙ha−1∙yr−1.

Relationship between microbial diversity and soil properties

The relationships between the α-diversity of soil microbes and soil properties were analyzed (Table 1). Chao1 and Shannon index values were significantly positively correlated with soil pH, and the Chao1 index was significantly negatively correlated with soil n class="Chemical">NO3− aclass="Chemical">nd class="Chemical">n class="Chemical">NH4+ concentrations. OTU number was also significantly positively correlated with pH but was negatively correlated with total N (TN), NO3−, and NH4+ concentrations. MBC was significantly negatively correlated with pH. Canonical correspondence analysis (CCA) demonstrated a correlation between primary soil properties and the microbial community; this showed that soil pH, C/N ratio, and available phosphorus (AP) and TN concentrations had the biggest impacts on the microbial community (Fig. 7).

Table 1

Pearson’s correlation coefficients between soil properties and soil microbial community indicators.

	OTUs	Chao1	Shannon	MBC
pH	0.844**	0.785**	0.588**	−0.553**
SOC	−0.222	−0.219	−0.172	0.169
TN	−0.475*	−0.403	−0.299	−0.079
C/N	0.348	0.262	0.14	0.095
AP	−0.163	−0.192	−0.164	−0.018
NO3−	−0.644**	−0.564**	−0.404	−0.018
NH4+	−0.661**	−0.502*	−0.392	−0.1

OTUs: operational taxonomic units (97% identity); MBC: microbial biomass carbon; SOC: soil organic carbon; TN: total nitrogen; AP: available phosphorus.

*P < 0.05, **P < 0.01.

Figure 7

Canonical correspondence analysis (CCA) of the relative abundances of dominant microbial and soil environmental factors.

Soil factors indicated in red text include pH, SOC (soil organic carbon), TN (total nitrogen), AP (available phosphorus), NH4 + , NO3−, and C/N ratio. Arrow lengths indicate the strength of the relationship between the soil property and the overall microbial community. CN0: conventional management with no added N; CN30: conventional management with 30 kg N∙ha−1∙yr−1; CN60: conventional management with 60 kg N∙ha−1∙yr−1; CN90: conventional management with 90 kg N∙ha−1∙yr−1; IN0: intensive management with no added N; IN30: intensive management with 30 kg N∙ha−1∙yr−1; IN60: intensive management with 60 kg N∙ha−1∙yr−1; IN90: intensive management with 90 kg N∙ha−1∙yr−1.

Discussion

Effect of management practices on soil microbial biomass and community diversity

In this study, we found significantly higher levels of n class="Disease">MBC iclass="Chemical">n the IM treatmeclass="Chemical">nts compared to the class="Chemical">n class="Disease">CM treatments without N addition (Figs 1 and 4), indicating that IM practices led to a significant increase in soil microbial biomass. This partially supports our first hypothesis, that IM practices increase soil MBC and diversity. IM practices such as fertilization may provide abundant nutrients for the growth of microorganisms and thus increase soil MBC3334. Fertilization can also increase the productivity of Moso bamboo35, resulting in increased litter and nutrient return, which may contribute to microbial growth. Similar results were observed by Yu et al.17.

The second part of this hypothesis, that IM increases microbial diversity, was not borne out by the results. In fact, we found that IM significantly reduced soil microbial diversity. Similarly, He et al.18 observed that after 15 years of IM, the abundance of N2-fixing bacteria decreased in n class="Species">Moso bamboo forests. A comparable result was also observed iclass="Chemical">n a class="Chemical">n class="Species">Castanea mollissima forest16. Some management practices, such as plowing and weeding, may reduce the diversity of the aboveground plant community32, alter soil conditions, and destroy original habitat, potentially creating an unfavorable environment for some microbes and reducing microbial diversity. In addition, the osmotic potential in the soil solution of the IM plots may have become toxic due to the introduction of additional ions via fertilizer36, further reducing microbial abundance.

Previous studies found that fertilization changed the microbial community’s structure and composition373839. Our results found that Actinobacteria were more prevalent with IM treatment, particularly with N addition. Actinobacteria can depolymerize the n class="Chemical">polyphenols iclass="Chemical">n litter iclass="Chemical">nto small, soluble molecules aclass="Chemical">nd thus play aclass="Chemical">n importaclass="Chemical">nt role iclass="Chemical">n the decompositioclass="Chemical">n of class="Chemical">n class="Chemical">lignin. The abundance of Actinobacteria therefore affects soil enzyme activity40. This result partially explains our previous finding that decomposition of leaf litter and lignin in Moso bamboo forests occurs more rapidly under IM than CM24. Acidobacteria, in contrast, were more prevalent in CM plots; these microbes are generally oligotrophic, consistent with their significantly lower abundances in the nutrient-rich rhizosphere and agricultural soils compared to that of bulk soil541. Crenarchaeota are a dominant group of microorganisms that govern NH4+ oxidation42. The decrease in the abundance of Crenarchaeota under IM that we observed may further affect the nitrifying process and N2O emissions in the soil of Moso bamboo forests.

Effect of nitrogen deposition on soil microbial biomass and community diversity

Our results demonstrated that a low amount of added N increased the soil n class="Disease">MBC, but wheclass="Chemical">n N additioclass="Chemical">n exceeded a certaiclass="Chemical">n threshold (30 kg N∙ha−1∙yr−1 for IM aclass="Chemical">nd 60 kg N∙ha−1∙yr−1 for class="Chemical">n class="Disease">CM), a sharp decrease in soil MBC occurred (Fig. 1). This partially supports our second hypothesis, that N addition increases soil MBC and microbial diversity, but it is clear that the relationship is dependent on the amount of N added. We previously observed a limited amount of N in the soil of Moso bamboo plantations under IM, even though fertilization provided some N43. It is also known that N deposition can increase soil microbial biomass in an N-limited region3144. Additional N input may directly increase available N and indirectly increase organic matter by enhancing plant productivity, thus promoting the growth of microorganisms26. However, excessive N input has negative effects on microbial activity. Many studies have reported that long-term nitrogen deposition reduced soil MBC2526. Nitrogen saturation induced by excess N input can decrease soil pH, leading to leaching of magnesium and calcium and mobilization of aluminum45. When this happens, microbes may become magnesium- or calcium-limited or suffer aluminum toxicity26, resulting in reduced microbial biomass. Our results indicate that the N input threshold for Moso bamboo plantations may be 60 kg N∙ha−1∙yr−1 (Fig. 1).

As for the effect of N addition on microbial diversity, we found that added N significantly decreased soil microbial community diversity, particularly under IM (Fig. 4). Similarly, Frey et al.29 observed that the diversity of the ectomycorrhizal fungal community was reduced under N addition (50 kg N∙ha−1∙yr−1) compared to that of control areas in the Harvard Forest. n class="Species">Compton et al.28 also observed that N additioclass="Chemical">n stroclass="Chemical">ngly iclass="Chemical">nflueclass="Chemical">nced the DNA profiles of the microbial commuclass="Chemical">nity. A meta-aclass="Chemical">nalysis revealed that the abuclass="Chemical">ndaclass="Chemical">nces of microbes decreased uclass="Chemical">nder N additioclass="Chemical">n aclass="Chemical">nd that these reductioclass="Chemical">ns were more evideclass="Chemical">nt wheclass="Chemical">n higher total amouclass="Chemical">nts of N were added aclass="Chemical">nd for loclass="Chemical">nger duratioclass="Chemical">ns26. Iclass="Chemical">n our study, the maiclass="Chemical">n reasoclass="Chemical">n for this reduced microbial diversity appears to be a reductioclass="Chemical">n iclass="Chemical">n soil pH iclass="Chemical">nduced by N fertilizatioclass="Chemical">n, particularly with IM treatmeclass="Chemical">nt. Previous studies have also fouclass="Chemical">nd that fertilizatioclass="Chemical">n treatmeclass="Chemical">nt reduces soil pH3946.

The N-induced reduction in microbial diversity may affect the structure of the microbial community47. n class="Disease">CM is class="Chemical">not coclass="Chemical">nducive to class="Chemical">nutrieclass="Chemical">nt accumulatioclass="Chemical">n, therefore, N additioclass="Chemical">n resulted iclass="Chemical">n a decrease iclass="Chemical">n the relative abuclass="Chemical">ndaclass="Chemical">nces of eutrophic taxa (iclass="Chemical">ncludiclass="Chemical">ng members of the Proteobacteria aclass="Chemical">nd Bacteroidetes) aclass="Chemical">nd aclass="Chemical">n iclass="Chemical">ncrease iclass="Chemical">n the abuclass="Chemical">ndaclass="Chemical">nces of oligotrophic taxa (maiclass="Chemical">nly Acidobacteria) (Fig. 3). Uclass="Chemical">nder IM, the relative abuclass="Chemical">ndaclass="Chemical">nces of Proteobacteria aclass="Chemical">nd Creclass="Chemical">narchaeota were higher iclass="Chemical">n the N90 treatmeclass="Chemical">nts thaclass="Chemical">n iclass="Chemical">n the coclass="Chemical">ntrol, whereas Acidobacteria aclass="Chemical">nd Verrucomicrobia exhibited the opposite treclass="Chemical">nd. Agaiclass="Chemical">n, this may be because the relative abuclass="Chemical">ndaclass="Chemical">nces of eutrophic taxa (iclass="Chemical">ncludiclass="Chemical">ng members of the Proteobacteria aclass="Chemical">nd Bacteroidetes) iclass="Chemical">ncreased, whereas those of oligotrophic taxa (maiclass="Chemical">nly Acidobacteria) decreased iclass="Chemical">n high-N plots48. Similar effects have beeclass="Chemical">n also observed iclass="Chemical">n previous studies2849. Acticlass="Chemical">nobacteria, aclass="Chemical">n importaclass="Chemical">nt microbe iclass="Chemical">n class="Chemical">n class="Chemical">lignin decomposition40, exhibited a greater relative abundance in the N30 treatment than in the N90 treatment under both IM and CM, thus providing a reasonable explanation for our previous observation that a low level of N addition facilitates the decomposition of Moso bamboo leaf litter, while a high level suppresses decomposition24.

Interactive effect of nitrogen deposition and management practices on soil microbial biomass and community diversity

Our third hypothesis was that the combination of management practices and N deposition would exert greater effects on soil microbial biomass and diversity than either practice independently. Two-way ANOVA demonstrated that both management practices and N deposition had significant positive effects on soil n class="Disease">MBC (Fig. 1 aclass="Chemical">nd Table 2) but sigclass="Chemical">nificaclass="Chemical">nt class="Chemical">negative effects oclass="Chemical">n soil microbial commuclass="Chemical">nity diversity (Fig. 4 aclass="Chemical">nd Table 3). These patterclass="Chemical">ns were true for the two variables, both aloclass="Chemical">ne aclass="Chemical">nd iclass="Chemical">n combiclass="Chemical">natioclass="Chemical">n, supporticlass="Chemical">ng our third hypothesis. PCoA aclass="Chemical">nd cluster aclass="Chemical">nalyses democlass="Chemical">nstrated that maclass="Chemical">nagemeclass="Chemical">nt practices had a greater impact oclass="Chemical">n soil microbial diversity thaclass="Chemical">n class="Chemical">n class="Chemical">nitrogen deposition (Figs 5 and 6). IM practices such as plowing, weeding, and fertilization alter the structural and physicochemical properties of soil more than the addition of nitrogen alone. Moreover, IM has been performed at the site for approximately 15 years, whereas the N addition experiment was only conducted for 2.5 years. The cumulative effects of IM practices on soil microbes were therefore greater than those of N addition owing to the short duration of the experiment.

Table 2

Two-way ANOVA indicating the effects of N addition and management type (intensive or conventional) on soil microbial biomass carbon.

Source of variation	SS	df	MS	F	P-value
N addition	2742792.2	3	914264.1	1574.6	<0.0001
Management type	861123.7	1	861123.7	1483.0	<0.0001
Interaction	1264365.6	3	421455.2	725.8	<0.0001
Within	9290.4	16	580.7
Total	4877572.0	23

SS: sum-of-squares; df: degrees of freedom; MS: mean square.

Table 3

Two-way ANOVA indicating the effects of N addition and management type (intensive or conventional) on the Chao 1 index of soil microbial community composition.

Source of variation	SS	df	MS	F	P-value
N addition	1295767.9	3	431922.6	9.4	0.0008
Management type	8005801.7	1	8005802	175.2	<0.0001
Interaction	840578.7	3	280192.9	6.1	0.006
Within	730924.5	16	45682.8
Total	10873072.9	23

SS: sum-of-squares; df: degrees of freedom; MS: mean square.

Previous research indicated that soil pH may strongly influence the soil microbial community composition and bacterial diversity. Soil pH is thought to be a good predictor of bacterial community composition2837395051525354. The results of this study demonstrated that bacterial diversity and community composition were mainly correlated with soil pH (Fig. 7 and Table 1). A strong negative relationship between bacterial diversity and soil pH was also observed by Lauber et al.55 and has been supported by other studies850. Our results also demonstrated that N, available P, and the C/N ratio are important parameters for shaping microbial community structure (Fig. 7). In accordance with this, previous studies observed that bacterial community structure was closely correlated with the n class="Chemical">NO3− coclass="Chemical">nceclass="Chemical">ntratioclass="Chemical">n395657. Additioclass="Chemical">nally, Zhao et al.58 observed that Betaproteobacteria aclass="Chemical">nd Deltaproteobacteria abuclass="Chemical">ndaclass="Chemical">nces were positively correlated with available soil P coclass="Chemical">nteclass="Chemical">nt. Chu et al.59 observed that the C/N ratio was sigclass="Chemical">nificaclass="Chemical">ntly correlated with the relative abuclass="Chemical">ndaclass="Chemical">nces of various domiclass="Chemical">naclass="Chemical">nt phyla iclass="Chemical">n class="Chemical">n class="Disease">Arctic tundra soils, consistent with the abundances and composition of soil bacterial communities in the Changbai Mountains60. These results indicate that there are significant correlations between the diversity and composition of microbial communities and soil properties. Therefore, management practices and N addition may affect microbial diversity through changes in soil parameters.

In the present study, measurement of the soil microbial biomass included both bacteria and fungi, whereas assessment of the V4 hypervariable region of the 16S rRNA gene only reflected changes in the bacterial community composition. A previous study revealed that N deposition also inhibits fungal community composition26. Therefore, future studies should explore the effects of N deposition and management practices on fungal diversity.

In conclusion, we found that while both intensive management and N addition increased soil microbial biomass and reduced microbial diversity on a n class="Species">Moso bamboo placlass="Chemical">ntatioclass="Chemical">n, maclass="Chemical">nagemeclass="Chemical">nt practices exerted a greater iclass="Chemical">nflueclass="Chemical">nce over regulatioclass="Chemical">n of these parameters thaclass="Chemical">n N depositioclass="Chemical">n. These results class="Chemical">not oclass="Chemical">nly have practical implicatioclass="Chemical">ns for the maclass="Chemical">nagemeclass="Chemical">nt of class="Chemical">n class="Species">Moso bamboo forests, but also serve to broaden our understanding of the effect of human interventions and N deposition on soil microbial communities.

Methods

Study site

The study site was located in Qingshan Town, Lin’an City (30°14′N, 119°42′E), Zhejiang Province, China. The area has a monsoonal subtropical climate with a mean annual precipitation of 1,420 mm and a mean annual temperature of 15.6 °C, ranging from 24 °C in July to 3 °C in January. The area receives an average of approximately 1,847 hours of sunshine per year and features an average of 230 frost-free days per year.

Experimental design and measurement

Detailed information on the experimental design can be found in Song et al.24. Briefly, n class="Disease">CM aclass="Chemical">nd IM class="Chemical">n class="Species">Moso bamboo plantations were established at the study site. The CM plantation was originally established in the late 1970 s, and the IM plantation was established based on the CM area in 2001. CM practices included only the selective and regular harvest of bamboo trunks and shoots, which also took place in the IM area24. The IM area included additional management practices, such as plowing, weeding with herbicides, and fertilization. Specifically, in September of each year, 450 kg∙ha−1 of compound fertilizer (15:6:20 N:P2O5:K2O) was manually and evenly scattered on the soil surface, followed by deep plowing to 0.3 m. There was a greater abundance and biomass of plant species in the CM area than in the IM area24.

For this study, we established 12 n class="Disease">CM plots aclass="Chemical">nd 12 IM plots iclass="Chemical">n November 2012. Each plot had aclass="Chemical">n area of 20 × 20 m aclass="Chemical">nd was surrouclass="Chemical">nded by a 20-mbuffer zoclass="Chemical">ne. Iclass="Chemical">n accordaclass="Chemical">nce with the local backgrouclass="Chemical">nd atmospheric N depositioclass="Chemical">n rate of 30–37 kg N∙ha−1∙yr−1 616263, N additioclass="Chemical">n treatmeclass="Chemical">nts iclass="Chemical">ncluded a low-N treatmeclass="Chemical">nt (30 kg N∙ha−1∙yr−1, N30), a medium-N treatmeclass="Chemical">nt (60 kg N∙ha−1∙yr−1, N60), aclass="Chemical">nd a high-N treatmeclass="Chemical">nt (90 kg N∙ha−1∙yr−1, N90). Three replicate plots per treatmeclass="Chemical">nt were raclass="Chemical">ndomly placed iclass="Chemical">n each maclass="Chemical">nagemeclass="Chemical">nt area. Begiclass="Chemical">nclass="Chemical">niclass="Chemical">ng iclass="Chemical">n Jaclass="Chemical">nuary 2013, followiclass="Chemical">ng a widely used method for simulaticlass="Chemical">ng N depositioclass="Chemical">n6465, class="Chemical">n class="Chemical">NH4NO3 was dissolved in 10 L water to the proper final N concentration and sprayed evenly on the forest floor with an electric sprayer at the beginning of every month. An equivalent amount of N-free water was sprayed on each control plot to control for the effect of the added water24.

Soil sampling

In July 2015, ten soil cores at a depth of 0–20 n class="Disease">cm were raclass="Chemical">ndomly collected from each plot aclass="Chemical">nd mixed together. Samples were sieved through 2-mm mesh to remove roots, placlass="Chemical">nt residues, aclass="Chemical">nd stoclass="Chemical">nes. A portioclass="Chemical">n of each soil sample was collected iclass="Chemical">n a 50-mL ceclass="Chemical">ntrifuge tube, placed iclass="Chemical">n a cooler, aclass="Chemical">nd traclass="Chemical">nsferred to the laboratory. The tubes were archived at −80 °C uclass="Chemical">ntil beiclass="Chemical">ng used for molecular aclass="Chemical">nalysis. The remaiclass="Chemical">niclass="Chemical">ng samples were used to measure the microbial biomass (from field-moist soil) aclass="Chemical">nd were theclass="Chemical">n air-dried to determiclass="Chemical">ne the physicochemical properties of the soil.

Analysis of soil microbial biomass and physiochemical properties

n class="Disease">MBC geclass="Chemical">nerally accouclass="Chemical">nts for 40–50% of the microbial biomass aclass="Chemical">nd 1–5% of soil class="Chemical">n class="Chemical">organic carbon, and is thus regarded as an important indicator of soil microbial biomass66. We therefore measured MBC to reflect the dynamics of soil microbial biomass. MBC was estimated using the chloroform fumigation-extraction method67. The pH of a soil-water (1:2.5 w/v) suspension was measured using a pH meter (FE20, Mettler Toledo, Switzerland) after shaking for 30 min. Soil organic carbon (SOC) and TN were determined using an elemental analyzer (Elementar Vario EL III, Germany). AP in the soil was extracted with sodium bicarbonate and determined using the molybdenum blue method68. The soil KCl-extractable NO3− and NH4+ were determined by extraction with 2 M KCl, followed by steam distillation and titration69.

DNA extraction and library construction

Microbial DNA was extracted from the soil samples (0.25 g of wet weight) using an Ezup Column Soil DNA Purification Kit (Sangon Biotech, Shanghai, China) according to the manufacturer’s instructions. Total DNA was evaluated on a 1.0% n class="Chemical">agarose gel, aclass="Chemical">nd the DNA coclass="Chemical">nceclass="Chemical">ntratioclass="Chemical">n aclass="Chemical">nd quality (A260/A280) of the extracts were estimated visually usiclass="Chemical">ng a Naclass="Chemical">noDrop ND-1000 UV-V spectrophotometer (Thermo Scieclass="Chemical">ntific, Rockwood, TN, USA).

16S rRNA gene amplification and sequencing

The V4 hypervariable region of the 16S rRNA gene was amplified from each soil sample using the PCR primers 515F (5′-GTGCCAGn class="Disease">CMGCCGCGGTAA-3′) aclass="Chemical">nd 806R (5′-GGACTACHVGGGTWTCTAAT-3′) aclass="Chemical">nd a sample taggiclass="Chemical">ng approach7071. PCRswere performed iclass="Chemical">n 30-μL reactioclass="Chemical">ns with 15 μL of Phusioclass="Chemical">n® High-Fidelity PCR Master Mix (New Eclass="Chemical">nglaclass="Chemical">nd BioLabs, Ipswich, MA, USA), 0.2 μM forward aclass="Chemical">nd reverse primers, aclass="Chemical">nd approximately 10 class="Chemical">ng of template DNA. Thermal cycliclass="Chemical">ng coclass="Chemical">nditioclass="Chemical">ns coclass="Chemical">nsisted of iclass="Chemical">nitial declass="Chemical">naturatioclass="Chemical">n at 98 °C for 1 miclass="Chemical">n; 30 cycles of declass="Chemical">naturatioclass="Chemical">n at 98 °C for 10 s, aclass="Chemical">nclass="Chemical">nealiclass="Chemical">ng at 50 °C for 30 s, aclass="Chemical">nd eloclass="Chemical">ngatioclass="Chemical">n at 72 °C for 30 s; aclass="Chemical">nd a ficlass="Chemical">nal exteclass="Chemical">nsioclass="Chemical">n at 72 °C for 5 miclass="Chemical">n. PCR products were purified with a Geclass="Chemical">ne JET Gel Extractioclass="Chemical">n Kit (Thermo Scieclass="Chemical">ntific) aclass="Chemical">nd combiclass="Chemical">ned iclass="Chemical">n equimolar ratios with the quaclass="Chemical">ntitative DNA biclass="Chemical">ndiclass="Chemical">ng method to create a DNA pool that was subsequeclass="Chemical">ntly used for sequeclass="Chemical">nciclass="Chemical">ng from the adaptor. High-throughput sequeclass="Chemical">nciclass="Chemical">ng of the 16S rRNA tag-eclass="Chemical">ncoded geclass="Chemical">ne was performed oclass="Chemical">n the Illumiclass="Chemical">na MiSeq platform at Novogeclass="Chemical">ne (Beijiclass="Chemical">ng, Chiclass="Chemical">na).

Bioinformatics analysis

Sequencing reads were assigned to each sample according to the unique barcode tags. Sequences were analyzed with the QIIME (Quantitative Insights Into Microbial Ecology) software package and UPARSE pipeline70. The reads were first filtered by QIIME quality filters. The default settings for Illumina processing in QIIME were used. Then, the UPARSE pipeline was used to discern OTUs with 97% identity. For each OTU, a representative sequence was selected and a taxonomic group was assigned using the Ribosomal Database Project (RDP) classifier. The species richness of each sample was estimated by rarefaction analysis; the Chao1 index of each library was determined as described previously72.

Statistical analysis

One-way analysis of variance (ANOVA) including post-hoc correction for multiple comparisons using the Bonferroni method was used to determine the statistical significance of differences in n class="Disease">MBC aclass="Chemical">nd the Chao1 iclass="Chemical">ndex for each sampliclass="Chemical">ng eveclass="Chemical">nt amoclass="Chemical">ng the four experimeclass="Chemical">ntal treatmeclass="Chemical">nts iclass="Chemical">n the two placlass="Chemical">ntatioclass="Chemical">n maclass="Chemical">nagemeclass="Chemical">nt schemes. Two-way ANOVA was performed to assess the combiclass="Chemical">ned effect of N depositioclass="Chemical">n aclass="Chemical">nd maclass="Chemical">nagemeclass="Chemical">nt practices oclass="Chemical">n microbial biomass aclass="Chemical">nd diversity. The Pearsoclass="Chemical">n correlatioclass="Chemical">n coefficieclass="Chemical">nt betweeclass="Chemical">n soil properties aclass="Chemical">nd α-diversity was also calculated. Aclass="Chemical">nalyses were coclass="Chemical">nducted usiclass="Chemical">ng SPSS (Statistical Package for the Social Scieclass="Chemical">nces) 18.0 for Wiclass="Chemical">ndows (SPSS Iclass="Chemical">nc., Chicago, Illiclass="Chemical">nois, USA).

QIIME was used to calculate the weighted UniFrac. PCoA and UPGMA (unweighted pair group method with arithmetic mean) clustering were conducted on the weighted UniFrac based on a protocol published previously73. In addition, a CCA was performed to identify the abiotic factors with the most impact on bacterial community composition, and these results were used to construct a soil property matrix for variation partitioning analysis in R v.2.8.1 with the vegan package.

Additional Information

How to cite this article: Li, Q. et al. class="Chemical">Nitrogen depositioclass="Chemical">n aclass="Chemical">nd maclass="Chemical">nagemeclass="Chemical">nt practices iclass="Chemical">ncrease soil microbial biomass class="Chemical">n class="Chemical">carbon but decrease diversity in Moso bamboo plantations. Sci. Rep. 6, 28235; doi: 10.1038/srep28235 (2016).