
==== Front
PeerJ
PeerJ
peerj
PeerJ
2167-8359
PeerJ Inc. San Diego, USA

17912
10.7717/peerj.17912
Biodiversity
Ecology
Microbiology
Freshwater Biology
Insights from single-strain and mixed culture experiments on the effects of heatwaves on freshwater flagellates
Boden Lisa lisa.boden@uni-due.de
1
Klagus Chantal 1
Boenigk Jens 12
1 Department Biodiversity, University of Duisburg–Essen, Essen, North Rhine Westphalia, Germany
2 Center for Water and Environmental Research, University of Duisburg–Essen, Essen, North Rhine Westphalia, Germany
Kormas Konstantinos
12 9 2024
2024
12 e1791220 6 2024
23 7 2024
©2024 Boden et al.
2024
Boden et al.
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

The increasing frequency and intensity of heatwaves driven by climate change significantly impact microbial communities in freshwater habitats, particularly eukaryotic microorganisms. Heterotrophic nanoflagellates are important bacterivorous grazers and play a crucial role in aquatic food webs, influencing the morphological and taxonomic structure of bacterial communities. This study investigates the responses of three flagellate taxa to heatwave conditions through single-strain and mixed culture experiments, highlighting the impact of both biotic and abiotic factors on functional redundancy between morphologically similar protist species under thermal stress. Our results indicate that temperature can significantly impact growth and community composition. However, density-dependent factors also had a significant impact. In sum, stabilizing effects due to functional redundancy may be pronounced as long as density-dependent factors play a minor role and can be overshadowed when flagellate abundances increase.

Chrysophyceae
Fluorescence in situ hybridization
Heatwave
Heterotrophic flagellate
Microbial food web
Abiotic stressors
Biotic interactions
Collaborative Research Center (CRC) RESISTGerman Research Foundation (DFG)- CRC1439/1 INST 20876/402-1 This study was performed within the Collaborative Research Center (CRC) RESIST and analyses were mainly done in Project A06, funded by the German Research Foundation (DFG) –CRC 1439/1; project number INST 20876/402–1. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
==== Body
pmcIntroduction

Microbial communities frequently experience and cope with natural fluctuations in abiotic environmental conditions (Ruokolainen et al., 2009). However, strong changes caused by climate change and anthropogenic activities put ecosystems under severe pressure (IPBES, 2019). Microorganisms have evolved mechanisms to withstand changing environmental conditions (Bernhardt et al., 2020), but tolerances to abiotic stressor vary strongly (Díaz-Almeyda et al., 2017; Mascarin et al., 2018; Jain & Saraf, 2021). Consequently, abiotic disturbances can result in lower abundance of individual species in the microbial community or total loss of species (Schulhof et al., 2020; Proesmans et al., 2022). Such changes in biodiversity can significantly impact ecosystem functions and food webs (Mooney et al., 2009; IPBES, 2019). However, the impact of losing individual species can be lowered by the presence of other species in the community that perform the same function, an attribute known as functional redundancy (Walker, 1992; Naeem, 1998). Functional redundancy is hypothesized to enhance the resilience and stability of ecosystems during the occurrence of stressors by maintaining ecological functions despite changes in taxon composition (Biggs et al., 2020; Chen et al., 2022). Here we address functional redundancy between cryptic species and its changing significance in the light of density-dependent and density-independent factors—an aspect about which, to our knowledge, little is known, particularly in freshwater habitats.

As a major abiotic stressor, increased temperatures threaten biodiversity across different habitats (Urrutia-Cordero et al., 2017; Maire et al., 2022; Lau et al., 2024). In addition to the long-term increase of the average global surface temperature, more frequent and intense heatwaves have been predicted as a consequence of climate change (Woolway et al., 2021; IPCC, 2021). Temperature thresholds and the duration of heatwaves are region-specific and cannot be universally defined (Xu et al., 2016). However, heatwaves are generally characterized by relatively short but intense increases in temperature, representing one of the most significant challenges in aquatic habitats (Sun & Arnott, 2022). Climate change scenarios predict that temperatures during heatwaves in freshwater bodies can increase by up to 4 °C for several days (Woolway et al., 2021). During the occurrence of a heatwave, organisms are rapidly pushed beyond the limit of their temperature tolerances, which may affect community composition and ecosystem functions more profoundly than a gradual long-term temperature increase (Vasseur et al., 2014; Stillman, 2019). Heatwaves have been shown to cause biodiversity loss in both marine and freshwater habitats (Brauko et al., 2020; Sabater et al., 2022). Particularly for eukaryotic microorganisms, heat stress has resulted in significant shifts in community composition (Hao et al., 2018; Thomson & Manoylov, 2019). Interestingly, it is assumed that short-term disturbances like heatwaves are succeeded by complete or partial recovery (Bender, Case & Gilpin, 1984; Harris et al., 2018).

Colorless heterotrophic nanoflagellates play an important role in regulating the quantity and biomass of bacterial populations and profoundly influence the morphological and taxonomic structure of bacterial communities (Sherr & Sherr, 2002). Consequently, alterations in both the quantity and quality of bacterivory, resulting from shifts in the flagellate community, can influence the overall structure of aquatic food webs. The Chyrosphytes are among the most important bacterivorous grazers in freshwater habitats (Finlay & Esteban, 1998). They are widely distributed across various habitats, primarily inhabiting freshwaters but also extending into terrestrial and marine ecosystems (Kristiansen & Škaloud, 2017). They further encompass a broad range of different morphologies (Škaloud, Kristiansen & Škaloudová, 2013; Skaloudova & Skaloud, 2013). The morphology of flagellates has been shown to correlate with various aspects of predator–prey interactions, including preferences for food size and feeding mechanisms (Boenigk & Arndt, 2000a; Boenigk & Arndt, 2000b). Consequently, flagellates exhibiting similar morphologies can feed on the same prey, thereby serving equivalent roles in the food web. Interestingly, molecular analysis demonstrated high molecular diversity behind distinct morphological forms (Grossmann et al., 2016). This cryptic diversity likely resulted from parallel evolution, wherein similar morphological forms evolved independently in different lineages (Graupner et al., 2018). The ecological significance of this phenomenon, however, remains unclear. In the order Ochromonadales, significant differences in temperature tolerance were observed among clades with similar morphologies (Boenigk et al., 2007; Nolte et al., 2010). Such clade-specific responses to changing environmental conditions can result in a shift of the community composition (Boden, Sieber & Boenigk, 2023). This suggest that cryptic diversity might disguise a turn-over of cryptic taxa that are differently adapted to abiotic factors, thereby stabilizing the structure of food webs and ecosystem functions during the occurrence of stressors despite significant changes in the taxon composition (Boden, Sieber & Boenigk, 2023).

Here, we investigate the effects of heatwaves on the growth of three cryptic taxa associated with Pedospumella, Spumella, and Poteriospumella to gain a better understanding of the ecological importance of cryptic diversity and its potential impact on the stability of ecosystems when facing abiotic disturbances. Given the variation in temperature tolerance, we anticipate that the three strains will respond differently to the temperature increase. Additionally, we expect the strains to grow differently in mixed cultures compared to the single-strain cultures because biotic interactions between the strains will influence the growth of each taxon. We predict more pronounced fluctuations in the total abundance of flagellates in single-strain cultures compared to mixed cultures because the presence of multiple taxa with varying temperature tolerances in mixed cultures may buffer the effects of the heatwave. Consequently, we also expect less fluctuation in the total abundance of prey bacteria in mixed cultures, as the presence of multiple flagellate strains feeding on the same prey should stabilize the predator–prey interaction and ensure functional redundancy in the face of environmental changes and potential shifts in the flagellate community.

Materials & Methods

Strains & cultivation

The strains Pedospumella encystans JBM/S11, Spumella rivalis AR4A6, and Poteriospumella lacustris JBM10 were previously isolated from soil or freshwater samples originating from different geographical locations (Table 1). The axenic strain JBM10 is routinely cultivated in NSY medium (nutrient broth, peptone from soybean, yeast extract; Hahn et al., 2003) in 100 ml Erlenmeyer flasks. Due to the lack of axenic strains associated with Peodospumella and Spumella, two xenic strains were used in this study. Both JBM/S11 and AR4A6 are routinely grown in inorganic basal medium (Hahn et al., 2003) in cell culture flasks (25 cm3 with filter screw cap; TTP Techno Plastic Products AG) using wheat grains as a food source. The compositions of the NSY medium and the IB medium are provided in Table S1. All strains are grown in a climate chamber (SANYO Electric Co. Ltd., Osaka, Japan) at 15 °C under a 14h:10 h light-dark cycle. The bacterial strain Linmohabitans spp. IID5 is routinely grown as an axenic culture in NSY medium in 100 ml Erlenmeyer flasks at room temperature (RT) with constant shaking (96 rpm, orbital shaker; LAUDA-Brinkmann, LP, Marlton, NJ, USA).

10.7717/peerj.17912/table-1 Table 1 Origin and affiliation of strains used in this study.

Strain	Species designation	18S clade	Geographical Origin	Habitat origin	Media	
JBMS11	Pedospumella encystans	C1	Austria, Mondsee, near Rauchhaus	soil	IB + grain of wheat	
AR4A6	Spumella rivalis	C2	Austria, River Fuschler Ache	freshwater	IB + grain of wheat	
JBM10	Poteriospumella lacustris	C3	Austria, Lake Mondsee	freshwater	NSY (3g/L)	

Experimental set-up with single-strain and mixed cultures

To create similar conditions for all three strains, all experiments were performed in IB medium. Therefore, flagellates were harvested by centrifugation for 10 min at 2,820 g at RT to remove the original medium and each resulting pellet was then resuspended in 50 ml IB medium. One ml aliquots of a Linmohabitans spp. IID5 culture were centrifuged for 10 min at 15,000 g at RT and, after removal of the supernatant, resuspended in IB medium. One milliliter of the washed bacteria was added to each flagellate culture as food source. These pre-cultures were then incubated at 23 °C with a 14h:10 h light-dark cycle for three days to acclimatize to the new conditions. The bacteria from both xenic strains were isolated through filtration (diameter 25 mm, pore size 0.45 µm; Millipore GTTP 02500, Millipore, Eschborn, Germany) and cultivated in NSY medium in a 100 ml Erlenmeyer flask at 15 °C with a 14h:10 h light-dark cycle. All experiments were conducted in triplicates in 75 ml of IB medium in 100 ml Erlenmeyer flasks. Total flagellate abundance in the pre-cultures at the end of the acclimatization phase was counted in Sedgewick–Rafter chambers using 1 ml of the Lugol–fixed subsamples to adjust the flagellate abundance to approximately 30,000 cells per ml in the single-strain cultures. In the mixed cultures, the three strains were mixed to achieve relative abundances of 20% Pedospumella, 30% Poteriospumella, and 50% Spumella, mimicking the relative abundances of the respective clades in natural environments during late spring (Weisse, 1997; Cleven & Weisse, 2001; Nolte et al., 2010). The concentration of food bacteria was set to 2 × 107 bacteria per ml, consisting predominately of Limnohabitans spp. IID5 in addition to 10,000 cells per ml of the previously cultivated bacteria from the xenic flagellate cultures to ensure uniformed bacterial backgrounds in all cultures. Further, wheat grains were preheated at 70 °C for one hour and one wheat grain was added to each culture as food source for the bacteria. All cultures were initially incubated at 23 °C with a 14h:10 h light-dark cycle for two days. The temperature was then increased to 27 °C for three days before being reduced back to 23 °C for a two-day recovery phase. While 23 °C provides optimal growth conditions for representatives of all three clades, it has been shown that growth rates generally decrease at temperatures above 25 °C, particularly for strains associated with Pedospumella and Spumella (Boenigk et al., 2007). The 4 °C temperature increase exposed the cultures to conditions that induced heat stress, while also ensuring the survival of the cultures. Four ml subsamples were collected daily and fixed with 2% paraformaldehyde (PFA) for 1 h at RT or overnight at 4 °C and filtered onto white polycarbonate filters (diameter 25 mm, pore size 0.2 µm; Millipore GTTP 02500; Millipore, Eschborn, Germany). The filters were air-dried and stored at −20 °C for further use.

Determination of cell counts and growth rates in single-strain and mixed cultures

Sections of the filters were cut out with a razor blade. To determine total flagellate abundance in single-strain cultures, the filter sections were stained with 4,6-diamidino-2-phenylindole (DAPI; 0.1 mg/ml) for 10 min at RT. To remove non-specific staining, the filters were washed in sterile distilled water for 1 min at RT, then dehydrated with 100% ethanol for 30 s at RT and air-dried. To determine the fraction of each strain in the mixed cultures, Pedospumella encystans, Spumella rivalis, and Poteriospumella lacustris were stained using fluorescently labelled clade-specific probes (O1C531, O2C613, O3C723; Boden, Sieber & Boenigk, 2023; 50 ng DNA/µL; Eurofins Genomics Germany GmbH, Ebersberg, Germany). Therefore, filter sections were incubated sequentially in 50%, 80%, and 100% ethanol solutions to dehydrate the samples. For each dehydration step, the filters were incubated for 3 min at RT. Then, the cells were hybridized with Cy3-labeled single-strand oligonucleotide probes as described by Glöckner et al. (1996). To optimize signal detection rates, the concentration of SDS was increased to 0.02% in the hybridization buffer, and the hybridization time was extended to 3 h at 46 °C. To determine total flagellate abundance in the mixed cultures, the filter sections were afterward counterstained with DAPI as described before. After staining, the filter sections were mounted in the non-hardening and non-bleaching mounting medium CitiFluorTM AF2 (Citifluor, Ltd., London, United Kingdom), and images were captured using the NIS-Elements BR imaging software (Nikon Corp., Tokyo, Japan) and a Nikon Eclipse 80i microscope (Nikon Corp., Tokyo, Japan) with a 100x objective and appropriate fluorescence filters to detect DAPI and Cy3 signals. DAPI and Cy3 signals were counted manually using the software package Fiji (Schindelin et al., 2012).

Data analysis

Growth rates (µ) were calculated for three distinct phases throughout the experiment using R (v4.3.3; R Core Team, 2024) and RStudio (v2023.12.1; RStudio Team, 2024) and the formula μ=lnCt− lnC0/t

where Ct is the total flagellate abundances at the end of the respective interval and C0 is the flagellate abundance at the beginning of the interval. All growth rates indicate growth per day. Specifically, growth rates were determined during the initial acclimatization phase from exponentially growing single strains as well as for all flagellates in the mixed cultures. Further, growth rates were assessed at the end of the heatwave, i.e., between days five and six, and during the recovery phase, i.e., between days seven and eight. Cell abundances, growth rates and the relative share of each taxa in the mixed cultures were plotted using R, RStudio and the R packages “readxl”, “gglot”, “cowplot” and “dplyr” (Wickham & Bryan, 2023; Wickham, 2016; Wilke, 2024; Wickham et al., 2023). Statistical data analysis was performed in R and RStudio. Growth rates were statistically compared among each other via Welch’s ANOVA and Games-Howell post-hoc test using the R-packages “readxl”, “dplyr” and “rstatix” (Wickham & Bryan, 2023; Wickham et al., 2023; Kassambara, 2023). Similarly, the relative shares of Pedospumella encystans, Spumella rivalis and Poteriospumella lacustris in the mixed cultures at different sampling days were statistically compared via Welch’s ANOVA and Games-Howell post-hoc test. Variability in bacterial abundances was statistically compared among all cultures via Brown-Forsythe test.

Results

After a lag phase between day one and day two, an increase in total flagellate abundance starting from day two was observed in all cultures (Fig. 1). Growth occurred between day two and day four in all cultures, with the exception of Poteriospumella single-strain experiments, where growth ceased by day three (Fig. 1). During this initial growth phase, hereafter referred to as “heatwave onset”, significant differences among the cultures were observed when comparing flagellate growth in the single-strain cultures and in the mixed cultures (Welch’s ANOVA; p = 0.04). Growth rates of the taxa were not significantly different when grown alone (Welch’s ANOVA; p = 1.00; Fig. 2). But growth rates significantly differed during the heatwave onset (Welch’s ANOVA; p = 0.004) between flagellates growing in single-strain culture and the same taxon growing in mixed culture: Growth rates for Spumella (Games-Howell test; p = 0.02) and Pedospumella (Games-Howell test; p = 0.02) were significantly lower in the mixed cultures compared to when they were grown individually. In contrast, growth rates of Poteriospumella did not differ between treatments (Games-Howell test; p = 1.00; Fig. 2). At the end of the heatwave total flagellate abundance decreased in all cultures (Fig. 1) and no significant differences were observed between the three taxa (Welch’s ANOVA, p = 0.98; Fig. 2). Similarly, no significant differences were observed when comparing growth of the three taxa in single-strain culture and in mixed culture (Games-Howell test; Pedospumella: p = 1.00, Spumella: p = 1.00, Poteriospumella: p = 1.00). During the recovery phase no significant differences were observed in total flagellate growth when comparing single-strain experiments and mixed cultures (Welch’s ANOVA; p = 0.37). Interestingly, the growth rates of Poteriospumella and Spumella in single-strain cultures were slightly but not significantly lower than in the mixed cultures (Games-Howell test; Poteriospumella: p = 0.68, Spumella: p = 0.64; Fig. 2). Replicates of the same treatment exhibited similar trends throughout the entire experiment (Fig. 2).

10.7717/peerj.17912/fig-1 Figure 1 Cell abundances under varying temperature conditions.

Total absolute abundance of flagellates (left) and bacteria (right) during phases of normal (blue) and increased (orange) temperature in single-stain cultures of (A) Pedospumella encystans, (B) Spumella rivalis, (C) Poteriospumella lacustris and (D) mixed cultures containing all three strains.

10.7717/peerj.17912/fig-2 Figure 2 Flagellate growth in single-strain and mixed cultures.

Growth rates [day−1] of Pedospumella encystans, Spumella rivalis and Poteriospumella lacustris in single-strain cultures and in mixed cultures as well as total flagellate growth in the mixed cultures during the onset of the heatwave, at the end of the heatwave and during the recovery phase. Negative growth rates correspond to mortality rates.

The relative share of all three taxa in the mixed culture (Fig. 3) changed significantly during the onset of the heatwave but was stable before and after: No significant changes were observed for any of the three taxa during the initial acclimatization phase, i.e., between day one and day two (Games-Howell test; Pedospumella: p = 1.00, Spumella: p = 0.20, Poteriospumella: p = 0.48). During the onset of the heatwave, i.e., between day two and three, the relative share of Spumella decreased significantly from 52% on day two to 34% on day three (Games-Howell test; p < 0.001). Similarly, the relative share of Pedospumella decreased significantly from 21% on day two to 11% on day three (Games-Howell test; p < 0.001). In contrast, the share of Poteriospumella increased significantly from 29% on day two to 55%, dominating the community on day three (Games-Howell test; p < 0.001). After this shift in community composition during onset of the heatwave, no further significant changes of the relative share were observed both during the heatwave (Games-Howell test; Pedospumella: p = 0.97, Spumella: p = 0.99, Poteriospumella: p = 0.96) and during the following recovery phase (Games-Howell test; Pedospumella: p = 0.95, Spumella: p = 1.00, Poteriospumella: p = 1.00). On day seven, Poteriospumella still dominated the community, constituting on average 56% of the total flagellate abundance, while Pedospumella and Spumella represented 12% and 34%, respectively, of all flagellates in the mixed cultures.

10.7717/peerj.17912/fig-3 Figure 3 Community composition in mixed cultures.

Relative shares of Pedospumella encystans, Poteriospumella lacustris and Spumella rivalis in the mixed cultures. The bars indicate the mean percentage of each strain in all mixed cultures on day 1, day 2, day 3, day 4, and day 7. Error bars represent standard deviations.

Furthermore, variability of the bacterial abundances over time differed significantly among the cultures (Brown-Forsythe test; p = 0.01). The relative variability in bacterial abundances was assessed by calculating the coefficient of variation (CV) for each culture. Pedospumella single-strain cultures exhibited the lowest CV at 39.1%, indicating the least variability in bacterial abundance throughout the experiment. In contrast, Spumella single-strain cultures showed the highest CV at 58.8%. Poteriospumella single-strain cultures and mixed cultures displayed similar CVs of 49.4% and 49.6%, respectively (Fig. S1).

Discussion

The impact of temperature shifts, e.g., due to global warming and heatwaves, has been the subject of numerous studies (Font et al., 2021; Woolway et al., 2021), revealing substantial shifts in protist communities (Hao et al., 2018; Thomson & Manoylov, 2019). We observed significant differences both among the different cultures and when comparing growth of each taxa in single-strain and mixed cultures (Fig. 2). While total flagellate growth was generally higher in single-strain experiments than in the mixed cultures, Poteriospumella exhibited the highest growth rate among all cultures. Interestingly, similar growth rates were observed for Poteriospumella in both single-strain and mixed cultures, contrasting with Pedospumella and Spumella, which displayed significantly lower growth rates in mixed cultures compared to when they were grown individually (Fig. 2). This is in accordance with previous research that identified Poteriospumella’s competitive advantage, particularly under elevated temperature conditions (Boden, Sieber & Boenigk, 2023). Additionally, we observed a shift in the community composition in the mixed cultures. At the beginning of the experiment, the community was dominated by Spumella but its relative share decreased significantly during the onset of the heatwave. In contrast, the relative abundance of Poteriospumella increased significantly, becoming the dominant taxa approximately 24 h after the temperature rise (Fig. 3). This shift in taxon composition in our mixed cultures is in accordance with seasonal environmental data that identified similar temperature-induced shifts in a natural lake, with Spumella typically dominating in winter and Poteriospumella in summer (Nolte et al., 2010). No shifts in community composition were observed in the mixed cultures at the end of the heatwave and during the recovery phase, i.e., the community composition remained unchanged despite changing temperature conditions (Fig. 3). This was unexpected as the influence of abiotic factors such as temperature, salinity or oxygen availability on flagellate abundance and community composition has been demonstrated in previous studies (Nolte et al., 2010; Princiotta & Sanders, 2017; Gran-Stadniczeñko et al., 2018). However, most of these studies were investigating the effects of abiotic factor in environmental samples, indicating that mock communities in laboratory experiments may be less dynamic than natural communities, possibly due to reduced competition in species-poor communities. At the end of the heatwave, the growth rates of Poteriospumella, both in single-strain and mixed cultures, were similar to those of Pedospumella and Spumella and no competitive advantage was identified for any of the three taxa during that later phase. The reasons why Poteriospumella failed to maintain its competitive advantage throughout the entire duration of the heatwave remain unclear. All cultures reached maximal abundance latest by day five, therefore potential reasons include density-dependent factors like the accumulation of toxic metabolic byproducts such as acids and alcohols, or oxygen depletion, which are common phenomena in cultures with high cell density and would affect all three taxa equally (Rouf et al., 2017). Additionally, the daily sampling resulted in a reduction of medium in all cultures which could further impact flagellate growth particularly during the latter half of the experiment. Beside abiotic factors, biotic interactions can also influence flagellate growth in the cultures. High cell densities can also lead to cannibalism among flagellates, particularly when alternative prey is scarce or competition is high (Martel & Flynn, 2008; Yang et al., 2020). However, we found no indications of cannibalism in our samples, likely because the concentration of prey bacteria remained above 5 ×106 cells per milliliter throughout the experiment, which was sufficiently high to prevent starvation (Boenigk et al., 2002). Given the consistently high concentration of prey bacteria, it is unlikely that nutrient depletion, a common cause of mortality in dense microbial cultures (Rouf et al., 2017), was responsible for the decline in total flagellate abundances at the end of the heatwave. Our results indicate that the abiotic and/or biotic density-dependent factors arising from the experimental setup overshadowed the effects of the temperature increase during the latter half of the heatwave. While heterotrophic nanoflagellates can reach abundances of several thousand cells per milliliter in natural habitats, their total flagellate abundances remain notably lower than the maximum abundances that were reached in this study (Boenigk & Arndt, 2002). Therefore, the response to the temperature increase that was observed during the onset of the heatwave may be a closer representation of natural habitats than towards the end of the heatwave. In the recovery phase, growth rates of Spumella and Poteriospumella were slightly but not significantly higher in the mixed cultures than in single-strain cultures, indicating that the lower total flagellate abundance during the recovery phase may reduce the impact of density-dependent factors, thereby allowing temperature-related effects to become dominant again.

As important players in the top-down control of bacterial populations, changes in both the total abundance and community composition of heterotrophic nanoflagellates significantly influence the dynamics of bacterial prey populations (Sherr & Sherr, 2002). In this study, variability in the bacterial abundances differed among the single strain-cultures, with the highest variation observed in Spumella single-strain cultures. Notably, within mixed cultures, bacterial abundances displayed similar variability than in the Poteriospumella single-strain cultures. The strong fluctuations observed in Spumella single-strain cultures were likely mitigated in the mixed cultures by the presence and dominance of Poteriospumella. These results further support the hypothesis that functional redundancy, demonstrated here by the simultaneous presence of taxa capable of feeding on the same prey but exhibiting different adaptations to abiotic factors, can stabilize habitats when confronted with changing environmental conditions. Interestingly, variability of bacterial abundance was lowest in Pedospumella single-strain cultures. However, Pedospumella showed low competitive strength in the presence of the other two taxa, thus failing to further stabilize bacterial populations in the mixed cultures. This highlights the importance of biotic interactions when assessing the impact of extreme climate events on natural microbial communities.

Conclusions

Our findings highlight the complexities of protist community dynamics under thermal stress, indicating that while temperature can drive significant changes in abundance and composition of flagellate communities, other factors such as density-dependent interactions and prey availability also significantly impact microbial communities. Our study indicates that stabilizing effects due to functional redundancy may be pronounced as long as density-dependent factors play a minor role. This study underscores the importance of considering both abiotic and biotic factors when assessing the impacts of climate change on microbial communities, as well as the potential for functional redundancy to stabilize ecosystems facing environmental fluctuations.

Supplemental Information

10.7717/peerj.17912/supp-1 Supplemental Information 1 Compositions of the IB and NSY media

10.7717/peerj.17912/supp-2 Supplemental Information 2 Variability of bacterial abundance

Coefficient of variation calculated from the bacterial abundances in single-strain cultures of Pedospumella encystans, Spumella rivalis, Poteriospumella lacustris and mixed cultures containing all three strains.

10.7717/peerj.17912/supp-3 Supplemental Information 3 Raw count data

The raw cell counts that were used to calculate growth rates and investigate community composition in the mixed cultures.

Additional Information and Declarations

Competing Interests

Author Contributions

Data Availability

The authors declare there are no competing interests.

Lisa Boden performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Chantal Klagus performed the experiments, analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

Jens Boenigk conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

The following information was supplied regarding data availability:

Raw data are available in the Supplemental Files.
==== Refs
References

Bender, Case & Gilpin (1984) Bender EA Case T Gilpin ME 1984 Perturbation experiments in community ecology: theory and practice Ecology 65 1 1 13 10.2307/1939452
Bernhardt et al. (2020) Bernhardt JR O’Connor MI Sunday JM Gonzalez A 2020 Life in fluctuating environments Philosophical Transactions of the Royal Society B: Biological Sciences 375 1814 20190454 10.1098/rstb.2019.0454
Biggs et al. (2020) Biggs CR Yeager LA Bolser DG Bonsell C Dichiera AM Hou Z Keyser SR Khursigara AJ Lu K Muth AF Negrete B Erisman BE 2020 Does functional redundancy affect ecological stability and resilience? A review and meta-analysis Ecosphere 11 7 e03184 10.1002/ecs2.3184
Boden, Sieber & Boenigk (2023) Boden L Sieber G Boenigk J 2023 Effects of stressors on growth and competition between different cryptic taxa affiliated with Ochromonadales (Chrysophyceae) Fottea 23 2 235 245 10.5507/fot.2023.003
Boenigk & Arndt (2000a) Boenigk J Arndt H 2000a Particle handling during interception feeding by four species of heterotrophic nanoflagellates Journal of Eukaryotic Microbiology 47 4 350 358 10.1111/j.1550-7408.2000.tb00060.x 11140448
Boenigk & Arndt (2000b) Boenigk J Arndt H 2000b Comparative studies on the feeding behavior of two heterotrophic nanoflagellates: the filter-feeding choanoflagellate Monosiga ovata and the raptorial-feeding kinetoplastid Rhynchomonas nasuta Aquatic Microbial Ecology 22 243 249 10.3354/ame022243
Boenigk & Arndt (2002) Boenigk J Arndt H 2002 Bacterivory by heterotrophic flagellates: community structure and feeding strategies Antonie Van Leeuwenhoek 81 1-4 465 80 10.1023/a:1020509305868 12448743
Boenigk et al. (2007) Boenigk J Jost S Stoeck T Garstecki T 2007 Differential thermal adaptation of clonal strains of a protist morphospecies originating from different climatic zones Environmental Microbiology 9 3 593 602 10.1111/j.1462-2920.2006.01175.x 17298360
Boenigk et al. (2002) Boenigk J Matz C Jürgens K Arndt H 2002 Food concentration-dependent regulation of food selectivity of interception-feeding bacterivorous nanoflagellates Aquatic Microbial Ecology 27 2 195 202 10.3354/ame027195
Brauko et al. (2020) Brauko KM Cabral A Costa NV Hayden J Dias CEP Leite ES Westphal RD Mueller CM Hall-Spencer JM Rodrigues RR Rörig LR Pagliosa PR Fonseca AL Alarcon OE Horta PA 2020 Marine heatwaves, sewage and eutrophication combine to trigger deoxygenation and biodiversity loss: a SW Atlantic case study Frontiers in Marine Science 7 590258 10.3389/fmars.2020.590258
Chen et al. (2022) Chen H Ma K Lu C Fu Q Qiu Y Zhao J Huang Y Yang Y Schadt CW Chen H 2022 Functional redundancy in soil microbial community based on metagenomics across the globe Frontiers in Microbiology 13 878978 10.3389/fmicb.2022.878978 35586865
Cleven & Weisse (2001) Cleven E Weisse T 2001 Seasonal succession and taxon-specific bacterial grazing rates of heterotrophic nanoflagellates in Lake Constance Aquatic Microbial Ecology 23 147 161 10.3354/AME023147
Díaz-Almeyda et al. (2017) Díaz-Almeyda EM Prada C Ohdera AH Moran H Civitello DJ Iglesias-Prieto R Carlo TA LaJeunesse TC Medina M 2017 Intraspecific and interspecific variation in thermotolerance and photoacclimation in Symbiodinium dinoflagellates Proceedings of the Royal Society B: Biological Sciences 284 1868 20171767 10.1098/rspb.2017.1767
Finlay & Esteban (1998) Finlay BJ Esteban GF 1998 Freshwater protozoa: biodiversity and ecological function Biodiversity and Conservation 7 9 1163 1186 10.1023/a:1008879616066
Font et al. (2021) Font RA Khamis K Milner AM Smith GHS Ledger ME 2021 Low flow and heatwaves alter ecosystem functioning in a stream mesocosm experiment Science of the Total Environment 777 146067 10.1016/j.scitotenv.2021.146067 33677285
Glöckner et al. (1996) Glöckner FO Amann R Alfreider A Pernthaler J Psenner R Trebesius K Schleifer KH 1996 An in situ hybridization protocol for detection and identification of planktonic bacteria Systematic and Applied Microbiology 19 3 403 406 10.1016/s0723-2020(96)80069-5
Gran-Stadniczeñko et al. (2018) Gran-Stadniczeñko S Egge E Hostyeva V Logares R Eikrem W Edvardsen B 2018 Protist diversity and seasonal dynamics in Skagerrak Plankton communities as revealed by metabarcoding and microscopy Journal of Eukaryotic Microbiology 66 3 494 513 10.1111/jeu.12700 30414334
Graupner et al. (2018) Graupner N Jensen M Bock C Marks S Rahmann S Beisser D Boenigk J 2018 Evolution of heterotrophy in chrysophytes as reflected by comparative transcriptomics FEMS Microbiology Ecology 94 4 fiy039 10.1093/femsec/fiy039 29518196
Grossmann et al. (2016) Grossmann L Bock C Schweikert M Boenigk J 2016 Small but manifold –hidden diversity in Spumella -like Flagellates Journal of Eukaryotic Microbiology 63 4 419 439 10.1111/jeu.12287 26662881
Hahn et al. (2003) Hahn MW Lünsdorf H Wu Q Schauer M Höfle MG Boenigk J Stadler P 2003 Isolation of novel ultramicrobacteria classified as actinobacteria from five freshwater habitats in Europe and Asia Applied and Environmental Microbiology 69 3 1442 1451 10.1128/aem.69.3.1442-1451.2003 12620827
Hao et al. (2018) Hao B Roejkjaer AF Wu H Cao Y Jeppesen E Li W 2018 Responses of primary producers in shallow lakes to elevated temperature: a mesocosm experiment during the growing season of Potamogeton crispus Aquatic Sciences 80 4 34 10.1007/s00027-018-0585-0
Harris et al. (2018) Harris RMB Beaumont LJ Vance TR Tozer CR Remenyi TA Perkins-Kirkpatrick SE Mitchell PJ Nicotra AB McGregor S Andrew NR Letnic M Kearney MR Wernberg T Hutley LB Chambers LE Fletcher MS Keatley MR Woodward CA Williamson G Duke NC Bowman DMJS 2018 Biological responses to the press and pulse of climate trends and extreme events Nature Climate Change 8 7 579 587 10.1038/s41558-018-0187-9
IPBES (2019) IPBES 2019 Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services IPBES secretariat Bonn Brondizio ES Settele J Díaz S Ngo HT 10.5281/zenodo.3831673
IPCC (2021) IPCC 2021 Summary for policymakers Climate Change 2021: The physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change https://www.ipcc.ch/report/ar6/wg1/
Jain & Saraf (2021) Jain R Saraf M 2021 Exploring the abiotic and biotic stress tolerance potential of Rhizobactera isolated from Cyamosis Journal of Advanced Scientific Research 12 03 190 194 10.55218/jasr.202112327
Kassambara (2023) Kassambara A 2023 rstatix: pipe-friendly framework for basic statistical tests R package version 0.7.2 https://CRAN.R-project.org/package=rstatix
Kristiansen & Škaloud (2017) Kristiansen J Škaloud P 2017 Chrysophyta Handbook of the Protists Springer International Publishing AG Cham, Switzerland Archibald JM Simpson AGB Slamovits CH 331 366
Lau et al. (2024) Lau GXH Fan HY Halim MA Cheah YK Wong CMVL 2024 Forecasting the impact of global warming on soil bacterial communities using simulated systems Malaysian Journal of Microbiology 20 3 397 405 10.21161/mjm.230323
Maire et al. (2022) Maire J Buerger P Chan W Deore P Dungan A Nitschke M Oppen M 2022 Effects of Ocean warming on the underexplored members of the coral microbiome Integrative and Comparative Biology 62 1700 1709 10.1093/icb/icac005 35259253
Martel & Flynn (2008) Martel C Flynn K 2008 Morphological controls on cannibalism in a planktonic marine phagotroph Protist 159 1 41 51 10.1016/J.PROTIS.2007.05.003 17768088
Mascarin et al. (2018) Mascarin G Pereira-Junior R Fernandes É Quintela E Dunlap C Arthurs S 2018 Phenotype responses to abiotic stresses, asexual reproduction and virulence among isolates of the entomopathogenic fungus Cordyceps javanica (Hypocreales: Cordycipitaceae) Microbiological Research 216 12 22 10.1016/j.micres.2018.08.002 30269851
Mooney et al. (2009) Mooney H Larigauderie A Cesario M Elmquist T Hoegh-Guldberg O Lavorel S Mace GM Palmer M Scholes R Yahara T 2009 Biodiversity, climate change, and ecosystem services Current Opinion in Environmental Sustainability 1 1 46 54 10.1016/j.cosust.2009.07.006
Naeem (1998) Naeem S 1998 Species redundancy and ecosystem reliability Conservation Biology 12 39 45 10.1111/j.1523-1739.1998.96379.x
Nolte et al. (2010) Nolte V Pandey RV Jost S Medinger R Ottenwälder B Boenigk J Schlötterer C 2010 Contrasting seasonal niche separation between rare and abundant taxa conceals the extent of protist diversity Molecular Ecology 19 14 2908 2915 10.1111/j.1365-294x.2010.04669.x 20609083
Princiotta & Sanders (2017) Princiotta SD Sanders RW 2017 Heterotrophic and mixotrophic nanoflagellates in a mesotrophic lake: abundance and grazing impacts across season and depth Limnology and Oceanography 62 2 632 644 10.1002/lno.10450
Proesmans et al. (2022) Proesmans W Andrews C Gray A Griffiths R Keith A Nielsen UN Spurgeon D Pywell R Emmett B Vanbergen AJ 2022 Long-term cattle grazing shifts the ecological state of forest soils Ecology and Evolution 12 4 e8786 10.1002/ece3.8786 35386880
R Core Team (2024) R Core Team 2024 R: a language and environment for statistical computing R Foundation for Statistical Computing Vienna https://www.R-project.org/
Rouf et al. (2017) Rouf A Kanojia V Naik HR Naseer B Qadri T 2017 An overview of microbial cell culture Journal of Pharmacognosy and Phytochemistry 6 6 1923 1928
RStudio Team (2024) RStudio Team 2024 RStudio: integrated development environment for R RStudio, PBC Boston http://www.rstudio.com/
Ruokolainen et al. (2009) Ruokolainen L Lindén A Kaitala V Fowler MS 2009 Ecological and evolutionary dynamics under coloured environmental variation Trends in Ecology & Evolution 24 10 555 563 10.1016/j.tree.2009.04.009 19699550
Sabater et al. (2022) Sabater S Freixa A Jiménez L López-Doval J Pace G Pascoal C Perujo N Craven D González-Trujillo JD 2022 Extreme weather events threaten biodiversity and functions of river ecosystems: evidence from a meta-analysis Biological Reviews 98 2 450 461 10.1111/brv.12914 36307907
Schindelin et al. (2012) Schindelin J Arganda-Carreras I Frise E Kaynig V Longair M Pietzsch T Preibisch S Rueden C Saalfeld S Schmid B Tinevez J-Y White DJ Hartenstein V Eliceiri K Tomancak P Cardona A 2012 Fiji: an open-source platform for biological-image analysis Nature Methods 9 7 676 682 10.1038/nmeth.2019 22743772
Schulhof et al. (2020) Schulhof MA Allen AE Allen EE Mladenov N McCrow JP Jones NT Blanton J Cavalheri HB Kaul D Symons CC Shurin JB 2020 Sierra Nevada mountain lake microbial communities are structured by temperature, resources and geographic location Molecular Ecology 29 11 2080 2093 10.1111/mec.15469 32578266
Sherr & Sherr (2002) Sherr EB Sherr BF 2002 Significant of predation by protists in aquatic microbial food webs Antonie Van Leewenhoek 81 293 308 10.1023/A:1020591307260
Skaloudova & Skaloud (2013) Skaloudova M Skaloud P 2013 A new species of Chrysosphaerella (Chrysophyceae: Chromulinales), Chrysosphaerella rotundata sp. nov. from Finland Phytotaxa 130 1 34 10.11646/phytotaxa.130.1.4
Stillman (2019) Stillman JH 2019 Heat waves, the new normal: Summertime temperature extremes will impact animals, ecosystems, and human communities Physiology 34 2 86 100 10.1152/physiol.00040.2018 30724132
Sun & Arnott (2022) Sun X Arnott SE 2022 Interactive effects of increased salinity and heatwaves on freshwater zooplankton communities in simultaneous and sequential treatments Freshwater Biology 67 9 1604 1617 10.1111/fwb.13964
Thomson & Manoylov (2019) Thomson AH Manoylov KM 2019 Algal community dynamics within the Savannah River Estuary, Georgia under anthropogenic stress Estuaries and Coasts 42 6 1459 1474 10.1007/s12237-019-00579-1
Urrutia-Cordero et al. (2017) Urrutia-Cordero P Ekvall MK Ratcovich J Soares M Wilken S Zhang H Hansson LA 2017 Phytoplankton diversity loss along a gradient of future warming and brownification in freshwater mesocosms Freshwater Biology 62 11 1869 1878 10.1111/fwb.13027
Vasseur et al. (2014) Vasseur DA De Long JP Gilbert B Greig HS Harley CDG McCann KS Savage V Tunney TD O’Connor MI 2014 Increased temperature variation poses a greater risk to species than climate warming Proceedings of the Royal Society B: Biological Sciences 281 1779 20132612 10.1098/rspb.2013.2612
Škaloud, Kristiansen & Škaloudová (2013) Škaloud P Kristiansen J Škaloudová M 2013 Developments in the taxonomy of silica-scaled chrysophytes –from morphological and ultrastructural to molecular approaches Nordic Journal of Botany 31 4 385 402 10.1111/j.1756-1051.2013.00119.x
Walker (1992) Walker BH 1992 Biodiversity and ecological redundancy Conservation Biology 6 1 18 23 10.1046/j.1523-1739.1992.610018.x
Weisse (1997) Weisse T 1997 Growth and production of heterotrophic nanoflagellates in a meso-eutrophic lake Journal of Plankton Research 19 703 722 10.1093/PLANKT/19.6.703
Wickham (2016) Wickham 2016 ggplot2: elegant graphics for data analysis Springer-Verlag New York
Wickham & Bryan (2023) Wickham H Bryan J 2023 readxl: read excel files R package version 1.4.3 https://CRAN.R-project.org/package=readxl
Wickham et al. (2023) Wickham H François R Henry L Müller K Vaughan D 2023 dplyr: a grammar of data manipulation https://CRAN.R-project.org/package=dplyr R package version 1.1.4
Wilke (2024) Wilke C 2024 cowplot: streamlined plot theme and plot annotations for ’ggplot2’ R package version 1.1.3 https://CRAN.R-project.org/package=cowplot
Woolway et al. (2021) Woolway R Jennings E Shatwell T Golub M Pierson D Maberly S 2021 Lake heatwaves under climate change Nature 589 402 407 10.1038/s41586-020-03119-1 33473224
Xu et al. (2016) Xu Z FitzGerald G Guo Y Jalaludin B Tong S 2016 Impact of heatwave on mortality under different heatwave definitions: a systematic review and meta-analysis Environment International 89–90 193 203 10.1016/j.envint.2016.02.007
Yang et al. (2020) Yang H Hu Z Shang L Deng Y Tang YZ 2020 A strain of the toxic dinoflagellate Karlodinium veneficum isolated from the East China Sea is an omnivorous phagotroph Harmful Algae 93 101775 10.1016/j.hal.2020.101775 32307067
