Next Article in Journal
Competition, Climate, and Size Effects on Radial Growth in an Old-Growth Hemlock Forest
Next Article in Special Issue
CO2 Transfer Characteristics of Calcareous Humid Subtropical Forest Soils and Associated Contributions to Carbon Source and Sink in Guilin, Southwest China
Previous Article in Journal
Ground Beetle (Coleoptera: Carabidae) Response to Harvest Residue Retention: Implications for Sustainable Forest Bioenergy Production
Previous Article in Special Issue
Effects of Plum Plantation Ages on Soil Organic Carbon Mineralization in the Karst Rocky Desertification Ecosystem of Southwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Land Use Patterns on Soil Bacterial Community Structure in the Karst Graben Basin of Yunnan Province, China

1
College of Environmental Science and Engineering, Guilin University of Technology, Guilin 541004, China
2
Key Laboratory of Karst Dynamics, MNR & GZAR, Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004, China
3
International Research Center on Karst under the Auspices of UNESCO, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Forests 2020, 11(1), 51; https://doi.org/10.3390/f11010051
Submission received: 21 November 2019 / Revised: 29 December 2019 / Accepted: 30 December 2019 / Published: 31 December 2019

Abstract

:
Land use patterns can change the structure of soil bacterial communities. However, there are few studies on the effects of land use patterns coupled with soil depth on soil bacterial communities in the karst graben basin of Yunnan province, China. Consequently, to reveal the structure of the soil bacterial community at different soil depths across land use changes in the graben basins of the Yunnan plateau, the relationship between soil bacterial communities and soil physicochemical properties was investigated for a given area containing woodland, shrubland, and grassland in Yunnan province by using next-generation sequencing technologies coupled with soil physicochemical analysis. Our results indicated that the total phosphorus (TP), available potassium (AK), exchangeable magnesium (E-Mg), and electrical conductivity (EC) in the grassland were significantly higher than those in the woodland and shrubland, yet the total nitrogen (TN) and soil organic carbon (SOC) in the woodland were higher than those in the shrubland and grassland. Proteobacteria, Verrucomicrobia, and Acidobacteria were the dominant bacteria, and their relative abundances were different in the three land use types. SOC, TN, and AK were the most important factors affecting soil bacterial communities. Land use exerts strong effects on the soil bacterial community structure in the soil’s surface layer, and the effects of land use attenuation decrease with soil depth. The nutrient content of the soil surface layer was higher than that of the deep layer, which was more suitable for the survival and reproduction of bacteria in the surface layer.

1. Introduction

The karst graben basin in the East Yunnan plateau is a special geomorphological formation due to the rifts of the plateau, which were uplifted at the same time [1,2]. This area is also the main area of the “two barriers and three belts” for China’s national ecological security. However, due to deforestation of this area over the past several decades, the karst ecosystem has seriously degenerated, thereby affecting soil quality, soil fertility, and ecological conditions, and resulting in abandoned bare land. To fight against this environmental problem, the “Green for Grain” program, or the Natural Forest Protection Project, was launched by the Chinese government in this region [3,4,5]. Accordingly, the size of the degenerated areas has shrunk with the revegetation process. However, little is known about how vegetation restoration types affect the soil bacterial community’s structure in the karst graben basin.
Vegetation restoration can affect soil microorganism communities, which can regulate the soil’s biogeochemical cycles and affect the stability of the soil’s ecosystem [6,7,8]. Although many studies have discussed soil microbial community structures with land use pattern changes worldwide, there are few studies on the influence of land use patterns on the soil bacterial community’s structure in karst graben basins. For example, Suleiman et al. compared the bacterial community in the original forest-covered area and grassland for eight years and found that the main bacterial community composition showed little difference [9]. Song et al. found that the number and composition of the soil microbial population in the karst peak-cluster depression were different in farmland, grassland, scrubland, forest plantation, secondary forest, and primary forest [10]. Ederson et al. found that there were significant differences in the community composition among crops, pastures, and agroforestry, as well as young secondary forest (up to 5 years old), old secondary forest (5 to 30 years old), and primary forest sites [11]. To better explore the influence of land use patterns on the soil bacterial community structure in the Luxi county in the Yunnan karst graben basin, the woodland, shrubland, and grassland in a given karst area were selected. Our objectives were to (i) gain insight into the difference of soil bacterial community structure with land use changes, (ii) inquire into the effects of land use patterns with soil depth, and (iii) identify the key factors in determining soil bacterial communities. Our findings provide a basis for understanding the influence of land use patterns on soil bacterial community structures in the karst graben basin of Yunnan province, China.

2. Materials and Methods

2.1. Study Sites

Luxi county (103°30′–104°03′ E, 24°15′–24°46′ N) is located in the north Honghe Hani and Yi Autonomous Prefecture in Yunnan province in the subtropical monsoon climate zone, which is rainy in summer and dry in winter. The rainy season extends from May to October, and the dry season runs from March to April. The average annual precipitation is 2026.5 mm. The average annual temperature is 16.3 °C. Moreover, the rocky desertification in Luxi country is serious.

2.2. Soil Sample Collection

Sampling occurred between the wet and dry season, January 2018. Three 20 × 20 m plots were established for each land use pattern. The minimum distance between plots was 400 m to avoid pseudoreplication. Soil samples were collected from the surface soil (0–10 cm), named the A layer, and the other samples were taken from the deep layer (10–20 cm), named B layer, with a split tube auger 5 cm in diameter. At each plot, six soil cores were selected in an S-shaped pattern to form one soil sample. A total of 18 soil samples from woodland (WL), shrubland (SL), and grassland (GL) were collected. Soil samples were named according to the soil sampling sites (such as SL) and soil layer (A: 0–10 cm) in that particular order (e.g., SLA). The three land use patterns were continuously managed for 15 years. All soil samples were transported to the laboratory immediately after collection in sterile plastic bags on dry ice and divided into two uniform portions. One portion was stored at −80 °C for DNA analysis, and the other part was air dried for physicochemical analysis.

2.3. Physicochemical Analysis

Soil moisture content (moisture), soil temperature (T), and EC were measured in situ by soil three-parameter tachometers (Delta-T Devices Ltd., Moisture Meter type HH2, UK). Soil organic carbon (SOC) was determined by wet digestion using the H2SO4 and K2Cr2O7 method [12]. Total phosphorus (TP) was determined using the molybdenum blue colorimetric method following HClO4 digestion [13]. Available potassium (AK) was extracted with ammonium acetate and analyzed using a flame photometer [14]. Total nitrogen (TN) was determined by an element analyzer [15]. Soil pH was determined by a 1:5 (m:V) soil/water ratio and measured by a corrected desktop pH meter of Maitre-Toledo S470-B Seven Excellence [16]. Exchangeable calcium (E-Ca) and exchangeable magnesium (E-Mg) were determined by ammonium acetate exchange-atomic absorption spectrophotometry [17].

2.4. DNA Extraction

DNA was extracted using the Powersoil DNA Isolation Kit (Mobio Laboratories, Inc., Carlsbad, CA, USA). For next-generation sequencing, the V3–V4 region of 16S rRNA genes was amplified using PCR primers 338F (ACTCCTACGGGAGGCAGCA) and 806R (GGACTACHVGGGTWTCTA AT). The PCR products targeting the V3–V4 region of the 16S rRNA genes were purified by using the TIANquick Maxi Purification Kit (TIANGEN Biotech (Beijing) Co., Ltd., China). Then, 16S rRNA gene sequencing was performed on the Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA, USA) at the MAGIENE (Guangzhou, China).

2.5. Bioinformatic Analysis and Statistical Analysis

According to the barcode sequence, the sequencing data were divided into different sample data, and the barcode sequence was truncated. After splicing each sample with FLASH 1.2.11 software [18], the Cutadapt 1.9.1 software was employed to truncate the sequence of PCR amplified primers and remove fragments (less than 200 bp) [19]. Using the SILVA SSURef 123 NR database as the reference database, chimeric sequences were removed by the UCHIME 4.2 software [20,21]. Afterwards, the processed sequences were clustered with the UCLUST 1.2.22 software in operational taxonomic units (OTUs) according to sequences with more than a 97% similarity, and the OTUs were classified by the UCLUST 1.2.22 software alignment against the most recent SILVA (123) database [22]. The QIIME 1.9.1 software was used to perform an alpha diversity (Chao1, Simpson, Shannon, and Observed OTUs) and beta diversity test on the OTU table [23].
SPSS 25 software (SPSS Inc., Chicago, IL, USA) for a one-way ANOVA was used to analyze the difference soil physicochemical properties and bacterial community structure and diversity in different land use patterns [24]. The least significant difference method was used for multiple comparisons, and the Pearson correlation coefficient method was used for a correlation analysis. The Origin 8.5 software was used to map the abundance of bacterial communities at the phyla level. The OTUs whose numbers were more than 0.5% of the total OTUs were defined as the dominant OTUs. RStudio 3.0.3 software was used to draw the heat map of the dominant OTUs and perform a distance-based redundancy analysis (db-RDA) between the soil bacteria and soil physicochemical parameters [25]. The similarity of the OTUs among samples was analyzed by using the Bray–Curtis and weighted UniFrac distance algorithm of principal coordinate analysis (PCoA) [26]. The network maps of dominant OTUs and soil physicochemical factors were drawn via the Gephi 0.9.2 software [27].

3. Results

3.1. Soil Physicochemical Parameters with Land Use Changes

It can be seen that the soil’s physicochemical properties were different in the woodland, shrubland, and grassland (Table 1). The TP, AK, E-Mg, EC, and T in grassland were higher than those in the woodland and shrubland. The TN and SOC in the woodland were higher than those in the shrubland and grassland. Moreover, the soil physicochemical properties, except for some physicochemical properties in woodland, decreased by increasing soil depth.

3.2. Soil Bacterial Community Structure and Diversity

The dominant phyla were different in the three land use types, as well as in the A and B layers, as seen in Figure 1. There were 13 dominant phyla in the A layer (Figure 1a). Proteobacteria, Verrucomicrobia, and Acidobacteria were the most abundant dominant bacteria. The total proportions of the three dominant bacteria in the A layer were woodland (80.22%), shrubland (80.07%), and grassland (59.52%); the Verrucomicrobia in grassland were significantly fewer than those in the woodland and shrubland; Acidobacteria, Actinobacteria, Bacteroidetes, and Chloroflexi in grassland were significantly higher than those in woodland and shrubland (Table S1). In the B layer, the total proportions of the three most abundant dominant bacteria were woodland (78.61%), shrubland (67.14%), and grassland (56.03%). Among the three different land use patterns, the composition and proportions of the other dominant phyla were different but not significant.
Moreover, the heat map clearly shows the distribution of the dominant OTUs (Figure 2). In the A layer (Figure 2a), OTU 410 (DA101 soil group) and OTU 24 (Candidatus Solibacter) were the dominant OTUs for the woodland; OTUs 238 and 30 (Candidatus Solibacter), OTUs 28 and 124 (Sphingomonas), OTUs 49, 69, and 27 (subgroup 6), and OTUs 58, 368, 26, and 47 (RB41) were the dominant OTUs for grassland; and OTUs 1, 31, 2312, 9, and 823 (DA101 soil group), and OTU 12 (Acidobacteria) were the dominant OTUs for shrubland. In the B layer (Figure 2b), OTUs 163, 823, 2, 1, and 772 (DA101 soil group), and OTUs 10 and 2858 (Acidobacteriaceae) were the dominant OTUs in shrubland. OTUs 65 and 16 (0319-6A21) and OTUs 53 and 11 (Acidobacteriaceae) were the dominant OTUs in woodland. OTU 21 (Candidatus Xiphinematobacter), OTU 71 (Gemmatimonadaceae), and OTU 54 (Pedomicrobium) were the dominant OTUs in grassland.
The alpha diversities of the soil bacterial communities were different between the three land use types. In the A layer, except for the Simpson index, the alpha diversity indices from grassland were significantly higher than those from shrubland and woodland (p < 0.05), as listed in Table 2. Moreover, there was no significant difference between shrubland and woodland (p > 0.05). In the B layer, the alpha diversity indices of the three land use patterns were not significantly different, except for the Simpson index. The alpha diversities decreased with an increase in soil depth, except for the Simpson index.
To investigate the effects of land use type and soil depth on soil bacterial communities, we examined the beta diversity (Figure 3). The shrubland and woodland can be accurately clustered together, which shows the similarity of their bacterial communities’ compositions. Compared with the A layer, the points in the B layer were more dispersed.

3.3. The Relationship between Soil Physicochemical Parameters and Soil Bacteria

To assess the influence of land uses on bacterial community composition, we performed a db-RDA on the dominant bacterial phyla. According to the db-RDA plot and Spearman’s correlations between the soil physicochemical parameters and dominant phyla in the land use types, the first two axes accounted for 2.86% of the variability of the bacterial community structure in the A layer. Further, the bacterial communities were positively correlated with six soil physicochemical properties (including SOC, TN, TP, AK, EC, and T) (Figure 4a). Verrucomicrobia were negatively correlated with TP, T, E-Mg, AK, and EC, and positively correlated with TN and SOC. Acidobacteria were positively correlated with EC, T, and E-Mg; Actinobacteria were positively correlated with AK, EC, T, and E-Mg; Bacteroidetes were positively correlated with E-Mg, AK, and EC; Chloroflexi were positively correlated with T, E-Mg, TP, AK, and EC; and Actinobacteria, Bacteroidetes, and Chloroflexi were negatively correlated with TN and SOC (Figure 5). In the B layer, db-RDA showed that the first two axes accounted for 2.86% of the variability of the bacterial community structure, and bacterial communities were positively correlated with SOC, TN, AK, and T (Figure 4b). Verrucomicrobia were positively correlated with SOC and negatively correlated with AK, T, and E-Mg. Acidobacteria were negatively correlated with TP; Chloroflexi were negatively correlated with TN and SOC, and positively correlated with T and E-Mg; Actinobacteria were positively correlated with AK, TP, T, E-Mg, and EC; and Bacteroidetes were positively correlated with AK and TP (Figure 5).
Heat maps showed that TP, T, E-Mg, AK, EC, TN, and SOC were significantly correlated with most of the dominant OTUs, whereas pH and E-Ca were significantly correlated with some OTUs (Figure 4c,d). In the A layer, OTUs 26, 47, 58, and 368 (RB41); OTUs 27 and 49 (subgroup 6); and OTUs 30 and 238 (Candidatus Solibacter) of Acidobacteria were significantly correlated with TP, T, E-Mg, AK, and EC (Figure 4c). Moreover, OTUs 28 and 124 (Sphingomonas), OTU 17 (Lysobacter), OTU 54 (Pedomicrobium), and OTU 43 (Steroidobacter) of Proteobacteria were also significantly correlated with TP, T, E-Mg, AK, and EC (Figure 4c). In the A layer, OTUs 1, 2, 9, 31, 52, 410, 823, and 2312 (DA101 soil group) and OTUs 114 and 143 (RB41) of Verrucomicrobia and Acidobacteria were significantly correlated with SOC and TN (Figure 4c). In the B layer, OTU 26 (RB41) and OTUs 27 and 69 (subgroup 6) of Acidobacteria were significantly correlated with TP, T, E-Mg, AK, and EC (Figure 4d), and OTUs 1, 2, 9, 31, 52, 410, 772, 823, and 2312 (DA101 soil group) of Verrucomicrobia were significantly correlated with SOC and TN (Figure 4d), as confirmed by their significant relationships in Tables S2 and S3. Although representing the relationship between soil physicochemical parameters and despite the observation that the most abundant OTUs were different in the three land uses, these tables, on the whole, show a certain regularity (Figure 4c,d and Figure S1). The OTU co-occurrence patterns and the relationships among soil bacterial communities and physicochemical parameters were also investigated. Overall, the ecological networks were markedly different among the bacterial groups at different soil depths. This network was comprised of highly connected genera and soil physicochemical properties, thereby forming a “small world” topology. The nodes in the network were assigned to 11 bacterial phyla (Figure 4e,f). Among these, three phyla (Proteobacteria, Acidobacteria, and Verrucomicrobia) were widely distributed, accounting for over 60% of all nodes. In the surface layer, OTUs 26, 47, and 58 (RB41) and OTUs 1, 2, 9, 31, 410, 823, and 2312 (DA101 soil group) of Acidobacteria and Verrucomicrobia were keystone taxa (Figure 4e). In the deep layer, OTUs 1, 2, 9, and 2312 of Acidobacteria were keystone taxa (Figure 4f). This study indicates that these OTUs may play key roles in maintaining the structure and function of soil bacterial communities. It can be seen that EC, E-Mg, SOC, AK, T, and TN were the most important soil physicochemical factors affecting the bacterial community composition in the A layer, whereas SOC, AK, E-Mg, and TN were the most important factors in the B layer. The networks of competition and cooperation were more complex among bacteria in the surface layer than those in the deep layer. The correlation between bacterial communities and soil physicochemical properties decreased with an increase of soil depth (Figure 4c–f).

4. Discussion

4.1. The Characteristics of the Soil Physicochemical Properties

Land use patterns determine the surface vegetation and soil management method, which in turn result in a difference between soil physicochemical properties [28,29]. In our study, TP, AK, E-Mg, EC, and the moisture in grassland were significantly higher than those in the woodland and shrubland. However, the TN and SOC in woodland were significantly higher than those in shrubland and grassland. It is well known that microorganisms are intimately involved in rock weathering and could use these elements as electron acceptors and nutrients [30]. Consequently, grassland, during the early stages of vegetation succession under the influence of microorganisms on rock weathering, had a high concentration of AK, E-Mg, and TP compared with shrubland and woodland. However, due to the deficient root systems and lower amount of litter in grassland, the SOC and TN were low in the grassland. Moreover, SOC and TN are elements of soil fertility and are closely related to ecosystem stability, environmental protection, and land use [31,32]. Compared with grassland, the litter in woodland was high, and the root system was dense. Moreover, TN and SOC were also positively correlated. Therefore, the contents of TN and SOC in woodland were higher than those in shrubland and grassland. Further, soil nutrient contents decreased gradually with increasing soil depth because the soil’s surface layer was the humus layer with rich litter, high nutrient content, a well-developed plant root system, good ventilation, and positive hydrothermal conditions [33,34,35,36].

4.2. Distribution of Bacterial Diversity Compositions

Bacterial communities from different land use types and different soil depths were highly diverse. In our study, soil bacterial diversity indices from grassland in the A layer and B layer were significantly higher than those in shrubland and woodland (p < 0.05). The diversity indices of soil bacterial communities in the three land use types decreased with increasing soil depth, which was consistent with the changes in soil physicochemical properties. Plant community richness usually increases during late successive stages. Moreover, increased plant community richness can significantly alter soil bacterial community composition and is negatively correlated with bacterial diversity [37]. An increased plant community can supply a diverse array of resources to the soil, thereby promoting coevolutionary niche differentiation and favoring nonantagonistic microbial communities in which antagonism plays little role in maintaining soil community diversity [38]. The decrease in bacterial diversity with increasing soil depth is due to the reduction in the nutrient substrates for bacteria [39]. This is consistent with the expected depth patterns of soil bacterial diversity found in other studies [40,41].
Beta diversity was used to describe the similarities and differences in community structure. The PCoA distance showed a pronounced influence of land use changes on soil bacterial community structure. In our study, bacterial communities from grassland were significantly different from those in woodland and shrubland. This shows that the rate of species replacement between grassland and shrubland was relatively fast, possibly due to the lack of soil nutrients in grassland and increased competition during successive stages, which intensifies the species replacement between communities [42].

4.3. Relationships of Bacterial Communities with Basic Soil Parameters

Soil physicochemical properties determine the living environment of bacteria, which affect soil bacterial communities [43]. Soil bacteria play an extremely important role in organic matter decomposition, soil nutrient cycling, and ecological environment improvement [44,45]. Changes in bacterial community compositions led us to evaluate the extent of different land uses. Then, in our study, 33 dominant phyla were found in the 18 soil samples. Proteobacteria, Verrucomicrobia, and Acrobacteria were the most abundant dominant soil bacteria from the three land use patterns. These phyla have also been found in different relative proportions in other ecosystems worldwide [46,47], suggesting that they may play a fundamental role in these environments.
Proteobacteria showed no significant difference in the three land use patterns, yet their proportion in surface soil (26.61%) was significantly higher than that in deep soil (20.95%), which could be related to their copiotrophic living attributes. Proteobacteria belong to aerobic heterotrophic and facultative trophic bacteria [48,49], which are able to live in soils with high organic matter content. As a result, with an increase in soil depth, the soil nutrients decreased, and the abundance of Proteobacteria decreased. Proteobacteria were not significantly correlated with other physicochemical factors, except for soil moisture (Figure 5). At the same soil depth, the difference in moisture among the three land use patterns was not significant, and the other physicochemical factors had no significant effect on the distribution of Proteobacteria because Proteobacteria contain many subgroups with different habits, which are distributed in different niches. The Proteobacteria’s structure is stable, and its anti-interference ability is strong. In the surface layer, OTUs from Proteobacteria were more prevalent in grassland (Figure 2a), and most of the OTUs were significantly correlated with TP, AK, E-Mg, EC, and T (Figure 4c). In the deep layer, Proteobacteria were not significantly different in their land use (Figure 2b). Further, the OTUs were not significantly correlated with the soil’s physicochemical properties (Figure 4d). These bacterial structural differences reflect the prominent changes in particular groups. The heat map revealed that the distribution of OTUs was strongly affected by different land uses and soil depths, and soil nutrients appear to determine the distribution and frequency of OTUs (Figure 4c,d). Among the most frequent Proteobacterial OTUs, six were classified at the genus level (Variibacter, Bradyrhizobium, Pedomicrobium, Steroidobacter, Sphingomonas, Lysobacter). Interestingly, some taxonomically related OTUs from land use, such as Proteobacteria-related OTUs (28 and 124, classified as part of the Sphingomonas group), exhibited similar occurrence patterns. Some recent advances have demonstrated that Sphingomonas have unique abilities in degrading refractory contaminants, serving as bacterial antagonists to phytopathogenic fungi [50]. The genus Bradyrhizobium (OTU 4) is one of several genera of nitrogen-fixing bacteria, which play an important role in agricultural productivity and global nitrogen cycling. Previous reports have found that Bradyrhizobium are dominant in forest soil [51]. However, in our study, Bradyrhizobium were found to be dominant in shrubland and woodland environments.
Verrucomicrobia were the most widely distributed and diverse phylum of bacteria in soil habitats [49]. The proportion of Verrucomicrobia in grassland was significantly lower than that in woodland and shrubland at different soil depths. Verrucomicrobia were distributed in various soil and aquatic habitats, including environments with poor nutrition, eutrophication, extreme pollution, and man-made habitats [52,53]. Verrucomicrobia may be more adaptable to oligotrophic environments in soils [54], which are widely distributed in rhizosphere and aggregate soils. They can use a variety of carbon compounds and may be closely related to the carbon cycle [55]. Our study found that Verrucomicrobia were significantly positively correlated with TN and SOC (Figure 5), and significantly negatively correlated with TP, AK, EC, and E-Mg. Because the grassland had more soil nutrients, the proportion of Verrucomicrobia in grassland was significantly smaller than that in woodland and shrubland. Moreover, it was found that there was a striking increase in verrucomicrobial abundances at different soil depths (Figure 1). The slow growth rate of Verrucomicrobia may allow them to exploit the sparse resources in subsurface soil [56]. In different soil layers, OTUs from Verrucomicrobia were enriched in the shrubland (Figure 2a,b). Moreover, most OTUs from Verrucomicrobia were significantly correlated with TN and SOC (Figure 4c,d). Among the most frequent Verrucomicrobia OTUs, only one (OTU 21) was classified at the genus level (Candidatus Xiphinematobacter). Previous studies exposed the difficulty of classifying Verrucomicrobia because the portion of culturable bacteria within the Verrucomicrobia is quite low [57].
Acidobacteria were widespread in all types of soils with high richness [58]. In the A layer, Acidobacteria were significantly positively correlated with EC and E-Mg (Figure 5). The proportion of Acidobacteria in grassland was significantly higher than that in woodland and shrubland. In the B layer, Acidobacteria had a significantly negative correlation with TP (Figure 5), though the proportion of Acidobacteria in the three land use patterns was not significantly different. Moreover, the dominant genera from Acidobacteria were the Blastocatella (OTU 37) and Candidatus Solibacter OTUs (21, 30, and 238). In the surface layer, OTUs from Acidobacteria were prevalent in the three land uses, reflecting the high adaptability of this group. In the deep layer, OTUs from Acidobacteria were more frequent in grassland and shrubland. In addition, most OTUs were significantly correlated with TP, AK, E-Mg, EC, and moisture. Accordingly, each abundant Acidobacteria OTU prevailed in a different niche, reflecting the high adaptability of this group. Acidobacteria are acidophilic bacteria, which can decompose animal and plant residues and cellulose to form organic carbon [59]. Consequently, they are more suitable for living in a low organic carbon environment. Our study found that Acidobacteria were negatively correlated with TOC. Moreover, there was more litter on the soil’s surface, making the surface more suitable for the survival of Acidobacteria.

5. Conclusions

Proteobacteria, Verrucomicrobia, and Acidobacteria were the dominant bacteria, but the abundance of bacterial communities was different. Verrucomicrobia were significantly positively correlated with TN and SOC and significantly negatively correlated with TP, AK, EC, and E-Mg. Acidobacteria were significantly correlated with EC, E-Mg, and TP. Different land use patterns have significant effects on the soil’s physicochemical properties. TP, AK, E-Mg, and EC in the grassland were significantly higher than those in the woodland and shrubland. TN and SOC in the woodland were higher than those in the shrubland and grassland. The soil physicochemical properties decreased with increasing soil depth, and SOC, TN, and AK were the most important physicochemical properties affecting the composition of the bacterial community. The diversity of bacterial communities in the three land use types decreased with an increase of soil depth, which was consistent with the trend of soil physicochemical properties. Different land use patterns and soil depths have significant effects on bacterial communities. Land use shapes strong effects in the soil’s bacterial community structure on the soil’s surface layer, and the effects of land use attenuation decrease with soil depth. Compared with deep soil, surface soil contains more nutrients that are more suitable for the growth and reproduction of bacteria. Our study provides a basis for understanding the influence of land use patterns and soil depths on the bacterial community structure in the karst graben basin of Yunnan province, China. Selecting a suitable land use according to the soil bacterial community structure of karst graben basins has great significance. We should also consider the impacts of soil bacterial communities on land use at different soil depths.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/11/1/51/s1: Table S1: Abundance of dominant phyla in the three land use types. Table S2: Spearman’s correlations showing the relationship between soil physicochemical parameters and dominant OTUs in the A layer. Table S3: Spearman’s correlations showing the relationship between the soil physicochemical parameters and dominant OTUs in the B layer. Figure S1: Heat map representing the relationship between the soil physicochemical parameters and the most abundant OTUs with abundances >0.5% of the three land uses in different soil layers: woodland in the surface layer (a), shrubland in the surface layer (b), grassland in the surface layer (c), woodland in the deep layer (d), shrubland in the deep layer (e), and grassland in the deep layer (f).

Author Contributions

J.Q. organized the available literature and data, and subsequently developed the original draft of the manuscript. G.L., Y.L., and H.W. collected the soil samples. Q.L. planned and designed the research. Q.L. and J.C. developed the original concept for the project, co-authored the proposal to fund the research, and edited all subsequent drafts of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Project supported by the National Key Research and Development Program of China “Ecological, environmental and geological differentiation of rocky desertification and its driving mechanism in the karst graben basin” (No. 2016YFC0502501).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yao, L.S. The formation mechanism and model of faulted karst basins in Yunnan province. Carsol. Sin. 1984, 3, 48–55. (In Chinese) [Google Scholar]
  2. Wang, Y.; Zhang, H.; Zhang, G.; Wang, B.; Peng, S.H.; He, R.S.; Zhou, C.Q. Zoning of environmental geology and functions in karst fault-depression basins. Carsol. Sin. 2017, 36, 283–295. [Google Scholar] [CrossRef]
  3. Deng, L.; Zhou, P.; Guan, S.; Li, R. Effects of the grain-for-green program on soil erosion in China. Int. J. Sediment. Res. 2012, 27, 120–127. [Google Scholar] [CrossRef]
  4. Yang, J. Yunnan Natural Forest Resources Conservation Project Phase II Started in an All-round Way. Yunnan For. 2011, 32, 24. (In Chinese) [Google Scholar]
  5. Bai, C.L. Deepening Reform and Accelerating Development Creating a New Situation of Natural Forest Protection Project Construction in Yunnan Province. Yunnan For. 2007, 28, 4–7. (In Chinese) [Google Scholar]
  6. Xu, H.J.; Wang, X.H.; Li, H.; Yao, H.Y.; Su, J.Q.; Zhu, Y.G. Biochar impacts soil microbial community composition and nitrogen cycling in an acidic soil planted with rape. Environ. Sci. Technol. 2014, 48, 9391–9399. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Du, B.H.; Jin, Z.G.; Li, Z.H.; Song, H.N.; Ding, Y.Q. Analysis of bacterial communities in rhizosphere soil of healthy and diseased cotton (Gossypium sp.) at different plant growth stages. Plant. Soil 2011, 339, 447–455. [Google Scholar] [CrossRef]
  8. McCann, K.S. The diversity-stability debate. Nature 2000, 405, 228–233. [Google Scholar] [CrossRef]
  9. Suleiman, A.K.A.; Manoeli, L.; Boldo, J.T.; Pereira, M.G.; Roesch, L.F.W. Shifts in soil bacterial community after eight years of land-use change. Syst. Appl. Microbiol. 2013, 36, 137–144. [Google Scholar] [CrossRef]
  10. Song, M.; Zou, D.S.; Du, H.; Peng, W.X.; Zeng, F.P.; Tan, Q.J.; Fan, F.J. Characteristics of soil microbial populations in depressions between karst hills under different land use patterns. Chin. J. Appl. Ecol. 2013, 24, 2471–2478. [Google Scholar] [CrossRef]
  11. Edarson, D.C.J.; Terence, L.M.; James, M.T.; Fatima, M.D.S.M. Changes in land use alter the structure of bacterial communities in Western Amazon soils. ISME J. 2009, 3, 1004–1011. [Google Scholar] [CrossRef]
  12. Bremner, J.M.; Jenkinson, D.S. Determination of organic carbon in soil. I. Oxidation by dichromate of organic matter in soil and plant materials. Eur. J. Soil Sci. 2010, 11, 394–402. [Google Scholar] [CrossRef]
  13. Parkinson, J.A.; Allen, S.E. A wet oxidation procedure suitable for determination of nitrogen and mineral nutrients in biological material. Commun. Soil Sci. Plant Anal. 1975, 6, 1–11. [Google Scholar] [CrossRef]
  14. Carson, P.L. Recommended potassium test. In Recommended Chemical Soil Test Procedures for the North Central Region; Dahnke, W.C., Ed.; North Dakota Agricultural Experiment Station: Fargo, ND, USA, 1980; Volume 499, pp. 17–18. [Google Scholar]
  15. Wilke, B.M. Determination of chemical and physical soil properties. In Monitoring and Assessing Soil Bioremediation; Springer: Heidelberg/Berlin, Germany, 2005. [Google Scholar] [CrossRef]
  16. Li, Q.; Hu, Q.; Zhang, C.; Müller, W.E.G.; Schröder, H.C.; Li, Z. The effect of toxicity of heavy metals contained in tailing sands on the organic carbon metabolic activity of soil microorganisms from different land use types in the karst region. Environ. Earth Sci. 2015, 74, 6747–6756. [Google Scholar] [CrossRef]
  17. Chan, C.O.; Yang, X.D.; Fu, Y.; Feng, Z.L.; Sha, L.Q.; Peter, C.; Zou, X.M. 16S rRNA gene analyses of bacterial community structures in the soils of evergreen broad-leaved forests in south-west China. FEMS Microbiol. Ecol. 2006, 58, 247–259. [Google Scholar] [CrossRef] [Green Version]
  18. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  19. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  20. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucl. Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  21. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [Green Version]
  22. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  23. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Goodrich, J.K.; Gordon, J.I.; Huttley, G.A.; et al. QIIME allows analysis of highthroughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zuur, A.F.; Ieno, E.N.; Smith, G.M. Analysing Ecological Data; Springer: New York, NY, USA, 2007. [Google Scholar]
  25. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014; Available online: http://www.R-project.org/ (accessed on 15 July 2019).
  26. Anderson, M.J.; Willis, T.J. Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology 2003, 84, 511–525. [Google Scholar] [CrossRef]
  27. Bastian, M.; Heymann, S.; Gephi, M.J. An open source software for exploring and manipulating networks. In Proceedings of the Third International Conference on Weblogs and Social Media, ICWSM 2009, San Jose, CA, USA, 17–20 May 2009. [Google Scholar] [CrossRef]
  28. Oberson, A.; Friesen, D.K.; Rao, I.M.; Bühler, S.; Frossard, E. Phosphorus Transformations in an Oxisol under contrasting land-use systems: The role of the soil microbial biomass. Plant Soil 2001, 237, 197–210. [Google Scholar] [CrossRef]
  29. Batlle-Aguilar, J.; Brovelli, A.; Porporato, A.; Barry, D.A. Modelling soil carbon and nitrogen cycles during land use change-A review. Agron. Sustain. Dev. 2011, 31, 251–274. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, Y.P.; Lu, M.X.; Zhang, X.W.; Sun, Q.B.; Liu, R.L.; Lian, B. Shift of the microbial communities from exposed sandstone rocks to forest soils during pedogenesis. Int. Biodeter. Biodegr. 2019, 140, 21–28. [Google Scholar] [CrossRef]
  31. Huang, X.F.; Zhou, Y.C.; Zhang, Z.M. Distribution characteristics of soil organic carbon under different land use in a karst rocky desertification area. J. Soil Water Conserv. 2017, 31, 215–221. [Google Scholar] [CrossRef]
  32. Tong, J.H.; Hu, Y.C.; Du, Z.L.; Zou, Y.Q.; Li, Y.Y. Effects of land use change on soil organic carbon and total nitrogen storage in karst immigration regions of Guanxi Province, China. Chin. J. Appl. Ecol. 2018, 29, 2890–2896. [Google Scholar] [CrossRef]
  33. Li, L.; Qin, F.C.; Jiang, L.N.; Yao, X.L. Vertical distribution of soil organic carbon content and its influenceing factors in Aaohan, Chifeng. Acta Ecol. Sin. 2019, 39, 345–354. [Google Scholar] [CrossRef]
  34. Chaplot, V.; Bouahom, B.; Valentin, C. Soil organic carbon stocks in Laos: Spatial variations and controlling factors. Glob. Chang. Biol. 2010, 16, 1380–1393. [Google Scholar] [CrossRef]
  35. Kuang, W.N.; Qian, J.Q.; Ma, Q.; Liu, Z.M. Vertical distribution of soil organic carbon content and its relation to root distribution in five desert shrub communities. Chin. J. Ecol. 2016, 35, 275–281. [Google Scholar] [CrossRef]
  36. Wang, D.; Geng, Z.C.; She, D.; He, W.X.; Hou, L. Soil organic carbon storage and vertical distribution of carbon and nitrogen across different forest types in the Qinling Mountains. Acta Ecol. Sin. 2015, 35, 5421–5429. [Google Scholar] [CrossRef] [Green Version]
  37. Daniel, C.S.; Matthew, G.B.; James, M.B.; Linda, L.K. Plant community richness and microbial interactions structure bacterial communities in soil. Ecology 2015, 96, 134–142. [Google Scholar] [CrossRef] [Green Version]
  38. Kinkel, L.L.; Bakker, M.G.; Schlatter, D.C. A coevolutionary framework for managing disease-suppressive soils. Annu. Rev. Phytopathol. 2011, 49, 47–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Zhang, C.; Li, J.; Wang, J.; Liua, G.B.; Wang, G.L.; Guo, L.; Peng, S.Z. Decreased temporary turnover of bacterial communities along soil depth gradient during a 35-year grazing exclusion period in a semiarid grassland. Geoderma 2019, 351, 49–58. [Google Scholar] [CrossRef]
  40. He, S.B.; Guo, L.X.; Niu, M.Y.; Miao, F.H.; Jiao, S.; Hu, T.M.; Long, M.X. Ecological diversity and cooccurrence patterns of bacterial community through soil profile in response to long-term switchgrass cultivation. Sci. Rep. 2017, 3608, 7. [Google Scholar] [CrossRef]
  41. Cheng, J.M.; Guanghua Jing, G.H.; Wei, L.; Jing, Z.B. Long-term grazing exclusion effects on vegetation characteristics, soil properties and bacterial communities in the semiarid grasslands of China. Ecol. Eng. 2016, 97, 170–178. [Google Scholar] [CrossRef]
  42. Wang, S.X.; Wang, X.A.; Guo, H. Change patterns of β-diversity in the succession process of plant communities on Loess Plateau of Northwest China. Chin. J. Ecol. 2013, 32, 1135–1140. [Google Scholar] [CrossRef]
  43. Tian, Q.; Taniguchi, T.; Shi, W.Y.; Li, G.Q.; Yamanaka, N.; Du, S. Land-use types and soil chemical properties influence soil microbial communities in the semiarid Loess Plateau region in China. Sci. Rep. 2017, 7, 45289. [Google Scholar] [CrossRef] [Green Version]
  44. Bru, D.; Ramette, A.; Saby, N.P.A.; Dequiedt, S.; Ranjard, L.; Jolivet, C.; Arrouays, D.; Philippot, L. Determinants of the distribution of nitrogen-cycling microbial communities at the landscape scale. ISME J. 2011, 5, 532–542. [Google Scholar] [CrossRef] [Green Version]
  45. Adria, L.F.; Craig, C.S.; Donald, L.W.; Christopher, S.; Trevor, J.G.; Michael, J.S. Associations between soil bacterial community structure and nutrient cycling functions in long-term organic farm soils following cover crop and organic fertilizer amendment. Sci. Total Environ. 2016, 566, 949–959. [Google Scholar] [CrossRef] [Green Version]
  46. Manuel, D.B.; Angela, M.O.; Tess, E.B.; Alberto, B.G.; David, J.E.; Richard, D.B.; Fernando, T.M.; Brajesh, K.S.; Noah, F. A global atlas of the dominant bacteria found in soil. Science 2018, 359, 320–325. [Google Scholar] [CrossRef] [Green Version]
  47. Zhou, J.; Guan, D.W.; Zhou, B.K.; Zhao, B.S.; Ma, M.C.; Qin, J.; Jiang, X.; Chen, S.F.; Cao, F.M.; Shen, D.L.; et al. Influence of 34-years of fertilization on bacterial communities in an intensively cultivated black soil in northeast China. Soil Biol. Biochem. 2015, 90, 42–51. [Google Scholar] [CrossRef]
  48. Griffiths, B.S.; Philippot, L. Insights into the resistance and resilience of the soil microbial community. FEMS Microbiol. Rev. 2013, 37, 112–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Anna, K.; Jorge, L.M.R.; Eiko, E.K.; Patrick, S.G.C.; Johannes, A.V.V.; George, A.K. Phylogenetic and metagenomic analysis of Verrucomicrobia in former agricultural grassland soil. FEMS Microbiol. Ecol. 2010, 71, 23–33. [Google Scholar] [CrossRef]
  50. David, C.W.; Susan, D.S.; David, B.R. The genus Sphingomonas: Physiology and ecology. Curr. Opin. Biotechnol. 1996, 7, 301–306. [Google Scholar] [CrossRef]
  51. David, V.; Kendra, R.M.; Erick, C.; Cameron, R.S.; Steven, J.H.; William, W.M. Non-symbiotic Bradyrhizobium ecotypes dominateNorth American forest soils. ISME J. 2015, 9, 2435–2441. [Google Scholar] [CrossRef] [Green Version]
  52. Heinz, S.; Cheryl, J.; Jamest, S. The phylum Verrucomicrobia: A phylogenetically heterogeneous bacterial group. Prokaryotes 2006, 7, 881–896. [Google Scholar] [CrossRef]
  53. Wagner, M.; Horn, M. The Planctomycetes, Verrucomicrobia, Chlamydiae and sister phyla comprise a superphylum with biotechnological and medical relevance. Curr. Opin. Biotechnol. 2006, 17, 241–249. [Google Scholar] [CrossRef]
  54. Noll, M.; Diethart, M.; Frenzel, P.; Manigee, D.; Liesack, W. Succession of bacterial community structure and diversity in a paddy soil oxygen gradient. Environ. Microbiol. 2005, 7, 382–395. [Google Scholar] [CrossRef]
  55. Fierer, N.; Ladau, J.; Clementei, J.C.; Leff, J.W.; Owens, S.M.; Katherine, S.P.; Knight, R.; Gilbert, J.A.; McCulley, R.L. Reconstructing the microbial diversity and function of pre-agricultural tallgrass prairie soils in the United States. Science 2013, 342, 621–624. [Google Scholar] [CrossRef] [Green Version]
  56. Janssen, P.H.; Yates, P.S.; Grinton, B.E.; Taylor, P.M.; Sait, M. Improved culturability of soil bacteria and isolation in pure culture of novel members of the divisions Acidobacteria, Actinobacteria, Proteobacteria, and Verrucomicrobia. Appl. Environ. Microb. 2002, 68, 2391–2396. [Google Scholar] [CrossRef] [Green Version]
  57. Stefan, S.; Boyke, B.; Cathrin, S.; Peter, S.; Manfred, R.; Brian, J.T.; Hans, P.K. Characterization of the first cultured representative of Verrucomicrobia subdivision 5 indicates the proposal of a novel phylum. ISME J. 2016, 10, 2801–2816. [Google Scholar] [CrossRef] [Green Version]
  58. Wang, P.; Chen, B.; Zhang, H. High throughput sequencing analysis of bacterial communities in soils of a typical Poyang Lake wetland. Acta Ecol. Sin. 2017, 37, 1650–1658. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, Y.G.; Cong, J.; Lu, H.; Li, G.L.; Qu, Y.Y.; Su, X.J.; Zhou, J.Z.; Li, D.Q. Community structure and elevational diversity patterns of soil Acidobacteria. J. Environ. Sci. 2014, 26, 171–1724. [Google Scholar] [CrossRef]
Figure 1. Comparison of the average quantitative contribution of the sequences affiliated with different bacterial phyla from the A layer (a) and B layer (b). Sequences not classified under any known phylum are included as unassigned bacteria. In each soil sample, the bacterial phylum with relative frequency of less than 1% is included as others.
Figure 1. Comparison of the average quantitative contribution of the sequences affiliated with different bacterial phyla from the A layer (a) and B layer (b). Sequences not classified under any known phylum are included as unassigned bacteria. In each soil sample, the bacterial phylum with relative frequency of less than 1% is included as others.
Forests 11 00051 g001
Figure 2. Heat map illustrating the mean relative frequency of the 47 most abundant operational taxonomic units (OTUs) with abundances >0.5% in the A layer with different land uses (a). Heat map illustrating the mean relative frequency of the 45 most abundant OTUs with abundances >0.5% in the B layer with different land uses (b). Taxonomic assignment of the OTUs is provided at the lowest level of classification possible according to the SILVA 123 database (D1: phylum, D2: class, D3: order, D4: family, and D5: genus).
Figure 2. Heat map illustrating the mean relative frequency of the 47 most abundant operational taxonomic units (OTUs) with abundances >0.5% in the A layer with different land uses (a). Heat map illustrating the mean relative frequency of the 45 most abundant OTUs with abundances >0.5% in the B layer with different land uses (b). Taxonomic assignment of the OTUs is provided at the lowest level of classification possible according to the SILVA 123 database (D1: phylum, D2: class, D3: order, D4: family, and D5: genus).
Forests 11 00051 g002
Figure 3. Principal coordinate analysis (PCoA) plots of soil microbial community structures based on the Bray–Curtis and weighted UniFrac results under different land use types in the Yunnan karst graben basin. Bray–Curtis (a); weighted UniFrac (b).
Figure 3. Principal coordinate analysis (PCoA) plots of soil microbial community structures based on the Bray–Curtis and weighted UniFrac results under different land use types in the Yunnan karst graben basin. Bray–Curtis (a); weighted UniFrac (b).
Forests 11 00051 g003
Figure 4. Distance-based redundancy analysis (db-RDA) plots revealing the relationship between the dominant phyla (mean relative frequency >1%) and soil physicochemical parameters in the surface layer (a) and deep layer (b). Heat map representing the relationship between the soil physicochemical parameters and the most abundant OTUs with abundances >0.5% in the surface layer (c) and deep layer (d). The networks revealing the co-occurring bacterial OTUs and soil properties in the surface layer (e) and deep layer (f). The co-occurring networks are colored at the phylum level, and the size of each node is proportional to the number of connections.
Figure 4. Distance-based redundancy analysis (db-RDA) plots revealing the relationship between the dominant phyla (mean relative frequency >1%) and soil physicochemical parameters in the surface layer (a) and deep layer (b). Heat map representing the relationship between the soil physicochemical parameters and the most abundant OTUs with abundances >0.5% in the surface layer (c) and deep layer (d). The networks revealing the co-occurring bacterial OTUs and soil properties in the surface layer (e) and deep layer (f). The co-occurring networks are colored at the phylum level, and the size of each node is proportional to the number of connections.
Forests 11 00051 g004
Figure 5. Spearman’s correlations showing the relationship between the soil physicochemical parameters and dominant phyla in the three land use types. Significance levels are denoted as follows: p < 0.05 (*) and p < 0.01 (**).
Figure 5. Spearman’s correlations showing the relationship between the soil physicochemical parameters and dominant phyla in the three land use types. Significance levels are denoted as follows: p < 0.05 (*) and p < 0.01 (**).
Forests 11 00051 g005
Table 1. The soil physicochemical properties in different land use patterns.
Table 1. The soil physicochemical properties in different land use patterns.
NameLand Use PatternTPTNSOCAKpHMoistureTECE-CaE-Mg
(g/kg)(g/kg)(g/kg)(g/kg)(%)(°C)(ms.m−1)(cmol/kg)(cmol/kg)
WLAwoodland0.828 ± 0.020 ab5.09 ± 0.24 a61.00 ± 1.96 a128.47 ± 5.78 b6.21 ± 0.11 b37.47 ± 5.04 a8.60 ± 0.21 c68.33 ± 4.10 b5.14 ± 1.15 a26.82 ± 0.65 c
SLAShrubland0.734 ± 0.025 b4.11 ± 0.16 b47.00 ± 2.5 b124.03 ± 14.05 b6.72 ± 0.09 a42.85 ± 2.56 a11.60 ± 1.12 b64.33 ± 2.03 b6.91 ± 0.09 a36.18 ± 3.48 b
GLAGrassland0.941 ± 0.059 a2.97 ± 0.14 c33.00 ± 1.47 c262.17 ± 44.37 a6.68 ± 0.08 a40.00 ± 1.42 a16.53 ± 0.74 a84.67 ± 2.03 a6.87 ± 0.08 a51.56 ± 2.32 a
WLBwoodland0.640 ± 0.037 a3.16 ± 0.49 a35.79 ± 5.32 a40.40 ± 2.14 b6.46 ± 0.10 a36.53 ± 2.98 a9.37 ± 0.09 c70.33 ± 6.64 ab6.64 ± 0.11 a29.21 ± 0.28 c
SLBShrubland0.595 ± 0.025 a2.32 ± 0.12 a24.53 ± 2.13 ab38.40 ± 3.98 b6.66 ± 0.08 a33.07 ± 1.91 a10.97 ± 0.33 b56.00 ± 3.00 b6.85 ± 0.08 a34.20 ± 1.02 b
GLBGrassland0.762 ± 0.084 a2.34 ± 0.37 a20.46 ± 3.71 b127.27 ± 15.21 a6.55 ± 0.03 a32.70 ± 1.16 a12.77 ± 0.34 a79.33 ± 5.70 a6.73 ± 0.03 a39.82 ± 1.05 a
Note: A = 0–10 cm (A layer), B = 10–20 cm (B layer); data are the means ± standard error (means ± SE); different lowercase letters (a, b and c) in the same column represent a significant difference from the different sample points in the same soil layer (p < 0.05). TP: total phosphorous, TN: total nitrogen, SOC: soil organic carbon, AK: available potassium, T: temperature, EC: electrical conductivity, E-Ca: exchangeable calcium, E-Mg: exchangeable magnesium.
Table 2. Alpha diversities of the soil bacterial communities with different land use patterns.
Table 2. Alpha diversities of the soil bacterial communities with different land use patterns.
NameLand Use PatternChao1ShannonSimpsonObserved OTUs
WLAWoodland1656 ± 69 b8.01 ± 0.23 b0.985 ± 0.004 a1162 ± 45 b
SLAShrubland1640 ± 70 b7.69 ± 0.17 b0.977 ± 0.004 b1103 ± 48 b
GLAGrassland1935 ± 75 a8.85 ± 0.10 a0.993 ± 0.001 a1409 ± 25 a
WLBWoodland1582 ± 89 a7.80 ± 0.31 a0.981 ± 0.005 b1115 ± 67 a
SLBShrubland1565 ± 62 a7.70 ± 0.07 a0.982 ± 0.002 a1030 ± 24 a
GLBGrassland1693 ± 119 a8.54 ± 0.37 a0.993 ± 0.003 a1212 ± 106 a
Note: A = 0–10 cm (A layer), B = 10–20 cm (B layer); data are the means ± standard error (means ± SE). Different lowercase letters (a and b) in the same column represent a significant difference from the different sample points in the same soil layer (p < 0.05).

Share and Cite

MDPI and ACS Style

Qiu, J.; Cao, J.; Lan, G.; Liang, Y.; Wang, H.; Li, Q. The Influence of Land Use Patterns on Soil Bacterial Community Structure in the Karst Graben Basin of Yunnan Province, China. Forests 2020, 11, 51. https://doi.org/10.3390/f11010051

AMA Style

Qiu J, Cao J, Lan G, Liang Y, Wang H, Li Q. The Influence of Land Use Patterns on Soil Bacterial Community Structure in the Karst Graben Basin of Yunnan Province, China. Forests. 2020; 11(1):51. https://doi.org/10.3390/f11010051

Chicago/Turabian Style

Qiu, Jiangmei, Jianhua Cao, Gaoyong Lan, Yueming Liang, Hua Wang, and Qiang Li. 2020. "The Influence of Land Use Patterns on Soil Bacterial Community Structure in the Karst Graben Basin of Yunnan Province, China" Forests 11, no. 1: 51. https://doi.org/10.3390/f11010051

APA Style

Qiu, J., Cao, J., Lan, G., Liang, Y., Wang, H., & Li, Q. (2020). The Influence of Land Use Patterns on Soil Bacterial Community Structure in the Karst Graben Basin of Yunnan Province, China. Forests, 11(1), 51. https://doi.org/10.3390/f11010051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop