Functional rarity and evenness are key facets of biodiversity to boost multifunctionality

Dynamic shifts of functional diversity through climate-resilient strategies and farmland restoration in a mountain protected area

Land-use land-cover (LULC) change contributes to major ecological impacts, particularly in areas undergoing land abandonment, inducing modifications on habitat structure and species distributions. Alternative land-use policies are potential solutions to alleviate the negative impacts of contemporary tendencies of LULC change on biodiversity. This work analyzes these tendencies in the Montesinho Natural Park (Portugal), an area representative of European abandoned mountain rural areas. We built ecological niche models for 226 species of vertebrates (amphibians, reptiles, birds, and mammals) and vascular plants, using a consensus modelling approach available in the R package 'biomod2'. We projected the models to contemporary (2018) and future (2050) LULC scenarios, under four scenarios aiming to secure relevant ecosystem services and biodiversity conservation for 2050: an afforestation and a rewilding scenario, focused on climate-smart management strategies, and a farmland and an agroforestry recovery scenario, based on re-establishing human traditional activities. We quantified the influences of these scenarios on biodiversity through species habitat suitability changes for 2018-2050. We analyzed how these management strategies could influence indices of functional diversity (functional richness, functional evenness and functional dispersion) within the park. Habitat suitability changes revealed complementary patterns among scenarios. Afforestation and rewilding scenarios benefited more species adapted to habitats with low human influence, such as forests and open woodlands. The highest functional richness and dispersion was predicted for rewilding scenarios, which could improve landscape restoration and provide opportunities for the expansion and recolonization of forest areas by native species. The recovery of traditional farming and agroforestry activities results in the lowest values of functional richness, but these strategies contribute to complex landscape matrices with diversified habitats and resources. Moreover, this strategy could offer opportunities for fire suppression and increase landscape fire resistance. An integrative approach reconciling rewilding initiatives with the recovery of extensive agricultural and agroforestry activities is potentially an harmonious strategy for supporting the provision of ecosystem services while securing biodiversity conservation and functional diversity within the natural park.

Aboveground and belowground biodiversity have complementary effects on ecosystem functions across global grasslands

Grasslands are integral to maintaining biodiversity and key ecosystem services and are under threat from climate change. Plant and soil microbial diversity, and their interactions, support the provision of multiple ecosystem functions (multifunctionality). However, it remains virtually unknown whether plant and soil microbial diversity explain a unique portion of total variation or shared contributions to supporting multifunctionality across global grasslands. Here, we combine results from a global survey of 101 grasslands with a novel microcosm study, controlling for both plant and soil microbial diversity to identify their individual and interactive contribution to support multifunctionality under aridity and experimental drought. We found that plant and soil microbial diversity independently predict a unique portion of total variation in above- and belowground functioning, suggesting that both types of biodiversity complement each other. Interactions between plant and soil microbial diversity positively impacted multifunctionality including primary production and nutrient storage. Our findings were also climate context dependent, since soil fungal diversity was positively associated with multifunctionality in less arid regions, while plant diversity was strongly and positively linked to multifunctionality in more arid regions. Our results highlight the need to conserve both above- and belowground diversity to sustain grassland multifunctionality in a drier world and indicate climate change may shift the relative contribution of plant and soil biodiversity to multifunctionality across global grasslands.

Identifying global conservation priorities for terrestrial vertebrates based on multiple dimensions of biodiversity

The Kunming-Montreal Global Biodiversity Framework of the Convention on Biological Diversity calls for an expansion of the current protected areas (PAs) to cover at least 30% of global land and water areas by 2030 (i.e., the 30×30 target). Efficient spatial planning for PA expansion is an urgent need for global conservation practice. A spatial prioritization framework considering multiple dimensions of biodiversity is critical for improving the efficiency of the spatial planning of PAs, yet it remains a challenge. We developed an index for the identification of priority areas based on functionally rare, evolutionarily distinct, and globally endangered species (FREDGE) and applied it to 21,536 terrestrial vertebrates. We determined species distributions, conservation status (global endangerment), molecular phylogenies (evolutionary distinctiveness), and life-history traits (functional rarity). Madagascar, Central America, and the Andes were of high priority for the conservation of multiple dimensions of terrestrial vertebrate biodiversity. However, 68.8% of grid cells in these priority areas had <17% of their area covered by PAs, and these priority areas were under intense anthropogenic and climate change threats. These results highlight the difficulties of conserving multiple dimensions of biodiversity. Our global analyses of the geographical patterns of multiple dimensions of terrestrial vertebrate biodiversity demonstrate the insufficiency of the conservation of different biodiversity dimensions, and our index, based on multiple dimensions of biodiversity, provides a useful tool for guiding future spatial prioritization of PA expansion to achieve the 30×30 target under serious pressures.

Soil enzyme profile analysis for indicating decomposer micro-food web

Highly diverse exoenzymes mediate the energy flow from substrates to the multitrophic microbiota within the soil decomposer micro-food web. Here, we used a "soil enzyme profile analysis" approach to establish a series of enzyme profile indices; those indices were hypothesized to reflect micro-food web features. We systematically evaluated the shifts in enzyme profile indices in relation to the micro-food web features in the restoration of an abandoned cropland to a natural area. We found that enzymatic C:N stoichiometry and decomposability index were significantly associated with substrate availability. Furthermore, the higher Shannon diversity index in the exoenzyme profile, especially for the C-degrading hydrolase, corresponded to a greater microbiota community diversity. The increased complexity and stability of the exoenzyme network reflected similar changes with the micro-food web networks. In addition, the gross activity of the enzyme profile as a parameter for soil multifunctionality, effectively predicted the substrate content, microbiota community size, diversity, and network complexity. Ultimately, the proposed enzymic channel index was closely associated with the traditional decomposition channel indices derived from microorganisms and nematodes. Our results showed that soil enzyme profile analysis reflected very well the decomposer food web features. Our study has important implications for projecting future climate change or anthropogenic disturbance impacts on soil decomposer micro-food web features by using soil enzyme profile analysis.

Shifts in bacterial traits under chronic nitrogen deposition align with soil processes in arbuscular, but not ectomycorrhizal-associated trees

Nitrogen (N) deposition increases soil carbon (C) storage by reducing microbial activity. These effects vary in soil beneath trees that associate with arbuscular (AM) and ectomycorrhizal (ECM) fungi. Variation in carbon C and N uptake traits among microbes may explain differences in soil nutrient cycling between mycorrhizal associations in response to high N loads, a mechanism not previously examined due to methodological limitations. Here, we used quantitative Stable Isotope Probing (qSIP) to measure bacterial C and N assimilation rates from an added organic compound, which we conceptualize as functional traits. As such, we applied a trait-based approach to explore whether variation in assimilation rates of bacterial taxa can inform shifts in soil function under chronic N deposition. We show taxon-specific and community-wide declines of bacterial C and N uptake under chronic N deposition in both AM and ECM soils. N deposition-induced reductions in microbial activity were mirrored by declines in soil organic matter mineralization rates in AM but not ECM soils. Our findings suggest C and N uptake traits of bacterial communities can predict C cycling feedbacks to N deposition in AM soils, but additional data, for instance on the traits of fungi, may be needed to connect microbial traits with soil C and N cycling in ECM systems. Our study also highlights the potential of employing qSIP in conjunction with trait-based analytical approaches to inform how ecological processes of microbial communities influence soil functioning.

Linking functional composition moments of the sub-Mediterranean ecotone with environmental drivers

Introduction

Functional trait-based approaches are extensively applied to the study of mechanisms governing community assembly along environmental gradients. These approaches have been classically based on studying differences in mean values among species, but there is increasing recognition that alternative metrics of trait distributions should be considered to decipher the mechanisms determining community assembly and species coexistence. Under this framework, the main aim of this study is to unravel the effects of environmental conditions as drivers of plant community assembly in sub-Mediterranean ecotones.

Methods

We set 60 plots in six plant communities of a sub-Mediterranean forest in Central Spain, and measured key above- and belowground functional traits in 411 individuals belonging to 19 species, along with abiotic variables. We calculated community-weighted mean (CWM), skewness (CWS) and kurtosis (CWK) of three plant dimensions, and used maximum likelihood techniques to analyze how variation in these functional community traits was driven by abiotic factors. Additionally, we estimated the relative contribution of intraspecific trait variability and species turnover to variation in CWM.

Results and discussion

The first three axes of variation of the principal component analyses were related to three main plant ecological dimensions: Leaf Economics Spectrum, Root Economics Spectrum and plant hydraulic architecture, respectively. Type of community was the most important factor determining differences in the functional structure among communities, as compared to the role of abiotic variables. We found strong differences among communities in their CWMs in line with their biogeographic origin (Eurosiberian vs Mediterranean), while differences in CWS and CWK indicate different trends in the functional structure among communities and the coexistence of different functional strategies, respectively. Moreover, changes in functional composition were primarily due to intraspecific variability.

Conclusion

We observed a high number of strategies in the forest with the different communities spreading along the acquisitive-conservative axis of resource-use, partly matching their Eurosiberian-Mediterranean nature, respectively. Intraspecific trait variability, rather than species turnover, stood as the most relevant factor when analyzing functional changes and assembly patterns among communities. Altogether, our data support the notion that ecotones are ecosystems where relatively minor environmental shifts may result in changes in plant and functional composition.

Plant biodiversity offsets negative effects of metals and metalloids soil multi-contamination on ecosystem multifunctionality

Despite increasing metals and metalloids (MM) human-driven soil contamination, how it simultaneously alters biodiversity and ecosystem functioning remains unknown. We used a wide gradient of a 170-year-old MM soil multi-contamination in Mediterranean scrublands to assess the effects of soil multi-contamination on multiple plant biodiversity facets, microbial communities and ecosystem multifunctionality (EMF). We found an overall positive effect of plant biodiversity on EMF mediated by microbial communities, and allowing offsetting the negative impacts of MM soil multi-contamination, especially on soil water holding capacity and nitrogen content. The diversity of distant plant lineages was the key facet promoting EMF by enhancing microbial communities, whereas the subordinate species richness altered EMF. By developing a holistic approach of these complex relationships between soil multi-contamination, plant biodiversity, microbial communities and ecosystem functioning, our results reveal the potential of plant biodiversity, and especially the diversity of evolutionary distant species, to offset the alteration of ecosystem functioning by MM soil multi-contamination. In this worldwide decade of ecosystems restoration, our study helps to identify relevant facets of plant biodiversity promoting contaminated ecosystem functioning, which is crucial to guide and optimize management efforts aiming to restore ecosystems and preserve human health.

Relationships of stomatal morphology to the environment across plant communities

The relationship between stomatal traits and environmental drivers across plant communities has important implications for ecosystem carbon and water fluxes, but it has remained unclear. Here, we measure the stomatal morphology of 4492 species-site combinations in 340 vegetation plots across China and calculate their community-weighted values for mean, variance, skewness, and kurtosis. We demonstrate a trade-off between stomatal density and size at the community level. The community-weighted mean and variance of stomatal density are mainly associated with precipitation, while that of stomatal size is mainly associated with temperature, and the skewness and kurtosis of stomatal traits are less related to climatic and soil variables. Beyond mean climate variables, stomatal trait moments also vary with climatic seasonality and extreme conditions. Our findings extend the knowledge of stomatal trait-environment relationships to the ecosystem scale, with applications in predicting future water and carbon cycles.

Direct and indirect effects of dominant plants on ecosystem multifunctionality

Biodiversity is essential for the provision of multiple ecosystem functions simultaneously (ecosystem multifunctionality EMF). Yet, it remains unclear whether and how dominant plant species impact EMF. Here, we aimed at disentangling the direct from indirect above- and belowground pathways by which dominant plant species influence EMF. We evaluated the effects of two dominant plant species (Dasiphora fruticosa, and the toxic perennial plant Ligularia virgaurea) with expected positive and negative impacts on the abiotic environment (soil water content and pH), surrounding biological communities (plant and nematode richness, biomass, and abundance in the vicinity), and on the EMF of alpine meadows, respectively. We found that the two dominant plants enhanced EMF, with a positive effect of L. virgaurea on EMF greater than that of D. fruticosa. We also observed that dominant plants impacted on EMF through changes in soil water content and pH (indirect abiotic effects), but not through changes in biodiversity of surrounding plants and nematodes (indirect biotic pathway). Our study suggests that dominant plants may play an important role in promoting EMF, thus expanding the pervasive mass-ratio hypothesis originally framed for individual functions, and could mitigate the negative impacts of vegetation changes on EMF in the alpine meadows.

How to improve scaling from traits to ecosystem processes

Scaling approaches in ecology assume that traits are the main attributes by which organisms influence ecosystem functioning. However, several recent empirical papers have found only weak links between traits and ecosystem functioning, questioning the usefulness of trait-based ecology (TBE). We argue that these studies often suffer from one or more widespread misconceptions. Specifically, these studies often (i) conflict with the conceptual foundations of TBE, (ii) lack theory- or hypothesis-driven selection and use of traits, (iii) tend to ignore intraspecific variation, and (iv) use experimental or study designs that are not well suited to make strong tests of TBE assumptions. Addressing these aspects could significantly improve our ability to scale from traits to ecosystem functioning.

Temporal change in functional rarity in marine fish assemblages

Recent research has uncovered rapid compositional and structural reorganization of ecological assemblages, with these changes particularly evident in marine ecosystems. However, the extent to which these ongoing changes in taxonomic diversity are a proxy for change in functional diversity is not well understood. Here we focus on trends in rarity to ask how taxonomic rarity and functional rarity covary over time. Our analysis, drawing on 30 years of scientific trawl data, reveals that the direction of temporal shifts in taxonomic rarity in two Scottish marine ecosystems is consistent with a null model of change in assemblage size (i.e. change in numbers of species and/or individuals). In both cases, however, functional rarity increases, as assemblages become larger, rather than showing the expected decrease. These results underline the importance of measuring both taxonomic and functional dimensions of diversity when assessing and interpreting biodiversity change.

Habitat loss predicts the functional extinction of fish from Amazonian streams during the Anthropocene

The evaluation of extinction risk has typically focused on individual species, although a shift to a focus on ecosystem functioning would appear to be an urgent priority for conservation planning, especially considering that a sixth mass extinction event has already begun. In the present study, we investigated how fish extinction driven by habitat loss may modify the functioning of freshwater Amazonian ecosystems. We sampled the fish and environmental conditions of 63 streams in the eastern Amazon and simulated extinction based on the vulnerability of the species to habitat loss, which is the principal threat to tropical biodiversity. The simulated extinction of vulnerable species led to a decrease in both the mean body size of the community and functional rarity and culminated in abrupt losses of ecosystem functions after 5% and 10% of extinction at local and regional scales. Our functional approach demonstrated the progressive loss of ecological functions in Amazon streams, which may collapse altogether following the extinction of functions related to protection against biological invasions, and associated alterations in nutrient cycling and water quality. We provide robust predictions on the modification of the ecosystem following the extinction of fish species, which is a major step toward the development of effective conservation measures that ensure the avoidance of the predicted processes, and help to prevent the loss of biodiversity and the potentially irreversible modifications to ecosystem functioning.

Size matters for linking traits to ecosystem multifunctionality

A priority research field addresses how to optimize diverse ecosystem services to people, including biodiversity support, regulatory, utilitarian and cultural services. This field may benefit from linking ecosystem services to the sizes of different body parts of organisms, with functional traits as the go-between. Using woody ecosystems to explore such linkages, we hypothesize that across stem diameter classes from trunk via branches to twigs, key wood and bark functional traits (especially those defining size-shape and resource economics spectra) vary both within individual trees and shrubs and across woody species, thereby together boosting ecosystem multifunctionality. While we focus on woody plants aboveground, we discuss promising extensions to belowground organs of trees and shrubs and analogs with other organisms, for example, vertebrate animals.

The biological sink of atmospheric H2 is more sensitive to spatial variation of microbial diversity than N2O and CO2 emissions in a winter cover crop field trial

The integration of winter cover crop (WCC) in culture rotations promotes multiple ecosystem services, but concomitant microbial diversity and functioning responses in soil have received less attention. A field trial was established to test the hypothesis that enhanced crop diversity with the integration of WCC in a conventional maize-soy rotation promotes microbial diversity and the biological sink of H₂ in soil, while reducing N₂O emissions to the atmosphere. Vicia villosa (hairy vetch), Avena sativa (oat), and Raphanus sativus (Daikon radish) were cultivated alone or in combinations and flux measurements were performed throughout two subsequent growing seasons. Soil acted as a net sink for H₂ and as a net source for CO₂ and N₂O. CO₂ flux was the most sensitive to WCC whereas a significant spatial variation was observed for H₂ flux with soil uptake rates observed in the most productive area two-fold greater than the baseline level. Sequencing and quantification of taxonomic and functional genes were integrated to explain variation in trace gas fluxes with compositional changes in soil microbial communities. Fungal communities were the most sensitive to WCC, but neither community abundance nor beta diversity were found to be indicative of fluxes. The alpha diversity of taxonomic and functional genes, expressed as the number of effective species, was integrated into composite variables extracted from multivariate analyses. Only the composite variable computed with the inverse Simpson's index displayed a reproducible pattern throughout both growing seasons, with functional genes and bacterial 16S rRNA gene defining the two most contrasting gradients. The composite variable was decoupled from WCC treatment and explained 19-20% spatial variation of H₂ fluxes. The coupling of composite alpha diversity metrics derived from multiple genes with soil processes warrants further investigations to implement novel indicators of soil health in response to changing management practices at the local scale.

Functionally distinct tree species support long-term productivity in extreme environments

Despite evidence of a positive effect of functional diversity on ecosystem productivity, the importance of functionally distinct species (i.e. species that display an original combination of traits) is poorly understood. To investigate how distinct species affect ecosystem productivity, we used a forest-gap model to simulate realistic temperate forest successions along an environmental gradient and measured ecosystem productivity at the end of the successional trajectories. We performed 10 560 simulations with different sets and numbers of species, bearing either distinct or indistinct functional traits, and compared them to random assemblages, to mimic the consequences of a regional loss of species. Long-term ecosystem productivity dropped when distinct species were lost first from the regional pool of species, under the harshest environmental conditions. On the contrary, productivity was more dependent on ordinary species in milder environments. Our findings show that species functional distinctiveness, integrating multiple trait dimensions, can capture species-specific effects on ecosystem productivity. In a context of an environmentally changing world, they highlight the need to investigate the role of distinct species in sustaining ecosystem processes, particularly in extreme environmental conditions.

Using geometry to rank evenness measures: Towards a deeper understanding of divergence

While recent work has established divergence as a key framework for understanding evenness, there is currently no research exploring how the families of measures within the divergence-based framework relate to each other. This paper uses geometry to show that, holding order and richness constant, the families of divergence-based evenness measures nest. This property allows them to be ranked based on their reactivity to changes in relatively even assemblages or changes in relatively uneven ones. We establish this ranking and explore how the distance-based measures relate to it for both order q = 2 and q = 1. We also derive a new family of distance-based measures that captures the angular distance between the vector of relative abundances and a perfectly even vector and is highly reactive to changes in even assemblages. Finally, we show that if we only require evenness to be a divergence, then any smooth, monotonically increasing function of diversity can be made into an evenness measure. A deeper understanding of how to measure evenness will require empirical or theoretical research that uncovers which kind of divergence best reflects the underlying concept.

Variation Patterns of Functional Trait Moments Along Geographical Gradients and Their Environmental Determinants in the Subtropical Evergreen Broadleaved Forests

Understanding how environmental change alters the composition of plant assemblages is a major challenge in the face of global climate change. Researches accounting for site-specific trait values within forest communities help bridge plant economics theory and functional biogeography to better evaluate and predict relationships between environment and ecosystem functioning. Here, by measuring six functional traits (specific leaf area, leaf dry matter content, leaf nitrogen, and phosphorus concentration, leaf nitrogen/phosphorus, wood density) for 292 woody plant species (48,680 individuals) from 250 established permanent forest dynamics plots in five locations across the subtropical evergreen broadleaved forests (SEBLF) in China, we quantified functional compositions of communities by calculating four trait moments, i.e., community-weighted mean, variance, skewness, and kurtosis. The geographical (latitudinal, longitudinal, and elevational) patterns of functional trait moments and their environmental drivers were examined. Results showed that functional trait moments shifted significantly along the geographical gradients, and trait moments varied in different ways across different gradients. Plants generally showed coordinated trait shifts toward more conservative growth strategies (lower specific leaf area, leaf N and P concentration while higher leaf nitrogen/phosphorus and wood density) along increasing latitude and longitude. However, trends opposite to the latitudinal and longitudinal patterns appeared in trait mean values along elevation. The three sets of environmental variables (climate, soil and topography) explained 35.0-69.0%, 21.0-56.0%, 14.0-31.0%, and 16.0-30.0% of the variations in mean, variance, skewness, and kurtosis across the six functional traits, respectively. Patterns of shifts in functional trait moments along geographical gradients in the subtropical region were mainly determined by the joint effects of climatic and edaphic conditions. Climate regimes, especially climate variability, were the strongest driving force, followed by soil nutrients, while topography played the least role. Moreover, the relationship of variance, skewness and kurtosis with climate and their geographical patterns suggested that rare phenotypes at edges of trait space were selected in harsher environments. Our study suggested that environmental filtering (especially climate variability) was the dominant process of functional assembly for forest communities in the subtropical region along geographical gradients.

Unveiling ecological assembly rules from commonalities in trait distributions

Deciphering the effect of neutral and deterministic processes on community assembly is critical to understand and predict diversity patterns. The information held in community trait distributions is commonly assumed as a signature of these processes, but empirical and modelling attempts have most often failed to untangle their confounding, sometimes opposing, impacts. Here, we simulated the assembly of trait distributions through stochastic (dispersal limitation) and/or deterministic scenarios (environmental filtering and niche differentiation). We characterized the shape of trait distributions using the skewness-kurtosis relationship. We identified commonalities in the co-variation between the skewness and the kurtosis of trait distributions with a unique signature for each simulated assembly scenario. Our findings were robust to variation in the composition of regional species pools, dispersal limitation and environmental conditions. While ecological communities can exhibit a high degree of idiosyncrasy, identification of commonalities across multiple communities can help to unveil ecological assembly rules in real-world ecosystems.

The dimensionality and structure of species trait spaces

Trait-based ecology aims to understand the processes that generate the overarching diversity of organismal traits and their influence on ecosystem functioning. Achieving this goal requires simplifying this complexity in synthetic axes defining a trait space and to cluster species based on their traits while identifying those with unique combinations of traits. However, so far, we know little about the dimensionality, the robustness to trait omission and the structure of these trait spaces. Here, we propose a unified framework and a synthesis across 30 trait datasets representing a broad variety of taxa, ecosystems and spatial scales to show that a common trade-off between trait space quality and operationality appears between three and six dimensions. The robustness to trait omission is generally low but highly variable among datasets. We also highlight invariant scaling relationships, whatever organismal complexity, between the number of clusters, the number of species in the dominant cluster and the number of unique species with total species richness. When species richness increases, the number of unique species saturates, whereas species tend to disproportionately pack in the richest cluster. Based on these results, we propose some rules of thumb to build species trait spaces and estimate subsequent functional diversity indices.