Effects of stand age and soil properties on soil bacterial and fungal community composition in Chinese pine plantations on the Loess Plateau

Abstract

The effects of Chinese pine (Pinus tabuliformis) on soil variables after afforestation have been established, but microbial community changes still need to be explored. Using high-throughput sequencing technology, we analyzed bacterial and fungal community composition and diversity in soils from three stands of different-aged, designated 12-year-old (PF1), 29-year-old (PF2), and 53-year-old (PF3), on a Chinese pine plantation and from a natural secondary forest (NSF) stand that was almost 80 years old. Abandoned farmland (BL) was also analyzed. Shannon index values of both bacterial and fungal community in PF1 were greater than those in PF2, PF3 and NSF. Proteobacteria had the lowest abundance in BL, and the abundance increased with stand age. The abundance of Actinobacteria was greater in BL and PF1 soils than those in other sites. Among fungal communities, the dominant taxa were Ascomycota in BL and PF1 and Basidiomycota in PF2, PF3 and NSF, which reflected the successional patterns of fungal communities during the development of Chinese pine plantations. Therefore, the diversity and dominant taxa of soil microbial community in stands 12 and 29 years of age appear to have undergone significant changes; afterward, the soil microbial community achieved a relatively stable state. Furthermore, the abundances of the most dominant bacterial and fungal communities correlated significantly with organic C, total N, C:N, available N, and available P, indicating the dependence of these microbes on soil nutrients. Overall, our findings suggest that the large changes in the soil microbial community structure of Chinese pine plantation forests may be attributed to the phyla present (e.g., Proteobacteria, Actinobacteria, Ascomycota and Basidiomycota) which were affected by soil carbon and nutrients in the Loess Plateau.

About the Authors:

Peng Dang

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing

Affiliations College of Forestry, Northwest A&F University, Yangling, China, State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, Yangling, China

Xuan Yu

Roles Conceptualization, Methodology, Writing – review & editing

Affiliation: College of Forestry, Northwest A&F University, Yangling, China

Hien Le

Roles Formal analysis

Affiliation: College of Forestry, Northwest A&F University, Yangling, China

Jinliang Liu

Roles Formal analysis

Affiliation: College of Forestry, Northwest A&F University, Yangling, China

Zhen Shen

Roles Investigation

Affiliation: College of Forestry, Northwest A&F University, Yangling, China

Zhong Zhao

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Writing – review & editing

* E-mail: [email protected]

Affiliations College of Forestry, Northwest A&F University, Yangling, China, State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, Yangling, China

ORCID http://orcid.org/0000-0002-7992-2512

Introduction

Soil bacterial and fungal communities play important roles in the decomposition of organic matter and provide available nutrients for plant growth [1, 2]. As such, these communities are dynamic components of terrestrial ecosystems and exhibit temporal and spatial variation [3, 4]. Previous studies have shown that vegetation changes with changes in land use [5, 6]. These land use changes directly affect litter, root and exudates that indirectly affect soil organic carbon (SOC) and other edaphic properties of the soil that eventually translate into alterations in ecological succession [7, 8]. Microbial communities are influenced by pH, SOC, and other soil factors [9, 10]. Additionally, soil organisms, nutrients and associated microbial communities are influenced by changes in land use or forest growth [5].

Microbial community structure and activity change during ecological restoration of abandoned agricultural soils [11, 12] and post-mining areas [13, 14] through plant secondary succession [15, 16]. Belowground properties, including plant diversity and richness, soil pH, total carbon, total nitrogen and other nutrient content, change simultaneously with changes in the aboveground structure. Furthermore, these indexes ultimately influence the structure and function of the soil microbial community [17–19]. In forest ecosystems, afforestation may alter the physicochemical and activity properties of the soil due to changes in the forest tree structure [20, 21]. Determining the soil microbial community dynamics that reflect changes in forest age and soil properties is important when evaluating the interplay between aboveground and belowground communities.

The Loess Plateau has an area of approximately 9.6×106 ha and lies in the upper-middle reaches of the Yellow River in China [22]. This region experiences serious soil erosion and drought. To accelerate ecological rehabilitation and improve ecological stability in this region, a project called “Grain for Green” has been implemented on the Loess Plateau. This project aims to improve the fragile ecosystem by converting farmland into forestland [23]. As the predominant pioneer tree species, and Chinese pine is widely planted for Loess Plateau afforestation due to its high tolerance to cold, drought, and poor soil quality [24, 25]. Previous studies of Chinese pine forests have focused primarily on the aboveground ecosystem, soil chemistry, enzyme activity and microorganisms across different stand ages or management strategies [26]. For example, a PCR-DGGE approach was used to characterize soil bacterial and fungal communities in the rhizosphere of 30-year-old Chinese pine stands on the Loess Plateau [27]. However, successional changes in soil bacterial and fungal communities, which are crucial for the assessment of ecological restoration in Chinese pine plantations, have never been studied.

The study aimed (1) to examine the bacterial and fungal community structure composition and diversity of a Chinese pine plantation in stands of different ages and (2) to reveal the relationship between soil physicochemical properties and microbial communities. We employed sequencing technology to identify and study microbial groups with relatively low abundance (<1%) [28], quantify and model successional dynamics, such as those associated with microbial community structure and diversity, and characterize relationship between microorganisms and soils [29, 30].

Materials and methods

Study area description

The study was conducted in the Chinese pine forests located in the Caijiachuan forest area of Huanglong County (35°28′–36°02′N, 109°38′–110°12′E), which is on the southern Loess Plateau in northwestern China (Fig 1). Permission was granted by the Huanglong Forestry Bureau to carry out our research. The area of the Loess Plateau is characterized by many hills and gullies, with an altitude of 1200 meters. The study area is in the transition zone between semi-humid and semi-arid. Its climate is warm and temperate with a mean annual temperature of 8.6°C. The annual precipitation is 611.8 mm, with the majority of rainfall falling between July and September. The frost-free season averages 158 days in length. The Loess Plateau receives an average of 2369.8 h of daylight annually. The soil type is Typic-Loessi Orthic Primosols [31].

[Figure omitted. See PDF.]

Fig 1. Map of the study sites located in the Caijiachuan forest area of Huanglong county on the Loess Plateau.

https://doi.org/10.1371/journal.pone.0186501.g001

Large areas of abandoned farmlands with low yields were converted into forest land to improve the fragile ecosystem on the Loess Plateau decades ago. As the dominant native tree species, Chinese pine was planted beginning in the 1960s. Before planting, the main crop species were maize (Zea mays subsp), millet (Setaria italica) and pea (Pisum sativum Linn), and the only moisture applied during the growing seasons was derived from the falling rain. In this study, we collected samples from three stands of different ages that were planted in 1962, 1986 and 2003. For comparison with the plantation samples, control samples were obtained from farmland left uncultured for three years and an 80-year-olds Chinese pine natural secondary forest.

Sampling

Sampling was conducted in August 2015. Stands of the following three ages from the Chinese pine plantation were examined: 12-years (PF1), 29-years (PF2) and 53-years (PF3). One abandoned farmland site and one natural secondary Chinese pine forest stand with an age of almost 80 years (NSF) were selected as controls. The basic characteristics of the different stand sites in the Chinese pine forest and abandoned farmland are described in Table 1.

[Figure omitted. See PDF.]

Table 1. The basic conditions of the forests and abandoned farmland sites.

https://doi.org/10.1371/journal.pone.0186501.t001

Three replicate plots (20 m × 20 m) with similar slopes, gradients and altitudes were established at each stand site. The distance between each sampling plot was approximately 100 m. For each triplicate plot, nine soil corns of topsoil (0~10 cm) were collected with an “S” shape after removing the litter layer and then mixed to form one sample. A stainless-steel corer with a 4.5-cm inner diameter was used to collect soil. Fifteen samples were collected in total (five stand sites × three plots). All samples were filtered to remove the roots, litter and stones and then transported to the laboratory on ice immediately. One aliquot of each sample was stored at -80°C for DNA extraction. The remaining aliquots were air-dried and used to analyze physicochemical properties.

Analyses of soil physicochemical properties

The soil pH was determined using an electrode pH meter (Sartorius PB-10, Germany) after shaking the soil water liquid suspension (1:5 wt/vol) [32]. The soil total carbon (SOC) was determined using an Elemental analyzer (Elementar, Germany). Total nitrogen (TN) content was assessed via the Kjeldahl method [33]. Ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3--N) levels were measured using an AA3 continuous flow analytical system with 1 M KCl extraction (AA3, Germany) [32]. Total phosphorus (TP) was measured by spectrophotometry after wet digestion with HClO4-H2SO4 [34]. Available phosphorus (AP) was measured using the colorimetric method with 0.5 M NaHCO3 extraction [35].

Soil DNA extraction

E.Z.N.A soil DNA kits (OMEGA, USA) were used to extract total soil microbial community DNA in triplicate per soil sample. Briefly, 0.5 g of soil was placed into tubes and processed with buffer containing detergent. Then, a heating-and-freezing process was conducted to precipitate soil contaminants such as protein and humic acid. Next, water was used to elute the pure DNA. Subsequently, the purified DNA was quantified using a spectrophotometer (Epoch, BioTek, USA).

PCR amplification and sequencing of 16S rRNA and 18S rRNA genes

The V3-V4 region of bacterial 16S rDNA genes was amplified by PCR using primers 338F (5’-barcode-ACTCCTACGGGAGGCAGCAG)-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’) [36]. The following PCR procedure was used: 95°C for 3min, followed by 27 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s with a final extension at 72°C for 10 min. The PCR reactions were performed in a 20 μL reaction volume (2 μL of 2.5 mM dNTPs, 4 μL of 5× FastPfu Buffer, 0.4 μL of FastPfu Polymerase, 0.8 μL of each primer at 5 μM, and 10 ng of template DNA). The primers for the fungal 18S rDNA genes were SSU0817F (5’-barcode-TTAGCATGGAATAATRRAATAGGA)-3’) and SSU1196R (5’-TCTGGACCTGGTGAGTTTCC-3’) [37]. The following PCR procedure was used: 95°C for 3 min, followed by 32 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s with a final extension at 72°C for 10 min. The PCR reactions were performed in a 20 μL reaction volume (4 μL of 5× FastPfu Buffer, 2 μL of 2.5 mM dNTPs, 0.4 μL of FastPfu Polymerase, 0.8 μL of each primer at 5 μM primer, and 10 ng of template DNA). Triplicate amplicons were extracted from 2% agarose gels and purified using a DNA Extraction Kit. Then, the DNA amplicons were combined in an equimolar ratio and analyzed via sequencing using the Illumina MiSeq platform.

Processing of sequencing data

Raw fastq files were demultiplexed and quality filtered using QIIME (version 1.17) as follows [38]. First, the 300 bp reads were truncated at any site that received an average quality score <20 over a 50 bp sliding window. Second, reads containing ambiguous characters and exact barcode matches were removed. Third, only sequences that had an overlap longer than 10 bp were assembled according to their overlapping sequence. Reads that could not be assembled were discarded. UPARSE (version 7.1 http://drive5.com/uparse/) was used to cluster the operational taxonomic units (OTUs) using a 97% similarity cutoff [39]. UCHIME was used to identify and remove chimeric sequences [40].

The taxonomy of each 16S rRNA gene sequence was analyzed with the RDP classifier (http://rdp.cme.msu.edu/) against the Silva reference database (http://www.arb-silva.de) using a confidence threshold of 70% (version 2013). The phylogenetic affiliation of each 18S rRNA gene sequence was analyzed with the RDP classifier (http://rdp.cme.msu.edu/) against the Silva (SSU117/119) 18S rRNA database using a confidence threshold of 70% [41].

Statistical analyses

All soil environmental factors (SOC, TN, TP, AP, NH4+-N, NO3--N and pH), microbial community diversity indexes (Shannon and Simpson), and the relative abundances of microbial communities were compared by performing one-way analysis (ANOVA). Multiple comparisons were carried out to compare the differences among sample sites according to the least significant difference (LSD). Principal coordinates analysis (PCoA) was used to assess the microbial community structures in the stands based on Bray-Curtis distances. Redundancy analysis (RDA) was applied to test the relationships between soil physicochemical properties and the major microbial groups. Analysis of variance was performed using SPSS version 11.5 software (SPSS Inc., Chicago, IL, USA). PCoA and RDA were performed using the R software package (V.3.4.0).

Results

Soil physicochemical properties

The SOC and TN concentrations and the C:N ratios in all forests samples increased significantly compared to those in the abandoned land samples(P<0.01). The concentrations of NO3--N in PF3 and NSF were significantly higher than those in BL, PF1 and PF2 (P<0.05). AP concentration and soil bulk density (BD) decreased with plantation age (Table 2).

[Figure omitted. See PDF.]

Table 2. Soil physicochemical characteristics of forest and abandoned land of different ages.

https://doi.org/10.1371/journal.pone.0186501.t002

Alpha and beta diversity of bacterial and fungal communities

The dataset comprised a total of 332,172 high-quality bacterial sequences after the elimination of chimeras, with an average of 22,145 sequences obtained from each soil sample. A total of 20,661 OTUs were identified at a 97% sequence similarity cutoff, with an average of 1,377 OTUs per sample. The rarefaction curves of all the soil samples are shown in S1A Fig. The Shannon index of PF1 was significantly higher than that of PF2, PF3, and NSF. However, the Simpson index was lower in PF1 than that in PF2, PF3, and NSF (Fig 2A).

[Figure omitted. See PDF.]

Fig 2.

Shannon index and Simpson index of bacterial (a) and fungal (b) communities in plantation stands of different ages and abandoned farmland.

https://doi.org/10.1371/journal.pone.0186501.g002

After the elimination of chimeras, a total of 477,045 fungal sequences were obtained from all soil samples. A total of 1,792 OTUs were identified at a 97% sequence similarity cutoff, with an average of 119 OTUs per sample. The rarefaction curves of all of the soil samples are shown in S1B Fig. The Shannon indexes of BL and PF1 were significantly higher than those of PF2, PF3, and NSF (P<0.05). However, the Simpson indexes of BL and PF1 were significantly lower than those of PF3 and NSF, respectively (P<0.05) (Fig 2B).

PCoA revealed a variation in the bacterial and fungal communities among all soil sample sites (Fig 3). For all the groups of bacterial communities, the first and second principal components explained 48.5% and 13.21% of the variance, respectively (Fig 3A). For fungal communities, the first and second principal components explained 66.23% and 13.09% of the variance, respectively (Fig 3B). Profiles of soil bacterial communities from the BL and PF1 sites tended to group together, while the PF3 and NSF sites also tended to group together. The fungal communities showed a trend similar to that of the bacterial communities.

[Figure omitted. See PDF.]

Fig 3.

PCoA for soil bacterial (a) and fungal (b) communities from forest sites (PF1, PF2, PF3, NSF) and abandoned land (BL).

https://doi.org/10.1371/journal.pone.0186501.g003

Composition of the bacterial and fungal communities

The dominant bacterial phyla were Proteobacteria (25.11%), Actinobacteria (23.77%), Acidobacteria (23.25%), Chloroflexi (14.32%), Gemmatimonadetes (4.06%), Nitrospirae (3.23%), Bacteroidetes (1.86%) and Planctomycetes (1.06%) (Fig 4). The relative abundances of the Proteobacteria subgroups (Alpha-, Beta-, Gamma-, and Delta-Proteobacteria) were 12.37%, 6.98%, 2.14%, and 4.14%, respectively (S2 Fig). The Alpha-proteobacteria sequences were dominated by Rhizobiales (6.87%), Sphingomonadales (2.40%), and Rhodospirillales (2.38%). MND1 (3.21%), Xanthomonadales (1.57%), and Syntrophobacterales (1.70%) were the most abundant groups within the Beta-proteobacteria, Gamma-proteobacteria, and Delta-proteobacteria, respectively (S3 Fig). The relative abundances of Proteobacteria and Bacteroidetes were significantly lower in the BL site than those in all other forest sites. Rhizobiales, Rhodospirillales, and Xanthomonadales all had the same variance. The relative abundance of Actinobacteria were significantly higher in the BL and PF1 sites than those in the older stands (PF2, PF3 and NSF) (Fig 4 and S4 Fig). The relative abundances of Nitrospirae in the NSF (5.05%) and PF3 (4.79%) sites were significantly higher than those in the PF2 (2.92%), PF1 (2.02%), and BL (1.38%) sites (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Mean relative abundances of dominant bacterial phyla across all different sites.

Different letters indicate signiFIcant differences (P<0.05) across the different sample sites. NS: not significant.

https://doi.org/10.1371/journal.pone.0186501.g004

Ascomycota and Basidiomycota sequences were the most abundant among the fungal communities in all soil samples, representing 51.29% and 45.06% of all sequences, respectively (Fig 5). The relative abundances of Ascomycota were significantly higher in the BL (80.19%) and PF1 (73.43%) sites than those in the PF2 (30.40%), PF3 (34.04%), and NSF (38.38%) sites. The relative abundances of Basidiomycota were significantly lower in the BL (10.94%) and PF1 (19.79%) sites than those in the PF2 (68.67%), PF3 (65.07%), and NSF (60.83%) sites. The remaining taxa included Glomeromycota (≈1.2%), Mucoromycotina (<1%), Chytridiomycota (<1%), and Zoopagales (<1%) (Fig 5). Of the Ascomycota, Sordariomycetes were the most abundant (16.58%), followed by Pezizomycetes (13.86%), Dothideomycetes (8.07%) and Eurotiomycetes (5.93%). Of the Basidiomycota, Agaricomycetes (42.21%) and Tremellomycetes (2.49%) were the most abundant groups (S5 Fig).

[Figure omitted. See PDF.]

Fig 5. Mean relative abundances of dominant fungal phyla across all different sites.

Different letters indicate significant differences (P<0.05) across the different sample sites.

https://doi.org/10.1371/journal.pone.0186501.g005

Correlation between microbial communities and soil physicochemical properties

The abundances of Delta- and Gamma-proteobacteria, Bacteroidetes, Planctomycetes, and Basidiomycota were correlated positively with the TOC, TN and C:N. In contrast, abundances of Actinobacteria and Ascomycota were correlated negatively with TOC, TN, and C: N. The abundances of Actinobacteria, Ascomycota and Glomeromycota were correlated positively with AP, while the abundances of Nitrospirae, Planctomycetes and Basidiomycota were correlated negatively with AP. The abundances of Delta- and Gamma-proteobacteria and Nitrospirae were correlated positively with NO3--N (S1 Table.). The first two axes explained 56.84% and 65.54% of the total variance for the bacterial and fungal communities, respectively. Based on these results, SOC, TN, C:N, AP and NO3--N showed close correlations with changes in the bacterial and fungal community composition (Figs 6 and 7). Soil bulk density (BD) was correlated negatively with fungal communities (Fig 7).

[Figure omitted. See PDF.]

Fig 6. RDA of abundant bacterial communities at the order level and soil chemical properties for soil samples from abandoned farmland and Chinese pine forests.

https://doi.org/10.1371/journal.pone.0186501.g006

[Figure omitted. See PDF.]

Fig 7. RDA of the dominant fungal communities at the order level and soil chemical properties for soil samples from abandoned land and Chinese pine forests.

https://doi.org/10.1371/journal.pone.0186501.g007

Discussion

Changes in microbial community diversity

The Shannon and Simpson indexes were determined to reveal differences in the complexity of soil microbial communities during the development of plantation stages. In our study, soil bacterial and fungal diversity was higher in PF1 than that in PF2, PF3 and NSF. However, the differences among PF2, PF3 and NSF were not significant. The soil microbial diversity of Chinese pine varied greatly among stands that were 12 and 29 years old during development of. However, after 29 years of age, the soil microbial diversity changed slightly. Plantation development affects soil microbes through several pathways [42–45]. For example, plantation trees cover directly affects the amount of light available and also affects the understory composition, which influences carbon allocation, the soil microbial community, nutrient cycling and water conservation of the soil. Hence, large differences in the plantation tree canopy between 12-year-old stands and other older stands may indirectly affect soil microbial diversity through the understory plants root system and litter substrates. Meanwhile, the major groups of Actinobacteria and Ascomycota demonstrated sharp decreases from the youngest 12-year-old plantation stands to the older stands undergoing succession. For example, most order-level fungal communities belonging to the Ascomycota phylum, significantly decreased after the 12-year-old plantation stage (S5 Fig). The decreases in these major groups, which were affected by soil nutrient and litter input, may contribute to the decline in microbial diversity observed with age in the study.

Clear differentiation in the PCoA plots illustrated that the bacterial and fungal communities at the BL and PF1 sites were similar to each other. Plantation stands in the later development stage (PF3) and NSF were similar to each other (Fig 3). The results indicate that afforestation plays an essential role in shaping microbial communities. With decades of plantation development, the soil microbial community structure of Chinese pine plantations has become similar to that of the NSF.

Changes in the composition of bacterial and fungal communities

The relative abundance of bacterial phyla in the present study was similar to that of other soil environments reported in previous studies [16, 46]. In our study, the relative abundance of Proteobacteria, in forest soils was significantly greater than that in BL and increased with forest age. Similar to the previous finding [47], Rhizobiales was the most abundant order of Proteobacteria among all sample sites in our study. Rhizobiales, which fixes nitrogen from the atmosphere, increased with forest age in our study site. Thus, N-cycling groups of bacteria (e.g., Rhizobiales) play vital roles in the process of soil restoration in arid regions [48]. As one of the major bacterial groups, the Actinobacteria play important roles in organic matter turnover and carbon cycling [49]. The negative association between Actinobacteria and Proteobacteria abundance in our study sites is similar to previous studies in grassland soils [50] and snowmelt in an alpine tundra soil when plant carbon inputs were greatest [51]. It is assumed that the Actinobacteria, compared with the less abundant Proteobacteria group, are spore-forming bacteria, meaning that they dominate under harsh and stressful soil conditions [52].

Basidiomycota and Ascomycota, which represent the main soil fungi decomposers [53, 54], were the two major fungal communities, accounting for more than 90% of total fungal phyla. The Agaricomycetes (Agaricales, Russulales), belonging to the phylum Basidiomycota, contain ectomycorrhizal fungi related to pine development. Agaricomycetes increased substantially after afforestation. The most common orders of mycorrhizal fungi, such as Russulales and Agaricales, are late-stage fungi that are always symbiotic with older and larger trees [55, 56]. This association may be related to the storage of lignocellulose organic matter in older forests, which is degraded by members of Basidiomycota [57–59]. Moreover, large overstory trees contribute to the construction of a sustainable environment for mycorrhizal fungi [60, 61]. This may explain why mycorrhizal fungi occupied a dominant position in the Chinese pine plantation stands 29-year and older in our study. Mycorrhizal fungi play important roles in the pine forest ecosystem by promoting the absorption and utilization of mineral elements in the roots of trees. After approximately 29 years of afforestation recovery, the Chinese pine plantation gradually formed a stable ecosystem, which provides a sustainable habitat environment for mycorrhizal fungi. However, the relative abundances of Ascomycota were greater in BL and PF1 than those in the other sites. Abundance was similar among the PF2, PF3 and NSF sites. Apparently, Ascomycota gradually withdraw from a dominant advantage when Chinese pine plantation growth progresses to later stages. We propose that shifts in the fungal community structure are mediated by the tree structure and plant cover after afforestation with Chinese pine trees, which is consistent with previous studies [62, 63].

Response of the microbial community to soil physicochemical properties

Previous studies showed that SOC played the most important role in microbial community structure, and N availability also had an impact on soil microbial communities [50, 64]. In our study, the relative abundances of Proteobacteria, Bacteroidetes and Nitrospirae were correlated positively with the SOC and TN concentrations. Similar results have been reported previously [65–67]. This finding indicates that these bacteria prefer enriched soil environments and belong to eutrophic groups. This finding supports previous reports that abundance of Proteobacteria and Bacteroidetes represent soils with abundant labile substrates, exhibiting copiotrophic attributes [16, 66]. As the dominant bacterial group of this study, Proteobacteria played a comparable functional role in long term plantation restoration in arid and semi-arid lands.

However, the relative abundance of Actinobacteria was correlated negatively with the SOC, TN, and C:N, similar to previous studies that Actinobacteria was linked with low carbon concentrations in temperate [49, 68, 69] and Arctic [70] soils with one exception [71]. In the exception [71], the abundance of Actinobacteria increased with nitrogen inputs, and suggested that the Actinobacteria may represent a copiotrophic group. Whether Actinobacteria populations belong to copiotrophic or oligotrophic groups remains an open question [65]. Actinobacteria may adopt an oligotrophic life strategy in the study sites environment. Meanwhile, abundance of Actinobacteria decreased markedly when croplands were converted to monoculture plantation [49]. In this latter case, the abundance of Actinobacteria correlated positively with the SOC mineralization rate. This study indicates that the conversion of cropland to forests decreases the amount of Actinobacteria abundance that are capable of decomposing more recalcitrant soil carbon. Actinobacteria may be well adapted to harsh environmental conditioning due to its metabolic versatility and widespread occurrence in pristine soil [72, 73].

Previous studies showed that NO3--N was a main soil-available element influenced by bacterial communities on the Tibetan Plateau [74]. Additionally, nitrogen addition and soil NO3--N content positively influenced the relative abundances of dominant bacterial groups [75, 76]. Zhang et al. [16] showed that NO3--N closely correlated with the relative abundances of bacterial communities during the natural succession of abandoned land on the loess plateau. Similarly, this result indicated that NO3--N was a dominant chemical factor that influences soil bacterial communities during plantation ecological processes on the Loess plateau. Thus, these dominant bacterial communities may be regulated by soil nutrients, particularly the soil carbon and nitrogen.

Basidiomycota were positively correlated with TOC, TN, C:N, and NO3--N. Values of TOC, TN and C:N ratio increased gradually with forest age. Basidiomycota play important roles in mediating the decomposition of low-quality lignified and aromatic substrates, which accumulate with increases in the tree cover and understory vegetation diversity of both coniferous and deciduous tree litter [77–79]. This finding explains why pine forest soils and hardwoods with high C:N ratios demonstrate a high prevalence of Basidiomycetes. The abundance of Ascomycota showed a positive correlation with AP. Meanwhile, Previous studies also considered the relationship between the relative abundance of Ascomycota and soil phosphorous content in soils. Therefore, Soil phosphorous is considered an important regulator of fungal communities in the soil [80].

Conclusion

Our study suggests that plantations influence soil properties and microbial communities. Microbial communities correlated with soil factors during the process of vegetation restoration. Obvious decreases in bacterial and fungal community diversity indexes occurred between 12 and 29 years after plantation. The relative abundance of Proteobacteria increased with age over 53 years of growth. However, the relative abundance of Actinobacteria decreased between 12 years and 29 years after plantation. Fungal communities transitioned from being Ascomycota dominant in the BL and 12-year-old plantation stand sites to Basidiomycota dominant after 29 years. These bacterial and fungal phyla reflected the patterns of belowground microorganisms during the process of plantation restoration in the Loess Plateau. Soil microbial community composition and diversity reached relatively stable states after 29 years of growth. Change in microbial communities correlated with changes in the availability of soil nutrients (SOC, TN, AP, C:N, and NO3--N).

Supporting information

[Figure omitted. See PDF.]

S1 Fig.

Rarefaction curves for soil bacterial (a) and fungal (b) communities at 97% sequence similarity in soil.

https://doi.org/10.1371/journal.pone.0186501.s001

(TIF)

S2 Fig. Mean relative abundances of Proteobacteria populations at the class level at each site.

https://doi.org/10.1371/journal.pone.0186501.s002

(TIF)

S3 Fig. Mean relative abundances of Proteobacteria populations at the order level at each site.

https://doi.org/10.1371/journal.pone.0186501.s003

(TIF)

S4 Fig. Mean relative abundances of Actinobacteria populations at the order level in each site.

https://doi.org/10.1371/journal.pone.0186501.s004

(TIF)

S5 Fig. Mean relative abundances of dominant fungal orders in each site.

https://doi.org/10.1371/journal.pone.0186501.s005

(TIF)

S1 Table. Correlations among dominant bacterial and fungal communities (phyla) and soil properties.

https://doi.org/10.1371/journal.pone.0186501.s006

(PDF)

Citation: Dang P, Yu X, Le H, Liu J, Shen Z, Zhao Z (2017) Effects of stand age and soil properties on soil bacterial and fungal community composition in Chinese pine plantations on the Loess Plateau. PLoS ONE12(10): e0186501. https://doi.org/10.1371/journal.pone.0186501

References

1. Allison V, Condron L, Peltzer D, Richardson S, Turner B. Changes in enzyme activities and soil microbial community composition along carbon and nutrient gradients at the Franz Josef chronosequence, New Zealand. Soil Biology and Biochemistry. 2007;39(7):1770–81.

2. Weber CF, Lockhart JS, Charaska E, Aho K, Lohse KA. Bacterial composition of soils in ponderosa pine and mixed conifer forests exposed to different wildfire burn severity. Soil Biology and Biochemistry. 2014;69:242–50.

3. Saetre P, Bååth E. Spatial variation and patterns of soil microbial community structure in a mixed spruce–birch stand. Soil Biology and Biochemistry. 2000;32(7):909–17.

4. Högberg MN, Briones MJ, Keel SG, Metcalfe DB, Campbell C, Midwood AJ, et al. Quantification of effects of season and nitrogen supply on tree below‐ground carbon transfer to ectomycorrhizal fungi and other soil organisms in a boreal pine forest. New Phytologist. 2010;187(2):485–93. pmid:20456043

5. Kasel S, Bennett LT, Tibbits J. Land use influences soil fungal community composition across central Victoria, south-eastern Australia. Soil Biology and Biochemistry. 2008;40(7):1724–32.

6. Lauber CL, Strickland MS, Bradford MA, Fierer N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biology and Biochemistry. 2008;40(9):2407–15.

7. Kasel S, Bennett LT. Land-use history, forest conversion, and soil organic carbon in pine plantations and native forests of south eastern Australia. Geoderma. 2007;137(3–4):401–13.

8. Malchair S, Carnol M. Microbial biomass and C and N transformations in forest floors under European beech, sessile oak, Norway spruce and Douglas-fir at four temperate forest sites. Soil Biology & Biochemistry. 2009;41(4):831–9.

9. Cao YS, Fu SL, Zou XM, Cao HL, Shao YH, Zhou LX. Soil microbial community composition under Eucalyptus plantations of different age in subtropical China. European Journal of Soil Biology. 2010;46(2):128–35.

10. Wen L, Lei PF, Xiang WH, Yan WD, Liu ST. Soil microbial biomass carbon and nitrogen in pure and mixed stands of Pinus massoniana and Cinnamomum camphora differing in stand age. Forest Ecology and Management. 2014;328:150–8.

11. Elhottová D, Szili-Kovács T, Tříska J. Soil microbial community of abandoned sand fields. Folia Microbiologica. 2002;47(4):435–40. pmid:12422524

12. Jia G-M, Zhang P-D, Wang G, Cao J, Han J-C, Huang Y-P. Relationship Between Microbial Community and Soil Properties During Natural Succession of Abandoned Agricultural Land. Pedosphere. 2010;20(3):352–60.

13. Hüttl RF, Weber E. Forest ecosystem development in post-mining landscapes: a case study of the Lusatian lignite district. Naturwissenschaften. 2001;88(8):322–9. pmid:11572012

14. Mummey DL, Stahl PD, Buyer JS. Soil microbiological properties 20 years after surface mine reclamation: spatial analysis of reclaimed and undisturbed sites. Soil Biology and Biochemistry. 2002;34(11):1717–25.

15. Kuramae EE, Gamper HA, Yergeau E, Piceno YM, Brodie EL, DeSantis TZ, et al. Microbial secondary succession in a chronosequence of chalk grasslands. ISME J. 2010;4(5):711–5. http://www.nature.com/ismej/journal/v4/n5/suppinfo/ismej201011s1.html. pmid:20164861

16. Zhang C, Liu G, Xue S, Wang G. Soil bacterial community dynamics reflect changes in plant community and soil properties during the secondary succession of abandoned farmland in the Loess Plateau. Soil Biology and Biochemistry. 2016;97:40–9.

17. Chen TH, Chiu CY, Tian GL. Seasonal dynamics of soil microbial biomass in coastal sand dune forest. Pedobiologia. 2005;49:645–53.

18. Tonon G, Boldreghini P, Gioacchini P. Seasonal changes in microbial nitrogen in an old broadleaf forest and in a neighbouring young plantation. Biology and Fertility of Soils. 2005;41(2):101–8.

19. Devi NB, Yadava PS. Seasonal dynamics in soil microbial biomass C, N and P in a mixed-oak forest ecosystem of Manipur, North-east India. Appl Soil Ecol. 2006;31(3):220–7.

20. Overby ST, Owen SM, Hart SC, Neary DG, Johnson NC. Soil microbial community resilience with tree thinning in a 40-year-old experimental ponderosa pine forest. Appl Soil Ecol. 2015;93:1–10.

21. Ren C, Zhao F, Kang D, Yang G, Han X, Tong X, et al. Linkages of C:N:P stoichiometry and bacterial community in soil following afforestation of former farmland. Forest Ecology and Management. 2016;376:59–66.

22. Zhang C, Liu G, Xue S, Wang G. Changes in rhizospheric microbial community structure and function during the natural recovery of abandoned cropland on the Loess Plateau, China. Ecological Engineering. 2015;75:161–71.

23. Chen LW, Wet ; Fu , Bojie ; Lu , Yihe . Soil and water conservation on the Loess Plateau in China: review and perspective. PROGRESS IN PHYSICAL GEOGRAPHY. 2007;31(4):389–403.

24. Zhang H, Tang M, Chen H, Tian Z, Xue Y, Feng Y. Communities of arbuscular mycorrhizal fungi and bacteria in the rhizosphere of Caragana korshinkii and Hippophae rhamnoides in Zhifanggou watershed. Plant and Soil. 2010;326(1–2):415.

25. Yuan BC, Yue DX. Soil Microbial and Enzymatic Activities Across a Chronosequence of Chinese Pine Plantation Development on the Loess Plateau of China. Pedosphere. 2012;22(1):1–12.

26. Zhang J-T. Succession analysis of plant communities in abandoned croplands in the eastern Loess Plateau of China. Journal of Arid Environments. 2005;63(2):458–74.

27. Yu HX, Wang CY, Tang M. Fungal and Bacterial Communities in the Rhizosphere of Pinus tabulaeformis Related to the Restoration of Plantations and Natural Secondary Forests in the Loess Plateau, Northwest China. Sci World J. 2013. pmid:24459438

28. Roesch LF, Fulthorpe RR, Riva A, Casella G, Hadwin AK, Kent AD, et al. Pyrosequencing enumerates and contrasts soil microbial diversity. The ISME journal. 2007;1(4):283–90. pmid:18043639

29. Lockwood JL, Powell RD, Nott MP, SL. P. Assembling ecological communities in time and space. Oikos. 1997;80:549–53.

30. Dini-Andreote F, Silva M, Triado-Margarit X, Casamayor EO, van Elsas JD, Salles JF. Dynamics of bacterial community succession in a salt marsh chronosequence: evidences for temporal niche partitioning. Isme Journal. 2014;8(10):1989–2001. pmid:24739625

31. An S, Mentler A, Mayer H, Blum WE. Soil aggregation, aggregate stability, organic carbon and nitrogen in different soil aggregate fractions under forest and shrub vegetation on the Loess Plateau, China. Catena. 2010;81(3):226–33.

32. Bao SD. Soil and Agricultural Chemistry Analysis. BeiJing: China Agriculture Press, Beijing, China (in Chinese); 2000.

33. Li H, Li J, He Y, Li S, Liang Z, Peng C, et al. Changes in carbon, nutrients and stoichiometric relations under different soil depths, plant tissues and ages in black locust plantations. Acta physiologiae plantarum. 2013;35(10):2951–64.

34. Parkinson J, Allen S. A wet oxidation procedure suitable for the determination of nitrogen and mineral nutrients in biological material. Communications in Soil Science & Plant Analysis. 1975;6(1):1–11.

35. Emteryd O. Chemical and physical analysis of inorganic nutrients in plant, soil, water and air: Swedish University of Agricultural Sciences. Department of Forest Site Research; 1989.

36. Dennis KL, Wang Y, Blatner NR, Wang S, Saadalla A, Trudeau E, et al. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by IL-10–producing T cells. Cancer research. 2013;73(19):5905–13. pmid:23955389

37. Borneman J, Hartin RJ. PCR primers that amplify fungal rRNA genes from environmental samples. Applied and environmental microbiology. 2000;66(10):4356–60. pmid:11010882

38. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nature methods. 2010;7(5):335–6. pmid:20383131

39. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature methods. 2013;10(10):996–8. pmid:23955772

40. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. pmid:21700674

41. Amato KR, Yeoman CJ, Kent A, Righini N, Carbonero F, Estrada A, et al. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. The ISME journal. 2013;7(7):1344–53. pmid:23486247

42. Mitchell RJ, Keith AM, Potts JM, Ross J, Reid E, Dawson LA. Overstory and understory vegetation interact to alter soil community composition and activity. Plant and soil. 2012;352(1–2):65–84.

43. Keith AM, Brooker RW, Osler GH, Chapman SJ, Burslem DF, van der Wal R. Strong impacts of belowground tree inputs on soil nematode trophic composition. Soil Biology and Biochemistry. 2009;41(6):1060–5.

44. Van Breemen N, Van Dijk H. Ecosystem effects of atmospheric deposition of nitrogen in the Netherlands. Environmental Pollution. 1988;54(3–4):249–74. pmid:15092521

45. Yin K, Zhang L, Chen D, Tian Y, Zhang F, Wen M, et al. Understory herb layer exerts strong controls on soil microbial communities in subtropical plantations. Scientific reports. 2016;6.

46. Peralta RM, Ahn C, Gillevet PM. Characterization of soil bacterial community structure and physicochemical properties in created and natural wetlands. Science of The Total Environment. 2013;443:725–32. pmid:23228718

47. Will C, Thurmer A, Wollherr A, Nacke H, Herold N, Schrumpf M, et al. Horizon-Specific Bacterial Community Composition of German Grassland Soils, as Revealed by Pyrosequencing-Based Analysis of 16S rRNA Genes. Applied and Environmental Microbiology. 2010;76(20):6751–9. pmid:20729324

48. Rojas C, Gutierrez RM, Bruns MA. Bacterial and eukaryal diversity in soils forming from acid mine drainage precipitates under reclaimed vegetation and biological crusts. Appl Soil Ecol. 2016;105:57–66.

49. Zheng J, Chen J, Pan G, Wang G, Liu X, Zhang X, et al. A long-term hybrid poplar plantation on cropland reduces soil organic carbon mineralization and shifts microbial community abundance and composition. Appl Soil Ecol. 2017;111:94–104.

50. Cheng J, Jing G, Wei L, Jing Z. Long-term grazing exclusion effects on vegetation characteristics, soil properties and bacterial communities in the semi-arid grasslands of China. Ecological Engineering. 2016;97:170–8.

51. Björk RG, Björkman MP, Andersson MX, Klemedtsson L. Temporal variation in soil microbial communities in Alpine tundra. Soil Biology and Biochemistry. 2008;40(1):266–8.

52. Teixeira LC, Peixoto RS, Cury JC, Sul WJ, Pellizari VH, Tiedje J, et al. Bacterial diversity in rhizosphere soil from Antarctic vascular plants of Admiralty Bay, maritime Antarctica. The ISME journal. 2010;4(8):989–1001. pmid:20357834

53. Vandenkoornhuyse P, Baldauf SL, Leyval C, Straczek J, Young JPW. Extensive fungal diversity in plant roots. Science. 2002;295(5562):2051–. pmid:11896270

54. Bastian F, Bouziri L, Nicolardot B, Ranjard L. Impact of wheat straw decomposition on successional patterns of soil microbial community structure. Soil Biology & Biochemistry. 2009;41(2):262–75.

55. Keizer PJ, Arnolds E. SUCCESSION OF ECTOMYCORRHIZAL FUNGI IN ROADSIDE VERGES PLANTED WITH COMMON OAK (QUERCUS-ROBUR L) IN DRENTHE, THE NETHERLANDS. Mycorrhiza. 1994;4(4):147–59.

56. Durall DM, Gamiet S, Simard SW, Kudrna L, Sakakibara SM. Effects of clearcut logging and tree species composition on the diversity and community composition of epigeous fruit bodies formed by ectomycorrhizal fungi. Canadian Journal of Botany-Revue Canadienne De Botanique. 2006;84(6):966–80.

57. Osono T, Takeda H. Fungal decomposition of Abies needle and Betula leaf litter. Mycologia. 2006;98(2):172–9. pmid:16894962

58. Yelle DJ, Ralph J, Lu FC, Hammel KE. Evidence for cleavage of lignin by a brown rot basidiomycete. Environmental Microbiology. 2008;10(7):1844–9. pmid:18363712

59. Ma AZ, Zhuang XL, Wu JM, Cui MM, Lv D, Liu CZ, et al. Ascomycota Members Dominate Fungal Communities during Straw Residue Decomposition in Arable Soil. Plos One. 2013;8(6). pmid:23840414

60. Jones MD, Durall DM, Cairney JWG. Ectomycorrhizal fungal communities in young forest stands regenerating after clearcut logging. New Phytologist. 2003;157(3):399–422.

61. Dove NC, Keeton WS. Structural Complexity Enhancement increases fungal species richness in northern hardwood forests. Fungal Ecology. 2015;13:181–92.

62. Ros M, Klammer S, Knapp B, Aichberger K, Insam H. Long-term effects of compost amendment of soil on functional and structural diversity and microbial activity. Soil Use and Management. 2006;22(2):209–18.

63. Macdonald CA, Thomas N, Robinson L, Tate KR, Ross DJ, Dando J, et al. Physiological, biochemical and molecular responses of the soil microbial community after afforestation of pastures with Pinus radiate. Soil Biology & Biochemistry. 2009;41:1642–51.

64. Zhang Y, Li Y, Wang L, Tang Y, Chen J, Hu Y, et al. Soil microbiological variability under different successional stages of the Chongming Dongtan wetland and its effect on soil organic carbon storage. Ecological engineering. 2013;52:308–15.

65. Fierer N, Bradford MA, Jackson RB. Toward an ecological classification of soil bacteria. Ecology. 2007;88(6):1354–64. pmid:17601128

66. Goldfarb KC, Karaoz U, Hanson CA, Santee CA, Bradford MA, Treseder KK, et al. Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Frontiers in microbiology. 2011;2.

67. Durrer A, Gumiere T, Taketani RG, da Costa DP, e Silva MdCP, Andreote FD. The drivers underlying biogeographical patterns of bacterial communities in soils under sugarcane cultivation. Appl Soil Ecol. 2016.

68. Griffiths B, Ritz K, Ebblewhite N, Dobson G. Soil microbial community structure: effects of substrate loading rates. Soil Biology and Biochemistry. 1998;31(1):145–53.

69. Fierer N, Schimel JP, Holden PA. Variations in microbial community composition through two soil depth profiles. Soil Biology and Biochemistry. 2003;35(1):167–76.

70. Bell TH, Yergeau E, Maynard C, Juck D, Whyte LG, Greer CW. Predictable bacterial composition and hydrocarbon degradation in Arctic soils following diesel and nutrient disturbance. The ISME journal. 2013;7(6):1200. pmid:23389106

71. Ramirez KS, Craine JM, Fierer N. Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Global Change Biology. 2012;18(6):1918–27.

72. Hamamura N, Olson SH, Ward DM, Inskeep WP. Microbial population dynamics associated with crude-oil biodegradation in diverse soils. Applied and Environmental Microbiology. 2006;72(9):6316–24. pmid:16957258

73. Quatrini P, Scaglione G, De Pasquale C, Riela S, Puglia A. Isolation of Gram‐positive n‐alkane degraders from a hydrocarbon‐contaminated Mediterranean shoreline. Journal of applied microbiology. 2008;104(1):251–9. pmid:17922832

74. Yuan Y, Si G, Wang J, Luo T, Zhang G. Bacterial community in alpine grasslands along an altitudinal gradient on the Tibetan Plateau. FEMS microbiology ecology. 2014;87(1):121–32. pmid:23991911

75. Yao MJ, Rui JP, Li JB, Dai YM, Bai YF, Hedenec P, et al. Rate-specific responses of prokaryotic diversity and structure to nitrogen deposition in the Leymus chinensis steppe. Soil Biology & Biochemistry. 2014;79:81–90.

76. Tian W, Wang L, Li Y, Zhuang K, Li G, Zhang J, et al. Responses of microbial activity, abundance, and community in wheat soil after three years of heavy fertilization with manure-based compost and inorganic nitrogen. Agriculture, Ecosystems & Environment. 2015;213:219–27.

77. Bardgett RD, McAlister E. The measurement of soil fungal: bacterial biomass ratios as an indicator of ecosystem self-regulation in temperate meadow grasslands. Biology and Fertility of Soils. 1999;29(3):282–90.

78. Six J, Frey SD, Thiet RK, Batten KM. Bacterial and fungal contributions to carbon sequestration in agroecosystems. Soil Science Society of America Journal. 2006;70(2):555–69.

79. Toljander JF, Eberhardt U, Toljander YK, Paul LR, Taylor AF. Species composition of an ectomycorrhizal fungal community along a local nutrient gradient in a boreal forest. New Phytologist. 2006;170(4):873–84. pmid:16684245

80. Lauber CL, Ramirez KS, Aanderud Z, Lennon J, Fierer N. Temporal variability in soil microbial communities across land-use types. Isme Journal. 2013;7(8):1641–50. pmid:23552625