Seasonal variations of microbial community structure, assembly processes, and influencing factors in karst river

The physicochemical properties and microbial communities have significant annual and seasonal changes in karst aquifers. To explore the changes of microbial community and their relationships with environmental factors, water samples were collected from a typical karst river. Microbial communities in winter (Jan-2017 and Jan-2019) were stable with high similarity in spite of the 2 years sampling interval, but the microbial communities in Aug-2017 was different from that in Aug-2018. In four sampling times, there were 275 shared genera, whose average relative abundance ranging from 89.04 to 96.27%. The winter and summer specific genera were mainly from the recharge of tributary site K6 and discharge of waste water treatment plant (K2 and K3), respectively. The deterministic processes had a more significant effect on the microbial community assembly in winter than that in summer, which was affected by environmental pressure from pollution. Furthermore, antibiotics and inorganic nitrogen pollution affected element cycles of nitrogen and sulfur indirectly through microbial ecological modules in karst river, and the denitrification and desulfurization processes were potentially inhibited. These findings contributed to understand the changes and its assembly mechanism of microbial community, as well as the feedback to environment in polluted karst river.

1. Introduction

Approximately 15% of the continents are karst terrain (Yuan, 1997). The water cycle was rapid in the karst system because of its underground-surface double-layered structure and high hydraulic conductivity (Torres et al., 2018). Due to the uneven distribution of karst aquifer fissures and tubes, rapid discharge, and significant seasonal flow velocity variation, the hydrological processes and hydrochemical characteristics exhibit substantial spatial and temporal variability (He et al., 2018; Sun et al., 2019, 2021). Karst water provided drinking water for about 25% of the global population (Dodgen et al., 2017). With the strengthening human activities, however, pollutants such as nitrogen, phosphorus, heavy metals, and organic matter (Wu et al., 2018; Huang F. et al., 2019; Qin et al., 2020, 2021) are continuously introduced into karst aquifers, which have altered the karst ecosystem and element cycles.

As the most active component of the aquatic ecosystem, microorganisms acquire nutrients and energy to survive in an oligotrophic environment via numerous metabolic pathways, which are crucial to the elemental cycling of the karst ecosystem (Webster et al., 2018; Xue et al., 2020; Zhu et al., 2022). According to several studies, Proteobacterium, Bacteroides, and Firmicutes were the most common phyla in karst aquifers (Li et al., 2017; Danczak et al., 2018; Tang et al., 2019; Li Q. et al., 2020), and there were shifts in both geographical and seasonal composition of microbial communities (Opalički Slabe et al., 2021). Deterministic and stochastic processes jointly govern the assembly of microbial communities (Chave, 2004). Previous studies have demonstrated that the relationship between these processes varies across different spatial and temporal scales due to the intensity of environmental change and the threshold of microbial tolerance (Liu et al., 2019). In karst aquifer, the relatively steady hydrodynamic circumstances during the dry season resulted in the distinct microbial communities in individual aquifers because of their distinct physicochemical conditions (Shabarova, 2013; Shabarova et al., 2014), whereas the rapid recharge during the rainy season enabled the passive migration of isolated species in densely connected aquifers, reducing the environmental selection on microbial communities and resulting in their convergence in karst aquifer (Yan et al., 2020). In addition to the aquifer connectivity, different hydrochemical conditions affected the composition and distribution of “permanent inhabitants” and specific species in different seasons in microbial communities. In the unexplored karst aquifer, pH, temperature (Yun et al., 2016; Li H. et al., 2021), salinity (Park et al., 2021; Fang et al., 2022), and organic matter (Iker et al., 2010) were the most impacted factors on the microbial communities. With increasing human activities, the constant input of nitrogen, phosphorus and emerging pollutants such as antibiotics had disrupted the physicochemical properties, the structure and interactions of microbial communities, as well as their carbon and nitrogen cycles (Zhang et al., 2013; Xi et al., 2015; Mahana et al., 2016; Zhernakova et al., 2016; Subirats et al., 2019; Roose-Amsaleg et al., 2021). Not only had antibiotics affected the microbial community structure, but also influenced their metabolic mechanisms, such as carbon sequestration, nitrification and denitrification (Gonzalez-Martinez et al., 2014; Xiong et al., 2015; Li H. et al., 2020; Slipko et al., 2021). Different bacteria acquired resistance genes variously under the selection pressure caused by antibiotics, which regulated their metabolic strategies and the efficiency of element cycles (Ming et al., 2020; Yang et al., 2020; Zhang et al., 2022). However, karst aquifers with intense human activities were often influenced by more than one pollutant (Zhu et al., 2019), which would have more complex effects on the community structure and function of microbial community.

The elements cycles were driven by biogeochemical processes in karst system, which would further influence the role of karst system in global climate change (Mahler et al., 2021). However, the microbial assembly and ecological roles in polluted karst rivers at annual and seasonal scales, as well as their feedback to environment were less studied. In this study, we focused on: (1) the physicochemical properties and pollution characteristics of typical karst river in dry and rainy seasons; (2) the variation and assembly mechanisms of microbial communities in polluted karst river; and (3) the response of microbial community composition and functions to environmental factors, especially the pollutants from human activities.

2. Materials and methods

2.1. Study area and sample collection

The study area, Kaiyang, is situated in Guizhou Province, southwest China, in a typical subtropical monsoon climate zone with an annual average temperature of 14.4°C and precipitation primarily concentrated from June to September. There were two karst rivers from west to east, the south river and the north river, with three underground entrances (K4, K8, K9) and three underground exits (K5, K7, and KX) (Figure 1). A total 40 samples were collected. Each 10 samples were collected in January 2017 (Jan-2017), August 2017 (Aug-2017) (Xiang et al., 2020), August 2018 (Aug-2018) and January 2019 (Jan-2019) (Zhang et al., 2021), respectively.

FIGURE 1

Figure 1. The sampling sites in study area.

Except for K9 and KX, all sampling sites were in the south river. K1 to K4 were located in the surface river from upstream; K1 was located in a reservoir; K2 and K3 were situated upstream and downstream of a waste water treatment plant (WWTP), respectively; and K4 was situated in a resort. There were agricultural fields and a small amount of livestock activity between K5 and K8. K5 was an outlet of underground river, with evident domestic wastes. K8 was an underground entrance. K6 and K7 were located in the upstream tributaries of K5 and K8, respectively. K6 was a surface river site with domestic wastes, and K7 was the outlet of an underground river. K9 was a sinkhole in the scenic region of the north river, while KX was the confluence of the north and south rivers.

2.2. Hydrochemistry monitoring and analyses

Temperature (T) and pH value of karst river water samples were determined by pH meter (pH30, CLEAN, CA, United States). Electronic conductivity (EC) and dissolved oxygen (DO) were determined by EC meter (CON, CLEAN, CA, United States) and DO meter (CON30, CLEAN, CA, United States), respectively. The total organic carbon (TOC) was determined by TOC analyzer (Shimadzu, Japan). HCO₃⁻ was tested by potentiometric titrator (Metrohm 877 Titrino plus, Swiss), and Cl⁻, SO₄²⁻ were measured by ion chromatograph (DIONEX ICS-900, Sunnyvale, CA, United States). Ca²⁺, Mg²⁺, K⁺, and Na⁺ were measured by ICP-OES (SPECTRO Blue Sop, Germany). And the spectrophotometer UV-1800 (Shimadzu, Japan) was used to determine dissolved inorganic nitrogen (DIN) - NH₄⁺, NO₃⁻ and NO₂⁻. Following our previous studies (Huang H. et al., 2019), auto-solid phase extraction (SPE) instrument (Auto SPE-06C, Reeko Instrument, TX, United States) with an Oasis HLB SPE column (6 mL, 500 mg, Waters, MA, United States) was used to extract antibiotics from water samples of karst river. Samples were collected and immediately chilled in 4°C and transported via overnight express to the laboratory for further test. The sampling process was followed our previous study (Xiang et al., 2020).

2.3. DNA extraction and 16S rRNA gene sequencing

The bacteria in each sample were filtered through a 0.22 μm Millipore GSWP membrane for 1.0 l water and then suspended in 10 mL of physiological saline for 6 h. The mentioned suspensions were separated and kept for DNA extraction using a PowerSoil DNA Extraction Kit (MoBio Laboratories, Carlsbad, CA, United States). The extracted DNA was stored at −20°C for further analysis and at −80°C for permanent preservation. Using a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, United States) and agarose gel electrophoresis (Bio-Rad, Hercules, CA, United States), the quantity and quality of the extracted DNA were established.

The primers 338F and 806R (338F: 5′-barcode-ACTCCTAC GGGAGGCAGCAG-3′ and 806R: 5′-barcode-GGACTACHVGG GTWTCTAAT-3′) were used to amplify the V3-V4 regions of the bacterial 16S ribosomal RNA gene (16S rRNA gene). The PCR amplification method was referred to our previous study (Zhang et al., 2020). Raw fastq files were demultiplexed, then quality-filtered using QIIME (Version 1.17). All clusters were subsampled by the minimum reads (19,551 reads), and operational taxonomic units (OTUs) were clustered with a 97% similarity cutoff by Uparse.¹ Mothur² was used to compute the alpha diversity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional profiling was predicted using PICRUSt2.³ The clustered sequence datasets in this study were deposited in the NCBI Genbank with the accession number PRJNA908735 and PRJNA936825.

2.4. Statistical analysis

The data were normalized except for pH value. Spearman correlation analysis and its significance test among each factor were completed by “corrplot” package.⁴ The ordinary least squares (OLS) analysis among ammonium, nitrate, sulfate, and other environmental factors were calculated by “car” package.⁵ Principal co-ordinates analysis (PCoA) based on Bray-Curtis distance and variance partitioning analysis (VPA) were both performed by “vegan” package.⁶ Venn diagram at genus level was plotted by “vennDiagram” package. The Kruskal-Wallis test (p < 0.05) was completed by “nparcomp” package, and corrected p-values were calculated using Dunn’s test. The microbial community assembly processes were analyzed by “iCAMP” package.⁷ And the heatmap figure was plotted by “pheatmap” package.⁸ All the above data analysis were implemented in R.⁹ The co-occurrence network of bacteria at genus level was visualized with Gephi 0.9.2,¹⁰ and ecological modules were identified using the Louvain algorithm. The saturation index of gypsum (SIg) and saturation index of calcite (SIc) were calculated by PHREEQC (Version 3.6.2). Structural equation model (SEM) was performed by SPSS AMOS (Version 24.0) to discuss the environmental factors effect on ecological module and element metabolic module.

3. Results and discussion

3.1. Characteristics of physicochemical properties of karst river

The pH value of karst river water was ranging from 7.70 to 8.82, and the temperature was ranging from 4.0°C to 26.1°C (Supplementary Table 1). DO ranged from 1.88 to 9.93 mg/L, and TOC ranged from 1.09 to 30.63 mg/L. Ca²⁺ and HCO₃⁻ were the major cation and anion due to the abundant carbonate rocks in the study area (Sun et al., 2019). For carbonate rock weathering involving only carbonic acid, the equivalent ratio of [Ca²⁺+Mg²⁺]/[HCO₃⁻] is typically around 1. But the average value of our study was 1.23, while the mean value of [Ca²⁺+Mg²⁺]/[HCO₃⁻ + NO₃⁻ + SO₄²⁻] was 0.93, indicating that the carbonate weathering was influenced by both carbonic acid and sulfuric/nitric acids (Martin, 2017; Sun et al., 2021).

The concentration of NO₃⁻ varied from 0.74 to 45.82 mg/L, with an average value of 18.23 mg/L, and the concentration of SO₄²⁻ varied from 20.14 to 74.25 mg/L. There were severe NH₄⁺ and NO₂⁻ pollution, with NH₄⁺ amounting from 0.02 to 45.35 mg/L and NO₂⁻ amounting from 0.01 to 9.51 mg/L. The concentrations of both NH₄⁺ and NO₃⁻ showed seasonal variations with higher in winter and lower in summer, mainly influenced by the dilution effect of precipitation recharge in rainy season. Meanwhile, the NH₄⁺ were higher in upstream, particularly from K1 to K4, while NO₃⁻ concentrated downstream from K7 to KX, due to nitrification and the NO₃⁻ accumulation for its high solubility (Li X. et al., 2021; Li H. et al., 2021). The negative correlations between NH₄⁺ and NO₃⁻ also supported this hypothesis (Eqs 1, 2):

$\begin{array}{l}{\text{NH}}_4^{+}=-0.35{\text{NO}}_3^{-}+0.17\text{TOC}-0.31\mathrm{T}\\ +0.20\text{DO}+0.50\text{Antibiotics}\\ +\mathrm{a}1\left(\mathrm{R}2=0.56,p<0.01,\mathrm{a}1<0.0001\right)\end{array}\text{ (Equation 1)}$

$\begin{array}{l}{\text{NO}}_3^{-}=-0.42{\text{NH}}_4^{+}+0.54{\text{NO}}_2^{-}-0.60\mathrm{T}+0.27\text{pH}\\ +\mathrm{a}2\mathrm{\hspace{0.33em}}\left({\mathrm{R}}^2=0.48,p<0.01,{\mathrm{a}}^2<0.0001\right).\end{array}\text{ (Equation 2)}$

As shown in Eq. 1, The concentration of NH₄⁺ was positively correlated with antibiotics, as well as its correlation coefficients with macrolides, lincomycin, tetracyclines were 0.78, 0.71 and 0.66 (p < 0.01), respectively. The concentrations of NH₄⁺ and antibiotics showed a similar spatial variation trend, displayed higher upstream. As one of the most prevalent forms of nitrogen pollution in WWTP, NH₄⁺ originated primarily from wastewater plant discharges, the same source of antibiotics.

SO₄²⁻ is mainly derived from both natural (dissolution of soluble sulfate, oxidation of sulfide minerals, and atmospheric precipitation) and anthropogenic inputs(agricultural fertilizers, domestic sewage, industrial wastewater, and mine waste water)(Li et al., 2018). In the study area, the [Mg²⁺]/[Ca²⁺] values ranged from 0.39 to 1.14 with average 0.65, which indicated that the chemical composition is mainly influenced by low-magnesium minerals such as calcite and gypsum. Furthermore, the SIg values were − 2.43 ~ −1.75 while SIc were − 0.79 ~ 1.28, implied that the dissolution of gypsum were the main source of sulfate. Additional, SO₄²⁻ were positively correlated with antibiotics (R = 0.66, p<0.01), which implied that SO₄²⁻ also originated from WWTP, which was consistent with its concentration spatial distribution of higher in upstream.

A total of 28 antibiotics were detected, with total concentrations ranging from 9.13 to 1411.33 ng/L (Supplementary Table 2). The mean concentrations of four sampling times from high to low were Jan-2019 (average value of 514.15 ng/L), Jan-2017 (472.60 ng/L), Aug-2017 (392.44 ng/L), and Aug-2018 (199.37 ng/L), with a seasonal trend of higher in winter and lower in summer (Figure 2), due to the dilution effect of precipitation recharge, which was consistent with our previous study (Huang F. et al., 2019; Huang H. et al., 2019). The maximum value for each time was frequently recorded at K2 and K3, suggesting that WWTP effluent discharge was the major source of antibiotics in karst river, which indicated the anthropogenic inputs. Antibiotic demonstrated a spatial changes that fluctuates and decreases along the stream, which was affected by the adsorption of soil and sediment and the degradation of bacteria (Zhao et al., 2012). Moreover, antibiotic concentrations in tributaries K6 and K7 were lower than those in the mainstream.

FIGURE 2

Figure 2. The concentrations of main antibiotics in karst river at four sampling times (Jan-2017, Aug-2017, Aug-2018 and Jan-2019). (A–F) The concentrations of quinolones, lincomycin, macrolides, tetracyclines, sulfonamides and Chloramphenicol.

Quinolones (QNs), lincomycin (LIN), and macrolides (MLs) were the major antibiotics, accounting on average for 43.45, 34.58, and 13.99% of total concentration, respectively, which were consistent with their frequent use in southwest China (Sun et al., 2015; Zeng et al., 2019). QNs ranged from 3.10 to 994.13 ng/L (Figure 2A) with a 95% detection frequency. QNs were persistent pollutants due to their high chemical stability and prolonged half-life (Andriole, 1999). Ofloxacin and norfloxacin, as the maximum content subclass of QNs in most WWTP (Dinh et al., 2017), also dominated in karst river, ranging from 4.60 to 486.00 ng/L and 4.09 to 204.00 ng/L, respectively. The detection frequency of LIN was 95%, and the concentration was 4.77–227.53 ng/L (Figure 2B). The detection frequency of MLs was 88%, with concentration ranging from 0.83 to 205.13 ng/L (Figure 2C). Erythromycin and roxithromycin were the main subclasses of MLs. While antibiotics were transferred underground by continuous precipitation recharge in summer, their natural degradation was weakened due to weak light and then accumulated (Liu et al., 2021). However, the high connectivity of karst aquifers enabled rapid hydration and a more quickly response to contaminants.

The concentrations of tetracyclines (TCs), chloramphenicols (CAP), and sulfonamides (SAs) were ranged from 3.43 to 177.80 ng/L (Figure 2D), 1.70 to 76.20 ng/L (Figure 2E), and 0.65 to 33.56 ng/L (Figure 2F), respectively. The detection rate of TCs was 60%, with tetracycline and oxytetracycline being the most commonly used antibiotics. The detection frequency of CAP was 53%. Both TCs and CAP were higher in winter. The detection rate of SAs was 75%, with sulfapyridine being the most common subclass. The maximum content and the most subclasses of SAs were observed in Aug-2018. Only samples from Aug-2018 contained sulfadimethoxine and sulfadimidin, both of which were widely applied in veterinary medicine. Increased livestock vaccination in study area in 2018 may result in an increase use of veterinary drugs.¹¹ As a notable veterinary drug, enrofloxacin was also detected only in Aug-2018.

3.2. The structure of bacterial community in karst river

A total of 3,955,799 optimum sequences were obtained in 40 samples (Supplementary Table 3). After subsampling by the minimum sequences, 5,242 OTUs were obtained. The species richness in Aug-2017 was the lowest. Nonetheless, the Shannon index had no significant difference at four sampling times (Supplementary Table 3).

Proteobacteria (relative abundance ranged from 36.03 to 61.15%), Bacteroidetes (18.97 to 38.03%), Cyanobacteria (6.14 to 13.85%), and five other phyla with relative abundance greater than 5% were the dominant phyla. Flavobacterium, Arcobacter, norank_c_Cyanobacteria, Pseudomonas and Acinetobacter were prevalent in a total 945 genera, with the relative abundance of 12.68, 5.78, 5.18, 4.41 and 4.35%, respectively.

PC1 and PC2 explained 22.65% of the variance in community structure together (Figure 3A). The bacterial communities of Jan-2017 and Jan-2019 were similar. The bacterial community structures of Aug-2017 and Aug-2018 were divided from each other, which apparently differed from samples taken in winter.

FIGURE 3

Figure 3. Differences of microbial communities in karst rivers among four sampling times (Jan-2017, Aug-2017, Aug-2018 and Jan-2019). (A) Principal co-ordinates analysis (PCoA) of microbial community structure based on Bray-Curtis distance at genus level. (B) Venn diagram of microbial community at genus level. (C) Significantly Different of genus among two winter, Aug-2017 and Aug-2018. Kruskal-Wallis test with p-values were *p < 0.05, **p < 0.01, ***p < 0.001.

The relative abundance of 275 shared genera varied from 89.04 to 96.27% (Figure 3B). Most of them showed obvious seasonal characteristics. Seven genera were mainly distributed in winter, including Flavobacterium, Pseudomonas and hgcI_clade. In previous study, we observed strong correlations between ARGs and Flavobacterium, Pseudomonas (Zhang et al., 2021), which were more likely to survive in winter with higher antibiotic concentrations. Six genera were enriched in summer, including Dechloromonas, Hydrogenophaga and Cloacibacterium. The relative abundance of unclassified_f__Comamonadaceae gradually increased while Limnohabitans showed no noticeable trend in this study. Nine out of the nineteen dominant genera (relative abundance above 1%) (Supplementary Figure 1) exhibited statistically significant differences (p < 0.05) (Figure 3C). Flavobacterium, Pseudomonas, hgcI_clade and Albirhodobacter were considerably enriched in winter with average relative abundances varying from 8.15 to 28.06%, 6.26% to 8.08, 3.53 to 4.54% and 1.49 to 1.62%, respectively. Acinetobacter, Malikia, Dechloromonas and Hydrogenophaga were highly expressed in Aug-2017 with average relative abundances of 12.37, 11.90, 5.24 and 1.94%, respectively. And the relative abundance of norank_c_Cyanobacteria (7.15%) was the highest in Aug-2018.

Compared with shared genera, the richness and abundance of specific genera were lower. Forty seven genera only presented in winter, with relative abundances of 0.23 and 0.24% in Jan-2017 and Jan-2019, including Candidatus_Azambacteria, Empedobacter, etc. The winter genera were mainly distributed at K4 and K6 in Jan-2017 and at K6 in Jan-2019. As a sampling site in the tributary, there might be a relatively higher flow at K6 that introduce a great number of exotic species in winter. Forty six genera were only found in Aug-2017, with an average relative abundance of 1.13%, and distributed mainly at K2 and K3, such as Macromonas, Candidatus_Aquiluna. Forty nine genera such as Sporacetigenium and Desulfitobacterium were only found in Aug-2018, with an average relative abundance of 1.73% and primarily enriched at K8.

Five clusters (MOD0-MOD4) were obtained through co-occurrence network analysis of genera with relative abundance above 0.1% (Figure 4A). Twenty one genera were included in MOD0, primarily Flavobacterium, Malikia, Dechloromonas, etc., each of which exhibited distinct seasonal differences. Thus, there was no discernible seasonal trend in MOD0 but a peak value in Jan-2019 (Figure 4C), where Flavobacterium concentrated. Twenty three genera of MOD1 were mainly included Limnohabitans and Planktothrix, which concentrated in summer. Twenty five genera were included in MOD2, such as Arcobactor, Pseudomonas, Cyanobacteria. MOD3 was consisted of 17 genera, such as hgcI_clade and Hydrogenophaga, etc., which were abundant in winter. Prosthecobacter, Nitromonas and Nitrospira were dominant in MOD4 among the total 24 genera.

FIGURE 4

Figure 4. Interactions of microbial community in karst river. (A) Co-occurrence network of genus with relative abundance above 0.1%. Genus with relative abundance above 0.1% were selected to calculate their correlations. Each node represented a genus, the nodes were colored by module, and the connection between nodes represented a significant correlation (Spearman, p < 0.05, |r| > 0.6). (B) Connections of modules. The correlation of each genus in every module, negative relations were represented in green edges and positive relations were represented in orange edges. (C) Stacked histogram of five modules in Jan-2017, Aug-2017, Aug-2018, and Jan-2019. Relative abundance of each module was represented in different color.

The co-occurrence network of major microorganisms helped to reveal functional group structure and potential interactions between bacteria (Curtis et al., 2002; Cardona et al., 2016). There were 165 positive and 68 negative edges among five MODs (Figure 4B), with a modularity of 0.41, suggesting a strongly associated community among the major communities which were mainly positive synergies. MOD0, MOD3 and MOD4 were strongly connected with other MODs with 100, 98 and 65 positive edges, respectively. There were more negative edges between MOD1 and MOD4, indicating a negatively competitive relationship among genera. Genera involved in carbon, nitrogen and sulfur metabolism were widely distributed in each module. Flavobacterium, Dechloromonas, Massilia in MOD0 participated in denitrification, nitrogen reduction and sulfate reduction (Tanikawa et al., 2018; Chang et al., 2019). The diversity of metabolic function makes it most closely related to other MODs. MOD3 were primarily denitrifying bacteria, such as hgcI_clade (Herrmann et al., 2015; Luo et al., 2018). Conversely, ecological modules with different functions often exhibit a tendency to compete for limiting nutrients (Tate, 2020). This competitive relationship may also exist between MOD1 and MOD4. MOD1 was dominated by Limnohabitans, Planktothrix, and other genera with both carbon and nitrogen metabolism function, whereas the MOD4 was ruled by nitrifying bacteria such as Prosthecobacter, Nitrosomonas, and Nitrospira. In addition, Hydrogenophaga and Limnohabitans in MOD1 and MOD4 participate in sulfur metabolism (Herrmann et al., 2015; Luo et al., 2018). Conjecturally, bacterial communities participated in nitrogen and sulfur cycles were potentially affected by the pollution of antibiotics and nitrogen in karst river.

3.3. Assembly mechanisms of bacterial community in karst river

The bacterial community assembly processes were consisted of deterministic processes (heterogeneous selection and homogeneous selection) and stochastic processes (dispersal limitation, homogenizing dispersal and drift) (Stegen et al., 2013). Deterministic processes represented the influence of environmental factors such as physical and chemical indicators on microbial community, while stochastic processes were related to spatial variation. In general, deterministic processes were dominant in small spatial scale while stochastic processes dominated in a larger spatial scale (Dini-Andreote et al., 2015). However, frequent anthropogenic impacts and intense variation in karst hydrodynamic conditions produced high variance and increased the influence of stochastic processes.

In Figure 5, stochastic processes had more advantage in Aug-2017 (56.96%) and Aug-2018 (50.66%), and primarily affected by dispersal limitation and drift effects. Stochastic processes played a major role in microbial community assembly with less environmental pressure (Dini-Andreote et al., 2015). Lower antibiotic concentrations in summer reduced the survival pressure for less resistant bacteria, which increased their relative abundance. The type of antibiotics in Aug-2017 and Aug-2018 were remarkable different, and their additive, synergistic and antagonistic effects led to greater differences in community structure (De Liguoro et al., 2018; Liang et al., 2022). The rapid surface-underground exchange and the drastic environmental changes in karst river resulted in short adaptation and rapid turnover in summer, which benefited distinguishable bacterial communities as well. Deterministic processes enhanced in winter, and homogeneous selection was predominant, with 49.51 and 68.14% in Jan-2017 and Jan-2019, respectively, which was caused by the survival pressure of high concentration of antibiotics and DIN. Therefore, we suggested that the environmental selection caused by the multiple pressure of antibiotics and DIN created a more stable bacterial community with antibiotic resistance in winter. Consequently, despite the two-year sampling interval between Jan-2017 and Jan-2019, the structure of bacterial community showed a highly similarity, while the stronger survival pressure from pollutions led to a decline in the relative abundance of specific genera. The proportion of dispersal limitation in the stochastic processes decreased gradually, while the proportion of homogenizing dispersal increased, indicating that the composition of microbial community tended to be more similar over time.

FIGURE 5

Figure 5. Assembly processes of bacteria community in karst river at four sampling times (Jan-2017, Aug-2017, Aug-2018, and Jan-2019) in karst river. Deterministic processes were composed of heterogeneous selection and homogeneous selection, and stochastic processes were composed of dispersal limitation, homogenizing dispersal, drift and other effects.

3.4. Effects of environmental factors on the functions of bacterial community in karst river

Variations in microbial community were accompanied by alterations in community function (Zhong et al., 2018). The results of KEGG functional profiling revealed that the average relative abundance of 16 functional genes exceeded 1%, which were mainly sub-functions of metabolism, environmental information processing, and genetic information processing and so on (Figure 6). Among them, the average relative abundance of energy metabolism ranged from 5.58 to 5.79% and was higher in summer, possibly due to nutrient limitation (Cline and Zak, 2015; Jiang et al., 2022). In winter, the copious nutrients conditions made nutrients more accessible to microorganisms, and hence the corresponding metabolism gene was less prevalent. However, nutrient availability was restricted by dilution effects in summer, and microorganisms maintained their growth by increasing the metabolic function gene (Jiang et al., 2022).

FIGURE 6

Figure 6. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional profiling at level-2 at four sampling time (Jan-2017, Aug-2017, Aug-2018 and Jan-2019).

From the results of function profiling, 13 carbon fixation pathways, 5 nitrogen metabolism pathways, and 2 sulfur metabolism pathways were annotated (Figure 7; Supplementary Table 4). The carbon fixation pathways mainly included 7 carbon fixation pathways in prokaryotes (average relative abundance ranged from 5.98 to 6.58%) and 8 pathways in photosynthetic organisms (from 3.01 to 3.26%), such as reductive citrate cycle (M00173), dicarboxylate-hydroxybutyrate cycle (M00374), and reductive pentose phosphate cycle (Calvin cycle) (M00165), etc. The nitrogen metabolism pathways included nitrogen fixation (M00175), nitrification (M00528), denitrification (M00529), and 2 nitrate reduction to ammonium pathways (M00530 and M00531), with average relative abundances from 0.43 to 0.81%. The sulfur metabolism pathways included 2 sulfate reduction pathways (M00176 and M00596), with relative abundance ranging from 0.55 to 0.66%. But all the metabolism pathways above showed no clear seasonal variations.

FIGURE 7

Figure 7. The stacked histogram of carbon, nitrogen, sulfur metabolic modules of Kyoto Encyclopedia of Genes and Genomes (KEGG) functional profiling at four sampling time (Jan-2017, Aug-2017, Aug-2018 and Jan-2019).

SEM showed that MOD0, MOD3 and MOD4 contributed significantly to nitrogen and sulfur metabolism (Figure 8A). MOD0 was mainly composed of denitrification and nitrate reduction bacteria, so the denitrification was significantly influenced directly by MOD0, with path coefficient value of 0.95, and the nitrate reduction was not significantly influenced directly by MOD0. MOD3 also directly contributed significantly to nitrate reduction and denitrification, with coefficients of 0.84 and 0.44, respectively. Meanwhile, MOD4 played an important role in sulfate reduction.

FIGURE 8

Figure 8. The environmental impacts on bacterial communities and their metabolic function in karst river. (A) Structure equation model (SEM) among environmental factors, modules and metabolic functions. Antibiotics were including tetracyclines, macrolides, chloramphenicol and quinolones. Green lines, orange lines indicated negative and positive impacts, respectively, and the associated numbers were the path coefficient value with p-values *p < 0.05, **p < 0.01, ***p < 0.001. (B) Variance partitioning analysis (VPA) diagram. The impact of antibiotics, DIN and hydrochemical factors on microbial communities. DIN included ammonium, nitrate and nitrite. Antibiotics included quinolones, lincomycin, tetracyclines, macrolides, sulfonamides and chloramphenicol. Hydrochemistry included temperature, pH, electronic conductivity, dissolved oxygen and total organic carbon.

On the other hand, the microbial community function was indirectly affected by environmental factors through certain ecological modules. Consistent with the result of VPA (Figure 8B), the antibiotics and DIN had significant impacts on microbial communities (with 5.34 and 0.53% explanation, respectively). MOD3 and MOD4 were resistant and promoted by the antibiotics. Trichococcus and norank_f__Leptotrichiaceae in MOD3 and Acidovorax and norank_p_Saccharibacteria in MOD4 carried resistance genes of tetracycline and multidrugs (Zhang et al., 2021), providing a stronger survival advantage in winter with higher concentration of antibiotics. DIN showed negative influence to MOD0, MOD3 and MOD4. The high nitrate content would produce significant inhibition or toxicity to bacteria (Cua and Stein, 2011), as well as to sulfate-reducing bacteria (He et al., 2010; Zhou and Xing, 2021). Besides, the high nitrate suppression to sulfate-reducing bacteria was diminished by the denitrifier (Zhou and Xing, 2021), which also explained the strong correlation between MOD4 and MOD0.

Denitrification, nitrate reduction and sulfur reduction were indirectly promoted by antibiotics via MOD3 and MOD4, whereas they were inhibited by DIN. As the substrate of nitrogen metabolism, DIN also indirectly promoted the process of nitrification, denitrification and nitrate reduction. SEM results showed that the total effect value of environmental factors on nitrification and nitrate reduction was −2.53 and 20.39, respectively. We suggested that the antibiotics and DIN inhibited nitrification, but strongly promoted the nitrate reduction, which lead a further accumulation of NH₄⁺. Although the effects of antibiotics and inorganic nitrogen pollution on microbial ecological modules and their functions in karst rivers was indicated, the influence relationship of different ecological groups on environmental factors and the possible environmental utility was identified, the influence on metabolic process and the fluxes of nitrogen and sulfur still requires further researches such as long-term monitoring, isotope tests and other experimental methods.

4. Conclusion

The structure of microbial community with physicochemical properties altered considerably from winter to summer in karst ecosystem. Pollution such as antibiotics and inorganic nitrogen were introduced by anthropogenic activities from WWTP and agricultural discharges. Microbial communities in winter of 2 years interval were stable with higher similarity than those in summer of different year. Deterministic processes, homogeneous selection process specifically, occupied in winter due to the environmental pressure of high pollutants content, while stochastic processes had more advantages in summer with a significant extent of spatial heterogeneity due to intensive human activities and karst hydrodynamic condition. Pollution would indirectly affect the nitrogen and sulfur metabolism via microbial functional group. Furthermore, the environmental response and feedback of microbial communities’ element cycles needs more attentions and further measurable research.

Data availability statement

The datasets presented in this study can be found in the NCBI repository (https://www.ncbi.nlm.nih.gov/), accession numbers PRJNA908735 and PRJNA936825.

Author contributions

XG: contributed to the conceptualization, methodology, and writing. RH and BZ: participated in the formal analysis and discussion. CG: participated in the investigation. FL: participated in the topic and discussion of the whole study. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the National Natural Science Foundation of China (Grant no. 42172336 awarded to XG), Guangxi key R & D program support (Grand no. Guike AB22080070 awarded to FL) and the Fundamental Research Funds for Central Universities (Grant no. 2652019077 awarded to XG).

Conflict of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

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.2023.1133938/full#supplementary-material

Footnotes

1. ^http://www.drive5.com/uparse/

2. ^http://www.mothur.org/wiki/Classify.seqs

3. ^https://github.com/picrust/picrust2/

4. ^https://github.com/taiyun/corrplot

5. ^https://callr.r-lib.org/

6. ^https://github.com/vegandevs/vegan

7. ^https://github.com/DaliangNing/iCAMP1

8. ^https://cran.rstudio.com/package=pheatmap

9. ^https://www.r-project.org

10. ^https://gephi.org/

11. ^http://nynct.guizhou.gov.cn/xwzx/zwdt/201712/t2017120525041433.html

References

Andriole, V. T. (1999). The Future of the Quinolones. Drugs 58, 1–5. doi: 10.2165/00003495-199958002-00001

CrossRef Full Text | Google Scholar

Cardona, C., Weisenhorn, P., Henry, C., and Gilbert, J. A. (2016). Network-based metabolic analysis and microbial community modeling. Curr. Opin. Microbiol. 31, 124–131. doi: 10.1016/j.mib.2016.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, M., Wang, Y., Pan, Y., Zhang, K., Lyu, L., Wang, M., et al. (2019). Nitrogen removal from wastewater via simultaneous nitrification and denitrification using a biological folded non-aerated filter. Bioresour. Technol. 289:121696. doi: 10.1016/j.biortech.2019.121696

PubMed Abstract | CrossRef Full Text | Google Scholar

Chave, J. (2004). Neutral theory and community ecology: Neutral theory and community ecology. Ecol. Lett. 7, 241–253. doi: 10.1111/j.1461-0248.2003.00566.x

CrossRef Full Text | Google Scholar

Cline, L. C., and Zak, D. R. (2015). Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession. Ecology 96, 3374–3385. doi: 10.1890/15-0184.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cua, L. S., and Stein, L. Y. (2011). Effects of nitrite on ammonia-oxidizing activity and gene regulation in three ammonia-oxidizing bacteria: Effects of nitrite on ammonia-oxidizing bacteria. FEMS Microbiol. Lett. 319, 169–175. doi: 10.1111/j.1574-6968.2011.02277.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Curtis, T. P., Sloan, W. T., and Scannell, J. W. (2002). Estimating prokaryotic diversity and its limits. Proc. Natl. Acad. Sci. U. S. A. 99, 10494–10499. doi: 10.1073/pnas.142680199

PubMed Abstract | CrossRef Full Text | Google Scholar

Danczak, R. E., Johnston, M. D., Kenah, C., Slattery, M., and Wilkins, M. J. (2018). Microbialcommunity cohesion mediates community turnover in unperturbed aquifers. mSystems, 3, e00066–e00018. doi: 10.1128/mSystems.00066-18

PubMed Abstract | CrossRef Full Text | Google Scholar

De Liguoro, M., Riga, A., and Fariselli, P. (2018). Synergistic toxicity of some sulfonamide mixtures on Daphnia magna. Ecotoxicol. Environ. Saf. 164, 84–91. doi: 10.1016/j.ecoenv.2018.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinh, Q., Moreau Guigon, E., Labadie, P., Alliot, F., Teil, M. J., Blanchard, M., et al. (2017). Fate of antibiotics from hospital and domestic sources in a sewage network. Sci. Total Environ. 575, 758–766. doi: 10.1016/j.scitotenv.2016.09.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Dini-Andreote, F., Stegen, J. C., van Elsas, J. D., and Salles, J. F. (2015). Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. Proc. Natl. Acad. Sci. U. S. A. 112, E1326–E1332. doi: 10.1073/pnas.1414261112

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodgen, L. K., Kelly, W. R., Panno, S. V., Taylor, S. J., Armstrong, D. L., Wiles, K. N., et al. (2017). Characterizing pharmaceutical, personal care product, and hormone contamination in a karst aquifer of southwestern Illinois, USA, using water quality and stream flow parameters. Sci. Total Environ. 578, 281–289. doi: 10.1016/j.scitotenv.2016.10.103

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Y., Liu, J., Yang, J., Wu, G., Hua, Z., Dong, H., et al. (2022). Compositional and metabolic responses of autotrophic microbial community to salinity in lacustrine. Environments 7:e0033522. doi: 10.1128/msystems.00335-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Martinez, A., Rodriguez-Sanchez, A., Martinez-Toledo, M. V., Garcia-Ruiz, M.-J., Hontoria, E., Osorio-Robles, F., et al. (2014). Effect of ciprofloxacin antibiotic on the partial-nitritation process and bacterial community structure of a submerged biofilter. Sci. Total Environ. 476-477, 276–287. doi: 10.1016/j.scitotenv.2014.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Q., He, Z., Joyner, D. C., Joachimiak, M., Price, M. N., Yang, Z. K., et al. (2010). Impact of elevated nitrate on sulfate-reducing bacteria: a comparative study of Desulfovibrio vulgaris. ISME J. 4, 1386–1397. doi: 10.1038/ismej.2010.59

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Z., Liang, H., Yang, C., Huang, F., and Zeng, X. (2018). Temporal–spatial evolution of the hydrologic drought characteristics of the karst drainage basins in South China. Int. J. Appl. Earth Obs. Geoinform. 64, 22–30. doi: 10.1016/j.jag.2017.08.010

CrossRef Full Text | Google Scholar

Herrmann, M., Rusznyák, A., Akob, D. M., Schulze, I., Opitz, S., Totsche, K. U., et al. (2015). Large Fractions of CO ₂ -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

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H., Chen, Z., Wang, T., Xiang, C., Zhang, L., Zhou, G., et al. (2019). Nitrate distribution and dynamics as indicators to characterize karst groundwater flow in a mined mineral deposit in southwestern China. Hydrogeol. J. 27, 2077–2089. doi: 10.1007/s10040-019-01987-0

CrossRef Full Text | Google Scholar

Huang, F., Zou, S., Deng, D., Lang, H., and Liu, F. (2019). Antibiotics in a typical karst river system in China: Spatiotemporal variation and environmental risks. Sci. Total Environ. 650, 1348–1355. doi: 10.1016/j.scitotenv.2018.09.131

PubMed Abstract | CrossRef Full Text | Google Scholar

Iker, B. C., Kambesis, P., Oehrle, S. A., Groves, C., and Barton, H. A. (2010). Microbial Atrazine Breakdown in a Karst Groundwater System and Its Effect on Ecosystem Energetics. J. Environ. Qual. 39, 509–518. doi: 10.2134/jeq2009.0048

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Sun, X., Feng, J., Zhu, S., and Shui, W. (2022). Metagenomic analysis reveals the different characteristics of microbial communities inside and outside the karst tiankeng. BMC Microbiol. 22:115. doi: 10.1186/s12866-022-02513-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Gao, D., Hou, L., Qian, W., Liu, M., Zeng, H., et al. (2021). Nitrogen loads alter the N2 production between denitrification and anammox in Min River Estuary, a highly impacted estuary in southeast China. Environ. Pollut. 277:116757. doi: 10.1016/j.envpol.2021.116757

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Huang, Y., Xin, S., and Li, Z. (2020). Comparative analysis of bacterioplankton assemblages from two subtropical karst reservoirs of southwestern China with contrasting trophic status. Sci. Rep. 10:22296. doi: 10.1038/s41598-020-78459-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Song, A., Peng, W., Jin, Z., Müller, W. E. G., and Wang, X. (2017). Contribution of aerobic anoxygenic phototrophic bacteria to total organic carbon pool in aquatic system of subtropical karst catchments, Southwest China: evidence from hydrochemical and microbiological study. FEMS Microbiol. Ecol. 93, fix065. doi: 10.1093/femsec/fix065

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Song, H.-L., Xu, H., Lu, Y., Zhang, S., Yang, Y.-L., et al. (2020). Effect of the coexposure of sulfadiazine, ciprofloxacin and zinc on the fate of antibiotic resistance genes, bacterial communities and functions in three-dimensional biofilm-electrode reactors. Bioresour. Technol. 296:122290. doi: 10.1016/j.biortech.2019.122290

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Wu, P., Han, Z., Zha, X., Ye, H., and Qin, Y. (2018). Effects of mining activities on evolution of water quality of karst waters in Midwestern Guizhou, China: evidences from hydrochemistry and isotopic composition. Environ. Sci. Pollut. Res. 25, 1220–1230. doi: 10.1007/s11356-017-0488-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Yang, S., Semenov, M. V., Yao, F., Ye, J., Bu, R., et al. (2021). Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community. Glob. Change Biol. 27, 2763–2779. doi: 10.1111/gcb.15593

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, D., Hu, Y., Huang, R., Cheng, J., and Chen, Y. (2022). Effects of various antibiotics on aerobic nitrogen removal and antibiotic degradation performance: Mechanism, degradation pathways, and microbial community evolution. J. Hazard. Mater. 422:126818. doi: 10.1016/j.jhazmat.2021.126818

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Meng, Z., Liu, X., and Zhang, X.-H. (2019). Microbial assembly, interaction, functioning, activity and diversification: a review derived from community compositional data. Mar. Life Sci. Technol. 1, 112–128. doi: 10.1007/s42995-019-00004-3

CrossRef Full Text | Google Scholar

Liu, X., Wang, Z., Zhang, L., Fan, W., Yang, C., Li, E., et al. (2021). Inconsistent seasonal variation of antibiotics between surface water and groundwater in the Jianghan Plain: Risks and linkage to land uses. J. Environ. Sci. 109, 102–113. doi: 10.1016/j.jes.2021.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Tan, X., Liu, K., and Lin, W. (2018). Survey of sulfur-oxidizing bacterial community in the Pearl River water using soxB, sqr, and dsrA as molecular biomarkers. 3 Biotech 8:73. doi: 10.1007/s13205-017-1077-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahana, D., Trent, C. M., Kurtz, Z. D., Bokulich, N. A., Battaglia, T., Chung, J., et al. (2016). Antibiotic perturbation of the murine gut microbiome enhances the adiposity, insulin resistance, and liver disease associated with high-fat diet. Genome Med. 8:48. doi: 10.1186/s13073-016-0297-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahler, B. J., Jiang, Y., Pu, J., and Martin, J. B. (2021). Editorial: Advances in hydrology and the water environment in the karst critical zone under the impacts of climate change and anthropogenic activities. J. Hydrol. 595:125982. doi: 10.1016/j.jhydrol.2021.125982

CrossRef Full Text | Google Scholar

Martin, J. B. (2017). Carbonate minerals in the global carbon cycle. Chem. Geol. 449, 58–72. doi: 10.1016/j.chemgeo.2016.11.029

CrossRef Full Text | Google Scholar

Ming, X., Groves, C., Wu, X., Chang, L., Zheng, Y., and Yang, P. (2020). Nitrate migration and transformations in groundwater quantified by dual nitrate isotopes and hydrochemistry in a karst World Heritage site. Sci. Total Environ. 735:138907. doi: 10.1016/j.scitotenv.2020.138907

PubMed Abstract | CrossRef Full Text | Google Scholar

Opalički Slabe, M., Danevčič, T., Hug, K., Fillinger, L., Mandic-Mulec, I., Griebler, C., et al. (2021). Key drivers of microbial abundance, activity, and diversity in karst spring waters across an altitudinal gradient in Slovenia. Aquat. Microb. Ecol. 86, 99–114. doi: 10.3354/ame01956

CrossRef Full Text | Google Scholar

Park, Y., Yu, J., Nguyen, V. K., Park, S., Kim, J., and Lee, T. (2021). Understanding complete ammonium removal mechanism in single-chamber microbial fuel cells based on microbial ecology. Sci. Total Environ. 764:144231. doi: 10.1016/j.scitotenv.2020.144231

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, W., Han, D., Song, X., and Liu, S. (2021). Sources and migration of heavy metals in a karst water system under the threats of an abandoned Pb–Zn mine. Southwest China. Environ. Pollut. 277:116774. doi: 10.1016/j.envpol.2021.116774

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, Y., Hao, F., Zhang, D., Lang, Y., and Wang, F. (2020). Accumulation of organic carbon in a large canyon reservoir in Karstic area, Southwest China. Environ. Sci. Pollut. Res. 27, 25163–25172. doi: 10.1007/s11356-020-08724-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Roose-Amsaleg, C., David, V., Alliot, F., Guigon, E., Crouzet, O., and Laverman, A. M. (2021). Synergetic effect of antibiotic mixtures on soil bacterial N2O-reducing communities. Environ. Chem. Lett. 19, 1873–1878. doi: 10.1007/s10311-020-01117-3

CrossRef Full Text | Google Scholar

Shabarova, T. (2013). Life in subsurface pools: Insight into microbial diversity and dynamics in the endokarst environment zur Erlangung der naturwissenschaftlichen Doktorwürde

Google Scholar

Shabarova, T., Villiger, J., Morenkov, O., Niggemann, J., Dittmar, T., and Pernthaler, J. (2014). Bacterial community structure and dissolved organic matter in repeatedly flooded subsurface karst water pools. FEMS Microbiol. Ecol. 89, 111–126. doi: 10.1111/1574-6941.12339

PubMed Abstract | CrossRef Full Text | Google Scholar

Slipko, K., Marano, R. B., Cytryn, E., Merkus, V., Wögerbauer, M., Krampe, J., et al. (2021). Effects of subinhibitory quinolone concentrations on functionality, microbial community composition, and abundance of antibiotic resistant bacteria and qnrS in activated sludge. J. Environ. Chem. Eng. 9:104783. doi: 10.1016/j.jece.2020.104783

CrossRef Full Text | Google Scholar

Stegen, J. C., Lin, X., Fredrickson, J. K., Chen, X., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Subirats, J., Di Cesare, A., Varela Della Giustina, S., Fiorentino, A., Eckert, E. M., Rodriguez-Mozaz, S., et al. (2019). High-quality treated wastewater causes remarkable changes in natural microbial communities and intI1 gene abundance. Water Res. 167:114895. doi: 10.1016/j.watres.2019.114895

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, P., He, S., Yu, S., Pu, J., Yuan, Y., and Zhang, C. (2021). Dynamics in riverine inorganic and organic carbon based on carbonate weathering coupled with aquatic photosynthesis in a karst catchment. Southwest China. Water Res. 189:116658. doi: 10.1016/j.watres.2020.116658

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, P., He, S., Yuan, Y., Yu, S., and Zhang, C. (2019). Effects of aquatic phototrophs on seasonal hydrochemical, inorganic, and organic carbon variations in a typical karst basin, Southwest China. Environ. Sci. Pollut. Res. 26, 32836–32851. doi: 10.1007/s11356-019-06374-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, J., Luo, Q., Wang, D., and Wang, Z. (2015). Occurrences of pharmaceuticals in drinking water sources of major river watersheds. China. Ecotoxicol. Environ. Saf. 117, 132–140. doi: 10.1016/j.ecoenv.2015.03.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, J., Tang, X., Qin, Y., He, Q., Yi, Y., and Ji, Z. (2019). Karst rocky desertification progress: Soil calcium as a possible driving force. Sci. Total Environ. 649, 1250–1259. doi: 10.1016/j.scitotenv.2018.08.242

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanikawa, D., Nakamura, Y., Tokuzawa, H., Hirakata, Y., Hatamoto, M., and Yamaguchi, T. (2018). Effluent treatment in an aquaponics-based closed aquaculture system with single-stage nitrification–denitrification using a down-flow hanging sponge reactor. Int. Biodeterior. Biodegrad. 132, 268–273. doi: 10.1016/j.ibiod.2018.04.016

CrossRef Full Text | Google Scholar

Tate, R. L. (2020). Soil Microbiology, 1st New York: Wiley.

Google Scholar

Torres, N. I., Yu, X., Padilla, I. Y., Macchiavelli, R. E., Ghasemizadeh, R., Kaeli, D., et al. (2018). The influence of hydrogeological and anthropogenic variables on phthalate contamination in eogenetic karst groundwater systems. Environ. Pollut. 237, 298–307. doi: 10.1016/j.envpol.2018.01.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Webster, K. D., Drobniak, A., Etiope, G., Mastalerz, M., Sauer, P. E., and Schimmelmann, A. (2018). Subterranean karst environments as a global sink for atmospheric methane. Earth Planet. Sci. Lett. 485, 9–18. doi: 10.1016/j.epsl.2017.12.025

CrossRef Full Text | Google Scholar

Wu, X., Wu, L., Liu, Y., Zhang, P., Li, Q., Zhou, J., et al. (2018). Microbial Interactions With Dissolved Organic Matter Drive Carbon Dynamics and Community Succession. Front. Microbiol. 9:1234. doi: 10.3389/fmicb.2018.01234

PubMed Abstract | CrossRef Full Text | Google Scholar

Xi, X., Wang, M., Chen, Y., Yu, S., Hong, Y., Ma, J., et al. (2015). Adaption of the microbial community to continuous exposures of multiple residual antibiotics in sediments from a salt-water aquacultural farm. J. Hazard. Mater. 290, 96–105. doi: 10.1016/j.jhazmat.2015.02.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, S., Wang, X., Ma, W., Liu, X., Zhang, B., Huang, F., et al. (2020). Response of microbial communities of karst river water to antibiotics and microbial source tracking for antibiotics. Sci. Total Environ. 706:135730. doi: 10.1016/j.scitotenv.2019.135730

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, W., Sun, Y., Zhang, T., Ding, X., Li, Y., Wang, M., et al. (2015). Antibiotics, Antibiotic Resistance Genes, and Bacterial Community Composition in Fresh Water Aquaculture Environment in China. Microb. Ecol. 70, 425–432. doi: 10.1007/s00248-015-0583-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, Y., Tian, J., Quine, T. A., Powlson, D., Xing, K., Yang, L., et al. (2020). The persistence of bacterial diversity and ecosystem multifunctionality along a disturbance intensity gradient in karst soil. Sci. Total Environ. 748:142381. doi: 10.1016/j.scitotenv.2020.142381

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, L., Herrmann, M., Kampe, B., Lehmann, R., Totsche, K. U., and Küsel, K. (2020). Environmental selection shapes the formation of near-surface groundwater microbiomes. Water Res. 170:115341. doi: 10.1016/j.watres.2019.115341

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Song, X., Hallerman, E., and Huang, Z. (2020). Microbial community structure and nitrogen removal responses of an aerobic denitrification biofilm system exposed to tetracycline. Aquaculture 529:735665. doi: 10.1016/j.aquaculture.2020.735665

CrossRef Full Text | Google Scholar

Yuan, D. X. (1997). Rock desertification in the subtropical karst of South China. Ztschrift Geomorphol. 108:81–90.

Google Scholar

Yun, Y., Wang, H., Man, B., Xiang, X., Zhou, J., Qiu, X., et al. (2016). The Relationship between pH and Bacterial Communities in a Single Karst Ecosystem and Its Implication for Soil Acidification. Front. Microbiol, 7: 1955. doi: 10.3389/fmicb.2016.01955

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, S., Xu, Z., Wang, X., Liu, W., Qian, L., Chen, X., et al. (2019). Time series analysis of antibacterial usage and bacterial resistance in China: observations from a tertiary hospital from 2014 to 2018. Infect. Drug Resist. 12, 2683–2691. doi: 10.2147/IDR.S220183

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Ji, T., Jiang, Y., Zheng, C., Yang, H., and Liu, Q. (2022). Long-term effects of three compound probiotics on water quality, growth performances, microbiota distributions and resistance to Aeromonas veronii in crucian carp Carassius auratus gibelio. Fish Shellfish Immunol. 120, 233–241. doi: 10.1016/j.fsi.2021.11.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B., Li, Y., Xiang, S., Yan, Y., Yang, R., Lin, M., et al. (2020). Sediment Microbial Communities and Their Potential Role as Environmental Pollution Indicators in Xuande Atoll. South China Sea. Front. Microbiol. 11:1011. doi: 10.3389/fmicb.2020.01011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B., Qin, S., Guan, X., Jiang, K., Jiang, M., and Liu, F. (2021). Distribution of Antibiotic Resistance Genes in Karst River and Its Ecological Risk. Water Res. 203:117507. doi: 10.1016/j.watres.2021.117507

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Xie, J., Liu, M., Tian, Z., He, Z., van Nostrand, J. D., et al. (2013). Microbial community functional structure in response to antibiotics in pharmaceutical wastewater treatment systems. Water Res. 47, 6298–6308. doi: 10.1016/j.watres.2013.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Gu, X., Gao, S., Geng, J., and Wang, X. (2012). Adsorption of tetracycline (TC) onto montmorillonite: Cations and humic acid effects. Geoderma 183-184, 12–18. doi: 10.1016/j.geoderma.2012.03.004

CrossRef Full Text | Google Scholar

Zhernakova, A., Kurilshikov, A., Bonder, M. J., Tigchelaar, E., Schirmer, M., Vatanen, T., et al. (2016). Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 352, 565–569. doi: 10.1126/science.aad3369

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, Y., Yang, Q., Fu, G., Xu, Y., Cheng, Y., Chen, C., et al. (2018). Denitrifying microbial community with the ability to bromate reduction in a rotating biofilm-electrode reactor. J. Hazard. Mater. 342, 150–157. doi: 10.1016/j.jhazmat.2017.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., and Xing, J. (2021). Haloalkaliphilic denitrifiers-dependent sulfate-reducing bacteria thrive in nitrate-enriched environments. Water Res. 201:117354. doi: 10.1016/j.watres.2021.117354

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H., Jiang, C., and Liu, S. (2022). Microbial roles in cave biogeochemical cycling. Front. Microbiol. 13:950005. doi: 10.3389/fmicb.2022.950005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Wang, J., Hu, M., and Jia, L. (2019). Geographical detection of groundwater pollution vulnerability and hazard in karst areas of Guangxi Province. China Environ. Pollut. 245, 627–633. doi: 10.1016/j.envpol.2018.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar