ORIGINAL RESEARCH article

Front. Plant Sci., 31 August 2026

Sec. Functional Plant Ecology

Volume 17 - 2026 | https://doi.org/10.3389/fpls.2026.1899078

Root-fungal interactions upon resource foraging in heterogeneous soils

  • 1. Department of Botany, Faculty of Science, Charles University, Prague, Czechia

  • 2. Laboratory of Fungal Biology, Institute of Microbiology of the Czech Academy of Sciences, Prague, Czechia

Abstract

Introduction:

Plants respond to soil resource heterogeneity by foraging with their roots. The root foraging precision of plant species varies; and some species (especially those with low precision) may rely on symbiosis with mycorrhizal fungi to access nutrient-rich patches, thereby improving nutrient acquisition. However, the way herbaceous plants and arbuscular mycorrhizal fungi (AMF) interact during nutrient foraging is poorly understood.

Methods:

We tested nine mycorrhizal herbaceous plant species in a greenhouse experiment by exposing them to soil with heterogeneous or homogeneous nutrient distribution and with or without AMF inoculation. Root-foraging precision, shoot and root biomass production, and root-associated AMF gene copy number were quantified using plant biomass measurements and qPCR analyses of AMF.

Results:

Root foraging precision (in terms of roots placement) differed between species. Although foraging was evident in the heterogeneous nutrient treatment, it was largely unaffected by mycorrhizal inoculation. The exception was Melilotus officinalis which showed reduced root foraging precision with increasing root-associated AMF gene copy number. Plant biomass responses to AMF were species-specific, with overall root and shoot biomass declining slightly with increasing root-associated AMF gene copy number.

Discussion:

These results suggest that AMF inoculation generally does not alter root foraging precision across species, supporting the view that foraging precision is a species-specific trait independent of the collaboration gradient in the root economic spectrum. The balance between root foraging and mycorrhizal symbiosis may depend on other factors such as nutrient stoichiometry and forms. Future studies should thus examine root foraging, extraradical hyphal development and nutrient transfer under varying levels of nutrient limitation and spatial heterogeneity to better understand plants’ nutrient acquisition strategies.

1 Introduction

To acquire a sufficient amount of nutrients, plants often form various associations with other organisms, amongst the most important being mycorrhizal symbiosis. Over 90% of vascular plants establish some type of mycorrhizal association, most often with arbuscular mycorrhizal fungi (AMF) (Brundrett and Tedersoo, 2018). Plant roots and mycorrhizal hyphae cooperate to different degrees to acquire nutrients, forming a continuum of collaboration strategies that range from ‘do-it-yourself’ to ‘outsourcing’ approaches to nutrient uptake (Brundrett, 2009; Bergmann et al., 2020). The degree of cooperation depends on the characteristics of the two partners, as well as on environmental conditions (Johnson, 2010; Unger et al., 2016; Zhang et al., 2019; Šmilauer et al., 2020). In heterogeneous soil environments, plants may forage for nutrients using either their roots or mycorrhizal hyphae (Drew, 1975; John et al., 1983; Cavagnaro et al., 2005; Shi et al., 2011). Foraging precision appears to be a species-specific trait of both plants and fungi (Campbell et al., 1991; Cavagnaro et al., 2005; Kembel and Cahill, 2005; Weiser et al., 2016; Šmilauer et al., 2020) and may, therefore, be distributed along the collaboration continuum. This has been demonstrated for trees (Cheng et al., 2016) but remains still poorly understood for herbaceous species. Furthermore, despite the importance of mycorrhizal fungi for plant nutrient acquisition, root foraging experiments with herbaceous species are often conducted without controlling for mycorrhizal colonisation of plants and frequently in the absence of mycorrhiza altogether.

In a heterogeneous soil environment, plants’ association with AMF may help compensate for limited root foraging abilities (Cui and Caldwell, 1996b). Compared to roots, hyphae of AMF rapidly colonise soil, may better compete with other microorganisms for nutrients from nutrient-rich patches and can exploit smaller pores (Friese and Allen, 1991; Allen, 2007). The cost of producing new hyphae is lower than that of producing new roots (Tibbett, 2000; Hodge, 2006). The balance between the do-it-yourself strategy and outsourcing of nutrient acquisition to AMF (Bergmann et al., 2020) supports the existence of a root–hyphae foraging trade-off (but see Stiblíková et al., 2023). Indeed, a trade-off in foraging strategies between trees and AMF or ectomycorrhizal fungi has been reported (e.g. Eissenstat et al., 2015; Liu et al., 2015; Chen et al., 2016; Cheng et al., 2016). AMF-associated tree species tend to be more precise root foragers, whereas ectomycorrhizal fungi-associated tree species are less precise foragers and rely more heavily on fungal-mediated nutrient acquisition. However, this trade-off is poorly understood for herbaceous species.

Only a few studies have directly examined root foraging precision of herbaceous plants under the influence of mycorrhizal associations, with results showing either no change in root foraging after the addition of AMF (Farley and Fitter, 1999b; Hodge, 2001; Wijesinghe et al., 2001) or a reduction in root foraging precision in their presence (Cui and Caldwell, 1996a; Felderer et al., 2013). These contrasting results may reflect differences in mycorrhizal colonisation rates, which were often quite low or measured only as the presence of arbuscules in roots (Farley and Fitter, 1999b; Wijesinghe et al., 2001). Molecular quantification of AMF in roots and/or soil substrate, which is commonly available today, provides an additional approach for characterising plant–AMF associations (Felderer et al., 2013). Successful inoculation does not necessarily result in effective mycorrhizal colonization, as its establishment depends on multiple factors such as host-fungus compatibility and environmental conditions (e.g. nutrient and water availability [Mason et al., 2000; Klironomos, 2003; Lekberg and Koide, 2008; Frew, 2023]). Therefore, quantifying associated AMF gene copies using quantitative real-time PCR (qPCR) provides a more informative approach for evaluating mycorrhizal effects, allowing them to be analysed not only as a categorical factor (inoculated vs. non-inoculated) but also in relation to the actual amount of fungal DNA present in plant roots. Differences in the nature of nutrient patches used (organic, single-nutrient, or complex-nutrient patches) also complicate comparisons between experiments. In our previous study on root foraging and mycorrhizal symbiosis, no clear pattern of a trade-off along the collaboration gradient was observed, although the data were compiled from multiple databases rather than from a single experiment (Stiblíková et al., 2023). From the available evidence, the degree of plant cooperation with mycorrhizal fungi during root foraging for nutrients appears to be species-specific; however, additional evidence is needed to confirm this.

To further investigate root–fungal interactions during root foraging for heterogeneous nutrient sources, the present study linked the root allocation patterns of nine grassland species in heterogeneous nutrient environments with the presence of AMF in their roots. We examined how plant root foraging ability and biomass were affected by both the presence and the number of root-associated AMF gene copies in the central part of the root system. Plant species were grown in pots with either heterogeneous or homogeneous nutrient distributions and with or without mycorrhizal inoculum. Root dry biomass served as a proxy for plant foraging precision, and root-associated AMF nuclear large ribosomal subunit (LSU) copy numbers in the central part of the root system were quantified using qPCR. Our null hypothesis was that plant root foraging precision was unaffected by the presence of root-associated AMF LSU copy number. Our alternative hypothesis was that root foraging precision would be affected by the presence and number of root-associated AMF LSU copies, with higher root-associated AMF LSU copy number being associated with lower root foraging precision. Because we did not measure the distribution of extraradical hyphae, nutrient uptake by fungal hyphae, or nutrient transfer to plants from hyphae (e.g., using isotopes), our analyses were limited to testing associations between root foraging precision and root-associated AMF gene copies rather than the mechanism underlying root and hyphal foraging.

2 Materials and methods

2.1 Plant species selection, substrate, and mycorrhizal inoculum

We selected nine herbaceous plant species occurring in meadow habitats in Czechia: Centaurea jacea L. (Asteraceae), Lotus corniculatus L. and Melilotus officinalis (L.) Lam. (Fabaceae), Lythrum salicaria L. (Lythraceae), Linaria vulgaris Mill. (Plantaginaceae), Anthoxanthum odoratum L. and Phleum phleoides (L.) H.Karst. (Poaceae), Ranunculus acris L. (Ranunculaceae), and Potentilla recta L. (Rosaceae). All selected species form symbiosis with AMF. Their ability to establish mycorrhizal symbiosis with the inoculum used in this study was tested beforehand in a pre-experiment including 24 species (see Supplementary Table 1 in the Supporting Information). We considered the relative representation of different plant families according to their presence in the flora of Czechia and the different rates of root foraging precision established in our previous experiments (Weiser et al., 2016). We obtained plant seeds from a commercial supplier (Planta Naturalis Ltd., Markvartice, Czechia, www.plantanaturalis.com).

For the substrate in this experiment, we used a mixture of autoclaved (at 100–120 °C/6 h) zeolite, autoclaved sand, gamma-sterilised soil (min. 25 kGy), and mycorrhizal/non-mycorrhizal inoculum at the proportions of 43%, 43%, 9%, and 5%, respectively. The mycorrhizal inoculum consisted of spores of three fungal species – Funneliformis mosseae (T.H.Nicolson & Gerd.) C.Walker & A. Schuessler, Entrophospora claroidea (N.C.Schenck & G.S.Sm.) Błaszk., Niezgoda, B.T.Goto & Magurno, and Rhizophagus irregularis (Błaszk., Wubet, Renker & Buscot) C.Walker & A.Schüsler – produced separately, and mixed together with chopped fragments of leek roots (Allium porrum L.; less than 2 cm in length) colonised by these fungal species, and substrate from the leek rhizosphere containing extraradical mycelial fragments of the respective fungi. The non-mycorrhizal inoculum contained leek root fragments and substrate prepared in the same way but without AMF propagules.

2.2 Experiment

In the experiment, we applied two types of treatment – mycorrhizal or non-mycorrhizal inoculum and heterogeneous or homogeneous nutrient distribution – in a full-factorial design. Each of the four treatment combinations was replicated 10 times for each species, resulting in a total of 360 pots. For the nutrient treatments, we used a drip irrigator with flow-powered proportional dosing pumps (Dosatron, D25RE2) that added fertiliser to the water. Fertilised water was then delivered by the drippers to the pots at a constant amount at the same time. Each pot was supplied with two drippers placed at the opposite sides of the pot approximately 1 cm from the edge (Supplementary Figure 1B in the Supporting Information). In the heterogeneous treatment, a dripper supplied water with fertiliser to one side of the pot, whilst a second dripper supplied water without fertiliser to the other side. In the homogeneous treatment, each pot received water with fertiliser from both drippers. To create heterogeneous and homogeneous nutrient treatments, we changed the proportion of the fertiliser in the drippers. We used the recommended fertiliser concentrations for adult plants (0.1% v/v solution; Wuxal Super, NPK 8:8:6 + micronutrients, Aglukon; Supplementary Table 2 in the Supporting Information). The fertiliser was mixed with water (10:1) in storage barrels and then diluted further to obtain a 2% solution for the heterogeneous treatment and a 1% solution for the homogeneous treatment. Thus, both treatments received the same total amount of fertiliser. We placed the pots on water-levelled perforated tables to allow excess water to drain freely. We irrigated the pots twice a day for 10 minutes. The method used to create nutrient treatments is a modified version of Campbell and Grime (1989) nutrient-dripping method. The non-retentive substrate in the pots ensured that the main flow of the water with fertiliser was vertical, thus creating heterogeneity with only limited mixing in a horizontal direction without using a physical barrier.

Before the experiment, we sterilised the seeds with a solution of sodium hypochlorite (5% w/w) and water at a ratio of 1:10 for 10 minutes and rinsed them with demineralised water. We also sterilised the equipment and pots by washing them with 60% ethyl alcohol and wiping them with paper towels. On 22 April 2019, we sowed the seeds in Petri dishes containing white silica sand sterilised in an autoclave (120 °C/6 h). We placed the Petri dishes in a growth chamber (16:8 h day/night cycle with 20/15 °C temperature, respectively). After three weeks, we transplanted seedlings into small square pots (7 × 7 × 8 cm, 0.23 L) filled with the prepared substrate (sand, zeolite, soil, and mycorrhizal/non-mycorrhizal inoculum) and a nonwoven fabric at the bottom to prevent roots from growing out of the pot. Each of these small pots was perforated on two opposite sides with multiple 7-mm holes and placed inside a similar pot without perforations to keep the substrate in the dark. The perforated area was at least 60% of the pot wall surface (Supplementary Figure 1A in the Supporting Information). For each species, 40 small pots were used, half with mycorrhizal and half with non-mycorrhizal inoculum (in total 360 small pots). Seedlings grew in the small pots for six weeks in the greenhouse and were watered regularly with tap water. During this time, plants had enough time to establish symbiosis with mycorrhizal fungi.

On June 24, we transplanted the plants into round 3-L pots (TEKU Pöppelman, MCI 19) filled with the same substrate as the small pots. We removed each small, perforated pot from the non-perforated pot and placed it, with the plant, in the middle of a 3-L pot (Supplementary Figure 1B in the Supporting Information). We installed drip irrigation in the large pot and ensured that the perforated sides of the small pot faced the drippers. Roots and fungal hyphae could grow freely through the perforations of the small pot into the substrate in the large pot. We started the nutrient treatments as described above and let the plants grow for five weeks.

We harvested the plants after five weeks on 29 July. First, we removed the aboveground parts of the plant; we then cut out the small pot with a sharp knife (we cut off the roots that grew through the holes) and carefully removed it from the large pot. We harvested all the roots from the small pot by washing off the substrate. Then, we divided the rest of the substrate in the large pot by cutting it with a sharpened metal plate. The plane of the cut was vertical and perpendicular to the line connecting the positions of the two drippers, thus perpendicular to the nutrient gradient in the pots with the heterogeneous treatment. We thus obtained the nutrient-rich and nutrient-poor parts of the substrate, together with roots, from the heterogeneous treatment pots (or two similar parts from the homogeneous treatment pots). We harvested all roots from each half of the pot by washing off the substrate. Shoot and root biomass was dried at 65 °C for 48 hours and then weighed. The roots from the small central pots were subsequently analysed using qPCR. The experiments were performed in the greenhouse of the experimental garden of the Faculty of Science, Charles University (50.069N, 14.425E). The greenhouse was not equipped with automatic air conditioning; thus, temperature regulation was limited to ventilation, water misting, and shading with shade cloth. However, neither soil nor air temperature was measured during the experiment.

2.3 qPCR analyses

For each plant species, we performed qPCR analyses of dried root biomass collected from the small central pots (10 mycorrhizal and 10 non-mycorrhizal pots in either the heterogeneous or the homogeneous nutrient treatment; 40 pots per species). Dried roots were milled in a ball mill (MM 200, Retsch, Haan, Germany) to physically homogenise the root samples and increase the representativeness of subsamples. Then, we extracted DNA from 10–12.5 mg of milled root powder using the glassmilk method described by Gryndler et al. (2014). To assess DNA extraction efficiency, we added 1.8 × 1010 gene copies of an internal standard (a linearised plasmid carrying a fragment of cassava mosaic virus, GenBank accession AJ427910) to each sample before DNA extraction (Thonar et al., 2012). qPCR analyses were performed using Luna® Universal Probe qPCR Master Mix (New England Biolabs) on a Lightcycler 480 instrument (Roche). We followed the protocol described by Thonar et al. (2012), with species-specific markers for F. mosseae (marker ‘moss’), R. irregularis (marker ‘intra’), and E. claroidea (marker ‘clar’). All markers targeted the nuclear large ribosomal subunit (LSU) of the respective AMF taxa. We converted the qPCR results (Cq values) to root-associated LSU copy numbers of individual AMF taxa per unit weight (=mg) of roots, and we corrected them for extraction efficiency based on the recovery of the internal standard according to Thonar et al. (2012). One sample (belonging to the species L. salicaria, highlighted by red colour in Supplementary Table 3 in the Supporting Information) was excluded from the data analyses due to very low DNA recovery, most likely resulting from the small amount of available root material, which rendered the corrected qPCR estimate of mycorrhizal abundance unreliable.

2.4 Data preparation and analysis

In the analyses, the aboveground biomass was denoted ‘shoot biomass’. ‘Root system biomass’ referred to the biomass of roots from both halves of the large pot together with roots from the small central pot. ‘Whole plant biomass’ was the shoot and root system biomass combined. ‘Root:shoot (R:S) ratio’ was calculated from the shoot and root system biomass. All biomass data were log-transformed in all subsequent analyses to correct for non-normality. For calculation of root foraging precision, we used only root biomass from the separated halves of the substrate from the large pot. We calculated root foraging precision in each pot as the decadic logarithm of N/W, where N was the biomass of dried roots from the nutrient-rich half of a given pot and W was the biomass of dried roots from the nutrient-poor half of the same pot. Values above zero indicated that plants allocated more roots to the nutrient-rich half of the pot, whereas values below zero indicated that plants allocated more roots to the nutrient-poor half of the pot. In the case of homogeneous pots, we assigned N or W randomly to one of the two halves. Thirty-seven plants died or were lost during sampling, three samples contained missing qPCR values, and one sample was excluded after qPCR analysis due to very low DNA recovery; thus, we used 319 pots from the original 360 pots in the analyses (Supplementary Table 4 in the Supporting Information). ‘Root-associated AMF LSU copy number’ in small central pots was the gene copy number of LSU of individual AMF taxa (F. mosseae, E. claroidea, and R. irregularis) per unit weight of roots summed together. In the case of fungal presence, 98% of AMF-inoculated pots showed a positive presence of F. mosseae, 51% showed the presence of R. irregularis, and 23% showed the presence of E. claroidea (Supplementary Tables 3, 7 in the Supporting Information). As E. claroidea was not detected in 77% of the AMF-inoculated pots, we decided to analyse all three AMF taxa together. Although AMF contamination (or false positive detection) was encountered in some non-inoculated control pots (numbers were not zero), root-associated AMF LSU copy number was lower than in most successfully colonised inoculated pots (Supplementary Table 3 in the Supporting Information, Figure 1). For the variable ‘root-associated AMF LSU copy number’, we applied a cube-root transformation in all subsequent analyses to correct for non-normality of the data. This transformation reduced the values from millions of root-associated AMF LSU copy numbers to thousands in units on the transformed scale. We chose cube-root transformation because it effectively decreases right skewness and can also be applied to zeros, which were frequent in our dataset.

Figure 1

To test whether the inoculation was successful, we used a linear model to explore the effect of mycorrhizal treatment (mycorrhizal or non-mycorrhizal inoculation), in combination with species identity and nutrient treatment, on root-associated AMF LSU copy number of each pot. Then, we estimated the effects of mycorrhizal treatment (mycorrhizal/non-mycorrhizal pots or the root-associated AMF LSU copy number), species identity, and nutrient treatment (both individually and in combination) on plant biomass (root system biomass, shoot biomass, whole plant biomass, and R:S ratio) and root foraging precision using linear models. Furthermore, we explored root and shoot biomass, root foraging precision and mycorrhizal colonisation at the level of individual species. We determined partial correlations of root and shoot biomass with mycorrhizal treatment as a category and partial correlations of root foraging precision with root-associated AMF LSU copy number for individual species using Welch’s two-sample t-test. We used data from mycorrhizal and non-mycorrhizal pots for analyses of root-associated AMF LSU copy number even though there was quite a high number of zeros. The results of analysis of data only with the mycorrhizal treatment were almost the same (Figure 2, Supplementary Table 5 in the Supporting Information). We fit all models using R (R Core Team, 2024, version 4.4.0) in Rstudio (Posit Team, 2024, version 2024.4.1.748). All images were created using the ggplot2 package (Wickham, 2016) and the cowplot package (Wilke et al., 2024).

Figure 2

3 Results

For the analyses that included plant biomass as the response variable, we show results only for root biomass and shoot biomass as they were highly correlated with whole-plant biomass (Supplementary Table 6 in the Supporting Information). The results of all the analyses were similar.

3.1 Success of mycorrhizal inoculation

The mycorrhizal pots showed significantly higher root-associated AMF LSU copy number compared to non-mycorrhizal pots, although the numbers were highly variable in the mycorrhizal treatment and the level of difference between mycorrhizal and non-mycorrhizal pots depended also on the plant species (Figure 1, Table 1, Supplementary Table 3 in the Supporting Information). Root-associated AMF LSU copy number also differed between plant species irrespectively of the mycorrhizal inoculation treatment. The mycorrhizal inoculation was measured only using qPCR analyses; we did not measure microscopical colonization data, extraradical hyphal density or plant tissue nutrient concentration to evaluate the success of mycorrhizal inoculation.

Table 1

Source of variation (factor/factor interaction)DfSum SqMean SqF valuePr (>F)
Mycorrhiza15990159901234.723<0.001
Species81278215986.261<0.001
Nutrient186860.3380.561
Mycorrhiza × Species8884611064.333<0.001
Mycorrhiza × Nutrient15815812.2760.133
Species × Nutrient838754841.8980.06
Mycorrhiza × Species × Nutrient817882230.8760.537
Residuals28372222255

The effects of mycorrhizal inoculation (Mycorrhiza: inoculated vs. non-inoculated), plant species identity (Species) and nutrient treatment (Nutrient: homogeneous vs. heterogeneous) on the root-associated AMF LSU copy number of the three mycorrhizal fungal inoculants per unit weight of roots.

Significant results are in bold. Root-associated AMF LSU copy number was cube-root transformed for the statistical analysis.

3.2 Biomass and root foraging in different treatments

Species identity, nutrient treatment, and mycorrhizal treatment (presence/absence of AMF) in interaction with species identity significantly affected root (and shoot) biomass of some plant species (Table 2, Figure 3). The root biomass and shoot biomass of R. acris were significantly higher in the mycorrhizal treatment (root: t = 3.01, p < 0.01; shoot: t = 2.26, p < 0.05) compared to the non-mycorrhizal treatment. In contrast, root biomass, but not shoot biomass, of L. vulgaris was significantly lower in the mycorrhizal treatment (t = -1.95, p < 0.05). Root biomass was slightly higher in the heterogeneous nutrient treatment, but not significantly so, compared to the homogeneous nutrient treatment in the test of nutrient treatments only (t = -1.00, p > 0.05). R:S ratio differed significantly only between species; mycorrhizal treatment and nutrient treatment and any of the possible interactions had no effect on R:S ratio (Table 2).

Table 2

Response variableSource of variation (factor/factor interaction)DfSum SqMean SqF valuePr (>F)
Root biomassMycorrhiza10.010.040.0570.811
Species8371.846.4873.812<0.001
Nutrient12.82.824.4730.035
Mycorrhiza × Species8182.243.564<0.001
Mycorrhiza × Nutrient111.041.6540.199
Species × Nutrient80.60.070.1170.999
Mycorrhiza × Species × Nutrient86.10.771.220.287
Residuals282177.60.63
Shoot biomassMycorrhiza100.030.0480.827
Species8408.151.0186.545<0.001
Nutrient11.71.672.8370.093
Mycorrhiza × Species813.81.732.9370.004
Mycorrhiza × Nutrient10.80.851.4370.232
Species × Nutrient82.20.280.4690.878
Mycorrhiza × Species ×Nutrient86.50.811.3720.208
Residuals282166.20.59
R:S ratioMycorrhiza10.130.1250.7640.383
Species853.996.74941.19<0.001
Nutrient10.010.0060.0390.845
Mycorrhiza × Species80.810.1010.6160.765
Mycorrhiza × Nutrient10.150.1470.8950.345
Species × Nutrient80.80.10.6080.771
Mycorrhiza × Species × Nutrient81.510.1891.1530.328
Residuals28146.040.164
Root-foraging precisionMycorrhiza10.130.1310.6810.41
Species81.110.1390.7250.67
Nutrient17.257.25137.802<0.001
Mycorrhiza × Species81.670.2081.0870.373
Mycorrhiza × Nutrient10.010.0060.0310.861
Species × Nutrient83.440.432.2430.025
Mycorrhiza × Species × Nutrient81.550.1941.0130.426
Residuals28354.280.192

The effects of mycorrhizal inoculation (Mycorrhiza: inoculated vs. non-inoculated), plant species identity (Species), nutrient treatment (Nutrient: homogeneous vs. heterogeneous) and their interactions on root and shoot biomass, root:shoot (R:S) ratio and root foraging precision of nine plant species.

Significant results are in bold. Root and shoot biomass and R:S ratio were log-transformed for the statistical analyses.

Figure 3

Root foraging precision differed between heterogeneous and homogeneous nutrient treatments (it was higher in the heterogeneous treatment) and between species (the nutrient treatment × species identity interaction was significant) (Table 2, Supplementary Figure 2 in the Supporting Information). All species showed the same trend of higher root foraging precision in the heterogeneous nutrient treatment compared to the homogeneous treatment (significant results: C. jacea: t = -3.15, p < 0.01; L. corniculatus: t = -3.89, p < 0.001; L. salicaria: t = -2.45, p < 0.05; M. officinalis: t = -2.79, p < 0.05; P. recta: t = -6.18, p < 0.001). Mycorrhizal treatment had no effect on root foraging precision (Table 2).

However, some pots assigned to the mycorrhizal treatment were almost non-mycorrhizal (Figure 1, Supplementary Table 3 in the Supporting Information). Thus, it was reasonable to test the presence and abundance of mycorrhizal fungi using root-associated AMF LSU copy number. Mycorrhizal treatment (root-associated AMF LSU copy number) significantly affected root and shoot biomass (Table 3, Figure 2). Both root biomass and shoot biomass of plants declined with increasing root-associated AMF LSU copy number (root biomass: F1, 316 = 5.93, p < 0.05,  = 0.018; shoot biomass: F1,316 = 6.76, p < 0.01,  = 0.021). Results for R:S ratio were almost similar whether the root-associated AMF LSU copy number was used or whether AMF presence/absence was used in the analysis. However, the nutrient treatment × mycorrhiza interaction had a weak but significant effect when the root-associated AMF LSU copy number was used (Table 3). The R:S ratio of plants in the heterogeneous nutrient treatment was positively correlated with the root-associated AMF LSU copy number, whereas for plants in the homogeneous treatment, the trend was the opposite, although these effects were not significant when homogeneous and heterogeneous nutrient treatments were analysed separately. Moreover, the significant effect of the nutrient treatment × mycorrhiza interaction on the R:S ratio disappeared when we used only mycorrhizal inoculated pots for the analyses (Supplementary Table 5 in the Supporting Information).

Table 3

Response variableSource of variation (factor/factor interaction)DfSum SqMean SqF valuePr (>F)
Root biomassMycorrhiza110.610.6415.561<0.001
Species8360.945.1165.948<0.001
Nutrient12.82.814.1070.044
Mycorrhiza × Species86.80.851.2490.271
Mycorrhiza × Nutrient11.21.21.7570.186
Species × Nutrient80.40.050.0721.000
Mycorrhiza × Species × Nutrient82.30.290.4210.908
Residuals282192.90.68
Shoot biomassMycorrhiza112.612.5619.704<0.001
Species839549.3877.475<0.001
Nutrient11.71.732.7070.101
Mycorrhiza × Species85.40.671.0540.396
Mycorrhiza × Nutrient10.30.270.4290.513
Species × Nutrient81.80.230.3550.943
Mycorrhiza × Species × Nutrient82.90.360.5690.803
Residuals282179.70.64
R:S ratioMycorrhiza10.010.0050.0330.856
Species854.096.76242.62<0.001
Nutrient100.0040.0260.873
Mycorrhiza × Species82.050.2571.6180.119
Mycorrhiza × Nutrient10.80.7985.030.026
Species × Nutrient80.690.0860.540.826
Mycorrhiza × Species × Nutrient81.20.1510.9490.476
Residuals28144.580.159
Root-foraging precisionMycorrhiza10.190.1881.0260.312
Species81.090.1370.7470.65
Nutrient17.217.21239.351<0.001
Mycorrhiza × Species84.340.5422.9590.003
Mycorrhiza × Nutrient1000.0020.967
Species × Nutrient83.470.4342.3680.018
Mycorrhiza × Species × Nutrient81.280.160.8730.539
Residuals28351.860.183

Linear model analyses of the effects of mycorrhizal root colonisation (Mycorrhiza: root-associated AMF LSU copy number), plant species identity (Species) and nutrient treatment (Nutrient: homogeneous vs. heterogeneous) on the root and shoot biomass, root:shoot (R:S) ratio and root foraging precision of plants.

Mycorrhizal colonisation of roots was expressed as the root-associated AMF LSU copy number of each of the three mycorrhizal inoculants per unit weight of roots. Significant results are in bold. Root and shoot biomass and R:S ratio were log-transformed. The root-associated AMF LSU copy number was cube-root transformed.

Root foraging was affected by the mycorrhiza × species identity interaction when the root-associated AMF LSU copy number was used in the analysis of mycorrhizal treatment (Table 3, Figure 4). This effect was only marginally significant when we used only mycorrhizal inoculated pots for the analyses (Supplementary Table 5 in the Supporting Information) The root foraging precision of M. officinalis was significantly lower in pots with a higher root-associated AMF LSU copy number (F1,10 = 12.44, p < 0.01,  = 0.554). Other species showed no significant effect.

Figure 4

4 Discussion

Most vascular plants form some type of mycorrhizal symbiosis, the majority of the cases being the arbuscular type of mycorrhiza (Brundrett and Tedersoo, 2018). By exchanging nutrients and water for carbon products of photosynthesis, the symbiotic fungi affect the growth and nutrient acquisition strategies of the plant. One such strategy is plant roots’ precise foraging for heterogeneously distributed nutrients (Campbell et al., 1991; Kembel and Cahill, 2005). However, previous experimental evidence of mycorrhizal fungi affecting a plant’s ability to forage precisely for nutrients is contradictory (Cui and Caldwell, 1996a; Farley and Fitter, 1999b; Hodge, 2001; Wijesinghe et al., 2001; Felderer et al., 2013). In this study, we addressed the question of whether presence of mycorrhizal inoculum and root-associated AMF LSU copy number affected plant foraging precision by testing nine plant species in heterogeneous soil with or without mycorrhizal fungi.

4.1 Root foraging precision

In this study, plant species foraged for soil nutrients with varying degrees of precision as was also shown in previous studies (Campbell et al., 1991; Kembel and Cahill, 2005; Weiser et al., 2016). However, mycorrhizal fungi did not generally affect the root foraging precision of our species selection. The only species affected by mycorrhizal fungi was M. officinalis; it showed decreased root foraging precision with higher root-associated AMF LSU copy number in the central part of the root system. In previous studies, the effect of mycorrhizal fungi on root foraging precision was variable: there was either no effect of fungi on root foraging (Farley and Fitter, 1999b; Hodge, 2001; Wijesinghe et al., 2001) or plants showed reduced root foraging precision in association with mycorrhizal fungi (Cui and Caldwell, 1996a; Felderer et al., 2013). The observed reduction in root foraging precision of M. officinalis with higher root-associated AMF LSU copy number in the central part of the root system corresponds to the findings of studies on Lotus japonicus (Felderer et al., 2013) and Agropyron desertorum (Cui and Caldwell, 1996a). These studies showed a decrease in root foraging only in phosphorus (P) patches, whilst we investigated foraging for a mixture of nutrients (nitrogen [N], P, potassium [K] + micronutrients). As a leguminous species, M. officinalis is limited mainly by P availability in the soil, since it obtains most of its N through its association with symbiotic N2-fixing bacteria. Therefore, even though we tested foraging responses to all three major nutrients (N, P, and K), M. officinalis could have been limited only (or mainly) by the P. Consequently, if mycorrhizal fungi contributed to P acquisition of M. officinalis, the plant may have had sufficient amounts of all essential nutrients and could reduce its root foraging activity. In contrast, other species in the study might have remained limited by N and thus continued to forage even when associated with mycorrhizal fungi. However, we did not directly measure nodulation by nitrogen-fixing bacteria, nitrogen fixation activity, plant tissue N:P ratios, or responses to isolated N, P, or K patches here. Moreover, the second leguminous species in the study, L. corniculatus, did not exhibit the same trend as M. officinalis; it did not show a change in root foraging precision after mycorrhizal inoculation. Thus, it remains unclear whether the pattern observed in M. officinalis reflected P limitation resulting from the association with nitrogen-fixing bacteria, and this hypothesis requires further investigation.

These contrasting results raise a question that has not yet been fully answered: Which nutrients are plant species actually foraging for? Previous studies have shown that some species forage with varying degrees of precision for major nutrients, particularly N, P and K (Drew, 1975; Campbell et al., 1991), as well as specifically for ammonium (Drew, 1975), calcium and magnesium (Skálová et al., 2025). The use of a nutrient mixture in our experiment helped to overcome this uncertainty; however, it came at the cost of a more complex interpretation of the results, as it did not allow us to disentangle foraging responses to individual nutrients.

In soil, N is usually more abundant and mobile than P (Farley and Fitter, 1999a). Therefore, plants are often able to acquire N effectively on their own, without the assistance of AMF, which they rely on primarily for P uptake (Smith et al., 2011). One possible explanation for our results is that, in nutrient patches containing a mixture of nutrients, plants may have responded primarily to N availability, whilst AMF may have been more strongly associated with P acquisition. This might explain why we did not observe a reduction in root foraging precisions, except in M. officinalis, despite the observed values of root-associated AMF LSU copy number. However, this remains unclear, as it is not supported by direct measurements of extraradical hyphal development, fungal nutrient uptake, or nutrient transfer to plants. Another possible explanation for the absence of an effect of mycorrhizal inoculation on root foraging precision is the generally lower level of root-associated AMF LSU copy number, similarly to Farley and Fitter (1999b).

This study partially supports previous research (Stiblíková et al., 2023) that linked root-foraging precision to the root economic spectrum but did not find a relationship between foraging precision and mycorrhizal association or other root traits. Under optimal experimental conditions, root-foraging precision appears to be a species-specific trait that could theoretically be positioned along the collaboration gradient of the root economic spectrum (Bergmann et al., 2020). However, a theoretical study involving 123 species (Stiblíková et al., 2023) suggested that root-foraging precision constitutes an independent axis within this spectrum; our results may thus support this conclusion. In our study, when the fungal partner was removed from reportedly mycorrhizal plant species (Guerrero-Ramirez et al., 2021), root-foraging precision generally remained unaffected; only plant biomass was altered. Clearly, additional factors, such as nutrient type or form, could play a role in the allocation between root and fungi.

The balance between allocation of resources (primarily the carbon) to roots or fungal partner may play a role in the natural distribution of plant species over the course of succession as it depends on nutrient availability. Plants are more likely to rely on mycorrhizal hyphae to gain nutrients in nutrient-poor environments, whilst they tend to use roots in nutrient-rich sites (Hoepfner et al., 2015a). Indeed, it is more advantageous for plants to acquire nutrients using mycorrhizal fungi than their own fine roots in sand ecosystems (Hoepfner et al., 2015b). In our study, the nutrient supply was constant and high, which could have led to plants not requiring mycorrhizal symbionts for sufficient nutrient uptake; thus, different results may be obtained under nutrient deficiency.

4.2 Plant biomass allocation

Although our main focus was on the effects of mycorrhizal inoculation on root foraging precision, we also analysed plant biomass allocation to better understand plant responses associated with AMF. Biomass declined or did not change as root-associated AMF LSU copy number increased. The effect was similar for both aboveground and belowground biomass. Mycorrhizal symbiosis should not be regarded unequivocally as a mutualistic relationship. Plants’ growth responses may range from positive to neutral or negative, thus the relationship is a continuum from parasitism to mutualism (Johnson et al., 1997). Whether mycorrhizal fungi help plants or parasitise them also depends on the environment and local adaptation (Klironomos, 2003), as well as nutrient and water availability (Mason et al., 2000; Lekberg and Koide, 2008; Frew, 2023), both of which was quite high in our experiment. At least one fungal species, Rhizophagus irregularis, has a more mutualistic effect on plant P acquisition under drought conditions (Püschel et al., 2021). We used this species in our study, but it was grown under sufficient moisture.

Interestingly, plant investment in roots compare to shoots in our study did not change under the mycorrhizal or nutrient treatment. The root-associated AMF LSU copy number affected the R:S ratio only in interaction with nutrient treatments; plus, this effect was not apparent in the analyses only with pots to which mycorrhizal inoculum was added – thus without zeroes (corresponding to plants without detectable root-associated AMF LSU copy number) in non-mycorrhizal inoculum. The expected relationship with mutualistic fungi, i.e. higher biomass investment to shoots when the plant has a mycorrhizal partner that helps it to acquire nutrients, occurred in the homogeneous treatment. The effect in the heterogeneous treatment was the opposite, i.e. plants invested more biomass in roots when the root-associated AMF LSU copy number per unit weight of roots was higher, but this effect was rather weak and possibly spurious. We assessed mycorrhizal inoculation using the root-associated AMF LSU copy number per unit root mass in the central part of the root system as a proxy for the presence and abundance of AMF DNA. It does not necessarily indicate mycorrhizal activity outside of this zone and mycorrhiza-mediated nutrient transfer to the host plant. We did not measure arbuscule abundance or extraradical hyphal density; thus, we cannot determine whether differences in root-associated AMF LSU copy number reflected differences in the functional activity of the mycorrhizal fungi.

The reduction in plant biomass in the presence of mycorrhizal fungi in our foraging experiment was also observed by Wijesinghe et al. (2001). In their experiment, the mycorrhizal abundance was fairly low; thus, the effects of fungi would have hardly led to any prominent benefits. They observed higher biomass investment in shoots in the presence of mycorrhiza, but only in species with the highest level of mycorrhizal colonisation. Unlike our study, Wijesinghe et al. (2001) also measured biomass investment in flowers and flowering shoots compared to that in roots. A previous foraging study using the same fungal species as ours (Felderer et al., 2013) observed an increase in plant biomass in the mycorrhizal treatment, but the level of root-associated AMF LSU copy number was much higher compared to that in our experiment. Other previous studies of root foraging under mycorrhizal and non-mycorrhizal treatments also found no change in the R:S ratio in response to the combined mycorrhizal and nutrient treatment.

4.3 Limitations

The extent of mycorrhizal inoculation of plant species differed in our study despite giving plant roots and fungal hyphae sufficient time to establish mycorrhizal symbiosis. The number of pots with detectable root-associated AMF LSU copy number ranged from one to all inoculated pots (max. 19, Supplementary Tables 3 and 7 in the Supporting Information) for both plant and fungal species. The root-associated AMF LSU copy number in pots within the same treatment also differed widely. All plant species used in the study are naturally mycorrhizal but their mycorrhizal responsiveness differs (Soudzilovskaia et al., 2020) and, consequently, so may their belowground carbon allocation trade-off, i.e. the allocation to either roots or fungal hyphae, which depends on it (Unger et al., 2016). Moreover, the fungal symbionts we provided might not have suited all plant species, although we used a mixture of three fungal species which imitates natural conditions better than the use of a single fungal species. Usually, most AMF can enter symbiosis with most host species, although specific plant–fungus combinations are probably more profitable than others and affect the growth of plants and fungi differently (Bever et al., 1996; Klironomos, 2003).

Beyond natural factors, the experimental water and temperature regime might have contributed to the relatively low values of detected root-associated AMF LSU copy number, although its influence was not measured and tested. Frequent irrigation (a condition necessary for stable nutrient patch formation) and occasionally high temperatures might have contributed to reduced AMF establishment (Mason et al., 2000; Lekberg and Koide, 2008; Püschel et al., 2021; Frew, 2023; but see Stock et al., 2021). In addition to the suboptimal conditions for the development of mycorrhizal symbiosis, another reason for rather low values of root-associated AMF LSU copy number could be the overall favourable conditions for plants even without mycorrhizal support (enough water and soluble nutrients); thus, the benefits of the symbiosis may have been limited. Absence of strong N and/or P limitation in the substrate may thus have led to decreased allocation by plants to AMF (Johnson et al., 2003), although it should be noted that we did not directly measure nutrient concentrations in either the substrate or the plants. The soil N:P ratio may affect plant investment in AMF, although this may be species-specific depending on the plants’ mycorrhizal responsiveness (Unger et al., 2016).

4.4 Conclusion

Our results suggest that root-associated AMF LSU copy number did not consistently affect the root foraging precision of the nine herbaceous plant species we investigated. Although one species (M. officinalis) showed reduced root foraging precision with increasing root-associated AMF LSU copy number in the central part of the root system, the remaining species were unaffected, suggesting that the role of mycorrhizal inoculation in root nutrient foraging is neither universal nor particularly common. These results support the view that root foraging precision is a species-specific trait that is not necessarily aligned with the collaboration gradient of the root economic spectrum but instead may represent an independent axis within it (Stiblíková et al., 2023). The removal of the fungal partner did not consistently and strongly affect root foraging precision, indicating that factors such as nutrient type rather than the position of the plant species in the collaboration gradient may influence the observed relationship between root foraging precision and root-associated AMF LSU copy number.

However, our results might also have been affected by variability in root-associated AMF LSU copy number. The growing conditions might have shifted the mycorrhizal symbiosis from mutualistic towards a more parasitic interaction, as indicated by reduced plant biomass at higher levels of root-associated AMF LSU copy number. The relationship between root foraging precision and root-associated AMF LSU copy number likely depends on several factors, including the amount, distribution, and type of nutrients (Hoepfner et al., 2015a; Wang et al., 2016). Future research should investigate how root foraging responses are associated with extraradical hyphal development and nutrient transfer under varying degrees of nutrient limitation, with a focus on the spatial separation of individual nutrients (especially the N and P). A better grasp of these dynamics is essential for predicting plant nutrient acquisition strategies across diverse environmental contexts which helps to understand the root economic spectrum of plants.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

PS: Writing – review & editing, Conceptualization, Investigation, Writing – original draft, Visualization, Data curation, Project administration, Formal Analysis, Methodology. JJ: Methodology, Data curation, Conceptualization, Project administration, Investigation, Funding acquisition, Writing – review & editing, Resources. MR: Data curation, Methodology, Writing – review & editing. MK: Methodology, Writing – review & editing, Data curation. MB: Writing – review & editing, Methodology. JŠ: Writing – review & editing, Methodology. MW: Investigation, Funding acquisition, Resources, Writing – review & editing, Methodology, Formal Analysis, Project administration, Supervision, Conceptualization.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research was supported by the Czech Science Foundation No. 19-0630S.

Acknowledgments

We acknowledge and thank the colleagues that helped us with the realization of the experiments. We thank the language editor D Large, PhD for improving the language and style of the manuscript.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fpls.2026.1899078/full#supplementary-material

References

  • 1

    AllenM. F. (2007). Mycorrhizal fungi: Highways for water and nutrients in arid soils. Vadose Zone J.6, 291297. doi: 10.2136/vzj2006.0068

  • 2

    BergmannJ.WeigeltA.PlasF.LaughlinD. C.KuyperT. W.Guerrero-RamirezN.et al. (2020). The fungal collaboration gradient dominates the root economics space in plants. Sci. Adv.6, eaba3756. doi: 10.1126/sciadv.aba3756

  • 3

    BeverJ. D.MortonJ. B.AntonovicsJ.SchultzP. A. (1996). Host-dependent sporulation and species diversity of arbuscular mycorrhizal fungi in a mown grassland. J. Ecol.84, 7182. doi: 10.2307/2261701

  • 4

    BrundrettM. C. (2009). Mycorrhizal associations and other means of nutrition of vascular plants: Understanding the global diversity of host plants by resolving conflicting information and developing reliable means of diagnosis. Plant Soil320, 3777. doi: 10.1007/s11104-008-9877-9

  • 5

    BrundrettM. C.TedersooL. (2018). Evolutionary history of mycorrhizal symbioses and global host plant diversity. New Phytol.220, 11081115. doi: 10.1111/nph.14976

  • 6

    CampbellB. D.GrimeJ. P. (1989). A new method of exposing developing root systems to controlled patchiness in mineral nutrient supply. Ann. Bot.63, 395400. doi: 10.1093/oxfordjournals.aob.a087758

  • 7

    CampbellB. D.GrimeJ. P.MackeyJ. M. L. (1991). A trade-off between scale and precision in resource foraging. Oecologia87, 532538. doi: 10.1007/BF00320417

  • 8

    CavagnaroT. R.SmithF. A.SmithS. E.JakobsenI. (2005). Functional diversity in arbuscular mycorrhizas: Exploitation of soil patches with different phosphate enrichment differs among fungal species. Plant Cell Environ.28, 642650. doi: 10.1111/j.1365-3040.2005.01310.x

  • 9

    ChenW.KoideR. T.AdamsT. S.DeForestJ. L.ChengL.EissenstatD. M. (2016). Root morphology and mycorrhizal symbioses together shape nutrient foraging strategies of temperate trees. Proc. Natl. Acad. Sci. U.S.A.113, 87418746. doi: 10.1073/pnas.1601006113

  • 10

    ChengL.ChenW.AdamsT. S.WeiX.LiL.McCormackM. L.et al. (2016). Mycorrhizal fungi and roots are complementary in foraging within nutrient patches. Ecology97, 28152823. doi: 10.1002/ecy.1514

  • 11

    CuiM.CaldwellM. M. (1996a). Facilitation of plant phosphate acquisition by arbuscular mycorrhizas from enriched soil patches. I. Roots and hyphae exploiting the same soil volume. New Phytol.133, 453460. doi: 10.1111/j.1469-8137.1996.tb01912.x

  • 12

    CuiM. Y.CaldwellM. M. (1996b). Facilitation of plant phosphate acquisition by arbuscular mycorrhizas from enriched soil patches .2. Hyphae exploiting root-free soil. New Phytol.133, 461467. doi: 10.1111/j.1469-8137.1996.tb01913.x

  • 13

    DrewM. C. (1975). Comparison of the effects of a localised supply of phosphate, nitrate, ammonium and potassium on the growth of the seminal root system, and the shoot, in barley. New Phytol.75, 479490. doi: 10.1111/j.1469-8137.1975.tb01409.x

  • 14

    EissenstatD.M KucharskiJ.ZadwornyM.S. AdamsT.KoideR. (2015). Linking root traits to nutrient foraging in arbuscular mycorrhizal trees in a temperate forest. New Phytol.208, 114–124. doi: 10.1111/nph.13451

  • 15

    FarleyR. A.FitterA. H. (1999a). Temporal and spatial variation in soil resources in a deciduous woodland. J. Ecol.87, 688696. doi: 10.1046/j.1365-2745.1999.00390.x

  • 16

    FarleyR. A.FitterA. H. (1999b). The responses of seven co-occurring woodland herbaceous perennials to localized nutrient-rich patches. J. Ecol.87, 849859. doi: 10.1046/j.1365-2745.1999.00396.x

  • 17

    FeldererB.JansaJ.SchulinR. (2013). Interaction between root growth allocation and mycorrhizal fungi in soil with patchy P distribution. Plant Soil373, 569582. doi: 10.1007/s11104-013-1818-6

  • 18

    FrewA. (2023). Water availability alters the community structure of arbuscular mycorrhizal fungi and determines plant mycorrhizal benefit. Plants People Planet5, 683689. doi: 10.1002/ppp3.10372

  • 19

    FrieseC. F.AllenM. F. (1991). The spread of VA mycorrhizal fungal hyphae in the soil: Inoculum types and external hyphal architecture. Mycologia83, 409418. doi: 10.2307/3760351

  • 20

    GryndlerM.ČernáL.BukovskáP.HršelováH.JansaJ. (2014). Tuber aestivum association with non-host roots. Mycorrhiza24, 603610. doi: 10.1007/s00572-014-0580-9

  • 21

    Guerrero-RamirezN. R.MommerL.FreschetG. T.IversenC. M.McCormackM. L.KattgeJ.et al. (2021). Global root traits (GRooT) database. Glob. Ecol. Biogeogr.30, 2537. doi: 10.1111/geb.13179

  • 22

    HodgeA. (2001). Arbuscular mycorrhizal fungi influence decomposition of, but not plant nutrient capture from, glycine patches in soil. New Phytol.151, 725734. doi: 10.1046/j.0028-646x.2001.00200.x

  • 23

    HodgeA. (2006). Plastic plants and patchy soils. J. Exp. Bot.57, 401411. doi: 10.1093/jxb/eri280

  • 24

    HoepfnerI.BeyschlagW.BartelheimerM.WernerC.UngerS. (2015a). Role of mycorrhization and nutrient availability in competitive interactions between the grassland species Plantago lanceolata and Hieracium pilosella. Plant Ecol.216, 887899. doi: 10.1007/s11258-015-0476-6

  • 25

    HoepfnerI.FriedeM.UngerS.BeyschlagW. (2015b). Potential advantages of highly mycotrophic foraging for the establishment of early successional pioneer plants on sand. Funct. Plant Biol.42, 95104. doi: 10.1071/FP14097

  • 26

    JohnT. V. S.ColemanD. C.ReidC. P. P. (1983). Growth and spatial distribution of nutrient-absorbing organs: Selective exploitation of soil heterogeneity. Plant Soil71, 487493. doi: 10.1007/BF02182690

  • 27

    JohnsonN. C. (2010). Resource stoichiometry elucidates the structure and function of arbuscular mycorrhizas across scales. New Phytol.185, 631647. doi: 10.1111/j.1469-8137.2009.03110.x

  • 28

    JohnsonN. C.GrahamJ.-H.SmithF. A. (1997). Functioning of mycorrhizal associations along the mutualism–parasitism continuum*. New Phytol.135, 575585. doi: 10.1046/j.1469-8137.1997.00729.x

  • 29

    JohnsonN. C.RowlandD. L.CorkidiL.Egerton-WarburtonL. M.AllenE. B. (2003). Nitrogen enrichment alters mycorrhizal allocation at five mesic to semiarid grasslands. Ecology84, 18951908. doi: 10.1890/0012-9658(2003)084[1895:NEAMAA]2.0.CO;2

  • 30

    KembelS.CahillJ. (2005). Plant phenotypic plasticity belowground: A phylogenetic perspective on root foraging trade‐offs. Am. Nat.166, 216–230. doi: 10.1086/431287

  • 31

    KlironomosJ. N. (2003). Variation in plant response to native and exotic arbuscular mycorrhizal fungi. Ecology84, 22922301. doi: 10.1890/02-0413

  • 32

    LekbergY.KoideR. T. (2008). Effect of soil moisture and temperature during fallow on survival of contrasting isolates of arbuscular mycorrhizal fungi. Botany86, 11171124. doi: 10.1139/B08-077

  • 33

    LiuB.LiH.ZhuB.KoideR. T.EissenstatD. M.GuoD. (2015). Complementarity in nutrient foraging strategies of absorptive fine roots and arbuscular mycorrhizal fungi across 14 coexisting subtropical tree species. New Phytol.208, 125136. doi: 10.1111/nph.13434

  • 34

    MasonP. A.IbrahimK.InglebyK.MunroR. C.WilsonJ. (2000). Mycorrhizal development and growth of inoculated Eucalyptus globulus (Labill.) seedlings in wet and dry conditions in the glasshouse. For. Ecol. Manage.128, 269277. doi: 10.1016/S0378-1127(99)00155-3

  • 35

    Posit Team (2024). RStudio: Integrated Development Environment for R (Boston, MA: Posit Software, PBC). Available online at: http://www.posit.co/.

  • 36

    PüschelD.BitterlichM.RydlováJ.JansaJ. (2021). Drought accentuates the role of mycorrhiza in phosphorus uptake. Soil Biol. Biochem.157, 108243. doi: 10.1016/j.soilbio.2021.108243

  • 37

    R Core Team (2024). R: A Language and Environment for Statistical Computing (Vienna, Austria: R Foundation for Statistical Computing). Available online at: https://www.R-project.org/.

  • 38

    ShiZ.WangF.ZhangC.YangZ. (2011). Exploitation of phosphorus patches with different phosphorus enrichment by three arbuscular mycorrhizal fungi. J. Plant Nutr.34, 10961106. doi: 10.1080/01904167.2011.558154

  • 39

    SkálováH.PánkováK.StiblíkováP.KřivohlavýF.HadincováV.TylováE.et al. (2025). Root foraging for calcium and magnesium in common meadow plant species. Plant Soil. 517, 1605–21. doi: 10.1007/s11104-025-07941-2

  • 40

    ŠmilauerP.ŠmilauerováM.KotilínekM.KošnarJ. (2020). Foraging speed and precision of arbuscular mycorrhizal fungi under field conditions: An experimental approach. Mol. Ecol.29, 15741587. doi: 10.1111/mec.15425

  • 41

    SmithS. E.JakobsenI.GrønlundM.SmithF. A. (2011). Roles of arbuscular mycorrhizas in plant phosphorus nutrition: Interactions between pathways of phosphorus uptake in arbuscular mycorrhizal roots have important implications for understanding and manipulating plant phosphorus acquisition. Plant Physiol.156, 10501057. doi: 10.1104/pp.111.174581

  • 42

    SoudzilovskaiaN. A.VaessenS.BarceloM.HeJ.RahimlouS.AbarenkovK.et al. (2020). FungalRoot: Global online database of plant mycorrhizal associations. New Phytol.227, 955966. doi: 10.1111/nph.16569

  • 43

    StiblíkováP.KlimešA.CahillJ. F.KoubekT.WeiserM. (2023). Interspecific differences in root foraging precision cannot be directly inferred from species’ mycorrhizal status or fine root economics. Oikos. 2023 (1), e08995. doi: 10.1111/oik.08995

  • 44

    StockS. C.KoesterM.BoyJ.GodoyR.NajeraF.MatusF.et al. (2021). Plant carbon investment in fine roots and arbuscular mycorrhizal fungi: A cross-biome study on nutrient acquisition strategies. Sci. Total Environ.781, 146748. doi: 10.1016/j.scitotenv.2021.146748

  • 45

    ThonarC.ErbA.JansaJ. (2012). Real-time PCR to quantify composition of arbuscular mycorrhizal fungal communities—marker design, verification, calibration and field validation. Mol. Ecol. Resour.12, 219232. doi: 10.1111/j.1755-0998.2011.03086.x

  • 46

    TibbettM. (2000). Roots, foraging and the exploitation of soil nutrient patches: The role of mycorrhizal symbiosis. Funct. Ecol.14, 397399. doi: 10.1046/j.1365-2435.2000.00417.x

  • 47

    UngerS.FriedeM.HundackerJ.VolkmarK.BeyschlagW. (2016). Allocation trade-off between root and mycorrhizal surface defines nitrogen and phosphorus relations in 13 grassland species. Plant Soil407, 279292. doi: 10.1007/s11104-016-2994-y

  • 48

    WangG.LiX.ChristieP.ZhangJ.LiX. (2016). Response of arbuscular mycorrhizal fungi to soil phosphorus patches depends on context. Crop Pasture Sci.67, 11161125. doi: 10.1071/CP16168

  • 49

    WeiserM.KoubekT.HerbenT. (2016). Root foraging performance and life-history traits. Front. Plant Sci.7, 779. doi: 10.3389/fpls.2016.00779

  • 50

    WickhamH. (2016). Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. Available online at: https://ggplot2.tidyverse.org.

  • 51

    WijesingheD. K.JohnE. A.BeurskensS.HutchingsM. J. (2001). Root system size and precision in nutrient foraging: responses to spatial pattern of nutrient supply in six herbaceous species. J. Ecol.89, 972983. doi: 10.1046/j.0022-0477.2001.00618.x

  • 52

    WilkeC.FoxS. J.BatesT.ManaloK.LangB.BarrettM.et al. (2024). wilkelab/cowplot: 1.1.3. Zenodo. doi: 10.5281/zenodo.10553544

  • 53

    ZhangJ.WangP.XueK.HaoY.WangY.CuiX. (2019). Trait complementarity between fine roots of Stipa purpurea and their associated arbuscular mycorrhizal fungi along a precipitation gradient in Tibetan alpine steppe. J. Mountain Sci.16, 542547. doi: 10.1007/s11629-018-5311-9

Summary

Keywords

mycorrhizal symbiosis, nutrient availability, root foraging precision, root–fungal cooperation, soil heterogeneity, herbs, quantitative real-time PCR (qPCR)

Citation

Stiblíková P, Jansa J, Rozmoš M, Kotianová M, Brindzák M, Šašek J and Weiser M (2026) Root-fungal interactions upon resource foraging in heterogeneous soils. Front. Plant Sci. 17:1899078. doi: 10.3389/fpls.2026.1899078

Received

03 June 2026

Revised

07 August 2026

Accepted

18 August 2026

Published

31 August 2026

Volume

17 - 2026

Edited by

Julio Javier Diez Casero, University of Valladolid, Spain

Reviewed by

Khuram Shehzad Khan, China Agricultural University, China

Shuang Liang, Zhejiang Academy of Agricultural Sciences, China

Updates

Copyright

*Correspondence: Pavlína Stiblíková,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics