- 1Phycology Working Group, Faculty of Biology, University of Duisburg-Essen, Essen, Germany
- 2University of Belgrade, Institute of Chemistry, Technology and Metallurgy, National Institute of the Republic of Serbia, Belgrade, Serbia
- 3Aquatic Ecosystem Research, Faculty of Biology, University of Duisburg-Essen, Essen, Germany
Introduction: Microalgae form an essential group of benthic organisms that respond swiftly to environmental changes. They are widely used as bioindicators of anthropogenic stressors in freshwater ecosystems. We aimed to assess the responses of microalgae communities to multiple environmental stressors in the Kinzig River catchment, home to a long-term ecological monitoring site, in Germany.
Methods: We used a photosynthetic biomass proxy alongside community composition of diatoms assessed by digital light microscopy, and of microalgae by 18S-V9 amplicon sequencing, to characterise microalgae at 19 sampling sites scattered across the catchment.
Results: Our results revealed significant effects of physical and chemical factors on microalgae biomass and community compositions. We found that conductivity, water temperature and pH were the most important factors affecting microalgae community composition, as observed in both microscopy and amplicon analysis. In addition to these three variables, the effect of total phosphate on all microalgae, together with water discharge on the diatom (Bacillariophyta) communities, as assessed by amplicon analysis, may reveal taxon-specific variations in the ecological responses of different microalgal groups.
Discussion: Our results highlighted the complex relationship between various environmental variables and microalgae biomass and community composition. Further investigations, involving the collection of time series data, are required to fully understand the underlying biotic and abiotic parameters that influence these microalgae communities.
1 Introduction
The complex assemblages of microorganisms colonizing all submerged substrates also referred to as benthic microalgae biofilms, periphyton, or microphytobenthos (MPB), play an essential role in the biogeochemical cycles, primary production, and sediment stabilisation (B-Béres et al., 2023). MPB assemblages are mainly dominated by diatoms (Bacillariophyta), but also other eukaryotes, including other microalgae such as Chlorophyta, Cyanobacteria and Euglenophyta, as well as other mixo- and heterotrophic protists, and fungi, may be present in smaller proportions (Underwood, 2010; Lemley et al., 2016; Dalu et al., 2020). At first glance, these tiny organisms may seem insignificant, but their presence and importance in rivers should not be underestimated. Through photosynthesis, these microscopic organisms act as microforests, converting sunlight into energy, and producing oxygen and organic matter. They contribute substantially to the global carbon and nitrogen cycles and provide essential ecosystem services that sustain diverse forms of life (Castro-Català et al., 2020; B-Béres et al., 2023). The environment and other organisms with which MPB biofilms interact (such as grazers) affect their species richness, community composition, and productivity (Birk et al., 2020). The ability of MPB biofilms to grow in different habitats and conditions also explains the compositional variability and complexity observed in their assemblages.
Diatoms are secondary endosymbiotic eukaryotic algae belonging to the phylum Stramenopiles. Their glassy, siliceous cell walls known as frustules have been the focus of MPB community studies, mainly by observing their morphology under a light or a scanning electron microscope (Morales et al., 2001; Smol and Stoermer, 2010). Due to their short life cycles and rapid response to abiotic and biotic changes, MPB, particularly diatoms, are valuable bioindicators that can provide important insights into the health of aquatic ecosystems (Smol and Stoermer, 2010). Therefore, since 2000, the European Commission Water Framework Directive (European Commission, 2000) recommends the use of Biological Quality Elements (BQE), including diatoms, for the ecological assessment of surface waters.
Generally, studies of MPB biofilms have focused on the study of diatoms using light microscopy to analyse their community compositions. As diatoms are often the main component of these benthic communities, this approach has provided valuable insights into their ecological dynamics because they swiftly respond to both present and past environmental effects (Round, 1991; Passy, 2007; Soininen et al., 2016; Viso and Blanco, 2023). Because of their well-preserved silica cell walls, diatoms are considered powerful indicators of environmental change, including acidification, eutrophication, and climate change, both in neolimnological and paleolimnological studies (Dixit et al., 1992; van Dam et al., 1994). The robust frustules of diatoms form a stable ecological archive, enabling long-term studies and re-examination of materials embedded in permanent slides. In contrast to the degradation susceptibility of samples for amplicon sequencing, diatom valves are valuable for reconstructing past environmental changes. Several studies have examined the responses of diatom communities to various anthropogenic disturbances (Walsh and Wepener, 2009; Hlúbiková et al., 2014; Teittinen et al., 2015; Stenger-Kovács et al., 2020). Diatom responses to a range of stressors, including heavy metal pollution (Lavoie et al., 2018), habitat fragmentation (Szczepocka et al., 2021) and climate change (Tornés et al., 2022) are well documented. This wealth of knowledge has greatly improved our understanding of diatom ecology and its resilience or vulnerability under different stressful conditions.
However, traditional identification and counting of diatoms using light microscopy offers a somewhat biased perspective. Often overlooked is the diversity and functional potential of the broader range of other groups of organism that compose MPB biofilms, which include not only diatoms but also other microalgae, mixotrophic and heterotrophic protists and even fungi. This narrow focus on diatoms may lead us to underestimate the true ecological complexity and functional diversity of MPB biofilms. Additionally, it is important to study not only the community composition but also the functional traits (Stenger-Kovács et al., 2020). Traits related to photosynthesis are particularly significant as they underpin productivity and biogeochemical cycling in freshwater ecosystems. A crucial aspect that remains underexplored is the potential differential impact of anthropogenic activities on microphytobenthos communities. As we continue to reshape the planet in the Anthropocene, understanding the nuanced effects of multiple anthropogenic stressors has become increasingly important. Human activities introduce different types of anthropogenic pressures in river catchments. These human-induced stressors have various effects on microphytobenthos communities (Feld et al., 2016; Birk et al., 2020). For example, in urban areas, with high population densities, extensive infrastructure and industrial activities, cocktail of point source pollutants may enter in freshwater ecosystems, that include wastewater treatment plant effluents, channel erosion, nutrients load, and inflow of contaminants such as heavy metals and other chemicals (Trábert et al., 2020). This is in contrast with the rural landscapes dominated by agricultural land use, which is characterized typically by more diffuse pollution loads, such as the influx of nutrients, agrochemicals, pesticides and fine sediments, as well as hydromorphological alterations (Feld et al., 2016; Schürings et al., 2022). Yet, the effects of these multiple anthropogenic stressors (chemical and hydro-morphological) on microphytobenthos communities remain largely unknown.
In this study, we investigated how stream microphytobenthos communities respond to environmental variables in the Kinzig catchment, which is part of the Rhine-Main-Observatory (RMO), a long-term ecological research site (LTER) in central Germany. We focused on the primary environmental correlates of periphyton distribution in this river network. Our objectives were to: (1) assess the dynamics of photosynthetic biomass as a result of water physical and chemical conditions; (2) examine microalgae community composition and diversity metrics (abundances, richness, etc.) in relation to environmental variables and associated changes in ecological status; (3) identify the potential environmental factors that influence community composition in response to environmental conditions and ecological status, as compared to the macrozoobenthos studies of the LTER (Nguyen et al., 2023). We postulated that photosynthetic biomass and composition of MPB communities in the Kinzig catchment correlate with environmental factors. Higher biomass and reduced diversity will be associated with increased nutrient availability and eutrophication in the system. Similarly, hydrology (e.g. high flow velocity and depth, producing cloudy and low-light waters), will result also in a reduction in algal biomass and change in community composition. We seek to understand the resilience and adaptability of the MPB community by contrasting their responses to multiple anthropogenic stressors. In the face of numerous human pressures, this knowledge may help develop more potent conservation strategies for these essential freshwater ecosystem components. In addition to shedding light on the wider effects of multiple stressors on freshwater ecosystems, our study advances our understanding of MPB ecology.
2 Materials and methods
2.1 Study stie
The Kinzig catchment area (1058 km2) is a long-term ecosystem research site (LTER) in the Rhine-Main Observatory (RMO) (Nguyen et al., 2023). This low-mountain river system flows through different types of land use introducing various stressors. We are currently unaware of any ongoing restoration initiative in this river catchment, although land use in the study area has changed over time. Germany is actively working to reduce nutrient emissions in all catchment areas, in addition to ongoing efforts to address climate change (Nguyen et al., 2023). These factors could alter the structural composition of the microphytobenthos community over time. This study was part of the Collaborative Research Centre Multilevel Response to Stressor Increase and Release in Stream Ecosystems (CRC RESIST, https://sfb-resist.de/), which investigates the effects of global change in rivers and streams using field studies, mesocosm field and lab experiments and modelling approaches on several organism groups and scopes. We collected 19 samples in spring 2021 from ten tributaries selected within the Kinzig catchment (Figure 1). Sites were sampled based on adjacent land use.
2.2 Sample collection, preparation and processing
In spring 2021, nineteen samples were collected from various streams across the Kinzig catchment and examined using a range of methods, including a recently developed digital version of diatom light microscopy for improved resolution of community composition (Burfeid-Castellanos et al., 2022), 18S-V9 amplicon sequencing for more comprehensive taxonomic characterisation of MPB communities (Amaral-Zettler et al., 2009) and in situ chlorophyll fluorescence measurements as a proxy for overall and group-specific photosynthetic biomass estimation (Kahlert and McKie, 2014). At least five separate chlorophyll (Chl) fluorescence measurements were made in situ on each sampled rock using a BenthoTorch (bbe Moldaenke GmbH, Schwentinental, Germany). The measurements were made using factory-programmed settings, and total Chl-a concentrations were given as the sum of biomass values for the three photosynthetic algal groups, namely cyanobacteria, green algae and diatoms. Biofilm samples were collected according to the recommendations and guidelines for sampling benthic diatoms from European rivers (Kelly et al., 1998; CEN. UNE-EN, 2014). A new, clean toothbrush was used to scrape approximately 20 cm2 from each of the five cobble stones sampled at each site. Approximately 20 mL of water from the river was used to collect the scrapings from the 5 cobble stones. As these samples were collected for both digital microscopy and DNA extraction, absolute ethanol (99%) was added, resulting in an approximate final concentration of approximately 75%, transported on ice and then frozen at -20°C once in the laboratory until they were divided for preparation for further analysis by digital microscopy and molecular methods. Water temperature, electrical conductivity and pH were measured in situ using a PCE-PHD 1 pH meter (Meschede, Germany). Other physical and chemical properties of the water including dissolved oxygen (DO), oxygen saturation, average velocity, average depth, discharge and water level were also measured in situ during biofilm sampling. In the laboratory, the Deutsches Institut für Normung (DIN) ISO norms were followed to determine the content of ortho-phosphate by photometry using a test kit, total phosphate (TP) via ICP-MS (PerkinElmer), ammonium (NH4+) by photometry according to DIN 38406-5 (E5 – 1), nitrite by photometry according to DIN EN 26 777, nitrate, chloride and sulfate via ion chromatography, total nitrogen (TN) was calculated from nitrite, nitrate and ammonium values.
2.3 Diatom analysis by digital microscopy
We chose the digital microscopy workflow because it represents a good replacement for the “traditional” light microscopy analysis of diatom community composition. It combines high-resolution slide scanning with web-based virtual slide annotation. As Burfeid-Castellanos et al. (2022) compared the two techniques and found that while the digital microscopy workflow sometimes provided finer taxonomic resolution, it generally produced results comparable to traditional light microscopy. Digital microscopy makes quality control, taxonomic updates, multi-expert collaboration, and even taxonomic intercalibration much easier to accomplish. It also improves reproducibility and transparency. Sample preparation and processing for digital microscopy followed the procedures described by Burfeid-Castellanos and coworkers (Burfeid-Castellanos et al., 2022). After five prewash cycles, the biofilm samples were digested with hydrogen peroxide and hydrochloric acid on a hot plate (Taylor et al., 2007). After seven rinsing cycles, centrifugation at 1200 rpm for 3 minutes (Eppendorf Centrifuge 5427 R; Eppendorf, Hamburg, Germany) and decanting and replenishing with distilled water, a 400 µL aliquot of each cleaned sample was pipetted onto a coverslip (15 x 15 mm, #1,5), allowed to dry and then mounted onto permanent glass slides using Naphrax resin (refractive index = 1.71, Biologie-Bedarf Thorns, Deggendorf, Germany). To produce virtual slides, an area of approximately 5 x 5 mm2 on the part of the permanent glass slide with the highest density of evenly distributed diatom valves was scanned using a VS200 slide scanner (EVIDENT, Tokyo, Japan) using a high-resolution oil immersion objective (UPlanXApo 60x/1.42) with the ASW 3.1 software (Olympus Soft Imaging Solutions GmbH, Münster, Germany).
These virtual slides were analysed for diatom counting and identification with the BIIGLE 2.0 web browser-based image annotation platform (Langenkämper et al., 2017; Zurowietz and Nattkemper, 2021) as described previously (Burfeid-Castellanos et al., 2022). Using a combination of general (Bey and Ector, 2013; Cantonati et al., 2017; Spaulding et al., 2022) and specific taxonomic literature (Cox, 1995; Levkov et al., 2013; Trobajo et al., 2013), at least 400 diatom valves were counted and identified at the highest level of possible taxonomic resolution, mostly at species level. For each sample, a list of the taxa and their relative abundances was generated to calculate the IPS (“Indice de Polluosensibilité Spécifique”) (Prygiel and Coste, 1993), using OMNIDIA software (Lecointe et al., 1993) and to run multivariate statistical analyses using R open source software environment for statistical analyses (R Core Team, 2023).
2.4 Molecular analyses
2.4.1 DNA extraction and sequencing
Sample preparation and DNA extraction were performed following the silica bead-based extraction protocol described by Buchner (Buchner, 2022a; 2022b). Each sample was processed in 2 replicates. Genomic DNA amplification of the V9 hypervariable region of the small subunit ribosomal RNA genes (ca. 130 bp) was performed using the universal specific forward primer 1389F (5′-TTGTACACACCGCCC-3′) and the eukaryotic specific reverse primer 1510R (5′-CCTTCYGCAGGTTCACCTAC-3′) (Amaral-Zettler et al., 2009), following the protocol of Stoeck et al (Stoeck et al., 2010). We selected the 18S-V9 region as a marker to obtain a broader taxonomic resolution of the eukaryotic organisms present in the biofilm. The resulting library was sequenced on an Illumina MiSeq at CeGat GmbH (Tübingen, Germany).
2.4.2 Bioinformatics analysis
The Natrix2 workflow (Deep et al., 2023) was used to perform bioinformatics analyses on microphytobenthos Illumina amplicon sequencing data. The operational taxonomic units (OTUs) variant of the workflow was implemented using the SWARM v3.0.0 (Mahé et al., 2015) clustering algorithm. Paired-end reads were assembled using the simple Bayesian algorithm in PANDAseq v2.11 (Masella et al., 2012). Primers were trimmed with cutadapt v3.2 (Martin, 2011), and then filtered with a PANDAseq threshold score of 0.9, minimum length of 77, and maximum length 196 nucleotides. Sequence dereplication was performed using the CD-HIT v4.8.1 (Fu et al., 2012) algorithm at 100% similarity and chimeric sequences were identified and removed using VSEARCH v2.15.2 (Rognes et al., 2016). The split-sample approach (Lange et al., 2015) was applied to reduce erroneous sequences without strict abundance cut-offs. Statistical analysis of the amplicon sequencing data resulting from the split-sample process was conducted using the AmpliconDuo v1.1 R package (Lange et al., 2015). OTUs were generated by clustering sequences with Swarm v3.0.0 (Mahé et al., 2015). Then, mothur v1.40.5 (Schloss et al., 2009) was used to align OTUs against the protist ribosomal reference database (PR2) v.4.14.0 (Guillou et al., 2013) to assign taxonomy. Finally, mumu (https://github.com/frederic-mahe/mumu), a C++ implementation of lulu (Frøslev et al., 2017), was used for post-clustering curation of the periphyton amplicon sequencing data. The replicates of sequenced samples were merged into one representative sample by calculating the sum of reads for each OTU and the negative controls with maximum reads were subtracted from the sum of both samples to standardize the proportion of sequences. This is important to reduce redundancy in the resulting dataset for subsequent statistical analyses. The resulting OTU table with their relative read abundances is used for subsequent multivariate statistical analyses.
2.5 Data analysis
We further examined the results from all above-described analyses using R (v4.2.2, R Core Team, 2023). Data preparation and exploratory data analysis were conducted using various packages including ‘tidyverse’ (v2.0.0, Wickham et al., 2019), ‘corrplot’ (v0.92, Wei and Simko, 2021), ‘factoextra’ (v1.0.7, Kassambara and Mundt, 2020), ‘ggstatsplot’ (v0.12.0, Patil, 2021). Correlation analysis was performed to determine the relationship between photosynthetic biomass and environmental variables, visualized using ‘corrplot’ package. The ‘vegan’ package (v2.6-4, Oksanen et al., 2022) was used to conduct multivariate statistical analyses. Principal Components Analysis (PCA) was conducted after scaling the variables to have a mean of zero and a standard deviation of one, to show the correlation between environmental variables and sampling sites throughout the catchment. Classical diversity indices (species richness and Shannon) were calculated using the diversity() function of the vegan package. We used the Hellinger transformed abundance matrix for species composition and OTU reads (Legendre and Gallagher, 2001). PERMANOVA was conducted to determine which environmental variables influenced diatoms and microalgae OTU assemblages. Furthermore, to characterise the influence of environmental variables on diatom species composition and microalgae OTU assemblages, we conducted the Redundancy Analysis (RDA). To simplify the RDA model, we performed a stepwise selection of statistically important environmental variables using ‘ordiR2step()’ function. Only variables with significant correlations were visualized in the graphs.
3 Results
3.1 Environmental variables analysis
The physical and chemical parameters measured in this study are summarized in Table 1. The pH value of the water was neutral to relatively alkaline (7.03 to 9.55), with low to moderate electrical conductivity (108 to 862 µS/cm) (More details are given in Supplementary Table 1). Water temperature varied between 6 and 12.9°C, mean 8.69°C. The PCA performed on the matrix of physical and chemical variables showed that the first two axes explained 74.6% of the total variation in the data, so the first axis explained 58.2% of the variation and the second axis explained 16.4% of the variation (Figure 2). The most significant gradients in the environmental data were represented by chloride, conductivity, sulfate, nitrate and total nitrogen, which were the main contributors to the first PCA axis, while total phosphate and average depth were the main contributor to the second PCA axis (Figure 2). The PCA ordination diagram showed two main groups of sampling locations. On the left side, the locations are distributed along a hydrological gradient (average depth vs. average velocity) that can be interpreted to roughly reflect altitude or upstream – downstream position in the catchment. The right side of the diagram shows three downstream sampling locations (FAL1, KRE1 and KRE2) offset from the others by increased nutrients and salinity (Figure 2; Table 1). This means that while our sampling locations seem to capture well the catchment-wide hydrological gradient, it does not resolve the nutrient and salinization gradient that appears close to orthogonal to the latter.
Figure 2. PCA biplot showing the relationship between the physical and chemical parameters and sampling sites. The arrows and colour gradient show the contribution of different variables to the principal component.
3.2 Photosynthetic biomass
We found statistically significant differences between photosynthetic pigment fluorescence attributed to mean concentrations of cyanobacteria, green algae, diatoms and total algae between different streams in the catchment (ANOVA, all p values < 0.05). The mean biomass of cyanobacteria was lower in Fallbach and higher in Steinau. The mean green algae biomass was lower in Orb and higher in Steinau. The mean biomass of diatoms was lower in Krebsbach and higher in Bracht. Finally, the total algae concentration was lower in Krebsbach and higher in Steinau. The correlation matrix of the physical and chemical environmental variables and mean photosynthetic biomass of different microalgae groups showed that pH was positively correlated with total and cyanobacterial fluorescence; sulfate strongly negatively correlated with total, cyanobacterial, and diatom (but not with green algal) chlorophyll signal. Total and green algal pigment fluorescence was also negatively correlated with average depth, as well as cyanobacterial and green algal signal with average velocity. Besides, water temperature and ammonia both showed a negative correlation with the diatom signal (Figure 3).
Figure 3. Correlation matrix showing the relationship between physical and chemical variables, and photosynthetic biomass of different microalgal groups. The variables in the matrix are ordered such that the correlated ones are next to each other. Positive correlations are displayed in blue and negative correlations in red colour. Colour intensity and the size of the circle are proportional to the correlation coefficients. Correlation coefficients with a p-value > 0.05 are considered insignificant and their fields remain blank. T, Total microalgae concentration; C, Cyanobacteria; G, Green algae; D, Diatoms.
3.3 Diatom community composition: morphological analysis
A total of 142 diatom taxa belonging to 46 genera were identified by digital microscopy. Diatom community composition was dominated by Navicula lanceolate (NLAN), Navicula gregaria (NGRE), Nitzschia dissipata (NDIS), Achnanthidium minutissimum (ADMI) and Fragilaria deformis (FDEF) (Figure 4; Supplementary Table 2). Species richness of individual samples ranged from 21 to 56 (mean 37.3) and Shannon diversity index from 1.26 to 3.25 (mean 2.56). Both species richness and Shannon were lower in KRE1 and higher in GRU2. PERMANOVA common model including all predictors showed that conductivity, pH, and water temperature had significant effects (p < 0.05) on diatom community composition. All other factors had no statistically significant effects on diatom community composition (p > 0.05, Table 2).
Figure 4. Relative abundance of dominant diatom taxa identified by digital light microscopy. Four letter code based on the OMNIDIA short code, explained in Supplementary Table 2.
Table 2. Results of PERMANOVA partitioning testing the effects of environmental variables on microalgae communities.
The RDA model explained a variation of 27% on the first axis and 17.8% on the second. The environmental factors pointed out by this analysis to show the most statistically significant association (p < 0.05) with the differentiation of diatom communities were pH, water temperature, oxygen saturation (oxygen_100%), conductivity, average velocity, average depth, ammonia, chloride and sulfate. Of these, only pH, water temperature and conductivity showed a strong effect, while the effects of all the other factors was weak. The RDA ordination diagram shows sites in the upper left quadrant that are largely influenced by pH and oxygen saturation, and were associated with Nitzschia soratensis (NSTS) and Fragilaria rinoi (FRIO). Sites plotted in the lower right quadrant, were mainly influenced by ammonia, sulfate, chloride, conductivity and average depth. These sites were associated with taxa such as Surirella brebissonii var. kuetzingii (SBKU) and Navicula gregaria (NGRE) (Figure 5). The ecological status index of diatoms indicated by the IPS ranged from 10.5 to 15.9 (mean 14.3), indicating poor to good ecological status. Only one site, KRE2 had poor ecological status (Figure 5). The main component of the environmental gradient correlated with community composition was that of water temperature vs. pH and oxygen saturation, while the second axis was most strongly correlated with water depth, salinity components (chloride, sulfate and conductivity), and ammonia concentration. Interestingly, the magnitude of differentiation in diatom community composition appeared lower than the separation in environmental variables between the three high-nutrient, high-sulfate locations (FAL1, KRE1, KRE2) observed in the above PCA.
Figure 5. Redundancy analysis (RDA) plot presenting the relationships between the analysed diatom assemblages in the samples based on digital light microscopy and environmental variables. ADMI, Achnanthidium minutissimum; FDEF, Fragilaria deformis; FRIO, Fragilaria rinoi; GCLF, Gomphonema calcifugum; GOLI, Gomphonema olivaceum; NDIS, Nitzschia dissipata var. dissipata; NDME, Nitzschia dissipata var. media; NLAN,Navicula lanceolata; NSOC, Nitzschia sociabilis; NSTS, Nitzschia soratensis; SBKU, Surirella brebissonii var. kuetzingii.
3.4 Microalgae community composition: 18S-V9 amplicon sequencing
A total of 21,293,230 reads were obtained after running the Natrix Bioinformatics pipeline on 18S-V9 amplicon sequencing data from the 19 sampling sites (average 1,120,696 reads per sample, minimum 533,819 reads and maximum 1,604,630 reads). These reads were clustered into 7294 OTUs from various protist groups and assigned taxonomy up to the highest possible level. Microalgae OTU richness ranged from 803 to 2295 (mean 1460), Shannon diversity from 6.69 to 7.74 (mean 7.24). The lowest alpha diversity values were recovered in STE1 and the highest in GRU2. Diatoms contributed most of the MPB OTUs, making up the majority of reads (13,322,277 or about 79.3% of total read counts; Figure 6A; Supplementary Table 2). We found 1400 different diatom OTUs, with OTU richness ranging from 140 to 392 (mean 247.59), Shannon diversity ranging from 1.52 to 2.66 (mean 2.16). The lowest value of diatom OTU richness was recorded in BIE2 and the highest in FAL1. The diatom OTUs with the highest read counts were assigned to Navicula radiosa (NRAD), Navicula gregaria (NGRE), Nitzschia dissipata (NDIS), Surirella brebissoni (SBRE), Melosira varians (MVAR) and Asterionella formosa (AFOR) (Figure 6B). Many of these taxa were also observed during the digital light microscopy analysis, except for N. radiosa and A. formosa (Supplementary Table 2).
Figure 6. Relative reads abundance of (A) dominant OTU classes in the samples and (B) taxa associated with dominant diatom (Bacillariophyta) OTUs. Four letter code based on the OMNIDIA short code, explained in Supplementary Table 2.
PERMANOVA common model including all predictors run on the full microalgae OTU dataset revealed that water pH, temperature, conductivity had significant effect on community composition (p < 0.05, Table 2). Total phosphate had a weak effect on microalgae community. Similarly, the same model run on the diatom subset of OTU reads revealed significant effects of water pH, temperature, conductivity, discharge and total phosphate (p < 0.05, Table 2).
The RDA model examining the relationship between the diatom subset of OTU reads and environmental variables showed a lower explanatory power than the one for the diatom microscopy data set, with a fitted variation of 20.5% explained by the first axis and 17% by the second axis (Figure 7A). The major explanatory gradient here was a salinization – depth gradient that played a slightly lower role in explaining the diatom microscopy results, with the pH vs. water temperature gradient occurring here along the second axis (Figure 7A). In the upper left quadrant of the RDA ordination diagram, sites such as OBR1 and OBR2 were influenced by flow velocity, while pH was the most important predictor at sites such as KIN2 and KIN3 in the lower left quadrant. Sites BIE1, BIE2 and GRU1 in the upper right quadrant were mostly associated with higher water temperature. Finally, FAL1, KIN1 and KIN2 in the lower right quadrant were mainly influenced by chloride, conductivity and depth. The ecological status index (IPS) calculated from the diatom amplicon data set was between 11.4 and 18.1 (mean value 15.6), indicating poor to very good water quality (one has to note that this index was not calibrated for amplicon data sets). Only site KRE2 had a poor ecological status (Figure 7A). In this RDA, the high nutrient, salinized set of samples (FAL1, KRE1, KRE2) appeared with the highest RDA1 scores, but not so strongly offset from the rest of the samples as in the environmental PCA.
Figure 7. Redundancy analysis (RDA) plots presenting the relationships between the analysed OTUs assemblages and environmental variables: (A) diatom (Bacillariophyta) subset of OTUs (ecological status according to IPS index) and (B) the whole microalgae OTU dataset (small blue dots show significant OTUs). AFOR, Asterionella formosa; EOMI, Eolimna minima; GAFF, Gomphonema affine var. affine; GMIC, Gomphonema micropus var. micropus; MVAR, Melosira varians; NGRE, Navicula gregaria; NDIS, Nitzschia dissipata var. dissipata; SBRE, Surirella brebissonii; SULN, Synedra ulna.
Similarly, the RDA model explained a fitted variation of 19.8% on the first axis and 17.4% on the second axis when analysing the relationship between the entire microalgae OTU read dataset and physical and chemical variables (Figure 7B). The main explanatory environmental gradients were highly similar to the diatom OTU data set (with opposite signs along the first axis). Major environmental factors contributing to these variations (p < 0.05) included pH, temperature, conductivity, average velocity, average depth, total nitrogen and chloride. The separation of FAL1, KRE1 and KRE2 from the rest of the samples was clearer in this case than for both above diatom community data sets. (Figure 7B). Sites located in the Kinzig Stream (KIN1 – KIN4), plotted in the lower left quadrant, were mainly influenced by pH, while sites plotted in the lower right quadrant (BRA1, BRA2, SAL1, SAL2, STE1 and STE2) were associated with flow velocity (Figure 7B).
4 Discussion
4.1 Microalgae community composition in response to environmental variables
The effort to understand the effects of anthropogenic stressors on microalgal communities within a river catchment remains at the forefront of research in aquatic ecosystems. The compelling reason for this increasing interest is the widespread recognition that anthropogenic activities may have differential impacts on freshwater ecosystems. This study aimed to identify the main environmental drivers of variability in photosynthetic biomass and microalgal community compositions in the anthropogenically stressed Kinzig catchment. Alongside a set of candidate environmental predictors, we used for this purpose data sets capturing total and group-specific photosynthetic biomass and community composition of diatoms (characterized by microscopy vs. by amplicon sequencing) as well as of all microalgal groups (amplicon sequencing).
The PCA ordination of environmental predictors highlighted a hydrological gradient (depth and discharge vs. velocity) well resolved by our samples. Close to orthogonal to that, the PCA indicated a water chemistry gradient (nutrient and salts, the latter mainly sulphate, but also chloride) captured by two clusters of sites (the high nutrient, high salinity locations FAL1, KRE1 and KRE2 vs. all others). In addition, total phosphate concentration showed an aspect of variability somewhat decoupled from both of these gradients (Figure 2).
In terms of photosynthetic biomass (as estimated using a BenthoTorch), the depth gradient was slightly, but significantly negatively associated with total and green algal chlorophyll concentration (Figure 3), indicating a slightly decreasing trend in green algal biomass content with increasing depth which could be due to light limitations (Virtanen and Soininen, 2012; Sabater et al., 2016). Total, cyanobacterial and diatom chlorophyll signals were strongly negatively correlated with sulfate, also being the main driver of the salinity gradient in the system. Thus, whereas biomass of green algae seems to respond more strongly along the hydrological gradient, complementary to this, the chlorophyll biomass of diatoms and cyanobacteria appears to be driven mainly by the negative effect of sulfate (presumably more due to its osmotic effect than to specific toxicity). Conductivity, widely acknowledged as a key factor influencing periphytic communities, particularly diatoms, is underscored by several studies (Soininen et al., 2004; Virtanen and Soininen, 2012; Stenger-Kovács et al., 2020). It mainly reflects the overall ionic concentrations and closely correlates with water pH, another important factor noted in this study, as besides these two main environmental gradients, pH was positively correlated with the total and cyanobacterial biomass estimates; and diatom biomass estimates negatively with water temperature (as well as ammonium that can be seen as a component of the above water chemistry gradient).
We used two statistical approaches (PERMANOVA and RDA) to identify the main environmental drivers of community variability in three data sets: diatoms as assessed by microscopy; diatoms assessed by amplicon sequencing; and all microalgal taxa by amplicon sequencing. Common to all PERMANOVA and RDA results for all three data sets was the identification of a strong role of conductivity in driving patterns of community variability in both diatoms and all microalgal taxa (Table 2; Figures 5, 7), reflecting the main water chemistry gradient described. Both PERMANOVA and RDA highlighted a strong response of diatom and microalgae communities to pH and water temperature besides conductivity in the case of both the amplicon and the microscopy data sets.
In the case of the microalgal amplicon data, discharge, total phosphate concentrations showed significant explanatory power of community composition besides conductivity, water temperature and pH, indicating possible taxon-specific differences in ecological responses of different microalgal groups. Cyanobacteria, which were found to respond to some environmental gradient in terms of their pigment signatures, were not captured by our community data sets at all, leaving their specific responses open.
In addition, RDA showed that flow velocity was a significant factor for microalgae community composition. This is consistent with earlier research, several studies have reported that stream flow velocity is a major factor influencing the heterogeneity in the composition of microphytobenthic communities and biodiversity in streams (Boulêtreau et al., 2010; Townsend et al., 2012; Breuer et al., 2017). Benthic microalgae are directly exposed to increased flow velocity and material movement carried by water current. Furthermore, gas bubbles produced through metabolic processes like photosynthesis or respiration can further disrupt and upset the stability or integrity of the microphytobenthic communities under high flow velocity conditions (Boulêtreau et al., 2010). The microalgae communities tend to be highly specialized at high flow velocity (Soininen, 2005). For example, attached taxa (like Achnantidium minutissimum) can tolerate higher flow velocity, it is also among the first colonizers. Flow velocity influences the availability of other resources, such as light, suspended particles and nutrients, as well as other abiotic factors (e.g. temperature) (Townsend et al., 2012; Sabater et al., 2016; Breuer et al., 2017), thus having a significant impact on microalgae community composition.
Judged qualitatively by the spread of points in the RDA, as well as by the number of significant environmental correlations revealed, or also by the number of OTUs being almost an order of magnitude higher than the number of taxa observed microscopically, the amplicon data set provided a better resolution of the ecological signals in diatom communities than the microscopy data set (Zimmermann et al., 2015). Despite the lower diatom taxonomic resolution of the 18S-V9 marker compared to other diatom specific markers like the rbcl and 18S-V4, we found more diatom OTUs (1400) than morphospecies identified by morphology (142). While some OTUs may represent the same species, the difference in taxonomic resolution between the two methods appears to be minimal. Since closely related species often have identical 18S rDNA sequences, it could be expected that even this signal could possibly be improved by the use of more variable markers (Zimmermann et al., 2014). Furthermore, morphologically defined species do not always correspond to true biological species. Both molecular and morphological methods offer different but imperfect perspectives on the structure of biological communities and should be considered as complementary approaches (Minerovic et al., 2020). By the same criteria, with slightly different main drivers revealed in the PERMANOVA, the ecological resolution of the microalgal amplicon data set seems to have been the highest. This might be expected simply due to the broader taxonomic coverage.
The main gradients explaining community variability in the RDA ordinations were highly similar across the three data sets (with most variation in depicted patterns corresponding to rotation and mirroring of highly similar configurations). Common to all three data sets, these ordinations differed from what might have been expected from the environmental PCA by showing a strong alignment between the effects of the above described main hydrological and water chemistry gradients, with water temperature and pH showing strong effects, close to orthogonal to this main axis of variability. This is in line with the PERMANOVA results, and could indicate that less prominent gradients (as identified by a PCA) can also have strong ecological effects. As to the alignment of the main hydrological and water chemistry axes of variability, it is possible that more extensive sampling could differentiate between their effects upon community composition more clearly, this aspect will need to be revisited in the future. These findings agree with other studies in which water physical and chemical variables were main drivers of changes in microalgae assemblages, often with shift from more sensitive species to more tolerant species composition (Wu et al., 2011; Soininen et al., 2004; Virtanen and Soininen, 2012; Teittinen et al., 2015; Breuer et al., 2017).
The community ordinations (RDAs) also delivered somewhat unexpected results with respect to the FAL1, KRE1, KRE2 group of samples. These appeared with a substantial offset from all other samples in the environmental PCA, indicating a large difference in terms of both nutrient and salt concentrations. Although they appeared in relatively close vicinity of each other also in the RDA plots, their outlier position compared to all other samples was less pronounced than in the environmental PCA. Also in this respect, one can observe a pattern of increasing resolution in the rank order of diatom microscopy; diatom amplicon; microalgal amplicon data set.
Our results are consistent with (Nguyen et al., 2023), who found that the macroinvertebrate communities in the Kinzig catchment followed an environmental gradient ranging from sites with good, moderate to poor water quality. For example, the majority of species in the macroinvertebrate communities found at good water quality sites were grazers that preferred clear water with coarser substrates, and oligotrophic nutrient content. On the other hand, species that were foragers or active filter feeders and had a preference for turbid water with finer substrates, brackish and eutrophic/mesotrophic conditions were more likely to be found in poor water quality sites (Nguyen et al., 2023).
Our IPS results from microscopy analysis showed poor to good ecological status, while amplicon analysis showed poor to high (very good) ecological status at the sampled sites throughout the catchment (Figures 5, 7A). This could be explained by the discrepancies in community composition observed between the two methods. These discrepancies may be due to biases in both methods. For example, microscopy may ignore some barely visible taxa, while amplicon sequencing still miss many taxa in the PR2 reference library used in the present study (Zimmermann et al., 2014). Diatom ecological status indices such as the IPS are calculated on the notion that the ecological quality of a site is represented by the average ecological value of its taxa, weighted by their relative abundances (Stevenson et al., 2010). As a result, high abundance of a taxa strongly influences these indices, whereas low abundance of a taxa only slightly affects the final value (Tapolczai et al., 2024). For example, when our microscopic observations revealed that Navicula lanceolata was the most dominant taxa in most of the samples, amplicon analysis showed that the most abundant OTUs were associated with Navicula radiosa. And, the latter had a higher sensitivity to the IPS (S = 5.0) than the former (S = 3.8). The variations in water quality classes between the two methods may also be due to overestimation of the relative reads abundance of large diatom taxa in amplicon sequencing analysis. For example, larger diatoms like Navicula lanceolata and N. radiosa have more gene copies than smaller ones like Achnanthidium minutissimum (Vasselon et al., 2018; Minerovic et al., 2020). However, our results were consistent in the case of site KRE2, which had a poor ecological status as demonstrated by both microscopy and amplicon analysis. The substrate at this site was largely sediment and it is one of the three sites discriminated by the PCA, located at the lowest altitude, which followed a nutrients and salinity gradient.
4.2 Taxonomic inconsistencies between diatom microscopy and amplicon sequencing
Our results showed that microalgae communities characterized by different methods used in this study were influenced by various environmental factors, which agrees with previous studies (Virtanen and Soininen, 2012; Teittinen et al., 2015; Stenger-Kovács et al., 2020). Digital microscopy analysis revealed that Navicula lanceolata was the most dominant diatom taxon in the majority of our sampling sites, except KIN3 and KIN4. This taxon was reported to be the most dominant in an anthropogenically modified river in Poland, reaching over 80% of the total number of counted valve in spring (Noga et al., 2013). The authors attributed this dominance to the sampling season rather than anthropogenic influences. It was also identified as one of the taxa that had a higher proportion in agricultural land, with a higher content of nutrients and ions, such as in the Carpathian Basin (Stenger-Kovács et al., 2020). It occurred massively in slightly eutrophic boreal streams (Virtanen and Soininen, 2012). Due to its broad trophic preference for electrolyte-rich water and cooler temperatures, N. lanceolata grows rapidly during the winter and spring months (Noga et al., 2013; Schröder et al., 2015; Cantonati et al., 2017). This may explain why this species was highly represented in our samples collected in the first weeks of spring 2021. Therefore, the dominance of this taxon in our samples could also be related to the sampling season.
In contrast, amplicon sequencing analysis showed that Navicula radiosa was the most dominant taxon in most of our samples. We did not identify this taxon in microscopic analysis, instead we found N. lanceolata to be the most abundant taxon. Given the common taxonomic grouping of these two diatom species in the genus Navicula, they certainly have some similar morphological features that could lead to confusion in identification. However, N. lanceolata has obtusely rounded ends and a roundish central area (Cantonati et al., 2017). The observed discrepancies in our results are likely due to misidentification in the reference database, suggesting that the name N. radiosa may have been incorrectly assigned to OTUs belonging to N. lanceolata (Zimmermann et al., 2014). Other taxa detected in high abundances through molecular analysis in this study, include Navicula gregaria and Surirella brebisonii. These two taxa were also identified in light microscopy analysis. They grow abundantly in various environmental conditions, from mesotrophic to eutrophic waters. For example, N. gregaria has been reported as one of the indicators of nutrient enrichment in boreal alkaline and eutrophic basins (Virtanen and Soininen, 2012). We suspect that the presence of this taxon in large quantities in our samples may be related to nutrients and salt concentration in the rive catchment.
We found 1400 different diatom OTUs, about ten times the number of taxa we detected with morphological observations (142), demonstrating the broader taxonomic coverage of amplicon sequencing analysis. Ribosomal DNA markers can be used to amplify several groups of organisms, including prokaryotes and eukaryotes. Most amplicon sequencing studies of eukaryotic communities have employed 18S markers targeting individual phylogenetic lineages (Minerovic et al., 2020). For example, fungus-specific ITS primers (Blaalid et al., 2013) and diatom-specific 18S-V4 primers (Zimmermann et al., 2015) are widely used, and reference libraries have been created to support their application to metabarcoding (Zimmermann et al., 2014). The V9 region was chosen in the present study due to its potential to simultaneously characterize multiple groups of eukaryotic microalgal organisms (Minerovic et al., 2020). This genetic marker cannot detect cyanobacteria because it targets the 18S rRNA gene, which is specific to eukaryotic organisms. Cyanobacteria are prokaryotes (bacteria) and therefore possess the 16S rRNA gene instead of the 18S rRNA gene (Salis et al., 2017). Thus, 16S markers should be used for amplicon sequencing of cyanobacteria. Furthermore, the reference sequences from the 18S-V9 region for many freshwater diatom taxa are still missing from the PR2 reference database used in the present study. This shows that a lot of work is still needed to complete the reference databases in this research area. Even the well-curated and updated diatom reference database as Diat.barcode (Rimet et al., 2019) is still incomplete, which can lead to both non-identification and misidentification of taxa (Zimmermann et al., 2014; Bailet et al., 2019). In addition, it is not clear whether the issues of cryptic and pseudocryptic species and infraspecific or intragenomic variation are covered by rDNA reference databases (Minerovic et al., 2020).
Our results showed that only a small number of diatom taxa could be identified through both morphological analysis and amplicon sequencing to the lowest taxonomic level. The reference database used for taxonomic assignment in our amplicon analysis is likely incomplete, which could explain the large differences in taxonomic composition between the two methods. Since diatoms are unicellular organisms, variables such as the number of genomes and the number of gene copies per genome influence the number of copies of a gene (Vasselon et al., 2018). This implies that gene copy number could be approximated by cell biovolume. From the smallest to the largest diatom species, there is variation in gene copy number that can have a major effect on metabarcoding quantification (Vasselon et al., 2018). This could explain why highly abundant small-size taxa observed in microscopy such as Achnanthidium minutissimum were not as dominant in amplicon analysis.
In conclusion, our study aimed to identify the key environmental variables that influence microalgal biomass and community composition in the Kinzig catchment. Our results showed that depth gradient slightly negatively correlated with total and green algal biomass. Sulfate showed a strong negative correlation with the total, cyanobacterial, and diatom biomass. Total and cyanobacterial biomass were positively associated with pH, while diatom biomass was negatively correlated with ammonium content and water temperature. PERMANOVA showed that conductivity, water temperature and pH were the most important factors influencing microalgae community composition, as observed in both microscopy and amplicon analysis. Together with these three variables, total phosphate in all microalgae OTU and water discharge in diatom (Bacillariophyta) OTU datasets may suggest taxon-specific variations in the ecological responses of various microalgae groups. Our results highlighted the sensitivity of eukaryotic algae to physical and chemical variables, demonstrating their potential as good indicators of environmental gradients in the Kinzig catchment. These results suggest the complex relationship between multiple environmental factors and microalgal biomass and community composition. Further investigations, including collection of time series data, are required to fully understand the underlying biotic and abiotic factors driving these microalgal communities.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
NM: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. MD: Data curation, Writing – review & editing, Formal analysis. MK: Investigation, Methodology, Software, Writing – review & editing. DV: Data curation, Writing – review & editing, Formal analysis. DB: Investigation, Methodology, Software, Writing – review & editing. AB-C: Investigation, Methodology, Supervision, Validation, Writing – review & editing, Formal analysis. BB: Conceptualization, Funding acquisition, Resources, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was performed within the framework of the Collaborative Research Centre (CRC) 1439 RESIST (Multilevel Response to Stressor Increase and Decrease in Stream Ecosystems; www.sfb-resist.de) funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation; CRC 1439/1, project number: 426547801).
Acknowledgments
We want to thank Marzena Spyra and all the students who helped in the field with sample collection and in the laboratory with sample processing. The editor and two anonymous reviewers are thanked for their comments that helped improve this manuscript. We acknowledge support by the Open Access Publication Fund of the University of Duisburg-Essen.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2024.1422925/full#supplementary-material.
Supplementary Table 1 | Physical and chemical parameters analysed in this study.
Supplementary Table 2 | List of taxa identified by microscopy (_MIC) or amplicon sequencing (_DNA) with their 4 letters OMINIDIA codes. The results are in % of valves for MIC and % of reads for DNA.
References
Amaral-Zettler L. A., McCliment E. A., Ducklow H. W., Huse S. M. (2009). A method for studying protistan diversity using massively parallel sequencing of V9 hypervariable regions of small-subunit ribosomal RNA genes. PloS One 4, e6372. doi: 10.1371/journal.pone.0006372
Bailet B., Bouchez A., Franc A., Frigerio J.-M., Keck F., Karjalainen S.-M., et al. (2019). Molecular versus morphological data for benthic diatoms biomonitoring in Northern Europe freshwater and consequences for ecological status. MBMG 3, 21–35. doi: 10.3897/mbmg.3.34002
B-Béres V., Stenger-Kovács C., Buczkó K., Padisák J., Selmeczy G. B., Lengyel E., et al. (2023). Ecosystem services provided by freshwater and marine diatoms. Hydrobiologia 850, 2707–2733. doi: 10.1007/s10750-022-04984-9
Bey M.-Y., Ector L. (2013). Atlas des Diatomées des Cours d’eau de la Région Rhône-Alpes (Préfet de la Région Auvergne-Rhône-Alpes: Lyon, France: Direction régionale de l'environnement, de l'aménagement et du logement).
Birk S., Chapman D., Carvalho L., Spears B. M., Andersen H. E., Argillier C., et al. (2020). Impacts of multiple stressors on freshwater biota across spatial scales and ecosystems. Nat. Ecol. Evol. 4, 1060–1068. doi: 10.1038/s41559-020-1216-4
Blaalid R., Kumar S., Nilsson R. H., Abarenkov K., Kirk P. M., Kauserud H. (2013). ITS1 versus ITS2 as DNA metabarcodes for fungi. Mol. Ecol. Resour. 13, 218–224. doi: 10.1111/1755-0998.12065
Boulêtreau S., Sellali M., Elosegi A., Nicaise Y., Bercovitz Y., Moulin F., et al. (2010). Temporal dynamics of river biofilm in constant flows: A case study in a riverside laboratory flume. Int. Rev. Hydrobiology 95, 156–170. doi: 10.1002/iroh.200911203
Breuer F., Janz P., Farrelly E., Ebke K.-P. (2017). Environmental and structural factors influencing algal communities in small streams and ditches in central Germany. J. Freshw. Ecol. 32, 65–83. doi: 10.1080/02705060.2016.1241954
Buchner D. (2022a). Sample preparation and lysis of homogenized malaise trap samples v1. doi: 10.17504/protocols.io.dm6gpjrmjgzp/v1 (Accessed June 25, 2023)
Buchner D. (2022b). Guanidine-based DNA extraction with silica-coated beads or silica spin columns v2. doi: 10.17504/protocols.io.eq2ly73mmlx9/v2 (Accessed June 25, 2023)
Burfeid-Castellanos A. M., Kloster M., Beszteri S., Postel U., Spyra M., Zurowietz M., et al. (2022). A digital light microscopic method for diatom surveys using embedded acid-cleaned samples. Water 14, 3332. doi: 10.3390/w14203332
Cantonati M., Kelly M. G., Lange-Bertalot H. (Eds.) (2017). Freshwater benthic diatoms of Central Europe: Over 800 common species used in ecological assessment; engl. ed. with updated taxonomy and added species (Schmitten-Oberreifenberg: Koeltz Botanical Books).
Castro-Català N., Dolédec S., Kalogianni E., Skoulikidis N. T., Paunovic M., Vasiljević B., et al. (2020). Unravelling the effects of multiple stressors on diatom and macroinvertebrate communities in European river basins using structural and functional approaches. Sci. Total Environ. 742, 140543. doi: 10.1016/j.scitotenv.2020.140543
CEN. UNE-EN (2014). Water quality—Guidance for the routine sampling and preparation of benthic diatoms from rivers and lakes Vol. 13946 (Brussels, Belgium: CEN-CENELEC Management Centre), 2014.
Cox E. J. (1995). Studies on the diatiom genus Navicula Bory. VII. The identity and typification of Navicula gregaria Donkin, N. Cryptocephala Kütz. and related taxa. Diatom Res. 10, 91–111. doi: 10.1080/0269249X.1995.9705330
Dalu T., Magoro M. L., Naidoo L. S., Wasserman R. J., Human L. R., Adams J. B., et al. (2020). Microphytobenthos diversity and community structure across different micro-estuaries and micro-outlets: Effects of environmental variables on community structure. Environ. pollut. 260, 114097. doi: 10.1016/j.envpol.2020.114097
Deep A., Bludau D., Welzel M., Clemens S., Heider D., Boenigk J., et al. (2023). Natrix2 – Improved amplicon workflow with novel Oxford Nanopore Technologies support and enhancements in clustering, classification and taxonomic databases. MBMG 7, 263–271. doi: 10.3897/mbmg.7.109389
Dixit S. S., Smol J. P., Kingston J. C., Charles D. F. (1992). Diatoms: powerful indicators of environmental change. Environ. Sci. Technol. 26, 22–33. doi: 10.1021/es00025a002
European Commission (2000). Directive 2000/60/EC of the European parliament and of the council of 23 October 2000 establishing a framework for community action in the field of water policy. Off. J. Eur. Communities. L327, 1-72.22.12.2000. Available online at: https://eur-lex.europa.eu/eli/dir/2000/60/oj (Accessed August 5, 2021).
Feld C. K., Birk S., Eme D., Gerisch M., Hering D., Kernan M., et al. (2016). Disentangling the effects of land use and geo-climatic factors on diversity in European freshwater ecosystems. Ecol. Indic. 60, 71–83. doi: 10.1016/j.ecolind.2015.06.024
Frøslev T. G., Kjøller R., Bruun H. H., Ejrnæs R., Brunbjerg A. K., Pietroni C., et al. (2017). Algorithm for post-clustering curation of DNA amplicon data yields reliable biodiversity estimates. Nat. Commun. 8, 1188. doi: 10.1038/s41467-017-01312-x
Fu L., Niu B., Zhu Z., Wu S., Li W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Guillou L., Bachar D., Audic S., Bass D., Berney C., Bittner L., et al. (2013). The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D597–D604. doi: 10.1093/nar/gks1160
Hlúbiková D., Novais M. H., Dohet A., Hoffmann L., Ector L. (2014). Effect of riparian vegetation on diatom assemblages in headwater streams under different land uses. Sci. Total Environ. 475, 234–247. doi: 10.1016/j.scitotenv.2013.06.004
Kahlert M., McKie B. G. (2014). Comparing new and conventional methods to estimate benthic algal biomass and composition in freshwaters. Environ. Sci. Process Impacts 16, 2627–2634. doi: 10.1039/C4EM00326H
Kassambara A., Mundt F. (2020). factoextra: Extract and visualise the reults of multivariate data analyses. Available online at: https://doi.org/10.32614/CRAN.package.factoextra (Accessed September 9, 2023)
Kelly M. G., Cazaubon A., Coring E., Dell'Uomo A., Ector L., Goldsmith B., et al. (1998). Recommendations for the routine sampling of diatoms for water quality assessments in Europe. J. Appl. Phycol 10, 215–224. doi: 10.1023/A:1008033201227
Lange A., Jost S., Heider D., Bock C., Budeus B., Schilling E., et al. (2015). AmpliconDuo: A split-sample filtering protocol for high-throughput amplicon sequencing of microbial communities. PloS One 10, e0141590. doi: 10.1371/journal.pone.0141590
Langenkämper D., Zurowietz M., Schoening T., Nattkemper T. W. (2017). BIIGLE 2.0 - browsing and annotating large marine image collections. Front. Mar. Sci. 4. doi: 10.3389/fmars.2017.00083
Lavoie I., Morin S., Laderriere V., Fortin C. (2018). Freshwater diatoms as indicators of combined long-term mining and urban stressors in junction creek (Ontario, Canada). Environments 5, 30. doi: 10.3390/environments5020030
Lecointe C., Coste M., Prygiel J. (1993). Omnidia": software for taxonomy, calculation of diatom indices and inventories management. Hydrobiologia 269-270, 509–513. doi: 10.1007/BF00028048
Legendre P., Gallagher E. D. (2001). Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280. doi: 10.1007/s004420100716
Lemley D. A., Adams J. B., Bate G. C. (2016). A review of microalgae as indicators in South African estuaries. South Afr. J. Bot. 107, 12–20. doi: 10.1016/j.sajb.2016.04.008
Levkov Z., Metzeltin D., Pavlov A. (Eds.) (2013). “Luticola and luticolopsis,” in Diatoms of the european inland waters and comparable habitats (Koeltz Botanical Books, Oberreifenberg, Germany), 698.
Mahé F., Rognes T., Quince C., de Vargas C., Dunthorn M. (2015). Swarm v2: highly-scalable and high-resolution amplicon clustering. PeerJ 3, e1420. doi: 10.7717/peerj.1420
Martin M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10. doi: 10.14806/ej.17.1.200
Masella A. P., Bartram A. K., Truszkowski J. M., Brown D. G., Neufeld J. D. (2012). PANDAseq: paired-end assembler for illumina sequences. BMC Bioinf. 13, 31. doi: 10.1186/1471-2105-13-31
Minerovic A. D., Potapova M. G., Sales C. M., Price J. R., Enache M. D. (2020). 18S-V9 DNA metabarcoding detects the effect of water-quality impairment on stream biofilm eukaryotic assemblages. Ecol. Indic. 113, 106225. doi: 10.1016/j.ecolind.2020.106225
Morales E. A., Siver P. A., Trainor F. R. (2001). Identification of diatoms (Bacillariophyceae) during ecological assessments: Comparison between Light Microscopy and Scanning Electron Microscopy techniques. ansp 151, 95–103. doi: 10.1635/0097-3157(2001)151[0095:IODBDE]2.0.CO;2
Nguyen H. H., Kiesel J., Peters K., Hering D., Sinclair J. S., Haase P. (2023). Stream macroinvertebrate community metrics consistently respond to a spatiotemporal disturbance gradient but composition is more context-dependent. Landsc Ecol. 38, 3133–3151. doi: 10.1007/s10980-023-01769-w
Noga T., Stanek-Tarkowska J., Peszek Ł., Pajączek A., Kowalska S. (2013). Use of diatoms to assess water quality of anthropologically modified Matysówka stream. J. Ecol. Eng. 14, 1–11. doi: 10.5604/2081139X.1043166
Oksanen J., Simpson G., Blanchet F., Kindt R., Legendre P., Minchin P., et al. (2022). vegan: Community Ecology Package. Available online at: https://cran.r-project.org/web/packages/vegan/index.html (Accessed November 2, 2022)
Passy S. I. (2007). Diatom ecological guilds display distinct and predictable behavior along nutrient and disturbance gradients in running waters. Aquat. Botany. 86, 171–178. doi: 10.1016/j.aquabot.2006.09.018
Patil I. (2021). Visualizations with statistical details: The 'ggstatsplot' approach. JOSS 6, 3167. doi: 10.21105/joss.03167
Prygiel J., Coste M. (1993). Utilisation des indices diatomiques pour la mesure de la qualité des eaux du bassin Artois-Picardie: bilan et perspectives. Ann. Limnol. - Int. J. Lim. 29, 255–267. doi: 10.1051/limn/1993021
Rimet F., Gusev E., Kahlert M., Kelly M. G., Kulikovskiy M., Maltsev Y., et al. (2019). Diat.barcode, an open-access curated barcode library for diatoms. Sci. Rep. 9, 15116. doi: 10.1038/s41598-019-51500-6
Rognes T., Flouri T., Nichols B., Quince C., Mahé F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584. doi: 10.7717/peerj.2584
Round F. E. (1991). Diatoms in river water-monitoring studies. J. Apllied Phycology. 3, 129–145. doi: 10.1007/BF00003695
Sabater S., Artigas J., Corcoll N., Proia L., Timoner X., Tornés E. (2016). “Ecophysiology of river algae,” in River algae [recurso electrónico]. Ed. Necchi JR O. (Springer, Cham, Cham), 197–217.
Salis R. K., Bruder A., Piggott J. J., Summerfield T. C., Matthaei C. D. (2017). High-throughput amplicon sequencing and stream benthic bacteria: identifying the best taxonomic level for multiple-stressor research. Sci. Rep. 7, 1–12. doi: 10.1038/srep44657
Schloss P. D., Westcott S. L., Ryabin T., Hall J. R., Hartmann M., Hollister E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schröder M., Sondermann M., Sures B., Hering D. (2015). Effects of salinity gradients on benthic invertebrate and diatom communities in a German lowland river. Ecol. Indic. 57, 236–248. doi: 10.1016/j.ecolind.2015.04.038
Schürings C., Feld C. K., Kail J., Hering D. (2022). Effects of agricultural land use on river biota: a meta-analysis. Environ. Sci. Eur. 34, 1–13. doi: 10.1186/s12302-022-00706-z
Smol J. P., Stoermer E. F. (2010). The diatoms: applications for the environmental and earth sciences. Cambridge University Press, Cambridge.
Soininen J., Jamoneau A., Rosebery J., Passy S. I. (2016). Global patterns and drivers of species and trait composition in diatoms. Global Ecol. Biogeogr. 25, 940–950. doi: 10.1111/geb.12452
Soininen J., Paavola R., Muotka T. (2004). Benthic diatom communities in boreal streams: community structure in relation to environmental and spatial gradients. Ecography 27, 330–342. doi: 10.1111/j.0906-7590.2004.03749.x
Spaulding S. A., Potapova M. G., Bishop I. W., Lee S. S., Gasperak T. S., Jovanoska E., et al. (2022). Diatoms.org: supporting taxonomists, connecting communities. Diatom Res. 36, 291–304. doi: 10.1080/0269249X.2021.2006790
Stenger-Kovács C., Lengyel E., Sebestyén V., Szabó B. (2020). Effects of land use on streams: traditional and functional analyses of benthic diatoms. Hydrobiologia 847, 2933–2946. doi: 10.1007/s10750-020-04294-y
Stevenson R. J., Pan Y., van Dam H. (2010). “Assessing environmental conditions in rivers and streams with diatoms,” in The diatoms: Applications for the environment and earth sciences. Eds. Smol J. P., Stoermer E. F. (Cambridge University Press, Cambridge), 57–85.
Stoeck T., Bass D., Nebel M., Christen R., Jones M. D. M., Breiner H.-W., et al. (2010). Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol. Ecol. 19 Suppl 1, 21–31. doi: 10.1111/j.1365-294X.2009.04480.x
Szczepocka E., Nowicka-Krawczyk P., Olszyński R. M., Żelazna-Wieczorek J. (2021). Response of diatom assemblages to the disruption of the running water continuum in urban areas, and its consequences on bioassessment. PeerJ 9, e12457. doi: 10.7717/peerj.12457
Tapolczai K., Chonova T., Fidlerová D., Makovinská J., Mora D., Weigand A., et al. (2024). Molecular metrics to monitor ecological status of large rivers: Implementation of diatom DNA metabarcoding in the Joint Danube Survey 4. Ecol. Indic. 160, 111883. doi: 10.1016/j.ecolind.2024.111883
Taylor J. C., Harding W. R., Archibald C. G. (2007). A methods manual for the collection, preparation and analysis of diatom samples WRC Report TT 281/07. Water Research Commision, Pretoria, South Africa.
Teittinen A., Taka M., Ruth O., Soininen J. (2015). Variation in stream diatom communities in relation to water quality and catchment variables in a boreal, urbanized region. Sci. Total Environ. 530-531, 279–289. doi: 10.1016/j.scitotenv.2015.05.101
Tornés E., Alández-Rodríguez J., Corrochano A., Nolla-Querol P., Trapote M. C., Sabater S. (2022). Impacts of climate change on stream benthic diatoms—a nation-wide perspective of reference conditions. Hydrobiologia 849, 1821–1837. doi: 10.1007/s10750-022-04829-5
Townsend S. A., Garcia E. A., Douglas M. M. (2012). The response of benthic algal biomass to nutrient addition over a range of current speeds in an oligotrophic river. Freshw. Sci. 31, 1233–1243. doi: 10.1899/11-163.1
Trábert Z., Duleba M., Bíró T., Dobosy P., Földi A., Hidas A., et al. (2020). Effect of land use on the benthic diatom community of the danube river in the region of budapest. Water 12, 479. doi: 10.3390/w12020479
Trobajo R., Rovira L., Ector L., Wetzel C. E., KELLY M., Mann D. G. (2013). Morphology and identity of some ecologically important small Nitzschia species. Diatom Res. 28, 37–59. doi: 10.1080/0269249X.2012.734531
Underwood G. J. C. (2010). Microphytobenthos and phytoplankton in the Severn estuary, UK: present situation and possible consequences of a tidal energy barrage. Mar. pollut. Bull. 61, 83–91. doi: 10.1016/j.marpolbul.2009.12.015
van Dam H., Mertens A., Sinkeldam J. (1994). A coded checklist and ecological indicator values of freshwater diatoms from The Netherlands. Netherlands Journal of Aquatic Ecology. 28, 117–133. doi: 10.1007/bf02334251
Vasselon V., Bouchez A., Rimet F., Jacquet S., Trobajo R., Corniquel M., et al. (2018). Avoiding quantification bias in metabarcoding: Application of a cell biovolume correction factor in diatom molecular biomonitoring. Methods Ecol. Evol. 9, 1060–1069. doi: 10.1111/2041-210X.12960
Virtanen L., Soininen J. (2012). The roles of environment and space in shaping stream diatom communities. Eur. J. Phycology 47, 160–168. doi: 10.1080/09670262.2012.682610
Viso R., Blanco S. (2023). River diatoms reflect better past than current environmental conditions. Water 15, 333. doi: 10.3390/w15020333
Walsh G. W., Wepener V. (2009). The influence of land use on water quality and diatom community structures in urban and agriculturally stressed rivers. WSA 35, 579–594. doi: 10.4314/wsa.v35i5.49184
Wei T., Simko V. (2021). R package ’corrplot’: Visualization of a Correlation Matrix. Available online at: https://github.com/taiyun/corrplot (Accessed November 15, 2023).
Wickham H., Averick M., Bryan J., Chang W., McGowan L., François R., et al. (2019). Welcome to the tidyverse. JOSS 4, 1686. doi: 10.21105/joss.01686
Wu N., Schmalz B., Fohrer N. (2011). Distribution of phytoplankton in a German lowland river in relation to environmental factors. J. Plankton Res. 33, 807–820. doi: 10.1093/plankt/fbq139
Zimmermann J., Abarca N., Enke N., Skibbe O., Kusber W.-H., Jahn R. (2014). Taxonomic reference libraries for environmental barcoding: a best practice example from diatom research. PloS One 9, e108793. doi: 10.1371/journal.pone.0108793
Zimmermann J., Glöckner G., Jahn R., Enke N., Gemeinholzer B. (2015). Metabarcoding vs. morphological identification to assess diatom diversity in environmental studies. Mol. Ecol. Resour 15, 526–542. doi: 10.1111/1755-0998.12336
Keywords: microphytobenthos, biofilm, Bacillariophyta, environmental factors, digital microscopy, amplicon sequencing, photosynthetic biomass
Citation: Mayombo NAS, Dani M, Kloster M, Vidaković D, Buchner D, Burfeid-Castellanos AM and Beszteri B (2024) Assessment of microphytobenthos communities in the Kinzig catchment using photosynthesis-related traits, digital light microscopy and 18S-V9 amplicon sequencing. Front. Ecol. Evol. 12:1422925. doi: 10.3389/fevo.2024.1422925
Received: 24 April 2024; Accepted: 22 August 2024;
Published: 11 September 2024.
Edited by:
Manel Leira, University of Santiago de Compostela, SpainReviewed by:
Salomé Almeida, University of Aveiro, PortugalGabor Varbiro, Hungarian Academy of Sciences, Hungary
Copyright © 2024 Mayombo, Dani, Kloster, Vidaković, Buchner, Burfeid-Castellanos and Beszteri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ntambwe Albert Serge Mayombo, ntambwe.mayombo@uni-due.de
†ORCID: Ntambwe Albert Serge Mayombo, orcid.org/0000-0001-6418-1419
Michael Kloster, orcid.org/0000-0001-9244-4925
Danijela Vidakovic, orcid.org/0000-0002-8696-7743
Dominik Buchner, orcid.org/0000-0002-8499-5863
Andrea M. Burfeid-Castellanos, orcid.org/0000-0002-8556-3157
Bánk Beszteri, orcid.org/0000-0002-6852-1588