Integrated biogeography of planktonic and sedimentary bacterial communities in the Yangtze River
- Tang Liu†1,
- An Ni Zhang†2,
- Jiawen Wang†1,
- Shufeng Liu1,
- Xiaotao Jiang2,
- Chenyuan Dang1,
- Tao Ma1,
- Sitong Liu1,
- Qian Chen3,
- Shuguang Xie1,
- Tong Zhang2Email author and
- Jinren Ni1, 4Email authorView ORCID ID profile
© The Author(s). 2018
Received: 4 April 2017
Accepted: 17 December 2017
Published: 19 January 2018
Bacterial communities are essential to the biogeochemical cycle in riverine ecosystems. However, little is presently known about the integrated biogeography of planktonic and sedimentary bacterial communities in large rivers.
This study provides the first spatiotemporal pattern of bacterial communities in the Yangtze River, the largest river in Asia with a catchment area of 1,800,000 km2. We find that sedimentary bacteria made larger contributions than planktonic bacteria to the bacterial diversity of the Yangzte River ecosystem with the sediment subgroup providing 98.8% of 38,906 operational taxonomic units (OTUs) observed in 280 samples of synchronous flowing water and sediment at 50 national monitoring stations covering a 4300 km reach. OTUs within the same phylum displayed uniform seasonal variations, and many phyla demonstrated autumn preference throughout the length of the river. Seasonal differences in bacterial communities were statistically significant in water, whereas bacterial communities in both water and sediment were geographically clustered according to five types of landforms: mountain, foothill, basin, foothill-mountain, and plain. Interestingly, the presence of two huge dams resulted in a drastic fall of bacterial taxa in sediment immediately downstream due to severe riverbed scouring. The integrity of the biogeography is satisfactorily interpreted by the combination of neutral and species sorting perspectives in meta-community theory for bacterial communities in flowing water and sediment.
Our study fills a gap in understanding of bacterial communities in one of the world’s largest river and highlights the importance of both planktonic and sedimentary communities to the integrity of bacterial biogeographic patterns in a river subject to varying natural and anthropogenic impacts.
Rivers play an important part in coupling biogeochemical cycles between continents and oceans . River water flow acts as a carrier of biotic and abiotic substances, whereas sediment serves as a sink or source in the cycling of nutrients. Previous studies on riverine biodiversity mostly focused on invertebrates or fish, and very limited reports were related to bacterial composition along river networks [2–5]. Bacteria hold key roles in microbial communities and contribute significantly to biogeochemical processes and the cycling of nutrients in river ecosystems [6–8]. Recent studies have shown that the biogeography patterns of bacterial communities in complex freshwater networks can be explained by their origins in upstream freshwater  and terrestrial sources . Planktonic bacterial taxa arise as the sum of multiple upstream sources of bacteria that grow in rainfall, lakes, groundwater, and soil. The resulting planktonic bacterial community is vulnerable to changes in its composition and structure. Sedimentary bacterial taxa develop from long-term cumulative processes of sediment erosion and deposition under ambient conditions . The spatiotemporal distribution of planktonic and sedimentary bacterial communities in rivers could be quite different. Moreover, bacterial diversity is significantly altered by varying fluvial landforms and severe human interference along a large river. A better understanding of bacterial responses to the changing environment of river ecosystems is useful in the context of riverine cycles of nutrients, e.g., carbon and nitrogen, which are highly relevant to emission or sequestration of greenhouse gases [12–14].
Until recently, riverine bacterial communities have proved highly diverse and variable. Spatial and temporal variability of bacterioplankton composition in rivers has been demonstrated [2–5, 15–19]. Crump et al.  found that synchronous shifts in the bacterial communities of six Arctic rivers were strongly correlated with seasonal changes in the environment, suggesting that microbial communities may shift in predictable patterns from season to season. Fortunato et al.  detected spatial variability of bacterioplankton communities from the river to ocean that overwhelmed the seasonal trends. However, most of the previous studies considered relatively short river reaches or coastal areas [4, 17–19]. Although two recent investigations examined the variability and diversity of planktonic bacterial communities in longer river, e.g., the in Danube River  and Yenisei River , seasonal patterns in water and sediment bacterial community and the co-occurrence of different populations were still unknown.
An understanding of the biogeography of microorganisms is essential to reveal the specific role of unique species and their links to ecosystem functions that affect global processes . The microbial community assembly can be regulated by the local deterministic processes of environmental selection and interspecific competition (species sorting perspective) and by the regional stochastic processes of extinction, emigration, immigration, and speciation (neutral perspective). Species sorting emphasizes that the difference in local community composition lies along an environmental gradient, not a geographical one. This perspective has much in common with niche theory about niche separation [21–23]. The neutral perspective assumes all species are ecologically equivalent, and community dynamics can be derived from probabilities of immigration and emigration extinction and genetic drift . Various studies have indicated that species sorting or neutral processes are important in microbial community assembly [25–28], but few studies have examined their relative importance in structuring riverine bacterial communities [29, 30].
Full biogeographic patterns of bacterial communities in large rivers should include spatiotemporal patterns of both planktonic and sedimentary communities. For complex aqueous ecosystems in large rivers, a combination of neutral and species sorting perspectives is often needed to interpret bacterial communities. Herein, we present the first integrated biogeography of bacterial communities in water and sediment of the Yangtze River, which is interpreted in terms of landforms and human impacts along the river (Additional file 1: Figure S1).
Richness of bacterial communities
For 280 samples comprising 97 water and 183 sediment samples, a total of 19,733,498 high-quality bacterial 16S rRNA gene sequences and 38,906 operational taxonomic units (OTUs) were obtained by high-throughput sequencing. The rarefaction curves illustrated that the bacterial OTUs obtained by the current sequencing depth were sufficient to represent the microbial communities in water samples; whereas the curves of most of the sediment samples did not reach a plateau (Additional file 2: Figure S2). Good’s coverage values for sediment samples (0.97 ± 0.014, see Additional files 3: Table S1), indicating the primary sedimentary microbial communities were correctly represented by the current bacterial profiles. Moreover, sediment samples had higher OTU richness than water samples, although no significant seasonal difference was found in each type of sample. The sediment subgroup contained 98.8% of the total OTUs, 58.6% of which were unique to sediment. OTUs in water subgroup only accounted for 41.4% of the total OTUs. For sediment samples, the number of OTUs shared between spring and autumn accounted for 85.5% of sediment OTUs, whereas a much smaller portion, only 42.7% of OTUs, were found in the water samples.
Taxonomic composition and season-associated taxa
Of the 38,906 clustered bacteria OTUs, 81.8% (31,842) were assigned to the phylum level, followed by 77.8% (30,286), 63.2% (24,569), 34.6% (13,477), 10.8% (4189), and 0.7% (262) at lower taxonomy levels (i.e., class, order, family, genus, and species, respectively).
To examine the seasonal (autumn and spring) effect on the structure of the bacterial communities, the persistent bacterial OTUs were further clustered into autumn-associated OTUs or spring-associated OTUs, according to whether the OTUs were significantly (P < 0.05) more abundant in autumn or in spring. The results indicate that bacterial communities in both water and sediment samples were enriched during the autumn; the autumn-associated OTUs (red nodes in Fig. 1) accounted for 30.5% (648) and 47.8% (269) of the sediment and water samples, which is higher than the 12.6% (267) and 26.1% (147) of spring-associated OTUs (blue nodes in Fig. 1) in the sediment and water. The phylogenetic distribution of these season-associated OTUs demonstrates that the phyla usually harbored either autumn-associated OTUs or spring-associated OTUs, and seldom a mixture of both types (Fig. 1a, c). This indicates a consistent seasonal response of the season-associated OTUs within the same phylum. Accordingly, the main phyla in the sediment or water can be classified into three groups: (a) phyla dominated by autumn-associated OTUs, including Chloroflexi, Planctomycetes, Acidobacteria, and Actinobacteria (in sediment); (b) phyla dominated by spring-associated OTUs, including Bacteroidetes and Cyanobacteria (in sediment); and (c) phyla containing a mixture of both autumn-associated OTUs and spring-associated OTUs, including Proteobacteria, Cyanobacteria (in water), and Actinobacteria (in water). Overall, the season-associated OTUs within the same phylum generally displayed uniform seasonal preference, except for the super-phyla, Proteobacteria, and Cyanobacteria/Actinobacteria in water, with most phyla preferring autumn than spring. The effect of season on shaping the microbial community structure is consistent with the monitored results (Additional file 4: Figure S3) and co-occurrence of season-associated OTUs within the same phylum in network analysis (Additional file 5: Figure S4, Fig. 1b, d). In each phylum, positive association (red lines) dominates negative association (blue lines), indicating the uniform seasonality response of season-associated OTUs within a phylum. This uniform response is quantified by the network density value (d) (“Season-associated taxa analysis” section). Phyla that respond consistently are identified by having a higher density value implying a denser distribution of the association among OTUs. For example, Verrucomicrobia and Spirochaetes are the two most tightly correlated clusters (d = 0.433 and 0.455, respectively), whereas Acidobacteria in sediment samples and Actinobacteria and Bacteroidetes in water samples are more dispersed. The co-occurrence network analysis reveals almost exclusively or overwhelmingly positive correlations, consistent with the general pattern of dominant season-associated OTUs in most phyla (Fig. 1a, c).
In summary, the results show that microbial communities in Yangtze River are particularly sensitive to the season and are more prevalent in the autumn.
Biogeography patterns of bacterial communities
The Yangtze River flows through various landform types, including plateau, mountain, basin, foothill, and plain. Mainstream samples were used to study the effects of spatial variation, such as the river continuum and landform on bacterial communities. For water bacterial communities, NMDS (Additional file 8: Figure S6 (a) and (b)) gave similar results in both spring and autumn, with five separate groups: group 1 (between stations 1 and 2), group 2 (between stations 3 and 5), group 3 (from stations 6 to 9), group 4 (from stations 10 to 13), and group 5 (from stations 14 to 24), corresponding to the local landforms: mountain, foothill, basin, foothill-mountain, and plain. Moving window analysis (Additional file 9: Figure S7) was used to characterize the change rate of bacterial communities along the mainstream based on a comparison of results between two consecutive sampling sites. Higher change rates were always found at sampling sites where the landform type changed. Meanwhile, ANOSIM analysis (Additional file 10: Figure S8) further illustrated that taxonomic compositions of microbial communities significantly varied by landform type (P = 0.001). For the sediment samples, a similar clustering result (Additional file 8: Figure S6 (c) and (d)) was also obtained; this indicated that bacterial communities from the same landform tended to be similar to each other. The results of these analyses revealed that spatial variation in bacterial compositions across the samples could be partially attributed to the landform.
Influential factors on bacterial community compositions
Even though canonical correspondence analysis (CCA) of bacterial communities in water and sediment indicated weak correlation to environmental factors (Additional file 11: Table S3, and Additional file 12: Figure S9), water temperature was found to be the primary factor in structuring bacterial community assemblages in both water and sediment of the Yangtze River. In addition, DO (dissolved oxygen) influenced the bacterial community in water, and pH, Mdd (mean dendritic distance), and TN (total nitrogen) influenced the bacterial community in sediment.
Impacts of the large dams
Evidence of neutral and species sorting processes
The intensely turbulent flow in the Yangtze River enhances random dispersal of bacterial species in water over long stream-wise distances. However, sediment transport is likely to be much slower and restricted to short stream-wise distance. In fact, fluvial rivers are capable of self-adjustment  and eventually reach equilibrium state through long-term processes of sediment erosion and deposition  with accumulated sediment exerting an ambient influence on both abiotic and biotic diversity [33, 34]. In general, increased habitat heterogeneity in sedimentary environments should lead to more significant spatial differences in bacterial communities. By examining the factors that determine bacterial biogeography, it is possible to investigate whether neutral and/or species sorting processes control the assembly of microbial communities in water and sediment.
The neutral model did not interpret 100% of the community species distribution, indicating that other community assembly mechanisms were perhaps operating at the same time. There is some evidence of species sorting in the present results. For example, the LEfSe analysis indicated that different abundant species or clades (biomarkers) occurred in five landform types (Fig. 4). Canonical correspondence analysis (CCA) of the bacterial communities in water and sediment showed a weak dependence on environmental factors. A significant discrepancy in the abundance of OTUs (Fig. 5) was observed upstream and downstream of the Three Gorges Dam, with most OTUs having higher abundance upstream of the dam except for several OTUs belonging to the genera Anaerolinea and Flavobacterium (Fig. 5). This implies that each sampling site possesses specific species due to its particular environment. The non-random distributions of bacteria are ascribed to the heterogeneous environment that affects their natural habitats and nutrients, influencing selection.
A rapidly increasing number of studies based on high-throughput sequencing technologies have revealed a tremendous diversity of bacterial communities residing in the aquatic environment [38–40]. Most previous investigations on the variability and diversity of bacterial communities in rivers focused on a single dimension, i.e., either a long-term time series or a small-scale spatial dimension across environmental gradients [4, 19]. Here, we describe spatiotemporal patterns of lotic bacterial communities over a 4300 km river continuum for both water and sediment during spring and autumn seasons. Although changing slightly across the seasons, bacterial communities in sediment provide the main contribution to the bacterial diversity of the Yangtze River, and only 1.2% of the total OTUs are unique to water samples. The bacterial population fluctuation in water samples is higher than in sediment samples, as expected. Few previous studies have compared planktonic and sedimentary bacterial diversity in river reaches and coastal areas [15, 16]. Feng et al.  reported that sedimentary communities were taxonomically richer than their planktonic counterparts in the Yangtze River estuary. However, Feng et al.’s water samples displayed a slightly higher bacterial diversity than sediment samples for certain reaches of the Yellow River. In our study of the Yangtze River, taxonomic richness was generally higher in sediment and lower in water, presumably due to the lower concentration of suspended sediment in the Yangtze.
The huge datasets established for this study permit a deeper investigation into temporal and spatial patterns of bacterial communities in water and sediment. Figure 2 is a NMDS diagram from which it can be seen that planktonic communities vary according to sampling seasons, unlike the sedimentary communities. In both spring and autumn, the dominant phyla in the water ecosystem were Proteobacteria, Actinobacteria, Bacteroidetes, and Cyanobacteria. The abundances of Bacteroidetes and Cyanobacteria in spring were lower than in autumn, whereas Proteobacteria and Actinobacteria were more prevalent in autumn than in spring. Notably, the abundance of Cyanobacteria in the mainstream of the Yangtze was higher in spring, especially in the upstream reach. In autumn, the spatial change in diversity from one sampling station to the next along the Yangtze exhibits peaks at stations 15 and 16 (see Additional file 13: Figure S10) is most likely caused by inflow from Dongting Lake which is contaminated by Cyanobacteria. High levels of SAR11were found in certain water samples; SAR11 is well known to occur in the oceans [41, 42], and more recently, as a diverged freshwater clade in lake and river [43, 44].
The significant change in water temperature of the Yangtze from spring (average 11 °C) to autumn (average 21 °C) appears to have strongly influenced the composition of the indigenous bacterioplankton communities. It would be expected that environmental gradients in temperature and dissolved oxygen should affect the spatial or temporal variability of bacterial assemblages (Additional files 11 and 12). Even so, sedimentary microbial communities exhibited little difference between spring and autumn. Gibbons et al.  showed an insignificant correlation between seasons and community structure in a study of sedimentary bacterial communities in the Tongue River (Montana, USA). Although the connectivity of the Yangtze River networks is reflected in the distribution of the planktonic bacterial communities, the significant seasonal fluctuation in flow discharge causes the planktonic bacterial community to vary due to the introduction of exotic species originating from upstream freshwater sources and terrestrial species through the different hydrologic processes that occur in spring and autumn. It should also be noted that the different precipitation patterns during the wet and dry seasons cause more water-induced sediment erosion during autumn than in spring, and this directly affects the bacterial communities. Non-significant seasonal changes of bacterial communities were observed in sediment, which eventually reaches equilibrium state through long-term processes of sediment erosion and deposition.
Distance-decay relationships for spatial bacterial community similarity have been demonstrated previously for lakes , soil fungi , reservoirs , and soil protistan communities . In the Yangtze River, the mean dendritic distance is a parameter that indicates spatial bacterial distribution along the mainstream and its tributaries, whereas the cumulative dendritic length and site catchment area merely relate to the mainstream bacterial community assembly (Fig. 3, Additional file 7: Table S2). It appears that the transformation of ecological processes, such as dispersal and species sorting, follows the water travel time rather than the river network density  and plays a significant role in bacterial community assembly. In addition, bacterial communities originating in tributaries may modify mainstream bacterial communities . Our distance decay analysis revealed that the spatial distribution patterns of bacterial communities in water and sediment were quite different. One possible reason was that the stochastic process of turbulent flow tended to generate homogeneous environments for planktonic bacteria along the Yangtze River. However, sedimentary bacteria exhibited much lower dispersal capability than planktonic bacteria. Moreover, the bacterial communities in water and sediment responded differently to the changing seasons, generating different spatial distance-decay patterns [36, 49].
Soil type is a key factor in microbial community assembly but is not of overwhelming importance in our study. Following Liu et al. , we categorize land use in the study area in 2015 into three types: woodland-grassland, cropland-grassland, and cropland. Noting that it is difficult to distinguish the impact of soil type on microbial community structure, we follow Wang et al.  who found that land surface topography was the dominant factor influencing soil property variation (due to the effect of topography on runoff, drainage, microclimate, soil erosion, and, consequently, on soil formation). Landform type is a comprehensive parameter that describes the local soil type, land use, and soil nutrients. It was therefore deemed reasonable to explain the microbial community assembly according to landform types across the sample sites (Fig. 4, Additional file 8: Figure S6). To interpret the beta-diversity of bacterial assemblages by NMDS and ANOSIM, the Yangtze River can be divided into five reaches, according to the surrounding landform. Reach 1 (between stations 1 and 2) is situated upstream in the Hengduan Mountains where altitude is > 500 m. Reach 2 (between stations 3 and 5) is situated in the South Sichuan Foothill at the verge of the Sichuan Basin with altitude < 500 m. Reach 3 (from stations 6 to 9) is situated in the Sichuan Basin and receives inflow from the Min River. Reach 4 (from stations 9 to 13) lies in the Three Gorges where the landform comprises foothills and mountains. Reach 5 from Shashi to the estuary is located in the Middle-Lower Yangtze Plain (from stations 13 to 24). The ANOSIM test exhibited the lowest dissimilarity values between reaches 2 and 4 (Additional file 15: Table S5, and Additional file 16: Table S6) due to the matching landform types. The dissimilarity value within reach 1 was relatively higher because of the differences in altitude ranging from middle mountain (station 1, > 1000 m) to low mountain (station 2, < 1000 m). The different spatial dissimilarity values correspond to different sub-landform types. Hence, it may be concluded that the landform (a crucial factor in determining soil properties due to its effect on runoff, drainage, and soil erosion ) strongly affected spatial bacterial structures in the Yangtze River system.
In slowly varying systems such as soils  and lakes , environmental factors can be used to interpret differences between bacterial community assemblages among the sampling sites. In the Yangtze River which is a large dynamic river system, most bacterial species were found through canonical correspondence analysis to be primarily influenced by landform and large-scale human activity but less influenced by local environmental parameters (Additional file 12: Figure S9). In fact, as the landform changes along the Yangtze River so do the altitude and temperature.
The influence of large dams on bacterial structure was observed to occur mainly in the sediment samples. Additional file 8: Figure S6 (c) and (d) showed that sediment samples can be divided into two groups immediately upstream and downstream of the Three Gorges Dam and further subdivided according to spring and autumn. Construction of huge dams along the Yangtze River has caused reduced sediment load, drastic sediment erosion, and particle coarsening . As a result, the dams have significantly altered the downstream river hydrology and aquatic ecology , and the local environment physically, chemically, and biologically. Figure 5 indicates that abundance of OTUs is significantly lower immediately downstream of the dam, implying that the changes to habitats and nutrients caused by serious sediment scouring have made the downstream sediments unfavorable as a habitat for the microorganism to thrive. The higher richness of bacterial community in sediment could be expected upstream of dams due to accumulative sediment deposition, which might provide more suitable habitat for a bacterial community of higher diversity.
The neutral model provides close fits to the bacterial community structures in water and sediment, implying that the assembly of bacterial communities in the Yangtze River is close to a neutral process. The distance-decay relationship indicates that the bacterial community experiences more limited dispersal in sediment than water (Fig. 3). Moreover, Battin et al.  have shown that residence times at the streambed are longer than in the water column, and so microbes in the stream bed are expected to be more responsive to the local environment than those in the water column. It is likely that deviations between the neutral model and measured data resulted from differences in the dispersal and species sorting processes of the taxa . As a result, sedimentary taxa of higher deviations were not randomly distributed, but instead limited by dispersal processes and selected by an environment dominated by long-term sediment erosion and deposition (Fig. 6).
Given that it is difficult to account for all key environmental factors and that many environmental factors vary together with dispersal, it is almost impossible to quantify the relative contribution of species sorting . However, the species sorting process may have a relatively more powerful impact on bacterial community assembly in sediment than in water. This may be for two reasons: first, the neutral model fitted sedimentary taxa less well than planktonic taxa; and second, the heterogeneity of nutrient conditions at different sediment sites. Therefore, in the Yangtze River, the bacterial community in water was mostly shaped by the neutral processes of random dispersal because of the relatively continuous flow environment. Nevertheless, species sorting may prevail in the sedimentary community, limited by dispersal and perturbation process but led by competition and selection processes.
In general, it is possible that bacterial species in the Yangtze River randomly dispersed from their source meta-communities to the different local conditions along the length of the river. Restricted by differences in available resources, specific taxa were selected as competitive winners according to the local environment. Our results suggest that the interplay of neutral and species sorting can be used to interpret observed biogeography patterns, but their relative contributions might be spatially dependent on geographical location, local environmental changes, and anthropogenic impacts such as man-made dams.
Finally, our study demonstrates that bacterial communities in a single medium (water or sediment) do not necessarily provide a complete perspective of biogeographic patterns of bacterial community composition. This implies that a holistic understanding of both planktonic and sedimentary bacterial communities is of the utmost importance in any large river ecosystem subject to varying natural and anthropogenic impacts.
Our study revealed that a complete understanding of the integrity of bacterial geographic patterns in an entire riverine ecosystem requires knowledge of both planktonic and sedimentary communities. Based on an integrated investigation of the first spatiotemporal profiles of bacterial communities along a 4300 km continuum of the Yangtze River, we found that the sediment subgroup maintained 98.8% of the total OTUs, demonstrating that bacteria in sediment provide the greatest contribution by far to the microbial diversity of the river. Persistent OTUs responded uniformly to spring or autumn in most phyla, with autumn-associated OTUs dominant. In water, bacterial communities were most sensitive to season change. The spatial patterns of bacterial co-occurrence in both water and sediment were primarily controlled by landform, categorized here as mountain, foothill, basin, foothill-mountain, and plain. Moreover, a significant difference between sedimentary bacterial communities was found either side of the Three Gorges Dam and Xiluodu Dam. The diversity of bacterial OTUs in sediment was higher upstream of the dams than downstream due to cumulative sediment deposition, and certain OTUs presented a similar significant decline in sediment downstream of the dams largely due to severe sediment scouring. Overall, the biogeographic patterns of bacterial communities in water and sediment of the Yangtze River were satisfactorily described using meta-community theory, which considers a combination of neutral and species sorting perspectives.
The Yangtze River, the largest river in China and the third longest river in the world, originates in the Qinghai-Tibetan Plateau and flows eastward to the East China Sea. With a catchment area of 1,800,000 km2, the river drainage basin is characterized by a subtropical, warm, and wet climate. To investigate seasonal and spatial distributions of bacterial communities of Yangtze River, we carried out sampling campaigns in spring (March 2014) and autumn (October 2014). No extreme weather occurred during our sampling period. Paired water and surface sediment samples were measured almost synchronously (i.e., within 1 week) at 50 national monitoring stations along the mainstream and six major tributaries of the Yangtze. The sampling points covered a total river length of 4300 km from Shigu in the upper reach to Xuliujing near the river mouth (Additional file 1: Figure S1). Nevertheless, samples could not be taken at several monitoring sites due to their steep terrain and rapid flow conditions. Moreover, parallel sampling (1~4 samples) was undertaken at certain sites where heterogeneous sediment occurred. Noting the high velocity and turbulent flow in most reaches of the Yangtze River, parallel samples were taken at six monitoring sites. Further details of the sampling sites are listed in Additional file 3: Table S1 and Additional file 17: Table S7.
For water sampling at each site, 10-L water was collected two 5-L sterile PET bottles and immediately transported to an adjacent laboratory at low temperature 0~4 °C. Then, all water was filtered through 0.22-μm polycarbonate membranes (Millipore, USA) to capture microbial cells within 24 h. The filtered membranes were kept frozen at − 80 °C until DNA extraction. Surface sediment samples were collected where flow depth was about 0.5 m and sealed in 50-ml sterilized polypropylene tubes, placed in dry ice, and immediately transported to the laboratory. All sediment samples were stored at − 80 °C until DNA extraction.
DNA extraction, PCR amplification, and sequencing
Genomic DNA was extracted in duplicate using the FastDNA® SPIN Kit for Soil (MP Biomedicals, USA) according to the manufacturer’s protocols. The duplicate DNA extracts were finally mixed together for the following PCR amplification. The V4-V5 region of the bacterial 16S ribosomal RNA gene was amplified by PCR (95 °C for 2 min, followed by 25 cycles at 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, and a final extension at 72 °C for 5 min) using barcoded primers 515F (5′-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′) , where barcode is an eight-base sequence unique to each sample. PCR reactions were performed in triplicate 20 μL mixture containing 4 μL of 5 × FastPfu Buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.4 μL of FastPfu polymerase, and 10 ng of template DNA.
Amplicons were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) according to the manufacturer’s instructions and quantified using QuantiFluor™-ST (Promega, USA). Purified amplicons were pooled in equimolar amounts and sequenced using the strategies of PE250 (paired-end sequenced 250 × 2) on an Illumina MiSeq platform (Majorbio Company in Shanghai).
In addition, in order to monitor any contamination during the molecular workflow, negative filtration, DNA extraction, and PCR controls were applied to six negative control samples. After this step, no quantifiable DNA was detected in these negative controls, therefore, we did not analyze them further.
Sequences of bacterial 16S rRNA gene amplicons were quality-filtered using the QIIME v1.8  with the following criteria: (i) the 300 bp reads were truncated at any site receiving an average quality score < 20 over a 50-bp sliding window, discarding the truncated reads that were shorter than 50 bp; (ii) exact barcode matching, two nucleotide mismatch in primer matching, reads containing ambiguous characters were removed; (iii) pair-ended sequences that overlap longer than 10 bp were merged into a single sequence. Operational taxonomic units (OTUs) were clustered with 97% similarity cutoff using UPARSE (version 7.1), and chimeric sequences were identified and removed using UCHIME . The taxonomy of each 16S rRNA gene sequence was analyzed by RDP Classifier  against the Greengenes 16S rRNA database .
Alpha diversity was calculated using MOTHUR . Relative abundance was calculated for the OTUs in each sample, and then pairwise similarities among samples were computed using the unweighted UniFrac metrics . Unweighted UniFrac distance matrices were visualized using non-metric multidimensional scaling (NMDS). Analysis of similarity (ANOSIM) statistics was calculated to test the significance of differences among a priori sampling groups based on environmental parameters. Similarity matrices, NMDS, and ANOSIM statistics were carried out using R-3.2 with vegan package .
P < 0.05 (adjusted by false discovery rate) was considered as significant for all statistical tests unless indicated otherwise. OTUs with abundance lower than 0.01% were excluded for all statistics analysis. Singletons and doubletons, i.e., OTUs represented by one or two sequences, were not considered.
The whole datasets were split into two subgroups as water and sediment samples, which were individually analyzed. Various network methods were performed to visualize the results. The OTU distribution patterns in water and sediment samples were displayed across the taxonomic tree by directed networks using the prefuse layout algorithm using CYTOSCAPE 3.2.1 . OTUs with occupancy less than 30% were not included to avoid inconsistent trend caused by transient OTUs. An individual bipartite network was generated for each phylum to visualize the correlations among spring or autumn associated OTUs (P < 0.01) and OTUs shared by both seasons. Spearman’s rank correlation coefficient was calculated between each pair and was represented by edge length using edge-weighted spring-embedded layout algorithm in CYTOSCAPE.
The distance-decay patterns of bacterial community similarity were described by comparing three measurements determining geographic distance of site location considering the flowing direction and river network length. Site catchment area, cumulative dendritic distance, and mean dendritic distance of the sampling sites were calculated using the ArcGIS V10.2 software. Site catchment area reflected the drainage area upstream of a sampling site to that sampling site. Cumulative dendritic distance was calculated as the sum of the all paths upstream of a sampling site to that sampling site, reflecting the river network density. Mean dendritic distance represented the average length of these paths, indicative of the water residence times. Mantel tests were carried out to examine the Spearman’s rank correlation between the geographic distance matrix and the bacterial community similarity using Bray-Curtis distance matrices with 999 permutations, using the vegan package in R-3.2.
Canonical correspondence analysis (CCA) was performed to determine the effects of monitored environmental variables (temperature (T), total suspended solids (SS), dissolved oxygen (DO), electricity conductivity (EC), pH, chemical oxygen demand (COD), total nitrogen (TN), dissolved organic carbon (DOC), nitrate nitrogen (NO3-N), total phosphorus (TP), ammonium nitrogen (NH4-N), longitude (Long), latitude (Lat), altitude (A), and mean dendritic distance (Mdd)). The function “envfit” was run with 999 permutations to select the significant variables (P < 0.05). Then the significance testing was assessed by “permutest” function based on 999 permutations in R-3.2.
To estimate the potential role of neutral processes in shaping microbial community structure, the Sloan et al.’s  neutral model was fitting to describe the relationship between the observed occurrence frequency of OTUs (the proportion of local communities in which each OTU is detected) and their abundance (the mean relative abundance across all local communities) . The model is an adaptation of Hubbell’s neutral community model adjusted to bacterial populations analyzed with molecular tools . This model emphasizes the effects of stochastic dispersal and drift (birth-death immigration process) but ignores the ecological difference between species and their response to the surrounding environment. In this model, the random loss of an individual is immediately replaced by immigration from the meta-community, with probability m, or reproduction within the local community, with probability 1-m . The immigration rate was determined using non-linear least squares fitting in minpack.lm package of R-3.2. To further assess the deviations from the neutral model fitting, OTUs were subsequently sorted into three partitions depending on whether they occurred more frequently than “above” partition, less frequently than “below” partition, or within “neutral” partition the 95% confidence interval of the neutral model predictions.
Season-associated taxa analysis
OTUs with occurrence in more than 30% of all sediment or water samples were defined as persistent bacterial OTUs. An occupancy criterion was employed in order to generate the overall trend for taxonomic dendrograms (Fig. 1a, c) in water and sediment sample. Persistent bacterial OTUs that differed significantly between spring and autumn (P < 0.05) were further characterized as season-associated OTUs. According to their abundance in two seasons, season-associated OTUs were classified as autumn-associated OTUs (with significantly higher abundance in autumn) or spring-associated OTUs (with significantly higher abundance in spring). The network density value (d) was determined as the number of significant co-correlations divided by the number of all nodes, that is, a higher value represents a more intensive or dense response. Only persistent bacterial OTUs were displayed in taxonomic dendrograms (Fig. 1a, c) to avoid unstable associations and inconsistent trend caused by transient OTUs. To obtain comprehensive correlation among season-associated OTUs and non-associated OTUs, association networks were applied to both persistent and transient OTUs for each dominant phylum using edge-weighted spring-embedded layout algorithm (Fig. 1b, d).
We acknowledge the support from the Majorbio Company (Shanghai, China). We also thank Dr. Alistair Borthwick for his contribution to language editing.
This work was supported by the National Natural Science Foundation of China (Grant No.91647211 and 51721006).
Availability of data and materials
The datasets supporting the conclusions of this article were deposited into the NCBI’s Sequence Read Archive (SRA) database (Accession Number: SRP079372).
TL, ANZ, and JWW analyzed the data and prepared the draft of the manuscript. SL, XJ, CD, TM, STL, QC, and SX assisted in the data analysis and provided important suggestions for the manuscript preparation. JRN designed the research. TZ and JRN guided the data analysis and revised the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Aufdenkampe AK, Mayorga E, Raymond PA, Melack JM, Doney SC, Alin SR, et al. Riverine coupling of biogeochemical cycles between land, oceans, and atmosphere. Front Ecol Environ. 2011;9:53–60.View ArticleGoogle Scholar
- Staley C, Unno T, Gould TJ, Jarvis B, Phillips J, Cotner JB, et al. Application of Illumina next-generation sequencing to characterize the bacterial community of the Upper Mississippi River. J Appl Microbiol. 2013;115:1147–58.View ArticlePubMedGoogle Scholar
- Savio D, Sinclair L, Ijaz UZ, Parajka J, Reischer GH, Stadler P, et al. Bacterial diversity along a 2600 km river continuum. Environ Microbiol. 2015;12:4994–5007.View ArticleGoogle Scholar
- Kolmakova OV, Gladyshev MI, Rozanov AS, Peltek SE, Trusova MY. Spatial biodiversity of bacteria along the largest Arctic river determined by next-generation sequencing. FEMS Microbiol Ecol. 2014;89:442–50.View ArticlePubMedGoogle Scholar
- Read DS, Gweon HS, Bowes MJ, Newbold LK, Field D, Bailey MJ, et al. Catchment-scale biogeography of riverine bacterioplankton. ISME J. 2015;9:516–26.View ArticlePubMedGoogle Scholar
- Battin TJ, Luyssaert S, Kaplan LA, Aufdenkampe AK, Richter A, Tranvik LJ. The boundless carbon cycle. Nat Geosci. 2009;2:598–600.View ArticleGoogle Scholar
- Findlay S. Stream microbial ecology. J North Am Benthol Soc. 2010;29:170–81.View ArticleGoogle Scholar
- Madsen EL. Microorganisms and their roles in fundamental biogeochemical cycles. Curr Opin Biotechnol. 2011;22:456–64.View ArticlePubMedGoogle Scholar
- Hauptmann AL, Markussen TN, Stibal M, Olsen NS, Elberling B, Bælum J, et al. Upstream freshwater and terrestrial sources are differentially reflected in the bacterial community structure along a small Arctic River and its estuary. Front Microbiol. 2016;7:1474.View ArticlePubMedPubMed CentralGoogle Scholar
- Ruiz-González C, Niño-García JP, del Giorgio PA. Terrestrial origin of bacterial communities in complex boreal freshwater networks. Ecol Lett. 2015;18:1198–206.View ArticleGoogle Scholar
- Qian N, Zhang R, De Zhou Z. The fluvial process. Beijing: Science Press (in Chinese); 1987.Google Scholar
- Raymond PA, Hartmann J, Lauerwald R, Sobek S, McDonald C, Hoover M, et al. Global carbon dioxide emissions from inland waters. Nature. 2013;503:355–9.View ArticlePubMedGoogle Scholar
- Ward ND, Keil RG, Medeiros PM, Brito DC, Cunha AC, Dittmar T, et al. Degradation of terrestrially derived macromolecules in the Amazon River. Nat Geosci. 2013;6:530–3.View ArticleGoogle Scholar
- Borges AV, Darchambeau F, Teodoru CR, Marwick TR, Tamooh F, Geeraert N, et al. Globally significant greenhouse-gas emissions from African inland waters. Nat Geosci. 2015;8:637–42.View ArticleGoogle Scholar
- Feng BW, Li XR, Wang JH, Hu ZY, Meng H, Xiang LY, et al. Bacterial diversity of water and sediment in the Changjiang estuary and coastal area of the East China Sea. FEMS Microbiol Ecol. 2009;70:80–92.View ArticlePubMedGoogle Scholar
- Xia N, Xia X, Liu T, Hu L, Zhu B, Zhang X, et al. Characteristics of bacterial community in the water and surface sediment of the Yellow River, China, the largest turbid river in the world. J Soils Sediments. 2014;14:1894–904.View ArticleGoogle Scholar
- Crump BC, Peterson BJ, Raymond PA, RMW A, Rinehart A, JW MC, et al. Circumpolar synchrony in big river bacterioplankton. Proc Natl Acad Sci U S A. 2009;106:21208–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Nelson CE. Phenology of high-elevation pelagic bacteria: the roles of meteorologic variability, catchment inputs and thermal stratification in structuring communities. ISME J. 2009;3:13–30.View ArticlePubMedGoogle Scholar
- Fortunato CS, Herfort L, Zuber P, Baptista AM, Crump BC. Spatial variability overwhelms seasonal patterns in bacterioplankton communities across a river to ocean gradient. ISME J. 2012;6:554–63.View ArticlePubMedGoogle Scholar
- Cavender-Bares J, Kozak KH, Fine PVA, Kembel SW. The merging of community ecology and phylogenetic biology. Ecol Lett. 2009;12:693–715.View ArticlePubMedGoogle Scholar
- Leibold MA, Holyoak M, Mouquet N, Amarasekare P, Chase JM, Hoopes MF, et al. The metacommunity concept: a framework for multi-scale community ecology. Ecol Lett. 2004;7:601–13.View ArticleGoogle Scholar
- Logue JB, Lindström ES. Biogeography of bacterioplankton in inland waters. Fr Rev. 2008;1:99–114.Google Scholar
- Vanwonterghem I, Jensen PD, Dennis PG, Hugenholtz P, Rabaey K, Tyson GW. Deterministic processes guide long-term synchronised population dynamics in replicate anaerobic digesters. ISME J. 2014;8:2015–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Hubbell SP. The unified neutral theory of biodiversity and biogeography. Princeton: Princeton University Press; 2001.Google Scholar
- Battin TJ, Besemer K, Bengtsson MM, Romani AM, Packmann AI. The ecology and biogeochemistry of stream biofilms. Nat Rev Microbiol. 2016;14:251–63.View ArticlePubMedGoogle Scholar
- Logue JB, Lindström ES. Species sorting affects bacterioplankton community composition as determined by 16S rDNA and 16S rRNA fingerprints. ISME J. 2010;4:729–38.View ArticlePubMedGoogle Scholar
- Dumbrell AJ, Nelson M, Helgason T, Dytham C, Fitter AH. Relative roles of niche and neutral processes in structuring a soil microbial community. ISME J. 2010;4:337–45.View ArticlePubMedGoogle Scholar
- Niño-García JP, Ruiz-González C, del Giorgio PA. Interactions between hydrology and water chemistry shape bacterioplankton biogeography across boreal freshwater networks. ISME J. 2016;10:1755–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Langenheder S, Székely AJ. Species sorting and neutral processes are both important during the initial assembly of bacterial communities. ISME J. 2011;5:1086–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Van der Gucht K, Cottenie K, Muylaert K, Vloemans N, Cousin S, Declerck S, et al. The power of species sorting: local factors drive bacterial community composition over a wide range of spatial scales. Proc Natl Acad Sci. 2007;104:20404–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu H, Wang Z, Lu Y. Self-adjustment mechanism of bed structures under hydrology and sediment regimes. Hydrol Res. 2016;47:136–48.View ArticleGoogle Scholar
- Wang Z-Y, Li Y, He Y. Sediment budget of the Yangtze River. Water Resour Res. 2007;43:W04401.Google Scholar
- Zhang W, Yuan J, Han J, Huang C, Li M. Impact of the Three Gorges Dam on sediment deposition and erosion in the middle Yangtze River: a case study of the Shashi Reach. Hydrol Res. 2016;47:175–86.View ArticleGoogle Scholar
- Sloan WT, Lunn M, Woodcock S, Head IM, Nee S, Curtis TP. Quantifying the roles of immigration and chance in shaping prokaryote community structure. Environ Microbiol. 2006;8:732–40.View ArticlePubMedGoogle Scholar
- Logares R, Lindström ES, Langenheder S, Logue JB, Paterson H, Laybourn-Parry J, et al. Biogeography of bacterial communities exposed to progressive long-term environmental change. ISME J. 2013;7:937–48.View ArticlePubMedGoogle Scholar
- Burns AR, Stephens WZ, Stagaman K, Wong S, Rawls JF, Guillemin K, et al. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 2016;10:655–64.View ArticlePubMedGoogle Scholar
- Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, et al. Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc Natl Acad Sci U S A. 2006;103:12115–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Galand PE, Casamayora EO, Kirchmand DL, Lovejoye C. Ecology of the rare microbial biosphere of the Arctic Ocean. Proc Natl Acad Sci U S A. 2009;106:22427–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersson AF, Riemann L, Bertilsson S. Pyrosequencing reveals contrasting seasonal dynamics of taxa within Baltic Sea bacterioplankton communities. ISME J. 2010;4:171–81.View ArticlePubMedGoogle Scholar
- Biers EJ, Sun S, Howard EC. Prokaryotic genomes and diversity in surface ocean waters: interrogating the global ocean sampling metagenome. Appl Environ Microbiol. 2009;75:2221–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Rusch DB, Halpern AL, Sutton G, Heidelberg KB, Williamson S, Yooseph S, et al. The Sorcerer II Global Ocean Sampling expedition: Northwest Atlantic through eastern tropical Pacific. PLoS Biol. 2007;5:e77.View ArticlePubMedPubMed CentralGoogle Scholar
- Zwart G, Crump BC, Agterveld MPK, Hagen F, Han SK. Typical freshwater bacteria: an analysis of available 16S rRNA gene sequences from plankton of lakes and rivers. Aquat Microb Ecol. 2002;28:141–55.View ArticleGoogle Scholar
- Bahr M, Hobbie JE, Sogin ML. Bacterial diversity in an arctic lake: a freshwater SAR11 cluster. Aquat Microb Ecol. 1996;11:271–7.View ArticleGoogle Scholar
- Gibbons SM, Jones E, Bearquiver A, Blackwolf F, Roundstone W, Scott N, et al. Human and environmental impacts on river sediment microbial communities. PLoS One. 2014;9:e97435.View ArticlePubMedPubMed CentralGoogle Scholar
- Green JL, Holmes AJ, Westoby M, Oliver I, David B, Dangerfield M, et al. Spatial scaling of microbial eukaryote diversity. Nature. 2004;432:747–50.View ArticlePubMedGoogle Scholar
- Liu L, Yang J, Yu Z, Wilkinson DM. The biogeography of abundant and rare bacterioplankton in lakes and reservoirs of China. ISME J. 2015;9:2068–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Bates ST, Clemente JC, Flores GE, Walters WA, Parfrey LW, Knight R, et al. Global biogeography of highly diverse protistan communities in soil. ISME J. 2013;7:652–9.View ArticlePubMedGoogle Scholar
- Östman Ö, Drakare S, Kritzberg ES, Langenheder S, Logue JB, Lindström ES. Regional invariance among microbial communities. Ecol Lett. 2010;13:118–27.View ArticlePubMedGoogle Scholar
- Liu J, Kuang W, Zhang Z, Xu X, Qin Y, Ning J, et al. Spatiotemporal characteristics, patterns and causes of land use changes in China since the late 1980s. J Geogr Sci. 2014;24:195–210.View ArticleGoogle Scholar
- Wang J, Fu B, Qiu Y, Chen L. Soil nutrients in relation to land use and landscape position in the semi-arid small catchment on the loess plateau in China. J Arid Environ. 2001;48:537–50.View ArticleGoogle Scholar
- Reith F, Brugger J, Zammit CM, Gregg AL, Goldfarb KC, Andersen GL, et al. Influence of geogenic factors on microbial communities in metallogenic Australian soils. ISME J. 2012;6:2107–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Logue JB, Langenheder S, Andersson AF, Bertilsson S, Drakare S, Lanzén A, et al. Freshwater bacterioplankton richness in oligotrophic lakes depends on nutrient availability rather than on species-area relationships. ISME J. 2012;6:1127–36.View ArticlePubMedGoogle Scholar
- Dai Z, Liu JT. Impacts of large dams on downstream fluvial sedimentation: an example of the Three Gorges Dam (TGD) on the Changjiang (Yangtze River). J Hydrol. 2013;480:10–8.View ArticleGoogle Scholar
- Janzen T, Haegeman B, Etienne RS. A sampling formula for ecological communities with multiple dispersal syndromes. J Theor Biol. 2015;374:94–106.View ArticlePubMedGoogle Scholar
- Crump BC, Adams HE, Hobbie JE, Kling GW. Biogeography of bacterioplankton in lakes and streams of an arctic tundra catchment. Ecology. 2007;88:1365–78.View ArticlePubMedGoogle Scholar
- Shan J, Ji R, Yu Y, Xie Z, Yan X. Biochar, activated carbon, and carbon nanotubes have different effects on fate of 14C-catechol and microbial community in soil. Sci Rep. 2015;5:16000.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.View ArticlePubMedGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Env Microbiol. 2007;73:5261–7.View ArticleGoogle Scholar
- McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.View ArticlePubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. Vegan: community ecology package. R package version 2.4–0. 2016. https://cran.r-project.org/package=vegan. Accessed 15 Aug 2016.
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Sloan WT, Woodcock S, Lunn M, Head IM, Curtis TP. Modeling taxa-abundance distributions in microbial communities using environmental sequence data. Microb Ecol. 2007;53:443–55.View ArticlePubMedGoogle Scholar