Vesna Gulin
Barbara Vlaičević
Mirela Sertić Perić1
Fran Rebrina- 1Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
- 2Department of Biology, Josip Juraj Strossmayer University of Osijek, Osijek, Croatia
Tufa-depositing streams provide great microhabitat complexity and are therefore inhabited by various periphytic phagotrophic organisms such as ciliates and amoeboid protists. Recent removal of invasive plant species Ailanthus altissima (Mill.) Swinge from the Skradinski buk tufa barrier (Krka National Park, Croatia) resulted in changes in the barrier hydromorphology including the reactivation (revitalization) of dry streams. The objective of this study was to investigate: (1) the taxonomic and functional response of periphytic ciliates and amoeboid protists to stream revitalization by comparing taxonomic (i.e., abundance, species richness and diversity) and functional (i.e., functional diversity) metrics between revitalized (N) and control sites (C) during 1 and 2-months immersion period; (2) which environmental and (3) periphyton-associated factors shape the taxonomic and functional metrics and to what extent; (4) how duration of immersion affects taxonomic and functional metrics at revitalized sites. Our results showed that taxonomic and functional metrics of ciliates and amoeboid protists responded to the prevailing conditions characteristic of revitalized tufa-depositing streams: changing hydrology (occasional high flow or drought), soil drainage, and extensive inorganic matter, i.e., tufa deposition, although their responses were somewhat different. The two assemblages also showed different responses of taxonomic and functional metrics with respect to immersion duration: while the taxonomic and functional diversity of ciliates at N sites increased with longer immersion, indicating niche diversification, those of amoeboid protists hardly changed with time. Our results suggest that a comprehensive analysis of taxonomic and functional metrics of ciliates and amoeboid protists could be a good proxy for assessing revitalization of tufa-depositing streams. However, the temporal component should always be considered when conducting such studies, as the colonization processes of ciliates and amoeboid protists are quite complex, especially in tufa-depositing streams.
Introduction
Protists are the most diverse and widespread eukaryotes (de Vargas et al., 2015; Adl et al., 2019; Sieber et al., 2020). They perform important ecosystem functions and play a fundamental part in aquatic food webs (Weisse et al., 2016). Phagotrophic protists (those that utilize food through phagocytosis) have recently been shown to be even more diverse than their phototrophic counterparts (Singer et al., 2021). As the most important microbial predators in aquatic and soil ecosystems, phagotrophic protists transfer carbon and energy between micro- and macrofauna, significantly influencing ecosystem-level processes (Sherr and Sherr, 1984, 2002; Risse-Buhl et al., 2015; Geisen et al., 2017).
Due to their ubiquity, abundance and sensitivity, phagotrophic protists have been recognized as excellent bioindicators in a variety of environments (Payne, 2013). Ciliates (Ciliophora) and testate amoebae (Amoebozoa and Rhizaria) are very common groups of phagotrophic protists in freshwater and soil ecosystems (Singer et al., 2021). They are also the best-studied protists, partly because they are easy to isolate and have a distinctive morphology (Foissner, 1999; Todorov and Bankov, 2019). Naked amoebae, on the other hand, are understudied due to their smaller size and a lack of knowledge on their morphological features, although they dominate many protist communities, particularly in periphyton and benthos, as well as in soil (Coleman and Wall, 2015; Geisen, 2016; Singer et al., 2021).
Different protist groups may play opposing or complementary roles in ecosystem functioning. Thus, analyzing their functional diversity (FD) may enhance our understanding of various protist traits and their role in the ecosystem (Dumack et al., 2020). Functional diversity analysis is based on linking ecosystem processes to species’ functional traits, i.e., organisms’ characteristics that influence not only their fitness, but also their interactions with the environment (Nock et al., 2016). Trait-based approaches offer an opportunity to understand the relationship between microbial diversity and ecosystem functioning better than taxonomy-based approaches, especially since taxonomic and functional diversity are not always linearly correlated (Laliberté and Legendre, 2010). However, the trait-based approach has rarely been applied to protists (Weisse, 2017; Fiore-Donno et al., 2019).
Phagotrophic protists have already been recognized as good bioindicators of the karst environments’ (e.g., tufa-depositing lakes, rivers, and streams) quality (Primc-Habdija et al., 2001; Matoničkin Kepčija et al., 2011; Kulaš et al., 2021). A recent study by Gulin et al. (2021) showed that periphyton responds to revitalization of dry tufa-depositing streams, opening new possibilities for the application of periphytic organisms as proxies of the stream revitalization success. The study focused on evaluating the response of periphyton based on sampling natural substrates, i.e., stream substrate containing tufa particles, which allowed analysis not only of the abundance and diversity of periphytic taxa, but also of the granulometric and mineral constituents of the substrate itself. However, the method of sampling with a corer proved too aggressive for periphytic organisms, resulting in low abundance and diversity of periphyton. In contrast, the present study is an improvement, because it is based on evaluating the response by monitoring periphyton development on artificial substrates over time so that the periphyton can be studied directly without damaging it. It focuses on phagotrophic protists (ciliates and amoeboid protists) and analyses both their taxonomic and functional diversity. Since periphyton exhibits complex spatio-temporal colonization trends (Franco et al., 1998; Risse-Buhl and Küsel, 2009), 1 and 2-months immersion periods were covered.
Our first objective was to assess whether taxonomic (abundance, species richness and diversity) and functional (i.e., functional diversity) metrics of ciliates and amoeboid protists differed between control (C) sites (located within permanent streams) and revitalized (N) sites (located within reactivated streams). We expected higher abundance and diversity (taxonomic and functional) of ciliates and amoeboid protists at revitalized sites than at control sites. This was based on the fact that conditions that prevail at reactivated streams due to newly created waterways and soil drainage, represent a positive disturbance, which can improve the diversity of habitats, species and their survival strategies (Timoner et al., 2014). A better understanding of the effects of disturbance on functional diversity is critical because functional diversity is likely related to ecosystem resilience, i.e., the capacity of a system to absorb shocks, reorganize, and maintain the same structure and function (Walker, 1995).
The second objective was to determine the influence of several environmental factors (flow velocity, temperature, dissolved oxygen concentration, pH, conductivity, alkalinity, total water hardness, concentrations of nitrites, nitrates and orthophosphates, and chemical oxygen demand) on the taxonomic and functional metrics of ciliates and amoeboid protists. Some of these factors have already shown a significant effect on the assemblage metrics of periphytic taxa in revitalized tufa streams (Gulin et al., 2021). Based on previous findings, we expected several organic matter associated factors (e.g., chemical oxygen demand, orthophosphates) to increase the abundance and taxonomic diversity of phagotrophic protists while others (e.g., nitrites) were expected to yield the opposite (decreasing) effect (Gulin et al., 2021).
Our third objective was to estimate how periphyton-associated factors (organic and inorganic matter (tufa) content, chlorophyll a concentration) affect the taxonomic and functional metrics of ciliates and amoeboid protists. Based on the findings of Gulin et al. (2021), we expected that increased organic matter content (and subsequently chlorophyll a concentration) would lead to increased abundance of protists. On the other hand, we expected that increased inorganic matter (i.e., tufa) content would decrease the taxonomic and functional diversity of protists as intense tufa deposition can lead to detachment and sloughing (Pitois et al., 2001, 2003).
Finally, our fourth objective was to examine the extent to which duration of immersion in water influences changes in taxonomic and functional metrics of ciliates and amoeboid protists at revitalized sites. We expected to observe a trend of increasing values of taxonomic and functional metrics with longer immersion duration, for both ciliates and amoeboid protists from revitalized sites.
The links between ecosystem functioning and the taxonomic and functional diversity of phagotrophic protists could be fundamental to resource management and conservation planning, particularly in protected and sensitive areas such as national parks. The results of this research would make an important contribution to understanding fragile and complex tufa-depositing ecosystems.
Materials and Methods
Study Area
The study was conducted on a 1 ha experimental plot at the Skradinski buk tufa barrier-the longest and final tufa barrier in the watercourse of the Krka River, a karst river in the Dinaric region of Croatia, protected as a national park since 1985. Following an extensive aerial survey (Phantom 4) and detailed vegetation mapping, which determined that the invasive plant Ailanthus altissima (Mill.) Swinge (tree of heaven) dominated the area and was causing the dryness of the barrier due to its strong root system, individual specimens of the invasive tree were mechanically removed in August 2017 with permission from the Croatian Agency for Environment and Nature. Within 2 months of the removal, five streams in the experimental plot, that had previously dried up completely, were reactivated. The sampling design included seven sampling sites selected after removal of an invasive plant species: two control sites (C) in permanent streams where water was present before and after plant removal, C1 (15.966381, 43.805752)-control site where water had been present before and after the removal of A. altissima and displaying well-developed moss cover; C2 (15.966279, 43.805772)-control site without moss cover and five revitalized (N) sites: N1 (15.965449, 43.806438), N2 (15.965324, 43.806624), N3 (15.965246, 43.806541), N4 (15.965186, 43.806489), N5 (15.965538, 43.806225), representing newly reactivated streams (Figure 1).
Figure 1. (A) Experimental plot at Skradinski buk barrier with graphical presentation before and after the removal of A. altissima. Gray lines in A1 represent dry streams which were revitalized upon invasive plant removal. A2 shows seven sampling sites: two control sites located within permanent streams, C1 and C2, and five revitalized sites N1–N5. (B) Details from the experimental plot: 1—dry streams overgrown with A. altissima, 2—revitalized streams upon the removal of A. altissima, 3—3D printed Plexiglass carriers, 4—Plexiglass carriers with glass slides placed in stream bed, 5—a detail of a glass slide from N site with deposited inorganic matter, i.e., tufa, 6—well-developed moss cover from C1 site.
Sampling Design and Identification
On May 10, 2019, 14 3D-printed Plexiglass carriers were placed, two per each sampling site, each containing three glass slides (7.6 × 2.6 cm) to serve as artificial substrates for periphyton development (Figure 1). Prior to immersion into the stream, the slides were washed in detergent, 1 M hydrochloric acid and distilled water. The total effective surface area of each glass slide was 17.18 cm2 because it was partially covered by Plexiglass. Sampling was designed to immerse slides in the water for a period of 1 or 2 months, with three slides collected per sampling event at each site. Due to the high seasonality of the Krka River (Schöll et al., 2012), four sampling events were chosen for each immersion period, covering four different seasons. The dates on which slides were placed in the water and collected for the 1-month period were May 10–June 10, 2019 (spring), June 10–July 15, 2019 (summer), September 30–November 3, 2019 (autumn), February 2–March 8, 2020 (winter). The dates on which the slides were immersed for the 2-months period were placed and collected were June 10–August 15, 2019 (summer), September 30–December 2, 2019 (autumn), December 2, 2019–February 2, 2020 (winter), March 8–June 5, 2020 (spring, immersion period prolonged due to COVID-19 pandemic). A total of 168 slides were placed, 84 for the 1-month immersion period and another 84 for the 2-months immersion period. However, due to unpredictable events such as disturbance by visitors (sampling sites were located in close proximity to visitor trails and probably aroused curiosity) and seasonal drying of streams, only 55 slides were collected for the 1-month immersion period and 52 for the 2-months immersion period. More slides were collected from N sites due to the higher number of sampling sites in the first place. After collection, the slides were transferred to plastic containers filled with a small amount of ambient water and stored at 4°C in the dark. They were then examined directly under the microscope using Zeiss Axioimager A2 with DIC objectives and Axiocam 305 digital camera, within 24 h of sampling, except for few slides that were examined within 48 h. Ciliates and amoeboid protists were identified at species level using Zen 2.4 imaging software and relevant literature (Kahl, 1930; Page, 1976, 1991; Foissner et al., 1991, 1992, 1994, 1995; Foissner and Berger, 1996; Smirnov and Brown, 2004; Smirnov et al., 2011; Todorov and Bankov, 2019). For each taxon, 10–15 photomicrographs were taken and subjected to morphometric measurement (Figure 2). Additional video clips were used to record movements and distinguishing features.
Figure 2. Photomicrographs of different organisms observed during the study period. The organisms were examined directly under the microscope using Zeiss Axioimager A2 with DIC objectives and captured using Axiocam 305 digital camera. (A) Ciliate Acineta sp., (B) ciliate Lacrymaria sp., (C) ciliate Stentor sp., (D) various naked amoebas.
Environmental Factors
The following environmental factors were measured at each sampling site using the respective portable field meters: temperature (T) and dissolved oxygen concentration (DO) (oximeter OXI 96, WTW GmbH, Weilheim, Germany), pH (pH meter 330i, WTW GmbH, Weilheim, Germany), conductivity (Cond) (conductometer Sension 5, Hach, Loveland, Colorado, United States), and flow velocity (FV) (flow velocity meter P600, Dostmann electronic GmbH, Wertheim- Reicholzheim, Germany). An additional 1 L water sample was collected and stored at 4°C for subsequent laboratory analysis of the water. This included the analysis of alkalinity (Alk), total water hardness (TWH), concentrations of nitrite (N-NO2–), nitrate (N-NO3–) and orthophosphate (P-PO43–) (according to APHA, 1985), and total chemical oxygen demand (COD) using the standardized acidic potassium permanganate titrimetric method (Deutsches Institut für Normung, 1986).
Periphyton-Associated Factors (Organic/Inorganic Matter Content and Chlorophyll a Concentration)
After microscopic examination, periphyton from each slide was divided into two equal parts-one was used to determine the organic/inorganic content and the other was scraped to measure the chlorophyll a concentration (a proxy of primary productivity) (Falkowski and Raven, 1997). For the determination of organic matter content as ash-free dry mass, the samples were dried at 104°C to constant weight, then ashed at 400°C for 4 h and reweighed. The mass difference between the dried sample and the ashed sample was used to express the amount of organic matter while the mass difference between the ashed sample and the glass slide was used to express the amount of inorganic matter, i.e., deposited tufa. Values were expressed as the mass of organic/inorganic matter (mg) content per cm2 surface area. Chlorophyll a concentration was determined by the ethanol extraction method (Nusch, 1980). Values were expressed as mass of chlorophyll a (μg) per 8.59 cm2 surface area.
Statistical Analyses
Data Exploration and Visualization
The data were summarized and displayed using standard statistical measures (mean and standard deviation) and graphically presented to illustrate possible trends using ggplot2 package v. 3.3.5 (Wickham, 2016) in R v. 4.1.1 (R Core Team, 2021).
Species richness, Shannon and Simpson diversity indices were calculated independently for ciliates and amoeboid protists as a measure of alpha diversity using R vegan package v. 2.5.7 (Oksanen et al., 2020). Diversity indices were subsequently converted to the effective number of species (true diversity) following a procedure proposed by Jost (2006, 2007). Abundance was calculated independently for ciliates and amoeboid protists as the number of individuals per cm2 surface area for 1 and 2-months immersion periods.
Functional traits were assigned to each species examined in this study based on the relevant literature for ciliates and amoeboid protists (Foissner and Berger, 1996; Foissner et al., 2002; Adl et al., 2019; Fiore-Donno et al., 2019). Ciliates were classified into categories (functional groups) based on the following functional traits: food source, feeding strategy, ecosystem preference, habitat preference, motility, mode of locomotion and life form. Likewise, the amoeboid protists were assigned to the following functional categories (groups): food source, habitat, and morphology (the presence or absence of shell). A detailed overview of the criteria for classification can be found in Supplementary Tables 1, 2.
To minimize the effects of qualitative and quantitative selection of functional traits on the results of functional diversity measures that are shown to be potentially significant (Pakeman, 2014; Zihao et al., 2021), we focused on the robust measures of functional dispersion (FDis) and RaoQ quadratic diversity (RaoQ) to quantify the functional diversity of ciliates and amoeboid protists. FDis is defined as the weighted mean distance in multidimensional trait space of individual species to the centroid of all species, where the weights are the relative abundances of the species (Laliberté and Legendre, 2010; Cappelatti et al., 2020). RaoQ is defined as the sum of pairwise distances between species in multidimensional trait space weighted by their relative abundance (Ricotta and Moretti, 2011). By construction, the two functional dispersion indices are not influenced by species richness (Teittinen and Virta, 2021). FDis and RaoQ values were calculated in R using the FD package v. 1.0–12 (Laliberté and Legendre, 2010; Laliberté et al., 2014). Indices of functional diversity were calculated separately for ciliate and amoeboid protist assemblages at C and N sites for the 1 and 2-months immersion periods.
Data Analysis
Each of the four data sets (1 and 2-months immersion periods for ciliates and amoeboid protists), environmental data and data regarding periphyton-associated factors were tested for normality using Shapiro-Wilk’s test in R prior to further analysis. Since the data were not normally distributed (Shapiro-Wilk’s test, p < 0.05) and sphericity was violated (Mauchly’s test, p < 0.05), differences in environmental and periphyton-associated factors between C and N sites were tested using the analysis of similarity (ANOSIM), a non-parametric test for evaluating a dissimilarity matrix instead of raw data (Clarke and Warwick, 1994) from the R package vegan v. 2.5.7 (Oksanen et al., 2020).
Generalized linear mixed models (GLMMs) were constructed using SPSS Statistics ver. 28.0 (IBM Corp, 2021) to test for differences in ciliate and amoeboid protist assemblages between C and N sites for the 1 and 2-months immersion periods. The taxonomic metrics considered were: abundance, species richness, taxonomic (True) diversity derived from the Shannon and Simpson indices while the functional diversity indices included FDis and RaoQ. Only those environmental factors that showed significant differences (ANOSIM) between C and N sites were selected as inputs to the models. GLMMs with the variables “site,” “COD,” “TWH” and “nitrites” as fixed effects were constructed to evaluate the relationships between 1-month assemblages of ciliates and amoeboid protists and environmental factors. For the 2-months immersion period, the variables “site,” “DO,” “conductivity,” “pH,” “COD,” “alkalinity,” “TWH” and “orthophosphates” were selected as fixed effects. The variables “replicate” and “season” were included as random effects for both model variations, by recommendation of Jost (2007).
Additional GLMMs with periphyton-associated factors (organic and inorganic matter content, chlorophyll a concentration) as fixed effects were constructed for both ciliate and amoeboid protist assemblages for the 1 and 2-months immersion periods. Due to repeated measures sampling setup, first-order autoregressive (AR1) covariance structure of random effects was assumed in all models constructed within this study (Field, 2009).
Results
Environmental Factors
Values of all environmental factors measured in this research can be found in Supplementary Table 3. The COD values for the 1-month period at N sites were significantly lower (R = 0.45, p = 0.015) in comparison to C sites while the opposite was found for nitrites (R = 0.64, p = 0.007) (Figure 3).
Figure 3. Box plots showing selected environmental factors at control (C) and revitalized sites (N) for the 1 and 2-months immersion periods. Asterisk symbol (*) indicates statistically significant differences among estimated means (ANOSIM, p < 0.05). Upper and lower edges of the boxes are the first and third quartiles; the line inside the box represents the median; individual dots are outliers. Abbreviations: COD—total chemical oxygen demand, N-NO2–—nitrites, P-PO43–—orthophosphates.
For the 2-months immersion period, N sites had significantly lower values of COD (R = 0.52, p = 0.004), pH (R = 0.69, p = 0.015) and orthophosphates (R = 0.70, p = 0.001), while higher values of conductivity (R = 0.91, p = 0.004) in comparison to C sites (Figure 3).
Periphyton-Associated Factors
All periphyton-associated factors differed significantly between C and N sites for both immersion periods (R = 1.00, p < 0.001). Organic matter content was lower at N sites for both periods (Figure 4), but it could be noticed that N sites occasionally had extremely high values (30.83 mg per cm2 during the 1-month period and 7.22 mg per cm2 during the 2-months period), while C sites had no such extreme outliers. Inorganic matter content was lower at N sites than at the C sites during the 1-month period, but showed the opposite during the 2-months period, reaching record values of 143.60 mg per cm2. Chlorophyll a concentration at N sites was higher compared to C sites during the 1-month period and lower during the 2-months period. However, it could be noticed that N sites, unlike C sites, had some extremely high chlorophyll a concentrations during both periods.
Figure 4. Box plots showing periphyton-associated factors at control (C) and revitalized sites (N) for the 1 and 2-months immersion periods. Asterisk symbol (*) indicates statistically significant differences among estimated means (ANOSIM, p < 0.05). Upper and lower edges of the boxes are the first and third quartiles; the line inside the box represents the median; individual dots are outliers.
Taxonomic and Functional Metrics of Ciliate Assemblages
In this study, a total of 78 ciliate species were identified, most of them (48) at N sites. Ciliate taxonomic metrics (abundance, species richness, Shannon and Simpson derived True diversity indices) were higher at N sites than at C sites for both immersion periods (Figure 5). Although no significant differences were confirmed by GLMMs with environmental and periphyton-associated factors for the 1-month period, all taxonomic metrics showed significant differences between N sites and C sites for the 2-months period (Tables 1, 2). As for correlations between taxonomic metrics and environmental factors, no significant correlations were found for the 1-month period, but several were observed for the 2-months period: abundance was significantly positively correlated with COD, species richness was significantly negatively correlated with conductivity, while True diversity values (Shannon and Simpson) were significantly positively correlated with pH and orthophosphates but significantly negatively correlated with conductivity (Table 3). Periphyton-associated factors were found to have significant correlations with taxonomic metrics for both 1 and 2-months periods (Table 4). All taxonomic metrics except abundance were positively correlated with inorganic matter and negatively correlated with organic matter and chlorophyll a for the 1-month period. The opposite was found for the 2-months period when all taxonomic metrics were positively correlated with organic matter and chlorophyll a, while negatively correlated with inorganic matter. Significant correlations were found for species richness and True diversity (Shannon) values for both periods and abundance for the 2-months period.
Figure 5. Box plots showing ciliate taxonomic and functional metrics at control (C) and revitalized sites (N) for the 1 and 2-months immersion periods. Asterisk symbol (*) indicates statistically significant differences among estimated means (GLMM, p < 0.05). Upper and lower edges of the boxes are the first and third quartiles; the line inside the box represents the median; individual dots are outliers.
Table 1. GLMM (full model) output showing differences in ciliate assemblage metrics between control and revitalized sites (1 and 2-months immersion periods) with sites and environmental factors as fixed effects, and season and replicate as random effects.
Table 2. GLMM (full model) output showing differences in ciliate assemblage metrics between control and revitalized sites (1 and 2-months immersion periods) with sites and periphyton-associated factors as fixed effects, and season and replicate as random effects.
Table 3. GLMM (full model) output showing main effects of environmental factors (fixed effects) on taxonomic and functional metrics of ciliate assemblages, with season and replicate as random effects.
Table 4. GLMM (full model) output showing main effects of periphyton-associated factors (fixed effects) on taxonomic and functional metrics of ciliate assemblages, with season and replicate as random effects.
FDis values for ciliates were lower at N sites compared to C sites for both immersion periods, whereas RaoQ values were higher at N sites for the 1-month period but lower for the 2-months period (Figure 5). When tested with GLMMs including environmental and periphyton-associated factors, significant differences were confirmed only for FDis values for the 1-month period and RaoQ values for the 2-months period (Tables 1, 2). Both functional metrics for the 2-months period were negatively correlated with conductivity, while they were positively correlated with pH, COD, and orthophosphates, although the correlations were significant only for RaoQ (Table 3). As for periphyton-associated factors, FDis and RaoQ showed quite different correlations for the 1-month period: FDis was significantly positively correlated with organic matter and negatively correlated with inorganic matter, whereas RaoQ was significantly negatively correlated with organic matter and positively correlated with inorganic matter. Both functional metrics were then negatively correlated with chlorophyll a, while a significant correlation was found only for RaoQ. With regards to the 2-months period, both FDis and RaoQ were positively correlated with chlorophyll a and negatively correlated with organic and inorganic matter. Correlations with inorganic matter were the only found significant (Table 4).
The detailed GLMM report with all correlations between environmental factors and the 1 and 2-months taxonomic and functional metrics of ciliates can be found in Supplementary Table 4.
In terms of immersion duration, it could be noticed that all ciliate taxonomic metrics at N sites increased with longer immersion. The same could be observed for FDis while RaoQ values remained the same.
Taxonomic and Functional Metrics of Amoeboid Protists
A total of 52 species of amoeboid protists were identified in this study, most of which (41) belonged to N sites. Almost no differences were observed between N and C sites with respect to most taxonomic metrics for both immersion periods (Figure 6). Abundance was slightly lower at N sites for the 1-month period while True diversity Shannon values were significantly lower at N sites for the 2-months period, when tested with GLMMs that included environmental and periphyton-associated factors (Tables 5, 6). As for correlations with environmental factors, all taxonomic metrics for the 1-month period were positively correlated with COD, while the correlation with abundance was found significant. With regards to the 2-months period, all taxonomic metrics negatively correlated with COD and orthophosphates, while a positive correlation was found with pH and conductivity. Most of these correlations were significant (Table 7). As for correlations with periphyton-associated factors, all taxonomic metrics for the 1-month period were positively correlated with all three factors, while significant correlations were found between species richness and inorganic matter. For the 2-months period, there was a significant positive correlation between abundance and organic matter, while all taxonomic metrics were significantly negatively correlated with chlorophyll a (Table 8).
Figure 6. Box plots showing taxonomic and functional metrics of amoeboid protists at control (C) and revitalized sites (N) for the 1 and 2-months immersion periods. Asterix symbol (*) indicates statistically significant differences among estimated means (GLMM, p < 0.05). Upper and lower edges of the boxes are the first and third quartiles; the line inside the box represents the median; individual dots are outliers.
Table 5. GLMM (full model) output showing differences in amoeboid protist assemblage metrics between control and revitalized sites (1 and 2-months immersion periods) with sites and environmental factors as fixed effects, and season and replicate as random effects.
Table 6. GLMM (full model) output showing differences in amoeboid protist assemblage metrics between control and revitalized sites (1 and 2-months immersion periods) with sites and periphyton-associated factors as fixed effects, and season and replicate as random effects.
Table 7. GLMM (full model) output showing main effects of environmental factors (fixed effects) on taxonomic and functional metrics of amoeboid protist assemblages, with season and replicate as random effects.
Table 8. GLMM (full model) output showing main effects of periphyton-associated factors (fixed effects) on taxonomic and functional metrics of amoeboid protist assemblages, with season and replicate as random effects.
FDis and RaoQ values were higher at N sites compared with C sites for the 1-month period, whereas they were lower for the 2-months period (Figure 6). The only significant difference was found for the 1-month period. Regarding correlations with environmental factors, both functional indices were negatively correlated with COD for both immersion periods, while significant correlations were found only for the 2-months period. In addition, FDis was significantly positively correlated with pH and orthophosphates (Table 7). As for periphyton-associated factors, both functional indices were positively correlated with all periphyton-associated factors for the 1-month period, whereas they were negatively correlated with chlorophyll a for the 2-months period, although not significantly (Table 8).
The detailed GLMM report with all correlations between environmental factors and the 1 and 2-months taxonomic and functional metrics of amoeboid protists can be found in Supplementary Table 5.
With respect to immersion duration, species richness and True diversity (Shannon and Simpson) values at N sites were found to increase with longer immersion, while abundance decreased slightly. The functional indices remained more or less the same with longer immersion duration.
Discussion
Our results revealed that ciliates and amoeboid protists respond differently to stream revitalization at the assemblage level. While ciliate assemblages at N sites had higher abundance and taxonomic diversity compared to those at C sites for both immersion periods, taxonomic metrics of amoeboid assemblages differed little between sites, even with longer immersion. The two protist assemblages also differed in their functional response: while the functional diversity of ciliates was lower at N sites compared to C sites for both immersion periods, that of amoeboid protists was higher at N sites for the 1-month period but lower for the 2-months period, in comparison to C sites.
The lower functional diversity of ciliates at N sites compared to C sites for both immersion periods suggests that the level of disturbance in reactivated streams, in the form of occasionally changing hydrologic conditions (high flow or drought), was either too high or too low to promote functional diversity, judging by the intermediate disturbance hypothesis (IDH). IDH rests on the assumption that a maximum level of diversity is achieved when the assemblage is exposed to an intermediate level of disturbance (Connell, 1978). It is likely that in our case, disturbance filtered out incompatible suites of traits, so that only a subset of species with disturbance-resistant traits could colonize a disturbed habitat (Lavorel et al., 2008; Biswas and Mallik, 2010). Thus, compared to C sites, ciliate assemblages at N sites could be expected to consist mainly of prospective colonizers, with disturbance-intolerant species (r-strategists) absent or their share in the assemblage minimized, leading to increased similarity and functional redundancy within the assemblage (Segovia et al., 2016). On the other hand, it appears that amoeboid assemblages at N sites during the 1-month period consisted mainly of highly tolerant or even drought-resistant species (K-strategists) that would thrive under hydrologically extreme conditions (Carballeira and Pontevedra-Pombal, 2021), in contrast to the r-strategists that probably dominated the 2-months period, leading to differences in functional diversity between the two periods.
Our results suggested that the reasons for these different responses to stream revitalization in terms of taxonomic and functional metrics could also lie in the effects of environmental and periphyton-associated factors. While some factors had a similar effect on both assemblages, others affected the assemblages quite differently. The overall observed positive effect of COD, pH, and organic matter content on taxonomic metrics of both ciliates and amoeboid protists is consistent with previous findings by numerous authors who have worked with these protists (Cowling, 1994; Smith, 1996; Verhoeven, 2002; Bates et al., 2013; Bradford, 2016; Geisen, 2016; Fiore-Donno et al., 2019). The likely reason for this is the increased diversity and quantity of bacteria that are their main food source, as the deposited organic material supports attachment of bacteria that decompose both dissolved and particulate organic matter (Foissner, 2014). However, the dominance of bacteria as a food source could have led to functional homogeneity and explain why functional diversity of both ciliates and amoeboid protists was negatively affected by COD and/or organic matter.
Inorganic material, i.e., tufa, on the other hand, had a different effect. While it positively affected the taxonomic and functional metrics of both protist groups for the 1-month period, the same was not observed for the 2-months period. It appears that the higher inorganic matter content during the 2-months period, whose deposition was facilitated by slightly higher temperature, conductivity, and flow velocity at N sites (Supplementary Table 3), negatively affected the ciliate taxonomic and functional metrics, but did not cause any negative effects on those of the amoeboid assemblages. Although the presence of inorganic matter can increase microhabitat heterogeneity (Beisel et al., 2000; Bednar et al., 2017; Singer et al., 2021) and thus serve as an “inoculum” for periphyton supporting high abundance and diversity (Primc and Habdija, 1987; Wörner et al., 2000; Zimmermann-Timm, 2002), there appears to be a critical point at which inorganic matter content no longer promotes high abundance and taxonomic diversity of ciliates. It could be that ciliates are much more sensitive to burial and sloughing than amoeboid protists because they do not predominantly glide along the substrate as amoebae do but exhibit various modes of locomotion: many of them are free-swimmers or even move by jumping and rotating (Foissner and Berger, 1996; Esteban and Fenchel, 2021). Amoeboid protists (especially naked amoebae) are also smaller in size compared to ciliates, which allows them to savor the microhabitats and be very closely associated with the substrate (Preston, 2003).
The two assemblages also responded differently with respect to immersion duration: while ciliate taxonomic and functional metrics at N sites increased with longer immersion, the same was observed only for some taxonomic metrics of amoeboid protists (species richness and True diversity values), while the abundance decreased, and functional metrics remained more or less the same. Since chloropyhll a concentration showed significant correlations with most taxonomic and functional metrics of both assemblages and reached extreme values at times during the 2-months period at N sites, it is reasonable to assume that the different responses to immersion duration could be due to this. Although the observed orthophosphate levels during the 2-months period were still quite low and within the range characteristic of tufa-depositing streams (Primc-Habdija et al., 2005), they may have caused algal proliferation and subsequently higher chlorophyll a concentration. The dense algal coating of the substrate may have contributed to heterogeneity of food sources and microhabitats (Algarte et al., 2017) which together led to niche diversification, an increase in the share of disturbance-resistant ciliate species and subsequently an increase in taxonomic and functional diversity. Another possible explanation could be that changing hydrological conditions at N sites activate the dormant part of the microbial assemblages, leading to increased functional diversity over time (Velasco-González et al., 2020). On the other hand, the alginate covering the bacteria likely triggered a different response in the amoeboid assemblages, leading to reduced or even absent uptake of bacteria by amoeboid protists. This is likely due to the alginate coating increasing the overall size of the bacterial particles or altering their “flavor” (Heaton et al., 2001; Parry, 2004), resulting in reduced abundance of amoeboid protists during the 2-months immersion period.
The taxonomic and functional diversity of protists is a valuable indication of environmental conditions and ecosystem stability. Understanding the factors that structure protist assemblages is a prerequisite for using these organisms to predict environmental quality and future changes in ecosystem functioning. Protists have been used in environmental monitoring, but only to a limited extent because of the difficulties associated with sophisticated morphological identification methods. However, comprehensive studies that also consider their functional traits (Fournier et al., 2012) could make them more recognizable and increase the frequency of using protists as indicators of ecosystem functioning.
Conclusion
Linking taxonomic and functional metrics to environmental conditions can improve our understanding of biological processes, especially in ecosystems undergoing extensive changes, such as revitalized streams. Our results showed that taxonomic and functional metrics of ciliates and amoeboid protists responded to the prevailing conditions characteristic of revitalized tufa-depositing streams: changing hydrology (occasional high flow or drought), soil drainage, and extensive inorganic matter, i.e., tufa deposition, although their responses were somewhat different. While both protist assemblages benefited from COD and the presence of organic matter supported by favorable pH, as this provided their main food source–bacteria, ciliates were more sensitive to the presence of inorganic matter content than amoeboid protists. However, in terms of chlorophyll a concentration, ciliates appeared to benefit more from the algal coating as it resulted in heterogeneity of food sources and microhabitats, which increased their taxonomic and functional diversity, while it decreased the uptake of bacteria for amoeboid protists, negatively affecting their abundance. The two assemblages also showed different responses of taxonomic and functional metrics with respect to immersion duration: while the taxonomic and functional diversity of ciliates at N sites increased with longer immersion, indicating niche diversification, those of amoeboid protists hardly changed with time. The results presented suggest that a comprehensive analysis of taxonomic and functional metrics of ciliates and amoeboid protists could be a good proxy for assessing revitalization of tufa-depositing streams. However, the temporal component should always be considered when conducting such studies, as the colonization processes of ciliates and amoeboid protists are quite complex, especially in tufa-depositing streams.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
VG, MSP, and RMK contributed to conception and design of the study. VG and BV organized the trait database. VG and FR performed the statistical analysis. VG wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This research was conducted as a part of the project “Monitoring the recovery of periphyton and macrozoobenthos at tufa barriers in the Krka National Park after tree of heaven removal” 106-F19-00015, financed by Krka National Park, Croatia.
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.
Acknowledgments
We would like to thank Gordana Goreta who supervised the project and Nikolina Smolić who helped with the selection of sampling sites in the Krka National Park; Antun Beljak for his assistance in the field, Sunčica Bosak for loaning microscope equipment and Edin Cvolić for his help with illustrations.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.842395/full#supplementary-material
References
Adl, S. M., Bass, D., Lane, C. E., Lukeš, J., Schoch, C. L., Smirnov, A., et al. (2019). Revisions to the classification, nomenclature, and diversity of eukaryotes. J. Eukaryot. Microbiol. 66, 114–119. doi: 10.1111/jeu.12691
Algarte, V. M., Siqueira, T., Landeiro, V. L., Rodrigues, L., Bonecker, C. C., Rodrigues, L. C., et al. (2017). Main predictors of periphyton species richness depend on adherence strategy and cell size. PLoS One 12:e0181720. doi: 10.1371/journal.pone.0181720
APHA (1985). Standard Methods for the Examination of Water and Wastewater. Washington, DC: American Public Health Association.
Bates, S. T., Clemente, J. C., Flores, G. E., Walters, W. A., Parfrey, L. W., Knight, R., et al. (2013). Global biogeography of highly diverse protistan communities in soil. ISME J. 7, 652–659. doi: 10.1038/ismej.2012.147
Bednar, J. P., Trobej, M., Schagerl, M., and Waringer, J. (2017). Which factors shape macrozoobenthic communities in tufa springs? Results from Austrian meteogene travertine-depositing sites. Hydrobiologia 799, 293–307. doi: 10.1007/s10750-017-3228-x
Beisel, J. N., Usseglio-Polatera, P., and Moreteau, J. C. (2000). The spatial heterogeneity of a river bottom: a key factor determining macroinvertebrate communities. Hydrobiologia 422, 163–171. doi: 10.1023/A:1017094606335
Biswas, S. R., and Mallik, A. U. (2010). Disturbance effects on species diversity and functional diversity in riparian and upland plant communities. Ecology 91, 28–35. doi: 10.1890/08-0887.1
Bradford, M. A. (2016). Re-visioning soil food webs. Soil Biol. Biochem. 102, 1–3. doi: 10.1016/j.soilbio.2016.08.010
Cappelatti, L., Mauffrey, A. R. L., and Griffin, J. N. (2020). Functional diversity of habitat formers declines scale-dependently across an environmental stress gradient. Oecologia 194, 135–149. doi: 10.1007/s00442-020-04746-1
Carballeira, R., and Pontevedra-Pombal, X. (2021). Diversity of Testate Amoebae as an Indicator of the Conservation Status of Peatlands in Southwest Europe. Diversity 13:269. doi: 10.3390/d13060269
Clarke, K. R., and Warwick, R. M. (1994). Similarity-based testing for community pattern: the two-way layout with no replication. Mar. Biol. 118, 167–176. doi: 10.1007/BF00699231
Coleman, D. C., and Wall, D. H. (2015). “Soil fauna: occurrence, biodiversity, and roles in ecosystem function,” in Soil Microbiology, Ecology and Biochemistry, 4th Edn, ed. E. A. Paul (Boston: Academic Press), 111–149. doi: 10.1016/B978-0-12-415955-6.00005-0
Connell, J. H. (1978). Diversity in tropical rain forests and coral reefs: high diversity of trees and corals is maintained only in a nonequilibrium state. Science 199, 1302–1310. doi: 10.1126/science.199.4335.1302
Cowling, A. J. (1994). “Protozoan distribution and adaptation,” in Soil Protozoa, ed. J. F. Darbyshire (Wallingford: CAB International), 5–42.
de Vargas, C., Audic, S., Henry, N., Decelle, J., Mahé, F., Logares, R., et al. (2015). Eukaryotic plankton diversity in the sunlit ocean. Science 348:1261605. doi: 10.1126/science.1261605
Deutsches Institut für Normung (1986). Deutsche Einheitsverfahren zur Wasser-, Abwasser- und Schlammuntersuchung, 16th Edn. Weinheim: Verlag Chemie.
Dumack, K., Fiore-Donno, A. M., Bass, D., and Bonkowski, M. (2020). Making sense of environmental sequencing data: ecologically important functional traits of the protistan groups Cercozoa and Endomyxa (Rhizaria). Mol. Ecol. Resour. 20, 398–403. doi: 10.1111/1755-0998.13112
Esteban, G. F., and Fenchel, T. M. (2021). Ecology of Protozoa: The Biology of Free-living Phagotrophic Protists, 2nd Edn. Switzerland AG: Springer Nature.
Fiore-Donno, A. M., Richter-Heitmann, T., Degrune, F., Dumack, K., Regan, K. M., Marhan, S., et al. (2019). Functional traits and spatio-temporal structure of a major group of soil protists (Rhizaria: Cercozoa) in a temperate grassland. Front. Microbiol. 10:1332. doi: 10.3389/fmicb.2019.01332
Foissner, W. (1999). Protist diversity: estimates of the near-imponderable. Protist 150, 363–368. doi: 10.1016/S1434-4610(99)70037-4
Foissner, W. (2014). ”Protozoa” in Earth Systems and Environmental Sciences. Amsterdam: Elsevier, doi: 10.1016/B978-0-12-409548-9.09130-2
Foissner, W., Agatha, S., and Berger, H. (2002). Soil Ciliates (Protozoa, Ciliophora) from Namibia (Southwest Africa), with Emphasis on Two Contrasting Environments, the Etosha Region and the Namib Desert. Part I: Text and Line Drawings. Linz, Austria: Biologiezentrum des Oberösterreichisches Landesmuseum.
Foissner, W., and Berger, H. (1996). A user-friendly guide to the ciliates (Protozoa, Ciliophora) commonly used by hydrobiologists as bioindicators in rivers, lakes, and waste waters, with notes on their ecology. Freshw. Biol. 35, 375–482. doi: 10.1111/j.1365-2427.1996.tb01775.x
Foissner, W., Berger, H., and Kohmann, F. (1992). Taxonomische und ökologische Revision der Ciliaten des Saprobiensystems – Band II: Peritrichia, Heterotrichida, Odontostomatida. Deggendorf: Bayerisches Landesamt für Wasserwirtschaft.
Foissner, W., Berger, H., and Kohmann, F. (1994). Taxonomische und ökologische Revision der Ciliaten des Saprobiensystems – Band III: Hymenostomata, Prostomatida, Nassulida. Deggendorf: Bayerisches Landesamt für Wasserwirtschaft.
Foissner, W., Berger, H., Blatterer, H., and Kohmann, F. (1995). Taxonomische und ökologische Revision der Ciliaten des Saprobiensystems – Band IV: Gymnostomatea, Loxodes, Suctoria. Deggendorf: Bayerisches Landesamt für Wasserwirtschaft.
Foissner, W., Blatterer, H., Berger, H., and Kohmann, F. (1991). Taxonomische und ökologische Revision der Ciliaten des Saprobiensystems – Band I: Cyrtophorida, Oligotrichida, Hypotrichia, Colpodea. Deggendorf: Bayerisches Landesamt für Wasserwirtschaft.
Fournier, B., Malysheva, E., Mazei, Y., Moretti, M., and Mitchell, E. A. D. (2012). Toward the use of testate amoeba functional traits as indicator of floodplain restoration success. Eur. J. Soil Biol. 49, 85–91. doi: 10.1016/j.ejsobi.2011.05.008
Franco, C., Esteban, G. F., and Tellez, C. (1998). Colonization and succession of ciliated protozoa associated with submerged leaves in a river. Limnologica 28, 275–283.
Geisen, S. (2016). The bacterial-fungal energy channel concept challenged by enormous functional versatility of soil protists. Soil Biol. Biochem. 102, 22–25. doi: 10.1016/j.soilbio.2016.06.013
Geisen, S., Mitchell, E. A. D., Wilkinson, D. M., Adl, S., Bonkowski, M., Brown, M. W., et al. (2017). Soil protistology rebooted: 30 fundamental questions to start with. Soil Biol. Biochem. 111, 94–103. doi: 10.1016/j.soilbio.2017.04.001
Gulin, V., Matoničkin Kepčija, R., Sertić Perić, M., Felja, I., Fajković, H., and Križnjak, K. (2021). Environmental and periphyton response to stream revitalization – a pilot study from a tufa barrier. Ecol. Indic. 126:107629. doi: 10.1016/j.ecolind.2021.107629
Heaton, K., Drinkall, J., Minett, A., Hunt, A. P., and Parry, J. D. (2001). “Ameoboid grazing on surface-associated prey,” in Biofilm Community Interaction– Chance or Necessity?, eds P. Gilbert, D. Allison, M. Brading, J. Verran, and J. Walker (Cardiff: BioLine), 293–301.
Jost, L. (2007). Partitioning diversity into independent alpha and beta components. Ecology 88, 2427–2439. doi: 10.1890/06-1736.1
Kahl, A. (1930). “Urtiere oder Protozoa I: Wimpertiere oder Ciliata (Infusoria),” in Die Tierwelt Deutschlands, ed. F. Dahl (Jena: G. Fisher), 1–886.
Kulaš, A., Gulin, V., Matoničkin Kepčija, R., Žutinić, P., Sertić Perić, M., Orlić, S., et al. (2021). Ciliates (Alveolata, Ciliophora) as bioindicators of environmental pressure: a karstic river case. Ecol. Indic. 124:107430. doi: 10.1016/j.ecolind.2021.107430
Laliberté, E., and Legendre, P. (2010). A distance-based framework for measuring functional diversity from multiple traits. Ecology 91, 299–305. doi: 10.1890/08-2244.1
Laliberté, E., Legendre, P., and Shipley, B. (2014). FD: Measuring Functional Diversity from Multiple Traits, and Other Tools for Functional Ecology R Package Version 1.0-12. Available online at: https://cran.r-project.org/web/packages/FD/ (accessed February 2022).
Lavorel, S., Grigulis, K., McIntyre, S., Williams, N. S. G., Garden, D., Dorrough, J., et al. (2008). Assessing functional diversity in the field – methodology matters! Funct. Ecol. 22, 134–147. doi: 10.1111/j.1365-2435.2007.01339.x
Matoničkin Kepčija, R., Miliša, M., Sertić Perić, M., Matijić Cvjetović, M., and Primc-Habdija, B. (2011). Response of periphyton to nutrient addition in tufa-depositing environment. Aquat. Microb. Ecol. 65:183.195. doi: 10.3354/ame01545
Nock, C. A., Vogt, R. J., and Beisner, B. E. (2016). Functional Traits in: eLS., (Chichester, UK: John Wiley & Sons), 1–8. doi: 10.1002/9780470015902.a0026282
Nusch, E. A. (1980). Comparison of different methods for chlorophyll and phaeopigment determination. Arch. Hydrobiol. Beih. Ergebn. Limnol. 14, 14–36.
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020). Vegan: Community Ecology Package. R package version v. 2.5.7. Available online at: https://CRAN.R-project.org/package=vegan (accessed February 2022).
Page, F. C. (1976). An Illustrated Key to Freshwater and Soil Amoebae: with Notes on Cultivation and Ecology. Ambleside: Freshwater Biological Association.
Page, F. C. (1991). ““Nackte Rhizopoda”,” in Nackte Rhizopoda und Heliozoa, eds F. C. Page and F. J. Siemensma (Stuttgart: Gustav Fischer Verlag), 1–170.
Pakeman, R. J. (2014). Functional trait metrics are sensitive to the completeness of the species’ trait data? Meth. Ecol. Evol. 5, 9–15. doi: 10.1111/2041-210X.12136
Parry, J. D. (2004). “Protozoan grazing of freshwater biofilms,” in Advances in Applied Microbiology, eds A. I. Laskin, J. W. Bennet, and G. M. Gadd (New York: Elsevier), 167–196. doi: 10.1016/S0065-2164(04)54007-8
Payne, R. J. (2013). Seven reasons why protists make useful bioindicators. Acta Protozool. 52, 105–113. doi: 10.4467/16890027AP.13.0011.1108
Pitois, F., Jigorel, A., and Bertru, G. (2001). Colonization dynamics of an encrusting cyanobacterial mat in a hardwater river (Eaulne, France). Geomicrobiol. J. 18, 139–155. doi: 10.1080/01490450151143435
Pitois, F., Jigorel, A., and Bertru, G. (2003). Development of cyanobacterial build-up and evolution of river bed morphology in the chalk stream Eaulne (Upper-Normandy, France). Biodivers. Conserv. 12, 621–636. doi: 10.1023/A:1022412211915
Preston, T. M. (2003). The water-air interface: a microhabitat for amoebae. Eur. J. Protistol. 39, 385–389. doi: 10.1078/0932-4739-00008
Primc, B., and Habdija, I. (1987). Ciliated colonization of artificial substrates in different saprobic conditions in a running water. Acta Hydrochim. Hydrobiol. 15, 487–494. doi: 10.1002/aheh.19870150509
Primc-Habdija, B., Habdija, I., and Plenković-Moraj, A. (2001). Tufa deposition and periphyton overgrowth as factors affecting the ciliate community on travertine barriers in different current velocity conditions. Hydrobiologia 457, 87–96. doi: 10.1023/A:1012265206470
Primc-Habdija, B., Habdija, I., Matoničkin, R., and Špoljar, M. (2005). Development of ciliate community on artificial substrates associated with vertical gradients of environmental conditions in a karstic lake. Arch. Hydrobiol. 164, 513–527. doi: 10.1127/0003-9136/2005/0164-0513
R Core Team (2021). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Ricotta, C., and Moretti, M. (2011). CWM and Rao’s quadratic diversity: a unified framework for functional ecology. Oecologia 167, 181–188. doi: 10.1007/s00442-011-1965-5
Risse-Buhl, U., and Küsel, K. (2009). Colonization dynamics of biofilm associated ciliate morphotypes at different flow velocities. Eur. J. Protistol. 45, 67–74. doi: 10.1016/j.ejop.2008.08.001
Risse-Buhl, U., Schlief, J., and Mutz, M. (2015). Phagotrophic protists are a key component of microbial communities processing leaf litter under contrasting oxic conditions. Freshw. Biol. 60, 2310–2322. doi: 10.1111/fwb.12657
Schöll, F., Birk, S., and Böhmer, J. (2012). XGIG Large River Intercalibration Exercise – WFD Intercalibration Phase 2: Milestone 6 Report – BQE: Phytobenthos. Ispra: European Commission Directorate General JRC Joint Research Centre Institute of Environment and Sustainability, 73.
Segovia, B. T., Lansac-Tôha, F. M., de Meira, B. R., Cabral, A. F., Lansac-Tôha, F. A., and Velho, L. F. M. (2016). Anthropogenic disturbances influencing ciliate functional feeding groups in impacted tropical streams. Environ. Sci. Pollut. Res. 23, 20003–20016. doi: 10.1007/s11356-016-7185-0
Sherr, B. F., and Sherr, E. B. (1984). “Role of heterotrophic protozoa in carbon and energy flow in aquatic ecosystems,” in Current Perspectives in Microbial Ecology, eds M. J. Klug and C. A. Reddy (Washington D.C: American Society for Microbiology), 412–423.
Sherr, E. B., and Sherr, B. F. (2002). Significance of predation by protists in aquatic microbial food webs. Antonie Leeuwenhoek 81, 293–308. doi: 10.1023/A:1020591307260
Sieber, G., Beisser, D., Bock, C., and Boenigk, J. (2020). Protistan and fungal diversity in soils and freshwater lakes are substantially different. Sci. Rep. 10:20025. doi: 10.1038/s41598-020-77045-7
Singer, D., Seppey, C. V. W., Lentendu, G., Dunthorn, M., Bass, D., Belbahri, L., et al. (2021). Protist taxonomic and functional diversity in soil, freshwater and marine ecosystems. Environ. Int. 146:106262. doi: 10.1016/j.envint.2020.106262
Smirnov, A. V., and Brown, S. (2004). Guide to the methods of study and identification of soil gymnamoebae. Protistology 3, 148–190.
Smirnov, A. V., Chao, E., Nassonova, E. S., and Cavalier-Smith, T. (2011). A revised classification of naked lobose amoebae (Amoebozoa: Lobosa). Protist 162, 545–570. doi: 10.1016/j.protis.2011.04.004
Smith, H. G. (1996). Diversity of Antarctic terrestrial protozoa. Biodivers. Conserv. 5, 1379–1394. doi: 10.1007/BF00051984
Teittinen, A., and Virta, L. (2021). Exploring Multiple Aspects of Taxonomic and Functional Diversity in Microphytobenthic Communities: effects of Environmental Gradients and Temporal Changes. Front. Microbiol. 12:668993. doi: 10.3389/fmicb.2021.668993
Timoner, X., Acuña, V., Frampton, L., Pollard, P., Sabater, S., and Bunn, S. E. (2014). Biofilm functional responses to the rehydration of a dry intermittent stream. Hydrobiologia 727, 185–195. doi: 10.1007/s10750-013-1802-4
Todorov, M., and Bankov, N. (2019). An Atlas of Sphagnum-Dwelling Testate Amoebae in Bulgaria. Sofia: Pensoft Publishers.
Velasco-González, I., Sanchez-Jimenez, A., Singer, D., Murciano, A., Díez-Hermano, S., Lara, E., et al. (2020). Rain-fed granite rock basins accumulate a high diversity of dormant microbial eukaryotes. Microb. Ecol. 79, 882–897. doi: 10.1007/s00248-019-01463-y
Verhoeven, R. (2002). Ciliates in coastal dune soils of different stages of development. Eur. J. Soil Biol. 38, 187–191. doi: 10.1016/S1164-5563(02)01147-0
Walker, B. (1995). Conserving biological diversity through ecosystem resilience. Conserv. Biol. 9, 747–752. doi: 10.1046/j.1523-1739.1995.09040747.x
Weisse, T. (2017). Functional diversity of aquatic ciliates. Eur. J. Protistol. 61, 331–358. doi: 10.1016/j.ejop.2017.04.001
Weisse, T., Anderson, R., Arndt, H., Calbet, A., Hansen, P. J., and Montagnes, D. J. S. (2016). Functional ecology of aquatic phagotrophic protists – concepts, limitations, and perspectives. Eur. J. Protistol. 55, 50–74. doi: 10.1016/j.ejop.2016.03.003
Wörner, U., Zimmermann-Timm, H., and Kausch, H. (2000). Succession of protists on estuarine aggregates. Microb. Ecol. 40, 209–222. doi: 10.1007/s002480000051
Zihao, Z., Jihua, H., and Nianpeng, H. (2021). Predictability of Functional Diversity Depends on the Number of Traits. J. Resour. Ecol. 12, 332–345. doi: 10.5814/j.issn.1674-764x.2021.03.003
Keywords: eukaryotic, single-celled, freshwater, soil, phagotrophic, functional traits, protozoa
Citation: Gulin V, Vlaičević B, Sertić Perić M, Rebrina F and Matoničkin Kepčija R (2022) Taxonomic and Functional Metrics of Ciliates and Amoeboid Protists in Response to Stream Revitalization. Front. Microbiol. 13:842395. doi: 10.3389/fmicb.2022.842395
Received: 23 December 2021; Accepted: 08 March 2022;
Published: 01 April 2022.
Edited by:
Per Juel Hansen, University of Copenhagen, DenmarkReviewed by:
Ute Risse-Buhl, Helmholtz Centre for Environmental Research, Helmholtz Association of German Research Centres (HZ), GermanyIndranil Mukherjee, Academy of Sciences of the Czech Republic (ASCR), Czechia
George B. Mcmanus, University of Connecticut Groton, United States
Anja Scherwaß, University of Cologne, Germany
Copyright © 2022 Gulin, Vlaičević, Sertić Perić, Rebrina and Matoničkin Kepčija. 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: Vesna Gulin, vesna.gulin@biol.pmf.hr