What can intraspecific trait variability tell us about fungal communities and adaptations?

Analyses of species functional traits are suitable to better understand the coexistence of species in a given environment. Trait information can be applied to investigate diversity patterns along environmental gradients and subsequently to predict and mitigate threats associated with climate change and land use. Species traits are used to calculate community trait means, which can be related to environmental gradients. However, while species traits can provide insights into the mechanisms underlying community assembly, they can lead to erroneous inferences if mean trait values are used. An alternative is to incorporate intraspecific trait variability (ITV) into calculating the community trait means. This approach gains increasing acceptance in plant studies. For macrofungi, functional traits have recently been applied to examine their community ecology but, to our knowledge, ITV has yet to be incorporated within the framework of community trait means. Here, we present a conceptual summary of the use of ITV to investigate the community ecology of macrofungi, including the underlying ecological theory. Inferences regarding community trait means with or without the inclusion of ITV along environmental gradients are compared. Finally, an existing study is reconsidered to highlight the variety of possible outcomes when ITV is considered. We hope this Opinion will increase awareness of the potential for within-species trait variability and its importance for statistical inferences, interpretations, and predictions of the mechanisms structuring communities of macro- and other fungi.


Introduction
Climate change and land use dramatically affect the biodiversity within ecosystems (Fahrig 2003;Lindner et al. 2010). However, our restricted understanding of the mechanisms underlying species diversity patterns (Vellend 2010) limits our ability to predict and mitigate threats to biodiversity and, further, to ecosystem function (e.g., Urban et al. 2016). A primary goal of community ecology is to elucidate the mechanisms that structure communities along environmental gradients, such as species sorting, biotic interactions, or legacy effects (Vellend 2010). Analyses of functional traits can help to predict community change in response to environmental change across spatial scales. In contrast to fungal community ecology, the use of species traits have a long history in animal and plant community ecology (Grime 2006;Westoby and Wright 2006;Shipley et al. 2016). Trait-based approaches allow an understanding of environmental filters' relative importance and biotic interactions in community assembly (McGill et al. 2006). Today, a large set of tools and methods are available to analyze traits within a community ecological framework (D'Amen et al. 2017;Weiher 2011;de Bello et al. 2021). However, the inclusion of intraspecific trait variability is still rare and a growing body of literature shows that a reliance on means of species trait can lead to (i) an underestimation of intraspecific trait changes, (ii) incorrect inferences regarding the correlations of particular traits with environmental gradients, and (iii) incorrect inferences of assembly processes, e.g., in tests against null-model-based neutral assembly (Jung et al. This article is part of the Topical collection on Basidiomycote Mycology in honor of Franz Oberwinkler who passed away in March 2018 Section Editor: Dominik Begerow 2010; Violle et al. 2012;Bolnick et al. 2011;Fang et al. 2019;Spasojevic et al. 2016;Shipley et al. 2016). Based on a common and widely used method (community mean), we first illustrate how traits can enter a community analysis and demonstrate the effect the consideration of intraspecific trait variability (ITV) can have on inferences based on this method.

Community-fixed versus community-specific trait means
Although other methods exist to study the link between traits, communities, and the environment (Weiher 2011), the focus here is on the community trait mean as a simple and intuitive mathematical framework with wide applications. After the quantification of species occurrences on plots/sites within a standardized study design, a community matrix is digitalized. To calculate the community trait mean, species traits are assembled (Shipley et al. 2016). Species traits in this framework are typically obtained either from literature or from the measurement of traits during field work. Since the measurement of all individuals is labor-intensive, in plant studies, individuals are often only measured on a few plots to obtain a species trait mean (Shipley et al. 2016). This is especially the case for traits not available from the literature (Shipley et al. 2016). In either case, a fixed, site-independent trait mean of each species is used to calculate community trait means, the "community- where p i is the proportion (e.g., relative abundance, biomass) of the ith species in a community (site, plot, region, etc.), S is the number of species in that community, and x i is the mean trait (e.g., literature-based species trait) of the ith species (Lepš et al. 2011). Trait changes along an environmental gradient are thus studied via changes in the community-fixed mean along the gradient and any change in traits is then necessarily attributed to species turnover. However, besides species turnover, trait changes may also or instead reflect intraspecific trait variability (ITV, Fig. 1a) (Lepš et al. 2011), caused by local adaptation or phenotypic plasticity (Bolnick et al. 2011;Violle et al. 2012). To account for ITV, rather than quantifying the community-fixed means, changes in species turnover and ITV are quantified. The field observations and resulting community matrix nonetheless follow a similar approach. Specific, site-dependent traits are measured for individuals on each plot. Thus, in terms of calculating a site-independent community trait mean, a site-dependent community trait mean is calculated (hereafter "community-specific mean"), as , where x i_site is the site-dependent trait value of the ith species on the respective site (Lepš et al. 2011). This individual trait potentially differs for different sites, characterized by different environmental conditions. Note that based on data that allows calculating the community-specific mean, also the community-fixed mean can be obtained by averaging the site-dependent trait values for each species. This approach allows to compare community-fixed with communityspecific means along a gradient and thus to disentangle the relative effects of species turnover and additionally individual trait turnover (e.g., using the R package cati, Taudiere and Violle (2016)). Many trait-based studies have demonstrated the importance of including ITV in plant community ecology. First, a meta-analysis of 629 plant communities and 36 plant traits revealed that, compared to inter-species variability, the ITV was higher than expected (Siefert et al. 2015). Second, in several studies, the application of community-fixed andspecific means showed that trait change is a combined effect of species turnover and intraspecific trait turnover (Jung et al. 2010;Kichenin et al. 2013;Lepš et al. 2011). Moreover, community-specific means can change substantially, independent of a change in the community-fixed means (Kichenin et al. 2013;Jung et al. 2010). Jung et al. (2010), for example, sampled traits (e.g., specific leaf area or leaf dry matter content) of 20 dominant plant species along an experimental drought gradient (rainfall exclusion), focusing on the strength and direction of the responses of community-fixed and -specific means. They found that in addition to species turnover, there was a large contribution of the intraspecific trait turnover. For some traits (e.g., leaf carbon concentration), there was no effect of the community-fixed mean; however, the specific trait mean responded strongly to drought. The authors concluded that considering community-fixed but not community-specific means leads to an underestimation of responses and probably to erroneous inferences about assembly processes and predictions. As an example, community responses to climate change may not be understood in their full range (Lepš et al. 2011) (Fig. 1b). Therefore, ITV is increasingly being considered in plant community ecology (Jung et al. 2010;Kichenin et al. 2013;Moran et al. 2016;Chalmandrier et al. 2017;Niu et al. 2020), with a growing body of theoretical considerations (Bolnick et al. 2011;Spasojevic et al. 2016;Violle et al. 2012;Albert et al. 2011).
Trait-based approaches have been applied only recently in studies of fungal community ecology, mostly to explore environment-fungi relationships (Aguilar-Trigueros et al. 2015;Crowther et al. 2014;Zanne et al. 2020). Mycelium and spore traits have been studied in Glomeromycota (Munkvold et al. 2004;Ho et al. 2017;Giauque et al. 2019), and mycelium, spore, and fruit body traits in Agaricomycetes (Maynard et al. 2019;Bässler et al. 2016a;Abrego et al. 2016;Krah et al. 2019). Independent of the particular method employed to examine traits within a community ecology framework, to our knowledge, only community-fixed but not community-specific means were applied. We, therefore, briefly review the existing studies of functional traits of macrofungi (fruit body-forming fungi), before demonstrating the importance of incorporating ITV in analyses of community trait means along environmental gradients.
The use of functional traits in community ecology of macrofungi Functional traits are increasingly being used to examine the community ecology of macrofungi, but such studies are still rare (Dawson et al. 2018). The majority of studies investigated traits of the sexually reproductive organ (fruit body), e.g., mushroom and spore size. The advantage of focusing on mushroom and spore traits is that they are easily available because they form an important basis of species identification. Less frequent are studies using mycelium traits, but we are not aware of studies that considered both fruit body and mycelium traits. Investigations of fungal traits have mainly been conducted at a local to landscape scale. Studies addressing fruit body traits showed that the mean fruit body size of assemblages increases with the availability of resources (Bässler et al. 2016a) and that total fruit body biomass is mainly affected by interannual precipitation (Alday et al. 2017). Smaller spores were associated with open shrub-like vegetation and larger spores with forests (Crandall et al. 2020). Forest management was suggested to act as a habitat filter sorting woodinhabiting fungi based on their fruit body characteristics (Bässler et al. 2014), and fruit body size of perennial species was larger in disturbed sites (Bässler et al. 2016b). Resource depletion increased spore elongation of assemblages, Fig. 1 Graphical summary of the importance of intraspecific trait variability (ITV) in community ecology studies of macrofungi. a Trait changes along environmental gradients may result from species turnover, individual turnover, or both. We currently do not know the relative importance of species turnover vs. individual trait turnover in macrofungi, limiting our ability to predict and thus mitigate threats to biodiversity, e.g., caused by climate change. b Individual traits may be more extreme than species trait means, thus influencing inferences of assembly processes. For example, even if species shift, individual traits may compensate for lost species traits and thus the community-specific trait mean may not change. c As an example: a previous study used the community-fixed mean of fruit body color lightness, which changed together with the mean average temperature across Europe (Krah et al. 2019). However, since species trait means were used in that study, whether trait change due to species turnover is strengthened, weakened, or neutralized can only be determined by the consideration of ITV (numbers indicate scenarios that are discussed in the main text). Thus, predictions regarding trait changes in response to climate warming may be over-or underestimated when based on community-fixed means alone probably to overcome dispersal limitation (Halbwachs et al. 2017). A study of fruit body types showed differing responses along with a forest management gradient (Abrego et al. 2016). Another study examined multiple traits and found that they correlate with host preferences of wood-inhabiting fungi (Purhonen et al. 2020). In contrast to studies on the local to landscape scale, studies on the continental to global scale are less frequent. Across Norway, species that fruited earlier and more towards continental climates had larger spore size (Kauserud et al. 2011). A study that assessed a variety of reproductive traits (e.g., fruit body and spore size) found that the latitudinal abundance of forest fungi across Great Britain could be predicted by their functional traits (Gange et al. 2018). One study investigated mushroom size globally and found larger mushrooms in temperate and smaller mushrooms in the tropics and boreal zone (Bässler et al. 2021). This effect was attributed to differences in seasonal resource input between biomes or thermoregulation. In a macroclimatic European-scale study, fruit body color lightness was shown to structure the examined fungal communities and followed a mean temperature gradient (Krah et al. 2019). In detail, it was found that mushroom assemblages are darker in cold environments, probably due to thermoregulatory mechanisms (thermal melanism theory). However, there was no evidence of changes in the community-fixed mean color lightness over time , despite expectations that the communities would have already become lighter in response to climate warming ("thermophilization," see Zellweger et al. (2020)). In parallel with the publication of trait-based studies, trait handbooks and guidelines for trait measurement were developed that drew upon plant functional trait frameworks (Dawson et al. 2018;Aguilar-Trigueros et al. 2015;Zanne et al. 2020) and thereby highlighted the importance of studying fungi based on their functional traits (Pringle et al. 2015).
This brief summary of studies and results demonstrates the value of studying fungal traits in a geographic-environmental context as well as their competitive abilities. The important insights provided by the cited studies are, however, limited because we do not know the within-species responses.

How can ITV provide better insights into fungal community ecology?
The timeline of the above-summarized studies reveals that community ecologists, studying macrofungi, only recently became aware of the advantages of trait-based analyses, in contrast to the much longer history of this approach in plant science. However, as outlined above, interspecific, siteindependent traits are not able to resolve the full spectrum of trait change, i.e., intraspecific trait turnover (Shipley et al. 2016). Although preliminary and based on a small dataset a recent study investigating mushroom size variability found high interspecific variability, demonstrating the potential effect ITV may have on community analysis (Dawson and Jönsson 2020). Thus, we strongly encourage ecologists studying fungal communities to revisit their hypotheses regarding macrofungi, by taking into account ITV and including individual trait measurements in field inventories.
The above example of the community-fixed mean color lightness can be used to illustrate the limitations of approaches based on the fixed community means alone. This example can also be used to demonstrate how to account for ITV in future studies. For example, the color lightness of fungal fruit bodies should be determined for at least three individuals per site per species, such as by taking standardized photos or reflection spectra on all plots across a macroclimatic gradient. Then both a community-fixed (black line in Fig. 1c) and communityspecific trait mean are calculated as described above (blue lines in Fig. 1c).
Based on the inclusion of ITV, three plausible responses of the specific trait mean are proposed as an example, although more are theoretically possible: (1) Individuals in cold climates are darker than mean color trait values, resulting in a more pronounced response of the community-specific mean than of the community-fixed mean.
(2) Individual traits do not deviate substantially from the mean lightness trait values and thus, the responses based on the community-fixed and community-specific means are not significantly different. (3) Individual traits deviate strongly from the mean color trait values so that individuals are more light-colored in cold environments and more dark-colored in warm climates. This scenario would lead to a change of the sign of the communityspecific compared to the community-fixed mean (Fig. 1c). As outlined above, Krah et al. (2019) did not find a shift in community-fixed lightness over time, but they did observe a shift in community-fixed lightness with mean annual temperature. One possible explanation may be that species already responded via their intraspecific color lightness, however, not yet via species shifts over time.

Conclusion
Despite recent recognition of the value of trait-based approaches in community ecology studies of macrofungi (Dawson et al. 2018), the pure focus on species mean trait values to calculate community means may lead to erroneous inferences, interpretations, and predictions of, e.g., the effects of climate change on fungal communities. We, therefore, encourage ecologists studying fungal communities to revise their hypotheses to include ITV in the calculation of community means. Therefore, measuring trait values in situ from a large number of individuals will be necessary in most cases if the retrospective incorporation of ITV into analyses is not possible.
We suggest a combination of field surveys/experiments with a focus on intraspecific trait values that is complemented by eco-physiological laboratory experiments. Sampling in the laboratory allows traits to be evaluated under standardized conditions. A potential model species for laboratory experiments is the cultivatable and edible mushroom Agrocybe aegerita (Herzog et al. 2019). The above-noted absence of studies simultaneously addressing mycelium and fruit body traits suggests an exciting avenue, especially for future laboratory studies. Finally, although we have discussed the potential pitfalls of ignoring ITV in the context of macrofungal community ecology, the same argument can be made for fungal lineages without conspicuous fruit bodies and/or nonmorphological traits (e.g., molecular traits). The field of fungal community ecology would thus benefit from a unified theoretical concept that incorporates ITV into studies beyond those of macrofungi or morphological traits.
Author contribution FK had the initial idea of the study and wrote the first draft of the manuscript. CB contributed development of the concept and to the writing of the manuscript.
Funding Open access funding enabled and organized by Projekt DEAL.

Declarations
Availability of data and material (data transparency) Not applicable.
Code availability (software application or custom code) Not applicable.
Competing interests The authors declare that they have no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.