- 1Institute of Resource Biology and Biotechnology, Department of Biotechnology, College of Life Science and Technology, Huazhong University of Science and Technology, Wuhan, China
- 2College of Environmental Science and Engineering, Guilin University of Technology, Guilin, China
- 3Key Laboratory of Karst Dynamics, MNR & GZAR, Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin, China
- 4Key Laboratory of Molecular Biophysics, Ministry of Education, Wuhan, China
Karst wetlands are important in the global carbon and nitrogen cycles as well as in security of water resources. Huixian wetland (Guilin) is the largest natural karst wetland in China. In recent years, groundwater nitrogen pollution has increasingly affected the wetland ecosystem integrity due to anthropogenic activities. In this study, it was hypothesized that autotrophic microbial diversity is impacted with the advent of pollution, adversely affecting autotrophs in the carbon and nitrogen cycles. Autotrophic microbes have important roles in abating groundwater nitrogen pollution. Thus, it is of great significance to study the characteristics of autotrophic bacterial communities and their responses to environmental parameters in nitrogen-polluted karst groundwaters. The abundances of the Calvin–Benson cycle functional genes cbbL and cbbM as well as the autotrophic CO2-fixing bacterial communities were characterized in the karst groundwater samples with different levels of nitrogen pollution. The cbbM gene was generally more abundant than the cbbL gene in the groundwater samples. The cbbL gene abundance was significantly positively correlated with dissolved inorganic nitrogen (DIN) concentration (P < 0.01). In the autotrophic CO2-fixing bacterial communities, Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria of the phylum Proteobacteria were predominant. At the genus level, Rubrivivax and Methylibium were the dominant cbbL gene containing genera, while Halothiobacillus and Endothiovibrio were the dominant genera for the cbbM gene. The abundance of autotrophic CO2-fixing bacterial communities increased but their diversity decreased with the inflow of nitrogen into the karst groundwater system. The community structure of autotrophic CO2-fixing bacteria in the groundwaters was also significantly affected by environmental factors such as the carbonic anhydrase (CA) activity, dissolved inorganic carbon (DIC) concentration, temperature, and oxidation-reduction potential (ORP). Nitrogen inflow significantly changed the characteristics of autotrophic CO2-fixing bacterial communities in the karst groundwaters. Some key genera such as Nitrosospira and Thiobacillus were clearly abundant in the karst groundwaters with high nitrogen levels. Their respective roles in nitrification and denitrification impact nitrogen removal in this ecosystem. The findings in this study provide an important reference for biological abatement of nitrogen pollution in the karst groundwater system.
Introduction
The karst area of China accounts for 15.6% of the total karst area in the world and about one third of the national land area (Yuan and Cao, 2008). Karst groundwater in China accounts for more than one fourth of the total groundwater resources and is an important source of potable water and production water in karst areas (Yuan et al., 1993). These groundwaters are liable to be contaminated because of the special structures of the surface and subsurface double layers in karst areas. Sewage and agricultural non-point source pollution caused by fertilization and breeding are common sources of groundwater nitrogen pollution (Guo et al., 2006). Groundwater nitrogen pollution is a serious threat to water security. It is important, therefore, to solve problems of karst groundwater nitrogen pollution in endeavors to protect groundwater resources.
Autotrophic microbes play an important role in the treatment of groundwater nitrogen pollution (Li et al., 2020). The bedrock composition and dissolved inorganic carbon (DIC) produced by the unique water-rock-gas interaction in karst groundwater can provide the energy and carbon sources for autotrophs. Therefore, karst groundwater is a compatible environment for the growth of autotrophs. Autotrophic CO2-fixing microbes can assimilate CO2 as the only carbon source, and light or reduced inorganic compounds such as H2, H2S, S0, S2O32–, NH4+, NO2–, or Fe2+ as the energy source (Yuan et al., 2011). Therefore, the autotrophic CO2-fixing microbes play an important role in the carbon and nitrogen cycles in aquatic ecosystems. Some autotrophs can oxidize ammonia and nitrite or reduce nitrate to N2 gas and can play important roles in nitrogen removal from water bodies (Wang et al., 2020). Since autotrophs do not need an organic carbon source (Jakus et al., 2021), they have potential application prospects for nitrogen removal in freshwater and groundwater habitats including karst groundwater systems.
Currently, there are six major CO2-fixing pathways in autotrophs: the Calvin–Benson cycle, reductive tricarboxylic acid cycle, reductive acetyl-CoA, 3-hydroxypropionate bicycle, 3-hydroxypropionate/4-hydroxybutyrate cycle, and dicarboxylate/4-hydroxybutyrate cycle (Thauer, 2007; Berg, 2011; Hügler and Sievert, 2011; Jeoung et al., 2014). The Calvin-Benson cycle is the most important pathway for bacterial assimilation of CO2 (Berg, 2011). The key enzyme is the ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) that catalyzes the carboxylation of ribulose 1,5-bisphosphate to yield two molecules of 3-phosphoglycerate (Hügler and Sievert, 2011). The RubisCO Type I (encoded by gene cbbL) and Type II (encoded by gene cbbM) are the main and most studied RubisCO enzymes (Berg, 2011; Liu et al., 2017). Due to the high conservation of these two functional genes and the appropriate nucleotide length, cbbL and cbbM genes have been used as specific probes for molecular analysis of the autotrophic microbial communities in the environment (Wu et al., 2014; Wang et al., 2021). These genes have also been used to characterize the diversity of autotrophic bacteria and their responses to environmental conditions (Long et al., 2015; Wrighton et al., 2016; Alfreider et al., 2017). However, there are only few karst groundwater studies on the characterization of autotrophic microbial communities and their responses to environmental changes.
The Huixian wetland in Guilin, Guangxi Province is the most representative karst wetland in China and even in the global tropical, subtropical middle, and low altitude karst areas (Cai et al., 2012). It is the largest natural karst wetland of China. The groundwater in the Huixian karst wetland was affected by sewage pollution and non-point source pollution caused mostly by farming and agriculture (Peng et al., 2018). It was hypothesized that autotrophic microbial diversity is impacted with the advent of pollution, adversely affecting autotrophs in the carbon and nitrogen cycles of karst wetland system. The purpose of this study was to characterize the autotrophic CO2-fixing bacterial communities in the Huixian karst groundwaters with different nitrogen levels and their relationships with dissolved inorganic nitrogen species and other selected groundwater environmental factors. A karst groundwater vein with different nitrogen levels in the aquifer was selected in the Huixian karst wetland for this study. Of special interest were the differences in the community structures of CO2-fixing bacteria in the karst groundwaters with high and low nitrogen levels because these differences are poorly understood in karst groundwaters. The results of this study will provide a scientific and biological basis for removal of nitrogen pollution, and an important reference for the protection and utilization of karst groundwater resources.
Materials and methods
Study area
The study area is in the Huixian karst wetland, Fengjia Village, Huixian County, Guilin, Guangxi Province, China. The Huixian karst wetland belongs to the subtropical monsoon climate with an annual average temperature range of 16.5∼20.5°C. The temperature in the summer is significantly different from that in the winter. The temperature difference between the spring and autumn is small, and the temperature in the autumn is slightly higher than that in the spring. The average annual rainfall is about 1890.4 mm. The rainy season is mainly concentrated from April to July, and October through March next year is the dry season, with January being the driest and coldest period (Cai et al., 2012).
Sampling
A groundwater vein (110°13′31″ ∼ 110°15′10″E, 25°09′13″ ∼ 25°05′51″N) from Baitengdang in the north to Fengjia in the south was selected in the east of Huixian karst wetland, as shown in Figure 1. Villages and farmlands are widely distributed along this groundwater vein. The main villages it flows through are Baitengdang, Qin Village, Shang Village, Nan Village, Dalu Village, Mamian Village and Fengjia Village, of which Fengjia Village has a large area of farming and agriculture. Thirteen groundwater sampling sites (D1–D13) were selected from the north to south transect along this groundwater vein. They represented three types of karst groundwater, namely a karst spring (Spring), a shallow pumping well (Well), and a shallow underground river (Shallow groundwater). Sites D1, D3, D4, and D13 were karst springs, sites D2, D5, D6, D10, and D11 were shallow pumping wells, and sites D7, D8, D9, and D12 were shallow underground rivers. The lithology of the sampling sites was mainly limestone, and mixed limestone and dolomite. The land use patterns around the sampling sites were mainly shrubs and farmland. The hydrogeological information and land use patterns of various sampling sites are presented in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g001.jpg)
Figure 1. Schematic diagram of the study area and sampling sites in the Huixian karst groundwaters. (A) Location of the Huixian karst wetland; (B) Location of the sampling sites in the Huixian karst groundwaters.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-t001.jpg)
Table 1. Hydrogeological information and land use patterns at various sampling sites in the Huixian karst groundwaters, China.
Sampling was carried out in January (winter, 201801), April (spring, 201804), July (summer, 201807), and October (autumn, 201810). Before sampling, the sampler (a polyethylene plastic bucket) was first washed three times with purified water followed by three rinses with groundwater from the sampling site. At each site, 5–10 L groundwater samples were collected from three points, then combined and mixed, and stored in polyethylene plastic containers. The samples were filtered with 0.22 μm cellulose acetate membrane filters to entrap the microorganisms from the groundwater. The membrane filters were then placed in sterile petri dishes and then placed in a dry-ice box, brought back to the laboratory and stored at –80°C for total DNA extraction.
About 500 mL of filtered groundwater samples were acidified with 2 mL of 1 mol/L sulfuric acid and stored in 500 mL polyethylene bottles in an ice box, brought back to the laboratory and stored at 4°C for analysis of NH4+, NO2–, NO3–, SO42–, and dissolved organic carbon (DOC). For the analysis of total nitrogen (TN) and total organic carbon (TOC), about 500 mL unfiltered groundwater samples at each site were fixed with 2 mL of saturated HgCl2 solution, stored in 500 mL polyethylene bottles in an ice box, brought back to the laboratory, and stored at 4°C before analysis.
Determination of groundwater physicochemical properties
Groundwater physicochemical parameters, including pH, temperature (T), turbidity, conductivity, oxidation-reduction potential (ORP), dissolved oxygen (DO), Ca2+, and DIC (HCO3–) concentration were measured in situ. The pH was measured using a Mettler-Toledo Seven2GO S2 portable pH meter with an accuracy of ± 0.01 pH unit. The T and turbidity were measured using a Eureka Manta2 multi-parameter water quality monitor (Eureka Water Probes, Austin, TX, United States). The accuracy of the measurements was ±0.1°C and ±0.02 NTU, respectively. The conductivity was measured using a CT-3030 conductivity probe (Kedida, Shenzhen) with an accuracy of ± 0.01 mS/cm. The ORP was measured using a CT-8022 portable ORP pen (Kedida, Shenzhen) with an accuracy of ±0.1 mV. The DO was measured using a JPB-607A portable dissolved oxygen meter (Ramag, Shanghai) with an accuracy of ± 0.01 mg /L. The concentration of Ca2+ was measured using a Aquamerck Calcium Kit (Merck) with an accuracy of ±0.1 mg/L. The concentration of DIC was measured using an Aquamerck Alkalinity Kit (Merck) with an accuracy of ±0.05 mmol/L (Nzung’a et al., 2018).
The concentrations of TN, NH4+, NO2–, and NO3– of groundwater samples were determined according to the Methods for Monitoring and Analyzing Water and Waste water (2002). The concentration of dissolved inorganic nitrogen (DIN) is the sum of NH4+, NO2–, and NO3– concentrations. The TOC and DOC concentrations were determined using a total organic carbon analyzer (Multi N/C 2100, Analytik, Jena, Germany). The SO42– concentration was determined using an Ion Chromatography with an accuracy of ± 0.01 mg/L.
Measurement of groundwater carbonic anhydrase activity
The carbonic anhydrase (CA) activity of groundwater samples was determined from the rate of CO2 hydration by following the change of pH traced on a chart recorder as previously described (Li et al., 2005; Nzung’a et al., 2018). The basic physicochemical properties and CA activity at various sites of the Huixian karst groundwaters are presented in Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-t002.jpg)
Table 2. Basic physicochemical properties and CA activity at the sampling sites of the Huixian karst groundwater.
Extraction of total DNA of groundwater microorganisms
PowerWater DNA Isolation Kit (Qiagen, Hilden, Germany) were used to extract the total DNA of the biomass on the membrane filters according to the instruction of the manufacturer. The DNA was dissolved in 100 μL of sterile deionized water, and the quality was checked by 1% agarose gel electrophoresis. DNA concentration and purity were determined with a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States). The DNA samples were stored at –80°C for subsequent analysis.
Detection of cbbL and cbbM genes
The abundances of the 16S rRNA gene and the cbbL and cbbM genes of bacteria in the Huixian karst groundwater samples were analyzed by quantitative PCR (qPCR). The amplification primers are shown in Supplementary Table 1. The reaction mixture was 10 μL including 5 μL SYBR PremixEX TaqTM, 0.5 μL 10 mM forward primer, 0.5 μL 10 mM reverse primer, 1 μL template DNA, and 3 μL ddH2O. The reaction conditions were as follows: pre-denaturation at 95°C for 60 s; 40 cycles of denaturation at 95°C for 10 s, annealing at specific annealing temperature for 10 s, and extension at 72°C for 30 s; and melting at 95°C for 15 s, 60°C for 60 s, and 95°C for 15 s. The annealing temperatures for the 16S rRNA, cbbL, and cbbM genes were 50, 54, and 57°C, respectively.
High-throughput sequencing
The corresponding primers listed in Supplementary Table 1 were used to amplify the cbbL and cbbM genes by PCR. The PCR reaction system was 30 μL including 3 μL template DNA, 25 μL 2 × Premix Taq, 1 μL 10 mM forward primer, and 1 μL 10 mM reverse primer. The reaction conditions were as follows: pre-denaturation at 94°C for 5 min; 30 cycles of denaturation at 94°C for 30 s, annealing at 54°C (cbbL) or 57°C (cbbM) for 30 s, extension at 72°C for 30 s, and a final extension at 72°C for 10 min. The PCR amplification products were detected by 2% agarose gel electrophoresis, and the target fragments were cut and recovered using a gel recovery kit (Omega Bio-tek, Norcross, GA, United States). The recovered target gene fragments were sequenced on an Illumina MiSeq 300 platform by Personal Biotechnology Co., Ltd. (Shanghai, China).
Sequence analysis
The raw data obtained by paired-end sequencing were first spliced by FLASH and then filtered by the Quantitative Insights Into Microbial Ecology (QIIME, v1.8.0) pipeline (Caporaso et al., 2010). After removing low-quality bases and adaptor contamination sequences, the data filtering was completed. The chimeras were removed by UCHIME software to obtain effective sequences. The effective sequences were clustered into operational taxonomic units (OTUs) at 97% sequence similarity using the UCLUST method in QIIME software (Edgar, 2010). The Alpha diversity indexes (Chao 1, ACE, Shannon, and Simpson indexes) of the samples were calculated by QIIME software, and the OTU taxonomic classification was conducted by BLAST searching the representative sequences set against the NCBI NT Database using the best hit.
The raw sequences of CO2-fixing bacteria in the Huixian karst groundwaters were deposited in the NCBI Sequence Read Archive under the accession number PRJNA821389.
Statistical analysis
Differences in the abundances of the cbbL and cbbM genes among different seasons were analyzed by one-way analysis of variance (ANOVA) at P < 0.05 using SPSS 18.0. Relationships between the groundwater environmental parameters and functional gene abundances as well as diversity indexes were analyzed by the Pearson correlation using SPSS 18.0. Relationships between the groundwater environmental parameters and the relative abundance of bacterial genera were analyzed by the Spearman correlation heatmap using SPSS 18.0. The redundancy analysis (RDA) was performed to explore the relationship between the CO2-fixing bacterial community structure and groundwater environmental parameters using the R software package. The linear discriminant analysis (LDA) effect size (LEfSe) method was used to identify significant differences in the community structure of CO2-fixing bacteria in the karst groundwaters between the high (HN) and low (LN) nitrogen levels using the galaxy platform1. Histograms of linear discriminate analysis (LDA) and cladograms of phylogenetic distribution were prepared to present the differentially abundant carbon-fixing bacterial taxa in the karst groundwaters. Variance partitioning analysis (VPA) was performed to distinguish the relative contribution of water environmental factors associated with carbon (environmental carbon factors for short) and water environmental factors associated with nitrogen (environmental nitrogen factors) on the autotrophic CO2-fixing bacterial communities in the karst groundwaters using the “varpart” function in the vegan package. The environmental carbon factors included DIC, DOC, TOC, and CA, and the environmental nitrogen factors included DIN, TN, NH4+, NO2–, and NO3–.
Results
Variation of carbon- fixing genes abundance in the karst groundwaters
In the Huixian karst groundwaters, the abundance of the carbon-fixing gene cbbM was generally higher than that of cbbL. The abundance of the cbbM gene in all samples ranged from 9.76 × 105 to 2.56 × 1011 copies/L, while the cbbL gene ranged from 4.93 × 105 to 1.79 × 109 copies/L in four seasons (Figure 2). The abundances of the cbbL and cbbM genes were significantly different among the four seasons (P < 0.05), and they were also significantly higher in the winter and spring than in the summer and autumn (P < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g002.jpg)
Figure 2. Abundances of the cbbL (left, blue) and cbbM (right, orange) genes in the Huixian karst groundwaters. (A) Winter (201801); (B) Spring (201804); (C) Summer (201807); (D) Autumn (201810).
The proportion of potential CO2-fixing bacteria in the bacterial community in the karst groundwaters was estimated based on the abundances of the 16S rRNA, cbbL and cbbM genes. In all samples, the ratio of cbbL/16S rRNA gene ranged from 9.68 × 10–11 to 2.31 × 10–5, while the ratio of cbbM/16S rRNA ranged from 1.08 × 10–10 to 1.75 × 10–3. The proportion of the cbbL gene was generally lower than that of cbbM gene except at a few sampling sites (Figure 3). The results of the one-way ANOVA showed that there was no significant difference in the proportion of carbon-fixing genes among the seasons.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g003.jpg)
Figure 3. Ratios of gene copy numbers of cbbL (left, blue) or cbbM (right, orange) to 16S rRNA genes in the Huixian karst groundwaters. (A) Winter (201801); (B) Spring (201804); (C) Summer (201807); (D) Autumn (201810).
Effects of nitrogen and other environmental factors on the abundance of CO2-fixing genes
To investigate whether nitrogen levels influenced the abundance of CO2-fixing bacterial community in the karst groundwaters, Pearson correlation analysis was conducted between the abundance of CO2-fixing genes and DIN content. The results showed that the cbbL gene abundance was significantly positively correlated with the DIN content (P < 0.01). The higher the nitrogen level, the higher the cbbL gene abundance (Figure 4). Although the cbbM gene abundance was positively correlated with DIN content, the correlation was not significant (P = 0.077).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g004.jpg)
Figure 4. Relationship between the carbon-fixing genes and DIN concentration in the Huixian karst groundwaters. (A) cbbL gene; (B) cbbM gene.
The results of Pearson correlation analysis between the abundance of carbon-fixing genes and various environmental factors showed that the cbbL gene abundance was significantly positively correlated with ORP, DOC, and turbidity, while negatively correlated with groundwater T (P < 0.05). The cbbM gene abundance was significantly positively correlated with ORP, TOC and DOC, while negatively correlated with the CA activity (P < 0.05) (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-t003.jpg)
Table 3. Results of the correlation analysis between carbon-fixing gene abundances and various environmental factors in the Huixian karst groundwaters.
Community structure and diversity of the CO2-fixing bacteria in the karst groundwaters
The results showed that the physicochemical properties and the CO2-fixing gene abundances in the Huixian karst groundwaters were significantly different between summer (201807) and winter (201801) (Table 1 and Figure 2). Therefore, the summer and winter were selected as the representatives of rainy season and dry season, respectively. High-throughput sequencing analysis was carried out on the groundwater samples in these two seasons to investigate the community characteristics of CO2-fixing bacteria. The rarefaction curves of the cbbL and cbbM genes reached the saturation plateau (Supplementary Figure 1), suggesting that the sequencing depth basically covered all species in the samples. The observed species richness in the CO2-fixing bacterial community containing the cbbL and cbbM genes in all samples was within 1200. The richness (ACE, Chao1) and diversity indexes (Simpson and Shannon) of the cbbL gene-containing bacterial community were higher than those of the cbbM gene except for a few sampling sites (D1, D2, D10, and D12) in the summer (Supplementary Tables 2, 3). This indicated that the diversity of CO2-fixing bacterial community using RubisCO Type I was generally higher than that using RubisCO Type II for carbon assimilation in the Huixian karst groundwaters.
A total of 26 samples were obtained for sequencing from the Huixian karst groundwaters representing the summer and winter seasons. 645429 sequences and 85985 OTUs were obtained from the cbbL gene sequencing, belonging to 7 phyla, 11 classes, 17 orders, 30 families and 45 genera. 1111850 sequences and 114438 OTUs were obtained from the cbbM gene sequencing, belonging to 1 phylum, 4 classes, 8 orders, 15 families, and 23 genera.
The results of the cbbL gene sequencing showed that at the phylum level, Proteobacteria was dominant in the groundwater samples from most sites, and its average relative abundance was 95.97% with a maximum of 100%. In summer, the Cyanobacteria phylum was dominant with an average relative abundance of 79.48% in the water samples from sites D7, D8, and D9. At the class level, except for a few sampling sites (D7, D8, and D9) in summer, Betaproteobacteria was in the dominant position with an average relative abundance of 83% in most sampling sites, followed by Alphaproteobacteria with an average relative abundance of 3.7%. Many Cyanobacteria species could not be identified at the class level in the samples from sites D7, D8, and D9 in the summer (Figure 5A). At the genus level, the top 20 genera with the relative abundance were selected to construct the distribution graph of the relative abundance of species (Figure 5B). The top genera were Rubrivivax, Cyanobium, Methyllium, Hydrogenophaga, and Cupriavidus with a total relative abundance of 92.1%. Rubrivivax was the most abundant species in most groundwater samples, and was found in the groundwater samples from all sites in the winter and summer seasons, with the highest relative abundance of 83.1% and the average relative abundance of 49.1%. Cyanobium was the second and was found only in the groundwater samples from sites D1, D6 and D7 in the winter and from sites D1, D6, D7, D8, and D9 in the summer. Its relative abundance in the samples from site D9 in the summer was the highest (98.1%), and the overall average relative abundance was 12.2%. Methyllium was found at all sampling sites in the winter and summer seasons, with an average relative abundance of 11.8%. It was dominant at sites D12 and D13 in the winter, with the relative abundance of 59.0 and 50.2% respectively. Except for the samples from sites D7, D8, D9, and D12 in the summer, Hydrogenophaga was found in the samples from other sites, with an average relative abundance of 10.0%. Cupriavidus was found at all sampling sites in the winter and summer, with an average relative abundance of 9.0%. The relative abundance ranged from 0.1 to 30.3%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g005.jpg)
Figure 5. Relative abundance and composition of the cbbL gene-containing CO2-fixing bacterial community at the class level (A) and genus level (B) in the Huixian karst groundwaters.
The results of the cbbM gene sequencing showed that Proteobacteria was also the dominant bacterial phylum at all sampling sites in the summer and winter seasons (Figure 6). At the class level, Gammaproteobacteria, Betaproteobacteria, and Alphaproteobacteria were dominant in all groundwater samples with average relative abundances of 62.3, 26.6, and 10.8%, respectively (Figure 6A). At the genus level, the top 20 genera with relative abundance were selected to construct the distribution graph of the relative abundance of species (Figure 6B). Halothiobacillus, Endothiovibrio, Sulfuritalea, Magnetospira, Thiobacillus, Limnohabitans, and Rhodopseudomonas were among the abundant genera, with average relative abundances of 33.4, 24.4, 14.2, 5.2, 4.8, 4.2, and 3.0%, respectively. Their total relative abundance was up to 89.2%. Halothiobacillus had the highest average relative abundance and was found at all sampling sites. It was dominant at sites D4, D5, D6, and D7 in the winter and at sites D11 and D12 in the summer, and the highest relative abundance was 77.8%. Endothiovibrio was found at all sampling sites in the winter and summer. It was dominant at site D9 in the winter and at sites D3, D4, D5, and D13 in the summer, and its relative abundance at site D4 in the summer was the highest (80.6%). Sulfuratalea, Magnetospira, Thiobacillus, and Limnohabitans were found at all sampling sites in the winter and summer. The relative abundance of Limnohabitans in the groundwater samples from site D1 in winter was much higher than in other samples, and the relative abundance was 35.0%. The relative abundance of Rhodopseudomonas in the groundwater samples at site D3 in the winter was much higher than in other samples, and the relative abundance was 33.2%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g006.jpg)
Figure 6. Relative abundance and composition of the cbbM gene-containing CO2-fixing bacterial community at the class level (A) and genus level (B) in the Huixian karst groundwaters.
Thus, these results indicated that there were substantial temporal and spatial differences in the community structure of CO2-fixing bacteria containing the cbbL and cbbM genes in the Huixian karst groundwaters.
Distribution of CO2-fixing bacterial communities as affected by dissolved inorganic nitrogen levels
According to the DIN concentration, the groundwater samples in the winter and summer were divided into high (DIN > 3.0 mg/L) and low nitrogen levels. The results of LEfSe analysis showed that there were significant differences in the distribution of some CO2-fixing bacterial taxa in the groundwater samples between the high and low nitrogen levels. For the cbbL gene-containing CO2-fixing bacterial communities, some autotrophs such as Nitrosomonadales, Nitrosomonadaceae, Thiobacillaceae, Nitrosospira and Thiobacillus were relatively abundant in the groundwaters with high nitrogen level (Figure 7A). For the cbbM gene-containing CO2-fixing bacterial communities, Rhodobacterales, Rhodobacteraceae and Rhodobacter were significantly more abundant in the groundwaters with high nitrogen level, while Chromatiaceae was significantly more abundant with low nitrogen level (Figure 7B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g007.jpg)
Figure 7. Differentially abundant CO2-fixing bacterial taxa as assessed using a histogram of linear discriminate analysis (LDA) and cladogram of phylogenetic distribution with effect size measurements (LEfSe) in the Huixian karst groundwaters. (A) cbbL gene; (B) cbbM gene.
Effects of environmental factors on the structure and diversity of CO2-fixing bacterial community
Pearson correlation analysis was conducted between the richness indexes ACE and Chao1 as well as the diversity indexes Simpson and Shannon of the CO2-fixing bacterial community and environmental factors. The richness and diversity indexes of the cbbL gene-containing CO2-fixing bacterial community were significantly positively correlated with ORP in water (P < 0.05), while significantly negatively correlated with the concentration of inorganic nitrogen compounds in water (P < 0.05) (Table 4). The richness and diversity indexes of the cbbM gene-containing CO2-fixing bacterial community were significantly negatively correlated with NH4+, NO2–, and DO concentrations in water (P < 0.05) (Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-t004.jpg)
Table 4. Results of the correlation analysis between the community diversity of the cbbL gene-containing carbon-fixing bacteria and environmental factors in the Huixian karst groundwaters.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-t005.jpg)
Table 5. Results of the correlation analysis between the community diversity of the cbbM gene-containing carbon-fixing bacteria and environmental factors in the Huixian karst groundwaters.
The results of the RDA analysis showed that the main groundwater environmental factors affecting the community structure of the cbbL gene-containing CO2-fixing bacteria in the winter were CA (r2 = 0.5270, P = 0.025) and DIC (r2 = 0.5040, P = 0.032), while in the summer they were T (r2 = 0.4964, P = 0.030) and ORP (r2 = 0.4666, P = 0.036) (Figure 8). The main factors affecting the community structure of the cbbM gene-containing CO2-fixing bacteria in the summer were CA (r2 = 0.5401, P = 0.015) and DIC (r2 = 0.4891, P = 0.035) (Figure 9), whereas in the winter, there was no significant correlation between the environmental factors and the community structure. The distribution of various sites in the RDA diagram was scattered, indicating that the community structures of CO2-fixing bacteria at various sites varied due to the differences in environmental factors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g008.jpg)
Figure 8. Redundancy analysis (RDA) of the cbbL gene-containing CO2-fixing bacterial community structures with environmental factors in the Huixian karst groundwaters. (A) Winter (201801); (B) Summer (201807).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g009.jpg)
Figure 9. Redundancy analysis (RDA) of the cbbM gene-containing CO2-fixing bacterial community structures with environmental factors in the Huixian karst groundwaters. (A) Winter (201801); (B) Summer (201807).
The results of VPA showed that the environmental carbon and nitrogen factors explained 24 and 15% of the variation, respectively, in the cbbL gene-containing CO2-fixing bacterial community (Figure 10A). For the cbbM gene-containing CO2-fixing bacterial community, the environmental carbon and nitrogen factors explained 10 and 7% of the variation, respectively (Figure 10B). Thus, carbon and inorganic nitrogen compounds in the groundwater affected the variation in the CO2-fixing bacterial community.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g010.jpg)
Figure 10. Variance partitioning analysis (VPA) of CO2-fixing bacterial communities containing the cbbL (A) or cbbM (B) genes based on environmental carbon and nitrogen factors.
Spearman correlation heatmap analysis was conducted between the relative abundance of varied-fixing bacterial genera and the groundwater environmental factors (Figure 11). The relative abundance of more bacterial genera in the cbbL-gene containing bacterial community was significantly correlated with the concentration of inorganic nitrogen components in the groundwaters. Specifically, the relative abundances of Thiobacillus, Halorhodospira, Limnohabitans, Thiomonas, Rhodobacter, and Nitrosopira were significantly positively correlated with NH4+ (P < 0.05). Rhodopseudomonas was negatively correlated with NO2– (P < 0.05). Thiobacillus, Nitrosospira, and Sulfuricaulis were significantly positively correlated with NO3– (P < 0.05). Thiobacillus and Nitrosospira were significantly positively correlated with the concentration of dissolved inorganic nitrogen (DIN) (P < 0.05). For the cbbM-gene containing CO2-fixing bacterial community, there were fewer bacterial genera which were significantly correlated with the concentration of inorganic nitrogen compounds in the groundwaters. Phaeospirillum and Magnetovibrio were negatively correlated with NH4+ (P < 0.05). Thioalkalicoccus was negatively correlated with the concentration of inorganic nitrogen compounds (P < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/949208/fmicb-13-949208-HTML-r1/image_m/fmicb-13-949208-g011.jpg)
Figure 11. Heatmaps of the Spearman correlation between the distribution of the CO2-fixing bacterial communities containing the cbbL (A) or cbbM (B) genes and environmental factors in the Huixian karst groundwaters.
Other groundwater environmental factors were also significantly correlated with the relative abundance of many bacterial genera in the CO2-fixing bacterial community (P < 0.05). For example, DOC was significantly negatively correlated with Sinorhizobium (P < 0.05), and significantly positively correlated with Limnohabitans, Halorhodospira, Thiobacillus, and Nitrosospira (P < 0.05) in the cbbL-gene containing carbon-fixing bacterial community. For the cbbM-gene containing CO2-fixing bacterial community, DOC was significantly negatively correlated with Magnetovibrio and Magnetospirillum (P < 0.05), and significantly positively correlated with Acidithiobacillus (P < 0.05).
Discussion
Community structure and diversity of CO2-fixing bacteria in the Huixian karst groundwaters
The cbbL and cbbM genes were detected in all sampling sites of the Huixian karst groundwaters. The cbbL and cbbM CO2-fixing gene abundances were higher in the Huixian karst groundwaters than in freshwater lakes (Alfreider et al., 2017), seawaters (Fang et al., 2015), and pristine limestone aquifers (Herrmann et al., 2015). These results indicated that a large number of autotrophic CO2-fixing bacteria used the Calvin-Benson cycle for carbon fixation in the Huixian karst groundwaters. Autotrophic CO2-fixing bacterial communities play an important role in the CO2 fixation in aquatic environments (Zhao et al., 2020). The abundant CO2-fixing bacterial communities implied a higher carbon fixation potential in the Huixian karst groundwaters. In most samples from the Huixian karst groundwaters, the ratio of cbbM/16S rRNA was significantly higher than that of cbbL/16S rRNA, suggesting that the CO2-fixing bacterial communities using the RubisCO Type II for carbon assimilation accounted for a higher proportion in the Huixian karst groundwaters. Low O2 and high CO2 concentration are conducive to the growth of cbbM gene-containing CO2-fixing bacterial communities (Badger and Bek, 2008). Due to the special water-rock-gas interaction in the karst groundwaters, the concentration of CO2 dissolved in the water is relatively high, which may favor the growth of the cbbM gene-containing CO2-fixing bacterial community in this karst groundwater environment.
The diversity index of the CO2-fixing bacterial community containing the cbbL gene was higher than that containing the cbbM gene in most samples from the Huixian karst groundwaters. Thus, the diversity of CO2-fixing bacterial communities using RubisCO Type I for carbon assimilation was higher in the karst groundwaters. Microbial community diversity responds to changes in environmental physicochemical properties (Delgado-Baquerizo et al., 2016; Cavicchioli et al., 2019). A high level of microbial diversity is associated with ecosystem stability. In this study, DIN was significantly negatively correlated with the diversity index of the cbbL gene-containing bacterial community. The diversity index of cbbM gene-containing bacterial community was significantly negatively correlated with NH4+ and NO2–. These data indicated that nitrogen influx in the Huixian karst groundwaters decreased the diversity, compromising the stability of the CO2-fixing bacterial community.
The Proteobacteria and Cyanophyta were the main predominant phyla in the CO2-fixing bacterial community, and Proteobacteria occupied a high proportion in most sampling sites from the Huixian karst groundwaters. The Betaproteobacteria and Gammaproteobacteria of the phylum Proteobacteria were the dominant bacterial classes. These data were basically consistent with the reports (Ji et al., 2016) that the CO2-fixing bacterial communities in the surface water of northern South China Sea were dominated by Gamma- and Betaproteobacteria, Cyanobacteria, and Firmicutes.
Rubrivivax of the Betaproteobacteria was the dominant genus in the cbbL gene-containing CO2-fixing bacterial community in the Huixian karst groundwaters. Members in this genus can grow photoautotrophically under light conditions and chemoautotrophically under dark conditions, and they can be found in soils, lakes, and groundwaters (Kane et al., 2001; Gomila et al., 2010; Wu et al., 2015). Gao et al. (2019) reported that Rubrivivax was also associated with a denitrifying community. Zhao (2016) established that Rubrivivax gelatinosus was the main denitrifier in low ammonia-nitrogen water bodies and may play an important role in the removal of nitrogen pollution in water bodies. Members of Hydrogenophaga was also found in large numbers in the cbbL gene-containing CO2-fixing bacterial community in the present study. Hydrogenophaga belongs to hydrogen-oxidizing bacteria. As the chemoautotrophic bacteria with the fastest growth, hydrogen-oxidizing bacteria can obtain energy by oxidizing H2 and assimilate inorganic carbon source into organic carbon (Hu et al., 2009; Meng et al., 2010). Tan (2016) found that Halothiobacillus, Sulfuritalea, Magnetospira, Thiobacillus, Limnohabitans, and Rhodopseudomonas were the main autotrophic sulfur-oxidizers in Pearl River waters and possessed both carbon-fixing functional genes and sulfur-oxidizing functional genes. Thiobacillus and Halothiobacillus have been identified in the cbbM gene-containing CO2-fixing bacterial community in ammonia nitrogen-polluted groundwater (Alfreider et al., 2012). In the present study, a large number of both genera were also found in the cbbM gene-containing CO2-fixing bacterial community. Halothiobacillus accounted for a large proportion, and the highest relative abundance was 77.8%. Halothiobacillus is a type of chemoautotrophic sulfur-oxidizing bacteria, which carbon and energy sources required for growth and metabolism are completely derived from the fixation of CO2 and the oxidation of inorganic sulfur compounds (Shively et al., 1998; Kelly and Wood, 2000). Therefore, most cbbM gene-containing CO2-fixing bacterial communities in the Huixian karst groundwaters were related to sulfur oxidation. Herrmann et al. (2015) also found that a large fraction of bacteria in pristine limestone aquifers had the genetic potential for autotrophic CO2 fixation, with energy most likely provided by the oxidation of reduced sulfur and nitrogen compounds.
Effects of environmental factors on CO2-fixing bacterial communities in the Huixian karst groundwaters
The VPA revealed that environmental carbon and nitrogen factors greatly influenced the CO2-fixing bacterial communities in in the Huixian karst groundwaters. Carbon and nitrogen compounds can serve as energy sources (Yang and Wu, 2020) as well as for assimilation into cell biomass. DIC is abundant as HCO3– in the karst groundwaters (Table 2). Tang et al. (2018) found that DIC was positively correlated with the composition of the CO2-fixing bacterial community in lake surface sediments in the northern Qinghai-Tibetan Plateau, which was comparable to that of the present study. Extracellular CA catalyzes the conversion of HCO3– into CO2, efficiently transporting CO2 into cells, and helps maintain high substrate concentration for the RubisCO enzyme active center. The CA activity in the groundwaters reflects the extracellular CA activity of autotrophic microbes to some extent. The correlation between the composition and abundance of CO2-fixing bacteria and the groundwater CA activity is a captivating paradigm. The CA results in the present study were comparable to the positive correlation between species and abundance of microalgae and the water CA activity previously reported for the Huixian karst wetland watershed (Yan et al., 2020).
The increase of inorganic nutrients, especially of nitrogen compounds, can enhance the number of CO2-fixing bacteria in aquatic ecosystems (Teufel et al., 2017). Yue et al. (2019) reported that increases of organic carbon and nitrogen compounds led to increased abundance of CO2-fixing genes in freshwater and hypersaline lakes of the Tibetan Plateau. In the present study, DIN was significantly positively correlated with the cbbL gene abundance, suggesting that the inflow of agricultural nitrogen into the Huixian karst groundwaters stimulated the growth of CO2-fixing bacterial communities in the groundwater system. In addition, Kamble et al. (2014) found that the addition of organic carbon significantly promoted bacterial growth in some high pH and high salinity soils. In the present study, the abundances of the cbbL and cbbM genes were significantly positively correlated with DOC, suggesting that the karst groundwaters with higher DOC content contained more CO2-fixing bacteria.
Additionally, the community structure of CO2-fixing bacteria in the Huixian karst groundwaters had significant correlation with groundwater T and ORP. It is well established that water temperature and nutrients affect the community structure of CO2-fixing bacteria in aquatic environments (Yue et al., 2019). Temperature controls the biological activity and thus affects the rate of microbial metabolism (Gillooly et al., 2001; Brown et al., 2004). Nutrient availability and suitable temperature are among the main drivers of the formation of microbial communities in aquatic ecosystems (Gilbert et al., 2009; Cross et al., 2015). Moreover, Alfreider et al. (2012) found that the differences in redox condition led to differences in the distribution and diversity of autotrophic carbon-fixing bacteria. These results were comparable to those of the present study.
Differences in the distribution of CO2-fixing bacterial community in Huixian karst groundwaters with different nitrogen levels
This study compared the differences in the CO2-fixing bacterial community structure between high and low nitrogen levels through LEfSe analysis to further understand their key role in the nitrogen cycle process in the karst groundwaters. Nitrosospira and Thiobacillus (containing cbbL gene), and Rhodobacter (containing cbbM gene) were enriched in high nitrogen groundwaters. The Spearman correlation analysis showed that the relative abundances of Nitrosospira and Thiobacillus were significantly positively correlated with the concentrations of NH4+, NO3– and DIN (P < 0.05). These results suggested that Nitrosospira and Thiobacillus responded positively to changes in nitrogen in the karst groundwaters, and may actively participate in the carbon and nitrogen cycles in the karst groundwater. Oved et al. (2001) found that Nitrosospira-like populations were dominant in soils irrigated with fertilizer-amended water. Shaw et al. (2006) reported that Nitrosospira spp. could use NO2– to produce N2O via a nitrifier denitrification pathway. They suggested that this nitrifier-denitrification path may be a universal trait in ammonia-oxidizing bacteria of the Betaproteobacteria. Some Thiobacillus species are denitrifiers, such as Thiobacillus denitrificans, which is a sulfur-oxidizing, autotrophic denitrifier, and is one of the main species in the biological denitrification (Ye and Zuo, 2011). Some Rhodobacter species have been reported to be associated with the simultaneous nitrification and denitrification process (Moura et al., 2018). For nitrogen-polluted groundwater with oligoorganic carbon sources, autotrophic denitrification only requires inorganic carbon sources, thus avoiding the secondary pollution of water bodies by the addition of organic carbon sources, which is an effective and safe way to deal with nitrogen pollution (Capua et al., 2019). Autotrophic denitrifying bacteria have been applied in the removal of nitrogen pollution in groundwater systems. For example, Moon et al. (2008) removed nitrate via denitrification in groundwater samples by inoculating a bacterial consortium containing Thiobacillus denitrificans as a key member. The denitrifying community almost completely transformed nitrate (60 mg N L–1) to dinitrogen gas. The discovery of bacterial genera related to autotrophic denitrification in the Huixian karst groundwaters suggests that CO2-fixing bacteria are involved in nitrogen transformations and are important in the nitrogen removal in karst groundwater ecosystem.
Conclusion
In the Huixian karst groundwaters, the abundance of the cbbM gene was higher than that of the cbbL gene. The phylum Proteobacteria was dominant in the autotrophic CO2-fixing bacterial communities containing the cbbL and cbbM genes. Rubrivivax and Methylibium were the dominant genera in the cbbL gene-containing CO2-fixing bacterial communities, while Halothiobacillus and Endothiovibrio were dominant genera in the cbbM gene-containing bacterial communities. In addition to environmental carbon factors (DIC, DOC, TOC, and CA), the abundance and diversity of autotrophic CO2-fixing bacterial communities containing the cbbL and cbbM genes in the karst groundwaters were modulated by inorganic nitrogen compounds. There were significant negative correlations between DIN, NH4+, NO2– and the community diversity of autotrophic CO2-fixing bacteria. Thus, the inflow of nitrogen increased the abundance of autotrophic CO2-fixing bacteria but reduced the community diversity of autotrophs in the karst groundwaters. Some key autotrophic CO2-fixing bacterial genera such as Nitrosospira and Thiobacillus were enriched in the karst groundwaters with high nitrogen levels, and their respective roles in nitrification and denitrification affect nitrogen biotransformations in the karst groundwaters. These findings provide ongoing scientific basis for in-depth understanding of the microbial autotrophs and nitrogen cycling in karst groundwater systems, and provide an important reference for the protection and utilization of karst groundwater resources.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.
Author contributions
WL designed the experiments. XW, TS, and XP collected all samples. XW conducted the experiments and wrote the manuscript. XW and WL analyzed the data. WL, AC, YX, MZ, and LY revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by Guangxi Key Research and Development Project of China (Grant No. Guike AB21196026).
Acknowledgments
We thank the School of Environmental Science and Engineering, Huazhong University of Science and Technology for DOC, TOC, and SO42– analysis. We are grateful to Zaiwang Yang for groundwater sampling and physicochemical property analysis assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.949208/full#supplementary-material
Footnotes
References
Alfreider, A., Baumer, A., Bogensperger, T., Posch, T., Salcher, M. M., and Summerer, M. (2017). CO2 assimilation strategies in stratified lakes: diversity and distribution patterns of chemolithoautotrophs. Environ. Microbiol. 19, 2754–2768. doi: 10.1111/1462-2920.13786
Alfreider, A., Schirmer, M., and Vogt, C. (2012). Diversity and expression of different forms of RubisCO genes in polluted groundwater under different redox conditions. FEMS Microbiol. Ecol. 79, 649–660. doi: 10.1111/j.1574-6941.2011.01246.x
Badger, M. R., and Bek, E. J. (2008). Multiple RubisCO forms in proteobacteria: their functional significance in relation to CO2 acquisition by the CBB cycle. J. Exp. Bot. 59, 1525–1541. doi: 10.1093/jxb/erm297
Berg, I. A. (2011). Ecological aspects of the distribution of different autotrophic CO2 fixation pathways. Appl. Environ. Microbiol. 77, 1925–1936. doi: 10.1128/AEM.02473-10
Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M., and West, G. B. (2004). Toward a metabolic theory of ecology. Ecology 85, 1771–1789. doi: 10.1890/03-9000
Cai, D., Ma, Z., and Jiang, Z. C. (2012). Research on Huixian Karst Wetland Ecosystem. Beijing: Geological Publishing House.
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Capua, F. D., Pirozzi, Lens, P. N. L., and Esposito, G. (2019). Electron donors for autotrophic denitrification. Chem. Eng. J. 362, 922–937. doi: 10.1016/j.cej.2019.01.069
Cavicchioli, R., Ripple, W. J., Timmis, K. N., Azam, F., Bakken, L. R., Baylis, M., et al. (2019). Scientists’ warning to humanity: microorganisms and climate change. Nat. Rev. Microbiol. 17, 569–586. doi: 10.1038/s41579-019-0222-225
Cross, W. F., Hood, J. M., Benstead, J. P., Huryn, A. D., and Nelson, D. (2015). Interactions between temperature and nutrients across levels of ecological organization. Glob. Change Biol. 21, 1025–1040. doi: 10.1111/gcb.12809
Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 7:10541. doi: 10.1038/ncomms10541
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461.
Fang, F., Wang, L., Xi, X. F., Hu, J. J., Fu, X. H., Lu, B., et al. (2015). Potential carbon fixation capability of non-photosynthetic microbial community at different depth of the South China Sea and its response to different electron donors. Environ. Sci. 36, 1550–1556.
Gao, W. X., Yan, J. H., Du, H. Y., and Zhang, K. Q. (2019). Biological response of soil nirS- and nosZ- type denitrifying bacteria to dairy slurry irrigation in community structure and diversity. J. Agro-Environ. Sci. 38, 1089–1100.
Gilbert, J. A., Field, D., Swift, P., Newbold, L., Oliver, A., Smyth, T., et al. (2009). The seasonal structure of microbial communities in the Western English Channel. Environ. Microbiol. 11, 3132–3139. doi: 10.1111/j.1462-2920.2009.02017.x
Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M., and Charnov, E. L. (2001). Effects of size and temperature on metabolic rate. Science 293, 2248–2251. doi: 10.1126/science.1061967
Gomila, M., Pinhassi, J., Falsen, E., Moore, E. R. B., and Lalucat, J. (2010). Kinneretia asaccharophila gen. nov., sp. nov., isolated from a freshwater lake, a member of the Rubrivivax branch of the family Comamonadaceae. Int. J. Syst. Evol. Microbiol. 60, 809–814. doi: 10.1099/ijs.0.011478-0
Guo, F., Jiang, G. H., and Jiang, Z. C. (2006). Groundwater and environmental geological problems in different kinds of karst environment system in karst rock areas of South China. Geol. Sci. Technol. Inf. 25, 83–87.
Herrmann, M., Rusznyák, A., Akob, D. M., Schulze, I., Opitz, S., Totsche, K. U., et al. (2015). Large fractions of CO2-fixing microorganisms in pristine limestone aquifers appear to be involved in the oxidation of reduced sulfur and nitrogen compounds. Appl. Environ. Microbiol. 81, 2384–2394. doi: 10.1128/AEM.03269-14
Hu, J. J., Wang, L., Li, Y. L., Fu, X. H., Le, Y. Q., Xu, D. S., et al. (2009). Breeding, optimization and community structure analysis of non-photosynthetic CO2 assimilation microbial flora. Environ. Sci. 30, 2438–2444.
Hügler, M., and Sievert, S. M. (2011). Beyond the Calvin cycle: autotrophic carbon fixation in the ocean. Annu. Rev. Mar. Sci. 3, 261–289. doi: 10.1146/annurev-marine-120709-142712
Jakus, N., Blackwell, N., Osenbrück, K., Straub, D., Byrne, J. M., Wang, Z., et al. (2021). Nitrate removal by a novel lithoautotrophic nitrate-reducing, iron(II)-oxidizing culture enriched from a pyrite-rich limestone aquifer. Appl. Environ. Microbiol. 87:e00460-21. doi: 10.1128/AEM.00460-21
Jeoung, J.-H., Goetzel, S., Hennig, S. E., Fesseler, J., Wörman, C., Dendra, J., et al. (2014). The extended reductive acetyl-CoA pathway: ATPases in metal cluster maturation and reductive activation. Biol. Chem. 395, 545–558. doi: 10.1515/hsz-2013-0290
Ji, F., Ming, H., Li, H., Zan, S., Wang, J., Su, J., et al. (2016). Diversity of CO2 fixation gene in the surface waters of northern South China Sea in the Calvin cycle. Acta Sci. Circumst. 36, 4037–4043.
Kamble, P. N., Gaikwad, V. B., Kuchekar, S. R., and Bååth, E. (2014). Microbial growth, biomass, community structure and nutrient limitation in high pH and salinity soils from Pravaranagar (India). Eur. J. Soil Biol. 65, 87–95. doi: 10.1016/j.ejsobi.2014.10.005
Kane, S. R., Beller, H. R., Legler, T. C., Koester, C. J., Pinkart, H. C., Halden, R. U., et al. (2001). Aerobic biodegradation of methyl tert-butyl ether by aquifer bacteria from leaking underground storage tank sites. Appl. Environ. Microbiol. 67, 5824–5829. doi: 10.1128/AEM.67.12.5824-5829.2001
Kelly, D. P., and Wood, A. P. (2000). Reclassification of some species of Thiobacillus to the newly designated genera Acidithiobacillus gen. nov., Halothiobacillus gen. nov. and Thermithiobacillus gen. nov. Int. J. Syst. Evol. Microbiol. 50, 511–516. doi: 10.1099/00207713-50-2-511
Li, S., Zhang, Y., Yin, S., Wang, X., Liu, T., and Deng, Z. Q. (2020). Analysis of microbial community structure and degradation of ammonia nitrogen in groundwater in cold regions. Environ. Sci. Pollut. Res. 27, 44137–44147. doi: 10.1007/s11356-020-10318-w
Li, W., Yu, L. J., Yuan, D. X., Wu, Y., and Zeng, X. D. (2005). A study of the activity and ecological significance of carbonic anhydrase from soil and its microbes from different karst ecosystems of Southwest China. Plant Soil 272, 133–141. doi: 10.1007/s11104-004-4335-9
Liu, Y. Y., Wang, S., Li, S. Z., and Deng, Y. (2017). Advances in molecular ecology on microbial functional genes of carbon cycle. Microbiol. China 44, 1676–1689.
Long, X. E., Yao, H., Wang, J., Huang, Y., Singh, B. K., and Zhu, Y. G. (2015). Community structure and soil pH determine chemoautotrophic carbon dioxide fixation in drained paddy soils. Environ. Sci. Technol. 49, 7152–7160. doi: 10.1021/acs.est.5b00506
Meng, Y., Wang, W. W., Chen, X. D., Xiong, B. T., and Tang, M. (2010). Research progress on isolation, screening and plant growth-promoting mechanism of hydrogen-oxidizing bacteria. Microbiol. China 37, 1525–1532.
Methods for Monitoring and Analyzing Water and Wastewater (2002). Methods for Monitoring and Analyzing Water and Wastewater, 4th Edn. Beijing: China Environmental Science Press, 244–284.
Moon, H. S., Shin, D. Y., Nam, K., and Kim, J. Y. (2008). A long-term performance test on autotrophic denltrification column for applicafion a permeable reactive barrier. Chemosphere 73, 723–728. doi: 10.1016/j.chemosphere.2008.06.065
Moura, R. B., Santos, C. E. D., Okada, D. Y., Martins, T. H., Ferraz, J. A. D. N., and Damianovic, M. H. R. Z. (2018). Carbon-nitrogen removal in a structured-bed reactor (SBRRIA) treating sewage: operating conditions and metabolic perspectives. J. Environ. Manage 224, 19–28. doi: 10.1016/j.jenvman.2018.07.014
Nzung’a, S. O., Pan, W., Shen, T., Li, W., Qin, X., Wang, C., et al. (2018). Comparative study of carbonic anhydrase activity in waters among different geological eco-environments of Yangtze River basin and its ecological significance. J. Environ. Sci. 66, 173–181. doi: 10.1016/j.jes.2017.05.027
Oved, T., Shaviv, A., Goldrath, T., and Mandelbaum, R. T. (2001). Influence of effluent irrigation on community composition and function of ammonia-oxidizing bacteria in soil. Appl. Environ. Microbiol. 67, 3426–3433. doi: 10.1128/AEM.67.8.3426-3433.2001
Peng, C., Pan, X. D., Jiao, Y. J., Ren, K., and Zeng, J. (2018). Identification of nitrate pollution sources through various isotopic methods: a case study of the Huixian wetland. Environ. Sci. 39, 5410–5417. doi: 10.13227/j.hjkx.201803212
Shaw, L. J., Nicol, G. W., Smith, Z., Fear, J., Prosser, J. I., and Baggs, E. M. (2006). Nitrosospira spp. can produce nitrous oxide via a nitrifier denitrification pathway. Environ. Microbiol. 8, 214–222. doi: 10.1111/j.1462-2920.2005.00882.x
Shively, J. M., van Keulen, G., and Meijer, W. G. (1998). Something from almost nothing: carbon dioxide fixation in chemoautotrophs. Annu. Rev. Microbiol. 52, 191–230. doi: 10.1146/annurev.micro.52.1.191
Tan, X. Q. (2016). The Diversity and Sulfur Metabolic Pathway of Sulfur Oxidizing Bacteria in the Pearl River Water. Unpublished M.S thesis, China: South China University of Technology.
Tang, Y., Liu, Y., Yang, J., and Jiang, H. (2018). Gene diversity involved in Calvin pathway of carbon fixation and its response to environmental variables in surface sediments of the northern Qinghai-Tibetan Plateau lakes. Earth Sci. 43(Suppl.1), 19–30.
Teufel, A. G., Li, W., Kiss, A. J., and Morgan-Kiss, R. M. (2017). Impact of nitrogen and phosphorus on phytoplankton production and bacterial community structure in two stratified Antarctic lakes: a bioassay approach. Polar Biol. 40, 1007–1022. doi: 10.1007/s00300-016-2025-8
Thauer, R. K. (2007). A fifth pathway of carbon fixation. Science 318, 1732–1733. doi: 10.1126/science.1152209
Wang, H., Li, Y., Zhang, S., Li, D., Liu, X., Wang, W., et al. (2020). Effect of influent feeding pattern on municipal tailwater treatment during a sulfur-based denitrification constructed wetland. Bioresour. Technol. 315:123807. doi: 10.1016/j.biortech.2020.123807
Wang, X., Li, W., Xiao, Y., Cheng, A., Shen, T., Zhu, M., et al. (2021). Abundance and diversity of carbon-fixing bacterial communities in karst wetland soil ecosystems. Catena 204:105418. doi: 10.1016/j.catena.2021.105418
Wrighton, K. C., Castelle, C. J., Varaljay, V. A., Satagopan, S., Brown, C. T., Wilkins, M. J., et al. (2016). RubisCO of a nucleoside pathway known from Archaea is found in diverse uncultivated phyla in bacteria. ISME J. 10, 2702–2714. doi: 10.1038/ismej.2016.53
Wu, X., Ge, T., Wang, W., Yuan, H., Wegner, C. E., Zhu, Z., et al. (2015). Cropping systems modulate the rate and magnitude of soil microbial autotrophic CO2 fixation in soil. Front. Microbiol. 6:379. doi: 10.3389/fmicb.2015.00379
Wu, X. H., Jian, Y., Chen, X. J., Li, B. Z., Yuan, H. Z., Ge, T. D., et al. (2014). Molecular mechanism on carbon dioxide assimilation of autotrophic microorganism and carbon translocation in agricultural soils. Acta Ecol. Sin. 34, 701–709. doi: 10.5846/stxb201211051538
Yan, Z., Li, W., Shen, T., Wang, X., Zhu, M., and Yu, L. (2020). Aquatic microalgal and bacterial communities in a karst wetland watershed. J. Hydrol. 591:125573. doi: 10.1016/j.jhydrol.2020.125573
Yang, Y., and Wu, P. (2020). Soil bacterial community varies but fungal community stabilizes along five vertical climate zones. Catena 195:104841. doi: 10.1016/j.catena.2020.104841
Ye, L., and Zuo, S. (2011). Research advances in autotrophic denitrification by Thiobacillus denitrificans and its influencing factors. Environ. Sci. Manag. 36, 58–63. doi: 10.1016/j.scitotenv.2020.137657
Yuan, D., Zhu, D., Weng, J., Zhu, X., and Han, X. (1993). Chinese karstology. Beijing: Geological Publishing House.
Yuan, H. Z., Qin, H. L., Liu, S. L., Nie, S. A., Wei, W. X., and Wu, J. S. (2011). Advances in research of molecular ecology of carbon fixation microorganism. Sci. Agr. Sin. 44, 2951–2958.
Yue, L., Kong, W., Ji, M., Liu, J., and Morgan-Kiss, R. M. (2019). Community response of microbial primary producers to salinity is primarily driven by nutrients in lakes. Sci. Total Environ. 696:134001. doi: 10.1016/j.scitotenv.2019.134001
Zhao, P. W. (2016). The Ecological Respond of Microorganisms to Environment in Different Nitrogen Polluted Water. Unpublished M.S thesis, China: Huazhong Agricultural University.
Keywords: CO2-fixing bacterial community, cbbL gene, cbbM gene, karst groundwater, nitrogen pollution
Citation: Wang X, Li W, Cheng A, Shen T, Xiao Y, Zhu M, Pan X and Yu L (2022) Community characteristics of autotrophic CO2-fixing bacteria in karst wetland groundwaters with different nitrogen levels. Front. Microbiol. 13:949208. doi: 10.3389/fmicb.2022.949208
Received: 20 May 2022; Accepted: 26 July 2022;
Published: 15 August 2022.
Edited by:
Yong Yuan, Guangdong University of Technology, ChinaReviewed by:
Hui Tong, Guangdong Institute of Eco-environmental and Soil Sciences (CAS), ChinaJibing Li, State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry (CAS), China
Copyright © 2022 Wang, Li, Cheng, Shen, Xiao, Zhu, Pan and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Li, hulwei_009@163.com