# Soil Microbial Community Structure Evolution along Halophyte Succession in Bohai Bay Wetland

Authors: Mingyang Cong; Di Cao; Jingkuan Sun; Fuchen Shi
Journal: Journal of Chemistry
 (2014)

---

## Abstract

It is urgent to recover Bohai Bay costal wetland ecosystem because of covering a large area of severe saline-alkali soil. To explore the relationship between halophyte herbaceous succession and microbial community structure, we chose four local communities which played an important role in improving soil microenvironment. We performed phospholipid fatty acid analysis, measured soil parameters, and evaluated shifts of microbial community structure. Results showed that microbial community structure changed significantly along succession and bacteria community was dominant. Total phospholipid fatty acid content increased in different successional stages but decreased with depth, with similar variations in bacterial and fungal biomass. Soil organic carbon and especially total nitrogen were positively correlated with microbial biomass. Colonization of pioneering salt-tolerant plantsSuaeda glauca in saline-alkali bare land changed total soil microorganism content and composition. These results showed that belowground processes were strongly related with aboveground halophyte succession. Fungal/bacterial ratio, Gram-negative/Gram-positive bacteria ratio, total microbial biomass, and fungi and bacteria content could indicate the degree of succession stages in Bohai Bay wetland ecosystem. And also these findings demonstrated that microbial community biomass and composition evolved along with vegetation succession environmental variables.


---

## Body

## 1. Introduction

Vegetation, soil, and microbes depend on and restrict each other [1]. Microbes are sensitive to external environmental changes, and the most remarkable characteristic is the ability to change the community structure [2]. The role of soil microorganisms in the ecosystem is to decompose organic matter in the soil and promote formation of humus [3, 4], absorb, fix, and release nutrients, and improve and regulate nutrition in plants [5, 6]. Previous studies examining relationship between vegetation succession and microbial communities were conducted in a forest ecosystem [7], karst ecosystem [8], or volcanic desert ecosystem [9]. However, studies on examining an extreme saline-alkali wetland ecosystem are rare. Beidagang wetland is one of the most important wetlands located in Bohai Bay, which is officially included in “list of important wetlands in China.” Coastal wetlands are fragile and sensitive ecological areas with poor soil. On account of different climate and soil conditions from inland, the vegetation has unique characteristics. High salinity soil samples from four plots representing different successional stages were collected for this study. Main local halophyte in study site includedSuaeda glauca andPhragmites australis.Suaeda glauca is a kind of leafy succulent plant that can accumulate salt ions. Such a strong alkali euhalophyte can decrease the salt content of topsoil during growing season. On the contrary,Phragmites australis is a kind of salt-excluding plant which can reduce its own salt concentration through physiological structure of cell membrane in roots [10].A phospholipid fatty acid (PLFA) method was employed to examine the soil microbial community. PLFAs are important components of living microbe cell membrane and are synthesized by various microorganisms through different biochemical pathways. Thus, PLFAs can accurately express the biomass and composition of the soil microbial community [11].The aims of this study were to elucidate (1) the response of microorganism community structural composition to, along with environmental changes caused by, halophyte succession, and (2) main factors influencing soil microbial community structure along with succession. The results of this study may offer critical information for ecological restoration and soil quality improvement in severe saline-alkali wetlands.


## 2. Materials and Methods

### 2.1. Study Sites

It was located in the southeast of Tianjin near Bohai Bay, middle-east of China (117°27′E, 38°43′N) (Figure1 and Table 1). It was characterized by a warm semihumid continental monsoon climate. Annual mean temperature was approximately 12.2°C, with the coldest (January) and hottest (July) mean temperatures of 3.5°C and 26.2°C, respectively. Annual mean precipitation was 520–660 mm, and 75% was summer rainfall (July–September). The study area was covered by extreme saline-alkali soil, together with larger amount of evaporation and smaller amount of precipitation. The soil in this area was barren. Landform was a coastal floodplain with sparse vegetation.Table 1
Dominant species in plots.
Stage
Plant community
Average height (cm)
Coverage (%)
Biomass (g m−2)
Altitude (m)
I
Saline-alkali bare land
—
—
—
10
II
Suaeda glauca
16
74
1.8
±
0.4
10.3
III
Suaeda glauca + Phragmites australis
24
69
2.3
±
0.3
10.6
109
66
169
±
10
10.6
IV
Phragmites australis
136
83
192
±
7
11Figure 1
Locations of study and sample sites (SS) near the Bohai Bay. Details of the plots and abbreviation names of dominant species are listed in Table1.


### 2.2. Soil Sampling

We investigated four community types with three repetitions in June 2013. Table1 summarized the information of four successional stages, covering average height, coverage, biomass, and altitude of dominant species. The four typical communities examined were (I) saline-alkali bare land (BL), (II)Suaeda glauca community (SG), (III)Suaeda glauca +Phragmites australis community (SG + PA), and (IV)Phragmites australis community (PA). Soil samples were collected using multipoint mixed method with three vertical layers, S1 (0–10 cm), S2 (10–20 cm), and S3 (20–30 cm). We collected in two times on account of different analyses. One set was used to measure bulk density and moisture content by a cutting ring. Another set was divided into two parts: for PLFA analysis, kept at −20°C; for parameter analyses, dried in shade, homogenized, and sieved before preserved.


### 2.3. Soil Physicochemical Analyses

Soil organic carbon (SOC) and total nitrogen (TN) were measured by combustion using a macroelement analyzer (vario MACRO CN, Elementar Analysensysteme GmbH, Germany). Bulk density and moisture content were determined with samples dried in oven at 105°C for 48 h until a constant weight. Soil pH was determined with air-dried soil (soil : water, 1 : 2.5) using a glass electrode Sartorius PB-10 pH meter (Sartorius, Germany). Salinity was determined by drying-weighing method. Filtrate (soil : water, 1 : 5) was placed in oven at 105°C to constant weight.


### 2.4. Phospholipid Fatty Acid Analysis (PLFA)

PLFAs are found in living microbial cell membranes and can be used as markers of biomass and a fingerprint of the microorganism community structure [2, 12]. PLFA affords a quick and reliable method for estimating microbial biomass and interpreting variations in community structures [13]. Ester-linked method was used for fatty acid methyl ester (FAME) profile analyses [14]. The specific procedure to extracting phospholipids was divided into six steps: (1) 15 mL of 0.2 M KOH in methanol was added to a 35 mL centrifuge tube containing 3 g of soil. (2) The mixture was incubated at 37°C for 1 h, and during the process ester-linked fatty acids were released and methylated. Samples were vortexed every 10 min. (3) 3 mL of 1.0 M acetic acid and 10 mL of hexane were added to the contents in the tube. (4) FAMEs were partitioned into an organic phase by centrifugation at 4800 r for 10 min. (5) The hexane layer was transferred to a clean glass test tube and the hexane was evaporated under a stream of N2. (6) FAMEs were dissolved in 0.5 mL of 1 : 1 hexane : methyltert-butyl ether and transferred to a GC vial for analysis. FAMEs extracted were detected using a gas chromatograph-mass spectrometer (7890A GC/5975C MSD; Agilent Technologies Co., USA) fitted with an HP-5MS quartz capillary chromatographic column (0.25 mm × 30 m × 0.25 μm), with an injector temperature at 250°C. Internal standard for quantification was C19:0 (Sigma-Aldrich Co., USA), and external standard for qualification was GLC NESTLE 37 mix (Nu-Chek).Bardgett et al. [15] previously demonstrated that the content and composition of PLFA in the soil reflected the biomass and structure of soil microbial communities. Results were analyzed by the MIDI system (Microbial ID, Inc.). Standard nomenclature was employed to describe FAMEs. We used the PLFA nomenclature pioneered by Frostegard et al. [12] and the total content of PLFA to indicate the total microbial biomass. Twelve fatty acids (15:0, i15:0, a15:0, i16:0, 16:1ω7, 16:1ω9, 17:0, i17:0, a17:0, 18:1ω7, i19:0, and cy19:0) were chosen as indicators of bacterial PLFA [16]. Fungal PLFA was represented by 18:2ω6 [12]; anaerobic bacterial PLFA was represented by cy17:0 and cy19:0 [17]; aerobic bacterial PLFA was represented by 16:1ω7, 16:1ω7t, and 18:1ω7; Gram-negative bacterial PLFA and monounsaturated fatty acids were represented by cy17:0 and cy19:0 [18]; Gram-positive bacterial PLFAs were represented by five fatty acids (i15:0, a15:0, i16:0, i17:0, and a17:0) [19].


### 2.5. Statistical Analyses

Statistical analyses were conducted with SPSS 17.0 software for Windows. One-way analysis of variance (ANOVA) was used to test for significant differences among the four successional stages and among the three soil layers in each stage. We subdivided the data into three sets: vegetation data (average height, coverage, and biomass), soil data (SOC, TN, C/N ratio, pH, soil moisture content, bulk density, and salinity), and microbial community diversity. We used Duncan’s multiple range tests and defined significance asP<0.05 in both the vertical and horizontal gradients. We used Spearman’s rank correlations to examine the association between the microbial PLFA content and the physicochemical characteristics in soil (P<0.01). Redundancy discriminate analysis (RDA), a linear canonical community ordination method, was used to explore the vital soil factors leading to changes in the microbial community structure with Canoco for Windows 4.5. The Monte Carlo permutation test (number of permutations, 999) was conducted with P values to assess the significance of soil variables in accounting for the shifts in microbial community parameters.


## 2.1. Study Sites

It was located in the southeast of Tianjin near Bohai Bay, middle-east of China (117°27′E, 38°43′N) (Figure1 and Table 1). It was characterized by a warm semihumid continental monsoon climate. Annual mean temperature was approximately 12.2°C, with the coldest (January) and hottest (July) mean temperatures of 3.5°C and 26.2°C, respectively. Annual mean precipitation was 520–660 mm, and 75% was summer rainfall (July–September). The study area was covered by extreme saline-alkali soil, together with larger amount of evaporation and smaller amount of precipitation. The soil in this area was barren. Landform was a coastal floodplain with sparse vegetation.Table 1
Dominant species in plots.
Stage
Plant community
Average height (cm)
Coverage (%)
Biomass (g m−2)
Altitude (m)
I
Saline-alkali bare land
—
—
—
10
II
Suaeda glauca
16
74
1.8
±
0.4
10.3
III
Suaeda glauca + Phragmites australis
24
69
2.3
±
0.3
10.6
109
66
169
±
10
10.6
IV
Phragmites australis
136
83
192
±
7
11Figure 1
Locations of study and sample sites (SS) near the Bohai Bay. Details of the plots and abbreviation names of dominant species are listed in Table1.


## 2.2. Soil Sampling

We investigated four community types with three repetitions in June 2013. Table1 summarized the information of four successional stages, covering average height, coverage, biomass, and altitude of dominant species. The four typical communities examined were (I) saline-alkali bare land (BL), (II)Suaeda glauca community (SG), (III)Suaeda glauca +Phragmites australis community (SG + PA), and (IV)Phragmites australis community (PA). Soil samples were collected using multipoint mixed method with three vertical layers, S1 (0–10 cm), S2 (10–20 cm), and S3 (20–30 cm). We collected in two times on account of different analyses. One set was used to measure bulk density and moisture content by a cutting ring. Another set was divided into two parts: for PLFA analysis, kept at −20°C; for parameter analyses, dried in shade, homogenized, and sieved before preserved.


## 2.3. Soil Physicochemical Analyses

Soil organic carbon (SOC) and total nitrogen (TN) were measured by combustion using a macroelement analyzer (vario MACRO CN, Elementar Analysensysteme GmbH, Germany). Bulk density and moisture content were determined with samples dried in oven at 105°C for 48 h until a constant weight. Soil pH was determined with air-dried soil (soil : water, 1 : 2.5) using a glass electrode Sartorius PB-10 pH meter (Sartorius, Germany). Salinity was determined by drying-weighing method. Filtrate (soil : water, 1 : 5) was placed in oven at 105°C to constant weight.


## 2.4. Phospholipid Fatty Acid Analysis (PLFA)

PLFAs are found in living microbial cell membranes and can be used as markers of biomass and a fingerprint of the microorganism community structure [2, 12]. PLFA affords a quick and reliable method for estimating microbial biomass and interpreting variations in community structures [13]. Ester-linked method was used for fatty acid methyl ester (FAME) profile analyses [14]. The specific procedure to extracting phospholipids was divided into six steps: (1) 15 mL of 0.2 M KOH in methanol was added to a 35 mL centrifuge tube containing 3 g of soil. (2) The mixture was incubated at 37°C for 1 h, and during the process ester-linked fatty acids were released and methylated. Samples were vortexed every 10 min. (3) 3 mL of 1.0 M acetic acid and 10 mL of hexane were added to the contents in the tube. (4) FAMEs were partitioned into an organic phase by centrifugation at 4800 r for 10 min. (5) The hexane layer was transferred to a clean glass test tube and the hexane was evaporated under a stream of N2. (6) FAMEs were dissolved in 0.5 mL of 1 : 1 hexane : methyltert-butyl ether and transferred to a GC vial for analysis. FAMEs extracted were detected using a gas chromatograph-mass spectrometer (7890A GC/5975C MSD; Agilent Technologies Co., USA) fitted with an HP-5MS quartz capillary chromatographic column (0.25 mm × 30 m × 0.25 μm), with an injector temperature at 250°C. Internal standard for quantification was C19:0 (Sigma-Aldrich Co., USA), and external standard for qualification was GLC NESTLE 37 mix (Nu-Chek).Bardgett et al. [15] previously demonstrated that the content and composition of PLFA in the soil reflected the biomass and structure of soil microbial communities. Results were analyzed by the MIDI system (Microbial ID, Inc.). Standard nomenclature was employed to describe FAMEs. We used the PLFA nomenclature pioneered by Frostegard et al. [12] and the total content of PLFA to indicate the total microbial biomass. Twelve fatty acids (15:0, i15:0, a15:0, i16:0, 16:1ω7, 16:1ω9, 17:0, i17:0, a17:0, 18:1ω7, i19:0, and cy19:0) were chosen as indicators of bacterial PLFA [16]. Fungal PLFA was represented by 18:2ω6 [12]; anaerobic bacterial PLFA was represented by cy17:0 and cy19:0 [17]; aerobic bacterial PLFA was represented by 16:1ω7, 16:1ω7t, and 18:1ω7; Gram-negative bacterial PLFA and monounsaturated fatty acids were represented by cy17:0 and cy19:0 [18]; Gram-positive bacterial PLFAs were represented by five fatty acids (i15:0, a15:0, i16:0, i17:0, and a17:0) [19].


## 2.5. Statistical Analyses

Statistical analyses were conducted with SPSS 17.0 software for Windows. One-way analysis of variance (ANOVA) was used to test for significant differences among the four successional stages and among the three soil layers in each stage. We subdivided the data into three sets: vegetation data (average height, coverage, and biomass), soil data (SOC, TN, C/N ratio, pH, soil moisture content, bulk density, and salinity), and microbial community diversity. We used Duncan’s multiple range tests and defined significance asP<0.05 in both the vertical and horizontal gradients. We used Spearman’s rank correlations to examine the association between the microbial PLFA content and the physicochemical characteristics in soil (P<0.01). Redundancy discriminate analysis (RDA), a linear canonical community ordination method, was used to explore the vital soil factors leading to changes in the microbial community structure with Canoco for Windows 4.5. The Monte Carlo permutation test (number of permutations, 999) was conducted with P values to assess the significance of soil variables in accounting for the shifts in microbial community parameters.


## 3. Results

### 3.1. Soil Physicochemical Analyses

Percentage of SOC decreased vertically but increased along four stages (Table2). The maximum value of SOC in S1 of stage IV was 1.8–3 times that of the other stages, and stage I had the lowest value. Variation in TN followed the same pattern as SOC. Although C/N ratio declined with depth, it increased along succession. Soil pH ranged from 8.38 to 9.31 and increased with soil depth. Soil salinity values decreased markedly with deeper soil layers, especially in stage I, and also decreased along with succession stages. Moisture content reduced not only in vertical direction but also along succession. Value of bulk density was the lowest in S1 and decreased with each successional stage but fluctuated in S2 and S3.Table 2
Soil parameters in successional stages.
Stage
Layer
SOC (%)
TN (%)
C/N ratio
pH
Moisture content (g g−1)
Bulk density (g cm−3)
Salinity (%)
I
S1
1.33 (0.06)Aa
0.06 (0.00)Aa
22.10 (0.66)Aa
8.38 (0.05)Ae
18.78 (0.43)Aa
1.31 (0.06)Aa
5.22 (0.18)Aa
S2
1.29 (0.02)ABa
0.03 (0.01)Ba
29.36 (2.52)Ba
8.59 (0.03)Be
20.69 (0.35)Ba
1.47 (0.09)Aa
1.85 (0.45)Ba
S3
1.22 (0.04)Ba
0.03 (0.00)Ba
40.01 (3.05)Ca
8.71 (0.02)Cf
21.70 (0.35)Ca
1.47 (0.06)Aa
0.94 (0.69)Ba
II
S1
1.70 (0.13)Aab
0.08 (0.02)Aa
20.05 (1.62)Aa
8.90 (0.01)Aa
16.47 (0.59)Ab
1.23 (0.08)Aa
0.71 (0.16)Ac
S2
1.29 (0.03)Ba
0.04 (0.00)Ba
27.26 (0.53)Bb
8.97 (0.03)Bb
20.73 (0.36)Ba
1.40 (0.05)Ba
0.84 (0.03)Bb
S3
1.23 (0.01)Ba
0.03 (0.01)Bab
34.54 (0.51)Ca
9.06 (0.02)Cb
21.13 (0.23)Ba
1.30 (0.02)ABabc
0.78 (0.02)Ba
III
S1
2.29 (0.39)Ab
0.14 (0.04)Aa
16.08 (1.25)Ab
8.45 (0.02)Ad
15.50 (0.58)Ab
1.03 (0.12)Ab
0.67 (0.24)Ac
S2
1.33 (0.01)Ba
0.04 (0.00)Ba
35.43 (1.26)Bb
8.74 (0.02)Ba
16.22 (0.09)Ac
1.41 (0.08)Ba
0.48 (0.01)Abc
S3
1.24 (0.01)Ba
0.04 (0.00)Bbc
43.16 (2.58)Cb
8.77 (0.01)Ca
17.43 (0.41)Bb
1.38 (0.10)Bab
0.45 (0.00)Aa
IV
S1
4.12 (0.88)Ac
0.33 (0.08)Ab
12.72 (0.31)Ac
8.46 (0.01)Ad
14.58 (0.14)Aa
0.88 (0.08)Ac
0.37 (0.03)Ac
S2
1.46 (0.02)Bb
0.05 (0.01)Bb
27.61 (3.31)Bb
9.04 (0.06)Bd
15.93 (0.04)Bb
1.18 (0.09)Bc
0.29 (0.01)Bc
S3
1.45 (0.09)Bb
0.04 (0.01)Bc
31.89 (1.88)Cb
9.31 (0.01)Ce
17.08 (0.42)Cb
1.17 (0.15)Bc
0.22 (0.01)Ca
Numbers in parentheses are standard errors. Capital and small letters represent significant differences among soil depths within the same successional stage and among successional stages within the same soil depth, respectively (P<0.05).


### 3.2. PLFA Composition and Content

We detected 43 PLFAs, with the least in stage I (19), increasing in stages II and III (25 and 26, resp.) and reached the maximum in stage IV (39). As shown in Figure2(a) (P<0.05) total PLFA increased significantly from stages I to IV, while vertically values in S1 were significantly higher than S2 and S3 (P<0.05). Pattern for bacterial (Figure 2(b)) and fungal PLFA (Figure 2(c)) was similar to total PLFA. By contrast, ratio of fungal PLFA to bacterial PLFA (F/B ratio) showed the opposite pattern in S2 and S3, with the value increasing as layers are deepening (Figure 2(d)).Total PLFA (a), bacterial PLFA (b), fungal PLFA (c), and F/B ratio (d) in three layers (S1, S2, and S3) from four successional stages: (I) saline-alkali bare land, (II)Suaeda glauca, (III)Suaeda glauca +Phragmites australis, and (IV)Phragmites australis. Capital and small letters represent significant differences among soil depths within the same successional stage and among successional stages within the same soil depth, respectively (P<0.05).
(a)
Total PLFA
(b)
Bacterial PLFA
(c)
Fungal PLFA
(d)
F/B ratio


### 3.3. Correlation Analysis of PLFA Content and Soil Parameters

The results showed that all types of microbial communities were positively correlated with SOC and TN but negatively correlated with C/N ratio, soil moisture content, and bulk density (Table3). Although fungi showed no significant correlations with pH value, PLFA content of all other microbes displayed negative correlations with pH value. Salinity showed negative correlations only with bacteria and fungi communities but showed no significant correlations with other kinds of microbial PLFA content.Table 3
Spearman’s rank correlation analysis between soil microbial PLFA content and soil parameters.
Microbial community
SOC
TN
C/N ratio
Moisture content
Bulk density
Soil pH
Salinity
Total PLFA
0.724
*
*
0.766
*
*
−
0.731
*
*
−
0.799
*
*
−
0.674
*
*
−
0.461
*
*
ns
Bacterial PLFA
0.753
*
*
0.760
*
*
−
0.737
*
*
−
0.844
*
*
−
0.777
*
*
−
0.481
*
*
−
0.353
*
*
Fungal PLFA
0.720
*
*
0.625
*
*
−
0.587
*
*
−
0.821
*
*
−
0.716
*
*
ns
−
0.389
*
*
Gram-positive bacterial PLFA
0.612
*
*
0.733
*
*
−
0.735
*
*
−
0.661
*
*
−
0.593
*
*
−
0.356
*
*
ns
Gram-negative bacterial PLFA
0.809
*
*
0.824
*
*
−
0.803
*
*
−
0.813
*
*
−
0.781
*
*
−
0.509
*
*
ns
“+” denotes positive correlation, and “−” denotes negative correlation.Significant** at P<0.01. Significant* at P<0.05. “ns” means not significant.


### 3.4. Microbial Community Ratio Alterations

Along with vertical soil layers deepening, the proportion of fungal (Figure3(a)) and Gram-negative bacterial PLFA (Figure 3(b)) decreased, whereas bacterial and Gram-positive bacterial PLFA increased. In addition, there was a large variation between S1 and S2 (Figures 3(a) and 3(b)) but a small variation between S2 and S3. In successional gradient, the proportion of fungal (Figure 3(a)) and Gram-negative bacterial PLFA (Figure 3(b)) in all layers increased with each stage, whereas bacterial and Gram-positive bacterial PLFA decreased in all layers.Ratios of fungal to bacterial PLFA (a) and Gram-negative bacterial to Gram-positive bacterial PLFA (b) in three layers (S1, S2, and S3) from four successional stages: (I) saline-alkali bare land, (II)Suaeda glauca, (III)Suaeda glauca +Phragmites australis, and (IV)Phragmites australis.
(a)
(b)


### 3.5. RDA Analyses

RDA revealed the relationship between two ratios (F/B and Gram-negative bacteria to Gram-positive bacteria) and seven (SOC, TN, C/N ratio, pH, moisture content, bulk density, and salinity) key environmental factors. From Figure4, we found that both SOC and TN markedly affected F/B and Gram-negative/Gram-positive ratios (P<0.001), while moisture content, pH, and C/N ratio were negatively correlated with the two ratios, and the effect for pH was less than moisture content and C/N ratio. Besides, salinity and bulk density showed no significant effects on and were not main environmental factors to the two ratios.Figure 4
Redundancy discriminate analysis representing seven parameters (soil organic carbon, total nitrogen, C/N ratio, soil moisture, bulk density, pH, and salinity) and ratios for two microbial groups: fungal PLFA/bacterial PLFA and Gram-negative/Gram-positive bacteria.


## 3.1. Soil Physicochemical Analyses

Percentage of SOC decreased vertically but increased along four stages (Table2). The maximum value of SOC in S1 of stage IV was 1.8–3 times that of the other stages, and stage I had the lowest value. Variation in TN followed the same pattern as SOC. Although C/N ratio declined with depth, it increased along succession. Soil pH ranged from 8.38 to 9.31 and increased with soil depth. Soil salinity values decreased markedly with deeper soil layers, especially in stage I, and also decreased along with succession stages. Moisture content reduced not only in vertical direction but also along succession. Value of bulk density was the lowest in S1 and decreased with each successional stage but fluctuated in S2 and S3.Table 2
Soil parameters in successional stages.
Stage
Layer
SOC (%)
TN (%)
C/N ratio
pH
Moisture content (g g−1)
Bulk density (g cm−3)
Salinity (%)
I
S1
1.33 (0.06)Aa
0.06 (0.00)Aa
22.10 (0.66)Aa
8.38 (0.05)Ae
18.78 (0.43)Aa
1.31 (0.06)Aa
5.22 (0.18)Aa
S2
1.29 (0.02)ABa
0.03 (0.01)Ba
29.36 (2.52)Ba
8.59 (0.03)Be
20.69 (0.35)Ba
1.47 (0.09)Aa
1.85 (0.45)Ba
S3
1.22 (0.04)Ba
0.03 (0.00)Ba
40.01 (3.05)Ca
8.71 (0.02)Cf
21.70 (0.35)Ca
1.47 (0.06)Aa
0.94 (0.69)Ba
II
S1
1.70 (0.13)Aab
0.08 (0.02)Aa
20.05 (1.62)Aa
8.90 (0.01)Aa
16.47 (0.59)Ab
1.23 (0.08)Aa
0.71 (0.16)Ac
S2
1.29 (0.03)Ba
0.04 (0.00)Ba
27.26 (0.53)Bb
8.97 (0.03)Bb
20.73 (0.36)Ba
1.40 (0.05)Ba
0.84 (0.03)Bb
S3
1.23 (0.01)Ba
0.03 (0.01)Bab
34.54 (0.51)Ca
9.06 (0.02)Cb
21.13 (0.23)Ba
1.30 (0.02)ABabc
0.78 (0.02)Ba
III
S1
2.29 (0.39)Ab
0.14 (0.04)Aa
16.08 (1.25)Ab
8.45 (0.02)Ad
15.50 (0.58)Ab
1.03 (0.12)Ab
0.67 (0.24)Ac
S2
1.33 (0.01)Ba
0.04 (0.00)Ba
35.43 (1.26)Bb
8.74 (0.02)Ba
16.22 (0.09)Ac
1.41 (0.08)Ba
0.48 (0.01)Abc
S3
1.24 (0.01)Ba
0.04 (0.00)Bbc
43.16 (2.58)Cb
8.77 (0.01)Ca
17.43 (0.41)Bb
1.38 (0.10)Bab
0.45 (0.00)Aa
IV
S1
4.12 (0.88)Ac
0.33 (0.08)Ab
12.72 (0.31)Ac
8.46 (0.01)Ad
14.58 (0.14)Aa
0.88 (0.08)Ac
0.37 (0.03)Ac
S2
1.46 (0.02)Bb
0.05 (0.01)Bb
27.61 (3.31)Bb
9.04 (0.06)Bd
15.93 (0.04)Bb
1.18 (0.09)Bc
0.29 (0.01)Bc
S3
1.45 (0.09)Bb
0.04 (0.01)Bc
31.89 (1.88)Cb
9.31 (0.01)Ce
17.08 (0.42)Cb
1.17 (0.15)Bc
0.22 (0.01)Ca
Numbers in parentheses are standard errors. Capital and small letters represent significant differences among soil depths within the same successional stage and among successional stages within the same soil depth, respectively (P<0.05).


## 3.2. PLFA Composition and Content

We detected 43 PLFAs, with the least in stage I (19), increasing in stages II and III (25 and 26, resp.) and reached the maximum in stage IV (39). As shown in Figure2(a) (P<0.05) total PLFA increased significantly from stages I to IV, while vertically values in S1 were significantly higher than S2 and S3 (P<0.05). Pattern for bacterial (Figure 2(b)) and fungal PLFA (Figure 2(c)) was similar to total PLFA. By contrast, ratio of fungal PLFA to bacterial PLFA (F/B ratio) showed the opposite pattern in S2 and S3, with the value increasing as layers are deepening (Figure 2(d)).Total PLFA (a), bacterial PLFA (b), fungal PLFA (c), and F/B ratio (d) in three layers (S1, S2, and S3) from four successional stages: (I) saline-alkali bare land, (II)Suaeda glauca, (III)Suaeda glauca +Phragmites australis, and (IV)Phragmites australis. Capital and small letters represent significant differences among soil depths within the same successional stage and among successional stages within the same soil depth, respectively (P<0.05).
(a)
Total PLFA
(b)
Bacterial PLFA
(c)
Fungal PLFA
(d)
F/B ratio


## 3.3. Correlation Analysis of PLFA Content and Soil Parameters

The results showed that all types of microbial communities were positively correlated with SOC and TN but negatively correlated with C/N ratio, soil moisture content, and bulk density (Table3). Although fungi showed no significant correlations with pH value, PLFA content of all other microbes displayed negative correlations with pH value. Salinity showed negative correlations only with bacteria and fungi communities but showed no significant correlations with other kinds of microbial PLFA content.Table 3
Spearman’s rank correlation analysis between soil microbial PLFA content and soil parameters.
Microbial community
SOC
TN
C/N ratio
Moisture content
Bulk density
Soil pH
Salinity
Total PLFA
0.724
*
*
0.766
*
*
−
0.731
*
*
−
0.799
*
*
−
0.674
*
*
−
0.461
*
*
ns
Bacterial PLFA
0.753
*
*
0.760
*
*
−
0.737
*
*
−
0.844
*
*
−
0.777
*
*
−
0.481
*
*
−
0.353
*
*
Fungal PLFA
0.720
*
*
0.625
*
*
−
0.587
*
*
−
0.821
*
*
−
0.716
*
*
ns
−
0.389
*
*
Gram-positive bacterial PLFA
0.612
*
*
0.733
*
*
−
0.735
*
*
−
0.661
*
*
−
0.593
*
*
−
0.356
*
*
ns
Gram-negative bacterial PLFA
0.809
*
*
0.824
*
*
−
0.803
*
*
−
0.813
*
*
−
0.781
*
*
−
0.509
*
*
ns
“+” denotes positive correlation, and “−” denotes negative correlation.Significant** at P<0.01. Significant* at P<0.05. “ns” means not significant.


## 3.4. Microbial Community Ratio Alterations

Along with vertical soil layers deepening, the proportion of fungal (Figure3(a)) and Gram-negative bacterial PLFA (Figure 3(b)) decreased, whereas bacterial and Gram-positive bacterial PLFA increased. In addition, there was a large variation between S1 and S2 (Figures 3(a) and 3(b)) but a small variation between S2 and S3. In successional gradient, the proportion of fungal (Figure 3(a)) and Gram-negative bacterial PLFA (Figure 3(b)) in all layers increased with each stage, whereas bacterial and Gram-positive bacterial PLFA decreased in all layers.Ratios of fungal to bacterial PLFA (a) and Gram-negative bacterial to Gram-positive bacterial PLFA (b) in three layers (S1, S2, and S3) from four successional stages: (I) saline-alkali bare land, (II)Suaeda glauca, (III)Suaeda glauca +Phragmites australis, and (IV)Phragmites australis.
(a)
(b)


## 3.5. RDA Analyses

RDA revealed the relationship between two ratios (F/B and Gram-negative bacteria to Gram-positive bacteria) and seven (SOC, TN, C/N ratio, pH, moisture content, bulk density, and salinity) key environmental factors. From Figure4, we found that both SOC and TN markedly affected F/B and Gram-negative/Gram-positive ratios (P<0.001), while moisture content, pH, and C/N ratio were negatively correlated with the two ratios, and the effect for pH was less than moisture content and C/N ratio. Besides, salinity and bulk density showed no significant effects on and were not main environmental factors to the two ratios.Figure 4
Redundancy discriminate analysis representing seven parameters (soil organic carbon, total nitrogen, C/N ratio, soil moisture, bulk density, pH, and salinity) and ratios for two microbial groups: fungal PLFA/bacterial PLFA and Gram-negative/Gram-positive bacteria.


## 4. Discussion

Microbial biomass reflects the size of populations involved in regulating energy and nutrient cycle of soil [20, 21]. Results above suggested that the average height, coverage, and biomass of dominant species nearly increased along with successional stages and also for altitude (Table 1). In addition, soil parameters (Table 2) affected microbial biomass and structures (Figure 2). Spearman’s correlation analysis (Table 3) showed that all microbial biomass was positively correlated with SOC and TN (P<0.01), which suggested that SOC and TN levels can be represented by total PLFA, bacteria, or fungi content. Some researchers [22–25] have also reached the same conclusion in other ecosystems. All microbe biomass was negatively correlated with C/N ratio, soil moisture, and bulk density, indicating that appropriate permeability and lower soil moisture content were beneficial to microbes [26]. Cook and Papendick [27] reported that bacteria were active at high water potential levels but fungi at low levels. However, we found the same results to bacteria (Figure 2(b)) and the opposite result to fungi (Figure 2(c)), which may be due to microtopography (Table 1). It is generally recognized that fungi are more tolerant to drought; then the reason why the highest value of fungi PLFA content appeared in successional stage IV is that the plot was located at the highest altitude. The highest topography gave rise to the lowest water potential and the most nutritious soil rich in SOC and TN. The pH was related to all microbial biomass, except for fungi in our study. However, Bååth and Anderson [28] reported that microbial biomass, especially fungi, was negatively related with pH. Aciego Pietri and Brookes [29] and Wu et al. [30] reported that soil with higher pH was richer in Gram-negative than Gram-positive bacteria. In our study, total microbial, bacterial, and fungal biomass increased along with succession accompanying salinity decline. It indicated that higher salinity inhibited growth of microbes, which may be responsible for sparse vegetation in degraded wetlands [31]. Thus, soil structure, air and water permeability, and nutrient influenced microbial communities during vegetation succession [32] and then these factors affected vegetation development and succession. No consensus has been reached regarding the impact of vegetation on microbes in coastal ecosystems. Su et al. [33] reported that presence or absence ofTypha angustifolia did not influence biomass of bacteria and fungi in sediments of wetlands.Microbial biomass showed obviously vertical distribution and decreased with depth increase (Figure2). Microbes lived mainly in topsoil, which demonstrated that rhizosphere was active interface for mass exchange between soil and plants. Total microbial biomass (Figure 2(a)) increased along succession, suggesting that succession promoted quantity of microbial accumulation. Similar results were observed for bacteria (Figure 2(b)) and fungi (Figure 2(c)). F/B ratio (Figure 2(d)) peaked in S1 of stage IV suggesting that the higher the F/B ratio, the greater the stability of wetland ecosystem, and it meant thatPhragmites australis had reached the climatic successional stage. Therefore, F/B ratio reflected degree of succession. de Vries et al. [34] reached the same conclusion. Bacteria were dominant in extreme saline-alkali wetland because it showed higher biomass than fungi in each successional stage, which is consistent with the report by Li et al. [35]. In general, the greater the amount of fungi, the better the relative soil condition. Fungal biomass rose significantly in stages III and IV as a result of roots development byPhragmites australis, which reduced soil bulk density, improved permeability of topsoil, and benefited the survival of fungi.The proportions of fungi and Gram-negative bacteria increased when high carbon material was introduced into the soil, while the proportion of Gram-positive bacteria was reduced (Figure3). It could hold true in study by Griffiths et al. [36]. Gram-negative bacteria were dominant microbial group. The proportion of Gram-negative bacteria increased along with vegetation succession. By contrast, decrease of Gram-positive bacteria may be due to a unique mechanism of adapting acidic environments [37], or pH may have increased the effectiveness of sugars and amino acids to promote Gram-negative bacterial growth. RDA also indicated that (Figure 3) the major factors especially SOC and TN impacted on ratios of microbes (P<0.001).


## 5. Conclusions

Taken together, our study indicated that the evolution of soil microbial communities underground in extreme saline-alkali wetland was closely associated with the succession of vegetation aboveground. Microbe biomass and structure changed significantly not only along with vegetation succession, but also in vertical soil layers. The colonization of pioneerSuaeda glauca improved the soil condition, including the increase of SOM and TN as well as the decrease of soil density and salinity. It was not due to a single factor but the result of multiple cooperating factors. SOC and especially TN levels were the main factors affecting microbial communities. However, soil salinity had no effect on fungi/bacteria and Gram-negative/Gram-positive ratios. Microorganism responded to microenvironment changes through rhizosphere. Bacteria were dominant, and fungi played an important role in soil nutrient cycling with succession.



---

*Source: 491347-2014-11-09.xml*