Skip to main content

Molecular biogeography of planktonic and benthic diatoms in the Yangtze River

Abstract

Background

Diatoms are of great significance to primary productivity in oceans, yet little is known about their biogeographic distribution in oligotrophic rivers.

Results

With the help of metabarcoding analysis of 279 samples from the Yangtze River, we provided the first integral biogeographic pattern of planktonic and benthic diatoms over a 6030 km continuum along the world’s third largest river. Our study revealed spatial dissimilarity of diatoms under varying landforms, including plateau, mountain, foothill, basin, foothill-mountain, and plain regions, from the river source to the estuary. Environmental drivers of diatom communities were interpreted in terms of photosynthetically active radiation, temperature, channel slope and nutrients, and human interference. Typical benthic diatoms, such as Pinnularia, Paralia, and Aulacoseira, experienced considerable reduction in relative abundance downstream of the Three Gorges Dam and the Xiluodu Dam, two of the world’s largest dams.

Conclusions

Our study revealed that benthic diatoms are of particular significance in characterizing motile guild in riverine environments, which provides insights into diatom biogeography and biogeochemical cycles in large river ecosystems.

Background

Diatoms play a particularly important role in the biogeochemical cycle [1] of primary elements such as carbon, nitrogen, phosphorus, and silica, contributing about 20–25% of global primary production [2]. Diatoms are ubiquitous and diverse species of single-celled, eukaryotic, photosynthetic microorganisms [3], and are often the dominant primary producers in marine and freshwater ecosystems [4]. Therefore, diatoms in such ecosystems may be remarkably dissimilar either in phylogenetic composition or biogeographic distribution [5, 6]. Freshwater bodies typically consist of lentic (particularly lakes and wetlands) and lotic waters (including streams and rivers), which are often dominated respectively by planktonic algae and benthic species [7].

Accurate identification of diatoms depends on the reliability of the methods used. Morphological analysis requires extensive taxonomic expertise and may exhibit shortcomings in characterizing specific diatoms in rivers [8]. With the development of high-throughput sequencing (HTS) technology, DNA metabarcoding has become a rapid, accurate, and reliable method for diatom detection [9]. Various DNA barcoding studies have been successfully conducted, based on different maker genes, including COI [10], ITS [11], and 18S rDNA [9, 12]. Malviya et al. [13] provided a new estimate of diversity and distribution of marine planktonic diatoms based on the V9 region of eukaryotic 18S rDNA. As a result, the most widespread and diverse diatom genera are derived from 46 marine stations. Recently, the V4 region of 18S rDNA was proposed for diatom barcoding in studies of diatoms in river and deltaic systems [9, 12].

Comparing the numerous studies of diatoms and eutrophication in oceans [14, 15] and lakes [16, 17] to date, it is clear that the present understanding of diatoms is relatively poor for lotic and oligotrophic rivers [7]. In fact, previous reports on the dynamics of riverine diatoms have mostly focused on tributaries, small rivers, reaches, stations, and estuaries [18,19,20]. Many studies examined the diversity and composition of planktonic [19, 20] or benthic diatoms [18, 21, 22] based on morphological identification. For example, Centis et al. [20] investigated planktonic diatoms dominated by physical constraints at two stations of the River Adige, Italy. Liu et al. [18] investigated the community structure of benthic diatoms in the Dong River, one of the three main tributaries of the Pearl River, China. Although Kireta et al. [23] observed that both planktonic and periphytic diatoms could be used as bio-indicators of river conditions, little is known about the distinction between planktonic and benthic diatoms regarding their spatiotemporal distributions.

Biogeography studies aim to reveal the spatial and temporal distribution of biodiversity and provide insight into the mechanisms that generate and sustain diversity [24]. Spatial dispersal and environmental selection processes are regarded as essential drivers for the biogeographical pattern of bacterial community [25]. The former promotes movement of species and their establishment at a new location, whereas the latter alters the abundance and composition of species, according to the ability to survive and reproduce under local environmental conditions. A similar explanation has been proved to apply to the biogeographical pattern of planktonic or benthic diatom communities in small rivers using morphological analysis [26,27,28]. However, it remains unclear how the integrated spatiotemporal distributions of planktonic and benthic diatom communities are shaped by spatial dispersal and environmental selection processes in large rivers subject to complex natural and anthropogenic impacts.

To close the above gap, we implemented large-scaled synchronous monitoring of diatom communities at 62 hydrologic stations over a 6030 km continuum of the Yangtze River in China. Consequently, we provided the first molecular biogeographic pattern of both planktonic and benthic diatoms in the largest river in Asia (Fig. 1). Meanwhile, environmental drivers of diatom communities were interpreted in terms of photosynthetic active radiation, temperature, channel slope, and nutrient conditions under varying landforms.

Fig. 1
figure1

Flowchart of the study. a Two hundred seventy-nine water and sediment samples at 62 hydrologic stations in the Yangtze River covering the actual sinuous channel reach of length 6030 km (equivalent to 1.83 times the 3290 km straight line from start to end sampling sites). b Metabarcoding analysis provides insights into biogeographic pattern of diatoms along the mainstream of the Yangtze River, represented by the spatial distribution of Shannon diversity. c Interpretations on the biogeographic patterns of diatom communities, with main influencing factors such as photosynthetically active radiation (PAR), channel slope, and nutrients characterized by ratio of total nitrogen to total phosphate (TN:TP)

Results

Our study generated a total of 8,602,620 V4 18S rDNA reads from 279 samples. All sequencing reads were classified into 3947 operational taxonomic units (OTUs) at a 97% similarity threshold, with 3144 OTUs well matching 454 diatom species in our reference database. Rarefaction curves (Additional file 1: Figure S1) together with high values of Good’s coverage ranging from 0.9854 to 0.9992 illustrated that OTUs obtained by the current sequencing depth gave a reasonable representation of the diatom communities. The phylogeny tree, constructed by representative OTUs (accounting for > 90% sequence in all samples) and reference sequences (Additional file 1: Figure S2), further confirmed the accuracy of taxonomic assignment.

Alpha and beta diversity of diatom communities

Molecular barcoding based on high-throughput sequencing (HTS) provided a detailed diatom directory for the whole Yangtze River at different taxonomy levels, i.e., 4 classes, 37 orders, 60 families, and 152 genera.

HTS is of particular use in detecting nano-sized diatoms (2–20 μm) in the Yangtze River, confirming the presence of Fragilaria perminuta, Achnanthidium minutissimum, Achnanthidium saprophilum, Amphora pediculus, Fistulifera saprophila, Mayamaea permitis, Sellaphora seminulum, Encyonema minutum, Fragilaria famelica, Fragilaria rumpens, Gomphonema pumilum, Staurosirella pinnata, Planothidium frequentissimum, Craticula buderi, and Craticula molestiformis.

Six types of environmental samples were taken along the Yangtze River, including water and sediment samples from the river source region (i.e., water-plateau (12 samples) and sediment-plateau (12 samples)) and those from the mainstream in the non-plateau area (i.e., water-spring (38 samples), water-autumn (46 samples), sediment-spring (87 samples), and sediment-autumn (84 samples)). Planktonic diatoms exhibited the highest alpha-diversity (Chao1 and Shannon indices) and benthic diatoms the lowest richness (Chao1) in the plateau (Additional file 1: Figure S3). In the non-plateau area, no significant differences were observed in the alpha richness and diversity of diatom communities in the four sample types.

Non-metric multidimensional scaling (NMDS) analysis of the compositional dissimilarities between diatom communities demonstrated not only a clear spatial differentiation in diatoms between the plateau and the main body of the Yangtze, but also a division between planktonic and benthic groups (Additional file 1: Figure S4). Water and sediment samples between spring and autumn in the non-plateau area were used for further seasonal analysis. Seasonal difference in planktonic diatoms was found much more significant than in benthic diatoms, as further confirmed by an analysis of similarity (ANOSIM) test (Additional file 1: Figure S5). Moreover, one-way analysis of variance (one-way ANOVA) indicated that more planktonic diatoms (42.75 ± 13.98% relative abundance, primarily belonging to Cyclotella, Stephanodiscus, and Skeletonema) than benthic diatoms (16.58 ± 5.06% relative abundance, primarily belonging to Pinnularia and Stephanodiscus) exhibited significant seasonal sensitivity (Additional file 1: Figure S6).

Biogeographic patterns of diatom communities

A variety of diatom species have been found closely relevant to carbon export [29]. In the Yangtze River, the planktonic diatoms such as Asterionella formosa, Diatoma vulgare, Lindavia viaradiosa, Gomphonema pumilum, and Thalassiosira nordenskioeldii were significantly strongly associated with dissolved carbon dioxide (pCO2, see the “Methods” section) (Spearman r > 0.3, P < 0.05), while the benthic diatoms Asterionella formosa, Encyonema prostratum, Eucocconeis laevis, Fistulifera saprophila, and Nitzschia sigmoidea were highly correlated with pCO2 (Additional file 1: Figure S7).

Obvious difference in species composition was observed in planktonic and benthic diatoms. In the Yangtze River, diatoms mainly consisted of Coscinodiscophyceae, Fragilariophyceae, Bacillariophyceae, and Mediophycea. Planktonic diatoms were dominated by Coscinodiscophyceae (about 43.76% of the total number of sequences) and Mediophyceae (17.91%), while benthic diatoms were dominated by Bacillariophyceae (54.88%) and Coscinodiscophyceae (30.96%) (Additional file 1: Figure S8). Planktonic and benthic diatoms were not always consistent in dominant genera (top 20, the relative abundance ranged from 55.6 to 83.6%) (Fig. 2a). The dominant genera were found to be Cyclotella, Stephanodiscus, Pinnularia, and Paralia, represented 12.2, 8.6, 7.3, and 6.6% of total sequences, respectively, in water samples. Meanwhile, Navicula, Pinnularia, and Cyclotella became the dominant genera, represented 14.4, 9.1, and 6.9% of total sequences, respectively, in sediment samples, in which Navicula was dominant in either sediment-plateau (17.1%), sediment-autumn (13.1%), or sediment-spring samples (13.2%).

Fig. 2
figure2

a Circular visualization of dominant diatoms at genus level in six sample types. Inner circular diagram shows relative abundance of different diatom genus in six sample types. Only the dominant genus with a mean relative abundance of ≥ 1% in all samples is depicted. The width of ribbons for each diatoms is directly proportional to their relative abundance in each sample type. Similarly, different colored ribbons of different width for each sample type describe the distribution of different genera. b Representative diatoms genera in different landform types from the river source to mouth along the Yangtze River

Over the 6030 km continuum from river source to mouth, landform played a significant role in spatial differentiation of both planktonic and benthic diatom communities. Referring to a previous study on landform types in the Yangtze [30], planktonic diatoms were represented by Cymbella, Asterionella, Stephanodiscus, Melosira, Cyclotella, and Conticribra in the plateau, mountain, foothill, basin, foothill-mountain, and plain regions, respectively, while benthic diatoms were abundant by Cymbella, Navicula, Melosira, Conticribra, Cyclotella, and Surirella, respectively, in the corresponding regions (Fig. 2b).

A completed description on biogeographic pattern over a large river requires to identify the difference in diatom compositions among different types of samples and their spatiotemporal heterogeneity. Using the indicator species analysis, the diatoms that were responsible for the observed community differences among the six types of samples could be well identified (Additional file 1: Table S1). The number of indicator diatom species in the river ranged from 6 (sediment-spring) to 41 (water-plateau). Diatom communities in the plateau region were quite different from those in the non-plateau region of the Yangtze River, as evidenced by the higher percentage of top indicator species in water-plateau and sediment-plateau samples (Additional file 1: Figure S9). The average relative abundance of indicator species in the source area exceeded 40%, and planktonic indicator species contributed more reads than benthic indicator species. Furthermore, a number of indicator species belonging to Tabellariales and Hemiaulales occurred in water-plateau and sediment-plateau samples, respectively.

Diatom composition in terms of ecological guilds showed spatial dissimilarity in water and sediment samples (Additional file 1: Figure S10). Diatoms were divided into four ecological guilds according to their biological traits, including low-profile, high-profile, motile, and planktic guilds in terms of different responses to nutrients and dynamic disturbances [31,32,33] (see “Ecological guilds classification” section). Benthic diatoms in the motile guild prevailed at most stations along the whole river, whereas those in high-profile and planktic guilds dominated upstream and downstream reaches, respectively. In addition, planktonic diatoms in the planktic guild were predominant at most stations along the Yangtze River.

Environmental effects on diatom biogeography

Significant distance-decay in diatom similarity was observed along the geographical distance (Additional file 1: Figure S11), with a greater slope of the curve for water (slope = − 0.042) than for sediment (slope = − 0.038) using least squares linear regression. The partial Mantel test demonstrated that both geographical and environmental distances played important roles in constraining diatom composition and distribution (Additional file 1: Table S2). Canonical correspondence analysis (CCA) showed significant correlations between diatom communities and specific environmental and spatial factors such as water temperature, pH, suspended solids, and PCNM-1 (Additional file 1: Table S3). Variation partitioning of diatom composition showed that a greater percentage (14.6–21.2%) could be explained by a purely environmental component than that (3.4–6.0%) of the total variation by a pure spatial component (Additional file 1: Figure S12), and a minor portion (0.4–5.4%) explained by spatially structured environmental heterogeneity, leaving the majority of the total variation (68.7–79.0%) inexplicable. As a deterministic process, environmental selection played a critical role in the biogeography of planktonic and benthic diatoms. Although environmental differentiation seemed more important than spatial dispersing in shaping a diatom community, neither could fully explain the total variation in diatom composition. Among others, the typical environmental components such as photosynthetically active radiation, temperature, channel slope, and nutrients condition are essential to diatom community accompanied with the spatial dispersal.

Photosynthetically active radiation (PAR, 400–700 nm) is utilized by diatoms to synthesize biomass through photosynthesis [34]. Spatially, the annual-averaged PAR exhibits four stages along the Yangtze River [35], i.e., the highest in the upper reach located in Qinghai-Tibet Plateau region (above 32 mol m−2 d−1), the higher in the reach located in the Hengduan Mountains, the lowest in the reach located in the Sichuan Basin (below 23 mol m d−1), and the moderate in the lower reach (Fig. 3c, see the “Photosynthetically active radiation (PAR) divisions” section). For a better understanding of the spatial heterogeneity of both planktonic and benthic diatom communities, LefSe analysis was used considering its advantages in identifying differentially abundant taxa under different environmental conditions [36]. Consequently, preferred planktonic and benthic diatom species in different PAR regions were identified (Fig. 3a, b). For example, the Caloneis, Cymbella, Fistulifera, and Fragilaria genera preferred very-high PAR zones, the Papiliocellulus genus favored medium PAR regions, and Conticribra and Cyclotella showed a preference for low PAR habitats. Planktonic Cymatopleura and Navicula, and benthic Asterionella, Biddulphia, Diatoma, and Encyonema genera preferred to high PAR conditions. Moreover, water temperature is a key environmental factor in structuring diatom community assemblages through its influence on diatom size and growth rate [37] in the Yangtze River (Additional file 1: Table S3). Although the richness of planktonic diatoms seemed to fluctuate with PAR, the richness of benthic diatoms tended to rise with increasing temperature (Fig. 3c).

Fig. 3
figure3

LEfSe cladogram of planktonic (a) and benthic (b) diatom communities from four PAR regions. Diatom taxa with a mean relative abundance of ≥ 0.1% in all samples, assigned to kingdom (innermost), phylum, class, order, family, and genus (outermost), are used to determine taxa or clades most likely to explain differences between PAR regions. Differentially abundant taxa (biomarkers) are colored by their most abundant PAR regions, i.e., red, green, blue, and purple circles stand for biomarkers in regions of very high, high, medium, and low. Orange and blue circles display the average alpha-diversity (Chao1) of planktonic and benthic diatoms respectively in different photosynthetically active radiation (PAR) regions, and their sizes correspond to the Chao1 index (c)

Stream power, often simply characterized by the river channel slope or the product of channel slope and flow discharge (except in plateau regions) [38], is another important factor altering the spatial distribution of diatoms. The channel slope dramatically changes along the Yangtze River, primarily due to the basis of geology, climate and geomorphology. In the mountainous reaches (stations 1~2), the river channel slope can be higher than 400 × 10−5. In the upper reaches (stations 3~14), the channel slope has dropped sharply to 10–30 × 10−5. In middle-lower reaches (stations 15~24), the slope of the riverbed is nearly to zero as the channel widens and shallows in the estuarine region. In general, the varying channel slope along the Yangtze River could be simplified into three stages, steep slope in mountainous reaches, moderate slope in upper reaches, and mild slope in middle-lower reaches (Fig. 4c). During the wet season (autumn), the higher flow discharge weakens the correlation between planktonic community similarity and channel slope, although a stronger correlation between benthic community similarity and channel slope is maintained due to higher mobility of the streambed (Fig. 4a, b). In view of their relative abundance, planktonic diatoms were characterized by Psammothidium, Nitzschia, and Cymbella for steep slope environments, Papiliocellulus for moderate slope, Mayamaea, Pinnularia, and Surirella for mild slope environments. Benthic diatoms were represented by Cocconeis, Entomoneis, and Melosira for steep slope environments, Fallacia, Psammothidium, and Skeletonema for moderate slopes, and Actinocyclus, Aulacoseira, and Conticribra for mild slopes (Additional file 1: Figure S13). Furthermore, slope effects on diatoms might be identified in terms of ecological guilds. Regardless of the diatoms in an unspecified ecological guild, planktonic diatoms were dominant in the planktic guild. Interestingly, species in motile guild stably constituted the main component of benthic diatoms in the whole lotic river (Fig. 4c).

Fig. 4
figure4

Relationships between community similarity and river channel slope for water-spring (a), water-autumn (b), sediment-spring (c), and sediment-autumn (d) samples. Values of Mantel Spearman correlations (r) and probabilities (P) are also provided. Gray lines denote ordinary least squares linear regression fits across all samples. Spatial distributions of ecological guilds for different channel slopes are shown in (e)

Nutrient condition, represented by the total nitrogen to total phosphorus ratio (TN:TP), was considered as an indicator of ecosystem to support for algae biomass [39]. Since diatoms were reported to predominate in phytoplankton at high nitrogen to phosphorus ratio (N:P > 16) in water [40], we investigated the response of planktonic and benthic diatoms to TN:TP in the Yangtze River. Higher TN:TP was observed in water samples (13.8~45.63) than in sediment samples (0.004~0.65). The spatial distribution of TN:TP in water and sediment samples varied in different sampling stations. The lowest value of TN:TP in water samples was obtained in river reaches receiving large inflow from the Min River, and the highest value was obtained in both water and sediment samples near river estuary. Moreover, TN:TP exhibited a fluctuation between upstream and downstream of large dams (Xiluodu Dam and Three Gorges Dam). The alpha-diversity of planktonic diatoms could be partly interpreted in terms of annual-averaged dissolved TN:TP (data range from 2005 to 2014) in water-spring (Adj R2 = 0.54) and water-autumn (Adj R2 = 0.41) samples (Additional file 1: Figure S14a-b); however, benthic diatoms demonstrated an even weaker response to TN:TP (monitored data) (Additional file 1: Figure S14c-d). Furthermore, the monitored TN:TP was significantly related to community composition of diatoms in both water and sediment (P < 0.05, Additional file 1: Table S3). For relatively oligotrophic rivers like the Yangtze River, it appeared that TN was more important than TP as limiting nutrients to planktonic and benthic community variability (Additional file 1: Table S3).

Nevertheless, the aforementioned environmental factors are subject to change due to human interference. For the Yangtze River, one of the most prominent impacts on its ecology arises from the construction and operation of large dams. In the present study, one-way ANOVA analysis revealed that certain OTUs exhibited significant variations in distribution immediately upstream and downstream of the Three Gorges Dam (P < 0.01) and Xiluodu Dam (P < 0.05). Sudden drops in relative abundance of OTUs belonging to specific benthic species (such as Pinnularia, Paralia, and Aulacoseira) occurred downstream of the dams (Additional file 1: Figure S15). Furthermore, one-way ANOVA analysis also revealed significant difference in environmental variables (pH (P = 0.003), NH4-N (P < 0.001), TN (P < 0.001), and TN:TP (P = 0.01)) for sediment samples immediately before and after the Three Gorges Dam.

Discussion

In riverine ecosystems, both planktonic and benthic algae are the important components of primary producers. Previous studies mainly focused on species composition and richness of planktonic algae, and their responses to pollutions, disturbance, as well as physic-chemical factors in large rivers [18,19,20,21,22]. In this study, we implemented the synchronous sampling of water and sediment in the Yangtze River, which enables to reveal integrative structure of planktonic and benthic diatoms, and environmental drivers for their biogeographic patterns.

Taxonomic compositions of riverine diatoms in the Yangtze River are markedly different from those in lakes [41] and oceans [13] because of the distinct differences in salinity and hydraulic conditions. For example, Actinocyclus spp., Aulacoseira spp., Cyclotella spp., Fragilaria spp., and Synedra spp. predominated in diatom communities in Lake Kasumigaura, Japan [41]. Moreover, Chaetoceros, followed by Fragilariopsis, Thalassiosira, and Corethron were reported to be the most abundant genera in the global oceans [13]. In marine ecosystems, diatoms were abundant in nutrient-rich coastal zones particularly at high latitudes. In the riverine ecosystem, diatoms of glacier origin from the Qinghai-Tibetan plateau were significantly different to those in the main body of the Yangtze River (Fig. 2, Additional file 1: Figure S4 and Figure S9). In order of abundance, the most common riverine diatoms in water and sediment of the Yangtze River were Cyclotella, Navicula, Pinnularia, Stephanodiscus, and Cymbella genera.

Synchronous water and sediment sampling along the Yangtze River makes it possible to compare the biogeography of planktonic and benthic diatoms. On the one hand, although the diversity of planktonic and benthic diatoms was very similar in non-plateau area (Additional file 1: Figure S3), their community composition was remarkably different. Planktonic diatoms were dominated by Cyclotella (average relative abundance: 17.66%), Stephanodiscus (12.81%), Pinnularia (7.53%), Paralia (7.34%), and Skeletonema (4.90%), whereas benthic diatoms were dominated by Navicula (13.12%), Cyclotella (10.33%), Pinnularia (10.12%), Surirella (7.10%), and Stephanodiscus (6.53%). On the other hand, seasonal differences were more evident in planktonic diatoms than benthic diatoms (Additional file 1: Figure S4-S6). One-way ANOVA analysis confirmed that water temperature of the Yangtze River was significantly different between spring and autumn (P < 0.01), with the water temperature (average 21 °C) in autumn being more conducive than that (average 11 °C) in spring for planktonic diatom growth [42]. Moreover, seasonal fluctuations in water discharge appear to affect the community structure of planktonic diatoms, owing to the introduction of diatom species from the upstream freshwater source and to different hydrologic processes in spring and autumn [43]. Meanwhile, the weak seasonal difference of diatom communities in sediment may be ascribed to the majority of benthic diatoms (Additional file 1: Figure S6) that weakly respond to seasonal changes and reached a state of relative equilibrium through long-term sediment erosion and deposition processes [30].

Interaction between planktonic and benthic diatoms has been one of the major concerns for large river ecosystems. In the Yangtze River, such interactions have several consequences. First, the community composition of planktonic community was significantly correlated with that of benthic diatoms in paired water and sediment samples (in spring: Spearman r = 0.3556, P = 0.001; in autumn: Spearman r = 0.1902, P = 0.006). Typical benthic diatoms (e.g., Nitzschia and Navicula) were found in high abundance in the water column. The local interactions could cause benthic and planktonic habitats to become coupled through migration of algal cells, meaning that phytoplankton can be derived from benthic diatoms, and sinking planktonic algae can become benthic algae [44]. Second, the richness of planktonic diatoms appeared to fluctuate with PAR, whereas the richness of benthic diatoms tended to change with temperature. Nutrient level (TN:TP) had different effects on the diversity and variation of planktonic and benthic diatoms. These phenomena further explain the local interactions in terms of light, temperature, and nutrient competition [45]. Third, given the “River continuum concept” [46] and “Continuous discontinuity concept” [47], the dominance of benthic or planktonic algae changes with the natural riverine gradient (e.g., channel slope), and interactions are invariably interrupted by anthropogenic disturbances such as dam construction and pollutants discharge. Damming has been shown to cause various changes in hydraulic conditions [48,49,50], such as sediment erosion, deepening river channel, coarse riverbed, and high-suspended particles, which makes the downstream sediments become a habitat with limited light, difficult to adhere, and unsuitable for benthic diatoms to thrive (Additional file 1: Figure S15). Furthermore, it not only weakens the competitive interactions for light and nutrients between benthic and planktonic diatoms, but also reduces their migration rate. Namely, as river channel deepens, the number of planktonic diatoms from detached benthic algae and benthic diatoms from sinking planktonic algae will reduce, resulting a dominant position of primary production by planktonic algae [45]. It can be further confirmed, that is, contrary to the planktonic diatoms, the average richness of benthic diatoms immediately downstream of Three Gorges Dam is significantly lower than that of upstream reaches. Other factors affecting interactions, such as river velocity, turbulent diffusion, algal sinking, and grazing remain to be further investigated [51].

The Yangtze River flows through a variety of landforms such as plateau, mountain, foothill, basin, foothill-mountain, and plain regions along the river continuum [30]. Planktonic and benthic diatoms demonstrate spatial dissimilarities in different types of landforms corresponding to varying soil property, temperature, altitude, light, and nutrients. For example, the Qinghai-Tibet Plateau is characterized with lower temperature, higher altitude, and the greatest PAR, while the basin (mostly located in Sichuan province) receives the lowest light resource. Further going down to the Middle-Lower Yangtze Plain, there exists a general decrease of altitude and increase of temperature and nutrient levels (nitrogen and phosphorus). The differentially abundant diatoms in both water and sediments in the six landform types could also be convinced by previous studies on special diatoms in different landforms such as mountain [52] and plateau [53]. Similar findings have also been reported recently on bacterial communities. For example, Liu et al. [30] confirmed that landform type is of significance to bacterial structures in the Yangtze River. Moreover, Chen et al. [54] indicated that anammox bacterial abundance and alpha diversity are spatially subject to landform variations.

It remains a practical challenge to examine the effects of environmental selection and spatial dispersal on determining algae communities in river ecosystems. In the Yangtze River, environmental selection appears to explain at least three times more of the community variance (either in planktonic or benthic diatoms) than spatial dispersal process (Additional file 1: Figure S12), which is in line with previous findings on planktonic [26] and benthic diatoms [28] in streams and rivers. Keck et al. reported the importance of both environment and dispersal-related processes in controlling planktonic diatom community structure in stream and rivers, Sweden [26]. Soininen et al. concluded that local environmental factors were even more important than spatial factors in explaining benthic diatom distributions in boreal streams, Finland [28]. Overall, environmental selection is the dominant driving force on the biogeographical pattern of planktonic and benthic diatom communities in the Yangtze River.

The environmental-based selection process driving the biogeographic pattern of diatoms is influenced by PAR, temperature, channel slope, and nutrient level. The community structure of benthic diatoms is a typical consequence driven by environmental selection in the Yangtze River. First, it has been reported that sufficient PAR drives the growth and production of diatoms [55], but excess PAR can affect various cellular processes and reduce the growth or viability of diatoms [56]. In the present study, specific diatoms were identified as preferring to different levels of light intensity (Fig. 3), indicating that diatoms possess diverse light-regulatory mechanisms and adaptive responses [57]. Although both light and temperature are essential for diatom growth, planktonic and benthic diatoms exhibit different preferences for PAR and water temperature. In the surface oceans, planktonic diatoms have been shown to be replaced by small phytoplankton, causing decreased primary production and carbon export, due to global warming [58]. In the Yangtze River, we observed that the richness of benthic diatoms tended to rise with increasing temperature (Fig. 3c). Moreover, the composition of benthic diatoms appeared to be strongly influenced by temperature (Additional file 1: Table S3), suggesting that benthic diatoms are more sensitive to temperature changes than planktonic diatoms in river ecosystem.

Second, as a primary driver of stream power that shapes the spatial distribution of diatoms, the channel slope not only affects competitive and succession processes among species but also alters nutritional utilization strategies and hence the production and growth of diatoms [59], leading to diatoms preferring to different flow conditions [60]. In benthic diatoms, over 40% abundance was taken by motile-guild species capable of moving fast and choosing the best microhabitat under certain circumstances (Fig. 4c). In other words, the motile guild can respond quickly to environmental changes, and benthic diatoms dominated by motile guild are more appropriate to reflect the influence of environmental heterogeneity along the river.

Third, the weak relationship between annual-averaged TN:TP and alpha-diversity of planktonic diatoms along the Yangtze River might be resulted from the selection of simpler indicator for nutrient level in a relatively oligotrophic river (Additional file 1: Figure S14). In fact, other environmental factors such as hydrodynamic condition, light, temperature, and biological predation may also play considerable roles over the full year. Moreover, species-specific responses to nitrogen and phosphorus in the production and growth of diatoms would be more complex. For example, Nitzschia palea [61] is phosphorus limited, but Chaetoceros calcitrans [62] is efficient in nitrogen assimilation, regardless the unclear mechanism on diatom utilization of nitrogen and phosphorus in different forms [63, 64].

Finally, damming in rivers is a typical anthropogenic perturbation which could profoundly modify material fluxes and biogeochemical cycles of downstream [48, 49]. Kunz et al [48] reported that sediment, carbon, nitrogen, and phosphorus were trapped by the reservoir immediately upstream of the Itezhi-Tezhi Dam, increasing the N:P ratio downstream of the dam. High flow downstream of dams has caused severe erosion of the riverbed and led to coarsening of bed materials [49]. Changes in water level also affect the light intensity and temperature to benthic diatoms. In short, a large dam disrupts the hydraulic gradient, nutrient conditions, light availability, and temperature in rivers, resulting in local changes to the environment in which benthic diatoms thrive (Additional file 1: Figure S15).

Conclusions

This study provided the first molecular biogeographic pattern of both planktonic and benthic diatoms over a continuum of 6030 km in the Yangtze River. Temporally, planktonic diatom communities demonstrated significant seasonal difference in the mainstream, and spatially dominant diatoms in water and sediment varied with landforms, such as the plateau, mountain, foothill, basin, foothill-mountain, and plain, from the river source to the estuary of the Yangtze. Comparing with the spatial dispersal process, the environmental selection process was the major driver for diatom biogeography, which could be further interpreted in terms of photosynthetically active radiation, hydraulic slope, nutrients, and human activities (the Three Gorges Dam and the Xiluodu Dam). Our study revealed that benthic diatoms represented by motile species in ecological guilds are typical consequences driven by environmental selection in a lotic-oligotrophic river, which highlights the particular importance of benthic diatoms in understanding biogeochemical cycles in world’s large river ecosystems.

Methods

Sample collection

The Yangtze River is the longest river situated wholly in Asia and the third longest in the world, with a drainage basin of 1.8 million km2. The river is over 6300 km long, has its source in the Qinghai-Tibet Plateau, and flows eastwards into the East China Sea near Shanghai. Over its length, the Yangtze River experiences great changes in landform type and hydrological regime and supports more than 588 million people [65].

In March (spring) and October (autumn) of 2014, water and sediment samples were synchronously (i.e., within 1 week) collected for planktonic and benthic diatom identification at 50 national monitoring stations (the non-plateau area) along the mainstream and six major tributaries of the Yangtze River. Finally, 84 water samples (i.e., 38 water-spring samples and 46 water-autumn samples) and 171 sediment samples (i.e., 87 sediment-spring samples and 84 sediment-autumn samples) were obtained for studying the seasonal and spatial distribution of diatom communities in the non-plateau area. To investigate the difference of diatoms in the plateau and the non-plateau areas, we implemented an additional synchronous monitoring along the source river of the Yangtze located at the Qinghai-Tibet Plateau in July 2017 (the ideal month for sampling in the plateau), collected 12 water and 12 sediment samples respectively at 12 sites. Although the most ideal sampling should be the synchronous monitoring in both plateau and non-plateau areas in the same year, the above remedial sampling is helpful considering insignificant inter-annual variations of water quality, riverine habitat, and aquatic organism in the plateau in recent years [66, 67]. Except for a very few samples missed due to restrictions of steep terrain and rapid flow as described in a previous study [30], up to four parallel samples were collected in most cases. At each sampling site, 10 L of well-mixed water was collected and then filtered onto 0.22 μm polycarbonate membranes (Millipore, USA) within 24 h. Filter membranes and sediment samples were stored in the laboratory at – 80 °C until further analysis took place.

Measurement of environmental variables

Environmental parameters including water temperature, chemical oxygen demand (COD), suspended solids (SS), dissolved oxygen (DO), pH, ammonium nitrogen (NH4-N), nitrate nitrogen (NO3-N), total nitrogen (TN), total phosphorus (TP), and dissolved organic carbon (DOC) were measured for water samples according to Environmental Quality Standards for Surface Water (GB3838-2002) recommended by the Ministry of Ecology and Environment of China. For sediment samples, pH, total organic carbon (TOC), NH4-N, NO3-N, TN, and TP contents were measured as described by Zhu et al. [68]. For each sampling sites, longitude, latitude, and altitude were recorded by a handheld GPS (Magellan, USA).

Dissolved carbon dioxide (pCO2) in water was measured using a headspace equilibration technique [69]. In short, 75 mL of river water was collected in a 100 mL polypropylene syringe. Air bubbles drawn in with the sample were eliminated by tapping and draining water while the syringe is pointing upwards. Dissolved CO2 was extracted by transferring 25 ml of ultra-high purity nitrogen at the field sites and then equilibrated with the headspace in the sample syringe by vigorous shake for 5 min. After equilibration, the headspace gas was immediately transferred to a pre-vacuum glass storage vial equipped with chlorobutyl septa. Finally, the CO2 partial pressure was measured using a gas chromatograph equipped with a thermal conductivity detector.

The information on channel slope in the Yangtze was sourced from Chen et al. [70]. Nutrients condition was represented by the atomic ratio of nitrogen to phosphorus. For water samples, we utilized dissolved total nitrogen (TN) and total phosphorus (TP) to calculate annual-averaged TN:TP for 2005 to 2014 and monitored TN:TP observed at 50 stations for 2014. For sediment samples, we utilized TN and TP to calculate monitored TN:TP observed at 50 stations for 2014.

DNA extraction, PCR amplification, and sequencing

DNA was extracted in triplicate using the FastDNA® SPIN Kit for Soil (MP Biomedicals, USA) following the manufacturer’s instructions. The triplicate DNA extracts were mixed together for later PCR amplification. Amplification of the V4 region of the 18S rDNA was performed by polymerase chain reaction (PCR) using barcoded primers DIV4for (5′-GCGGTAATTCCAGCTCCAATAG-3′) and DIV4rev3 (5′-CTCTGACAATGGAATACGAATA-3′) [12], where barcode is an eight-base sequence unique to each sample. Amplification was conducted under the following conditions: initial denaturation at 94 °C for 2 min, then 32 cycles of denaturation at 94 °C for 45 s, annealing at 50 °C for 45 s, elongation at 72 °C for 60 s, and final extension at 72 °C for 10 min, 10 °C until halted by user. PCR mixtures (20 μL volume) were prepared in triplicate contained 2 μL of 10 × buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.2 μL of rTaq polymerase, 0.2 μL of BSA, and 1 μL of 10 ng DNA sample. Amplicons were purified using the AxyPrep DNA Gel Extraction Kit (Axygen Bioscience, Union City, CA, USA) according to the manufacturer’s instructions and quantified using QuantiFluorTM-ST (Promega, USA). Adaptor was ligated onto the amplicons for the library construction. Afterwards, sample libraries were pooled in equimolar amounts and sequenced on Illumina MiSeq 2 × 250 PE platform (Majorbio Company, Shanghai, China).

Three negative control samples were used to monitor any contamination during the molecular workflow, negative filtration, DNA extraction, and PCR controls; however, no quantifiable DNA was detected for further analysis.

Bioinformatics analysis

Sequences of diatom 18S rDNA were quality-filtered using QIIME [71] as follows: (i) minimum sequence length of 300 bp, and minimum threshold quality score of Q20; (ii) maximum mismatches of 2 for matching the primer; any reads with ambiguous bases were removed; and (iii) merged pair-ended sequences that overlapped longer than 10 bp into a single sequence. UCHIME was used to remove chimeric sequences and UPARSE was used to cluster operational taxonomic units (OTUs) with 97% similarity cutoff [72].

We built a reference database of 18S rRNA reads composed of 4573 unique diatom sequences. First, we extracted all diatom sequences of 18S rRNA reads from GenBank (http://www.ncbi.nlm.nih.gov/). Second, short reads (less than 100 nucleotides) were refused access to the reference database, and redundant reads were eliminated by cd-hit to increase the taxonomy identification accuracy. Third, sequence alignment was performed by Mafft (ver 7.310) [73], then the sequences were analyzed to construct an approximately-maximum-likelihood phylogenetic tree using FastTree (ver 2.1.10) [74], and any incorrect reads discarded. Finally, a total of 4573 unique sequences were retained in our reference database.

To identify taxonomically OTUs obtained in this study against known diatom species, the BLASTN [75] program was applied to align clean 18S rRNA reads to the corrected diatom database. Those OTUs with the best BLAST hit scores, not only an e value ≤ 10−5 but also identity ≥ 80% with respect to the reference sequence were firstly selected. Then, the selected OTUs were checked by means of the phylogenetic tree, and only OTUs with correct taxonomical assignment were retained for further analysis. Clean reads were further assigned to known diatom species based on our reference database.

To estimate the community structure for each site, the Mothur program [76] was used to normalize all data sets with respect to the least-well-represented data set (11049 sequences). Alpha diversity indices (chao1, Shannon, and Goods coverage) were calculated using QIIME.

Statistical analysis

Diatom species that characterize each sample group were identified with indicator species analysis using labdsv and indval packages in R software [77]. Indicator values were calculated based on the relative frequency and relative average abundance of a given species in six types of environmental samples. Species with indicator value ≥ 0.3 and p value ≤ 0.01 were defined as indicator species at class, order, family, and genus levels. Nonmetric multidimensional scaling (NMDS) was performed to visualize the dissimilarity of different samples based on Bray-Curtis similarity matrices. Analysis of similarity (ANOSIM) was conducted to test the significance of differences among a priori sampling groups based on environmental parameters. NMDS and ANOSIM statistics were carried out using the vegan package in R. The linear discriminant analysis effect size (LEfSe) [36] was used with Kruskal-Wallis and Wilcoxon tests to discover high-dimensional biomarker and explain taxa difference at different environment conditions of PAR or channel slope. The LEfSe biomarker detection was performed in QIIME [71] using the logarithmic LDA threshold > 4 and the statistical parameters of P < 0.05. One-way analysis of variance (one-way ANOVA) was carried out to test significance of group differences using the vegan package in R.

Distance-decay patterns of diatom community similarity were described by considering geographical distance and environmental distance from the site location to river mouth among sample sites. Mantel tests were used to examine the Spearman’s rank correlation between geographical and environmental distance and diatom community similarity using Bray-Curtis distance matrices with 999 permutations in R. The geographical distance of each sampling site was calculated using ArcGIS V10.3 software. The environmental distance matrix (normalized Euclidean distance) was generated with a normalized combination of environmental variables such as water temperature, COD, SS, DO, pH, NH4-N, NO3-N, TN, TP, and DOC for water samples as well as TOC, pH, NH4-N, NO3-N, TN, and TP for sediment samples. The rate of distance-decay of diatom communities was calculated as the slope of ordinary least-squares regression line fitted to the relationship between geographic distance and community similarity. Partial Mantel tests were conducted to assess the pure effects of geographical distance (controlling for environmental distance) and environmental distance (controlling for geographic distance) on diatom community similarity with 9999 permutations.

A set of spatial variables was generated through the use of principal coordinates of neighbor matrices (PCNM) analysis based on the longitude and latitude coordinates of each sampling site [78]. The function “envfit” was run with 999 permutations to select significant environmental variables (P < 0.05). Significance testing was then assessed using the “permutest” function based on 999 permutations in R, while canonical correspondence analysis (CCA) was performed to determine the effects of selected environmental and spatial variables on diatom communities (Additional file 1: Table S3). Partial canonical correspondence analysis (pCCA) was performed to decompose the total variation in diatom community into a pure environmental component, a pure spatial component, a spatially structured environmental component, and residual variation.

Ecological guild classification

Based on their ecological characteristics, diatom species are classified into four ecological guilds (low profile, high profile, motile, and planktic guilds) [31,32,33], which are expected to respond in different ways to nutrients’ conditions and physical disturbances. A low-profile guild is defined as having high reproduction rate, low nutrient and light availability, and slow-moving diatoms. A high-profile guild possesses characteristics of high resource availability and low disturbance. A motile ecological guild has the ability to move fast and choose the best microhabitat in a given circumstance. A planktic guild adapts to lentic environments and resists sedimentation. We extended these guilds by adding supplemented classifications used in other studies [79, 80].

Photosynthetically active radiation (PAR) divisions

Solar radiation with wavelengths (400–700 nm), called photosynthetically active radiation (PAR), is able to be utilized by plants and algae through photosynthesis to convert light energy into biomass [34]. Monteith reported the linear correlation between net primary production (NPP) and PAR absorbed by green foliage [81]. Zhu et al. [35] also suggested that the spatial distribution of annual-averaged PAR is complex and inhomogeneous across China, using calculated PAR spatial data for the period 1961–2007 provided by China Meteorological Administration.

Zhu et al. [35] calculated and spatialized PAR using data simulation method [82] based on three climatic datasets, i.e., daily sunshine duration data at 740 weather stations across China for 1961–2007 and global radiation data at 122 radiation stations across China for 1961–2000 from China Meteorological Administration, and PAR observatory data at 36 field stations across China for 2004–2007 from Chinese Ecosystem Research Network. Then, PAR along the Yangtze River could be further derived. Thus, we define four zones of PAR intensity in different regions across the Yangtze River basin as follows:

  1. I.

    Very high, PAR > 32 mol m−2 d−1;

  2. II.

    High, 26 < PAR ≤ 32 mol m−2 d−1;

  3. III.

    Medium, 23 < PAR ≤ 26 mol m−2 d−1;

  4. IV.

    Low, PAR ≤ 23 mol m−2 d−1.

Availability of data and materials

Complete datasets supporting the findings of this article are in the NCBI Sequence Read Archive (SRA) database (Accession Number: SRP153344).

Abbreviations

OTUs:

Operational taxonomic units

HTS:

High-throughput sequencing

NMDS:

Nonmetric multidimensional scaling

ANOSIM:

Analysis of similarity

One-way ANOVA:

One-way analysis of variance

LefSe:

Linear discriminant analysis effect size

PAR:

Photosynthetically active radiation

TN:TP:

Ratio of total nitrogen to total phosphate

References

  1. 1.

    Armbrust EV. The life of diatoms in the world’s oceans. Nature. 2009;459:185–92.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Field CB, Behrenfeld MJ, Randerson JT, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science. 1998;281:237–40.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Sims PA, Mann DG, Medlin LK. Evolution of the diatoms: insights from fossil, biological and molecular data. Phycologia. 2006;45:361–402.

    Article  Google Scholar 

  4. 4.

    Mann DG, Droop SJM. Biodiversity, biogeography and conservation of diatoms. Hydrobiologia. 1996;336:19–32.

    Article  Google Scholar 

  5. 5.

    Oeding S, Taffs KH. Developing a regional diatom index for assessment and monitoring of freshwater streams in sub-tropical Australia. Ecol Indic. 2017;80:135–46.

    CAS  Article  Google Scholar 

  6. 6.

    Cermeño P, Falkowski PG. Controls on diatom biogeography in the ocean. Science. 2009;325:1539–41.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  7. 7.

    Orlando Necchi JR. River Algae. Necchi JR O, editor. Cham: Springer International Publishing; 2016.

  8. 8.

    Kelly MG. Use of the trophic diatom index to monitor eutrophication in rivers. Water Res. 1998;32:236–42.

    CAS  Article  Google Scholar 

  9. 9.

    Zimmermann J, Glöckner G, Jahn R, Enke N, Gemeinholzer B. Metabarcoding vs. morphological identification to assess diatom diversity in environmental studies. Mol Ecol Resour. 2015;15:526–42.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Evans KM, Wortley AH, Mann DG. An assessment of potential diatom “Barcode” genes (cox1, rbcL, 18S and ITS rDNA) and their effectiveness in determining relationships in Sellaphora (Bacillariophyta). Protist. 2007;158:349–64.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Moniz MBJ, Kaczmarska I. Barcoding of diatoms: nuclear encoded ITS revisited. Protist. 2010;161:7–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Visco JA, Apotheloz-Perret-Gentil L, Cordonier A, Esling P, Pillet L, Pawlowski J. Environmental monitoring: inferring the diatom index from next-generation sequencing data. Environ Sci Technol. 2015;49:7597–605.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Malviya S, Scalco E, Audic S, Vincent F, Veluchamy A, Poulain J, et al. Insights into global diatom distribution and diversity in the world’s ocean. Proc Natl Acad Sci U S A. 2016;113:E1516–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Boyd PW, Watson AJ, Law CS, Abraham ER, Trull T, Murdoch R, et al. A mesoscale phytoplankton bloom in the polar Southern Ocean stimulated by iron fertilization. Nature. 2000;407:695–702.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Tsuda A, Takeda S, Saito H, Nishioka J, Nojiri Y, Kudo I, et al. A mesoscale iron enrichment in the western subarctic pacific induces a large centric diatom bloom. Science. 2003;300:958–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Heinsalu A, Alliksaar T, Leeben A, Nõges T. Sediment diatom assemblages and composition of pore-water dissolved organic matter reflect recent eutrophication history of Lake Peipsi (Estonia/Russia). Hydrobiologia. 2007;584:133–43.

    CAS  Article  Google Scholar 

  17. 17.

    Kitner M, Poulíčková A. Littoral diatoms as indicators for the eutrophication of shallow lakes. Hydrobiologia. 2003;506–509:519–24.

    Article  Google Scholar 

  18. 18.

    Liu J, Soininen J, Han B-P, Declerck SAJ. Effects of connectivity, dispersal directionality and functional traits on the metacommunity structure of river benthic diatoms. J Biogeogr. 2013;40:2238–48.

    Article  Google Scholar 

  19. 19.

    Kheiri S, Solak CN, Edlund MB, Spaulding S, Nejadsattari T, Asri Y, et al. Biodiversity of diatoms in the Karaj River in the Central Alborz, Iran. Diatom Res. 2018;33:355–80.

    Article  Google Scholar 

  20. 20.

    Centis B, Tolotti M, Salmaso N. Structure of the diatom community of the river Adige (North-Eastern Italy) along a hydrological gradient. Hydrobiologia. 2010;639:37–42.

    CAS  Article  Google Scholar 

  21. 21.

    Rimet F. Benthic diatom assemblages and their correspondence with ecoregional classifications: case study of rivers in North-Eastern France. Hydrobiologia. 2009;636:137–51.

    Article  Google Scholar 

  22. 22.

    Dalu T, Wasserman RJ, Magoro ML, Mwedzi T, Froneman PW, Weyl OLF. Variation partitioning of benthic diatom community matrices: effects of multiple variables on benthic diatom communities in an Austral temperate river system. Sci Total Environ. 2017;601–602:73–82.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  23. 23.

    Kireta AR, Reavie ED, Sgro GV, Angradi TR, Bolgrien DW, Hill BH, et al. Planktonic and periphytic diatoms as indicators of stress on great rivers of the United States: testing water quality and disturbance models. Ecol Indic. 2012;13:222–31.

    Article  Google Scholar 

  24. 24.

    Martiny JBH, Bohannan BJM, Brown JH, Colwell RK, Fuhrman JA, Green JL, et al. Microbial biogeography: putting microorganisms on the map. Nat Rev Microbiol. 2006;4:102–12.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Hanson CA, Fuhrman JA, Horner-Devine MC, Martiny JBH. Beyond biogeographic patterns: processes shaping the microbial landscape. Nat Rev Microbiol. 2012;10:497–506.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Keck F, Franc A, Kahlert M. Disentangling the processes driving the biogeography of freshwater diatoms: a multiscale approach. J Biogeogr. 2018;45:1582–92.

    Article  Google Scholar 

  27. 27.

    Heino J, Bini LM, Karjalainen SM, Mykrä H, Soininen J, Vieira LCG, et al. Geographical patterns of micro-organismal community structure: are diatoms ubiquitously distributed across boreal streams? Oikos. 2010;119:129–37.

    Article  Google Scholar 

  28. 28.

    Soininen J. Determinants of benthic diatom community structure in boreal streams: the role of environmental and spatial factors at different scales. Int Rev Hydrobiol. 2004;89:139–50.

    Article  Google Scholar 

  29. 29.

    Tréguer P, Bowler C, Moriceau B, Dutkiewicz S, Gehlen M, Aumont O, et al. Influence of diatom diversity on the ocean biological carbon pump. Nat Geosci. 2018;11:27–37.

    Article  CAS  Google Scholar 

  30. 30.

    Liu T, Zhang AN, Wang J, Liu S, Jiang X, Dang C, et al. Integrated biogeography of planktonic and sedimentary bacterial communities in the Yangtze River. Microbiome. 2018;6:1–14.

  31. 31.

    Passy SI. Diatom ecological guilds display distinct and predictable behavior along nutrient and disturbance gradients in running waters. Aquat Bot. 2007;86:171–8.

    Article  Google Scholar 

  32. 32.

    Rimet F, Bouchez A. Life-forms, cell-sizes and ecological guilds of diatoms in European rivers. Knowl Manag Aquat Ecosyst. 2012;406:01.

  33. 33.

    B-Béres V, Lukács Á, Török P, Kókai Z, Novák Z, T-Krasznai E, et al. Combined eco-morphological functional groups are reliable indicators of colonisation processes of benthic diatom assemblages in a lowland stream. Ecol Indic. 2016;64:31–8.

    Article  Google Scholar 

  34. 34.

    Jacovides CP, Timvios FS, Papaioannou G, Asimakopoulos DN, Theofilou CM. Ratio of PAR to broadband solar radiation measured in Cyprus. Agric For Meteorol. 2004;121:135–40.

    Article  Google Scholar 

  35. 35.

    Zhu X, He H, Liu M, Yu G, Sun X, Gao Y. Spatio-temporal variation of photosynthetically active radiation in China in recent 50 years. J Geogr Sci. 2010;20:803–17.

    Article  Google Scholar 

  36. 36.

    Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Montagnes DJS, Franklin DJ. Effect of temperature on diatom volume, growth rate, and carbon and nitrogen content: reconsidering some paradigms. Limnol Oceanogr. 2001;46:2008–18.

    CAS  Article  Google Scholar 

  38. 38.

    Yu G-A, Brierley G, Huang HQ, Wang Z, Blue B, Ma Y. An environmental gradient of vegetative controls upon channel planform in the source region of the Yangtze and Yellow Rivers. Catena. 2014;119:143–53.

    Article  Google Scholar 

  39. 39.

    Redfield AC, Ketchum BH, Richards FA. The influence of organisms on the composition of seawater. Sea. New York: Wiley-Interscience; 1963.

  40. 40.

    Mikaelyan AS, Pautova LA, Chasovnikov VK, Mosharov SA, Silkin VA. Alternation of diatoms and coccolithophores in the north-eastern Black Sea: a response to nutrient changes. Hydrobiologia. 2015;755:89–105.

    CAS  Article  Google Scholar 

  41. 41.

    Arai H, Fukushima T. Impacts of long-term increase in silicon concentration on diatom blooms in Lake Kasumigaura, Japan. Ann Limnol. 2014;50:335–46.

    Article  Google Scholar 

  42. 42.

    Cohn SA, Farrell JF, Munro JD, Ragland RL, Weitzell RE, Wibisono BL. The effect of temperature and mixed species composition on diatom motility and adhesion. Diatom Res. 2003;18:225–43.

    Article  Google Scholar 

  43. 43.

    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.

  44. 44.

    Tekwani N, Majdi N, Mialet B, Tornés E, Urrea-Clos G, Buffan-Dubau E, et al. Contribution of epilithic diatoms to benthic-pelagic coupling in a temperate river. Aquat Microb Ecol. 2013;69:47–57.

    Article  Google Scholar 

  45. 45.

    Jäger CG, Diehl S. Resource competition across habitat boundaries: asymmetric interactions between benthic and pelagic producers. Ecol Monogr. 2014;84:287–302.

    Article  Google Scholar 

  46. 46.

    Vannote RL, Minshall GW, Cummins KW, Sedell JR, Cushing CE. The river continuum concept. Can J Fish Aquat Sci. 1980;37:130–7.

    Article  Google Scholar 

  47. 47.

    Ward JV, Stanford JA. The serial discontinuity concept of lotic ecoystems. In: Fontaine, T.D., Bartell SM, editor. Dynamics of Lotic Ecosystems. Ann Arbor Science Publications. Ann Arbor, Michigan; 1983.

  48. 48.

    Kunz MJ, Wüest A, Wehrli B, Landert J, Senn DB. Impact of a large tropical reservoir on riverine transport of sediment, carbon, and nutrients to downstream wetlands. Water Resour Res. 2011;47:W12531.

    Article  CAS  Google Scholar 

  49. 49.

    Graf WL. Downstream hydrologic and geomorphic effects of large dams on American rivers. Geomorphology. 2006;79:336–60.

    Article  Google Scholar 

  50. 50.

    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.

    CAS  Article  Google Scholar 

  51. 51.

    Jäger CG, Borchardt D. Longitudinal patterns and response lengths of algae in riverine ecosystems: a model analysis emphasising benthic-pelagic interactions. J Theor Biol. 2018;442:66–78.

    PubMed  Article  PubMed Central  Google Scholar 

  52. 52.

    Gillett N, Pan Y, Parker C. Should only live diatoms be used in the bioassessment of small mountain streams? Hydrobiologia. 2009;620:135–47.

    Article  Google Scholar 

  53. 53.

    Smucker NJ, Vis ML. Spatial factors contribute to benthic diatom structure in streams across spatial scales: considerations for biomonitoring. Ecol Indic. 2011;11:1191–203.

    Article  Google Scholar 

  54. 54.

    Chen L, Liu S, Chen Q, Zhu G, Wu X, Wang J, et al. Anammox response to natural and anthropogenic impacts over the Yangtze River. Sci Total Environ. 2019;665:171–80.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Domingues RB, Anselmo TP, Barbosa AB, Sommer U, Galvão HM. Light as a driver of phytoplankton growth and production in the freshwater tidal zone of a turbid estuary. Estuar Coast Shelf Sci. 2011;91:526–35.

    CAS  Article  Google Scholar 

  56. 56.

    Buma AGJ, Engelen AH, Gieskes WWC. Wavelength-dependent induction of thymine dimers and growth rate reduction in the marine diatom Cyclotella sp. exposed to ultraviolet radiation. Mar Ecol Prog Ser. 1997;153:91–7.

    CAS  Article  Google Scholar 

  57. 57.

    Depauw FA, Rogato A, D’Alcalá MR, Falciatore A. Exploring the molecular basis of responses to light in marine diatoms. J Exp Bot. 2012;63:1575–91.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  58. 58.

    Bopp L, Aumont O, Cadule P, Alvain S, Gehlen M. Response of diatoms distribution to global warming and potential implications : a global model study. Geophys Res Lett. 2005;32:1–4.

    Article  CAS  Google Scholar 

  59. 59.

    Peters F, Arin L, Marrasé C, Berdalet E, Sala MM. Effects of small-scale turbulence on the growth of two diatoms of different size in a phosphorus-limited medium. J Mar Syst. 2006;61:134–48.

    CAS  Article  Google Scholar 

  60. 60.

    Mitrovic SM, Chessman BC, Davie A, Avery EL, Ryan N. Development of blooms of Cyclotella meneghiniana and Nitzschia spp. (Bacillariophyceae) in a shallow river and estimation of effective suppression flows. Hydrobiologia. 2008;596:173–85.

    Article  Google Scholar 

  61. 61.

    Han L, Xu B, Qi F, Chen Z. Effect of nitrogen/phosphorus concentration on algal organic matter generation of the diatom Nitzschia palea: total indicators and spectroscopic characterization. J Environ Sci (China). 2016;47:130–42.

    Article  Google Scholar 

  62. 62.

    Tantanasarit C, Englande AJ, Babel S. Nitrogen, phosphorus and silicon uptake kinetics by marine diatom Chaetoceros calcitrans under high nutrient concentrations. J Exp Mar Biol Ecol. 2013;446:67–75.

    CAS  Article  Google Scholar 

  63. 63.

    Zeng Q, Qin L, Bao L, Li Y, Li X. Critical nutrient thresholds needed to control eutrophication and synergistic interactions between phosphorus and different nitrogen sources. Environ Sci Pollut Res. 2016;23:21008–19.

    CAS  Article  Google Scholar 

  64. 64.

    Cotner JB Jr, Wetzel RG. Uptake of dissolved inorganic and organic phosphorus compounds by phytoplankton and bacterioplankton. Limnol Oceanogr. 1992;37:232–43.

    CAS  Article  Google Scholar 

  65. 65.

    Qu S, Wang L, Lin A, Zhu H, Yuan M. What drives the vegetation restoration in Yangtze River basin, China: climate change or anthropogenic factors? Ecol Indic. 2018;90:438–50.

    Article  Google Scholar 

  66. 66.

    Huang Z, Liu Y, Zhao W, Zhao L. Discussion on recent spatial-temporal distribution of water quality in Changjiang River source area. J Yangtze River Sci Res Inst. 2016;33:46–50 (in Chinese).

    Google Scholar 

  67. 67.

    Zhao L, Li W, Lin L, Guo W, Zhao W, Tang X, et al. Field investigation on river hydrochemical characteristics and larval and juvenile fish in the source region of the Yangtze River. Water (Switzerland). 2019;11:1342.

    Article  Google Scholar 

  68. 68.

    Zhu G, Wang S, Wang Y, Wang C, Risgaard-Petersen N, Jetten MS, et al. Anaerobic ammonia oxidation in a fertilized paddy soil. ISME J. 2011;5:1905–12.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  69. 69.

    Beaulieu JJ, Shuster WD, Rebholz JA. Controls on gas transfer velocities in a large river. J Geophys Res Biogeosci. 2012;117:G02007.

    Article  CAS  Google Scholar 

  70. 70.

    Chen Z, Li J, Shen H, Zhanghua W. Yangtze River of China: historical analysis of discharge variability and sediment flux. Geomorphology. 2001;41:77–91.

    Article  Google Scholar 

  71. 71.

    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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.

    CAS  PubMed  Article  Google Scholar 

  73. 73.

    Yamada KD, Tomii K, Katoh K. Application of the MAFFT sequence alignment program to large data - reexamination of the usefulness of chained guide trees. Bioinformatics. 2016;32:3246–51.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Price MN, Dehal PS, Arkin AP. FastTree 2 - approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  75. 75.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: anew generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Dufrtne M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr. 1997;67:345–66.

    Google Scholar 

  78. 78.

    Borcard D, Legendre P. All-scale spatial analysis of ecological data by means of principal coordinates of neighbor matrices. Ecol Model. 2002;153:51–68.

    Article  Google Scholar 

  79. 79.

    Berthon V, Bouchez A, Rimet F. Using diatom life-forms and ecological guilds to assess organic pollution and trophic level in rivers: a case study of rivers in south-eastern France. Hydrobiologia. 2011;673:259–71.

    CAS  Article  Google Scholar 

  80. 80.

    B-Béres V, Török P, Kókai Z, Lukács Á, T-Krasznai E, Tóthmérész B, et al. Ecological background of diatom functional groups: comparability of classification systems. Ecol Indic. 2017;82:183–8.

    Article  Google Scholar 

  81. 81.

    Monteith JL. Solar-radiation and productivity in tropical ecosystems. J Appl Ecol. 1972;9:747–66.

    Article  Google Scholar 

  82. 82.

    Lin Z-H, Mo X-G, Li H-X, Li H-B. Comparison of three spatial interpolation methods for climate variables in China. Dili Xuebao/Acta Geogr Sin. 2002;57:47–54 (in Chinese).

    Google Scholar 

Download references

Acknowledgements

Supports from BIC-ESAT and Majorbio Company (Shanghai, China) are gratefully acknowledged.

Funding

Financial support from the National Natural Science Foundation of China (No. 91647211 and 51721006) is appreciated.

Author information

Affiliations

Authors

Contributions

JRN designed the research. JWW and QXL performed the research. JWW, QXL, JRN, and AGLB wrote the paper. XZ and YL contributed new ideas and information. All of the authors contributed to interpretation of the findings and approved the final version of the manuscript.

Corresponding author

Correspondence to Jinren Ni.

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 diatom richness for 279 samples. Figure S2. Phylogenetic distribution of reference sequence and abundant OTUs. Figure S3. Alpha diversity index per diatom community for each sample type. Figure S4. Nonmetric multidimensional scaling diagram of Bray-Curtis dissimilarities between diatom communities for all samples. Figure S5. ANOSIM statistics concerning differences in diatom communities within and between sample types. Figure S6. Diatom genera exhibiting significant seasonal differences in water (a) and sediment (b) samples. Figure S7. Spearman relationships between diatoms and dissolved carbon dioxide. Figure S8. Biogeographical distribution of diatoms at class level throughout the mainstream of the Yangtze for: (a) water-spring, (b) water-autumn, (c) sediment-spring, and (d) sediment-autumn samples. Figure S9. Proportion of indicator diatoms in each sample type obtained using indicator taxa analysis at Class (a), Order (b), Family (c), and Genus (d) levels. Figure S10. Biogeographical distribution of ecological guilds throughout the mainstream of the Yangtze River. Figure S11. Distance-decay relationships between community similarity and geographic distance for water (a) and sediment (b) samples. Figure S12. Variation in community composition explained by environmental, spatial, and spatially structured environmental components. Figure S13. LEfSe cladogram of planktonic (a) and benthic (b) diatom communities for the three channel slope regions. Figure S14. Relationships between Shannon diversity and TN:TP for water-spring (a), water-autumn (b), sediment-spring (c), and sediment-autumn (d) samples. Distance relationship of TN:TP for sampling sites along the mainstream (e). Figure S15. Significant differences in abundance of benthic diatoms upstream and downstream of Xiluodu Dam (a) and Three Gorges Dam (b). Table S1. Numbers of Indicator species and Top Indicator species across sample sites. Table S2. Partial Mantel test for Spearman correlations between community similarity and geographic and environmental distances. Table S3. Effects of selected factors on diatom community using canonical correspondence analysis.

Rights and permissions

Open Access This 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, J., Liu, Q., Zhao, X. et al. Molecular biogeography of planktonic and benthic diatoms in the Yangtze River. Microbiome 7, 153 (2019). https://doi.org/10.1186/s40168-019-0771-x

Download citation

Keywords

  • Diatoms
  • Biogeography
  • Water
  • Sediment
  • Landform
  • Environmental driver
  • Human interference
  • Yangtze River