Skip to main content

The bacterial community composition and its environmental drivers in the rivers around eutrophic Chaohu Lake, China



Bacterial community play a key role in environmental and ecological processes in river ecosystems. Rivers are used as receiving body for treated and untreated urban wastewaters that brings high loads of sewage and excrement bacteria. However, little is known about the bacterial community structure and functional files in the rivers around the eutrophic Chaohu Lake, the fifth largest freshwater lake in China, has been subjected to severe eutrophication and cyanobacterial blooms over the past few decades. Therefore, understanding the taxonomic and functional compositions of bacterial communities in the river will contribute to understanding aquatic microbial ecology. The main aims were to (1) examine the structure of bacterial communities and functional profiles in this system; (2) find the environmental factors of bacterial community variations.


We studied 88 sites at rivers in the Chaohu Lake basin, and determined bacterial communities using Illumina Miseq sequencing of the 16 S rRNA gene, and predicted functional profiles using PICRUSt2. A total of 3,390,497 bacterial 16 S rRNA gene sequences were obtained, representing 17 phyla, and 424 genera; The dominant phyla present in all samples were Bacteroidetes (1.4-82.50 %), followed by Proteobacteria (12.6–97.30 %), Actinobacteria (0.1–17.20 %). Flavobacterium was the most numerous genera, and accounted for 0.12–80.34 % of assigned 16 S reads, followed by Acinetobacter (0.33–49.28 %). Other dominant bacterial genera including Massilia (0.06–25.40 %), Psychrobacter (0-36.23 %), Chryseobacterium (0.01–22.86 %), Brevundimonas (0.01–12.82 %), Pseudomonas (0-59.73 %), Duganella (0.08–23.37 %), Unidentified Micrococcaceae (0-8.49 %). The functional profiles of the bacterial populations indicated an relation with many human diseases, including infectious diseases. Overall results, using the β diversity measures, coupled with heatmap and RDA showed that there were spatial variations in the bacterial community composition at river sites, and Chemical oxygen demand (CODMn) and (NH4+ )were the dominant environmental drivers affecting the bacterial community variance.


The high proportion of the opportunistic pathogens (Acinetobacter, Massilia, Brevundimonas) indicated that the discharge of sewage without adequate treatment into the rivers around Chaohu Lake. We propose that these bacteria could be more effective bioindicators for long-term sewage monitoring in eutrophic lakes.

Peer Review reports


Rivers are the primary receiver of organic matter and nutrients from terrestrial ecosystems, and play a key role in biogeochemical cycles in aquatic ecosystems [1]. Although many studies emphasize the importance of river ecosystem services, most river worldwide have suffered serious deterioration, caused primarily by rapid industrialization and urbanization. Therefore, investigation of the overall ecological condition of the river ecosystems is of paramount importance.

In watersheds, lakes and their input rivers are highly linked in multiple ways, yet microbial diversity in river is less commonly studied than lake ecosystems [2]. As a crucial constituent of the river ecosystem, microbes are widely distributed in water column and are diverse in terms of numbers of species. They play a key role in the mineralization of organic matter, and biogeochemical processes [3, 4]. Although many studies have focused on bacterial community composition (BCC) in the water body of lakes [5,6,7,8,9,10], BCC in the input rivers around lakes has not been examined in as much detail. Futhermore, the bacterioplankton assemblage composition in lakes could be affected by the input of allochthonous bacteria [11]. In addition, rivers are also often used as receiving body for treated and untreated urban wastewaters [12], which brings high loads of sewage and excrement bacteria [13]. These bacterial genera usually include waterborne pathogens which are a danger to human health [14,15,16]. Nevertheless, the taxonomic and functional compositions of bacterial communities, and the influencing factors in the river around lakes have been largely ignored and the investigations are crying needed.

Our attention has been drawn to Chaohu Lake, the fifth largest freshwater lake in China, which is located in the downstream of the Yangtze River, which will serves many social, economic, and ecological purposes in the drainage basin [17, 18]. The western lake region receives major inflows, including the Nanfei and Shiwuli rivers (both have sewage outfalls), the Hangbu, and the Pai river. These western rivers account for almost 60 % of the total runoff volume contributed annually to the lake. The eastern lake region connects to the Yuxi river, an only outflowing river, which is the only channel connecting the eastern region to the Yangtze River. However, in the past decades, there has been increasing industrial and agricultural pollution and other strong human activities, causing serious deterioration in water quality of the lake, and increasing coverage and duration of cyanobacterial blooms [19]. Extensive research has been directed to the causes and mechanisms of eutrophication in Chaohu Lake [20]. However, no studies have characterized the bacterial community composition and functional profiles in the rivers around Chaohu Lake, nor has community structure been correlated with environmental factors.

In this study, excepted for investigating the composition of the bacterial community in the rivers around Chaohu Lake, we further studied which one is the environmental determinant of variations in the river bacterial communities. This work could decipher the spatial distribution patterns of BCC and its functional profiles in rivers around Chaohu Lake, and will be significant to understand the microbial ecology of the rivers and assessing ecological risk, as well as provide a scientific basis for the ameliorating pollution of the freshwater lake.


Environmental characterization

The means, and maximum and minimum values, for the 11 environmental parameters measured at the 88 sampling sites in the rivers around the Chaohu Lake basin are summarized in Table 1. During the sampling period, the water temperature varies from 6 to 9 °C. Measurements related to trophic status typically varied greatly among sites. For example, TN concentration ranged from 0.71 mg L− 1 at C45 to 18.80 mg L− 1 at C25 (mean = 4.00 mg L− 1), TP concentration ranged from 0.03 mg L− 1 at C45 to 3.00 mg L− 1 at C25 (mean = 0.22 mg L− 1) and CODMn concentration ranged from 1.43 mg L− 1 at C81 to 16.71 mg L− 1 at C25 (mean = 5.02 mg L− 1). Concentrations of TN, TP and CODMn and other nutrients in the water column showed high levels of pollution in most of the rivers around Chaohu Lake. All water samples had a pH greater than 7, with the highest value of 9.66, indicating alkaline conditions in the study area (Table 1).

Table 1 Mean, and maximum and minimum values, for 12 physicochemical parameters for the 88 sampling sites in the rivers on 15 February 2018. Abbreviations are, Temp water temperature, TN total nitrogen, DTN dissolved total nitrogen, NH4+ ammonium, NO2 nitrite, NO3 nitrate, TP total phosphorus, DTP dissolved total phosphorus, PO4+ orthophosphate, CODMn Chemical Oxygen Demand, BA Bacterial abundance

Patterns of bacterial abundance and diversity

Bacterial abundance varied greatly among the 88 sites in the rivers around the Chaohu Lake basin (Fig. 1). The bacterial abundance was lowest at site C64 with an abundance of only 1.52 × 104 cells mL− 1, and highest in site C25 with an abundance of 3.86 × 107 cells mL− 1; across all sites the mean abundance was 1.55 × 106 cells mL− 1 (Table 1; Fig. 1).

Fig. 1
figure 1

Spatial distribution of bacterial abundance at the 88 sites in the rivers around Chaohu Lake

Based on the identity level of 97 %, the 3,390,497 high quality sequence reads were classified into different OTUs after quality control. Among the 3,390,497 sequence reads, 125,901 OTUs were classified at the phylum level (Fig. 2). Rarefaction curves suggesting that the sequencing effort was sufficient to capture the community diversity (Fig. S1).

Fig. 2
figure 2

Spatial distribution of the OTUs and bacterial α-diversity, represented by Shannon, Richness and Pielou indices, of the 88 sites in the rivers around Chaohu Lake

The bacterial α-diversity patterns, including the Richness, Shannon, and Pielou indexes were distinct between the 88 sites along the rivers (Fig. 2). The Shannon index ranged from 3.471 to 5.860, and was highest in sample C6. The Pielou index ranged from 0.02515 to 0.1378, and was highest in sample C6. The Richness index ranged from 79.74 to 308.30, with the highest value in sample C70.

We investigated the β-diversity patterns of the bacterial communities by employing non-metric multidimensional scaling (NMDS). Notably, these groups were scattered on the NMDS plot, suggesting lower similarity of bacterial community compositions among the 88 sites (Fig. 3). This result was also supported by the hierarchical cluster analysis, indicateing that the 88 samples were separated by site condition (Fig. S2).

Fig. 3
figure 3

Non-metric multi-dimensional scaling (NMDS) plot based on the UniFrac weighted distance calculated from read numbers among the 88 sites around Chaohu Lake

Phylogenetic composition of bacterial community

To visualize the bacterial community compositions in the rivers of the Chaohu Lake basin, the profiles of all taxa at all 88 sites were plotted, and are show in Fig. 4 (phylum level) and Fig. 5 (genus level). A total of 17 different phyla were observed, and the dominant bacterial phyla (those with ≥ 5 % relative abundance in any sample, average values of 88 sites) belonged to Bacteroidetes (51.6 %, average relative abundance), Proteobacteria (38.3 %), and Actinobacteria (5.5 %), which together accounted for 95.4 % of the bacterial sequences; the other phyla accounted for a low fraction of the average relative abundance (4.6 %). Among Proteobacteria, the most abundance class was Gammaproteobacteria, comprising 15.8 % of the total sequences, followed by Betaproteobacteria (13.5 %), and Alphaproteobacteria (9.0 %).

Fig. 4
figure 4

Profile of bacterial community composition, at the phylum level, in the 88 sites of the rivers around Chaohu Lake

Fig. 5
figure 5

Profile of bacterial community composition, at the genus level, in the 88 sites of the rivers around Chaohu Lake

All reads were classified into 424 genus-level taxonomic groups, and 20 of them accounted for 80.89 % of all reads (Fig. 5). Among the 20 abundant genera, there were two genera of Bacteroidetes (48.3 %), followed by six Proteobacteria genera comprising α- (one genus), β- (two genera) and γ- subdivisions (three genera), and one genus of Actinobacteria. Deep taxonomic analyses showed that prominent genera (those with ≥ 1.4 % relative abundance in any sample, average values of 88 sites) consisted of Flavobacterium (45.6 %, average relative abundance), Acinetobacter (8.8 %), Massilia (4.6 %), Psychrobacter (3.2 %), Chryseobacterium (2.7 %), Brevundimonas (2.6 %), Pseudomonas (2.3 %), Duganella (2.1 %), and unidentified Micrococcaceae (1.4 %), which together accounted for 73.3 % of the bacterial sequences; other genera accounted for only 26.7 %. Heatmaps of the 50 most abundant bacterial genera based on the relative abundances indicated that more differences of bacterial community composition were observed between each site at the genus level (Fig. 6).

Fig. 6
figure 6

Heatmap profile showing 50 representative, predominant, 16 S rRNA gene-based sequences, classified at the genus level, using the Silva rRNA database

Functional profiles of the bacterial community

The results of PICRUSt program based on the KEGG classification showed that the predominant predicted bacterial functions were those related to, metabolism, genetic information processing, environmental information processing, human diseases and organismal systems (Fig. 7). Due to the importance of water quality for human health, we further targeted bacterial functional classes related to human diseases. We found that predicated functions were related to six kinds of diseases, cardiovascular, infectious, immune, metabolic, neurodegenerative, and cancer (Fig. 7). Involvement with infectious diseases was the most dominant (Fig. 7).

Fig. 7
figure 7

The relative abundance of predicted functions of bacterial communities in the rivers around Chaohu Lake, calculated using PICRUSt2 grouped into level-2 functional categories

Linking bacterial community to environmental parameters

To assess the correlation of community structure with environmental parameters, we performed a redundancy analysis (RDA) biplot of the BCC of 88 sites and 11 physicochemical parameters (temp, pH, TN, DTN, NH4+, NO2, NO3N, TP, DTP, PO4+, and CODMn) (Fig. 8). The plot demonstrated that CODMn and NH4+ played a significant role in the spatial changes of bacterial communities of the rivers (Monte Carlo test P < 0.05). The first RDA dimension explained 8.9 % of the variation of bacterial communities, and the second explained 5.4 %.

Fig. 8
figure 8

Redundancy analysis ordination plots show the relationship between two key environmental variables (CODMn and NH4+) and bacterial community structure, as reviewed by assessing 16 S rRNA gene OTUs. The nominal variable ‘‘sampling site’’ is represented as centroids. Environmental variables that are significant are shown by * (P < 0.05), and *** if very significant (P < 0.01)


Dominant taxa groups in the BCC in the rivers

Bacteroidetes, Proteobacteria and Actinobacteria were the most abundant phyla in the rivers around the Chaohu Lake, China. The sum of the average relative abundances of these three phylum exceeded 95 % of the total community. Among them, Bacteroidetes occupied 51.6 %. These bacteria play key role in the degradation of protein, cellulose, pectin and chitin, which are typical components of the high molecular mass fraction of dissolved organic matter [7, 21, 22]. The high abundance of Bacteroidetes in the rivers may be related to the organic matter contamination and algal blooms from the study area [18, 23]. Previous study in the Tama River also showed that the soaring numbers of Bacteroidetes may be related to cyanobacterial blooms [7] and river pollution levels [24]. In terms of genera, Flavobacterium was the most dominant one (45.6 % among bacteria) in this group (Fig. 2), which has been related to harmful algal blooms because of its algicidal activity and antagonism to other bacteria [25]. Previous study also indicated that Flavobacterium was often found in high abundance in eutrophic and hypertrophic urban rivers [26].

Proteobacteria was another predominant phylum, with Gammaproteobacteria, Betaproteobacteria, and Alphaproteobacteria predominating in all samples. Proteobacteria is the largest phylum of bacteria, and is involved in a variety of biogeochemical processes in aquatic ecosystems [27, 28]. Deep taxonomic analyses revealed that the genera Acinetobacter, Massilia, and Brevundimonas were the main components of the Gammaproteobacteria, Betaproteobacteria, and Alphaproteobacteria respectively, accounting for 8.8 %,4.6 and 2.6 % of all reads. Examination of the scientific literature showed they are the opportunistic pathogens and infectious [14, 18, 29]. In this study, the high abundance of Acinetobacter, Massilia, Brevundimonas may be related to failure in sewage treatment processes that may be reflected in the high concentration of TN、TP and CODMn of the river water environment. Previous studies also showed that the occurrence of pathogenic bacteria in river following failure in sewage treatment processes may increase near large urban populations [30, 31].

Actinobacteria was the third dominant phylum in our study. Actinobacteria are widely distributed and are the most abundant phylum in freshwater ecosystems [32, 33]. Deep taxonomic analyses showed that the main component of the Actinobacteria was the genus “unclassified Micrococcaceae”, which accounted for 88.1 % of the relative abundance. In total, 1.4 % of them were “unclassified Micrococcaceae within the family Micrococcaceae”. This finding is of special local note, and may reflect a particular local food. In the processing of dry-fermented sausages, Micrococcaceae are the crucial microorganisms used as starter cultures, and are used in the preservation of meat products that avoid rancidness and develop the typical red colour owing to catalase and nitrate reductase [34]. The Chaohu Lake lies in the Anhui Province of China, where preserved bacon is local speciality food [35].

Functional profiles of bacteria in the rivers

Understanding the functional profiles of bacterial communities is of great importance because it may shed light on ecosystem processes and community assembly mechanisms [1]. In our present study, functional profiles using PICRUSt revealed that the bacteria in the rivers of the Chaohu lake were involved in many diverse pathways (Fig. 4), most of which were related to metabolism systems, genetic information processing and environmental information processing. In addition, a considerable group of bacteria were involved in human diseases, including infectious, neurogenerative, metabolic, immune system, cardiovascular diseases and cancer. Among these, the most dominant was involvement in infectious diseases. This finding is consistent with previous finding in the Apies River, in South Africa [36].

The occurrence, and increase of bacteria associated with human diseases in river water may be related to the release of effluent from wastewater treatment plants, ineffective septic tank systems, and storm water runoff [37]. Although human pathogenic bacteria often occur at low levels in the environment [38], the pollutants that enter these systems may settle into the bottom sediments, and pathogenic bacteria will gradually increase in numbers, and finally increase the risk of infections to humans and animals [36]. Therefore, effective governance of untreated sewage in an urban river is of paramount importance.

Response of the community structure of river bacteria to CODMn and NH4 +

We found that there were significant differences in bacterial community composition among the 88 sites in the rivers (Figs. 3 and 6). This agrees with other findings that there is variation largely in the taxonomic composition and spatial distribution of freshwater bacterial community among different rivers [39,40,41,42,43]. Our RDA results revealed that CODMn and NH4+ were the most significant determining factors related to the variances in the river bacterial community around the Chaohu Lake. These may relate strongly to the significant variations of CODMn and NH4+ concentrations at different rivers. A correlation between CODMn and NH4+ and bacterial community structure has been reported previously. In a 2008 study, Wei et al. found that bacterial community structures in Chaohu Lake were influenced significantly by the influent CODMn [9]. Previous studies in China have also revealed that CODMn was significantly related to the bacterial composition diversity in an urban river [44], and CODMn and NH4+ were significant structuring factors for bacterial community compositions in an urban lake [45].


In summary, using Illumina miseq sequencing, we explored bacterial community diversity, composition and functional profiles of 88 sites in the rivers around Chaohu Lake, China. The results showed that Bacteroidetes, Proteobacteria and Actinobacteria were the dominant phyla, and dominant genera included Flavobacterium, Acinetobacter, Massilia, Psychrobacter, Chryseobacterium, Brevundimonas, Pseudomonas, Duganella, and Unidentified Micrococcaceae. The functional profiles of the bacterial populations revealed an association with many human diseases, including infectious diseases. We also found site differences in the bacterial community structure in river water, CODMn and NH4+ were the main drivers regulating these variations. Our results indicated that the discharge of sewage without adequate treatment into the rivers around the Chaohu Lake owing to widespread occurrence of pathogenic bacteria (Acinetobacter, Massilia, Brevundimonas), and these bacteria could be more effective bioindicators for long-term sewage monitoring in eutrophic lakes. Therefore, the capacity of sewage treatment needs to be substantially strengthened around the Chaohu watershed to protect Chaohu Lake from further contamination.


Study area and sampling

Chaohu Lake (31°25′- 31°43′N, 117°16′- 117°51′E), is located in the center of Anhui Province, China, and in the downstream of the Yangtze River. The lake has a surface area of 760 km2 and can be divided into two regions, from the Zhongmiao Temple to Qitouzui Cape (Fig. 9), the western region is eutrophic to hypertrophic, and the eastern region is mesotrophic. The western region receives major inflows, including the Nanfei and Shiwuli rivers (both have sewage outfalls), the Hangbu, and the Pai river. These western rivers account for almost 60 % of the total runoff volume contributed annually to the lake. The eastern region connectes to the Yuxi river, which is the only channel connecting the eastern region to Yangtze River, permitting water exchange [18, 46].

Fig. 9
figure 9

Map showing the location of the 88 sites in the rivers around Chaohu Lake

On 15 February 2018 we carried out our field work at 88 sites of the rivers around the Lake Chaohu basin (Fig. 9). At each sites, surface water (top 50 cm) was collected with a 5 L Schindler sampler. For 16 S rRNA gene analysis, a subsample of water (500 mL) was pre-filtered in situ with polycarbonate membrane (0.2 μm pore-size and 47 mm diameter, Millipore) using a hand-driven vacuum pump. These filters were frozen at -80 oC until DNA extraction was performed. For enumeration of bacterial abundance, an additional subsample (46 mL) was transferred into an autoclaved tube containing 4 mL of prefiltered (pore size, 0.2 μm) glutaraldehyde (final concentration 2 % [v/v]). These samples were then stored in a refrigerator at 4 oC, until slides were prepared for enumeration of bacterial abundance. The remaining water samples were transported to the laboratory in dark cooling boxes, and processed 3–5 h after sampling within 4 h for immediate chemical analysis.

Measurement of environmental parameters and bacterial abundance

Water temperature (Temp) and pH were measured in situ using a multi-parameter water quality sonde (YSI 6600V2, Yellow Springs Instruments, USA). Chemical analyses of water samples for nine parameters (total nitrogen (TN), dissolved total nitrogen (DTN), ammonium (NH4+), nitrate (NO3), nitrite (NO2), total phosphorus (TP), dissolved total phosphorus (DTP), orthophosphate (PO4), and chemical oxygen demand (CODMn) were conducted in the laboratory according to standard methods. The abundance of bacteria in the water samples was determined by the 4’,6’-diamidino-2-phenylindole (DAPI)-combined epifluorescence direct counting method [47].

DNA extraction, PCR amplification and Illumina Miseq sequencing

The total DNA was extracted using proteinase K, sodium dodecy1 sulfate, and cetyltrimethyl ammonium bromide, follow by phenol-chloroform extraction and isopropanol precipitation [48]. Crude DNA extracts were then purified by the E.Z.N.A® cycle-Pure kit (Omega Bio-Tek).

The V4-V5 regions of the 16 S rRNA genes were amplified using the primers 515 F (GTGCCAGCMGCCGCGGTAA) and 907R (CCGTCAATTCMTTTRAGTTT). Polymerase chain reaction (PCR) amplification was performed in a 25 µL reaction mixture containing 5 µL of 5 × PCR buffer, 1 µL of each primer (10 µmol L− 1 each), 20 ng µL− 1 DNA template, and 0.25 µL of Taq polymerase (5 U µL− 1; Fermentas). The PCR cycling was performed in a thermocycler (Applied Biosystems Veriti Thermal Cycler) under the following conditions: 98 °C for 2 min, 28 cycles each at 98 °C for 15 s, 55 °C for 30 s and 72 °C for 30 s; with final extension at 72 °C for 5 min. Sequencing service was performed by an Illumina Miseq platform at Personal Biotechnology Co., Ltd (Shanghai, China).

Sequencing data processing

Sequence reads were processed by the Quantitative Insights Into Microbial Ecology (QIIME) v. 1.9.1 pipeline [49]. After quality controlling, bacterial phylotypes were assigned to operational taxonomic units (OTUs, 97 % cutoff) using the Uclust algorithm [50]. The longest sequence in each cluster was chosen as the representative sequence, which were annotated by the Silva rRNA database project (SILVA VERSION SSU11;

Predictive metagenome analysis

The metagenome functional content was predicted using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) software package (version 2.0.0, [51]. PICRUSt2-compatible OTU tables were made using the closed-reference OTU picking protocol in QIIME against the RDP database. The nearest sequenced taxon index was used as a measure to represent the novelty of bacteria within an OTU table, in respect of previous sequenced genomes. The obtained OTU table was normalized as the true abundance, and applied with default settings to obtain the predicted metagenomics table with Kyoto Encyclopedia of Genes.

Statistical analysis

Bacterial alpha-diversity was processed using the QIIME pipeline. Community differences between sites were visualized by non-metric multidimensional scaling (NMDS) that was performed using the R statistical program [52]. Heat maps of the most abundant bacterial genera were analyzed using a “pheatmap” package in the R environment. Correlations between environmental variables and bacterial communities were measured with ordination methods using the vegan package in the R environment.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. The raw data are available from the corresponding author on reasonable request. Sequence data of this project have been deposited in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI), with accession number SRP189003 (Persistent web link to datasets:



Bacterial community composition


Chemical oxygen demand


Dissolved total nitrogen


Dissolved total phosphorus


Kyoto Encyclopedia of Genes


National Center for Biotechnology Information

NH4 + :



Non-metric multidimensional scaling analysis

NO3 :


NO2 :



Operational Taxonomic Units


Polymerase chain reaction


Permutational Multivariate Analysis of Variance


Phylogenetic Investigation of Communities by Reconstruction of Unobserved States

PO4 + :



Quantitative Insights Into Microbial Ecology


Total nitrogen


Total phosphorus


  1. Wu HN, Li Y, Zhang WL, Wang C, Wang PF, Niu LH, Du JM, Gao Y. Bacterial community composition and function shift with the aggravation of water quality in a heavily polluted river. J Environ Manage. 2019;237:433–441.

    Article  CAS  PubMed  Google Scholar 

  2. Ren Z, Wang F, Qu XD, Elser JJ, Liu Y, Chu JM. Taxonomic and functional differences between microbial communities in Qinghai lake and its input streams. Front Microbiol. 2017;8:2319.

  3. Cotner JB, Biddanda BA. Small players, large role, microbial influence on biogeochemical processes in pelagic aquatic ecosystems. Ecosystems. 2002;5:105–121.

    Article  CAS  Google Scholar 

  4. Newton RJ, Bootsma MJ, Morrison HG, Sogin ML, Mclellan SL. A microbial signature approach to identify fecal pollution in the waters off an urbanized coast of lake Michigan. Microb Ecol. 2013;65 (4):1011–1023.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lindströ ES. Bacterioplankton Community Composition in Five Lakes Differing in Trophic Status and Humic Content. Microb Ecol. 2000;40(2):104–113.

    Article  Google Scholar 

  6. Gucht KVD, Sabbe K, Meester LD, Vloemans N, Zwart G, Gillis M, Vyverman W. Contrasting bacterioplankton community composition and seasonal dynamics in two neighbouring hypertrophic freshwater lakes. Environ Microbiol. 2001;3(11):680–690.

    Article  Google Scholar 

  7. Eiler A, Bertilsson S. Composition of freshwater bacterial communities associated with Cyanobacterial blooms in four Swedish lakes. Environ Microbiol. 2004;6:1228–1243.

    Article  PubMed  Google Scholar 

  8. Gucht KVD, Vandekerckhove T, Vloemans N, Cousin S, Muylaert K, Sabbe K, Gillis M, Declerk S, Meester LD, Vyverman W. Characterization of bacterial communities in four freshwater lakes differing in nutrient load and food web structure. FEMS Microbiol Ecol. 2005;53(2):205–220.

    Article  PubMed  CAS  Google Scholar 

  9. Wei CL, Bao SM, Zhu XY, Huang XM. Spatio-temporal variations of the bacterioplankton community composition in Chaohu Lake, China. Prog Nat Sci. 2008;18(9):1115–1122.

    Article  CAS  Google Scholar 

  10. Kurilkina MI, Zakharova YR, Galachyants YP, Petrova DP, Bukin YS, Domysheva VM, Blinov VV, Likhoshway YV. Bacterial community composition in the water column of the deepest freshwater lake baikal as determined by next-generation sequencing. FEMS Microbiol Ecol. 2016;92(7):fiw094.

    Article  PubMed  CAS  Google Scholar 

  11. Lindström ES, Bergström AK. Influence of inlet bacteria on bacterioplankton assemblage composition in lakes of different hydraulic retention time. Limnol Oceanogr. 2004;49(1):125–136.

    Article  Google Scholar 

  12. Abraham WR. Megacities as sources for pathogenic bacteria in rivers and their fate downstream. Int J Microbiol. 2011:798282.

  13. Shanks OC, Newton RJ, Kelty CA, Huse SM, Sogin ML, Mclellan SL. Comparison of the microbial community structures of untreated wastewaters from different geographic locales. Appl Environ Microbiol. 2013;79(9):2906–2913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Pandey PK, Kass PH, Soupir ML, Biswas S, Singh VP. Contamination of water resources by pathogenic bacteria. AMB Express. 2014;4(51):2–16.

    Google Scholar 

  15. Rochelle-Newall E, Nguyen TMH, Le TPQ, Sengtaheuanghoung O, Ribolzi O. A short review of fecal indicator bacteria in tropical aquatic ecosystems, knowledge gaps and future directions. Front Microbiol. 2015;6:308.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Wang P, Chen B, Yuan RQ, Li CQ, Li Y. Characteristics of aquatic bacterial community and the influencing factors in an urban river. Sci Total Environ. 2016;569–570:382–389.

    Article  PubMed  CAS  Google Scholar 

  17. Chen SJ, Luo XJ, Lin Z, Luo Y, Li KC, Peng XZ, Mai BX, Ran Y, Zeng EY. Time trends of polybrominated diphenyl ethers in sediment cores from the Pearl River estuary, south China. Environ Sci Technol. 2007;41(16):5595–5600.

    Article  CAS  PubMed  Google Scholar 

  18. Huo SL, Li CC, Xi BD, Yu ZQ, Yeager KM, Wu FC. Historical record of polychlorinated biphenyls (PCBs) and special occurrence of PCB 209 in a shallow fresh-water lake from eastern China. Chemosphere. 2017;184:832–840.

    Article  CAS  PubMed  Google Scholar 

  19. Cai YJ, Gong ZJ, Xie P. Community structure and spatiotemporal patterns of macrozoobenthos in Lake Chaohu (China). Aquat Biol. 2012;17:35–46.

    Article  Google Scholar 

  20. Zhang Y, Cheng L, Tolonen KE, Yin HB, Gao JF, Zhang ZM, Li KY, Cai YJ. Substrate degradation and nutrient enrichment structuring macroinvertebrate assemblages in agriculturally dominated lake Chaohu basins, China. Sci Total Environ. 2018;627:57–66.

    Article  CAS  PubMed  Google Scholar 

  21. Kirchman DL. The ecology of Cytophaga-Flavobacteria in aquatic environments. FEMS Microbiol Ecol. 2002;39:91–100.

    CAS  PubMed  Google Scholar 

  22. Eiler A, Bertilsson S. Flavobacteria blooms in four eutrophic lakes, Linking population dynamics of freshwater bacterioplankton to resource availability. Appl Environ Microbiol. 2007;73:3511–3518.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Araya R, Tani K, Takagi T, Yamaguchi N, Nasu M. 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. 2003;43 (1):111–119.

    Article  CAS  PubMed  Google Scholar 

  24. Ghai R, Rodriguez-Valera F, McMahon KD, Toyama D, Rinke R, de, Oliveira TCS, Wagner GJ, deMiranda FP, Henrique-Silva F. Metagenomics of the water column in the pristine upper course of the Amazon river. PLoS One. 2011;6:e23785.

  25. Roth PA, Mikulski CM, Doucette GJ. The influence of microbial interactions on the susceptibility of Karenia spp. to algicidal bacteria. Aquat Microb Ecol. 2008;50:251–259.

    Article  Google Scholar 

  26. Drury B, Rosi-Marshall E, Kelly JJ. Wastewater treatment effluent reduces the abundance and diversity of benthic bacterial communities in urban and suburban rivers. Appl Environ Microbiol. 2013;79:1897–1905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Wang L, Zhang J, Li HL, Yang H, Peng C, Peng ZS, Lu L. Shift in the microbial community composition of surface water and sediment along an urban river. Sci Total Environ. 2018;627:600–612.

    Article  CAS  PubMed  Google Scholar 

  28. Liu WX, He W, Wu JY, Qin N, He QS, Xu FL. Residues, bioaccumulations and biomagnification of perfluoroalkyl acids (PFAAs) in aquatic animals from Lake Chaohu, China. Environ Poll. 2018;240:607–614.

    Article  CAS  Google Scholar 

  29. Bergogne-Bérézin E, Towner KJ. Acinetobacter spp. as nosocomial pathogens, microbiological, clinical, and epidemiological features. Clin Microbiol Rev. 1996;9(2:148–165.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Lindquist D, Murrill D, Burran WP, Winans G, Janda JM, Probert W. Characteristics of Massilia timonae and Massilia timonae-like isolates from human patients, with an emended description of the species. J Clin Microbiol. 2003;41:192–196.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Brown SP, Cornforth DM, Mideo N. Evolution of virulence in opportunistic pathogens, generalism, plasticity, and control. Trends Microbiol. 2012;20(7):336–342.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Crump BC, Hobbie JE. Synchrony and seasonality in bacterioplankton communities of two temperate rivers. Limnol Oceanogr. 2005;50:1718–1729.

    Article  CAS  Google Scholar 

  33. Wilhelm SW, Lecleir GR, Bullerjahn GS, McKay RM, Saxton MA, Twiss MR, Bourbonniere RA. Seasonal changes in microbial community structure and activity imply winter production is linked to summer hypoxia in a large lake. FEMS Microbiol Ecol. 2014;87:475–485.

    Article  CAS  PubMed  Google Scholar 

  34. Alberto M, Blanca C, Aranda E, María JB, María GC. Characterization of Micrococcaceae from iberian dry-cured sausages. Meat Sci. 2007;75(4):696–708.

    Article  CAS  Google Scholar 

  35. Su CS, Wang YL, Min LJ, Yang YH, Bei YQ, Wei CG, Zhang SJ. Exchange of experience in meat processing technology–Cured meant products. Food Science (China). 1981;2(7):59–64.

    Google Scholar 

  36. Abia A, Alisoltani A, Keshri J, Ubombajaswa E. 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. 2017;326:616–617.

    Google Scholar 

  37. Hong H, Qiu J, Liang Y. Environmental factors influencing the distribution of total and fecal coliform bacteria in six water storage reservoirs in the Pearl River Delta region, China. J Environ Sci-China. 2010;22:663–668.

    Article  PubMed  Google Scholar 

  38. Straub TM, Chandler DP. Towards a unified system for detecting waterborne pathogens. J Microbiol Meth. 2003;53:185–197.

    Article  Google Scholar 

  39. Bai Y, Qi W, Liang J, Qu J. Using high-throughput sequencing to assess the impacts of treated and untreated wastewater discharge on prokaryotic communities in an urban river. Appl Environ Microbiol. 2014;98:1841–1851.

    CAS  Google Scholar 

  40. Garcia-Armisen T, Inceoglu O, Ouattara NK, Anzil A, Verbanck MA, Brion N, Servais P. Seasonal variations and resilience of bacterial communities in a sewage polluted urban river. PLoS One. 2014;9 (3):e92579.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Ouattara NK, Garcia-Armisen T, Anzil A, Brion N, Servais P. Impact of wastewater release on the faecal contamination of a small Urban River, the Zenne River in Brussels (Belgium). Water Air Soil Poll. 2014;225:2043.

    Article  CAS  Google Scholar 

  42. Ruiz-González C, Niño-García JP, Giorgio PA. Terrestrial origin of bacterial communities in complex boreal freshwater networks. Ecol Lett. 2015;18:1198–1206.

    Article  PubMed  Google Scholar 

  43. Ibekwe AM, Ma J, Murinda SE. Bacterial community composition and structure in an urban river impacted by different pollutant sources. Sci Total Environ. 2016;566–567:1176–1185.

    Article  PubMed  CAS  Google Scholar 

  44. Tang FY, Zhang QF, Zhou YJ, Chen LX, Xu JR, Cheng JR. Impacts of ecological restoration treatments on diversity of prokaryotic bacteria in urban river. Chinese J Environ Eng. 2016;10(1):445–454.

    CAS  Google Scholar 

  45. Zhang HH, Wang Y, Chen SN, Zhao ZF, Feng J, Zhang ZH, Lu KY, Jia JY. Water bacterial and fungal community compositions associated with urban lakes, xi’an, China. Int J Environ Res Pub Heal. 2018;15(3):469.

    Article  CAS  Google Scholar 

  46. Shang GP, Shang JC. Spatial and temporal variations of eutrophication in western Chaohu Lake, China. Environ Monit Assess. 2007;130(1–3):99–109.

    Article  CAS  PubMed  Google Scholar 

  47. Tang XM, Gao G, Qin BQ, Zhu LP, Chao JY, Wang JJ, Yang GJ. Characterization of bacterial communities associated with organic aggregates in a large, shallow, eutrophic freshwater lake (lake Taihu, China). Microb Ecol. 2009;58(2):307–322.

    Article  CAS  PubMed  Google Scholar 

  48. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning, A laboratory manual, pp. 25–28. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, 1989, New York, N.Y., USA.

  49. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–336.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. Uchime improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–2200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Thurber RLV, Knight R, Beiko RG, Huttenhower C. Predictive functinal profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan, Community Ecology Package, R Package Version 2.0–10, 2013, Available from:

Download references


We are also grateful to Hao Liu for his assistance in field sampling work.


This study was supported by the National Natural Science Foundation of China (Grant No. 41977322), the Major Science and Technology Program for Water Pollution Control and Treatment (2017ZX07603-001), the “One-Three-Five” Strategic Planning of Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences (Grant No. NIGLAS2017GH05, NIGLAS2017GH01, NIGLAS2017GH02), the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (QYZDJSSWDQC008), and the Key Program of the Chinese Academy of Sciences (ZDRW-ZS-2017-3-4).

Author information

Authors and Affiliations



KS, QX and GG conceived and designed the experiments. XY, ZW, XJ and YH performed sample collection and analyzed the data; KS, and XT wrote the paper. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Guang Gao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1. 

Rarefaction curves of the number of operational taxonomic units (OTUs) at 97% similarity boxplot for each of 88 samples. Figure S2. UPGMA result based on the unweighted Unifrac metric. The hierarchical clustering structure helps to determine the similarity of the bacterial communities between different samples.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shao, K., Yao, X., Wu, Z. et al. The bacterial community composition and its environmental drivers in the rivers around eutrophic Chaohu Lake, China. BMC Microbiol 21, 179 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: