Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 20 April 2021
Sec. Terrestrial Microbiology

Exploring Trait Trade-Offs for Fungal Decomposers in a Southern California Grassland

  • 1Ecology and Evolutionary Biology, University of California Irvine, Irvine, CA, United States
  • 2Department of Earth System Science, University of California Irvine, Irvine, CA, United States
  • 3Department of Microbiology and Plant Pathology, University of California, Riverside, CA, United States

Fungi are important decomposers in terrestrial ecosystems, so their responses to climate change might influence carbon (C) and nitrogen (N) dynamics. We investigated whether growth and activity of fungi under drought conditions were structured by trade-offs among traits in 15 fungal isolates from a Mediterranean Southern California grassland. We inoculated fungi onto sterilized litter that was incubated at three moisture levels (4, 27, and 50% water holding capacity, WHC). For each isolate, we characterized traits that described three potential lifestyles within the newly proposed “YAS” framework: growth yield, resource acquisition, and stress tolerance. Specifically, we measured fungal hyphal length per unit litter decomposition for growth yield; the potential activities of the extracellular enzymes cellobiohydrolase (CBH), β-glucosidase (BG), β-xylosidase (BX), and N-acetyl-β-D-glucosaminidase (NAG) for resource acquisition; and ability to grow in drought vs. higher moisture levels for drought stress tolerance. Although, we had hypothesized that evolutionary and physiological trade-offs would elicit negative relationships among traits, we found no supporting evidence for this hypothesis. Across isolates, growth yield, drought stress tolerance, and extracellular enzyme activities were not significantly related to each other. Thus, it is possible that drought-induced shifts in fungal community composition may not necessarily lead to changes in fungal biomass or decomposer ability in this arid grassland.

Introduction

We are just beginning to understand how fungal growth and activity might shift under climate change (Todd-Brown et al., 2014; Xiao et al., 2018; Maynard et al., 2019). Fungal communities may transform to favor certain groups more adapted to these new environmental conditions (Kivlin and Treseder, 2014; Boddy et al., 2016). However, which fungal groups will be proliferating is still unknown. This issue is important in predicting decomposition, and ultimately, soil carbon (C) storage (Wieder et al., 2013). Specifically in the southwestern United States, which is experiencing major changes in precipitation patterns (Seager et al., 2007; Gibson et al., 2020), understanding how fungi will respond to drought is critical. One way to predict fungal response to drought stress is to investigate trade-offs among fungal traits.

Evolutionary and physiological trade-offs among traits may structure microbial communities and their contributions to ecosystem function (Ho et al., 2013; Malik et al., 2020b). Thus, examining trade-offs among fungal traits may be a useful approach in which to assess fungal response to drought. For example, under stressful conditions like drought, allocation of finite resources within organisms might require investment in one function at the expense of another function (e.g., production of compounds to combat stress at the cost of deferred growth; Bennett and Lenski, 2007). Here, we determined what trade-offs exist among traits of decomposer fungi in a Southern California grassland. Knowledge of these trade-offs may improve parameterization of trait-based microbial decomposition models as well as understanding of fungal physiology in this changing ecosystem.

We examined traits associated with three fungal strategies: growth yield, resource acquisition, and stress tolerance. Potential trade-offs between these life history strategies are described in Malik et al. (2020b) as a microbial modification to Grime’s competitor-stress tolerator-ruderal framework. In this newly proposed “YAS” framework (for yield-acquisition-stress tolerance), they hypothesize that microbes trade off resource investments among high yield (microbial biomass per unit of resource consumed), resource acquisition, or stress tolerance traits. Fungi with high yield strategies favor investment in central metabolism and assimilatory pathways necessary for building cellular components (Malik et al., 2020b). In contrast, fungi that allocate energy toward resource acquisition produce extracellular enzymes to break down various compounds in their environment. Investment in stress tolerance might involve increased production of osmolytes or polysaccharides to protect against desiccation (Schimel et al., 2007) and more effort spent on maintenance and repair of cellular structures (Kultz, 2005; Fuchs and Mylonakis, 2009).

These trade-offs are predicted in theoretical models (Pfeiffer et al., 2001; Allison, 2014; Manzoni et al., 2014). Empirical evidence for these trade-offs is strong for some, but limited for others. Several studies document trade-offs between growth yield and stress tolerance (Sterne and McCarver, 1978; Killham and Firestone, 1984; Dijksterhuis and de Vries, 2006; Gasch, 2007; Tiemann and Billings, 2011; Crowther et al., 2014; Garcia et al., 2020; Malik et al., 2020a). In contrast, trade-offs between resource acquisition and the other traits are less commonly studied. However, there is some support for trade-offs between stress tolerance and extracellular enzyme production, particularly for yeast (Treseder and Lennon, 2015; Morrison et al., 2018; Romero-Olivares et al., 2019). Soil microbial communities also exhibit trade-offs between C use efficiency, which could be a proxy for growth yield, and extracellular enzyme activity, especially for C-associated enzyme activity (Malik et al., 2019). Similarly, there are a few examples of trade-offs between growth rate and resource acquisition for fungal and bacterial species (Rinkes et al., 2011; Ramin and Allison, 2019; Whitney, 2019). Connecting trade-offs among these traits may be useful in predicting fungal response to climate change in the southwestern United States.

We hypothesized that fungal taxa from Southern California grassland preferentially invest in either growth yield, stress tolerance, or resource acquisition at the expense of the other traits, resulting in negative relationships between traits. To test this hypothesis, we measured growth yield (as fungal hyphal length per unit litter mass loss), drought stress tolerance (growth in drought vs. higher moisture levels), and resource acquisition (potential activity of four extracellular enzymes) traits from 15 decomposer fungi. Specifically, we investigated trade-offs between these traits by comparing their relationships across taxa.

Materials and Methods

Experimental Set-Up

Microcosms were established with grassland litter and fungal isolates. We collected dry, standing litter from Loma Ridge National Landmark in Southern California (33° 44' 13.2'' N, 117° 42' 42.0'' W, 365 m elevation), which is located on the traditional territory of the Acjachemen and Kizh communities (KIZH Nation, n.d.; Haas, 1995). This Mediterranean grassland is dominated by exotic and native annual grasses and forbs, including Avena, Bromus, Lolium, Erodium, and Lupinus, and native perennial grass Nassella pulchra (Potts et al., 2012; Martiny et al., 2016). The study site has an annual mean temperature of 17°C and mean precipitation of 30 cm (Kimball et al., 2014).

We roughly ground the litter (~1–2 cm) with coffee grinders and placed 5 g into 120 ml, amber widemouthed jars (Thermo Scientific; Supplementary Figure S1), then autoclaved them for 90 min at 121°C. Next, we inoculated the jars with one of 15 fungal strains that were previously isolated from Loma Ridge in 2017. These fungal strains represented 12 different species (some of the strains were the same species) from both Ascomycota and Basidiomycota phyla (Supplementary Table S1). The strains were identified via Sanger Sequencing using the ITS1F/ITS4 primer sets (White et al., 1990; Gardes and Bruns, 1993) and query coverage was at least 96% for all isolates according to the BLAST results. Isolates were stored as plugs in sterile water at room temperature for approximately 6 months (Jones et al., 1991) before regrowing on potato dextrose agar (with ampicillin and gentamicin) for the experiment. We transferred the regrown fungi to potato dextrose broth (with ampicillin and gentamicin) with continuous shaking. After 1 week, we centrifuged and rinsed the hyphae. The hyphal pellet was broken up and diluted with sterile tap water (optical density of 0.2 ± 0.02). We added 0.5 ml of dilute fungal hyphae to each jar of sterilized litter.

After inoculation, we added sterile deionized water to each microcosm so that the percent water holding capacity (WHC) equaled 4, 27, or 50%, and then placed the jars in an incubator set at 20°C. As a reference, we added ~0.05 ml to reach 4% WHC and ~17 ml of water to reach 50% WHC to the 5 g of litter. These moisture levels spanned conditions in the field site (Parolari et al., 2015) and aimed to capture the limits of water stress in this system. Additionally, from a practical sense, this was the minimum dryness level that was possible after the addition of the dilute hyphae. Each isolate (plus an additional uninoculated control) × moisture combination was replicated twice, for a total of 96 microcosms (15 isolates and 1 control × 3 moisture levels × 2 replicates = 96 microcosms). We incubated the microcosms (sterile litter, fungi, and water) for 5 weeks, airing the jars out weekly in a laminar flow hood for approximately 1 min to allow for oxygenation, while minimizing water loss and contamination risk. After the 5 week incubation period, fungal hyphae had spread throughout the microcosm litter. At that time, we froze a subset of the litter at −20°C for fungal hyphal biomass measurements and the rest at −80°C for extracellular enzyme activity measurements.

Trait Measurements: Growth Yield

We measured fungal hyphal length and mass loss to represent the growth yield lifestyle “Y” lifestyle, sensu (Malik et al., 2020b). We calculated fungal hyphal length by staining and microscopy, using a procedure modified from Allison et al. (2013). In brief, 0.5 g of frozen litter was stirred in a 39.5 g/L sodium hexametaphosphate solution to extract the fungal hyphae. Subsamples (5 ml) of this solution were vacuum-pumped through a 0.2 μm nylon filter and stained with acid fuchsin. This was repeated twice for each sample. We mounted the filters on slides and dried them overnight at 60°C in a drying oven. Using an Axioplan 2 imaging microscope, we took five photos at random of each filter (10 photos total per sample) to be representative of the whole filter paper. We measured fungal hyphal length using AxioVision and calculated fungal hyphal length per gram litter using the method described in Shen et al. (2016). Mass loss was determined by calculating the difference in litter mass at the beginning and end of the 5 week incubation. To estimate growth yield, we divided fungal hyphal length (m per microcosm) by total litter mass loss from the microcosm.

Trait Measurements: Resource Acquisition

We measured the potential extracellular enzyme activity of four enzymes involved in decomposition to assess resource acquisition traits (“A” lifestyle, sensu, Malik et al., 2020b). The enzymes measured included cellobiohydrolase (CBH) and β-glucosidase (BG) to estimate potential for cellulose degradation, β-xylosidase (BX) to estimate potential for hemicellulose degradation, and N-acetyl-β-D-glucosaminidase (NAG) to estimate potential for chitin degradation. We chose these enzymes because they are highly relevant in degrading plant fibers (Sinsabaugh et al., 2002) and fungal cell walls (Wohl and McArthur, 2001). We prepared sample homogenates and conducted fluorometric enzyme assays following methods described in Alster et al. (2013). In brief, 0.2 g of litter was added to 75 ml of 25 mM maleate buffer (pH 6.0) and homogenized using a Polytron automated homogenizer. About 200 μl of homogenate was added per well to 96-well plates with eight replicates per sample per assay. Each well also contained 50 μl of a fluorescent substrate. Controls of the sample homogenate, substrate solution, blanks, and standards were also replicated in eight wells. After 1 h, 10 μl of 1.0 M NaOH solution was added to each well to stop the reaction. Fluorescence was measured at 365 nm excitation and 450 nm emission. Enzyme activity was calculated according to Equation 1 in Alster et al. (2013) using the average values from the samples and controls. To control for enzyme activity being a function of fungal growth, we standardized potential extracellular enzyme activities per microcosm by multiplying by total litter mass remaining at the end of the experiment.

Trait Measurements: Stress Tolerance

To quantify drought stress tolerance of the isolates (“S” lifestyle, sensu, Malik et al., 2020b), we characterized their response to the three moisture levels. We reasoned that if an isolate produced more hyphae under the driest treatment, compared with the wettest treatment, then that isolate was relatively stress tolerant. Conversely, if an isolate grew best at the wettest treatment, but declined under drier conditions, then that isolate was less stress tolerant. To capture these patterns, we calculated a moisture association index (MAI) that considered total fungal hyphal length in each microcosm at each moisture level:

MAI=Hyphal length4%×0.04+Hyphal length27%×0.27+Hyphal length50%×0.50Hyphal length4%+Hyphal length27%+Hyphal length50%×100%

where Hyphal length4% was the hyphal length of the isolate in the 4% moisture treatment; Hyphal length27% was the hyphal length in the 27% moisture treatment; and so on. In theory, MAI could range between 4 (if the isolate only grew at 4% moisture) and 50% (if the isolate only grew at 50% moisture). The larger an isolate’s MAI, the less drought stress tolerant it was. We used this indicator because it weights the effect of moisture rather than treating growth at all moisture levels equally. This therefore “penalizes” fungi that have a much higher growth rate at higher moisture levels. While it is possible that this index could overestimate stress tolerance if growth peaked at the intermediate moisture level, or underestimate stress tolerance if growth was the same at all moisture levels, we generally did not find this to be the case for our isolates (Supplementary Figure S2).

Statistical Analysis

First, we tested whether isolates differed in the traits related to the three lifestyles: high growth yield, resource acquisition, and stress tolerance. To determine whether isolates possessed different degrees of stress tolerance (S lifestyle), we asked whether they varied in their responses to moisture. We performed a fully factorial ANOVA, with growth yield as the dependent variable, and moisture treatment and isolate as independent variables. A significant interaction between moisture treatment and isolate would indicate that isolates differed in responses to moisture. In addition, a significant main effect of isolate on biomass would support differences among isolates in the Y lifestyle.

Next, we conducted a series of fully-factorial ANOVAs to assess whether isolates varied in extracellular enzyme activity (A lifestyle). Dependent variables were activities of CBH, BG, BX, or NAG. Independent variables were isolate and moisture level.

To test our hypothesis that fungal isolates would exhibit trade-offs, we checked for negative relationships between the Y, A, and S traits. Within each moisture treatment there were no non-linear relationships, so we used a linear model to test relationships between traits. Specifically, we conducted a series of hierarchical linear mixed regressions between pairs of traits, with moisture as a covariate. Because isolates that were more closely related to one another may have tended to display similar traits, we included phylum, class, order, and family as nested random factors. This approach is similar to applying a phylogenetic independent contrast (Ricklefs and Starck, 1996), while allowing us to account for the moisture covariate. Significant negative relationships between pairs of traits would support the hypothesis. Growth yield was ranked before this analysis, because those data were not normally distributed.

We also checked for relationships between each of the four extracellular enzyme activities and looked at how taxonomic variation affected each of the traits. We conducted hierarchical linear mixed regressions between pairs of enzymes, with moisture as a covariate and phylum, class, order, and family as nested random factors. We also ran hierarchical linear models, with each trait as the dependent variable; and phylum, class, order, and family as nested factors to determine the influence of taxonomy. In all cases, statistical tests were conducted with Systat 13 (SPSS, 2017).

Finally, as an additional measure to assess trade-offs among the traits simultaneously, we conducted a redundancy analysis (RDA) in R version 3.5.3 (R Core Team, 2019). We used the “vegan” package (Oksanen, 2020) to conduct the RDA with moisture and fungal species as factors and phylum, order, and family as a random hierarchy, and growth yield, hyphal length, mass loss, and the enzyme activities as the dependent variables. We conducted a permutation test, similar to an ANOVA, for the RDA to determine significance.

Results

Variation Among Fungal Isolates in Traits and Moisture Responses

Fungal isolates differed significantly in their responses to moisture, based on their growth yield at low vs. medium and high moisture levels (Table 1, moisture × isolate interaction: p = 0.032). Overall, average growth yield from all of the isolates peaked at the lowest moisture level (Figure 1). Isolates also varied in total hyphal length. Some isolates, like Coprinellus aff xanthothrix Isolate 52, appeared most abundant in the driest treatment; others, like Fusarium poae, appeared to peak in the wettest treatment (Supplementary Figure S2). Across isolates, the MAI ranged from 18 to 36%.

TABLE 1
www.frontiersin.org

Table 1. Statistical results from fully-factorial ANOVA.

FIGURE 1
www.frontiersin.org

Figure 1. Growth yield for each isolate and their relationship to incubation moisture level. Each line represents one isolate. Symbols are means ± 1SE of two replicates. Bars are averages of all isolates.

In contrast, isolates did not vary significantly in the degree to which extracellular enzyme activity was sensitive to moisture (Table 1, moisture × isolate interaction: p > 0.176 for all enzymes). Although, across all moisture levels, isolates differed significantly from one another in CBH activity and BG activity (Supplementary Figure S3). In addition, across all isolates, growth yield, CBH activity, and NAG activity differed among moisture regimes. Specifically, CBH activity peaked at intermediate moisture, while NAG activity increased linearly with moisture.

Several traits also varied with taxonomy. MAI varied significantly at the phylum and family levels (Supplementary Table S2, p < 0.01). BG activity varied significantly at the family level (Supplementary Figure S2, p = 0.015), while NAG activity varied significantly at the class level (Supplementary Table S2, p = 0.047). Growth yield, CBH activity, and BX activity did not vary with taxonomic rank.

Relationships Between Moisture Responses, Growth Yield, and Enzymes

In no case were MAI or growth yield significantly related to any other trait we measured (Table 2; Figure 2). Thus, we rejected our hypothesis – these traits were not negatively related to one another. We found some support for a trade-off between growth yield and the extracellular enzyme activities based on the RDA (Supplementary Figure S5), however, this analysis was not significant (p = 0.561). Additionally, we also examined relationships between extracellular enzyme activities. CBH activity was significantly and positively related to BG and NAG (Table 2; Figure 3, p < 0.001, R = 0.875 for BG and R = 0.507 for NAG). β-glucosidase and NAG were also significantly positively related to one another (Figure 3; p < 0.001, R = 0.474).

TABLE 2
www.frontiersin.org

Table 2. Statistical results examining pairwise relationships between fungal traits.

FIGURE 2
www.frontiersin.org

Figure 2. Correlations between moisture association index (MAI), growth yield, and extracellular enzyme activities of fungal isolates. Each point represents one fungal isolate (mean of two sampling replicates). Colors represent moisture level. p-values are for relationships between traits across all samples, with moisture as a covariate and phylum, class, order, family, and genus as nested random factors. Correlations between the ranked growth yield data and the other traits can be found in Supplementary Figure S4.

FIGURE 3
www.frontiersin.org

Figure 3. Relationships between extracellular enzyme activities of fungal isolates. Each point represents one fungal isolate (mean of two sampling replicates). Colors represent moisture level. p-values are for relationships between traits across all samples, with moisture as a covariate and phylum, class, order, family, and genus as nested random factors. Lines are best fit for significant relationships.

Discussion

In this study, we investigated if fungi from a Mediterranean grassland sorted according to the YAS framework (Malik et al., 2020b), in order to improve predictions of fungal response to climate change in southern California. We found minimal evidence for evolutionary or physiological trade-offs between traits associated with growth yield, resource acquisition, and stress tolerance. Specifically, there were no significant negative relationships among growth yield, drought stress tolerance, or activities of four extracellular enzymes that target cellulose-C, hemicellulose-C, and chitin-N (Figure 2). Perhaps investment in stress tolerance traits or extracellular enzyme production did not come at a noticeable cost to growth yield or to one another. Alternately, traits related to the YAS lifestyles may have all been important for maintaining fungal fitness in this ecosystem, so that fungi may prioritize investment in all of these traits at the cost of other traits not measured here. No matter the underlying mechanism, these results suggest that greater drought stress under climate change will not necessarily select for or against fungi with high growth yield or strong resource acquisition at this site.

For comparison, previous empirical studies at this field site have reported negative relationships between growth yield and stress tolerance, albeit in the whole microbial communities and bacterial isolates (Ramin and Allison, 2019; Malik et al., 2020a) as oppose to fungal isolates. In particular, the abundance of growth yield metabolites (ectoine and 5-oxo-proline abundances) are negatively related to stress tolerance metabolites (aspartic acid and adenosine) among whole microbial communities at this ecosystem (Malik et al., 2020a). One reason for this discrepancy may be the absence of competition from other microbes in our study (Anthony et al., 2020). If competition is driving fungal community dynamics, fungi may increase growth yield or resource acquisition to monopolize space and resources at the expense of stress tolerance (Ghoul and Mitri, 2016; Wood et al., 2018). Fungi could also produce toxins or other compounds to target competitors (Chao and Levin, 1981; Riley and Gordon, 1999) at the expense of growth yield or resource acquisition. It would be interesting to examine how these fungal isolates would fare with respect to trait trade-offs with the addition of competition. Yet, a lab-based study from this site also found a trade-off between fast growth rate and potential extracellular enzyme activity for bacterial isolates, even though isolates were grown in the absence of competition (Ramin and Allison, 2019). Thus, even if these fungal isolates did not appear to experience trade-offs between growth yield and resource acquisition, it seems that other members of the microbial community do.

Studies examining these relationships specifically using fungal isolates in other systems have also reported trade-offs between stress tolerance vs. growth, or between growth rate and resource acquisition (Gasch, 2007; Whitney, 2019). In a review of genomic expression in response to environmental stress for several strains of Ascomycota isolated from unspecified environments, stress resistance came at the expense of cellular growth (Gasch, 2007). Similarly, in a study of 10 saprotrophic fungal isolates from a temperate forest in Massachusetts, Whitney (2019) observed a negative relationship between resource acquisition and specific growth rate. These studies were conducted on liquid media or agar (instead of litter; Gasch, 2007; Whitney, 2019), which may partially explain our differing conclusions. However, perhaps trade-offs between these traits become significant under wetter or cooler conditions than are experienced by the fungi in our Southern Californian grassland (Wallenstein and Hall, 2012).

Fungal activity seemed limited at lower moisture levels because decomposition rates were slower (Supplementary Table S1). However, it remains possible that the moisture levels in our microcosms were not dry enough to cause sufficient drought stress despite large differences in litter water content (4–50% WHC). Nevertheless, there is evidence that these fungal strains are adapted to drought stress. We found that some isolates had their highest growth at the lowest moisture level (Supplementary Figure S2), and growth yield was largest in the lowest moisture level for nearly all of the fungal isolates (Figure 1). The higher growth yield at the lower moisture levels could be driven by greater hyphal growth to forage for resources (Jennings, 1987), compared to the higher moisture levels. Lower growth yield in the higher moisture levels could also be due to lower carbon use efficiency (i.e., more CO2 production; Manzoni et al., 2012b). Investigating these confounding factors may be warranted in future studies of fungal response to drought in this ecosystem.

In general, fungi tend to withstand drought better than bacteria (Yuste et al., 2011; De Vries et al., 2012; Manzoni et al., 2012a; Alster et al., 2013), due to their desiccation-resistant morphologies (Barnard et al., 2013; Treseder and Lennon, 2015) and more stable fungal networks (de Vries et al., 2018). In our study system, the ability of fungal isolates to maintain growth yield and extracellular enzyme production, even during drought, may allow fungi to maintain ecosystem-level decomposition rates where climate change increases drought severity. The knowledge gained in this study has implications for trait-based models of microbial decomposition, which often assume that trade-offs exist between growth yield and functions such as enzyme production (Allison, 2014; Wieder et al., 2014; Allison and Goulden, 2017). This assumption for soil fungi may require revisiting in this arid ecosystem. However, it is also possible that resource acquisition (e.g., by extracellular enzymes) may be used to support investment in drought tolerance or yield, especially in carbon-rich litter. For example, C gained from litter could be invested in C-rich polysaccharides to protect against desiccation. If this mechanism dominates, it could explain the lack of trade-offs observed in this study.

Production of extracellular cellulases (CBH and BG), hemicellulases (BX), and chitinases (NAG) were, in many cases, positively correlated with one another among fungal taxa (Figure 3). These results are supported by several other studies also noting positive relationships between potential activities of these extracellular enzymes (Keeler et al., 2009; Sinsabaugh et al., 2009; Waring et al., 2014; Talbot et al., 2015; Jian et al., 2016). These positive relationships may exist because multiple enzymes are needed to digest the available organic matter. For example, BG is needed to break down the products of CBH (Tsai et al., 2014). Given that these microcosms contained fungal isolates, instead of a community, the fungi would have had to produce both BG and CBH, as well as endocellulases, to hydrolyze cellulose into glucose (Tsai et al., 2014). The C-acquiring enzymes BG and CBH and N-acquiring NAG may also be produced synergistically in order to maintain stoichiometric ratios if C and N are co-limiting (Sinsabaugh et al., 2009; Waring et al., 2014).

In this study, we selected traits representative of fungal physiology with known consequences for soil decomposition dynamics (Treseder and Lennon, 2015; Malik et al., 2020b). Other traits, such as hyphal morphology, sporocarp production, spore size, or pathogenicity (Zanne et al., 2020), may elicit stronger trade-offs than those associated with the YAS lifestyles. Measurement of other extracellular enzymes, for example, leucine amino peptidase to capture peptide degradation (Talbot et al., 2015) or oxidative enzymes to capture lignin degradation (Alster et al., 2013), might also have elicited stronger trade-offs than what we observed here. A more targeted approach to measuring stress tolerance, for example, by measuring markers for synthesis of osmolyte production or biomolecular repair (Malik et al., 2020b), could also have yielded stronger negative relationships with the growth yield or resource acquisition traits. We acknowledge that this study examined a relatively small fungal collection with limited diversity, potentially biasing our results. Additionally, perhaps the ability to be cultured self-selects for a specific type of isolate, which could influence the observed trait responses. However, prevalence of Ascomycota among these strains is representative of the sequences dominating dried leaf litter from the Loma Ridge site (Matulich et al., 2015; Glassman et al., 2018). Regardless, it is difficult – if not impossible – with current technology to include here the full complement of fungal species in this ecosystem. Finally, we limit our interpretations to Southern California grasslands until trade-offs can be assessed in fungi from other dry ecosystems.

In conclusion, fungi isolated from this grassland exhibited no detectable trade-offs between growth yield, resource acquisition, and drought stress tolerance. Since the drought stress tolerance exhibited by a given fungal isolate was not significantly related to its ability to break down organic matter, we do not necessarily expect that shifts in fungal communities under drought would alter fungal decomposer ability. By including this information in trait-based ecosystem models, we may improve predictions of ecosystem function under climate change in Southern California.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

KT, SA, and AM developed the original concept and obtained the funding. CA designed and planned the experiment with supervision from KT. SG isolated the fungi. CA performed the experiment. CA and KT analyzed the data with input from SA. CA wrote the original manuscript draft with significant involvement from KT. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by Grants from NSF (DEB 1912525) and the Department of Energy Office of Biological and Environmental Research (DE-SC0016410).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank D. Nguyen, H. Samy, and M. Yang for their laboratory assistance, E. Morrison and L. A. Cat for their field help, and C. Weihe, J. Martiny, and T. Tsai for their support, particularly with the fungal collection. We also thank D. Banuelas, E. Esch, A. Ferreri, K. Gallego, H. Maughan, A. Kuhn, S. M. Kosina, T. R. Northen, and T. Swenson for providing comments on earlier versions of this manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.655987/full#supplementary-material

References

Allison, S. D. (2014). Modeling adaptation of carbon use efficiency in microbial communities. Front. Microbiol. 5:571. doi: 10.3389/fmicb.2014.00571

PubMed Abstract | CrossRef Full Text | Google Scholar

Allison, S. D., and Goulden, M. L. (2017). Consequences of drought tolerance traits for microbial decomposition in the DEMENT model. Soil Biol. Biochem. 107, 104–113. doi: 10.1016/j.soilbio.2017.01.001

CrossRef Full Text | Google Scholar

Allison, S. D., Lu, Y., Weihe, C., Goulden, M. L., Martiny, A. C., Treseder, K. K., et al. (2013). Microbial abundance and composition influence litter decomposition response to environmental change. Ecology 94, 714–725. doi: 10.1890/12-1243.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Alster, C. J., German, D. P., Lu, Y., and Allison, S. D. (2013). Microbial enzymatic responses to drought and to nitrogen addition in a southern California grassland. Soil Biol. Biochem. 64, 68–79. doi: 10.1016/j.soilbio.2013.03.034

CrossRef Full Text | Google Scholar

Anthony, M. A., Crowther, T. W., Maynard, D. S., van den Hoogen, J., and Averill, C. (2020). Distinct assembly processes and microbial communities constrain soil organic carbon formation. One Earth 2, 349–360. doi: 10.1016/j.oneear.2020.03.006

CrossRef Full Text | Google Scholar

Barnard, R. L., Osborne, C. A., and Firestone, M. K. (2013). Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 7, 2229–2241. doi: 10.1038/ismej.2013.104

PubMed Abstract | CrossRef Full Text | Google Scholar

Bennett, A. F., and Lenski, R. E. (2007). An experimental test of evolutionary trade-offs during temperature adaptation. Proc. Natl. Acad. Sci. U. S. A. 104, 8649–8654. doi: 10.1073/pnas.0702117104

PubMed Abstract | CrossRef Full Text | Google Scholar

Boddy, L., Büntgen, U., and Kauserud, H. (2016). Climate impacts on fungal community and trait dynamics. Fungal Ecol. 22, 17–25. doi: 10.1016/j.funeco.2016.03.005

CrossRef Full Text | Google Scholar

Chao, L., and Levin, B. R. (1981). Structured habitats and the evolution of anticompetitor toxins in bacteria. Proc. Natl. Acad. Sci. U. S. A. 78, 6324–6328. doi: 10.1073/pnas.78.10.6324

PubMed Abstract | CrossRef Full Text | Google Scholar

Crowther, T. W., Maynard, D. S., Crowther, T. R., Peccia, J., Smith, J. R., and Bradford, M. A. (2014). Untangling the fungal niche: the trait-based approach. Front. Microbiol. 5:579. doi: 10.3389/fmicb.2014.00579

PubMed Abstract | CrossRef Full Text | Google Scholar

de Vries, F. T., Griffiths, R. I., Bailey, M., Craig, H., Girlanda, M., Gweon, H. S., et al. (2018). Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9. doi: 10.1038/s41467-018-05516-7

PubMed Abstract | CrossRef Full Text | Google Scholar

De Vries, F. T., Liiri, M. E., Bjørnlund, L., Bowker, M. A., Christensen, S., Setälä, H. M., et al. (2012). Land use alters the resistance and resilience of soil food webs to drought. Nat. Clim. Chang. 2, 276–280. doi: 10.1038/nclimate1368

CrossRef Full Text | Google Scholar

Dijksterhuis, J., and de Vries, R. P. (2006). Compatible solutes and fungal development. Biochem. J. 399, e3–e5. doi: 10.1042/BJ20061229

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuchs, B. B., and Mylonakis, E. (2009). Our paths might cross: the role of the fungal cell wall integrity pathway in stress response and cross talk with other stress response pathways. Eukaryot. Cell 8, 1616–1625. doi: 10.1128/EC.00193-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, M. O., Templer, P. H., Sorensen, P. O., Sanders-demott, R., Groffman, P. M., and Bhatnagar, J. M. (2020). Soil microbes trade-off biogeochemical cycling for stress tolerance traits in response to year-round climate change. Front. Microbiol. 11:616. doi: 10.3389/fmicb.2020.00616

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardes, M., and Bruns, T. D. (1993). ITS primers with enhanced specificity for basidiomycetes‐application to the identification of mycorrhizae and rusts. Mol. Ecol. 2, 113–118. doi: 10.1111/j.1365-294x.1993.tb00005.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasch, A. P. (2007). Comparative genomics of the environmental stress response in ascomycete fungi. Yeast 24, 961–976. doi: 10.1002/yea.1512

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghoul, M., and Mitri, S. (2016). The ecology and evolution of microbial competition. Trends Microbiol. 24, 833–845. doi: 10.1016/j.tim.2016.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibson, P. B., Waliser, D. E., Guan, B., Deflorio, M. J., Ralph, F. M., and Swain, D. L. (2020). Ridging associated with drought across the western and southwestern United States: characteristics, trends, and predictability sources. J. Clim. 33, 2485–2508. doi: 10.1175/JCLI-D-19-0439.1

CrossRef Full Text | Google Scholar

Glassman, S. I., Weihe, C., Li, J., Albright, M. B. N., Looby, C. I., Martiny, A. C., et al. (2018). Decomposition responses to climate depend on microbial community composition. Proc. Natl. Acad. Sci. U. S. A. 115, 11994–11999. doi: 10.1073/pnas.1811269115

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, L. (1995). Conquests and historical identities in California, 1769–1936. Berkeley and Los Angeles, California: University of California Press.

Google Scholar

Ho, A., Kerckhof, F. M., Luke, C., Reim, A., Krause, S., Boon, N., et al. (2013). Conceptualizing functional traits and ecological characteristics of methane-oxidizing bacteria as life strategies. Environ. Microbiol. Rep. 5, 335–345. doi: 10.1111/j.1758-2229.2012.00370.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jennings, D. H. (1987). Translocation of solutes in fungi. Biol. Rev. 62, 215–243. doi: 10.1111/j.1469-185X.1987.tb00664.x

CrossRef Full Text | Google Scholar

Jian, S., Li, J., Chen, J., Wang, G., Mayes, M. A., Dzantor, K. E., et al. (2016). Soil extracellular enzyme activities, soil carbon and nitrogen storage under nitrogen fertilization: a meta-analysis. Soil Biol. Biochem. 101, 32–43. doi: 10.1016/j.soilbio.2016.07.003

CrossRef Full Text | Google Scholar

Jones, R., Sizmur, K., and Wildman, H. (1991). A miniaturised system for storage of fungal cultures in water. Top. Catal. 5, 184–186. doi: 10.1016/S0269-915X(09)80482-5

CrossRef Full Text | Google Scholar

Keeler, B., Hobbie, S., and Kellogg, L. (2009). Effects of long-term nitrogen addition on microbial enzyme activity in eight forested and grassland sites: implications for litter and soil organic matter decomposition. Ecosystems 12, 1–15. doi: 10.1007/s10021-008-9199-z

CrossRef Full Text | Google Scholar

Killham, K., and Firestone, M. K. (1984). Proline transport increases growth efficiency in salt-stressed Streptomyces griseus. Appl. Environ. Microbiol. 48, 239–241. doi: 10.1128/AEM.48.1.239-241.1984

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimball, S., Goulden, M. L., Suding, K. N., and Parker, S. (2014). Altered water and nitrogen input shifts succession in a southern California coastal sage community. Ecol. Appl. 24, 1390–1404. doi: 10.1890/13-1313.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kivlin, S. N., and Treseder, K. K. (2014). Soil extracellular enzyme activities correspond with abiotic factors more than fungal community composition. Biogeochemistry 117, 23–37. doi: 10.1007/s10533-013-9852-2

CrossRef Full Text | Google Scholar

KIZH Nation (n.d.). KIZH NATION (Pronounced Keech), Gabrieleño Band Of Mission Indians. Available at: https://gabrielenoindians.org/ (Accessed September 23, 2020).

Google Scholar

Kultz, D. (2005). Molecular and evolutionary basis of the cellular stress response. Annu. Rev. Physiol. 67, 225–257. doi: 10.1146/annurev.physiol.67.040403.103635

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, A. A., Brodie, E. L., Swenson, T., Northen, T. R., Weihe, C., Morrison, E. W., et al. (2020a). Drought and plant litter chemistry alter microbial gene expression and metabolite production. ISME J. 14, 2236–2247. doi: 10.1038/s41396-020-0683-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, A. A., Martiny, J. B. H., Brodie, E. L., Allison, S. D., and Martiny, A. C. (2020b). Defining trait-based microbial strategies with consequences for soil carbon cycling under climate change. ISME J. 14, 1–9. doi: 10.1038/s41396-019-0510-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, A. A., Puissant, J., Goodall, T., Allison, S. D., and Griffiths, R. I. (2019). Soil microbial communities with greater investment in resource acquisition have lower growth yield. Soil Biol. Biochem. 132, 36–39. doi: 10.1016/j.soilbio.2019.01.025

CrossRef Full Text | Google Scholar

Manzoni, S., Schaeffer, S. M., Katul, G., Porporato, A., and Schimel, J. P. (2014). A theoretical analysis of microbial eco-physiological and diffusion limitations to carbon cycling in drying soils. Soil Biol. Biochem. 73, 69–83. doi: 10.1016/j.soilbio.2014.02.008

CrossRef Full Text | Google Scholar

Manzoni, S., Schimel, J. P., and Porporato, A. (2012a). Responses of soil microbial communities to water stress: results from a meta-analysis. Ecology 93, 930–938. doi: 10.1890/11-0026.1

CrossRef Full Text | Google Scholar

Manzoni, S., Taylor, P., Richter, A., Porporato, A., and Ågren, G. I. (2012b). Environmental and stoichiometric controls on microbial carbon-use efficiency in soils. New Phytol. 196, 79–91. doi: 10.1111/j.1469-8137.2012.04225.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Martiny, J. B. H., Martiny, A. C., Weihe, C., Lu, Y., Berlemont, R., Brodie, E. L., et al. (2016). Microbial legacies alter decomposition in response to simulated global change. ISME J. 11, 490–499. doi: 10.1038/ismej.2016.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Matulich, K. L., Weihe, C., Allison, S. D., Amend, A. S., Berlemont, R., Goulden, M. L., et al. (2015). Temporal variation overshadows the response of leaf litter microbial communities to simulated global change. ISME J. 9, 2477–2489. doi: 10.1038/ismej.2015.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Maynard, D. S., Bradford, M. A., Covey, K. R., Lindner, D., Glaeser, J., Talbert, D. A., et al. (2019). Consistent trade-offs in fungal trait expression across broad spatial scales. Nat. Microbiol. 4, 846–853. doi: 10.1038/s41564-019-0361-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrison, E. W., Pringle, A., Diepen, L. T. A.Van, and Frey, S. D. (2018). Simulated nitrogen deposition favors stress-tolerant fungi with low potential for decomposition. Soil Biol. Biochem. 125, 75–85. doi: 10.1016/j.soilbio.2018.06.027

CrossRef Full Text | Google Scholar

Oksanen, J. (2020). Vegan: an introduction to ordination. vegan 2.5-7. Available at: http://doi.acm.org/10.1145/2037556.2037605%5Cnftp://ftp3.ie.freebsd.org/ pub/cran.r-project.org/web/packages/vegan/vignettes/intro-vegan.pdf (Accessed February 18, 2021).

Google Scholar

Parolari, A. J., Goulden, M. L., and Bras, R. L. (2015). Ecohydrological controls on grass and shrub above-ground net primary productivity in a seasonally dry climate. Ecohydrology 8, 1572–1583. doi: 10.1002/eco.1605

CrossRef Full Text | Google Scholar

Pfeiffer, T., Schuster, S., and Bonhoeffer, S. (2001). Cooperation and competition in the evolution of ATP-producing pathways. Science 292, 504–508. doi: 10.1126/science.1058079

PubMed Abstract | CrossRef Full Text | Google Scholar

Potts, D. L., Suding, K. N., Winston, G. C., Rocha, A. V., and Goulden, M. L. (2012). Ecological effects of experimental drought and prescribed fire in a southern California coastal grassland. J. Arid Environ. 81, 59–66. doi: 10.1016/j.jaridenv.2012.01.007

CrossRef Full Text | Google Scholar

Ramin, K. I., and Allison, S. D. (2019). Bacterial tradeoffs in growth rate and extracellular enzymes. Front. Microbiol. 10:2956. doi: 10.3389/fmicb.2019.02956

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2019). R: A Language and Environment for Statistical Computing. Available at: https://www.r-project.org/ (Accessed February 18, 2021).

Google Scholar

Ricklefs, R. E., and Starck, J. M. (1996). Applications of phylogenetically independent contrasts: a mixed progress report. Oikos 77, 167–172. doi: 10.2307/3545598

CrossRef Full Text | Google Scholar

Riley, M. A., and Gordon, D. M. (1999). The ecological role of bacteriocins in bacterial competition. Trends Microbiol. 7, 129–133. doi: 10.1016/S0966-842X(99)01459-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Rinkes, Z. L., Weintraub, M. N., Deforest, J. L., and Moorhead, D. L. (2011). Microbial substrate preference and community dynamics during decomposition of Acer saccharum. Fungal Ecol. 4, 396–407. doi: 10.1016/j.funeco.2011.01.004

CrossRef Full Text | Google Scholar

Romero-Olivares, A. L., Meléndrez-Carballo, G., and Lago-Lestón, A. (2019). Soil metatranscriptomes under long-term experimental warming and drying: fungi allocate resources to cell metabolic maintenance rather than decay. Front. Microbiol. 10:1914. doi: 10.3389/fmicb.2019.01914

PubMed Abstract | CrossRef Full Text | Google Scholar

Schimel, J., Balser, T. C., and Wallenstein, M. (2007). Microbial stress-response physiology and its implications for ecosystem function. Ecology 88, 1386–1394. doi: 10.1890/06-0219

PubMed Abstract | CrossRef Full Text | Google Scholar

Seager, R., Ting, M., Held, I., Kushnir, Y., Lu, J., Vecchi, G., et al. (2007). Model projections of an imminent transition to a more arid climate in southwestern North America. Science 316, 1181–1184. doi: 10.1126/science.1139601

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Q., Kirschbaum, M. U. F., Hedley, M. J., and Arbestain, M. C. (2016). Testing an alternative method for estimating the length of fungal hyphae using photomicrography and image processing. PLoS One 11, 1–12. doi: 10.1371/journal.pone.0157017.

CrossRef Full Text | Google Scholar

Sinsabaugh, R. L., Carreiro, M. M., and Alvarez, S. (2002). Enzyme and microbial dynamics of litter decomposition. Enzym. Environ. Act. Ecol. Appl. Marcel Dekker, New York, Basel, 249–265.

Google Scholar

Sinsabaugh, R. L., Hill, B. H., and Shah, J. J. F. (2009). Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature 462, 795–798. doi: 10.1038/nature08632

PubMed Abstract | CrossRef Full Text | Google Scholar

SPSS (2017). IBM SPSS Statistics for Windows, Version 25.0. Armonk, NY: IBM Corp.

Google Scholar

Sterne, R. E., and McCarver, T. H. (1978). Osmotic effects on radial growth rate and specific growth rate of three soil fungi. Can. J. Microbiol. 24, 1434–1437. doi: 10.1139/m78-231

PubMed Abstract | CrossRef Full Text | Google Scholar

Talbot, J. M., Martin, F., Kohler, A., Henrissat, B., and Peay, K. G. (2015). Functional guild classification predicts the enzymatic role of fungi in litter and soil biogeochemistry. Soil Biol. Biochem. 88, 441–456. doi: 10.1016/j.soilbio.2015.05.006

CrossRef Full Text | Google Scholar

Tiemann, L. K., and Billings, S. A. (2011). Changes in variability of soil moisture alter microbial community C and N resource use. Soil Biol. Biochem. 43, 1837–1847. doi: 10.1016/j.soilbio.2011.04.020

CrossRef Full Text | Google Scholar

Todd-Brown, K. E. O., Randerson, J. T., Hopkins, F., Arora, V., Hajima, T., Jones, C., et al. (2014). Changes in soil organic carbon storage predicted by earth system models during the 21st century. Biogeosciences 11, 2341–2356. doi: 10.5194/bg-11-2341-2014

CrossRef Full Text | Google Scholar

Treseder, K. K., and Lennon, J. T. (2015). Fungal traits that drive ecosystem dynamics on land. Microbiol. Mol. Biol. Rev. 79, 243–262. doi: 10.1128/MMBR.00001-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, C. T., Morales-Rodriguez, R., Sin, G., and Meyer, A. S. (2014). A dynamic model for cellulosic biomass hydrolysis: a comprehensive analysis and validation of hydrolysis and product inhibition mechanisms. Appl. Biochem. Biotechnol. 172, 2815–2837. doi: 10.1007/s12010-013-0717-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallenstein, M. D., and Hall, E. K. (2012). A trait-based framework for predicting when and where microbial adaptation to climate change will affect ecosystem functioning. Biogeochemistry 109, 35–47. doi: 10.1007/s10533-011-9641-8

CrossRef Full Text | Google Scholar

Waring, B. G., Weintraub, S. R., and Sinsabaugh, R. L. (2014). Ecoenzymatic stoichiometry of microbial nutrient acquisition in tropical soils. Biogeochemistry 117, 101–113. doi: 10.1007/s10533-013-9849-x

CrossRef Full Text | Google Scholar

White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics” in Pcr protocols: A guide to methods and applications. eds. M. A. Innis, D. H. Gelfand, J. J. Sninsky, and T. J. White (San Diego, California, USA; London, England, UK: Academic Press, Inc.), 315–322.

Google Scholar

Whitney, S. (2019). Fungal physiological response to temperature and nitrogen availability.

Google Scholar

Wieder, W. R., Bonan, G. B., and Allison, S. D. (2013). Global soil carbon projections are improved by modelling microbial processes. Nat. Clim. Chang. 3, 909–912. doi: 10.1038/nclimate1951

CrossRef Full Text | Google Scholar

Wieder, W. R., Grandy, A. S., Kallenbach, C. M., and Bonan, G. B. (2014). Integrating microbial physiology and physio-chemical principles in soils with the MIcrobial-MIneral carbon stabilization (MIMICS) model. Biogeosciences 11, 3899–3917. doi: 10.5194/bg-11-3899-2014

CrossRef Full Text | Google Scholar

Wohl, D. L., and McArthur, J. V. (2001). Aquatic actinomycete-fungal interactions and their effects on organic matter decomposition: a microcosm study. Microb. Ecol. 42, 446–457. doi: 10.1007/s00248-001-0005-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, J. L., Tang, C., and Franks, A. E. (2018). Competitive traits are more important than stress-tolerance traits in a cadmium-contaminated rhizosphere: a role for trait theory in microbial ecology. Front. Microbiol. 9:121. doi: 10.3389/fmicb.2018.00121

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, W., Chen, X., Jing, X., and Zhu, B. (2018). A meta-analysis of soil extracellular enzyme activities in response to global change. Soil Biol. Biochem. 123, 21–32. doi: 10.1016/j.soilbio.2018.05.001

CrossRef Full Text | Google Scholar

Yuste, J. C., Penuelas, J., Estiarte, M., Garcia-Mas, J., Mattana, S., Ogaya, R., et al. (2011). Drought-resistant fungi control soil organic matter decomposition and its response to temperature. Glob. Chang. Biol. 17, 1475–1486. doi: 10.1111/j.1365-2486.2010.02300.x

CrossRef Full Text | Google Scholar

Zanne, A. E., Abarenkov, K., Afkhami, M. E., Aguilar-Trigueros, C. A., Bates, S., Bhatnagar, J. M., et al. (2020). Fungal functional ecology: bringing a trait-based approach to plant-associated fungi. Biol. Rev. 95, 409–433. doi: 10.1111/brv.12570

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: fungal traits, extracellular enzymes, drought, YAS framework, litter decomposition

Citation: Alster CJ, Allison SD, Glassman SI, Martiny AC and Treseder KK (2021) Exploring Trait Trade-Offs for Fungal Decomposers in a Southern California Grassland. Front. Microbiol. 12:655987. doi: 10.3389/fmicb.2021.655987

Received: 19 January 2021; Accepted: 05 March 2021;
Published: 20 April 2021.

Edited by:

Manuel Delgado Baquerizo, University of Alicante, Spain

Reviewed by:

Carlos Aguilar-Trigueros, Freie Universität Berlin, Germany
Ashley Lang, Indiana University, United States

Copyright © 2021 Alster, Allison, Glassman, Martiny and Treseder. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Charlotte J. Alster, charlotte.alster@waikato.ac.nz

Present address: Charlotte J. Alster, School of Science, University of Waikato, Hamilton, New Zealand

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.