Original Research ARTICLE
Soil Bacterial Community Structure and Co-occurrence Pattern during Vegetation Restoration in Karst Rocky Desertification Area
- 1Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Fuyang, China
- 2Observation and Research Station for Rock Desert Ecosystem, Puding, China
Vegetation restoration has been widely used in karst rocky desertification (KRD) areas of southwestern China, but the response of microbial community to revegetation has not been well characterized. We investigated the diversity, structure, and co-occurrence patterns of bacterial communities in soils of five vegetation types (grassland, shrubbery, secondary forest, pure plantation and mixed plantation) in KRD area using high-throughput sequencing of the 16S rRNA gene. Bray-Curtis dissimilarity analysis revealed that 15 bacterial community samples were clustered into five groups that corresponded very well to the five vegetation types. Shannon diversity was positively correlated with pH and Ca2+ content but negatively correlated with organic carbon, total nitrogen, and soil moisture. Redundancy analysis indicated that soil pH, Ca2+ content, organic carbon, total nitrogen, and soil moisture jointly influenced bacterial community structure. Co-occurrence network analysis revealed non-random assembly patterns of bacterial composition in the soils. Bryobacter, GR-WP33-30, and Rhizomicrobium were identified as keystone genera in co-occurrence network. These results indicate that diverse soil physicochemical properties and potential interactions among taxa during vegetation restoration may jointly affect the bacterial community structure in KRD regions.
Southwestern China is one of the three distribution centers of Karst in the world. From the 1950s to the 1980s, Karst regions in China experienced severe soil erosion, vegetation loss, and desertification caused by intensive anthropogenic activity, referred to as karst rocky desertification (KRD; Wang et al., 2004; Jiang et al., 2014). Revegetation has been widely implemented in ecological conservation and restoration efforts because it improves soil nutrient conditions and other environmental factors (Crouzeilles et al., 2016; Feng et al., 2016). In recent decades, vegetation restoration has been used to overcome the ecological degradation in these regions (Qi et al., 2013). For example, a number of ecological projects have been implemented in southwestern China, such as the Yangtze River Shelter-Forest Project, the Pearl River Shelter-Forest Project and diversified water and soil conservation projects, all of which have been based on vegetation restoration (Wang et al., 2004; Jiang et al., 2014).
Soil microbial communities play an essential role in shaping the aboveground biodiversity and functioning of terrestrial ecosystems by driving biogeochemical processes and mediating nutrient turnover (Doran and Zeiss, 2000; Bardgett and van der Putten, 2014). Thus, microbial community structure is an essential indicator of the health and sustainability of an ecosystem (Harris, 2003; Zak et al., 2003; Lewis et al., 2010). Co-occurrence network analysis provides insight into the structure of microbial communities and the interactions among microorganisms (Barberán et al., 2012). However, soil microbial communities and co-occurrence patterns to vegetation restoration are largely unknown. This information provides invaluable reference to the appropriate management and conservation of the degraded karst ecosystem.
Previous studies showed that microbial community structure can be affected by soil properties such as pH, C, N and moisture (Yergeau et al., 2007; Lauber et al., 2009; Chong et al., 2010; Chu et al., 2011; Yu et al., 2012; Bakker et al., 2013; Schlatter et al., 2015), and plant factor such as vegetation type (Oh et al., 2012; Shi et al., 2015). In this study, we elucidate the relationships among vegetation type, soil physicochemical properties, and microbial community structure in the KRD area using high-throughput sequencing of the 16S rRNA gene. The main objectives of this study are to test the following hypotheses: (i) the bacterial communities differ among the principal vegetation types of the KRD areas, and (ii) vegetation-associated soil properties play significant impacts on bacterial community structure. Considering the critical role of microbial interactions in determining the soil microbial communities, co-occurrence patterns among bacterial communities were also explored in the KRD soils.
Materials and Methods
Study Area and Sampling
A field experiment was carried out in five types of vegetation at the San-do-qing Forestry Station (25°02′30″–25°58′22″N, E103°58′37″–104°49′48″), located in Fuyuan County, Yunnan Province. Fuyuan County is on the eastern Yunnan Karst plateau. This area has a northern subtropical monsoon climate. The mean annual temperature is 13.8°C, with a mean minimum temperature of 5.7°C in January and a mean maximum temperature of 19.8°C in July (Li et al., 2014). The main type of soil is red limestone. The KRD area covers 601.018 km2, accounting for 29.33% of the Karst region (Li et al., 2016). The five typical types of vegetation in the area are grassland, shrubbery, secondary forest, pure plantation, and mixed plantation. The grassland, shrubbery, and secondary forest are representative of three successional stages, having been abandoned as farmland 3, 12, and 45 years prior, respectively. The pure plantation and mixed plantation were afforested from grassland 25 years before sampling. The grassland is covered mainly by Chrysopogon orientalis, Fragaria vesca, and Dicranopteris dichotoma. The shrubbery is dominated by Pyracantha fortuneana, Corylus yunnanensis, and Myrica nana. The secondary forest is dominated by Quercus aquifolioides and Lithocarpus dealbatus. The pure plantation consists of Pinus armandi, whereas the mixed plantation is dominated by Pinus yunnanensis and Alnus ferdinandi-coburgii. The understory of all three types of forest is dominated by Hypericum monogynum and Myrica nana.
Three 20 × 20 m plots were established at each site in August 2015. The minimum distance between plots was 500 m to avoid pseudoreplication. Soil samples were collected from the surface soil (0–10 cm) with a stainless steel cutting ring (5 cm in diameter). Six soil cores were selected in an S-shaped pattern from each plot and mixed to form one replicate. All soils were transported to the laboratory immediately after collection in sterile plastic bags on dry ice and divided into two portions. One subset was stored at -80°C for DNA analysis, and the other subset was air dried for physicochemical analysis.
Soil Physicochemical and Biological Parameters
Soil organic carbon (SOC) content was determined using potassium dichromate oxidation (Nelson and Sommers, 1982). Total nitrogen (TN) content was estimated with a TOC analyzer (Multi N/C 3100 TOC, Analytik, Jena, Germany). Soil moisture (SM) was analyzed by weighing the soil and calculating the mass lost after oven drying at 105°C until weight was stable (24 h). Soil pH was determined with a soil-to-water ratio of 1:2.5 (w/v) using a pH meter (FE20, Mettler-Toledo Instruments, China). Ca2+ content was measured using an atomic absorption spectrophotometer (ICE 3500, Thermo Scientific, United States).
DNA Extraction and Purification
DNA was extracted directly from the soil samples using the Power Soil Extraction Kit (Mo Bio Laboratories, San Diego, CA, United States) according to the manufacturer’s instructions. The concentration and purity of the extracted DNA were measured using a Nanodrop 2000 spectrometer (Thermo Fisher Scientific, Wilmington, DE, United States). The soil DNA was stored at -20°C until use.
The V4–V5 variable regions of the bacterial 16S ribosomal RNA gene were amplified via polymerase chain reaction (PCR) using two universal eubacterial primer pairs, 515F (5′-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′; Xiong et al., 2012). The forward and reverse primers were tagged with adapter, pad, and linker sequences. Barcode sequences (10-mer) were added to the reverse primer to pool multiple samples into one run for sequencing. All primers were synthesized by Invitrogen Life Technologies (Shanghai, China). PCR amplification was conducted using TransGen AP221-02: TransStart Fastpfu DNA Polymerase (TransGen Biotech, Beijing, China) and performed in a GeneAmp 9700 thermal cycler (Applied Biosystems, Foster City, CA, United States). The reaction mixture included 4 μL 5× FastPfu buffer, 2 μL 2.5 mM dNTPs, 0.2 μL BSA, 0.8 μL each primer (5 μM), 10 ng template DNA, and H2O to a final volume of 20 μL. Thermal cycling conditions were as follows: 95°C for 3 min 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. All samples were amplified in triplicate. PCR amplification was detected using 2% agarose gel electrophoresis. The triplicate amplification products were pooled and purified using the AxyPrep DNA Gel Extraction Kit (AXYGEN, Union City, CA, United States) and then quantified using QuantiFluorTM-ST (Promega, United States).
Bacterial 16S rRNA Gene Sequencing and Processing
The V4–V5 region (515F–907R) of the bacterial 16S rRNA gene was sequenced on the Illumina Miseq PE 250 platform. Bacterial raw reads were deposited in NCBI Sequence Read Archive (SRA) under the submission ID SUB2918840 and BioProject ID PRJNA397824. Processing of the 16S rRNA–derived sequence inventories was performed using QIIME (quantitative insights into microbial ecology; Caporaso et al., 2010). Briefly, partial 16S rRNA bacterial sequences were filtered using Mothur version 1.22.2 (Schloss et al., 2009) with the inclusion criteria of mean quality score ≥20 and length ≥ 250bp. Sequences were assigned to samples by exact matches of 10 bp barcodes. Then the Uchime algorithm was used to detect chimeric sequences from a chimera-free reference database (Edgar et al., 2011) via the Usearch tool. All chimeras were removed prior to further analysis. Operational taxonomic units (OTUs) were clustered at the 97% similarity level using UPARSE version 7.1 (Edgar, 2010)1. Final OTUs were generated based on the clustering results, and taxonomic assignment was performed with the RDP 16S Classifier2 (Wang et al., 2007).
Statistical analysis of OTU richness via Good’s coverage, Chao1, and Shannon’s index was performed with Mothur (version 1.22.2; Schloss et al., 2009). One-way analysis of variance (ANOVA) followed by Duncan’s multiple range test (DMRT) was performed to assess the significance of the effects of vegetation type on soil properties and microbial community composition and diversity using SPSS version 17.0 (SPSS Inc., Chicago, IL, United States). Bray-Curtis dissimilarity values were calculated to reveal the relationships among samples based on bacterial community composition. Shared and unique OTUs among the five vegetation types were used to generate a Venn diagram. The 50 most abundant OTUs among the five vegetation types were analyzed using the hierarchical clustering software Cluster version 3.03 and visualized using Java TreeView version 18.104.22.168 Redundancy analysis (RDA) was performed with Canoco (version 4.5 for Windows; Ithaca, NY, United States) to determine which environmental variables best explained the assemblage’s variability. Forward selection was based on Monte Carlo permutation tests (permutations = 999). The ordination in the x- and y-axes and the length of the corresponding arrows indicated the importance of each physicochemical factor in explaining the taxon distribution across communities. The co-occurrence of OTUs in microbial communities across the five vegetation types was analyzed. To reduce network complexity and facilitate the identification of the core soil community, we selected OTUs with more than five sequences for further analysis (Barberán et al., 2012). Significant Spearman correlations (p < 0.01) were noted, and visualization of the co-occurrence network was conducted using the Fruchtermann-Feingold layout of the interactive platform Gephi version 0.9.0. Possible keystone genera were those that demonstrated high betweenness centrality values (Vick-Majors et al., 2014). The modular structure of the community was evaluated via the modularity index (Lambiotte et al., 2015).
Soil Physicochemical Properties Associated with the Five Vegetation Types
Soil pH decreased as the vegetation developed from grassland (5.6) and shrubbery (5.6) to secondary forest (5.27), with pure plantation (4.69) exhibiting the lowest value (Table 1). In addition, no significant difference in Ca2+ concentration was observed among vegetation types except for pure plantation, which had the lowest value (0.37). In general, TN, SOC, and SM increased significantly (p < 0.05) as the vegetation changed from grassland to forest. Furthermore, carbon to nitrogen C:N ratios were comparable among the natural succession vegetation types (12.99–14.61), although pure plantation and mixed plantation showed higher values (14.9 and 15.79, respectively).
Distribution of Taxa and Phylotypes
We obtained a total of 1,182,391 high-quality bacterial V4–V5 Illumina sequences and 3778 OTUs (at 3% evolutionary distance) after applying all quality filters. Almost all sequences (99.95%) were between 350 and 400 bp, with an average read length of 396.21 bp. The number of sequences obtained from each sample ranged from 59,788 to 108,809. Good’s coverage values ranged from 98.9 to 99.4% (Table 2), with the number of OTUs increasing sharply before reaching a plateau, which indicates that the number of bacterial sequences obtained represented the bacterial communities well, as the rarefaction curves tended toward saturation (Supplementary Figure S1). The bacteria were from 32 phyla, 76 classes, 170 orders, 312 families, and 506 genera. The dominant phyla (except for Proteobacteria, which were characterized at the class level) across all samples were Acidobacteria (21.73–57.08%), Actinobacteria (2.21–22.89%), Alphaprotebacteria (9.95–16.8%), Chloroflexi (0.94–10.26%), Planctomycetes (3.59–9.79%), Deltaproteobacteria (4.21–8.64%), Gammaproteobacteria (3.34–7.34%), and Betaproteobacteria (1.42–3.97%) and to a lesser degree Gemmatimonadetes (0.063–6.02%), Bacteroidetes (0.88–3.38%), Armatimonadetes (0.26–1.81%), Nitrospirae (0–1.78%), and Latescibacteria (0–1.73%; >1%), which together accounted for more than 95% of bacterial sequences from each of the vegetation types (Supplementary Figure S2).
TABLE 2. Characteristics of soil bacteria richness and diversity indices under different vegetation types.
Bacterial Diversity and Differences in Community Structure among the Vegetation Types
According to OTU diversity estimated by Shannon’s index, the greatest bacterial diversity was in grassland and shrubbery soils (average = 6.21 and 6.18, respectively) followed by secondary forest and mixed plantation (average = 5.56 and 5.58, respectively), whereas pure plantation showed the lowest bacterial diversity (average = 4.72; Table 2). These results indicate that vegetation restoration plays an important role in determining soil bacterial diversity.
The relative abundance of each bacterial taxonomic group varied among the five vegetation types (Supplementary Figure S2). It is remarkable that the relative abundance of bacterial phyla associated with grassland and shrub differed significantly from the other three vegetation types. For example, the Acidobacteria and Gammaproteobacteria phylotypes were less abundant in grassland (33.95 and 3.49%, respectively) and shrubbery (21.80 and 3.56%) than in the secondary forest (49.53 and 5.57%), mixed plantation (50.20 and 5.55%), and pure plantation (53.62 and 6.82%) sites. Conversely, some taxa decreased markedly in relative abundance from grassland and shrubbery to the three forest types, including Actinobacteria, Betaproteobacteria, Gemmatimonadetes, and Nitrospirae.
The Impact of Vegetation-Associated Soil Characteristics on Bacterial Community Composition and Diversity
Bray-Curtis dissimilarity analysis revealed that the 15 bacterial community samples clustered into five groups that corresponded very well to the five vegetation types (Figure 1). Clustering indicated that the grassland and shrubbery sites were closely related and secondary forest and mixed plantation shared a close relationship. A Venn diagram demonstrated that OTUs differed among the five vegetation types (Figure 2). The number of site-specific OTUs ranged from 56 (pure plantation) to 321 (shrubbery). In addition, a total of 624 OTUs were shared among all five vegetation types; these were defined as generalists. Generalist OTUs were composed of a number of bacterial groups, including Acidobacteria, Proteobacteria, Actinobacteria, and Bacteroidetes.
FIGURE 1. Clustering analysis of bacterial communities under five vegetation types based on OTU abundance-based Bray-Curtis similarity coefficients.
FIGURE 2. Venn diagram of exclusive and shared bacterial OTUs (at the 3% evolutionary distance) under five vegetation types.
To investigate the differences in the soil bacterial communities of the five vegetation types, we used heatmap analysis of the 50 most abundant OTUs, which highlighted their relative distributions and abundances (Supplementary Table S2). As shown in the heatmap (Figure 3), the abundance of these 50 dominant OTUs differed among the five vegetation types. The dominant OTUs in each vegetation type were also different. For example, the mixed plantation site was dominated by OTU2328 and OTU2824, whereas pure plantation was enriched by OTU686.
Redundancy analysis and a Monte Carlo permutation test were used to determine the relationships among six biochemical factors and bacterial community structure. All of the edaphic variables explained 81.2% of the variance, with axis 1 explaining 62.1% of the variance and axis 2 explaining another 11.3% (Figure 4). The major biochemical characteristics driving soil bacterial community composition were pH (F = 3.55, p = 0.001), soil moisture (F = 3.33, p = 0.001), Ca2+ (F = 2.90, p = 0.001), soil organic C (F = 2.724, p = 0.004), total N (F = 2.31, p = 0.004), and C/N (F = 2.051, p = 0.013).
FIGURE 4. Redundancy analysis to show correlation between the bacterial communities and physicochemical properties under five vegetation types.
Bacterial Co-occurrence Network Analysis
Across all 15 samples from five Karst vegetation types, correlation network analysis showed 219 strong positive correlations among 49 genera (R > 0.6, p ≤ 0.01; Figure 5). The network of positive correlations formed three distinct major (≥10 OTUs) modules of co-occurring taxa. The average path length between two nodes (APL) was 2.03 edges with a diameter of six edges. The clustering coefficient (CC) was 0.72 and the modularity index (MD) was 0.53, where MD > 0.4 suggests that the network has a modular structure.
FIGURE 5. Co-occurring network of bacterial communities across five vegetation types based on correlation analysis. The nodes in network (A) are colored by phylum. The nodes in network (B) are colored by modularity class. The connections stands for a strong (spearman’s ρ > 0.6) and significant (P < 0.01) correlations. The size of each node is proportional to the relative abundance of specific genus. The thickness of each edge is proportional to the ρ.
All genera in the network were assigned to bacteria phyla. Among these, Proteobacteria and Acidobacteria made up the two largest proportions, accounting for 38.78 and 20.41% of all nodes, respectively. When the distribution of nodes was modularized, all nodes were classified into three major modules (>10 nodes). Based on betweenness centrality scores, the top three genera identified were Bryobacter, GR-WP33-30, and Rhizomicrobium, which indicates the critical roles these microbes play as keystone taxa in the co-occurrence network.
Modules varied in their environmental associations, which indicates that taxa from the same module were ecologically similar (Supplementary Table S1). For instance, the abundance of genera in module I showed no significant correlations with any environmental factors, whereas the other two modules strongly correlated with such factors. In addition, the abundances of genera in the three major modules were negatively correlated with soil pH, nutrients, and water content, which suggests that these modules were composed primarily of oligotrophic, acidophilic, and drought-tolerant taxa.
Responses of Soil Environmental Properties to the Vegetation Types
Both vegetation succession and afforestation can enhance soil organic C and N dynamics by increasing SOM input and decreasing the decomposition rate (Grünzweig et al., 2007; Morris et al., 2007; Zhu et al., 2012; Cheng and An, 2015; Nadal-Romero et al., 2016). Furthermore, vegetation restoration can increase the water-holding capacity of the soil (Selma, 2014; Zhang et al., 2016). In the present study, TN, SOC, and SM content generally increased as plant communities changed from grass to shrub and forest, which suggests that vegetation restoration accelerates organic matter accumulation and improves soil moisture conditions in Karst ecosystems. By contrast, a decreasing trend in pH was observed as vegetation restoration progressed, in accordance with previous studies (Robertson and Vitousek, 1981; Schipper et al., 2001; Dmowska and Ilieva-Makulec, 2006; Holtkamp et al., 2008; Zhao et al., 2014). Compared to Quercus aquifolioides × Lithocarpus dealbatus secondary forest, Pinus armandii pure plantation and Pinus yunnanensis × Alnus ferdinandi-coburgii mixed plantation exhibited lower pH values, which may have been due to acidification caused by coniferous afforestation (Brand et al., 1986; Jönsson et al., 2003). The native plants in Karst calcareous soil are mainly calcicoles, which accumulate amount of calcium in their leaves (White and Broadley, 2003). In the present study, a lower calcium accumulation capacity in these two calcifuge Pinus foliage may have reduced the calcium content of leaf litter and consequently by a decrease of pH value in topsoil. Besides, various origins of root-mediated changes such as organic acid may also affect soil pH (Hinsinger et al., 2003). Tree species was the most important factor determining the C:N ratio in European forest soils (Cools et al., 2014). In this study, the Pinus armandii plantation and Pinus yunnanensis × Alnus ferdinandi-coburgii mixed plantation had higher soil C:N ratios than other vegetation types, perhaps because of the divergent chemical characteristics of conifers compared to other plant species (Laganière et al., 2013; Hume et al., 2016). These results suggest that both vegetation restoration stage and chemical characteristics of the plant species present impact soil properties.
Vegetation-Associated Edaphic Impacts on Bacterial Community Structure
The critical role of pH in shaping bacterial community structure is well characterized (Fierer and Jackson, 2006; Lauber et al., 2009; Chu et al., 2010; Griffiths et al., 2011; Ling et al., 2016). In the present study, soil pH was closely correlated with bacterial diversity across Karst revegetation sites. Moreover, soil microbial communities were influenced primarily by soil pH. The strong correlation between soil pH and microbial distribution could be due to the relatively narrow growth tolerances exhibited by most bacterial taxa. Indeed, each type of microorganism has an optimal pH value, and a slight change in pH might favor distinct bacterial taxa. Therefore, pH is a universal factor for predicting bacterial diversity and community structure (Fierer and Jackson, 2006; Singh et al., 2012). We observed a significant increase in the abundance of Acidobacteria with decreased pH, which has also been observed across terrestrial soil types (Eichorst et al., 2007). Calcium ion is implicated in a broad array of bacterial functions, including heat shock, pathogenicity, chemotaxis, differentiation, and the cell cycle (Norris et al., 1996). Ca2+ is closely related to bacterial community structure in typical Chinese forest soils (Xia et al., 2016). In this study, bacterial richness and diversity demonstrated significant positive correlations with Ca2+, which implies that the available calcium content acts as a determining factor in shaping soil bacterial populations and activity. In addition, soil moisture is an important factor driving microbial diversity across Karst vegetation types. Changes in soil water conditions affect oxygen content and substrate availability and consequently the microbial community (Yu et al., 2012; Banerjee et al., 2016). This finding is in agreement with previous studies in Antarctic soils (Yergeau et al., 2007; Chong et al., 2010), Canadian low Arctic tundra (Chu et al., 2011), the Yellow River Estuary in China (Yu et al., 2012), and Beilu River permafrost soils (Zhang et al., 2013).
In our study, soil bacterial diversity was negatively correlated with soil nutrient concentrations (SOC and TN). This pattern may indicate that low-resource environments in grassland and shrubbery lead to more unique niches, whereas relatively high-resource forestland habitats are less inclined to microbial niche differentiation, which has a major impact on microbial diversity (Bakker et al., 2013; Schlatter et al., 2015). Soil resource elemental stoichiometry plays an essential role in bacterial diversity and community composition (Högberg et al., 2007; Wan et al., 2015; Delgado-Baquerizo et al., 2017). In this study, significant relationships of C:N ratio with bacterial diversity and community structure were also observed, which indicates that variables associated with nitrogen transformations may be crucial determinants of bacterial community structure (Shi et al., 2015). Taken together, these results indicate that vegetation-associated soil properties play a vital role in determining bacterial community composition.
Co-occurrence Patterns of Bacteria during Vegetation Restoration
The most abundant phyla studied in the co-occurrence network were Acidobacteria and Proteobacteria, which indicates that these generalists are adapted to a variety of environments (Jiao et al., 2016). Moreover, the CC value of 0.72 for the bacterial co-occurrence network in this Karst region was higher than those reported for other ecosystems (Ju et al., 2014; Peng et al., 2014), which demonstrates stronger correlations in this ecological network. The fact that Bryobacter, GR-WP33-30, and Rhizomicrobium had the top three betweenness centrality values indicates the importance of these nodes in the co-occurrence network. Bryobacter, an aerobic chemo-organotrophic bacterium that utilizes various sugars, polysaccharides, and organic acids, plays an important role in the biogeochemical carbon cycle (Dedysh et al., 2016). Rhizomicrobium is a symbiotic mycorrhizophere bacterium and is crucial for nitrogen fixation (Ueki et al., 2010). GR-WP33-30 has also been identified as a keystone species in previous research, and future work is needed to better understand the role of this bacterium in co-occurrence networks (Ma et al., 2016). Thus, all three keystone taxa may play critical roles in ecological processes.
Microbial species interactions are crucial to the structure and dynamics of soil bacterial communities (Czárán et al., 2002; Prosser et al., 2007; Hibbing et al., 2010; Schlatter et al., 2015). Correlation-based network analysis has been successfully used to explore the co-occurrence patterns of microbial communities (Barberán et al., 2012; Purkamo et al., 2015; Chao et al., 2016; Ma et al., 2016; Jiao et al., 2016; Pérez-Valera et al., 2017). Network analysis showed non-random co-occurrence patterns in microbial communities. Note that the 10 strongest positive correlations were all between different phyla, which indicates that metabolic cooperation may play an important role in shaping species co-occurrence (Zelezniak et al., 2015). Moreover, the different correlations between the abundance of genera and soil characteristics among modules indicate specific ecological characteristics in these assemblies. These non-random assembly patterns of bacteria indicate the dominance of species interactions and environmental filtering in shaping community structure.
In the present study, we determined the effects of vegetation restoration and microbial interactions on the structure of bacterial communities in five vegetation types in KRD areas. Our work suggested that variations in soil physiochemical properties following revegetation led to shifts in the structure of the bacterial communities. Coniferous-broadleaved mixed plantation showed more similar bacterial community structure with secondary forest than coniferous plantation, indicating a better effect on ecological restoration of KRD land. Soil pH emerged as the major determinant of bacterial community characteristics. Bryobacter, GR-WP33-30, and Rhizomicrobium are the keystone taxa in KRD areas of southwestern China. Finally, we demonstrated non-random co-occurrence and modular patterns of bacterial communities. This information improves understanding of microbial responses to vegetation restoration in degraded karst regions.
In this study work, SL was responsible for experiment design and writing guidance, LX and HR were responsible for experiment performance and paper writing, XL and XY were responsible for experimental data processing and analysis.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the Fundamental Research Funds for the Central Non-profit Research Institution of CAF (RISF2014007), the National Key R&D Program of China (2016YFC0502605-3), the National Natural Science Foundation of China (31470632; 31600586) and the Fundamental Research Funds for the Central Non-profit Research Institution of CAF (CAFYBB2017MA001 and CAFYBB2016QA002).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.02377/full#supplementary-material
FIGURE S1 | Rarefaction analysis for each sample.
FIGURE S2 | The proportion of the bacterial phyla/subphyla under five vegetation types. The phylum Proteobacteria contains the four subphyla of Alphaproteobacteria, Betaproteobacteria, Deltaproteobacteria and Gammaproteobacteria. G, Grassland; S shrubbery; SF, secondary forest; MP, mixed plantation; PP, pure plantation.
- ^ http://drive5.com/uparse/
- ^ http://rdp.cme.msu.edu
- ^ http://bonsai.hgc.jp/mdehoon/software/cluster/index.html
- ^ http://sourceforge.net/projects/jtreeview.htm
Bakker, M. G., Otto-Hanson, L., Lange, A. J., Bradeen, J. M., and Kinkel, L. L. (2013). Plant monocultures produce more antagonistic soil Streptomyces communities than high-diversity plant communities. Soil Biol. Biochem. 65, 304–312. doi: 10.1016/j.soilbio.2013.06.007
Banerjee, S., Helgason, B., Wang, L., Winsley, T., Ferrari, B. C., and Siciliano, S. D. (2016). Legacy effects of soil moisture on microbial community structure and N2O emissions. Soil Biol. Biochem. 95, 40–50. doi: 10.1007/s11356-017-8799-6
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chao, Y., Liu, W., Chen, Y., Chen, W., Zhao, L., Ding, Q., et al. (2016). Structure, variation, and co-occurrence of soil microbial communities in abandoned sites of a rare earth elements mine. Environ. Sci. Technol. 50, 11481–11490. doi: 10.1021/acs.est.6b02284
Cheng, M., and An, S. (2015). Responses of soil nitrogen, phosphorous and organic matter to vegetation succession on the Loess Plateau of China. J. Arid Land 7, 216–223. doi: 10.1007/s40333-014-0043-3
Chong, C. W., Pearce, D. A., Convey, P., Tan, G. Y. A., Wong, R. C. S., and Tan, I. K. P. (2010). High levels of spatial heterogeneity in the biodiversity of soil prokaryotes on Signy Island, Antarctica. Soil Biol. Biochem. 42, 601–610. doi: 10.1016/j.soilbio.2009.12.009
Chu, H., Neufeld, J. D., Walker, V. K., and Grogan, P. (2011). The influence of vegetation type on the dominant soil bacteria, archaea, and fungi in a low arctic tundra landscape. Soil Biol. Biochem. 75, 1756–1765. doi: 10.2136/sssaj2011.0057
Chu, H. Y., Fierer, N., Lauber, C. L., Caporaso, J. G., Knight, R., and Grogan, P. (2010). Soil bacterial diversity in the Arctic is not fundamentally different from that found in other biomes. Environ. Microbiol. 12, 2998–3006. doi: 10.1111/j.1462-2920.2010.02277.x
Cools, N., Vesterdal, L., De Vos, B., Vanguelova, E., and Hansen, K. (2014). Tree species is the major factor explaining C:N ratios in European forest soils. For. Ecol. Manage. 311, 3–16. doi: 10.1016/j.foreco.2013.06.047
Crouzeilles, R., Curran, M., Ferreira, M. S., Lindenmayer, D. B., Grelle, C. E. V., and Benayas, G. M. R. (2016). A global meta-analysis on the ecological drivers of forest restoration success. Nat. Commun. 7:11666. doi: 10.1038/ncomms11666
Dedysh, S. N., Kulichevskaya, I. S., Huber, K. J., and Overmann, J. (2016). Defining the taxonomic status of described subdivision 3 Acidobacteria: the proposal of Bryobacteraceae fam. nov. Int. J. Syst. Evol. Microbiol. 67, 498–501. doi: 10.1099/ijsem.0.001687
Delgado-Baquerizo, M., Reich, P. B., Khachane, A. N., Campbell, C. D., Thomas, N., Freitag, T. E., et al. (2017). It is elemental: soil nutrient stoichiometry drives bacterial diversity. Environ. Microbiol. 19, 1176–1188. doi: 10.1111/1462-2920.13642
Dmowska, E., and Ilieva-Makulec, K. (2006). Secondary succession of nematodes in power plant ash dumps reclaimed by covering with turf. Eur. J. Soil Biol. 42, S164–S170. doi: 10.1016/j.ejsobi.2006.07.005
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Eichorst, S. A., Breznak, J. A., and Schmidt, T. M. (2007). Isolation and characterization of soil bacteria that define Terriglobus gen. nov., in the phylum Acidobacteria. Appl. Environ. Microbiol. 73, 2708–2717. doi: 10.1128/AEM.02140-06
Feng, X., Fu, B., Piao, S., Ciais, P., Zeng, Z., Lü, Y., et al. (2016). Revegetation in China’s Loess Plateau is approaching sustainable water resource limits. Nat. Clim. Chang. 6, 1019–1022. doi: 10.1038/nclimate3092
Griffiths, R. I., Thomson, B. C., James, P., Bell, T., Bailey, M., and Whiteley, A. S. (2011). The bacterial biogeography of British soils. Environ. Microbiol. 13, 1642–1654. doi: 10.1111/j.1462-2920.2011.02480.x
Grünzweig, J. M., Gelfand, I., Fried, Y., and Yakir, D. (2007). Biogeochemical factors contributing to enhanced carbon storage following afforestation of a semi-arid shrubland. Biogeosciences 4, 2111–2145. doi: 10.5194/bgd-4-2111-2007
Hinsinger, M., Plassard, C., Tang, C., and Jaillard, B. (2003). Origins of root-mediated pH changes in the rhizosphere and their responses to environmental constraints: a review. Plant Soil 248, 43–59. doi: 10.1023/A:1022371130939
Holtkamp, R., Kardol, P., Wal, A., Dekker, S. C., Putten, W. H., and Ruiter, P. C. (2008). Soil food web structure during ecosystem development after land abandonment. Appl. Soil Ecol. 39, 23–34. doi: 10.1016/j.apsoil.2007.11.002
Hume, A., Chen, H. Y. H., Taylor, A. R., Kayahara, G. J., and Man, R. (2016). Soil C:N:P dynamics during secondary succession following fire in the boreal forest of central Canada. For. Ecol. Manage. 369, 1–9. doi: 10.1016/j.foreco.2016.03.033
Jiao, S., Liu, Z., Lin, Y., Yang, J., Chen, W., and Wei, G. (2016). Bacterial communities in oil contaminated soils: Biogeography and co-occurrence patterns. Soil Biol. Biochem. 98, 64–73. doi: 10.1016/j.soilbio.2016.04.005
Jönsson, U., Rosengren, U., Thelin, G., and Nihlgård, B. (2003). Acidification-induced chemical changes in coniferous forest soils in southern Sweden 1988–1999. Environ. Pollut. 123, 75–83. doi: 10.1016/S0269-7491(02)00335-4
Ju, F., Xia, Y., Guo, F., Wang, Z., and Zhang, T. (2014). Taxonomic relatedness shapes bacterial assembly in activated sludge of globally distributed wastewater treatment plants. Environ. Microbiol. 16, 2421–2432. doi: 10.1111/1462-2920.12355
Laganière, J., Paré, D., Bergeron, Y., Chen, H. Y. H., Brassard, B. W., and Cavard, X. (2013). Stability of soil carbon stocks varies with forest composition in the Canadian boreal biome. Ecosystems 16, 852–865. doi: 10.1007/s10021-013-9658-z
Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. doi: 10.1128/AEM.00335-09
Lewis, D. E., White, J. R., Wafula, D., Athar, R., Dickerson, T., Williams, H. N., et al. (2010). Soil functional diversity analysis of a bauxite-mined restoration chronosequence. Microb. Ecol. 59, 710–723. doi: 10.1007/s00248-009-9621-x
Li, S., Birk, S., Xue, L., Ren, H., Chang, J., and Yao, X. (2016). Seasonal changes in the soil moisture distribution around bare rock outcrops within a karst rocky desertification area (Fuyuan County, Yunnan Province, China). Environ. Earth Sci. 75, 1482. doi: 10.1007/s12665-016-6290-1
Li, S., Ren, H. D., Xue, L., Chang, J., and Yao, X. H. (2014). Influence of bare rocks on surrounding soil moisture in the karst rocky desertification regions under drought conditions. Catena 116, 157–162. doi: 10.1016/j.catena.2013.12.013
Ling, N., Zhu, C., Xue, C., Chen, H., Duan, Y., Peng, C., et al. (2016). Insight into how organic amendments can shape the soil microbiome in long-term field experiments as revealed by network analysis. Soil Biol. Biochem. 99, 137–149. doi: 10.1016/j.soilbio.2016.05.005
Ma, B., Wang, H., Dsouza, M., Lou, J., He, Y., Dai, Z., et al. (2016). Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 10, 1891–1901. doi: 10.1038/ismej.2015.261
Morris, S. J., Bohm, S., Haile-Mariam, S., and Paul, E. A. (2007). Evaluation of carbon accrual in afforested agricultural soils. Glob. Change Biol. 13, 1145–1156. doi: 10.1111/j.1365-2486.2007.01359.x
Nadal-Romero, E., Cammeraat, E., Pérez-Cardiel, E., and Lasanta, T. (2016). Effects of secondary succession and afforestation practices on soil properties after cropland abandonment in humid Mediterranean mountain areas. Agric. Ecosyst. Environ. 228, 91–100. doi: 10.1016/j.agee.2016.05.003
Nelson, D. W., and Sommers, L. E. (1982). “Total carbon, organic carbon, and organic matter,” in Methods of Soil Analysis: Chemical and Microbiological Properties, eds A. L. Page, H. Millerr, and D. R. Keeney (Madison: America Society of Agronomy), 539–552.
Oh, Y. M., Kim, M., Lee-Cruz, L., Lai-Hoe, A., Go, R., Ainuddin, N., et al. (2012). Distinctive bacterial communities in the rhizo-plane of four tropical tree species. Microb. Ecol. 64, 1018–1027. doi: 10.1007/s00248-012-0082-2
Peng, X. X., Guo, F., Ju, F., and Zhang, T. (2014). Shifts in the microbial community, nitrifiers and denitrifiers in the biofilm in a full-scale rotating biological contactor. Environ. Sci. Technol. 48, 8044–8052. doi: 10.1021/es5017087
Pérez-Valera, E., Goberna, M., Faust, K., Raes, J., García, C., and Verdú, M. (2017). Fire modifies the phylogenetic structure of soil bacterial co-occurrence networks. Environ. Microbiol. 19, 317–327. doi: 10.1111/1462-2920.13609
Prosser, J. I., Bohannan, B. J., Curtis, T. P., Ellis, R. J., Firestone, M. K., and Freckleton, R. P. (2007). The role of ecological theory in microbial ecology. Nat. Rev. Microbiol. 5, 384–392. doi: 10.1038/nrmicro1643
Purkamo, L., Bomberg, M., Kietäväinen, R., Salavirta, H., Nyyssönen, M., Nuppunen-Puputti, M., et al. (2015). The keystone species of Precambrian deep bedrock biosphere belong to Burkholderiales and Clostridiales. Biogeosci. Discuss. 12, 18103–18150. doi: 10.5194/bgd-12-18103-2015
Qi, X., Wang, K., and Zhang, C. (2013). Effectiveness of ecological restoration projects in a Karst region of southwest China assessed using vegetation succession mapping. Ecol. Eng. 54, 245–253. doi: 10.1016/j.ecoleng.2013.01.002
Schipper, L. A., Degens, B. P., Sparling, G. P., and Duncan, L. C. (2001). Changes in microbial heterotrophic diversity along five plant successional sequences. Soil Biol. Biochem. 33, 2093–2103. doi: 10.1016/S0038-0717(01)00142-0
Schlatter, D. C., Bakker, M. G., Bradeen, J. M., and Kinkel, L. L. (2015). Plant community richness and microbial interactions structure bacterial communities in soil. Ecology 96, 134–142. doi: 10.1890/13-1648.1
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Shi, Y., Xiang, X., Shen, C., Chu, H., Neufeld, J. D., Walker, V. K., et al. (2015). Vegetation-associated impacts on arctic tundra bacterial and microeukaryotic communities. Appl. Environ. Microbiol. 81, 492–501. doi: 10.1128/AEM.03229-14
Singh, D., Takahashi, K., Kim, M., Chun, J., and Adams, J. M. (2012). A hump-backed trend in bacterial diversity with elevation on mount Fuji, Japan. Microb. Ecol. 63, 429–437. doi: 10.1007/s00248-011-9900-1
Ueki, A., Kodama, Y., Kaku, N., Shiromura, T., Satoh, A., Watanabe, K., et al. (2010). Rhizomicrobium palustre gen. nov., sp. nov., a facultatively anaerobic, fermentative stalked bacterium in the class Alphaproteobacteria isolated from rice plant roots. J. Gen. Appl. Microbiol. 56, 193–203. doi: 10.2323/jgam.56.193
Vick-Majors, T. J., Priscu, J. C. A., Amaral-Zettler, L. (2014). Modular community structure suggests metabolic plasticity during the transition to polar night in ice-covered Antarctic lakes. ISME J. 8, 778–789. doi: 10.1038/ismej.2013.190
Wan, X., Huang, Z., He, Z., Yu, Z., Wang, M., Davis, M. R., et al. (2015). Soil C:N ratio is the major determinant of soil microbial community structure in subtropical coniferous and broadleaf forest plantations. Plant Soil 387, 103–116. doi: 10.1007/s11104-014-2277-4
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, S. J., Liu, Q. M., and Zhang, D. F. (2004). Karst rocky desertification in southwestern China: geomorphology, landuse, impact and rehabilitation. Land Degrad. Dev. 15, 115–121. doi: 10.1002/ldr.592
Xia, Z., Bai, E., Wang, Q., Gao, D., Zhou, J., Jiang, P., et al. (2016). Biogeographic distribution patterns of bacteria in typical Chinese forest soils. Front. Microbiol. 7:1106. doi: 10.3389/fmicb.2016.01106
Xiong, J., Liu, Y., Lin, X., Zhang, H., Zeng, J., Hou, J., et al. (2012). Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau. Environ. Microbiol. 14, 2457–2466. doi: 10.1111/j.1462-2920.2012.02799.x
Yergeau, E., Bokhorst, S., Huiskes, A. H. L., Boschker, H. T. S., Aerts, R., and Kowalchuk, G. A. (2007). Size and structure of bacterial, fungal and nematode communities along an Antarctic environmental gradient. FEMS Microbiol. Ecol. 59, 436–451. doi: 10.1111/j.1574-6941.2006.00200.x
Yu, Y., Wang, H., Liu, J., Wang, Q., Shen, T., Guo, W., et al. (2012). Shifts in microbial community function and structure along the successional gradient of coastal wetlands in Yellow River Estuary. Eur. J. Soil Biol. 49, 12–21. doi: 10.1016/j.ejsobi.2011.08.006
Zak, D. R., Holmes, W. E., White, D. C., Peacock, A. D., and Tilman, D. (2003). Plant diversity, soil microbial communities, and ecosystem function: are there any links? Ecology 84, 2042–2050. doi: 10.1890/02-0433
Zelezniak, A., Andrejev, S., Ponomarova, O., Mende, D. R., Bork, P., and Patil, K. R. (2015). Metabolic dependencies drive species co-occurrence in diverse microbial communities. Proc. Natl. Acad. Sci. U.S.A. 112, 6449–6454. doi: 10.1073/pnas.1421834112
Zhang, X. F., Zhao, L., Xu, S. J. Jr., Liu, Y. Z., and Cheng, G. D. (2013). Soil moisture effect on bacterial and fungal community in Beilu River (Tibetan Plateau) permafrost soils with different vegetation types. J. Appl. Microbiol. 114, 1054–1065. doi: 10.1111/jam.12106
Zhang, Y. W., Deng, L., Yan, W. M., and Shangguan, Z. P. (2016). Interaction of soil water storage dynamics and long-term natural vegetation succession on the Loess Plateau, China. Catena 137, 52–60. doi: 10.1016/j.catena.2015.08.016
Zhao, J., Li, S., He, X., Liu, L., and Wang, K. (2014). The soil biota composition along a progressive succession of secondary vegetation in a karst area. PLOS ONE 9:e112436. doi: 10.1371/journal.pone.0112436
Keywords: bacterial communities, co-occurrence, vegetation restoration, karst, rocky desertification
Citation: Xue L, Ren H, Li S, Leng X and Yao X (2017) Soil Bacterial Community Structure and Co-occurrence Pattern during Vegetation Restoration in Karst Rocky Desertification Area. Front. Microbiol. 8:2377. doi: 10.3389/fmicb.2017.02377
Received: 04 June 2017; Accepted: 16 November 2017;
Published: 01 December 2017.
Edited by:Frank Rasche, University of Hohenheim, Germany
Reviewed by:Kim Yrjälä, University of Helsinki, Finland
Roberta Fulthorpe, University of Toronto Scarborough, Canada
Copyright © 2017 Xue, Ren, Li, Leng and Yao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sheng Li, email@example.com
†These authors have contributed equally to this work.