Opposing patterns of zooplankton diversity and functioning along a natural stress gradient: when the going gets tough, the tough get going

Contrasting response of microeukaryotic and bacterial communities to the interplay of seasonality and local stressors in shallow soda lakes

Seasonal environmental variation is a leading driver of microbial planktonic community assembly and interactions. However, departures from usual seasonal trends are often reported. To understand the role of local stressors in modifying seasonal succession, we sampled fortnightly, throughout three seasons, five nearby shallow soda lakes exposed to identical seasonal and meteorological changes. We characterised their microeukaryotic and bacterial communities by amplicon sequencing of the 16S and 18S rRNA gene, respectively. Biological interactions were inferred by analyses of synchronous and time-shifted interaction networks, and the keystone taxa of the communities were topologically identified. The lakes showed similar succession patterns during the study period with spring being characterised by the relevance of trophic interactions and a certain level of community stability followed by a more dynamic and variable summer-autumn period. Adaptation to general seasonal changes happened through shared core microbiome of the lakes. Stochastic events such as desiccation disrupted common network attributes and introduced shifts from the prevalent seasonal trajectory. Our results demonstrated that, despite being extreme and highly variable habitats, shallow soda lakes exhibit certain similarities in the seasonality of their planktonic communities, yet local stressors such as droughts instigate deviations from prevalent trends to a greater extent for microeukaryotic than for bacterial communities.

Environmental changes associated with drying climate are expected to affect functional groups of pro- and microeukaryotes differently in temporary saline waters

Temporary ponds are among the most sensitive aquatic habitats to climate change. Their microbial communities have crucial roles in food webs and biogeochemical cycling, yet how their communities are assembled along environmental gradients is still understudied. This study aimed to reveal the environmental drivers of diversity (OTU-based richness, evenness, and phylogenetic diversity) and community composition from a network of saline temporary ponds, soda pans, in two consecutive spring seasons characterized by contrasting weather conditions. We used DNA-based molecular methods to investigate microbial community composition. We tested the effect of environmental variables on the diversity of prokaryotic (Bacteria, Cyanobacteria) and microeukaryotic functional groups (ciliates, heterotrophic flagellates and nanoflagellates, fungi, phytoplankton) within and across the years. Conductivity and the concentration of total suspended solids and phosphorus were the most important environmental variables affecting diversity patterns in all functional groups. Environmental conditions were harsher and they also had a stronger impact on community composition in the dry spring. Our results imply that these conditions, which are becoming more frequent with climate change, have a negative effect on microbial diversity in temporary saline ponds. This eventually might translate into community-level shifts across trophic groups with changing local conditions with implications for ecosystem functioning.

Life-history omnivory in the fairy shrimp Branchinecta orientalis (Branchiopoda: Anostraca)

Very little is known about the feeding of naupliar and juvenile life stages of omnivorous fairy shrimps (Crustacea: Anostraca). Here, we aim to reveal whether the fairy shrimp Branchinecta orientalis is an ontogenetic omnivore and at which age and ontogenetic stage they gain the ability to feed on zooplankton. We assess how food uptake rates change with age until reaching maturity by providing algae (pico- and nanoplanktonic unicellular algae) and zooplankton (rotifers and copepod nauplii) as food in individual experiments. We found that the fairy shrimp B. orientalis started to feed on both types of algal prey immediately after hatching. Nanoplanktonic algae likely represented the most important food source until reaching maturity. Moreover, fairy shrimps started to feed on zooplankton already when they were 7 days old. Slow-moving rotifers gradually gained importance in the fairy shrimp diet with time. Our results reveal an ontogenetic change in the prey spectrum of fairy shrimp. The systematic shift towards omnivory likely affects both phyto- and zooplankton community composition, possibly contributing to temporal changes in food web dynamics in fairy shrimp habitats, and temporary ponds, which may warrant more detailed investigations in future studies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-022-05132-z.

The Structuring Effects of Salinity and Nutrient Status on Zooplankton Communities and Trophic Structure in Siberian Lakes

Many continental saline lakes are under the effects of salinity increase and anthropogenic eutrophication exacerbated by global change. The response of the food web to these drivers of change is not straightforward. To understand the consequences of salinity and eutrophication interactive effects on the food web, we studied the seasonal dynamics of zooplankton and phytoplankton and water quality parameters in 20 lakes of different salinity (from freshwater to hypersaline) and nutrient status (from oligotrophic to eutrophic) located in southern Siberia. We observed a pronounced bottom-up effect of nutrients, which induced an increase in the biomass of phytoplankton and zooplankton and a decline in water quality. A significant decrease in the species abundance of zooplankton was observed at a threshold salinity of 3 g L−1 and the disappearance of fish at 10 g L−1. The top-down effect induced by salinity manifested itself in an increase in the biomass of zooplankton with the disappearance of fish, and in the change of the size distribution of phytoplankton particles with an increase in the proportion of cladocerans in the zooplankton. Even though we observed that with the salinity increase the food web in saline lakes transformed from three-trophic to two-trophic without fish, we conclude that in the salinity range from 10 to 20–30 g L−1 this transition in most cases will not increase the ability of zooplankton to control phytoplankton. Interactive effects of salinity and eutrophication strongly depend on the size and depth of the lake, as deep stratified lakes tend to have a better water quality with lower biomasses of both phyto- and zooplankton. Thus, the salinity per se is not the driver of the decline in water clarity or the uncontrolled development of phytoplankton. Moreover, for deep lakes, salinity may be a factor affecting the stability of stratification, which mitigates the consequences of eutrophication. Thus, small shallow lakes will be the most vulnerable to the joint effect of salinity increase and eutrophication with the degradation of ecosystem functioning and water quality at moderate salinities of 3–20 g L−1.

Microbial communities of soda lakes and pans in the Carpathian Basin: a review

In this review, I would like to summarize the current knowledge on the microbiology of soda lakes and pans of the Carpathian Basin. First, the characteristic physical and chemical features of these sites are described. Most of the microbiological information presented deals with prokaryotes and algae, but protists and viruses are also mentioned. Planktonic bacterial communities are dominated by members of the phyla Actinobacteria, Bacteroidetes and Proteobacteria; small-sized trebouxiophycean green algae and Synechococcus/Cyanobium picocyanobacteria are the most important components of phytoplankton. Based on the current knowledge, it seems that mainly temperature, salinity, turbidity and grazing pressure regulate community composition and the abundance of individual microbial groups, but the external nutrient load from birds also has a significant impact on the ecological processes.

Key drivers structuring rotifer communities in ponds: insights into an agricultural landscape

Understanding the influence of environmental and spatial factors on the structure of aquatic communities remains a major challenge in community ecology. This study aims to identify main drivers of rotifer abundance and diversity in ponds embedded in an intensive agricultural landscape in Northeast Germany. We studied 42 ponds of glacial origin (kettle holes) covering a wide range of environmental parameters. The predominant factors structuring the rotifer metacommunity shifted from abiotic environmental filtering in spring to unstudied factors in autumn, while spatial factors remained less important. Fertilizer-driven salinization, internal nutrient recycling, primary productivity and sediment phosphorus release were the prevalent biogeochemical processes in the ponds. Both fertilizer-driven salinization and primary productivity negatively affected rotifer alpha diversity, and positively affected beta diversity. This impact was lower in forest ponds than in those surrounded by arable fields or grassland. However, rotifer diversity did not significantly differ among land-use categories. Our results indicate that the long-term impact of intensive agriculture in the region and the associated widespread eutrophication overrides the direct influence of land use on rotifer diversity but point to an indirect effect via fertilizer-driven salinization. Furthermore, this study highlights the role of ponds in enhancing regional biodiversity in agricultural landscapes.

Grazing pressure-induced shift in planktonic bacterial communities with the dominance of acIII-A1 actinobacterial lineage in soda pans

Astatic soda pans of the Pannonian Steppe are unique environments with respect to their multiple extreme physical and chemical characteristics (high daily water temperature fluctuation, high turbidity, alkaline pH, salinity, polyhumic organic carbon concentration, hypertrophic state and special ionic composition). However, little is known about the seasonal dynamics of the bacterial communities inhabiting these lakes and the role of environmental factors that have the main impact on their structure. Therefore, two soda pans were sampled monthly between April 2013 and July 2014 to reveal changes in the planktonic community. By late spring in both years, a sudden shift in the community structure was observed, the previous algae-associated bacterial communities had collapsed, resulting the highest ratio of Actinobacteria within the bacterioplankton (89%, with the dominance of acIII-A1 lineage) ever reported in the literature. Before these peaks, an extremely high abundance (> 10,000 individuum l-1) of microcrustaceans (Moina brachiata and Arctodiaptomus spinosus) was observed. OTU-based statistical approaches showed that in addition to algal blooms and water-level fluctuations, zooplankton densities had the strongest effect on the composition of bacterial communities. In these extreme environments, this implies a surprisingly strong, community-shaping top-down role of microcrustacean grazers.

Trait-based diatom functional diversity as an appropriate tool for understanding the effects of environmental changes in soda pans

Saline lakes, among the most seriously endangered ecosystems, are threatened due to climate change and human activities. One valuable feature of these environments is that they constitute areas of high biodiversity. Ecologists are, therefore, under great pressure to improve their understanding of the effects of natural and anthropogenic disturbances on the biodiversity of saline lakes. In this study, a total of 257 samples from 32 soda pans in Central Europe between 2006 and 2015 were examined. The effects of environmental variables and of geographical and limnoecological factors on functional diversity were analyzed. Furthermore, the explanatory power of the trait-based approach was assessed, and the applicability of the indices for biomonitoring purposes was determined. It was found that low habitat heterogeneity and harsh environments lead to the selection of a small number of suitable traits, and consequently, to a naturally low level of functional diversity. Anthropogenic activities enhance diversity at functional level due to the shift toward freshwater characteristics. On the regional scale, the effects of the region and status (natural, degraded, reconstructed) on diatom functional diversity were significant and more pronounced than that of the environmental and other limnoecological factors. The degree of variance found in functional diversity ascribed to environmental variables is five times greater in the case of the application of a trait-based approach, than when a taxonomic one is employed in the literature. Each of the tested functional diversity indices was sensitive to the most important environmental variables. Furthermore, these were type-specific and proved to be more complex indicators than taxonomic metrics. It is possible to suggest four functional diversity indices (FGR, FRic, FDis, and FDiv) which emphasize their independence from substrate and seasonal variations for ecological status assessment and conservation planning.

Trait-based community assembly of epiphytic diatoms in saline astatic ponds: a test of the stress-dominance hypothesis

The stress dominance hypothesis (SDH) postulates that strong environmental gradients drive trait convergence in communities over limiting similarity. Previous studies, conducted mostly with terrestrial plant communities, found controversial evidence for this prediction. We provide here the first test for SDH for epiphytic diatoms. We studied community assembly in diatom communities of astatic ponds. These water bodies serve as a good model system for testing SDH because they exhibit stress gradients of various environmental factors. Functional diversity of diatom communities was assessed based on four traits: (1) combined trait reflecting the trade-off between stress tolerance and competitive dominance, (2) cell size, (3) oxygen requirement and (4) N-uptake strategy. According to our results, salinity, pH and the width of the macrophyte belt appeared as significant predictors of the trait convergence/divergence patterns presumably acting through influencing the availability of carbon dioxide and turbidity. Lower trait diversity was found in turbid, more saline and more alkaline ponds and functional diversity was higher in transparent, less saline and less alkaline ponds. Overall, our results supported the stress dominance hypothesis. In habitats representing increased environmental stress, environmental filtering was the most important community assembly rule, while limiting similarity became dominant under more favourable conditions.

Dual bloom of green algae and purple bacteria in an extremely shallow soda pan

In April 2014, dual bloom of green algae and purple bacteria occurred in a shallow, alkaline soda pan (Kiskunság National Park, Hungary). The water was only 5 cm deep, in which an upper green layer was clearly separated from a near-sediment purple one. Based on microscopy and DNA-based identification, the upper was inhabited by a dense population of the planktonic green alga, Oocystis submarina Lagerheim, while the deeper layer was formed by purple, bacteriochlorophyll-containing bacteria, predominated by Thiorhodospira and Rhodobaca. Additional bacterial taxa with a presumed capability of anoxygenic phototrophic growth belonged to the genera Loktanella and Porphyrobacter. Comparing the bacterial community of the purple layer with a former blooming event in a nearby soda pan, similar functional but different taxonomic composition was revealed. Members from many dominant bacterial groups were successfully cultivated including potentially new species, which could be the result of the application of newly designed media.

Habitat loss over six decades accelerates regional and local biodiversity loss via changing landscape connectance

When habitats are lost, species are lost in the region as a result of the sampling process. However, it is less clear what happens to biodiversity in the habitats that remain. Some have argued that the main influence of habitat loss on biodiversity is simply due to the total amount of habitat being reduced, while others have argued that fragmentation leads to fewer species per site because of altered spatial connectance among extant habitats. Here, we use a unique data set on invertebrate species in ponds spanning six decades of habitat loss to show that both regional and local species richness declined, indicating that species loss is compounded by habitat loss via connectivity loss, and not a result of a sampling process or changes in local environmental conditions. Overall, our work provides some of the clearest evidence to date from a longitudinal study that habitat loss translates into species loss, even within the remaining habitats.

Unexpected consequences of bombing. Community level response of epiphytic diatoms to environmental stress in a saline bomb crater pond area

The spatial response of epiphytic diatom communities to environmental stress was studied in a moderately saline wetland area located in the plain of Danube-Tisza Interfluve, Hungary. The area is characterised by World War II bomb crater ponds and can be regarded as an excellent ecological model system where the dispersion of species is slightly limited by distance. To study the effect of environmental variables on the communities, canonical correspondence analysis was applied. Salinity, pH, total suspended solids, total phosphorous and depth proved to be significant environmental drivers in this analysis. The ecological status of the ponds was assessed with Ziemann’s halobity index, as the trophity-depending metric cannot be applied to these habitats (due to the naturally high phosphorus content). Ponds in “good” ecological status significantly differed from those appertaining to water quality category of “not-good” ecological status considering characteristic of natural astatic soda pans (e.g. salinity, pH, ammonium, total phosphorous concentration, nitrogen:phosphorous ratio and turbidity). The differences between epiphytic diatom communities inhabiting the ponds were detected using non-parametric multidimensional scaling. The samples formed three groups according to the types of ponds (“transparent”, “transitional” and “turbid”) based on the width of the macrophyte belt around them. Indicator species related to the ecological status of the ponds and diatom communities contributing to the separation of groups of ponds were identified. One of the indicator species differed from species already described. Light and scanning electron microscopy features and phylogenetic analyses based on three genes (18S and 28S rRNA genes, rbcL) proved that it was a new species of Nitzschia genus, closely related to Nitzschia frustulum and Nitzschia inconspicua. Therefore, description of a new species, Nitzschia reskoi Ács, Duleba, C.E.Wetzel & Ector is proposed. We concluded that the increasing abundance of Nitzschia reskoi was a signal of the degradation of the intermittent saline wetlands.

Trophic state (TSISD) and mixing type significantly influence pelagic zooplankton biodiversity in temperate lakes (NW Poland)

Background

Lake depth and the consequent mixing regime and thermal structure have profound effects on ecosystem functioning, because depth strongly affects the availability of nutrients, light, and oxygen. All these conditions influence patterns of zooplankton diversity. Zooplankton are a key component of the aquatic environment and are essential to maintaining natural processes in freshwater ecosystems. However, zooplankton biodiversity can be different regard to depth, mixing type and trophic state. Therefore, the aim of this study was to examine how depth and mixing regime affect zooplankton diversity in lakes. We also investigated the vertical distribution of diversity across a trophic gradient of lakes.

Methods

A total of 329 zooplankton samples from 79 temperate lakes (36 polymictic and 43 dimictic) were collected. The biodiversity of zooplankton was calculated using species richness (SR) and the Shannon index (SI). An index based on Secchi disc visibility was used to determine the trophic state index (TSI_SD) of lakes. The one-way ANOVA with Duncan's post hoc test were used to determine differences in zooplankton biodiversity between mictic lake types and thermal layers. To find the best predictors for zooplankton biodiversity a multiple stepwise regression was used. The rarefaction method was used to evaluate the impact of mixing types, thermal layers, and the TSI_SDon zooplankton biodiversity indices. A Sørensen similarity analysis and nonmetric multidimensional scaling (NMDS) were performed to describe the similarity patterns in species composition among lakes.

Results

We identified a total of 151 taxa from 36 polymictic and 43 dimictic lakes. Lake depth and the TSI_SD were significantly correlated with the biodiversity of lake zooplankton. The results of ANOVA and Duncan tests show that mictic type and thermal zones had a significant effect on zooplankton biodiversity. The rarefaction curve showed significant differences in zooplankton biodiversity, which was greater in lakes with lower trophic state. Ordination by NMDS showed clustering of different mictic types, thermal layers, and composition changes throughout the TSI_SDprofile. Moreover, we determined that polymictic lakes are more heterogeneous than dimictic lakes in regard to zooplankton similarities.

Discussion

Both mictic lake types were characterized by varying levels of zooplankton biodiversity, which is shaped by the communities' response to lake depth, thermal layers and TSI_SD values. The zooplankton SR and SI (during daylight hours) depends greatly on the mixing type. Lake type also indicates the importance of the metalimnion in shaping zooplankton biodiversity in dimictic lakes. In addition, data from NW Polish lakes indicated that the increase of the TSI_SD leads to taxonomic shifts and has a negative effect on the diversity of all groups of zooplankton.

Wartime scars or reservoirs of biodiversity? The value of bomb crater ponds in aquatic conservation

Considering the ongoing loss of aquatic habitats, anthropogenic ponds are gaining importance as substitute habitats. It is therefore important to assess their functioning in comparison to their natural precursors. Here we assess the biodiversity value of sodic bomb crater ponds by comparing their gamma diversity to that of natural reference habitats, astatic soda pans, and assess their importance on the landscape level by studying alpha and beta diversity. We studied aquatic organisms ranging from algae to vertebrates in a dense cluster of 54 sodic bomb crater ponds in Central Europe. Despite the overall small area of the pond cluster, gamma diversity was comparable to that found in surveys of natural habitats that encompassed much wider spatial and temporal scales. We also found a considerable number of species shared with reference habitats, indicating that these anthropogenic habitats function as important refuge sites for several species that are associated with the endangered soda pans. Moreover, we found a number of regionally or worldwide rare species. Among the components of beta diversity, species replacement dominated community assembly. Individual ponds contributed similarly to beta diversity in terms of replacement, being equally important for maintaining high gamma diversity and emphasising the role of the pond network rather than individual ponds. This pattern was seen in all studied groups. Bomb crater ponds therefore acted as important contributors to aquatic biodiversity. Considering the tremendous losses of ponds throughout Europe, anthropogenic ponds should be taken into consideration in nature conservation, especially when occurring in pond networks.

Multiple extreme environmental conditions of intermittent soda pans in the Carpathian Basin (Central Europe)

Soda lakes and pans represent saline ecosystems with unique chemical composition, occurring on all continents. The purpose of this study was to identify and characterise the main environmental gradients and trophic state that prevail in the soda pans (n=84) of the Carpathian Basin in Central Europe. Underwater light conditions, dissolved organic matter, phosphorus and chlorophyll a were investigated in 84 pans during 2009-2010. Besides, water temperature was measured hourly with an automatic sensor throughout one year in a selected pan. The pans were very shallow (median depth: 15 cm), and their extremely high turbidity (Secchi depth median: 3 cm, min: 0.5 cm) was caused by high concentrations of inorganic suspended solids (median: 0.4 g L^-1, max: 16 g L^-1), which was the dominant (>50%) contributing factor to the vertical attenuation coefficient in 67 pans (80%). All pans were polyhumic (median DOC: 47 mg L^-1), and total phosphorus concentration was also extremely high (median: 2 mg L^-1, max: 32 mg L^-1). The daily water temperature maximum (44 °C) and fluctuation maximum (28 °C) were extremely high during summertime. The combination of environmental boundaries: shallowness, daily water temperature fluctuation, intermittent hydroperiod, high turbidity, polyhumic organic carbon concentration, high alkalinity and hypertrophy represent a unique extreme aquatic ecosystem.

Wind dispersal results in a gradient of dispersal limitation and environmental match among discrete aquatic habitats

Directional dispersal by wind and other dispersal agents may generate spatial patterns in passively dispersing metacommunities which cannot be detected by classical eigenvector methods based on Euclidean distances. We analysed zooplankton communities (Rotifera, Cladocera, Copepoda) in a cluster of soda pans distributed over a short spatial scale of 18 km and tested explicitly for directional signals in their spatial configuration. The study area is exposed to a prevailing northwestern wind direction. By applying asymmetric eigenvector maps (AEM), we were able to identify corresponding directionality in the spatial structure of communities. Furthermore, the match between community composition and environmental conditions exhibited a spatial pattern consistent with the prevailing wind corridor, with best match found downwind the dominant wind direction. We also found that classical eigenvector methods based on Euclidean distances underestimated the role of spatial processes in our data. Our study furthermore shows that dispersal limitation may constrain community assembly in highly mobile organisms even at spatial scales below 5 km.