Tree seedling shade tolerance arises from interactions with microbes and is mediated by functional traits

1 Introduction

Tree seedling mortality responses to understory light availability are an important filter of mature tree species composition and drivers of forest community dynamics (Pacala et al., 1996). Thus, plant survival in low light, or shade tolerance (Shirley, 1943) is a central concept in forest ecology. The seedling establishment phase, a major demographic bottleneck, is also critical for future community dynamics (Gurevitch et al., 2020). However, seedling responses to shade are far more complex than responses to light availability alone (Valladares and Niinemets, 2008; Valladares et al., 2016). Understanding the functional traits and biotic and abiotic conditions that convey shade tolerance are key to a more mechanistic understanding of forest community dynamics.

A plant’s ability to tolerate low light conditions can be moderated by soil-borne microbes, like soil-borne pathogens and mycorrhizal fungi (Jiang et al., 2020). Soil-borne pathogens (including fungi, oomycetes, and bacteria) can cause high seedling mortality (Song and Corlett, 2022), especially in shade (O’Hanlon-Manners and Kotanen, 2004; McCarthy-Neumann and Kobe, 2008; McCarthy-Neumann and Ibáñez, 2013), where wetter conditions enhance microbe reproduction and dispersal (Reinhart et al., 2010; Hersh et al., 2012; Liu and He, 2019). Arbuscular mycorrhizal fungi (AMF) can provide water and nutrients in exchange for sugars (Wipf et al., 2019). However, they may parasitize tree seedlings and increase mortality in low light (Ibáñez and McCarthy-Neumann, 2016; Konvalinková and Jansa, 2016), despite mutualistic tendencies and greater abundance in high light (Bereau et al., 2000; Shi et al., 2014; Koorem et al., 2017).

AMF, pathogens, and light levels can interact to influence seedling mortality. In high light, AMF root colonization is greater (Ibáñez and McCarthy-Neumann, 2016; Konvalinková and Jansa, 2016; Koorem et al., 2017), which can reduce the growth of fungal pathogens, potentially by competing for root space (Borowicz, 2001). AMF may also indirectly ameliorate pathogen effects (Liang et al., 2015), by providing water and nutrients to the host plant (Graham, 2001) and inducing production of defensive traits (Pozo and Azcón-Aguilar, 2007; Zamioudis and Pieterse, 2012) that protect against pathogens (Azcón-Aguilar et al., 2002; Violle et al., 2012). Conversely, in low light, seedling mortality may increase, due to the combined carbon costs of maintaining the AMF mutualism and recovery from pathogen attack.

Seedlings could also experience higher mortality from soil-borne microbes due to shade-induced changes in defensive functional trait values, such as reduced phenolics (Ichihara and Yamaji, 2009) and lignin (Rogers et al., 2005; Falcioni et al., 2018). Additionally, carbon limitation in shade and lower stored nonstructural carbohydrates (NSC) may constrain recovery from disease (Kobe, 1997; Kobe et al., 2010). Functional trait values may not only differ within a species based on light availability, but may also vary among species in relation to shade tolerance (Imaji and Seiwa, 2010). Shade tolerant species are typically less vulnerable to mortality by soil-borne microbes than shade intolerant species (McCarthy-Neumann and Kobe, 2008; McCarthy-Neumann and Kobe, 2010; Cortois et al., 2016) at least partly because shade tolerant species allocate more carbon to chemical and physical defenses (Coley et al., 1985; Coley and Barone, 1996), and recovery (Kitajima, 1994; Myers and Kitajima, 2007; Poorter et al., 2010).

To examine the effects of light availability, soil-borne microbes, tree seedling functional traits, and their interactions on light-dependent seedling survival, we established an experiment to test the following hypotheses: 1) within species, decreased survival under low versus high light only occurs in the presence of soil microbes; 2) mycorrhizal colonization is lower in the <20 µm microbial filtrate where pathogens are likely to be the dominant microbial group; 3) as defensive traits, phenolics and lignin are induced to higher levels in soils where microbes are present (non-sterilized) and in higher light availability; 4) as a recovery trait, NSC is lower in soils where microbes are present and in low light availability; 5) across all species and when microbes are present, survival increases as phenolics, lignin, and NSC increase.

2 Materials and methods

We conducted a fully factorial blocked-design greenhouse experiment at the Michigan State University Tree Research Center in Lansing, Michigan, USA (42.7 °N, 84.5 °W). The experiment consisted of three species (Acer negundo, A. rubrum, and A. saccharum), four microbial communities (<20 µm, representing pathogenic microbes; 40-250 µm, representing AMF; combined filtrate [both <20 µm and 40-250 µm]; and sterilized combined filtrate) and two light levels (2% and 30% full sun, representing shade and light gap environments). Individual pots were set up on six different benches (three per light level), where all treatment combinations were represented. We planted 80 seedlings per treatment combination for a total of 1,920 seedlings. We monitored seedlings every three days for survival, and randomly selected subsets for trait measurements at three, six, and nine weeks.

2.1 Species selection

We selected three biogeographically widespread, co-occurring tree species within the genus Acer: saccharum, rubrum, and negundo. These species have similar seed sizes (Osunkoya et al., 1994), but vary in shade tolerance (Burns and Honkala, 1990; Niinemets and Valladares, 2006).

2.2 Light availability

We grew seedlings at two light levels (2% and 30% full sun). We created light treatments by covering six greenhouse benches (three per treatment) with an inner layer of black shade cloth and an outer layer of reflective knitted poly-aluminum shade cloth (BFG Supply, Burton, Ohio, USA). We confirmed light levels using PAR (photosynthetically active radiation) measurements at each bench with a LI-COR 250A quantum sensor (LI-COR, Lincoln, Nebraska, USA) on a uniformly overcast day.

2.3 Soil collection and preparation of soil inocula

Soils were collected from Alma College’s Ecological Field Station in Vestaburg, Michigan, USA (43.4°N, 84.9 °W), in a 100-ha mixed-hardwood forest stand containing a 3-ha subplot with mapped and tagged trees. In August 2016, we randomly selected three adult trees per species. We selected adults that were at least two crown diameters away from other study species to reduce potential cross-culturing of soil. We collected soil (top 15 cm) within 1 m of each focal tree stem, maintained as separate replicates throughout the remainder of the experiment (as recommended by Rinella and Reinhart, 2018). We prepared soil by dicing roots and sifting soil through a 1 cm mesh sieve, retaining all roots that passed through the 1 cm sieve, as they may harbor host-specific microbial communities. Soil samples were stored at 4 °C for up to two months before preparation of soil inocula filtrates.

We created four microbial communities from sifted field soil using a wet-sieving method (Klironomos, 2002; Callaway et al., 2011; Pizano et al., 2014; Liang et al., 2015; König et al., 2016). For each extraction, we agitated 50 g of soil in a blender with 250 mL of deionized water at high speed for 60 sec, then passed the slurry through three analytical sieves (250-, 40-, and 20 µm) using a high-pressure water hose, for a total volume of 800 mL. To minimize contamination between treatments, we cleaned the sieves ultrasonically for 5 min between each extraction. The 250 µm sieve collected larger roots and coarser soil. We floated material retained by the 40 µm sieve on the surface of a 60% sucrose solution and centrifuged it at 688 g for 20 min. We collected material in the water and at the water-sucrose interface on 47 mm Whatman no.1 filter paper, surface sterilized it with 10% NaOCl for 10 sec, and washed it with distilled water under a filtration vacuum. We divided each filter paper into eight equal pieces. We collected the filtrate that passed through the 20 µm sieve and separated it into eight 100 mL containers. In sum, we created three microbial communities based on filtrate size classes: <20 µm, 40-250 µm, and combined (containing both <20 µm and 40-250 µm filtrates).

To test for abiotic effects of conspecific cultured soils due to nutrients or allelopathy, we combined filtrates (20- and 40-250 µm) and sterilized by steam autoclaves (at 121 °C for two hours). We also quantified and compared AMF colonization in the sterilized and <20 µm filtrates against a control containing only filter paper and deionized water to test the effectiveness of the sterilization for reducing microbes. We found no AMF colonization in soils treated with the sterilized (t = 2.97, df =176, p = 0.003) or <20 µm filtrates (t = 3.79, df = 178, p < 0.001), relative to a distilled water control. Filtrates were kept refrigerated at 4 °C for up to 48 hours, before adding them to the greenhouse pots.

Figure 1 Difference in predicted survival (mean ± 90% credible interval) in (A) low versus high light and (B) relative to sterilized soil. For each species and microbial communities, at the end of 9-weeks. Negative values indicate decreased survival in treatment versus (A) high light or (B) sterilized soil. Statistically significant differences (90% CI do not overlap 0) are indicated by asterisks.

In the greenhouse, we filled 1,920 (655 cm³) deepots (Stuewe and Sons, Tangent, Oregon, USA) with sterilized commercial topsoil (Hammond Farms Landscape Supply, Lansing, Michigan, USA). To aid seedling germination, we also topped each pot with 2 cm of 85% sterilized commercial soil mix, containing peat moss, perlite, and vermiculite (Fafard 4P Mix, Sun Gro Horticulture, Agawan, Massachusetts, USA). In pilot trials of this experiment, we found over 50% seedling mortality in the first two weeks, without the addition of the commercial soil mix (Wood pers. observation). We steam sterilized the topsoil and soil mix by autoclaving twice at 121 °C for 2 hours, with a 48-hour incubation period between cycles.

Within 48 hours of wet sieving, we added microbe treatments to soil, keeping filtrate from each adult and species separate. To enable microbial communities to sporulate and AMF hyphae to establish, we cultured the soil with Allium as bait plants, before planting Acer seedlings (Klironomos et al., 1999; Al-Yahya’ei et al., 2011). In January 2017, each pot was planted with three germinating Allium tuberosum seedlings. After two months, we removed aboveground A. tuberosum seedling biomass.

One week after Allium removal, we planted Acer seedlings whose hypocotyls had emerged within the three previous days. We purchased seeds from Sheffield’s Seed Company (Locke, New York, USA). To minimize microbes from non-experimental soil sources, we surface sterilized the seeds with 0.6% NaOCl both prior to cold stratification and germinating in perlite. We watered the seedlings three times per week with 50 mL deionized water.

2.4 Survival and functional trait measurements

We grew the Acer seedlings for nine weeks. We recorded emergence and survival every three days and assigned date of death as the first census with total leaf and stem tissue necrosis. We harvested a random subset of 20 seedlings for each species and each treatment at three, six, and nine weeks, to quantify AMF colonization, phenolics, lignin, and NSC. We chose these times for our sub-harvests, because in a previous greenhouse experiment, mortality curves for tree seedlings subjected to soil-borne pathogens often increased at week three and peaked between four to six weeks after germination (McCarthy-Neumann and Ibáñez, 2012).

To quantify percent AMF colonization, we stained seedling roots with a 5% Shaeffer black ink in vinegar solution (Vierheilig et al., 1998) and counted fungal structures (e.g., arbuscules, coils, vesicles, hyphae) along 100 intersections under the microscope (McGonigle et al., 1990). To quantify phenolics, we analyzed hypocotyl samples, using a microplate-adapted colorimetric total phenolics assay with Folin-Ciocalteu reagent (Ainsworth and Gillespie, 2007). To quantify lignin, we analyzed root and stem samples with an Ankom 200 fiber analyzer (ANKOM Technologies, Macedon, NY, USA), using the acid-fiber detergent fiber filter technique. To quantify NSC, we analyzed stem samples, using a standardized enzyme method for sugar and starch extraction and quantification (Landhäusser et al., 2018).

2.5 Statistical analyses

To test the influence of light availability and microbial community on tree seedling survival, we used an individual based counting process in a Cox survival model (Burnham and Anderson, 2002; McCarthy-Neumann and Ibáñez, 2012). Data for each seedling and each time, were coded as 0 until the seedling was found dead. We used a count process to model the number of events (mortality, ) until the experiment ended at nine weeks. We modeled the likelihood as:

and the process as:

where parameters were estimated as a function of the hazard (h), which is the intrinsic rate of mortality due to individual age or time within the experiment), and of risk (μ), which is the extrinsic rate of mortality due to light availability and microbial community. Simulations (3 chains) were run until convergence of the parameters was ensured (25,000 iterations) and then run for another 50,000 iterations, from which the posterior parameter values (Supplementary Figures 1, 2) and predicted survival (Supplementary Figure 3) were estimated. Predicted survival values were used to assess whether there were differences in how species responded to microbe treatments and light. We then used predicted survival values and their associated uncertainty to test if there were differences in how species responded to low versus high light and in different microbe treatments. Differences that did not include zero in their 90% credible intervals were considered statistically significant (Kruschke, 2014).

We used linear mixed effects models to evaluate the effects of light availability and microbial community on seedling traits. We ran individual models for each species and trait, where light level and microbial community were treated as fixed effects, and harvest time, bench (nested within light level), and adult tree were treated as random effects.

We then used linear mixed effects models to assess the effects of traits on seedling survival in the three weeks following trait measurement (i.e., traits collected at three weeks were used to predict seedling survival from three to six weeks, and traits at six weeks were used to predict seedling survival from six to nine weeks). We treated AMF colonization, phenolics, lignin, and NSC as fixed effects, and included light level as a covariate. Species was treated as a random effect.

We evaluated relative support for each of our linear mixed effects models using multi-modal inference with corrected Akaike’s Information Criterion (Burnham and Anderson, 2002). Models with ΔAICc < 6 of the best-approximating model were considered plausible (Richards, 2008).

We performed all analyses in R 3.5.1 (R Core Team, 2020). We used the rjags package (Plummer, 2019) to fit survival models and to run predicted survival and contrast simulations. We used the built-in “lmer” function to fit linear mixed effects models and tested significance of main effects using the “Anova” function in the car package (Fox and Weisberg, 2019). Model selection for linear mixed effects models was determined with the “step” function in the lmerTest package (Kuznetsova et al., 2017). Post-hoc Tukey pairwise comparisons of significant main effects were made using the “emmeans” and “joint_tests” functions in the multcomp package (Hothorn et al., 2008; Lenth, 2020).

3 Results

3.1 Differences in seedling survival in low versus high light occurred only when microbes were present

Differences in seedling survival in low light versus high light, a measure of shade tolerance, appeared only in the presence of soil biota (Figure 1A). In low light, survival decreased up to 14% with the <20 µm filtrate, 12% with the 40-250 µm filtrate, and 17% with the combined filtrate. A. negundo experienced the largest decreases in survival (17% reduction in the combined filtrate), whereas A. saccharum had the lowest (11% reduction in the combined filtrate). Furthermore, A. saccharum survival decreased only in the presence of combined filtrate.

Figure 2 (A) percent colonization AMF (%), (B) phenolic content (nmol Gallic acid equivalents per mg dry extract), (C) percent dry mass lignin (%), and (D) percent dry mass NSC (%). Means within each panel not sharing a letter are statistically different by the Tukey test at alpha = 0.05. Best-fit model terms are overlayed on each panel.

Microbial treatments influenced seedling survival in only three of the nine cases in high light, but in eight of nine in low light (Figure 1B). At low light, the presence of soil biota relative to sterilized soil decreased survival 12-26% for all species and filtrates, except for A. saccharum with 40-250 µm filtrate. Additionally, negative soil biota effects were more common with the <20 µm filtrate, decreasing survival for all three species, regardless of light level.

3.2 Seedling functional trait values varied with light availability and microbial community

Across all treatments, percent root colonization by AMF was 11-18% greater in high than low light for A. saccharum and A. negundo (Figure 2A). There was no AMF colonization in the <20 µm filtrate, and colonization was similar in the 40-250 µm and combined filtrates. Overall, A. negundo (76%) had the highest AMF colonization, compared to A. saccharum (66%) and A. rubrum (45%).

Phenolic content (nmol Gallic acid equivalents per mg dry extract) increased with light availability across almost all microbe filtrate treatments (Figure 2B); the only exception was A. saccharum with the <20 µm filtrate. For A. rubrum and A. negundo, phenolic content was negligible in low light with the sterilized and combined filtrates. Overall, A. saccharum had the highest phenolic content (0.23 nmol/mg), compared to A. negundo (0.14 nmol/mg) and A. rubrum (0.06 nmol/mg).

Percent dry mass lignin was greater in high than low light but depended upon species and microbial filtrate (Figure 2C). For A. saccharum, lignin increased 5% across microbe filtrates, in high versus low light. For A. negundo, lignin increased 27% but only in the 40-250 µm filtrate. Across species, A. saccharum had the highest percent dry mass lignin (12%), compared to A. rubrum (8%) and A. negundo (6%).

Percent dry mass NSC was greater in high than low light for all species and was generally greatest with the 40-250 µm filtrate (Figure 2D). Additionally, for A. rubrum and A. negundo, NSC decreased 40-70% in the <20 µm and combined microbe filtrates, in low versus high light. Across species, A. saccharum had the highest dry mass NSC (12%), compared to A. negundo (11%) and A. rubrum (8%).

3.3 Seedling functional trait values are associated with survival

Survival increased with both phenolics (χ² = 5.93, df = 1, p = 0.015; Figure 3A) and NSC (χ² = 7.72, df = 1, p = 0.005; Figure 3B). However, there was no significant relationship between light availability and survival (p > 0.05), either alone or interacting with phenolics and NSC.

Figure 3 Percent seedling survival as a function of (A) phenolic content (nmol Gallic acid equivalents per mg dry extract) and (B) percent dry mass NSC (%). Each point represents a mean of trait values at a harvest time (3 or 6 weeks) and survival for seedlings in the 3 weeks after harvest.

4 Discussion

Our results support that survival of newly germinated tree seedlings in low versus high light, or “shade tolerance,” may be due to interactions between low light and soil-borne microbes and be mediated by defense and recovery functional traits. Previous studies have demonstrated that seedling functional traits are influenced by light and microbes and that functional traits can influence growth and survival at low light (Falster et al., 2018), but have not linked resources, traits, and survival, as in this study.

Across species, overall seedling survivorship and insensitivity to shading corresponded with shade tolerance categorizations. A. saccharum (shade tolerant) had the highest overall survival and was least sensitive to the microbial filtrates, compared to A. rubrum (intermediate) and A. negundo (intolerant). In our study, A. negundo had the highest AMF colonization of the three study species (Figure 2A) and had the largest decreases in low-light survival with added microbial filtrates (Figure 1B). This aligns with prior studies which have shown that shade intolerant species experience greater mortality from disease in shade (Augspurger, 1984a; Pizano et al., 2014) and higher growth when grown in high light (Xi et al., 2023). Similarly, shade tolerant species showed no significant growth responses to microbial filtrate treatments, in contrast to pioneer species that were more sensitive to the habitat from which soil microbial filtrates were collected (Pizano et al., 2017).

4.1 Soil-borne microbes explain variation in tree seedling survival responses to light

There is good support that the 20 µm filtrate in our study is primarily composed of fungi and bacteria. We found no evidence of mycorrhizal colonization in our <20 µm soil filtrate (Figure 2A), suggesting that this treatment is mainly composed of non-mycorrhizal fungi and bacteria (also consistent with Klironomos, 2002). Additionally, there is support that the 40 µm filtrate is primarily composed of the mycorrhizal community associated with soils cultured by conspecific adults. In a prior study that investigated this methodology, Wagg et al. (2014) found that soil passing through 250 µm sieves contained ~80% of the mycorrhizal community with an additional ~20% of the mycorrhizal community passing through the 50 µm sieves. In the Wagg et al. (2014) study, filtrate ≤25 µm in size, effectively had no nematodes, <10% mycorrhizal fungi, ~70% other fungi and ~90% bacteria of the original soil community. In addition, several other studies have used these filtrate size classes to isolate and investigate the roles of soil-borne pathogens (McCarthy-Neumann and Kobe, 2008; König et al., 2016) and AMF (Klironomos, 2002; Callaway et al., 2011; Liang et al., 2015; Pizano et al., 2017).

Differences in low versus high light survival appeared only when soil-borne microbes were present (Figure 1A), supporting hypothesis 1. With microbes present, survival decreased for all species, with the largest differences occurring in low light (Figure 1B). Our results are consistent with previous research demonstrating that tree seedlings have higher mortality in shade and that the major cause of seedling death arises from disease (Vaartaja, 1962; Augspurger, 1984a; Augspurger, 1984b). Our results are also consistent with Liang et al. (2015), who utilized the wet-sieving method and found that pathogens were associated with decreased biomass and survival, AMF were associated with increased biomass, and combined filtrate treatments canceled each other out for both biomass and survival; however, Liang et al. (2015) did not consider interactions with light availability.

We did not see increases in seedling survival when AMF were present, despite high percent colonization of seedling roots. This is in contrast to the study by Liang et al. (2015), in which they found higher seedling survival with AMF. However, benefits of AMF in our study may have manifested in growth (Gehring, 2003), which we did not measure. Young seedlings still relying on maternal seed reserves and high resource availability in the greenhouse may have diminished the importance of AMF (Kulmatiski and Kardol, 2008; Forero et al., 2019; Heinze et al., 2020). Moreover, AMF may have acted indirectly by enhancing production of phenolics and NSC, which were both positively correlated with AMF colonization (Supplementary Figure 4). We also did not see any negative relationship between AMF colonization and NSC in low light, which could have suggested that AMF act parasitically when photosynthates are limited.

4.2 Amounts of functional traits varied with light availability and soil-borne microbes

For all seedling species examined, phenolic content increased when microbes were present, supporting this aspect of hypothesis 3. Both AMF (Vierheilig, 2004; Whipps, 2004; Pozo and Azcón-Aguilar, 2007) and pathogens (Nicholson and Hammerschmidt, 1992; Witzell and Martín, 2008) can induce phenolics production. Seedling phenolics also consistently increased with light availability (Figure 2B), supporting this aspect of hypothesis 3 and suggesting alleviation of photosynthate constraints on chemical defense production (Ballaré, 2014).

Similarly, percent dry mass lignin increased with light availability and was highest in the 40-250 µm filtrate (Figure 2C), also supporting hypothesis 3. Higher lignin in the 40-250 µm filtrate and lower lignin in the <20 µm and combined filtrates for A. rubrum suggest that AMF have a positive effect while pathogens have a negative effect on lignin production.

In partial support of hypothesis 4, we found that seedling NSC decreased when pathogens were present (Figure 2D), consistent with NSC reserves acting as a carbon buffer after damage (McPherson and Williams, 1998; Gleason and Ares, 2004; Myers and Kitajima, 2007; Kobe et al., 2010). However, contrary to predictions, NSC increased when AMF were present, but only when not combined with the pathogen filtrate. The positive association between AMF colonization and NSC is consistent with Li et al., (2022).

4.2 Functional traits were associated with greater seedling survival

Partly supporting hypothesis 5, we found that both phenolics (chemical defense) and NSC (carbon buffer precluding recovery from damage) had positive associations with tree seedling survival. These results are consistent with previous studies that have speculated higher allocation to defensive traits increases survival of shade tolerant seedlings in low light conditions (Vaartaja, 1962; Augspurger, 1984a; Augspurger and Kelly, 1984; Kitajima, 1994; Alvarez-Clare and Kitajima, 2007).

We found no evidence of an association between lignin and survival. This was in contrast to previous studies that have posited that differences in lignin development impact seedling susceptibility to pathogens (Sattler and Funnell-Harris, 2013; Lee et al., 2019; Zhu et al., 2021). In this study, A. saccharum, the most shade tolerant species, had up to 50% greater lignin than the other two species (Figure 2D), but these differences in lignin values did not manifest in differences in survival.

4.3 Caveats and future research

There are several areas upon which future research could build on this work. We used sucrose-centrifugation to separate AMF spores from most other microbes and debris, and a bleach sterilization step to kill potential pathogens. While we were able to see high AMF colonization in the 40-250 µm and combined filtrates, and no colonization in the <20 µm or sterilized filtrates (Figure 2A), we cannot eliminate the possibility that we excluded some AMF in smaller filtrate sizes and included additional microbes and debris in the 40-250 µm filtrate. A more robust method would include further isolating pure AMF spores, as with Callaway (2011) and Pizano et al. (2017). Alternatively, by adding a genetic analysis of the microbial inoculum, we could have determined more accurately which microbial groups were present, which would enhance our understanding of the results. We recommend that future research utilize the wet-sieving method in conjunction with spore isolation and/or genetic analyses.

Likewise, by culturing the pots with Allium, we may have inadvertently increased AMF presence, disproportionate to pathogens, or increased the relative abundance of microbes that specialize with Allium, rather than Acer species. Although often thought of as generalists, AMF can show some host specificity (Yang et al., 2012; Kajihara et al., 2022), which could influence post-culturing microbial communities. Thus, we recommend that future experiments culture AMF and other microbe communities with the host species of interest (in this case, Acer spp.). In future experiments, we recommend using the host species as a bait plant in the culturing step.

While we investigated three co-occurring species within a single genus to make broader generalizations, subsequent studies trying to generalize these results should include more species across additional levels of shade tolerance. Also, effects of light availability could be caused by changes in microclimate (e.g., soil temperature and moisture) and not directly due to irradiance. Similarly, effects of light may be mediated by photoreceptors and jasmonates, not just assimilate availability through higher photosynthetic rates (Ballaré and Austin, 2019; Pierik and Ballaré, 2021).

Additionally, we utilized shade cloth in the greenhouse to create shaded conditions, not vegetation shade. Future studies should consider teasing apart these mechanisms in the field, rather than the greenhouse, to provide more realistic seedling responses. Furthermore, although we were interested in light availability and shade tolerance, other environmental variables, such as nutrient or water availability, also could influence seedling survivorship (McCarthy-Neumann and Kobe, 2019).

4.4 Implications for forest community dynamics

This study provides a needed first step in developing a mechanistic understanding of how soil-borne microbes impact seedling shade tolerance, explained through functional traits. Although fast-growing shade intolerant species may be expected to outcompete shade tolerant species in high light (Pacala et al., 1996), shade intolerant species can be limited by the negative interactive effects of soil-borne microbes at low light (McCarthy-Neumann and Kobe, 2010; Liu and He, 2019), restricting their recruitment niche to areas with higher light and fewer soil-borne microbes. In this paper, we have demonstrated the importance of interactions between soil-borne microbes and light availability in determining tree seedling survival. Furthermore, we have related both intra- and interspecific differences in survival to functional traits, supporting a more trait-based and mechanistic approach to understanding forest community dynamics.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Dryad Repository https://doi.org/10.5061/dryad.7d7wm380b.

Author contributions

All authors conceived the ideas and methodology. KW collected and analyzed the data and lead the writing of the manuscript. RK and SM-N contributed critically to the drafts. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Science Foundation (NSF DEB 145732), Michigan State University, and Alma College.

Acknowledgments

We would like to thank Dr. Inés Ibáñez at University of Michigan for her help with developing the survival models and Dr. Walter Carson for his comments on earlier versions of the manuscript.

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.

Publisher’s note

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.

Supplementary material

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

References