J. Microbiol. Biotechnol. 2020. 30(9): 1355–1366 https://doi.org/10.4014/jmb.2004.04023 Review

Sediment Bacterial Community Structure under the Influence of Different Domestic Sewage Types

Lei Zhang1*, Mengli Xu1, Xingchen Li1, Wenxuan Lu2, and Jing Li2

1School of Civil Engineering and Architecture, Chuzhou University, Chuzhou 239000, P.R. China 2Fisheries Research Institute, Anhui Academy of Sciences, Hefei 230001, P.R. China

Sediment bacterial communities are critical to the biogeochemical cycle in river ecosystems, but our understanding of the relationship between sediment bacterial communities and their specific input streams in rivers remains insufficient. In this study, we analyzed the sediment bacterial community structure in a local river receiving discharge of urban domestic sewage by applying Illumina MiSeq high-throughput sequencing. The results showed that the bacterial communities of sediments samples of different pollution types had similar dominant phyla, mainly Proteobacteria, Actinobacteria, Chloroflexi and Firmicutes, but their relative abundances were different. Moreover, there were great differences at the genus level. For example, the genus Bacillus showed statistically significant differences in the hotel site. The clustering of bacterial communities at various sites and the dominant families (i.e., Nocardioidaceae, and Sphingomonadaceae) observed in the residential quarter differed from other sites. This result suggested that environmentally induced species sorting greatly influenced the sediment bacterial community composition. The bacterial co- occurrence patterns showed that the river had a nonrandom modular structure. Microbial taxonomy from the same module had strong ecological links (such as the nitrogenium cycle and degradation of organic pollutants). Additionally, PICRUSt metabolic inference analysis showed the most important function of river bacterial communities under the influence of different types of domestic sewage was metabolism (e.g., genes related to xenobiotic degradation predominated in residential quarter samples). In general, our results emphasize that the adaptive changes and interactions in the bacterial community structure of river sediment represent responses to different exogenous pollution sources. Keywords: Urban river, domestic sewage, bacterial community, metabolic function, co-occurrence network

Introduction Received: April 10, 2020 Rivers are the major source of renewable water for freshwater ecosystems, and they have important ecological Accepted: June 22, 2020 and economic value [1-3]. Currently, with the rapid growth of urbanization and population, much domestic sewage is discharged into rivers [4, 5]. Untreated or only partially treated domestic sewage contains a large number First published online: June 24, 2020 of pollutants, including antibiotics [6, 7], pathogenic microorganisms [4, 8], organic pollutants [9], suspended solids [10] and nutrients [11]. The chemical degradation of these pollutants may cause the eutrophication of the *Corresponding author river, decrease biological diversity, and accelerate the spread of pathogens in natural water systems [5, 12, 13], Phone: +86-550-3511822 thereby exposing water users to serious health risks. A more serious issue is the slow accumulation of pollutants in Fax: +550-3511822 the water sediment. Once released, the pollutants can be suspended in the water, which will again affect the quality E-mail: [email protected] of the water column and increase the risk of human infection [14]. Therefore, a sensitive and robust indicator is needed to identify and monitor changes within river sediments under the influence of domestic sewage. Supplementary data for this Fortunately, sediment microbial properties, especially bacterial diversity and community structure, represent paper are available on-line only very promising candidate bioindicators of polluted rivers [15]. Sedimentary bacteria are highly responsive to at http://jmb.or.kr. environmental changes [16], and they are the first to interact with dissolved matter [17]. Domestic sewage pISSN 1017-7825 containing organic and inorganic pollutants significantly changes the basic parameters of river sediment, such as eISSN 1738-8872 the nitrogen [18], phosphorus [19], and dissolved oxygen [4] content, thus further reshaping the bacterial community structure. Copyright© 2020 by There has long been keen interest to investigate not only the physical and chemical properties of rivers (such as The Korean Society for heavy metals [20], hydrocarbons [21] or nutrients [22] ), but also the composition of microbial communities in Microbiology and extreme or complex river ecosystems, their structure, and the relationship between diversity and the environment Biotechnology [23-26]. However, the contribution of these complex pollution sources to river ecosystems is comprehensive [1, 62]. Therefore, exploring a single source of pollution can provide an accurate understanding of the relationship

September 2020 ⎪ Vol. 30 ⎪ No. 9 1356 Zhang et al.

between bacterial communities in rivers and their specific input streams. Studies have shown that domestic wastewater is a critical issue in public health and disease prevention worldwide [28]. In recent years, research examining the impact of domestic pollution on river biodiversity has mainly focused on the impact on aquatic plants [29, 30] or plankton [31] and has poorly addressed bacterial communities in river sediments. Abiotic environmental factors affect microbial community function and the ecological mechanisms of interspecific interactions, excluding the composition and structure of microbial communities [32]. Research has shown that the function of the microbial community can provide a new perspective for the study of ecosystem function changes under different environmental conditions, including those related to the biogeochemical cycles of C, N, and P, as well as biomass production [33-35]. In addition, it has been proposed that microbial communities usually have nonrandom co-occurrence patterns and a modular structure [36], which strongly suggests that biological interactions play an important role in controlling community aggregation and ecosystem functions. We studied the composition of bacterial communities in river sediment and predicted the potential functional consequences under the influence of different domestic sewage types by uniting the bacteria-environment relationship and bacterial-bacterial interactions. In this study, we used the Qingliu River and its four typical domestic pollution input streams as an example. Our specific goals were (1) to explore the reasons for the difference of bacterial diversity composition in river sediments affected by different domestic sewage and (2) to investigate the correlation between the bacterial diversity of river sediments and their habitats under different domestic sewage pollution. The results from this study can provide a reference for the rational management and conservation of river ecosystems. Materials and Methods Description of the Study Sites and Sample Collection The sampling sites were located in the Qingliu River (watershed area of 1,265 km2 and mainstream length of 84.1 km), a typical urban river in Chuzhou City, China. This river is surrounded by residential and commercial areas, throughout the entire city of Chuzhou from the northwest to the southeast (Fig. 1), and it finally enters the Yangtze River through the Chu River. Excessive domestic sewage input results in the deterioration of this urban river. Sediment pollution under the long-term influence of domestic sewage has become an important issue in river ecosystem [37-39]. In this study, four representative areas affected by typical domestic pollution from the school (XX), residential quarter (JM), hospital (YY), and hotel (JD), were chosen as the sampling sites (Fig. 1). On the morning of May 6, 2019, three replicate sediment samples were collected at each sampling site, and these subsamples were obtained at 3 m intervals. The 12 sampling sites are shown in Fig. 1.The surface (< 5 cm in depth) sediment at each site was collected using a Peterson stainless steel grab sampler. Large fragments were removed, and the samples were stored in sterile polyethylene zipper bags, sealed, and transported to the laboratory on ice within 4 h. The samples

Fig. 1. Map showing the location and distribution of sampling sites across the urban region of Qingliu River (Chuzhou section).

J. Microbiol. Biotechnol. Bacterial Communities of River Sediments 1357 were then homogenized completely, and subsamples were placed in 2.5 ml sterile centrifuge tubes. The remaining samples were kept at 4°C for immediate physicochemical analysis.

Analysis of Physicochemical Properties The temperature (T), dissolved oxygen (DO), and pH of the overlying water on the samples were measured with a multiparameter sonde (YSI 6600V2, USA) during sample collection. The total nitrogen (TN), total phosphorus (TP), and total organic carbon (TOC) contents were determined in the laboratory according to standard methods [38].

DNA Extraction A 0.25 g sediment sample was collected, and then genomic DNA from each sediment sample was extracted in triplicate using a DNA Isolation Kit (E.Z.N.A., Omega, USA) in accordance with the manufacturer’s instructions. Subsequently, the concentration and purity of genomic DNA were detected with an ultraviolet spectrophotometer, and verified by 1% agarose gel electrophoresis. The measured DNA sample was stored at -20°C for 16S rRNA gene PCR amplification and sequence analyses [40].

PCR Amplification and Pyrosequencing The V3-V4 region of the bacterial 16S rRNA gene was amplified using primers 338F (5′-ACTCCTACGGGA GGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). According to formal experimental conditions, all the samples were processed. The PCR products were mixed and examined by 2% agarose gel electrophoresis. The PCR products were purified using AMPure beads and finally entrusted to Majorbio Bio- pharm Technology Co., Ltd., Shanghai, China for sequencing.

Sequencing Data Optimization and OTU Clustering In this study, Illumina MiSeq sequencing was used to obtain double-end sequence data, which were optimized with fastp software. The details were as follows: First, pairwise reads were combined into a sequence following the overlapping relationships between PE reads, and the effects of the merge were filtered by quality control. Then, the samples were distinguished based on the barcode and primer sequences at both ends of the sequence. Obtaining an effective sequence and correcting the sequence direction were performed to optimize the data. The optimized sequences were clustered into Operational Taxonomic Units (OTUs) using Usearch (version 7.1 http://drive5.com/uparse/) [41]. The RDP classifier (version 2.2 http://sourceforge.net/projects/rdp-classifier/) Bayesian algorithm was used to classify the 97% similarity level representative OTU sequences to obtain the species classification information corresponding to each OTU [42].

Statistical Analysis Based on the OTU information, the alpha diversity index of differences in community composition (Chao 1, ACE, Shannon, Simpson) and a sequencing depth index (Good's coverage) were calculated using Mothur software (version v.1.30.1 http://www.mothur.org/wiki/Schloss_SOP#Alpha_diversity) [43], and rarefaction curve graphs were generated with the R language tool (https://cran.R-project.org/web/packages/Vegan/) [44]. Beta diversity was represented by Bray-Curtis distance matrices generated from the OTUs table and statistical analysis and mapping using principal coordinate analysis (PCoA) in R language [45, 46]. According to the results of the taxonomic analysis, the community structure composition at different classification levels (such as the phylum, genus, and OTU, etc.) could be obtained. Based on the Kruskal-Wallis ranking and test, using the state package of R and the scipy package of Python to generate histograms, a significant difference test was performed between groups [47]. A similarity analysis (ANOSIM) was performed with the Bray-Curtis similarity matrix of the original pyrophosphate/phosphate data using the R-language vegan software package [48]. The relationship between the bacterial community structure and environmental variables was assessed using the Redundancy discriminant analysis (RDA) vegan package in R language [49]. The prediction of bacterial function and metabolic pathways was performed using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) [35, 50]. Based on the correlation between phylogeny and functions, PICRUSt was used to obtain information at different pathway levels (levels 1–3) and a table of the abundance of each level [88]. The co–occurrence networks of the microbial communities were determined via Spearman correlation. Gephi was used for network visualization and modularization analysis [51]. The network topology characteristics determined by Gephi include modularity, average path length, and the topology properties of each note include the degree and betweenness centrality [52]. Specific steps of the analysis were based on an online analysis platform (https://www.i-sanger.com/). The original sequencing data have been submitted to the NCBI Sequence Read Archive and the project accession code is SRP241215. Results Physicochemical Properties of the Sediment Samples The physicochemical parameters at the different sampling sites are presented in Table S1. Overall, the pH value of the overlying water ranged from neutral to alkaline conditions (pH = 7.46 to 8.72). The total phosphorus (TP) content of the river sediments did not significantly differ under the different types of domestic sewage pollution. However, there was a large difference in total organic carbon (TOC), ranging from 0.49 to 4.73 mg/g and the samples from the hotel site (JDS) has the lowest content. In addition, the total nitrogen (TN) content of JDS was

September 2020 ⎪ Vol. 30 ⎪ No. 9 1358 Zhang et al.

significantly less than the other samples. Notably, the contents of total phosphorus (TP), total nitrogen (TN), and total organic carbon (TOC) were higher in the school site sample (XXS) than the other samples.

Diversity and Composition of the Bacterial Communities In our study, the diversity and phylogenetic structure of the bacterial communities in the sediment samples were analyzed using Illumina MiSeq high-throughput sequencing. A total of 452,025 effective 16S rRNA gene sequences were acquired after removing low-quality sequences and chimeras. Of these, after the clustering and alignment of 12 sediment sludge samples, 1619 OTUs were obtained according to the standard of 0.97 similarities. The average coverage of all samples was 99.01 ± 0.1% (mean ± s.d.) (Table S2), indicating that the sequencing depth was sufficient for community analysis. All rarefaction curves verged toward saturation (Fig. 2A), indicating a high degree of confidence in the bacterial community structure. As seen in Table S2, the total number of OTUs in the four samples differed. The average total numbers of OTUs in the hospital samples (YYS) and hotel samples (JDS) were 1094 and 1056, respectively. The samples from the school site (XXS) and residential quarter site (JMS) had 920 and 838 OTUs, respectively. The results showed that YYS had the highest species richness and that JMS had the lowest species richness. The change of Chao1 and ACE index was consistent with the total number of OTUs (Fig. 2C and Table S2), which confirmed the foregone conclusion. The Shannon and Simpson indices were used to estimate the microbial diversity in the samples [1]. The Shannon index values for the YY and JD samples were higher than those for the XXS and JMS samples (YYS and JDS were 6.71 and 6.67, respectively, and XXS and JMS were 6.30 and 5.9, respectively) (Table S2). The Simpson index for JMS was the largest (0.0123), and that for YYS was the smallest (0.0039); the values for XXS and JDS were 0.0078 and 0.0061, respectively (Table S2). Consequently, the YYS and JDS had a high, and XXS and JMS, a low species community diversity. The PCoA was used to observe the composition of the various bacterial communities (Fig. 3). PC1 and PC2 accounted for 20.79% and 45.92% of the variation, respectively. The XX, YY, and JD samples were gathered on the right side of the figure, and JMS was distributed in the middle. Generally, the closer two samples were plotted, the more similar was their composition. It is worth noting that JMS was far from other samples in the graph, indicating a large difference in bacterial community composition from the other samples. In addition, ANOSIM indicated that the bacterial community structure of river sediment significantly differed among the different types of domestic sewage pollution (ANOSIM R = 1, p = 0.001). The sequences of the samples were sorted at the phylum classification level, as shown in Fig. 4A. Thirteen different phyla and unclassified bacterial members were detected in the urban river sediment samples. Proteobacteria, Actinobacteria, Chloroflexi, Firmicutes, Acidobacteria, and Cyanobacteria (relative abundance > 1% at least in a sample) were the dominant bacteria in all samples, accounting for 83.98% to 93.89% of the total sequences. Proteobacteria was the most abundant phylum in the XXS, YYS and JMS, with relative proportions of 48.5%, 28.74%, and 33.89%, respectively. This dominant bacterial taxon was the second highest in JMS, accounting for 28.74% of the total. Actinobacteria had the highest content in JMS, accounting for 36.18% of the total community. In addition, compared with the other samples, Firmicutes had the lowest content in JMS, accounting for 1.67% of the total community. It is worth noting that Proteobacteria, Actinobacteria,

Fig. 2. (A) Rarefaction curves of OTUs clustered at 97% sequence identity across the five samples. Properties of the bacterial communities of the twelve river sediments. (B) Diversity, expressed as the Shannon Index: Higher values indicate higher diversity. (C) Richness, indicated by Chao's richness estimator Chao 1: higher values indicate higher diversity.

J. Microbiol. Biotechnol. Bacterial Communities of River Sediments 1359

Fig. 3. PCoA analysis of sedimentary bacterial community diversity.

Firmicutes, Bacteroidetes, Gemmatimonadetes, and Nitrospirae were identified in the four sample types (p < 0.05) (Fig. 4C). To further reveal the structure of the bacterial communities in the river sediments, the relative abundance and classification of OTUs were analyzed at the family and genus levels (Figs. S1 and 4B). When comparing the species composition of a particular bacterial family at all sampling locations at the taxonomic level, the difference in JMS was more pronounced. Nocardioidaceae, Sphingomonadaceae, Intrasporangiaceae, and Xanthomonadaceae were the dominant families (Fig. S1). At the genus level, Dechloromonas, Nitrospira, Mycobacterium, Ilumatobacter, Bacillus and Pseudarthrobacter were the richest bacterial taxa in the sediment samples (Fig. 4B). These bacteria were also the common genera in XXS, YYS, JMS, and JDS. Notably, Nocardia (3.83%), Pseudarthrobacter (2.97%), Marmoricola (7.26%), and Nocardioides (4.05%) were unique genera in JMS. Dechloromonas was a unique species in XXS, with a relative proportion of 6.87%. In addition, Sphingomonas was the main genus in JMS, with a relative abundance of 7.59%. The genus Agromyces was found in XXS (0.88%), YYS (0.32%), JMS (0.48%), and JDS (0.28%). We observed that 10 of the 15 genera showed statistically significant differences among the four groups (Fig. 5A). We found statistically significant differences in bacillus among the four groups, but not between the two

Fig. 4. Community composition of bacteria in the phylum level (A) and the genus level (B) in sediment samples. (C) Comparison of bacterial abundances in river sediments under the influence of domestic sewage at the phylum level.

September 2020 ⎪ Vol. 30 ⎪ No. 9 1360 Zhang et al.

Fig. 5. Comparison of bacterial abundances in river sediments under the influence of domestic sewage at the genus level (A). One-way ANOVA bar plot on genus level for Sphingomonas (B), Dechloromonas (C) and Marmoricola (D), and Bacillus (E) * stands for 0.01 b p ≤ 0.05, ** stands for 0.001 b p ≤ 0.01 and *** stands for p ≤ 0.001.

groups (Fig. 5E). Marmoricola (Fig. 5D), Dechloromonas (Fig. 5C), and Sphingomonas (Fig. 5B) showed statistically significant differences among the four sample types and the two groups.

Co–Occurrence Networks Considering the nonrandom community aggregation model of bacteria under different levels of domestic sewage pollution, a network interface was built to exhibit the topological and taxonomic characteristics of the bacterial co–occurrence patterns in river sediment (Fig. 6). In accordance with the Spearman correlation analysis, the sediment bacterial network in the river was composed of 265 nodes (i.e., OTUs) and 446 edges, and there was a strong correlation between nodes and nodes in pairs (Spearman’s ρ ≥ 0.7, p < 0.05) (Fig. 6). In addition, the analysis showed that the topological characteristics of the bacterial communities in river sediments could be distinguished among different levels of domestic pollution. The networks exhibited higher values of modularity (MD, 0.635; Values> 0.4 indicate that the network is modular) and the Clustering coefficient (CC) than the global network of all sampled sites, suggesting that the different pollution sites had obvious co–occurrence patterns and modular structures. Modularity analysis showed that the whole network could be divided into six main modules (Fig. 6A). In these modules, species interactions were more frequent and intense than in the rest of the community. The OTUs from modules II and VI had the highest relative abundances in YYS, while those from modules I and IV tended to occur in XXS. The OTUs of modules II and VI were relatively high in abundance in YYS, while those from modules I and

J. Microbiol. Biotechnol. Bacterial Communities of River Sediments 1361

Fig. 6. Spearman’s correlation-based network analysis between OTUs for the inflow rivers samples. The size of each node is proportional to the number of connections (i.e., degree), and the nodes are colored according to different types of modularity classes (A) and phylum (B), respectively. Colored nodes (except light gray) represent the six major groups and light gray represents all modules except the six main major modules. Orange and blue edges indicate positive and negative correlations.

IV tended to be high in abundance in XXS. However, the classification relationship was obviously a critical factor in determining the module structure in the network. The nodes in modules I and IV were mostly Proteobacteria. The nodes in module II belonged to Proteobacteria, Firmicutes and Chloroflexi. The nodes in module III mostly belonged to Actinobacteria and Proteobacteria. The nodes in module V were mainly Cyanobacteria. Unclassified Anaerolineaceae dominated module VI (16.67%), while unclassified Cyanobacteria dominated module V (15.79%). According to the betweenness centrality scores, Variibacter, Parasegetibacter, unclassified Caldilineaceae, unclassified Nitrosomonadaceae and unclassified Geodermatophilaceae were identified as the top five genera. This finding suggested that these bacteria played a key role as keystone taxa within the co-occurrence network.

Environmental Factor Analysis Environmental factors are an important driver of the structure of the microbial community [53]. RDA was performed at the phylum level to identify the relationships between the composition of the bacterial community and environmental parameters, as shown in Fig. 7. RDA1 and RDA2 accounted for 74.91% and 21.74% of the interpretations, respectively. Notably, TP (p = 0.001) and pH (p < 0.005) were the environmental variables affecting the environmental relationships of bacterial communities, and TN (p > 0.05) and TOC (p > 0.05) had weaker effects. Furthermore, Actinobacteria, Bacteroidetes and Acidobacteria were positively correlated to pH, TN, and TOC and negatively correlated with TP. Interestingly, Proteobacteria and Cyanobacteria were positively related to all environmental factors. The most abundant phylum, Actinobacteria, was negatively related to pH, TOC, and TP but had little relationship with TN.

Fig. 7. RDA analysis of bacterial community species (phylum) in sediments and environmental factors.

September 2020 ⎪ Vol. 30 ⎪ No. 9 1362 Zhang et al.

Fig. 8. The comparison of bacterial community functions predicted by PICRUSt at level 3 related to the degradation and metabolism of xenobiotics between different groups with the influence of domestic sewage.

PICRUSt Functional Predictive Analysis To understand the function of river sediment bacteria under domestic sewage pollution, this study used PICRUSt software for prediction by comparing the relative abundance of the four groups in functional categories (level 1, Fig. S2), subcategories (level 2, Fig. S3), and individual pathways (level 3, Fig. 8). The results showed that the major functional gene families of all samples were involved in Metabolism, Environmental Information Processing, Genetic Information Processing, Cellular Processes, Human Diseases and Organismal Systems. It is worth noting that there was a high abundance of genes related to metabolism in bacterial communities of XXS, YYS, JMS and JDS (Fig. S2). Therefore, to further determine the functional differences, we compared the relative abundance of genes related to metabolism at level -2. As shown in Fig. S3, the genes involved in amino acid and carbohydrate metabolism were present in all samples with the highest abundance. In addition, the samples from residential quarter (JMS) had a relatively high abundance of genes related to xenobiotics biodegradation and metabolism, lipid metabolism and terpenoid and polyketide metabolism (Fig. S3). In general, rivers are the main recipients of pollutants and xenobiotics, wherein most aquatic organisms and bacteria are exposed to these xenobiotics [1, 54]. Among the metabolic genes and pathways associated with xenobiotics biodegradation at level- 3, the genes related to benzoate and aminobenzoate biodegradation were the most abundant genes in JMS and XXS (Fig. 8). Samples with an abundance of predictive genes related to Bisphenol degradation were found in JM (Fig. 8). Discussion Freshwater ecosystems are particularly vulnerable to the surrounding environment [1, 27, 55]. Fluctuations caused by various pollutants accelerate the energy flow and material circulation in the river ecosystem, which may have a direct impact on the microbial community in the river. Our research showed significant differences in the bacterial diversity of river sediments under different types of domestic sewage pollution (Table S2). A relatively low bacterial diversity was detected in samples XXS and JMS, potentially due to the massive enrichment of some species, which are very suitable for conditions of sewage pollution in schools and residential areas. In addition, JMS showed the lowest species diversity (Table S2). The main reason may be related substances (such as household chemicals, fungicide and pharmaceutical residues) in the complex living environment of the residential area, which inhibits the activity of certain bacteria in the river [87, 88]. PCoA (Principal co-ordinates analysis) showed that the bacterial community structure of YYS (the samples from the hospital) and JDS (the samples from the hotel) had more similarities (Figs. 3 and 4A). This phenomenon may mean that many taxa are generalists, able to utilize multiple nutrients [17], and/or able to tolerate severe environmental disturbances [59]. Based on the above results, the bacterial community is driven by various environmental factors, which indicates that the bacterial communities in river sediments may have developed phylogenetic diversity to cope with pollutants from different sources. The bacterial communities in our study included aquatic microorganisms and microorganisms of human origin, and the bacterial compounds from different presumptive sources also differed among sampling locations.

J. Microbiol. Biotechnol. Bacterial Communities of River Sediments 1363

By analyzing sediments under the influence of domestic sewage in this study, Proteobacteria, Actinobacteria, Chloroflexi, Firmicutes, Bacteroidetes, and Acidobacteria were found to be the most abundant phyla (percent of sequences). Consistent with other studies, these members are usually detected in freshwater sediment environments [60-62]. Although most bacterial groups are present in sediments, a comprehensive comparison revealed differences among the dominant bacteria in the samples (Fig. 4A). Historically, Actinobacteria were considered to be the primary active resident taxon in soils [63]. However, according to the data in this study, the main phylum of gram-negative bacteria, Proteobacteria, was the most abundant phylum involved in various biogeochemical processes in aquatic ecosystems (e.g., carbon and nitrogen cycles) [38]. In addition, RDA showed that the Proteobacteria in the river sediment was directly proportional to TN and TP (Fig. 7), which may have been caused by the high nutrient content in the sewage. This observation indicated that a high abundance of Proteobacteria is related to high nutrient utilization [56, 57]. Among all the classification results, the main sedimentary bacterial populations in the samples from the residential quarter (JMS) were the most significantly different (Fig. 4A), mainly including taxa related to nutrient contamination and other anthropogenic disturbances. The families Nocardioidaceae, Sphingomonadaceae, and Intrasporangiaceae were dominant at the JM site. Sphingomonadaceae are often found in habitats polluted by recalcitrant aromatic compounds from natural or anthropogenic sources [64] (Fig. S1). Specifically, these populations have the ability to degrade recalcitrant high-molecular-weight compounds [65]. Members of Chitinophagaceae can degrade chitin and cellulose [66]. Nocardioides species are chemoorganotrophs and utilize a variety of carbon and nitrogen sources [67]. Therefore, we conclude that the relative advantages of these high-abundance organisms reflect that the carbon cycle plays an important role at the JM site. In addition, nutrient conditions may select certain bacterial species, resulting in differences in the distribution of bacterial populations. The genus Dechloromonas was the most dominant group at the XX and YY sites. This genus is commonly found in freshwater sediment environments and can oxidize aromatic compounds such as toluene and chlorobenzoate [48]. Hospitals represent an irrefutable source of the release of numerous compounds into the aquatic environment due to laboratory activities or the release of pharmaceuticals into wastewater [70]. Therefore, the relative dominance of the genus Dechloromonas at the XX and YY sites may reflect the presence of high levels of organics. Acidobacteria and Bacteroidetes also showed good properties in the samples from the school and hotel (XXS and YYS, respectively) (Fig. 4A). Acidobacteria are renowned for their resistance to contaminants such as petroleum compounds [71, 72]. Bacteroides are gram-negative heterotrophic bacteria that are commonly found in freshwater ecosystems. They are ubiquitous in organic-rich systems and are considered to be adept at the degradation of high-molecular-weight organics in different ecological environments [73-75]. The RDA showed that Acidobacteria were positively correlated with YYS, Bacteroidetes were positively correlated with XXS, and the XXS and YYS samples were rich in TOC (Table S2). Therefore, the presence of Acidobacteria and Bacteroidetes further confirmed the prevalence of organic compound waste in areas where YY and XX deposits were located. These observations are consistent with other studies, showing changes in bacterial communities caused by different types of wastewater discharge [17, 38]. Network analysis can help explain the basic network of microbe-microbe interactions in addition to the changing bacterial community diversity and composition. The analysis of habitat networks in this study provides a visualization of the mutuality of bacterial genera and valuable supplementary information for the study of river sediment communities under the influence of domestic sewage (Fig. 6A). According to the betweenness centrality scores, the genera Variibacter, Parasegetibacter, unclassified Caldilineaceae, unclassified Nitrosomonadaceae and unclassified Geodermatophilaceae were identified as the top five, indicating that they may play a key role in maintaining the structure and function of sediment ecological communities in rivers. The genus Variibacter is a strictly aerobic bacterial taxon that provides essential information for understanding bacterial activities associated with the growth of plants [76]. The genera Parasegetibacter and unclassified Caldilineaceae have potential nitrification and phosphorus accumulation abilities [39, 77]. Unclassified Nitrosomonadaceae are closely related to N cycling. Unclassified Geodermatophilaceae members have properties that inhibit the production of reactive oxygen species and play a key role in extreme environments where nutrient supplies are scarce [78]. By integrating the classification allocation of OTU nodes and network structure, we found that closely related categories tended to interconnect and gather together (Fig. 6B), which can be explained by the overlapping strong niches in closely associated species or synergistic relations. In our study, the nodes were divided into four main modules, and different modules responded to different functions [79]. A strong ecological connection of co- occurrence cluster correlations was observed. In module IV, certain bacteria were related to the degradation of organic carbonate. For example, unclassified Anaerolineaceae, which are dominant in YYS, play a key role in the initial activation of long-chain alkane biodegradation [80], and can anaerobically decompose carbohydrates via fermentation [81]. Therefore, the large presence of unclassified Anaerolineaceae in module IV indicated that the organic matter in YY controlled these bacterial activities. In module II, the genus Bacillus was relatively abundant in JDS (Fig. 4B). Urban sewage, especially domestic sewage, has been confirmed to contain a large number of pathogenic microorganisms, including pathogenic bacteria and enteroviruses [82, 83]. In this study, the genus Bacillus was abundant in JDS. These bacteria are potentially pathogenic and can be found in fresh and aging feces [84]. There is no doubt that the organic and even toxic substances in JD wastewater have a significant impact on the role of the genus Bacillus in the modular structure. Conversely, the major taxa in module II likely participate in biogeochemical N cycles. Nitrospira species are the dominant nitrite oxidizers within freshwater sediments [85]. The genus Paenisporosarcina is associated with the function of nitrate reduction [10], and Comamonadaceae are known denitrifiers [86]. Therefore, the complexity of community structure, as well as the co-occurrence patterns

September 2020 ⎪ Vol. 30 ⎪ No. 9 1364 Zhang et al.

of bacterial communities in river ecosystems affected by domestic sewage, seems to be nonrandom and functionally driven. Based on PICRUSt prediction, this study revealed the potential functional characteristics of bacterial communities in river sediments caused by domestic pollution. A relationship between xenobiotic degradation genes and xenobiotic biodegradation rates has been suggested. Additionally, the biodegradable genes could be used as indicators for the existence of xenobiotics and their metabolites [89, 90]. In this study, genes related to benzoate and aminobenzoate biodegradation were the most abundant xenobiotic degradation genes in all the samples. As the key mediator in PAH biodegradation, there are a large number of genes related to benzoate and aminobenzoate biodegradation, indicating that the areas polluted by domestic sewage contain persistent organic pollutants with high contents from both JMS and XXS. This functional response of the bacterial community may accelerate bioremediation in polluted areas. Moreover, genes related to bisphenol degradation were detected in all river sediment samples (Fig. 8), and especially those from JMS with the highest value (8.19%), indicating that organic toxic chemicals may be contained in untreated wastewater from JMS. Bisphenol has been shown to exert acute toxicity toward aquatic organisms and to disrupt their endocrine function [91, 92]. Therefore, direct discharge of domestic sewage into rivers may pose a high risk to the bacterial ecology. In conclusion, this study revealed the structure, diversity, co-occurrence patterns, and metabolic functions of bacteria in the middle layer of urban rivers under the influence of different types of domestic sewage. Research has shown that different types of sewage in cities affect the bacterial community composition in rivers, and the difference in bacterial community composition is particularly significant in JMS areas. Moreover, we also found evidence of the same high abundance of bacterial communities in different locations (such as the high abundance of the genus Bacillus in JD and YY), as nutrient factors may affect the bacterial community assemblage within the river sediment. The occurrence of nonrandom co-occurrence and ecological function-driven modular patterns provides a new perspective of microbial assembly in terms of the effects of pollution on the sediment bacterial community. Our results indicate that TP is a more important factor affecting the distribution of bacteria in river sediments under the influence of domestic sewage. This study provides insights into our understanding and monitoring of pollution in the aquatic environment and contributes to the management of river environments. Acknowledgments This work was supported by the Key Program in the Youth Elite Support Plan in Universities of Anhui Province (gxyqZD2020046), the Key Research and development program of Anhui Province (202004i07020010) and the Anhui Academy of Agricultural Sciences Research and Application Innovation Team Project of Aquaculture Water Environment (2020YL043). Conflict of Interest The authors have no financial conflicts of interest to declare. References 1. Ibekwe AM, Ma J, Murinda SE. 2016. Bacterial community composition and structure in an Urban River impacted by different pollutant sources. Sci. Total Environ. 566: 1176-1185. 2. Araya R, Tani K, Takagi T, Yamaguchi N, Nasu M. 2003. Bacterial activity and community composition in stream water and biofilm from an urban river determined by fluorescent in situ hybridization and DGGE analysis. FEMS Microbiol. Ecol. 43: 111-119. 3. Koubar M, Houhou J, Hoteit M, Al-Iskandarani M. 2016. A comparative study of the physicochemical and bacterial characteristics of the three drinking water sources distributed in Haret-Hreik city, Lebanon. 4. Gao Y, Wang C, Zhang W, Di P, Yi N, Chen C. 2017. Vertical and horizontal assemblage patterns of bacterial communities in a eutrophic river receiving domestic wastewater in southeast China. Environ. Pollut. 230: 469-478. 5. Sheng Y, Qu Y, Ding C, Sun Q, Mortimer RJ. 2013. A combined application of different engineering and biological techniques to remediate a heavily polluted river. Ecol. Eng. 57: 1-7. 6. 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. 7. Chen J, Ying G-G, Wei X-D, Liu Y-S, Liu S-S, Hu L-X, et al. 2016. Removal of antibiotics and antibiotic resistance genes from domestic sewage by constructed wetlands: effect of flow configuration and plant species. Sci. Total Environ. 571: 974-982. 8. Kim M, Guerra P, Theocharides M, Barclay K, Smyth S, Alaee M. 2013. Parameters affecting the occurrence and removal of polybrominated diphenyl ethers in twenty Canadian wastewater treatment plants. Water Res. 47: 2213-2221. 9. Wang J-Z, Chen T-H, Zhu C-Z, Peng S-C. 2014. Trace organic pollutants in sediments from Huaihe River, China: Evaluation of sources and ecological risk. J. Hydrol. 512: 463-469. 10. Wang Z, Yang J, Zhang D, Zhou J, Zhang C, Su X, et al. 2014. Composition and structure of microbial communities associated with different domestic sewage outfalls. Genet. Mol. Res. 13: 7542-7552. 11. Nsenga Kumwimba M, Dzakpasu M, Zhu B, Wang T, Ilunga L, Kavidia Muyembe D. 2017. Nutrient removal in a trapezoidal vegetated drainage ditch used to treat primary domestic sewage in a small catchment of the upper Yangtze River. Water Environ. J. 31: 72-79. 12. He D, Zhang K, Tang J, Cui X, Sun Y. 2018. Using fecal sterols to assess dynamics of sewage input in sediments along a human- impacted river-estuary system in eastern China. Sci. Total Environ. 636: 787-797. 13. Liang Z, Siegert M, Fang W, Sun Y, Jiang F, Lu H, et al. 2018. Blackening and odorization of urban rivers: a bio-geochemical process. FEMS Microbiol. Ecol. 94: 10.1093/femsec/fix180. 14. Wildi W, Dominik J, Loizeau JL, Thomas RL, Favarger PY, Haller L, et al. 2004. River, reservoir and lake sediment contamination by heavy metals downstream from urban areas of Switzerland. Lakes Reservoir. Res. Manage. 9: 75-87. 15. Hu Y, Bai C, Cai J, Dai J, Shao K, Tang X, et al. 2018. Co-occurrence Network reveals the higher fragmentation of the bacterial community in Kaidu River than its tributaries in Northwestern China. Microbes Environ. 33: 127-134. 16. Abia ALK, Alisoltani A, Keshri J, Ubomba-Jaswa E. 2018. Metagenomic analysis of the bacterial communities and their functional profiles in water and sediments of the Apies River, South Africa, as a function of land use. Sci. Total Environ. 616: 326-334.

J. Microbiol. Biotechnol. Bacterial Communities of River Sediments 1365

17. Roberto AA, Van Gray JB, Leff LG. 2018. Sediment bacteria in an urban stream: spatiotemporal patterns in community composition. Water Res. 134: 353-369. 18. Chen ZF, Ying GG, Liu YS, Zhang QQ, Zhao JL, Liu SS, et al. 2014. Triclosan as a surrogate for household biocides: an investigation into biocides in aquatic environments of a highly urbanized region. Water Res. 58: 269-279. 19. Yadav R, Goyal B, Sharma R, Dubey S, Minhas P. 2002. Post-irrigation impact of domestic sewage effluent on composition of soils, crops and ground water—A case study. Environ. Int. 28: 481-486. 20. Virsek MK, Hubad B, Lapanje A. 2013. Mercury induced community tolerance in microbial biofilms is related to pollution gradients in a long-term polluted river. Aquat. Toxicol. 144: 208-217. 21. Scholes L, Faulkner H, Tapsell S, Downward S. 2008. Urban rivers as pollutant sinks and sources: a public health concern for recreational river users? Water Air Soil Poll. 8: 543-553. 22. Yang X, Huang S, Wu Q, Zhang R, Liu G. 2013. Diversity and vertical distributions of sediment bacteria in an urban river contaminated by nutrients and heavy metals. Front. Environ. Sci. Eng. 7: 851-859. 23. Moitra M, Leff LG. 2015. Bacterial community composition and function along a river to reservoir transition. Hydrobiologia 747: 201-215. 24. Cébron A, Coci M, Garnier J, Laanbroek HJ. 2004. Denaturing gradient gel electrophoretic analysis of ammonia-oxidizing bacterial community structure in the lower Seine River: impact of Paris wastewater effluents. Appl. Environ. Microbiol. 70: 6726-6737. 25. Köchling T, Sanz JL, Galdino L, Florencio L, Kato MT. 2017. Impact of pollution on the microbial diversity of a tropical river in an urbanized region of northeastern Brazil. Int. Microbiol. 20: 11-24. 26. Staley C, Gould TJ, Wang P, Phillips J, Cotner JB, Sadowsky MJ. 2014. Bacterial community structure is indicative of chemical inputs in the Upper Mississippi River. Front. Microbiol. 5: 524. 27. Feng B, Li X, Wang J, et al. 2009. Bacterial diversity of water and sediment in the Changjiang estuary and coastal area of the East China Sea. FEMS Microbiol. Ecol. 70: 236-248. 28. Mara, D.D. 2004. Domestic wastewater treatment in developing countries. Earthscan Publications, London. 29. Chati R. 2004. Effect of domestic sewage on fresh water bivalve molluses from manjara river. 30. Avery EL. 1970. Effects of domestic sewage on aquatic insects and salmonids of the East Gallatin River, Montana. Water Res. 4: 165- 177. 31. Kushwaha VB, Agrahari M. 2014. Effect of domestic sewage on zooplankton community in River Rapti at Gorakhpur, India. World J. Zool. 9: 86-92. 32. Wu H, Li Y, Zhang W, Wang C, Wang P, Niu L, et al. 2019. Bacterial community composition and function shift with the aggravation of water quality in a heavily polluted river. J. Environ. Manage 237: 433-441. 33. Kuang J, Huang L, He Z, Chen L, Hua Z, Jia P, et al. 2016. Predicting taxonomic and functional structure of microbial communities in acid mine drainage. ISME J. 10: 1527. 34. Su Z, Dai T, Tang Y, Tao Y, Huang B, Mu Q, et al. 2018. Sediment bacterial community structures and their predicted functions implied the impacts from natural processes and anthropogenic activities in coastal area. Mar. Pollut. Bull. 131: 481-495. 35. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. 2013. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31: 814-821. 36. Barberán A, Bates ST, Casamayor EO, Fierer N. 2014. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 8: 952. 37. Fernandes H, Jungles MK, Hoffmann H, Antonio RV, Costa RH. 2013. Full-scale sequencing batch reactor (SBR) for domestic wastewater: performance and diversity of microbial communities. Bioresour. Technol. 132: 262-268. 38. Wang L, Zhang J, Li H, Yang H, Peng C, Peng Z, et al. 2018. Shift in the microbial community composition of surface water and sediment along an urban river. Sci. Total Environ. 627: 600-612. 39. Kim SJ, Ahn JH, Weon HY, Hong SB, Seok SJ, Kwon SW. 2015. Parasegetibacter terrae sp. nov., isolated from paddy soil and emended description of the genus Parasegetibacter. Int. J. Syst. Evol. Microbiol. 65: 113-116. 40. Keshri J, Ram AP, Sime-Ngando T. 2018. Distinctive patterns in the taxonomical resolution of bacterioplankton in the sediment and pore waters of contrasted freshwater lakes. Microb. Ecol. 75: 662-673. 41. Edgar RC. 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26: 2460-2461. 42. Zhang L, Shen Z, Fang W, Gao G. 2019. Composition of bacterial communities in municipal wastewater treatment plant. Sci. Total Environ. 689: 1181-1191. 43. Shao J, He Y, Zhang H, Chen A, Lei M, Chen J, et al. 2016. Silica fertilization and nano-MnO 2 amendment on bacterial community composition in high arsenic paddy soils. Appl. Microbiol. Biotechnol. 100: 2429-2437. 44. Team, R. D. C. 2009. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.Available online: http://www. R-project. org (accessed on Feb. 10, 2015). 45. Bray JR, Curtis JT. 1957. An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27: 325-349. 46. Ibarbalz FM, Figuerola EL, Erijman L. 2013. Industrial activated sludge exhibit unique bacterial community composition at high taxonomic ranks. Water Res. 47: 3854-3864. 47. Fan L, Wang Z, Chen M, Qu Y, Li J, Zhou A, et al. 2019. Microbiota comparison of Pacific white shrimp intestine and sediment at freshwater and marine cultured environment. Sci. Total Environ. 657: 1194-1204. 48. Keshri J, Ram ASP, Colombet J, Perriere F, Thouvenot A, Sime-Ngando T. 2017. Differential impact of lytic viruses on the taxonomical resolution of freshwater bacterioplankton community structure. Water Res. 124: 129-138. 49. Liu J, Liu X, Wang M, Qiao Y, Zheng Y, Zhang X-H. 2015. Bacterial and archaeal communities in sediments of the north Chinese marginal seas. Microb. Ecol. 70: 105-117. 50. Dong Z, Hong M, Hu H, Wang Y, Zhang D, Wang K. 2018. Effect of excess nitrogen loading on the metabolic potential of the bacteria community in oligotrophic coastal water. Acta Sci. Circumst. 38: 457-466. 51. Bordenave S, Goñi-Urriza MS, Caumette P, Duran R. 2007. Effects of heavy fuel oil on the bacterial community structure of a pristine microbial mat. Appl. Environ. Microbiol. 73: 6089-6097. 52. Li Y, Wu H, Shen Y, Wang C, Wang P, Zhang W, et al. 2019. Statistical determination of crucial taxa indicative of pollution gradients in sediments of Lake Taihu, China. Environ. Pollut. 246: 753-762. 53. Haller L, Tonolla M, Zopfi J, Peduzzi R, Wildi W, Pote J. 2011. Composition of bacterial and archaeal communities in freshwater sediments with different contamination levels (Lake Geneva, Switzerland). Water Res. 45: 1213-1228. 54. Mark Ibekwe A, Leddy MB, Bold RM, Graves AK. 2012. Bacterial community composition in low-flowing river water with different sources of pollutants. FEMS Microbiol. Ecol. 79: 155-166. 55. Schmeller DS, Loyau A, Bao K, Brack W, Chatzinotas A, Vleeschouwer, FD, et al. 2018. People, pollution and pathogens – Global change impacts in mountain freshwater ecosystems[J]. Sci. Total Environ. 622: 756-763. 56. MacDonald DD, Ingersoll CG. 2002. A guidance manual to support the assessment of contaminated sediments in freshwater ecosystems. Volume I: An ecosystem-based framework for assessing and managing contaminated sediments. Chicago (IL): U.S. Environmental Protection Agency Great Lakes National Program Office. 57. Withers P, Jarvie H. 2008. Delivery and cycling of phosphorus in rivers: a review. Sci. Total Environ. 400: 379-395.

September 2020 ⎪ Vol. 30 ⎪ No. 9 1366 Zhang et al.

58. Rui L, Xiao D. 2017. Influences of sediment characteristics and overlying water quality on sediment bacterial communities in a seasonal sandy river. China Environ.l Sci. 37: 4342-4352. 59. Wittebolle L, Marzorati M, Clement L, Balloi A, Daffonchio D, Heylen K, et al. 2009. Initial community evenness favours functionality under selective stress. Nature 458: 623-626. 60. Li D, Yang M, Li Z, Qi R, He J, Liu H. 2008. Change of bacterial communities in sediments along Songhua River in Northeastern China after a nitrobenzene pollution event. FEMS Microbiol. Ecol. 65: 494-503. 61. Tamaki H, Sekiguchi Y, Hanada S, Nakamura K, Nomura N, Matsumura M, et al. 2005. Comparative analysis of bacterial diversity in freshwater sediment of a shallow eutrophic lake by molecular and improved cultivation-based techniques. Appl. Environ. Microbiol. 71: 2162-2169. 62. Guo X, Lu D, Niu Z, et al. 2018. Bacterial community structure in response to environmental impacts in the intertidal sediments along the Yangtze Estuary, China. Mar Pollut Bull, 126: 141-149. 63. Sun W, Xia C, Xu M, Guo J, Sun G. 2017. Seasonality affects the diversity and composition of bacterioplankton communities in dongjiang river, a drinking water source of Hong Kong. Front. Microbiol. 8: 1644. 64. Glaeser SP, Kämpfer P. 2014. The family sphingomonadaceae. In: Rosenberg E, DeLong E, Lory S, Stackebrandt E, Thompson F (Eds.), The Prokaryotes - Alphaproteobacteria and Betaproteobacteria. Springer, Berlin Heidelberg, pp. 641-707. 65. Glaeser SP, Grossart HP, Glaeser J. 2010. Singlet oxygen, a neglected but important environmental factor: short‐term and long‐term effects on bacterioplankton composition in a humic lake. Environ. Microbiol. 12: 3124-3136. 66. Rosenberg E. 2014. The Family Chitinophagaceae, pp. 493-495. Ed. Springer Berlin Heidelberg. 67. Tóth E, Borsodi A. 2014. The Family Nocardioidaceae, pp. 651-694. The Prokaryotes: Actinobacteria, Ed. Springer-Verlag Berlin Heidelberg 68. Alvarado P, Huang Y, Wang J, Garrido I, Leiva S. 2018. Phylogeny and bioactivity of epiphytic Gram-positive bacteria isolated from three co-occurring antarctic macroalgae. Antonie Van Leeuwenhoek 111: 1543-1555. 69. von Wintersdorff CJ, Penders J, van Niekerk JM, Mills ND, Majumder S, van Alphen LB, et al. 2016. Dissemination of antimicrobial resistance in microbial ecosystems through horizontal gene transfer. Front. Microbiol. 7: 173. 70. Jolibois B, Guerbet M. 2005. Hospital wastewater genotoxicity. Annal. Occup. Hyg. 50: 189-196. 71. Abed RM, Safi NM, Köster J, De Beer D, El-Nahhal Y, Rullkötter J, et al. 2002. Microbial diversity of a heavily polluted microbial mat and its community changes following degradation of petroleum compounds. Appl. Environ. Microbiol. 68: 1674-1683. 72. Sánchez-Peinado MM, González-López J, Martínez-Toledo MV, Pozo C, Rodelas B. 2010. Influence of linear alkylbenzene sulfonate (LAS) on the structure of Alphaproteobacteria, Actinobacteria, and Acidobacteria communities in a soil microcosm. Environ. Sci. Pollut. Res. 17: 779-790. 73. Obernosterer I, Catala P, Lebaron P, West NJ. 2011. Distinct bacterial groups contribute to carbon cycling during a naturally iron fertilized phytoplankton bloom in the Southern Ocean. Limnol. Oceanogr. 56: 2391-2401. 74. Wang SY, Sudduth EB, Wallenstein MD, Wright JP, Bernhardt ES. 2011. Watershed urbanization alters the composition and function of stream bacterial communities. PLoS One 6: e22972. 75. Thomas F, Hehemann JH, Rebuffet E, Czjzek M, Michel G. 2011. Environmental and gut bacteroidetes: the food connection. Front. Microbiol. 2: 93. 76. Lee JS, Lee KC, Kim KK, Lee B. 2016. Complete genome sequence of the Variibacter gotjawalensis GJW-30T from soil of lava forest, Gotjawal. J. Biotechnol. 218: 64-65. 77. Zhang B, Xu X, Zhu L. 2017. Structure and function of the microbial consortia of activated sludge in typical municipal wastewater treatment plants in winter. Sci. Rep. 7: 17930. 78. Normand P, Daffonchio D, Gtari M. 2014. The family geodermatophilaceae. The Prokaryotes 2014: 361-379. 79. Song W, Qi R, Zhao L, Xue N, Wang L, Yang Y. 2019. Bacterial community rather than metals shaping metal resistance genes in water, sediment and biofilm in lakes from arid northwestern China. Environ. Pollut. 254: 113041. 80. Liang B, Wang L-Y, Zhou Z, Mbadinga SM, Zhou L, Liu J-F, et al. 2016. High frequency of Thermodesulfovibrio spp. and Anaerolineaceae in association with Methanoculleus spp. in a long-term incubation of n-alkanes-degrading methanogenic enrichment culture. Front. Microbiol. 7: 1431. 81. Wang C, Liu S, Zhang Y, Liu B, He F, Xu D, et al. 2018. Bacterial communities and their predicted functions explain the sediment nitrogen changes along with submerged macrophyte restoration. Microb. Ecol. 76: 625-636. 82. Quiñónez-Dìaz MdJ, Karpiscak MM, Ellman ED, Gerba CP. 2001. Removal of pathogenic and indicator microorganisms by a constructed wetland receiving untreated domestic wastewater. J. Environ. Sci. Heal A. 36: 1311-1320. 83. Sagy M, Kott Y. 1990. Efficiency of rotating biological contactor in removing pathogenic bacteria from domestic sewage. Water Res. 24: 1125-1128. 84. Wong K, Shaw TI, Oladeinde A, Glenn TC, Oakley B, Molina M. 2016. Rapid microbiome changes in freshly deposited cow feces under field conditions. Front. Microbiol. 7: 500. 85. Daims H, Maixner F, Lücker S, Stoecker K, Hace K, Wagner M. 2006. Ecophysiology and niche differentiation of Nitrospira-like bacteria, the key nitrite oxidizers in wastewater treatment plants. Water Sci. Technol. 54: 21-27. 86. Wang J, Chu L. 2016. Biological nitrate removal from water and wastewater by solid-phase denitrification process. Biotechnol. Adv. 34: 1103-1112. 87. Moeder M, Carranzadiaz O, Lopezangulo G, Avina R, Chavez-Duran, A, Jomaa S, et al. 2017. Potential of vegetated ditches to manage organic pollutants derived from agricultural runoff and domestic sewage: A case study in Sinaloa (Mexico). Sci. Total Environ. 598: 1106-1115. 88. Chen Z, Ying G, Liu Y, Zhang Q-Q, Zhao JL, Liu SS, et al. 2014. Triclosan as a surrogate for household biocides: An investigation into biocides in aquatic environments of a highly urbanized region. Water Res. 58: 269-279. 89. Kawahigashi H, Hirose S, Ohkawa H, Ohkawa Y. 2008. Transgenic rice plants expressing human P450 genes involved in xenobiotic metabolism for phytoremediation. J. Mol. Microbiol. Biotechnol. 15: 212-219. 90. Guo G, Tian F, Ding K, Wang L, Liu T, Yang F. 2017. Effect of a bacterial consortium on the degradation of polycyclic aromatic hydrocarbons and bacterial community composition in Chinese soils. Int. Biodeter. Biodegr. 123: 56-62. 91. Chang BV, Liu JH, Liao CS. 2014. Aerobic degradation of bisphenol-A and its derivatives in river sediment. Environ. Technol. 35: 416- 424. 92. Kang J-H, Kondo F. 2002. Bisphenol A degradation by bacteria isolated from river water. Arch. Environ. Contam. Toxicol. 43: 0265- 0269.

J. Microbiol. Biotechnol.