Eutrophication increases the similarity of cyanobacterial community features in lakes and reservoirs.

Eutrophication of inland waters is a mostly anthropogenic phenomenon impacting aquatic biodiversity worldwide, and might change biotic community structure and ecosystem functions. However, little is known about the patterns of cyanobacterial community variations and changes both on alpha and beta diversity levels in response to eutrophication. Here, we investigated cyanobacterial communities sampled at 140 sites from 59 lakes and reservoirs along a strong eutrophication gradient in eastern China through using CPC-IGS and 16S rRNA gene amplicon sequencing. We found that taxonomic diversity increased, but phylogenetic diversity decreased significantly along the eutrophication gradient. Both niche width and niche overlap of cyanobacteria significantly decreased from low- to high-nutrient waterbodies. Cyanobacterial community distance-decay relationship became weaker from mesotrophic to hypereutrophic waterbodies, while ecological uniqueness (i.e., local contributions to beta diversity) tended to increase in high-nutrient waterbodies. Latitude and longitude were more important in shaping cyanobacterial community structure than other environmental variables. These findings suggest that eutrophication affects alpha and beta diversity of cyanobacterial communities, leading to increasingly similar community structures in lakes and reservoirs with a higher level of eutrophication. Our work highlights how cyanobacterial communities respond to anthropogenic eutrophication and calls for an urgent need to develop conservation and management strategies to control lake eutrophication and protect freshwater biodiversity.

Local eukaryotic and bacterial stream community assembly is shaped by regional land use effects.

With anticipated expansion of agricultural areas for food production and increasing intensity of pressures stemming from land-use, it is critical to better understand how species respond to land-use change. This is particularly true for microbial communities which provide key ecosystem functions and display fastest responses to environmental change. However, regional land-use effects on local environmental conditions are often neglected, and, hence, underestimated when investigating community responses. Here we show that the effects stemming from agricultural and forested land use are strongest reflected in water conductivity, pH and phosphorus concentration, shaping microbial communities and their assembly processes. Using a joint species distribution modelling framework with community data based on metabarcoding, we quantify the contribution of land-use types in determining local environmental variables and uncover the impact of both, land-use, and local environment, on microbial stream communities. We found that community assembly is closely linked to land-use type but that the local environment strongly mediates the effects of land-use, resulting in systematic variation of taxon responses to environmental conditions, depending on their domain (bacteria vs. eukaryote) and trophic mode (autotrophy vs. heterotrophy). Given that regional land-use type strongly shapes local environments, it is paramount to consider its key role in shaping local stream communities.

European river typologies fail to capture diatom, fish, and macrophyte community composition.

Typology systems are frequently used in applied and fundamental ecology and are relevant for environmental monitoring and conservation. They aggregate ecosystems into discrete types based on biotic and abiotic variables, assuming that ecosystems of the same type are more alike than ecosystems of different types with regard to a specific property of interest. We evaluated whether this assumption is met by the Broad River Types (BRT), a recently proposed European river typology system, that classifies river segments based on abiotic variables, when it is used to group biological communities. We compiled data on the community composition of diatoms, fishes, and aquatic macrophytes throughout Europe and evaluated whether the composition is more similar in site groups with the same river type than in site groups of different river types using analysis of similarities, classification strength, typical species analysis, and the area under zeta diversity decline curves. We compared the performance of the BRT with those of four region-based typology systems, namely, Illies Freshwater Ecoregions, the Biogeographic Regions, the Freshwater Ecoregions of the World, and the Environmental Zones, as well as spatial autocorrelation (SA) classifications. All typology systems received low scores from most evaluation methods, relative to predefined thresholds and the SA classifications. The BRT often scored lowest of all typology systems. Within each typology system, community composition overlapped considerably between site groups defined by the types of the systems. The overlap tended to be the lowest for fishes and between Illies Freshwater Ecoregions. In conclusion, we found that existing broad-scale river typology systems fail to delineate site groups with distinct and compositionally homogeneous communities of diatoms, fishes, and macrophytes. A way to improve the fit between typology systems and biological communities might be to combine segment-based and region-based typology systems to simultaneously account for local environmental variation and historical distribution patterns, thus potentially improving the utility of broad-scale typology systems for freshwater biota.

Seasonal variation of phytoplankton community assembly processes in Tibetan Plateau floodplain.

Uncovering the mechanisms underlying phytoplankton community assembly remains a major challenge in freshwater ecology. The roles of environmental filtering and spatial processes in shaping phytoplankton metacommunity in Tibetan floodplain ecosystems under various hydrological conditions are still unclear. Here, multivariate statistics and a null model approach were used to compare the spatiotemporal patterns and assembly processes of phytoplankton communities in the river-oxbow lake system of Tibetan Plateau floodplain between non-flood and flood periods. The results showed that phytoplankton communities had significant seasonal and habitat variations, with the seasonal variations being more remarkable. Phytoplankton density, biomass, and alpha diversity were distinctly lower in the flood than non-flood period. The habitat differences (rivers vs. oxbow lakes) in phytoplankton community were less pronounced during the flood than non-flood period, most likely due to the increased hydrological connectivity. There was a significant distance-decay relationship only in lotic phytoplankton communities, and such relationship was stronger in the non-flood than flood period. Variation partitioning and PER-SIMPER analysis showed that the relative role of environmental filtering and spatial processes affecting phytoplankton assemblages varied across hydrological periods, with environmental filtering dominating in the non-flood period and spatial processes in the flood period. These results suggest that the flow regime plays a key role in balancing environmental and spatial factors in shaping phytoplankton communities. This study contributes to a deeper understanding of ecological phenomena in highland floodplains and provides a theoretical basis for floodplain ecosystem maintenance and ecological health management.

On the shape and origins of the freshwater species-area relationship.

The species-area relationship (SAR) has over a 150-year-long history in ecology, but how its shape and origins vary across scales and organisms remains incompletely understood. This is the first subcontinental freshwater study to examine both these properties of the SAR in a spatially explicit way across major organismal groups (diatoms, insects, and fish) that differ in body size and dispersal capacity. First, to describe the SAR shape, we evaluated the fit of three commonly used models, logarithmic, power, and Michaelis-Menten. Second, we proposed a hierarchical framework to explain the variability in the SAR shape, captured by the parameters of the SAR model. According to this framework, scale and species group were the top predictors of the SAR shape, climatic factors (heterogeneity and median conditions) represented the second predictor level, and metacommunity properties (intraspecific spatial aggregation, γ-diversity, and species abundance distribution) the third predictor level. We calculated the SAR as a sample-based rarefaction curve using 60 streams within landscape windows (scales) in the United States, ranging from 160,000 to 6,760,000 km2 . First, we found that all models provided good fits (R2 ≥ 0.93), but the frequency of the best-fitting model was strongly dependent on organism, scale, and metacommunity properties. The Michaelis-Menten model was most common in fish, at the largest scales, and at the highest levels of intraspecific spatial aggregation. The power model was most frequent in diatoms and insects, at smaller scales, and in metacommunities with the lowest evenness. The logarithmic model fit best exclusively at the smallest scales and in species-poor metacommunities, primarily fish. Second, we tested our framework with the parameters of the most broadly used SAR model, the log-log form of the power model, using a structural equation model. This model supported our framework and revealed that the SAR slope was best predicted by scale- and organism-dependent metacommunity properties, particularly spatial aggregation, whereas the intercept responded most strongly to species group and γ-diversity. Future research should investigate from the perspective of our framework how shifts in metacommunity properties due to climate change may alter the SAR.

How does invasion degree shape alpha and beta diversity of freshwater fish at a regional scale?

Freshwater ecosystems appear more vulnerable to biodiversity loss due to several anthropogenic disturbances and freshwater fish are particularly vulnerable to these impacts. We aimed to (1) identify the contribution of land use, spatial variables, and invasion degree in determining freshwater fish alpha (i.e., species richness) and beta (i.e., local contributions to beta diversity, LCBD) diversity, evaluating also the relationship between invasion degree and nestedness ( ß nes) and turnover ( ß sim) components of beta diversity. (2) Investigate the relationship between alpha diversity and LCBD, under the hypothesis that alpha diversity and LCBD correlate negatively and (3) investigate the relationship between species contributions to beta diversity (SCBD) and species occurrence, hypothesizing that non-native species show a lower contribution to beta diversity. The linear mixed models and the partition of R 2 retained the invasion degree as the most important variables explaining alpha and beta diversity, having a positive relationship with both diversity components. Furthermore, land use related to human impacts had a positive influence on alpha diversity, whereas it showed a negative effect on LCBD. Regression model further showed that invasion degree related positively with ß sim, but negatively with ß nes, suggesting that non-native species were involved in the replacement of native species in the fish community. Alpha diversity and LCBD showed a weak positive correlation, meaning that sites with low species richness have higher LCBD. SCBD scaled positively with species occurrence highlighting that rarer species contribute less to SCBD. Finally, native and exotic species contributed similarly to beta diversity. These results suggest that invasion degree plays a central role in shaping alpha and beta diversity in stream fish, more than land use features reflecting habitat alteration or other geospatial variables. Furthermore, it is important to evaluate separately the native and the non-native components of biotic communities to identify linkages between invasion dynamics and biodiversity loss.

The scale-dependence of spatial distribution of reservoir plankton communities in subtropical and tropical China.

Distance-decay relationships (DDRs) represent a very useful approach to describing the spatial distribution of biological communities. However, plankton DDR patterns and community assembly mechanisms are still poorly understood at different spatial scales in reservoir ecosystems. We collected phytoplankton, zooplankton and water samples in 24 reservoirs from subtropical and tropical China from July to August 2018. We examined DDR patterns across three distinct spatial scales, i.e., within-reservoir, within-drainage (but between reservoirs) and between drainages. We tested whether the rate of change (i.e., slope) of DDRs is consistent across different spatial scales. We assessed the relative importance of spatial and environmental variables in shaping the community distribution of plankton and quantitatively distinguished the community assembly mechanisms. We observed significant DDR curves in phytoplankton and zooplankton communities, in which slopes of the DDRs were steepest at the smallest spatial scale. Both spatial and environmental factors had significant impacts on DDR and dispersal assembly was a slightly stronger process in reservoir phytoplankton and zooplankton community assembly than niche-based process. We conclude that DDRs of reservoir phytoplankton and zooplankton vary with spatial scale. Our data shed light on how spatial and environmental variables contribute to plankton community assembly together. However, we revealed that dispersal process contributes to the biogeography of reservoir plankton slightly more strongly than environmental filtering. Collectively, this study enhances the understanding of plankton biogeography and distribution at multiple spatial scales.

Microbial and Environmental Processes Shape the Link between Organic Matter Functional Traits and Composition.

Dissolved organic matter (DOM) is a large and complex mixture of molecules that fuels microbial metabolism and regulates biogeochemical cycles. Individual DOM molecules have unique functional traits, but how their assemblages vary deterministically under global change remains poorly understood. Here, we examine DOM and associated bacteria in 300 aquatic microcosms deployed on mountainsides that span contrasting temperatures and nutrient gradients. Based on molecular trait dimensions of reactivity and activity, we partition the DOM composition into labile-active, recalcitrant-active, recalcitrant-inactive, and labile-inactive fractions and quantify the relative influences of deterministic and stochastic processes governing the assembly of each. At both subtropical and subarctic study sites, the assembly of labile or recalcitrant molecules in active fractions is primarily governed by deterministic processes, while stochastic processes are more important for the assembly of molecules within inactive fractions. Surprisingly, the importance of deterministic selection increases with global change gradients for recalcitrant molecules in both active and inactive fractions, and this trend is paralleled by changes in the deterministic assembly of microbial communities and environmental filtering, respectively. Together, our results highlight the shift in focus from potential reactivity to realized activity and indicate that active and inactive fractions of DOM assemblages are structured by contrasting processes, and their recalcitrant components are consistently sensitive to global change. Our study partitions the DOM molecular composition across functional traits and links DOM with microbes via a shared ecological framework of assembly processes. This integrated approach opens new avenues to understand the assembly and turnover of organic carbon in a changing world.

Ecological networks of dissolved organic matter and microorganisms under global change.

Microbes regulate the composition and turnover of organic matter. Here we developed a framework called Energy-Diversity-Trait integrative Analysis to quantify how dissolved organic matter and microbes interact along global change drivers of temperature and nutrient enrichment. Negative and positive interactions suggest decomposition and production processes of organic matter, respectively. We applied this framework to manipulative field experiments on mountainsides in subarctic and subtropical climates. In both climates, negative interactions of bipartite networks were more specialized than positive interactions, showing fewer interactions between chemical molecules and bacterial taxa. Nutrient enrichment promoted specialization of positive interactions, but decreased specialization of negative interactions, indicating that organic matter was more vulnerable to decomposition by a greater range of bacteria, particularly at warmer temperatures in the subtropical climate. These two global change drivers influenced specialization of negative interactions most strongly via molecular traits, while molecular traits and bacterial diversity similarly affected specialization of positive interactions.

Studying biodiversity-ecosystem function relationships in experimental microcosms among islands.

Ecological studies on islands have provided fundamental insights into the mechanisms underlying biodiversity of larger organisms, but we know little about the factors affecting island microbial biodiversity and ecosystem function. We conducted a field experiment on five Baltic Sea islands where we placed aquatic microcosms with different levels of salinity mimicking environmental stress and allowed diatoms to colonize the microcosms via the air. Using structural equation models (SEM), we investigated the interconnections among environmental and dispersal-related factors, diatom biodiversity, and ecosystem productivity (represented by chlorophyll a concentration). We also tested whether the body size structure of the community influences productivity together with biodiversity. In SEMs, we found no relationship between species richness or evenness and productivity. However, productivity increased with increasing mean body size of species in the communities. The effects of environmental stress on both biodiversity and ecosystem productivity were highlighted as species richness and evenness declined, whereas productivity increased at the highest salinity levels. In addition to salinity, wind exposure affected both biodiversity metrics and productivity. This study provides new insights into microbial community assembly in a field experimental setting and the relationship between biodiversity and ecosystem function. Our results indicate that salinity presents a strong abiotic filter, leading to communities that may be species poor, yet comprise salinity-tolerant and relatively productive species at high salinity. Our findings also emphasize the importance of mean community body size in mediating the effects of environmental conditions on productivity and suggest that this trait should be considered more broadly in biodiversity-ecosystem function studies.

Embracing mountain microbiome and ecosystem functions under global change.

Mountains are pivotal to maintaining habitat heterogeneity, global biodiversity, ecosystem functions and services to humans. They have provided classic model natural systems for plant and animal diversity gradient studies for over 250 years. In the recent decade, the exploration of microorganisms on mountainsides has also achieved substantial progress. Here, we review the literature on microbial diversity across taxonomic groups and ecosystem types on global mountains. Microbial community shows climatic zonation with orderly successions along elevational gradients, which are largely consistent with traditional climatic hypotheses. However, elevational patterns are complicated for species richness without general rules in terrestrial and aquatic environments and are driven mainly by deterministic processes caused by abiotic and biotic factors. We see a major shift from documenting patterns of biodiversity towards identifying the mechanisms that shape microbial biogeographical patterns and how these patterns vary under global change by the inclusion of novel ecological theories, frameworks and approaches. We thus propose key questions and cutting-edge perspectives to advance future research in mountain microbial biogeography by focusing on biodiversity hypotheses, incorporating meta-ecosystem framework and novel key drivers, adapting recently developed approaches in trait-based ecology and manipulative field experiments, disentangling biodiversity-ecosystem functioning relationships and finally modelling and predicting their global change responses.

Pathways for cross-boundary effects of biodiversity on ecosystem functioning.

The biodiversity-ecosystem functioning concept asserts that processes in ecosystems are markedly influenced by species richness and other facets of biodiversity. However, biodiversity-ecosystem functioning studies have been largely restricted to single ecosystems, ignoring the importance of functional links - such as the exchange of matter, energy, and organisms - between coupled ecosystems. Here we present a basic concept and outline three pathways of cross-boundary biodiversity effects on ecosystem processes and propose an agenda to assess such effects, focusing on terrestrial-aquatic linkages to illustrate the case. This cross-boundary perspective of biodiversity-ecosystem functioning relationships presents a promising frontier for biodiversity and ecosystem science with repercussions for the conservation, restoration, and management of biodiversity and ecosystems from local to landscape scales.

The application of Uniform Manifold Approximation and Projection (UMAP) for unconstrained ordination and classification of biological indicators in aquatic ecology.

The analysis of community structure in studies of freshwater ecology often requires the application of dimensionality reduction to process multivariate data. A high number of dimensions (number of taxa/environmental parameters × number of samples), nonlinear relationships, outliers, and high variability usually hinder the visualization and interpretation of multivariate datasets. Here, we proposed a new statistical design using Uniform Manifold Approximation and Projection (UMAP), and community partitioning using Louvain algorithms, to ordinate and classify the structure of aquatic biota in two-dimensional space. We present this approach with a demonstration of five previously published datasets for diatoms, macrophytes, chironomids (larval and subfossil), and fish. Principal Component Analysis (PCA) and Ward's clustering were also used to assess the comparability of the UMAP approach compared to traditional approaches for ordination and classification. The ordination of sampling sites in 2-dimensional space showed a much denser, and easier to interpret, grouping using the UMAP approach in comparison to PCA. The classification of community structure using the Louvain algorithm in UMAP ordinal space showed a high classification strength for data with a high number of dimensions than the cluster patterns obtained with the use of a Ward's algorithm in PCA. Environmental gradients, presented via heat maps, were overlayed with the ordination patterns of aquatic communities, confirming that the ordinations obtained by UMAP were ecologically meaningful. This is the first study that has applied a UMAP approach with classification using Louvain algorithms on ecological datasets. We show that the performance of local and global structures, as well as the number of clusters determined by the algorithm, make this approach more powerful than traditional approaches.

Anthropogenic land‒use impacts on the size structure of macroinvertebrate assemblages are jointly modulated by local conditions and spatial processes.

Body size descriptors and associated resemblance measurements may provide useful tools for forecasting ecological responses to increasing anthropogenic land‒use disturbances. Yet, the influences of agriculture and urbanisation on the size structure of biotic assemblages have seldom been investigated in running waters. Using a comprehensive dataset on stream macroinvertebrates from 21 river basins across Western Finland, we assessed whether the structure of assemblages via changes in taxonomic composition and body size distributions responded predictably to anthropogenic land‒use impacts. Specifically, we applied a combination of resemblance measurements based on cumulative abundance profiles and spatially constrained null models to understand faunal impairment by agricultural and urban development, and the most likely mechanisms underlying the observed shifts in assemblage size structure. Anthropogenically impacted stream sites showed less variation in assemblage composition and size distributions compared with least‒disturbed sites, with strong declines in internal variation also occurring for the transition from near‒pristine to moderately impacted landscapes. These results were consistent whether based on species‒level or genus‒level data. Variation in assemblage size structure seemed to be more predictable than taxonomic composition, supporting the notion that resemblance measurements based on body size distributions can represent an improvement to more traditional approaches based on taxonomic identities alone. In addition, we showed that macroinvertebrate assemblages resulted from effects of land‒use degradation mediated through local conditions and spurious spatial structures in the distribution of anthropogenic activities across the landscape. Overall, our findings suggest that existing water policies and agri‒environment schemes should be guided not only by understanding the individual effects of agricultural and urban development on taxonomic composition at a given stream site. Rather, we should also acknowledge the size structure of stream assemblages and whether concomitant changes in local conditions and the non‒random distribution of human infrastructures are likely to mitigate or accelerate these effects.

Taxonomic and functional diversity covary in rock pool microalgal communities despite their different drivers.

Local biodiversity has traditionally been estimated with taxonomic diversity metrics such as species richness. Recently, the concept of biodiversity has been extended beyond species identity by ecological traits determining the functional role of a species in a community. This interspecific functional diversity typically responds more strongly to local environmental variation compared with taxonomic diversity, while taxonomic diversity may mirror more strongly dispersal processes compared with functional metrics. Several trait-based indices have been developed to measure functional diversity for various organisms and habitat types, but studies of their applicability on aquatic microbial communities have been underrepresented. We examined the drivers and covariance of taxonomic and functional diversity among diatom rock pool communities on the Baltic Sea coast. We quantified three taxonomic (species richness, Shannon's diversity, and Pielou's evenness) and three functional (functional richness, evenness, and divergence) diversity indices and determined abiotic factors best explaining variation in these indices by generalized linear mixed models. The six diversity indices were highly collinear except functional evenness, which merely correlated significantly with taxonomic evenness. All diversity indices were always explained by water conductivity and temperature-sampling month interaction. Taxonomic diversity was further consistently explained by pool distance to the sea, and functional richness and divergence by pool location. The explained variance in regression models did not markedly differ between taxonomic and functional metrics. Our findings do not clearly support the superiority of neither set of diversity indices in explaining coastal microbial diversity, but rather highlight the general overlap among the indices. However, as individual metrics may be driven by different factors, the greatest advantage in assessing biodiversity is nevertheless probably achieved with a simultaneous application of the taxonomic and functional diversity metrics.

Disentangling the relative roles of natural and anthropogenic-induced stressors in shaping benthic ciliate diversity in a heavily disturbed bay.

Coastal areas are facing biodiversity loss and degradation of habitats due to intensified human activities. However, our understanding of the relative contribution of natural gradients and human induced disturbance to biodiversity is limited. Here, we investigated the response of three facets of alpha and beta diversity of benthic ciliates to environmental gradients in a highly disturbed estuarine bay in China. We used linear regression and distance-based redundancy analysis to determine the key driving factors for biodiversity. Variation partitioning was further used to examine the relative influence of natural gradients and anthropogenic disturbances on ciliate communities. Our results revealed that ciliate alpha diversity and functional composition remained similar despite notable variation in species composition along salinity gradient. Sediment grain size, together with heavy metals were the strongest determinants in shaping both alpha and beta diversity. After controlling for the effect of natural factors, heavy metals still had significant impacts on beta diversity. Human induced nitrogen enrichment was positively correlated with algivorous functional group with possible impacts on benthic food webs. These results suggest that beta diversity is overall more sensitive to anthropogenic stressors than alpha diversity and give important insights into the role of anthropogenic disturbance on coastal diversity, being also useful for developing ecosystem protection and conservation strategies.

Are Bacterio- and Phytoplankton Community Compositions Related in Lakes Differing in Their Cyanobacteria Contribution and Physico-Chemical Properties?

Bacterioplankton community composition has become the center of research attention in recent years. Bacteria associated with toxic cyanobacteria blooms have attracted considerable interest. However, little is known about the environmental factors driving the bacteria community, including the impact of invasive cyanobacteria. Therefore, our aim has been to determine the relationships between heterotrophic bacteria and phytoplankton community composition across 24 Polish lakes with different contributions of cyanobacteria including the invasive species Raphidiopsis raciborskii. This analysis revealed that cyanobacteria were present in 16 lakes, while R. raciborskii occurred in 14 lakes. Our results show that bacteria communities differed between lakes dominated by cyanobacteria and lakes with minor contributions of cyanobacteria but did not differ between lakes with R. raciborskii and other lakes. Physical factors, including water and Secchi depth, were the major drivers of bacteria and phytoplankton community composition. However, in lakes dominated by cyanobacteria, bacterial community composition was also influenced by biotic factors such as the amount of R. raciborskii, chlorophyll-a and total phytoplankton biomass. Thus, our study provides novel evidence on the influence of environmental factors and R. raciborskii on lake bacteria communities.

Does trait-based joint species distribution modelling reveal the signature of competition in stream macroinvertebrate communities?

The occupancy and abundance of species are jointly driven by local factors, such as environmental characteristics and biotic interactions, and regional-scale factors, such as dispersal and climate. Recently, it has been shown that biotic interactions shape species occupancies and abundances beyond local extents. However, for small ectothermic animals, particularly for those occurring in freshwater environments, the importance of biotic interactions remains understudied. Species-to-species associations from joint species distribution models (i.e. species associations while controlling for environmental characteristics) are increasingly used to draw hypotheses of which species possibly show biotic interactions. We studied whether species-to-species associations from joint species distribution models show signs of competition using a hypothesis testing framework in stream macroinvertebrate communities at regional extent. We sampled aquatic macroinvertebrates from 105 stream sites in western Finland encompassing a latitudinal gradient of c. 500 kilometres. We hypothesized that if competition drives these associations (H1) functionally, similar species are mostly negatively associated, whereas functionally dissimilar species show random associations. We further hypothesized that the relationship between functional dissimilarity and the strength of association is more pronounced (H2) for abundances rather than occupancies, (H3) at small grain (i.e. stream site) rather than at large grain (i.e. river basin), and (H4) among species having weak dispersal ability than among species with high dispersal ability. Stream macroinvertebrates showed both negative and positive species-to-species associations while controlling for habitat characteristics. However, the negative associations were mostly at large grain (river basin) rather than at small grain (stream site), in occupancy rather than abundance, and not related to species functional dissimilarity or to their dispersal ability. Thus, all our hypotheses considering possible competition (H1-H4) were rejected. Competition does not appear to be a major driving force of stream macroinvertebrate communities at the spatial grain sizes considered. The observed positive associations in occupancy at small grain (stream site) may be attributed to species' similar microhabitat preferences, whereas at large grain (river basin), they may stem from metacommunity dynamics. Our results highlight that species traits were necessary to interpret whether or not species-to-species associations from joint species distribution models resulted from biotic interactions.

Cross-taxon congruence of aquatic microbial communities across geological ages in Iceland: Stochastic and deterministic processes.

Biotic groups usually have nonrandom cross-taxon relationships in their biodiversity or compositions across sites, but it is poorly known how such congruence varies across long-term ecosystem development, and what are the ecological processes underlying biodiversity patterns. Here, we examined the cross-taxon congruence in diversity and compositions of bacteria, fungi and diatoms in streams from four regions with different geological ages in Iceland, and further studied their community assembly processes. Bacteria and fungi showed contrasting trends in alpha and gamma diversities across geological ages, while their beta diversity patterns were consistent, being the lowest in the oldest region. The three taxonomic groups had the strongest cross-taxon congruence of beta diversity in the oldest region, while the weakest for intermediate-aged regions. Although environmental variables played important roles in shaping their congruence, biotic interaction had nonnegligible influences. Deterministic processes, being dominant for bacteria and fungi, had the highest relative influence in intermediate-aged regions, whereas diatoms showed higher stochasticity. We proposed a four-phase conceptual model to show how the balance of deterministic and stochastic processes changes across geological ages. Taken together, our results provide an advanced understanding of cross-taxon congruence and community assembly processes for aquatic communities over long-term periods of geological age.

Biodiversity Loss Threatens the Current Functional Similarity of Beta Diversity in Benthic Diatom Communities.

The global biodiversity loss has increased the need to understand the effects of decreasing diversity, but our knowledge on how species loss will affect the functioning of communities and ecosystems is still very limited. Here, the levels of taxonomic and functional beta diversity and the effect of species loss on functional beta diversity were investigated in an estuary that provides a naturally steep environmental gradient. The study was conducted using diatoms that are among the most important microorganisms in all aquatic ecosystems and globally account for 40% of marine primary production. Along the estuary, the taxonomic beta diversity of diatom communities was high (Bray-Curtis taxonomic similarity 0.044) and strongly controlled by the environment, particularly wind exposure, salinity, and temperature. In contrast, the functional beta diversity was low (Bray-Curtis functional similarity 0.658) and much less controlled by the environment. Thus, the diatom communities stayed functionally almost similar despite large changes in species composition and environment. This may indicate that, through high taxonomic diversity and redundancy in functions, microorganisms provide an insurance effect against environmental change. However, when studying the effect of decreasing species richness on functional similarity of communities, simulated species loss to 45% of the current species richness decreased functional similarity significantly. This suggests that decreasing species richness may increase variability and reduce the stability and resilience of communities. These results highlight the importance of high taxonomic biodiversity for the stable functioning of benthic communities.