Journal of Oceanology and Limnology   2021, Vol. 39 issue(3): 880-891     PDF       
http://dx.doi.org/10.1007/s00343-020-0073-y
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

XIA Pinhua, ZHANG Jian, LIU Jinbo, YU Lifei
Shifts of sediment bacterial community and respiration along a successional gradient in a typical karst plateau lake wetland (China)
Journal of Oceanology and Limnology, 39(3): 880-891
http://dx.doi.org/10.1007/s00343-020-0073-y

Article History

Received Feb. 12, 2020
accepted in principle Apr. 28, 2020
accepted for publication Jul. 2, 2020
Shifts of sediment bacterial community and respiration along a successional gradient in a typical karst plateau lake wetland (China)
Pinhua XIA1,2, Jian ZHANG2,3, Jinbo LIU4, Lifei YU2,3     
1 Guizhou Key Laboratory for Mountainous Environmental Information and Ecological Protection, Guizhou Normal University, Guiyang 550001, China;
2 College of Life Science, Guizhou University, Guiyang 550025, China;
3 Collaborative Innovation Center for Mountain Ecology & Agro-Bioengineering, Guizhou University, Guiyang 550025, China;
4 Department of Hepatobiliary Surgery, the Affiliated Hospital of Southwest Medical University, Luzhou 646000, China
Abstract: Bacteria are important regulators of carbon cycling in lakes and are central to sediment ecosystem processes. However, the sediment microbial communities and their respiratory responses to the lake wetland succession are poorly understood. In this study, we collected sediment samples from four different succession points (the Potamogeton lucens zone, the Scirpus tabernaemontani zone, the Scirpus triqueter zone, and the Juncus effusus zone) in the Caohai Wetland of the Guizhou Plateau (China). The bacterial communities at these succession points were studied using a high-throughput sequencing approach. The sediment microbial respiration (SR) was measured using static chambers in the field and basal respiration (BR) was determined in the laboratory. The results show that the dominant bacterial taxa in the sediment was Proteobacteria (34.7%), Chloroflexi (17.8%), Bacteroidetes (7.3%), Acidobacteria (6.6%), and Cyanobacteria (6.1%). Principal coordinate analysis showed that the microbial community structure differs significantly at different sampling points along the successional gradient, indicating that the bacterial community structure is sensitive to the lake wetland succession. Different hydrological regimes and soil characteristics such as NH4+-N, Fe2+, Mn2+, and sediment organic carbon (SOC) content may be important factors responsible for the differences in the sediment microbial characteristics of the different successional stages in the Caohai wetland. Additionally, it was found that the SR increased significantly from the P. lucens zone to the J. effusus zone, but BR had the opposite response. The shifts in the bacterial community structure along the successional gradient may be the main reason for the observed differences in sediment respiration.
Keywords: sediment bacteria    microbial respiration    lake wetland    Caohai Wetland    
1 INTRODUCTION

Wetlands, the transitional zone between aquatic and terrestrial environments, have important ecological functions, and play an important role in carbon cycling. In wetlands, the long-term flooding and hypoxia lead to an imbalance between net primary production and decomposition. Therefore, wetlands represent a major sink for atmospheric carbon (C) (Schulze and Freibauer, 2005). Although wetlands occupy less than 4%–6% of the terrestrial land surface worldwide, they contain approximately one-third of the world's carbon (15×1014 kg) stored in terrestrial soil reservoirs (Kraigher et al., 2006; Kayranli et al., 2010).

Climate change is expected to increase the frequency and intensity of droughts in the future (Cook et al., 2014; Chen et al., 2018). Droughts in wetlands caused by climate change and human activity have been widely observed in recent years (Liu et al., 2006; Wu et al., 2017). It is estimated that more than half of the world's wetlands will disappear in the next century. Both drought and drainage can result in significant reductions in water levels, leading to wetland succession. Wetland succession occurs because of lake filling. The succession includes various stages, specifically, the submerged plant stage, the erect plant stage, and the wet herbaceous stage. These different succession stages have different hydrological conditions, soil characteristics, and vegetation communities (Xia et al., 2015).

Wetland succession often threatens carbon stores, as increasing decomposition leads these wetlands to serve as a carbon source rather than a carbon sink. Drought-induced carbon loss in wetlands at high latitudes has been observed (Fenner and Freeman, 2011; Planas-Clarke et al., 2020). Falling water levels result in increased oxygen availability in the soil, which potentially causes more rapid decomposition of organic matter. It is well known that soil respiration in wetlands under anoxic conditions is much lower than that in oxic environments (Chen et al., 2018). However, some studies have found that drought in wetlands may lead to unchanged or reduced CO2 emission in some wetlands, thus increasing the carbon levels in these environments (Laiho, 2006; Noormets et al., 2010). Thus, the long-term outcome of sediment carbon remains unclear. Hence, it remains a significant pressing global problem to understand wetland succession effects and their inherent feedback mechanisms related to carbon decomposition (Wang et al., 2015). Because of changes in the climate or surrounding land use, wetland succession can result in altered local environmental conditions that can in turn directly affect the carbon cycle (Laiho, 2006).

Bacteria are important regulators of carbon cycling in lakes and are central to the processes of sediment ecosystem processes. Several studies have shed light on the link between wetland soil respiration and microbial community composition (Jaatinen et al., 2008; Mäkiranta et al., 2009), but the effects of wetland success to the sediment bacterial communities and respiration are still unknown. Caohai is a complete and a typical karst plateau lake wetland system, located in the Guizhou Plateau within the Changjiang (Yangtze) River watershed area. As the water level declined, the Caohai Lake showed succession from a shallow lake (submerged plant stage) to swamp (erect plant stage) to swampy meadow (wet herbaceous stage). In this study, we investigated the bacterial communities in sediment samples along the succession gradient using the high-throughput sequencing of 16S rRNA genes. Sediment microbial respiration (SR) was measured using the static chambers in the field, and sediment basal respiration (BR) was determined in the laboratory. While both SR and BR are measures of soil microbial respiration. SR reflects the in-situ soil respiration, while BR reflects the optimal microbial respiration. In other words, BR can be recognized as a measure of the aerobic decomposition potential of a substrate, that is, how much substrate can be utilized under the optimal oxygen and temperature conditions (Jaatinen et al., 2008). We aimed to: 1) characterize the changes in the bacterial community structure along the succession gradient; 2) understand how sediment respiration responds to the wetland succession and the relation with the microbial community structures; and 3) identify the main factors influencing microbial communities and sediment respiration.

2 MATERIAL AND METHOD 2.1 Site description

The Caohai Wetland (26°47′32″N–26°52′52″N, 104°10′16″E–104°20′40″E) is located on the Guizhou Plateau, in southwestern China, at an altitude of 2 170 m above sea level and covers an area of 120 km2, with water covering 25 km2 of the area (Xia et al., 2020). Caohai is subject to a subtropical climate, with an average annual temperature of 10.5 ℃ and a relative humidity of 79%. The average annual precipitation is 950.9 mm, with the majority (70%) of rainfall occurring from May to September (Sun et al., 2020). The unique geographical and climatic conditions contribute to a rich fauna and flora and the Caohai Wetland was designated as the National Wetland Nature Reserve in 1992.

In the 1970s, due to artificial drainage for grain farming, the water area of the Caohai Lake has decreased sharply from 45 km2 to 5 km2 but recovered to more than 20 km2 after the 1980s. Because of the impact of regional climate change and human activities, the eco-hydrological situation of the Caohai Lake changed significantly over the years, with serious degradation of the wetland resulting in a clear hydrological gradient where the vegetation shifts from the submerged to the mesophyte succession. Four sample points (CH1, CH2, CH3, and CH4) were chosen based on the vegetation distribution. CH1 was located in the open water area of the lake, where Potamogeton lucens is the dominant vegetation species at the deepest water level. In CH2, Scirpus tabernaemontani is the dominant vegetation species at the sub-deepest water level. In CH3, the vegetation is dominated by Scirpus triqueter and shows seasonal submersion. In CH4, the driest area that is only occasionally submerged, J. effusus dominates (Xia et al., 2015).

2.2 Sediment sampling and experimental design

Sediment samples were collected in May 2014. In the P. lucens (CH1) and S. tabernaemontani (CH2) zones, undisturbed sediment core samples were collected using a plexiglas corer; while in the S. triqueter (CH3) and J. effusus (CH4) zones, samples comprising the top 10 cm of sediment were collected using a soil auger; Triplicate sampling plots were set in each site (Fig. 1), and five sediment cores with a depth of 0 to 10 cm (excluding surface litter) were taken randomly from each plot, and mixed into one composite sample in a sterile plastic bag, transported to the laboratory on ice and kept at -80 ℃. In addition, the surface material in five plots (0.1 m2), equally spaced along each sites, was clipped, and collected for aboveground plant biomass measurement.

Fig.1 The sketch map of study areas at the Caohai Lake CH1: Potamogeton lucens zone; CH2: Scirpus tabernaemontani zone; CH3: Scirpus triqueter zone; CH4: Juncus effusus zone. Triplicate sampling plots were set in each site.

In the laboratory, samples from each plot were homogenized for each plot and sieved through 2-mm mesh to remove visible stones, grass, and roots. A portion of the fresh sediment was used immediately for determination of sediment respiration and total water moisture content, and another portion was airdried and stored at room temperature until measurement of the other sediment characteristics. The remainder of the fresh sediment was stored at -80 ℃ for subsequent DNA extraction.

2.3 Sediment physicochemical properties analysis

Sediment physicochemical properties were measured using standard methods. Briefly, the sediment pH was measured by a pH meter (Sartorius PB-10, Germany), using a sediment to water ratio of 1꞉2.5. The Oxidation-Reduction Potential (ORP) was determined by an ORP analyzer (HAD-FJA-6, China). Ammonium (NH4+-N) and sediment nitrate (NO3--N) were extracted with 1-mol/L KCl and their concentrations measured using an Automated Discrete Analyzer (AQ2+, SEAL Analytical Inc., United Kingdom) (Guo et al., 2015). The sediment water content (SWC) was gravimetrically determined after drying at 105 ℃ for 12 h. The sediment organic carbon (SOC) was measured with the dichromate digestion method, and the dissolved organic carbon (DOC) using a TOC analyzer (Vario TOC, Elementar, Germany). The sediment bulk density (BD) was determined by the measurements of the cutting ring. The concentrations of Fe2+ and Mn2+ were analyzed using Atomic Absorption Spectroscopy (Zeenit700p, Analytik Jena, Germany). The water level (WL) was measured with a ruler. The plant biomass (PB) was the dry weight of the aboveground part of the plant.

2.4 Sediment respiration measurements

The SR in the CH1 and CH2 zones was measured using static floating opaque chambers (50 cm×50 cm ×50 cm). The Field SR in the CH3 and CH4 zones was measured using the same chambers as used in the CH1 zone; in addition, aluminum collars were installed in the study sites prior to study initiation. At four time points, at 30-min intervals, a 40-mL gas sample was collected from each of the chambers using a polypropylene syringe.

The sediment BR was measured in the laboratory. Twenty grams of original fresh sediment samples were incubated in 250-mL serum bottles at 25 ℃ for 24 h. The concentration of CO2 was measured via Gas Chromatography (GC-14C, Shimadzu) using a 10-m×2-mm stainless steel column and an Flame Ionization Detector (FID). The inlet, column, and detector temperatures were 50 ℃, 50 ℃, and 90 ℃, respectively. H2 served as a carrier gas with 40 mL/min flow rate (Tang et al., 2011a).

2.5 DNA extraction, library preparation, and sequencing

A Fast DNA spin kit for soil (MP Biomedical, CA, USA) was used to isolate gDNA from the sediment samples using the manufacturer's protocols. The 515F and 907R primers (Xiong et al., 2012) were then used for amplifying the bacterial 16S rRNA V4-V5 region, and each sample was given an 8-nucleotide barcode during amplification. The thermocycler settings were: 95 ℃ for 2 min, 25 cycles of 95 ℃ for 30 s, 55 ℃ for 30 s and 72 ℃ for 5 min. Three PCR reactions were conducted per sample, with the total volume of 20-μL volume containing 10-ng template DNA, 2-μL dNTPs (2.5 mmol/L), 4 μL of 5×FastPfu Buffer (Transgen, Biotech, China), 0.4 μL of FastPfu Polymerase (Transgen), and 0.8-μL each of forward and reverse primers (5 μmol/L).

After amplification, 2% agarose gel electrophoresis was performed and an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, CA, USA) was used to isolate the resultant amplicons. Next, QuantiFluorTM-ST (Promega, WI, USA) was used to quantify sample nucleic acid content. Equimolar amounts of all samples were then pooled and subjected to standard paired-end sequencing (2×250 bp) with an Illumina MiSeq platform (Illumina, San Diego, CA, USA).

2.6 Bioinformatic analysis

QIIME (v1.17) was used for the de-multiplexing and quality filtering of raw FASTQ files based on the following criteria: (ⅰ) truncation of reads was performed at any site with an average quality score of < 20 using a sliding 50-bp window, with resultant reads < 50-bp long being discarded; (ⅱ) any reads that contained only barcodes, ambiguous reads, or two or more nucleotide mismatches were discarded; and (ⅲ) only sequences with >10 bp of overlap were assembled based upon overlapping sequences. In addition, any reads which were not successfully assembled were discarded.

A 97% similarity score served as a cut-off criterion for operational taxonomic unit (OTU) clustering, which was performed using UPARSE v7.1, with UCHIME used to remove chimeric sequences (Edgar et al., 2011; Rui et al., 2015). Ribosomal Database Project (RDP) Classifiers were used to assign 16S rRNA sequence taxonomy against the SILVA (SSU115) database with a 70% confidence threshold (Amato et al., 2013).

2.7 Statistical analysis

To correct for sampling effects, we normalized all samples to an even number of sequences per sample (5604) for further analysis (including phylogenetic distance, Chao1 richness and Shannon's diversity index, and observed species). Principal coordinates analysis (PCoA) was used tor assessing structural changes in bacterial communities based upon BrayCurtis distance, using R (http://cran.r-project.org/web/packages/vegan/index.html) for all analyses. Variance partitioning analysis (VPA) and canonical correspondence analysis (CCA) were conducted with CANOCO5 (Microcomputer Power, Ithaca, NY) in order to assess how different environmental variables (WL, SWC, SOC, NH4+-N, NO3--N, Fe2+, Mn2+, ORP, and PB) contributed to observed phenotypes. Oneway analyses of variance (ANOVAs) were used for comparing relative OTU abundance between samples, with Pearson's correlations were used for comparisons of diversity, species abundance, soil properties, and environmental factors with SPSS 18.0 used for analysis.

3 RESULT 3.1 Changes in sediment properties among the different successional gradient

The basic physicochemical characteristics for the different successional gradients are shown in Table 1. The hydrological regime (SWC and WL) showed a significant declining trend along the wetland succession (from CH1 to CH4), while the ORP was significantly increased (P < 0.05), and the sediment BD increases firstly and then decreases. Most of the physical and chemical properties (SWC, NH4+-N, SOC, DOC, and Fe2+) of the sediment showed decreased in value from the CH1 zone to the CH4 zone. In addition, the Sediment Organic Carbon (SOC), which was highest (30.93%±7.62%) in CH1 decreased to its lowest level (3.34%±0.72%) in CH3. However, the PB and pH values were similar among the different points in the successional gradient, with the exception of the CH3 zone.

Table 1 The sampling stations and sediment physicochemical parameters with standard deviations for each plot
3.2 Changes in the overall bacterial communities along the successional gradient

In total, 98 292 reads were obtained on Illumina sequencing. The average read length was 397 bp, and the samples contained 5 604 to 10 165 reads. A total of 845 OTUs (each OTU containing more than five reads) were identified in these samples based on a 16S rRNA sequence similarity threshold of 97%. The number of OTUs, Chao1 richness index and Shannon index suggested that both richness and diversity were the highest for the CH2 zone samples. Simpson index was the highest in CH3 zone and lowest in the CH2 zone (Table 2). Overall, there was no significant change in the diversity index between the CH1 zone and the CH4 zone, although higher values were seen in the CH2 zone (Table 2).

Table 2 Bacteria diversity indices calculated based on 97% similarity of 16S rRNA gene sequences

The similarity of the bacterial communities in the 12 samples were evaluated via PCoA and the PERMANOVA test. The bacterial communities in these 12 samples clustered well into four groups (Fig. 2): (ⅰ) group Ⅰ contained the three samples from the CH1 zone with the highest water level; (ⅱ) group Ⅱ contained the three samples from the CH2 zone, in which water level was lower than that in the CH1 zone; (ⅲ) group Ⅲ contained the three samples from the CH3 zone, which is often subject to dry to wet cycling; (ⅳ) group Ⅳ contained the three samples from the CH4 which is dry and inundated occasionally. The PERMANOVA test showed that the bacterial community structures differ significantly in the four different successional stages (R2=0.860, P=0.001). These results indicate that the bacterial community structure is sensitive to wetland succession.

Fig.2 Principal coordinate analysis of soil bacterial communities using weighted UniFrac Triplicate samples from the same sites are shown as "S# -1" to "S# -3".
3.3 Changes in the dominant bacterial groups among the different stages of the successional gradient

The relative abundances of the dominant bacterial phyla changed significantly along the successional gradient (Fig. 3). The core bacterial phyla were defined as those detected at least two sites and with relative abundances higher than 5% in at least one site. The average abundance of core bacterial phyla was, in descending order, Proteobacteria (34.7%), Chloroflexi (17.8%), Bacteroidetes (7.3%), Acidobacteria (6.6%), and Cyanobacteria (6.1%). The relative abundances of Acidobacteria, Chloroflexi and Actinobacteria increased, whereas the Cyanobacteria decreased along the successional gradient. However, Proteobacteria and "Unclassified taxa" first increased and then decreased along the successional gradient. Bacteroidetes decreased from the CH1 zone to the CH3 zone, but increased slightly in the CH4 zone. These results indicate that the relative abundance of the core bacterial phyla differed along the different successional gradient.

Fig.3 Changes in the dominant bacterial groups (i.e., the percentage compared to the total number of reads) at the phylum level along the degradation gradient in Caohai Lake
3.4 The relationship between environmental parameters and bacterial communities

Canonical correspondence analysis (CCA) was performed to discern the possible relationships between bacterial community structure and environmental parameters (Fig. 4a). The first axis separated the hydrologically distinct locations from one another. The two axes explained 61.35% of the total variation in sediment bacterial communities, suggesting a remarkable correlation between microbial community structure and environmental factors. Variation partitioning analysis showed that sediment properties (NH4+-N, Fe2+, Mn2+, and SOC) could individually explain 52.9% of the bacterial community variation, while hydrological regime (WL, SWC) and PB individually explained 11.2% and 3.0%, respectively, leaving 32.9% of the variations unexplained (Fig. 4b). Hydrological regime alone explained 11.2% of the variation, but most of its contribution was shared with sediment properties (32.5%). This indicates that the changes in hydrological regime and sediment properties along the successional gradient are the major factors driving the shifts in the bacterial communities, and that the impact of plant biomass on the bacterial community structure is marginal.

Fig.4 Canonical correspondence analysis (a) and corresponding variation partitioning analyses tested the correlation between bacterial community structure and environmental factors (b) In Fig. 4b, the values 11.2%, 3.0%, and 52.9% represent the alone contributions of hydrological conditions, plant biomass and sediment physicochemical factors to community changes, respectively. 32.5% and 1.6% represent the contributions of hydrological conditions, sediment physicochemical and plant biomass interactions to community change, respectively, -0.3% represents the contributions of sediment physicochemical and plant biomass interactions to community change. In addition, -1.6% represents the contributions of hydrological conditions, plant biomass and sediment physicochemical interactions to community change.

The relative abundances of 27 out of 37 different bacterial phyla correlated significantly with the hydrological conditions (WL, SWC) and sediment properties (NH4+-N, Fe2+, Mn2+, and SOC) (Supplementary Table S1), and 11 of them correlated significantly with both water level and soil physicochemical properties; these included Acidobacteria, Actinobacteria, Chloroflexi, and Gemmatimonadetes. However, 17 phyla, including Armatimonadetes, Elusimicrobia, and Fibrobacteres, were significantly correlated with only soil properties and none of the phyla correlated with plant biomass (Supplementary Table S1). This indicated that different phyla have different responses to the changes in wetland environmental conditions.

3.5 Sediment respiration along the successional gradient

SR and BR showed the different patterns along the successional gradient (Fig. 5). The SR varied greatly among locations (P < 0.05): it was the highest (0.85 g CO2/(m2∙h)) in the CH4 zone and the lowest (-0.02 g CO2/(m2∙h)) in the CH1 zone. BR rates showed the opposite pattern. BR significantly decreased along the successional gradient. The highest rates were observed in the CH1 zone and the lowest rates in the CH4 zone. The CH1 zone represents anaerobic conditions, where the progress of decomposition is affected in the sediment. In the CH4 zone, the sediment is exposed to aerobic decomposition, thus the decomposition rate is likely to be higher in this site.

Fig.5 The changes of the sediment respiration (SR) and basal respiration (BR) rates along the succession gradient in Caohai Lake

SR was significantly negatively correlated with water level, SWC, NH4+-N, SOC, DOC, and Fe2+ (P < 0.01), and positively correlated with BD. BR was positively correlated with WL, SWC, NH4+-N, SOC, DOC, and Fe2+ (P < 0.01) and negatively correlated with BD (P < 0.01) (Table 3). The relative abundances of Acidobacteria, Actinobacteria, Armatimonadetes, Gemmatimonadetes, Verrucomicrobia and Cyanobacteria showed a positive correlation with the SR. However, only two phyla (Cyanobacteria and Lentisphaerae) were positively correlated with the BR, with Chloroflexi (P < 0.05) and Actinobactereria (P < 0.05) showing a negative correlation (Supplementary Table S1).

Table 3 Pearson's correlation analysis of hydrological parameters with sediment respiration and basal respiration along the successional gradient
4 DISCUSSION

Although there have been many reports on the study of microbial communities in wetland sediments, the effects of succession on these communities and the wetland carbon stocks are relatively unknown. In this study, the results showed that there are significant differences in the hydrological regime, plant community species composition and sediment physics and chemistry at different wetland successional stages. The microbial communities are also vary from the CH1 zone to the CH4 zone. All these factors correlate significantly with the bacterial community composition. These results suggest that changing environmental factors affect bacterial community composition in the sediments along the wetland succession and that these compositional differences may account for the observed change of respiration, leading to a loss of sediment carbon storage.

4.1 Effect of hydrological regime and soil characteristics on the bacterial community structure

Bacterial communities vary significantly from the CH1 zone to the CH4 zone (Fig. 2). These results are consistent with those of a previous study showing variations in the bacterial community composition along the successional stages of the Chongming Dongtan Wetland (Zhang et al., 2013). Tang et al. (2011b) have reported that microbial community structural shifts are linked to changes in the sediment microenvironment. It has also been found that soil microbial communities can be influenced by a wide range of variables, including sediment physicochemical characteristics and climate warming (Balasooriya et al., 2008; Bond-Lamberty and Thomson, 2010; Burns et al., 2015; Jeanbille et al., 2016; Bulseco et al., 2020). The sediment ORP has been considered to be a fundamental indicator of such changes because many sediment biochemical processes are closely associated with the oxygen concentration in sediment, which is, in turn, related to the water level and vegetation. In this study, it was found that different water levels led to different microenvironments and thus to altered microbial community composition along the successional gradient.

We observed significant relationships between the hydrological regime, sediment physicochemical properties, and microbial community composition. Most of the bacterial phyla that were significantly associated with the water level were also significantly correlated with the soil physicochemical properties. The hydrological regime not only directly affects the microbial community (Ma et al., 2018) but also indirectly affects the microbial community structure by altering the soil physicochemical properties. Therefore, sediment characteristics and hydrology are important factors driving the shifts in soil microbial characteristics of the different parts of the successional gradient in the Caohai Wetland. Sediment properties and hydrology have been observed to significantly affect the microbial community composition in floodplains (Mayr et al., 2020). Vegetation also has an important impact on the spatial variation of sediment microbial communities (Jiang et al., 2013). Plant type and sediment physicochemical properties are the two primary factors driving the composition of the soil microbiome due to complex interactions between plants, sediments, and microbes (Garbeva et al., 2004; Mitchell et al., 2010; Yin et al., 2020). However, in this study, no significant relationship between the bacterial relative abundance and plant biomass was found, indicating that the impacts of plant biomass are marginal. This does not, however, mean that the vegetation is not an important factor affecting the bacterial communities in sediments in the Caohai Wetland, but indicates that the influence of plants on the bacterial community is not related to the plant biomass. The role of plants in the soil microbial community construction should be further investigated in the Caohai Wetland. Above all, different hydrological regimes and sediment characteristics such as NH4+-N, Fe2+, Mn2+, and SOC may be important factors resulting in the differences in sediment microbial characteristics of different successional stages in the Caohai Wetland.

4.2 Effects of wetland succession on SR and BR and their relationships with bacterial community

SR and BR showed different patterns in relation to the different stages of the successional gradient (Fig. 5). SR clearly increased from the CH1 zone to the CH4 zone, whereas BR showed the opposite trend. SR was the sediment respiration measured in situ. Along the succession, the water depth gradually decreased and the ORP gradually increased from the CH1 zone to the CH4 zone, indicating that the sediment environment had changed from anaerobic to aerobic. Therefore, the decomposition of organic matter also changed from anaerobic to aerobic decomposition, resulting in a gradual increase in SR. BR can be recognized as a measure of the aerobic decomposition potential of a substrate, that is, how much substrate can be utilized under the optimum oxygen and temperature conditions (Jaatinen et al., 2008). BR was higher in the P. lucens-dominated zone CH1, because these areas represent anaerobic conditions, where the progress of decomposition is impeded in situ. However, in the J. effusus-dominated zone CH4, the surface soil has been exposed to aerobic decomposition, and thus the decomposition potential is likely to be higher. BR had a significant positive correlation with SOC (Table 3), suggesting that SOC might be an important factor that influences BR. Sediment substrate (carbon) availability might be a key factor that influences the BR (Jaatinen et al., 2008; Schiller et al., 2019). The soil carbon content decreased from the CH1 zone to the CH4 zone, resulting in a lower rate of BR.

Variation of SR at different stages along the successional gradient may be related to the shifts in the relative abundance of specific microbial populations. The shifts of microbial community structure have been shown to be linked to SR changes in wetland sediments, and the differences in sediment microbial community structure in different wetland sediments are important reasons for the different SR values (Li et al., 2011). Different microbial community composition may be a major reason for differences in SR along the successional gradient (Jaatinen et al., 2008). Proteobacteria are a group of aerobic or facultative bacteria that often have highly versatile successional capacities (Yang et al., 2005; Sato et al., 2009), which might contribute to the decomposition of organic matter owing to its high abundance in Caohai Wetland. The phylum Actinobacteria is another major taxon that is capable of organic matter decomposition (Tang et al., 2011b), and the relative abundance of Actinobacteria increased dramatically along the successional gradient, which suggests it might be a key taxon driving the changes of SR in the Caohai Wetland. The relative abundance of Anaerolineae and Acidobacteria also increased along the successional gradient, which might also contribute to the change of SR. Substrate is a strong determinant of the microbial community (Chen et al., 2008; Eilers et al., 2010), as the changes that take place in both the carbon quantity and quality of the wetland succession are likely to affect the soil communities. Therefore, the composition and community-level physiological profiles of the soil microbial communities change along the wetland succession, which, in turn, affect sediment microbial respiration (Tang et al., 2011a).

Overall, differences in the hydrology and vegetation of different stages of the successional gradient could lead to differences in the soil physical and chemical characteristics, including SWC, SOC, NH4+-N and ORP, which might consequently lead to different soil microbial community compositions. This explains why the bacteria with high successional ability, such as Acidobacteria, Actinobacteria and Chloroflexi, were found to be more abundant in CH4 zone with the lowest water level, which may be an important factor contributing to the higher SR in this site.

5 CONCLUSION

These results support the following conclusions: (1) the wetland succession is accompanied by changes in bacterial community structure; (2) the presence of distinct microbial communities in wetland sediment samples may be the main reason for the observed differences in rates of decomposition, which, in turn, affect the sediment carbon cycle; (3) hydrological conditions and sediment properties were the major factors driving the shifts in the bacterial communities. Overall, this study indicates that the wetland succession has significant effects on the composition of soil microbial communities and that these effects can alter the carbon cycle of the lake wetland.

6 DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available from the corresponding author.

Electronic supplementary material

Supplementary material (Supplementary Table S1) is available in the online version of this article at https://doi.org/10.1007/s00343-020-0073-y.

References
Amato K R, Yeoman C J, Kent A, Righini N, Carbonero F, Estrada A, Gaskins H R, Stumpf R M, Yildirim S, Torralba M, Gillis M, Wilson B A, Nelson K E, White B A, Leigh S R. 2013. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. The International Society for Microbial Ecology Journal, 7: 1 344-1 353. DOI:10.1038/ismej.2013.16
Balasooriya W K, Denef K, Peters J, Verhoest N E C, Boeckx P. 2008. Vegetation composition and soil microbial community structural changes along a wetland hydrological gradient. Hydrology and Earth System Sciences Discussions, 12(1): 277-291. DOI:10.5194/hess-12-277-2008
Bond-Lamberty B, Thomson A. 2010. Temperature-associated increases in the global soil respiration record. Nature, 464(7288): 579-582. DOI:10.1038/nature08930
Bulseco A N, Vineis J H, Murphy A E, Spivak A C, Giblin A E, Tucker J, Bowen J L. 2020. Metagenomics coupled with biogeochemical rates measurements provide evidence that nitrate addition stimulates respiration in salt marsh sediments. Limnology and Oceanography, 65(S1): S321-S339.
Burns J H, Anacker B L, Strauss S Y, Burke D J. 2015. Soil microbial community variation correlates most strongly with plant species identity, followed by soil chemistry, spatial location and plant genus. AOB Plants, 7(1): plv030. DOI:10.1093/aobpla/plv030
Chen H Y, Zou J Y, Cui J, Nie M, Fang C M. 2018. Wetland drying increases the temperature sensitivity of soil respiration. Soil Biology and Biochemistry, 120: 24-27. DOI:10.1016/j.soilbio.2018.01.035
Chen J, Xie H J, Zhuang X L, Zhuang G Q, Bai Z H, Zhuang H X. 2008. Substrate-induced changes in microbial community-level physiological profiles and their application to discriminate soil microbial communities. Journal of Environmental Sciences, 20(6): 725-731. DOI:10.1016/s1001-0742(08)62119-1
Cook B I, Smerdon J E, Seager R, Coats S. 2014. Global warming and 21st century drying. Climate Dynamics, 43(9): 2 607-2 627. DOI:10.1007/s00382-014-2075-y
Edgar R C, Haas B J, Clemente J C, Quince C, Knight R. 2011. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics, 27: 2 194-2 200. DOI:10.1093/bioinformatics/btr381
Eilers K G, Lauber C L, Knight R, Fierer N. 2010. Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biology and Biochemistry, 42(6): 896-903. DOI:10.1016/j.soilbio.2010.02.003
Fenner N, Freeman C. 2011. Drought-induced carbon loss in peatlands. Nature Geoscience, 4(12): 895-900. DOI:10.1038/ngeo1323
Garbeva P V, Van Veen J A, Van Elsas J D. 2004. Microbial diversity in soil: selection of microbial populations by plant and soil type and implications for disease suppressiveness. Annual Review of Phytopathology, 42: 243-270. DOI:10.1146/annurev.phyto.42.012604.135455
Guo G X, Kong W D, Liu J B, Zhao J X, Du H D, Zhang X Z, Xia P H. 2015. Diversity and distribution of autotrophic microbial community along environmental gradients in grassland soils on the Tibetan Plateau. Applied Microbiology and Biotechnology, 99(20): 8 765-8 776. DOI:10.1007/s00253-015-6723-x
Jaatinen K, Laiho R, Vuorenmaa A, Del Castillo U, Minkkinen K, Pennanen T, Penttilä T, Fritze H. 2008. Responses of aerobic microbial communities and soil respiration to water-level drawdown in a northern boreal fen. Environmental Microbiology, 10(2): 339-353. DOI:10.1111/j.1462-2920.2007.01455.x
Jeanbille M, Buée M, Bach C, Cébron A, Frey-Klett P, Turpault M P, Uroz S. 2016. Soil parameters drive the structure, diversity and metabolic potentials of the bacterial communities across temperate beech forest soil sequences. Microbial Ecology, 71(2): 482-493. DOI:10.1007/s00248-015-0669-5
Jiang X T, Peng X, Deng G H, Sheng H F, Wang Y, Zhou H W, Tam N F Y. 2013. Illumina sequencing of 16S rRNA tag revealed spatial variations of bacterial communities in a mangrove wetland. Microbial Ecology, 66(1): 96-104. DOI:10.1007/s00248-013-0238-8
Kayranli B, Scholz M, Mustafa A, Hedmark Å. 2010. Carbon storage and fluxes within freshwater wetlands: a critical review. Wetlands, 30(1): 111-124. DOI:10.1007/s13157-009-0003-4
Kraigher B, Stres B, Hacin J, Ausec L, Mahne I, Van Elsas J D, Mandic-Mulec I. 2006. Microbial activity and community structure in two drained fen soils in the Ljubljana Marsh. Soil Biology and Biochemistry, 38(9): 2 762-2 771. DOI:10.1016/j.soilbio.2006.04.031
Laiho R. 2006. Decomposition in peatlands: reconciling seemingly contrasting results on the impacts of lowered water levels. Soil Biology and Biochemistry, 38(8): 2 011-2 024. DOI:10.1016/j.soilbio.2006.02.017
Li Y L, Wang L, Zhang W Q, Wang H L, Fu X H, Le Y Q. 2011. The variability of soil microbial community composition of different types of tidal wetland in Chongming Dongtan and its effect on soil microbial respiration. Ecological Engineering, 37(9): 1 276-1 282. DOI:10.1016/j.ecoleng.2011.03.024
Liu C L, Xie G D, Huang H Q. 2006. Shrinking and drying up of Baiyangdian Lake wetland: a natural or human cause?. Chinese Geographical Science, 16(4): 314-319. DOI:10.1007/s11769-006-0314-9
Ma Y T, Li J Q, Wu J, Kong Z Y, Feinstein L M, Ding X, Ge G, Wu L. 2018. Bacterial and fungal community composition and functional activity associated with lake wetland water level gradients. Scientific Reports, 8: 760. DOI:10.1038/s41598-018-19153-z
Mäkiranta P, Laiho R, Fritze H, Hytönen J, Laine J, Minkkinen K. 2009. Indirect regulation of heterotrophic peat soil respiration by water level via microbial community structure and temperature sensitivity. Soil Biology and Biochemistry, 41(4): 695-703. DOI:10.1016/j.soilbio.2009.01.004
Mayr M J, Besemer K, Sieczko A, Demeter K, Peduzzi P. 2020. Bacterial community composition and function along spatiotemporal connectivity gradients in the Danube floodplain (Vienna, Austria). Aquatic Sciences, 82: 28. DOI:10.1007/s00027-020-0700-x
Mitchell R J, Hester A J, Campbell C D, Chapman S J, Cameron C M, Hewison R L, Potts J M. 2010. Is vegetation composition or soil chemistry the best predictor of the soil microbial community?. Plant and Soil, 333(1): 417-430. DOI:10.1007/s11104-010-0357-7
Noormets A, Gavazzi M J, Mcnulty S G, Domec J C, Sun G, King J S, Chen J Q. 2010. Response of carbon fluxes to drought in a coastal plain loblolly pine forest. Global Change Biology, 16(1): 272-287. DOI:10.1111/j.1365-2486.2009.01928.x
Planas-Clarke A M, Chimner R A, Hribljan J A, Lilleskov E A, Fuentealba B. 2020. The effect of water table levels and short-term ditch restoration on mountain peatland carbon cycling in the Cordillera Blanca, Peru. Wetlands Ecology and Management, 28(1): 51-69. DOI:10.1007/s11273-019-09694-z
Rui J P, Li J B, Wang S P, An J X, Liu W T, Lin Q Y, Yang Y F, He Z L, Li X Z. 2015. Responses of bacterial communities to simulated climate changes in alpine meadow soil of the qinghai-tibet plateau. Applied & Environmental Microbiology, 81: 6 070-6 077. DOI:10.1128/AEM.02161-15
Sato K, Kato Y, Taguchi G, Nogawa M, Yokota A, Shimosaka M. 2009. Chitiniphilus shinanonensis gen. nov., sp. nov., a novel chitin-degrading bacterium belonging to Betaproteobacteria. The Journal of General and Applied Microbiology, 55(2): 147-153. DOI:10.2323/jgam.55.147
Schiller V, Datry D, Corti T, Foulquier R, Tockner A, Marce K. 2019. Sediment respiration pulses in intermittent rivers and ephemeral streams. Global Biogeochemical Cycles, 33(10): 1 251-1 263. DOI:10.1029/2019GB006276
Schulze E D, Freibauer A. 2005. Carbon unlocked from soils. Nature, 437(7056): 205-206. DOI:10.1038/437205a
Sun R G, Yang J, Xia P H, Wu S L, Lin T, Yi Y. 2020. Contamination features and ecological risks of heavy metals in the farmland along shoreline of Caohai plateau wetland, China. Chemosphere, 254: 126828. DOI:10.1016/j.chemosphere.2020.126828
Tang Y S, Wang L, Jia J W, Fu X H, Le Y Q, Chen X Z, Sun Y. 2011a. Response of soil microbial community in Jiuduansha wetland to different successional stages and its implications for soil microbial respiration and carbon turnover. Soil Biology and Biochemistry, 43(3): 638-646. DOI:10.1016/j.soilbio.2010.11.035
Tang Y S, Wang L, Jia J W, Li Y L, Zhang W Q, Wang H L, Sun Y. 2011b. Response of soil microbial respiration of tidal wetlands in the Yangtze River Estuary to different artificial disturbances. Ecological Engineering, 37(11): 1 638-1 646. DOI:10.1016/j.ecoleng.2011.06.004
Wang H J, Richardson C J, Ho M. 2015. Dual controls on carbon loss during drought in peatlands. Nature Climate Change, 5(6): 584-587. DOI:10.1038/nclimate2643
Wu X W, Cao R, Wei X, Xi X Q, Shi P L, Eisenhauer N, Sun S C. 2017. Soil drainage facilitates earthworm invasion and subsequent carbon loss from peatland soil. Journal of Applied Ecology, 54(5): 1 291-1 300. DOI:10.1111/1365-2664.12894
Xia P H, Kou Y Z, Yu L F. 2015. Carbon metabolic soil microbial community in Caohai Karst Plateau degraded wetland: a case study in southwest China. Acta Scientiae Circumstantiae, 35(8): 2 549-2 555. (in Chinese with English abstract) DOI:10.13671/j.hjkxxb.2014.1010
Xia P H, Yan D B, Sun R G, Song X, Lin T, Yi Y. 2020. Community composition and correlations between bacteria and algae within epiphytic biofilms on submerged macrophytes in a plateau lake, southwest China. Science of the Total Environment, 727: 138398. DOI:10.1016/j.scitotenv.2020.138398
Xiong J B, Liu Y Q, Lin X G, Zhang H Y, Zeng J, Hou J Z, Yang Y P, Yao T D, Knight R, Chu H Y. 2012. Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau. Environmental Microbiology, 14(9): 2 457-2 466. DOI:10.1111/j.1462-2920.2012.02799.x
Yang H C, Im W T, An D S, Park W S, Kim I S, Lee S T. 2005. Silvimonas terrae gen. nov., sp. nov., a novel chitin-degrading facultative anaerobe belonging to the 'Betaproteobacteria'. International Journal of Systematic and Evolutionary Microbiology, 55(6): 2 329-2 332. DOI:10.1099/ijs.0.63837-0
Yin X J, Lu J, Wang Y C, Liu G L, Hua Y M, Wan X Q, Zhao J W, Zhu D W. 2020. The abundance of nirS-type denitrifiers and anammox bacteria in rhizospheres was affected by the organic acids secreted from roots of submerged macrophytes. Chemosphere, 240: 124903. DOI:10.1016/j.chemosphere.2019.124903
Zhang Y N, Li Y L, Wang L, Tang Y S, Chen J H, Hu Y, Fu X H, Le Y Q. 2013. Soil microbiological variability under different successional stages of the Chongming Dongtan wetland and its effect on soil organic carbon storage. Ecological Engineering, 52: 308-315. DOI:10.1016/j.ecoleng.2012.10.002