Research ArticlePLANT SCIENCES

The fungal collaboration gradient dominates the root economics space in plants

See allHide authors and affiliations

Science Advances  01 Jul 2020:
Vol. 6, no. 27, eaba3756
DOI: 10.1126/sciadv.aba3756


Plant economics run on carbon and nutrients instead of money. Leaf strategies aboveground span an economic spectrum from “live fast and die young” to “slow and steady,” but the economy defined by root strategies belowground remains unclear. Here, we take a holistic view of the belowground economy and show that root-mycorrhizal collaboration can short circuit a one-dimensional economic spectrum, providing an entire space of economic possibilities. Root trait data from 1810 species across the globe confirm a classical fast-slow “conservation” gradient but show that most variation is explained by an orthogonal “collaboration” gradient, ranging from “do-it-yourself” resource uptake to “outsourcing” of resource uptake to mycorrhizal fungi. This broadened “root economics space” provides a solid foundation for predictive understanding of belowground responses to changing environmental conditions.


The diversity of plant traits across the globe shapes ecosystem functioning (1). Seeking general patterns, ecologists have used economic theory to explain trait variation in leaves as the aboveground plant organs for resource acquisition by photosynthesis (13). Aboveground plant strategies thereby fall along a “leaf economics spectrum” (2) from cheaply constructed but short-lived leaves optimized for “fast” resource acquisition to more expensive but persistent leaves with a “slower” rate of return over longer time scale.

Fine roots acquire resources from the soil and are often considered the belowground equivalent of leaves (4). Therefore, fine-root trait variation has been hypothesized to follow a similar one-dimensional spectrum (1, 5). At one side of this spectrum, plants with a “fast” belowground resource acquisition strategy are expected to construct long, narrow-diameter roots with minimal biomass investment but high metabolic rates (1, 4, 6). At the opposite side of the spectrum, plants with a “slow” strategy are expected to achieve longer life span and prolonged return on investment by constructing thicker-diameter, denser roots (4, 7).

However, mixed empirical results caused ecologists to question whether variation in root traits can be adequately explained by a one-dimensional “fast-slow” economics spectrum (1, 5, 810). Instead, when taken together, earlier results reveal that root trait variation might be driven by multiple evolutionary pressures (8, 1114). Here, we aim to settle this debate by presenting a new conceptual framework of root economics that better captures the complexity of belowground resource acquisition strategies. First, we integrated existing knowledge to build a conceptual understanding of the covariation among four key root traits (Table 1 and Fig. 1). Second, we tested our conceptual model against root traits of 1810 plant species from a global database. Third, we investigated generality of the concept across all biomes of the world, different plant growth forms, and symbiotic partnerships. All analyses were phylogenetically informed using fine-root trait data from the Global Root Trait (GRooT) database (15).

Table 1 Rationale of the conceptual framework of root trait correlations depicted in Fig. 1.

Expected correlations are based on mathematical and ecological rationale and empirical support from the literature. De facto correlations (see also fig. S1) are phylogenetically informed correlation coefficients of species subsets with the respective trait coverage. D, root diameter; SRL, specific root length; RTD, root tissue density; N, root nitrogen content; CF, cortex fraction.

View this table:
Fig. 1 Conceptual framework of the root economics space.

On the basis of this concept, we hypothesize (i) a collaboration gradient ranging from do-it-yourself soil exploration by high specific root length (SRL) to outsourcing by investing carbon into the mycorrhizal partner and hence extraradical hypheae, which requires a large cortex fraction (CF) and root diameter (D) and (ii) a conservation gradient ranging from roots with high root tissue density (RTD) that show a slow resource return on investment but are long-lived and well-protected, to fast roots with a high nitrogen content (N) and metabolic rate for fast resource return on investment but a short life span. Arrows indicate negative correlations between the single traits (see Table 1).

The currency of root economics is the carbon input required to construct fine roots that explore the soil for resource acquisition. Specific root length (SRL)—the root length per unit mass—therefore reflects the rate of return per unit of investment and is a function of both root diameter (D) and root tissue density (RTD)—the root mass per unit of root volume—followingSRL=4/(π×D2×RTD)Although this equation (6) is a simplification when sampling heterogeneous fine-root populations (16), it implies that SRL increases with decreasing D and/or RTD. Besides efficient soil exploration, plants have to maintain a high metabolic rate to assure fast resource acquisition leading to high nitrogen content (N) in the fine roots (1, 17). While strong negative relationships between SRL and D (9, 14, 1820) and between RTD and N (9, 14, 19) have been observed, the relationships between SRL and RTD (19, 21, 22) and between D and N (10) have been less clear. In fact, observations across a wide range of species suggest that plants can construct roots with many combinations of SRL and RTD (9, 14), indicating complex trait interactions inconsistent with a one-dimensional root economics spectrum (811, 14).

A growing body of literature (8, 1114) indicates that this root trait complexity may result from the range of belowground resource uptake strategies. In contrast to aboveground photosynthesis, which is solely conducted by plant organs, belowground many species have the ability to outsource resource acquisition. This gradient of plant collaboration strategies ranges from “do-it-yourself” acquisition by cheap roots for efficient soil exploration to “outsourcing” acquisition via the investment of carbon in a mycorrhizal fungal partner for the return of limiting resources. However, these outsourcing strategies have consequences for root traits. This is particularly true for arbuscular mycorrhizal fungi (AMF) because plants must increase their root cortical area, and hence their D, to provide the intraradical habitat for their fungal partner (19, 23, 24). This is generalizable for plant symbiosis with AMF, the most widespread type of mycorrhizal fungi (24) and also well documented for ectomycorrhizal (EM) fungi (25). Using this body of literature as our foundation, we developed an overarching concept of root economics based on the understanding that plants can optimize resource uptake by investing carbon either in thin roots that efficiently explore the soil themselves (14) or in a mycorrhizal fungal partner, which often requires a thick root for efficient symbiosis (Fig. 1).

This conceptualized collaboration gradient from do-it-yourself to outsourcing challenges the traditional spectrum of root economics that assumes D to increase with RTD for tissue conservation. Both scaling laws and empirical data (22) show that as D increases, root cortex area increases at a faster rate than stele area such that D scales positively with the cortex fraction (CF) (19) [although patterns can vary among growth forms (10)]. The parenchymatous cortical tissue has a lower carbon content and dry weight than the stele tissue, which transports nutrients and water through lignified cells (26, 27). Thus, CF and RTD will be negatively correlated (Table 1). Furthermore, since D and CF are closely positively correlated and increase in unison with mycorrhizal symbiosis, D should be negatively correlated with RTD. These relationships contradict the assumption of a one-dimensional root economics spectrum, where plants with a slow strategy are expected to construct roots that are both thick and dense, and advocate for a multidimensional space of root trait variation.


By testing pairwise correlations of all traits, we confirmed the bivariate relationships underlying our new concept of a belowground economics trait space with two main dimensions (Table 1). The strongest negative correlation was found between SRL and D (R = −0.70), representing the “collaboration” gradient from do-it-yourself to outsourcing. We also found a negative correlation between RTD and root N (R = −0.26) as observed in previous studies (9, 14, 19), which corresponds to a “conservation” gradient, representing the traditional trade-off between fast and slow return on investment (Fig. 1).

On a subset of 748 species with complete information on the four main root traits (SRL, D, RTD, and root N), we could confirm these two distinct and largely independent gradients in a phylogenetically informed principal component analysis (PCA) where the first two axes encompass a plane with a cumulative explanatory power of 77% of all root trait variation. Henceforth, we refer to these gradients as the main dimensions of the “root economics space” (Fig. 2A). The first PCA axis (44% of total trait variation) represents a gradient from SRL to D, confirming our conceptualized collaboration gradient and suggesting that it actually represents the main source of root trait variation. The second PCA axis (33% of total trait variation) represents the conservation gradient from root N to RTD (table S1).

Fig. 2 The root economics space.

Phylogenetically informed principal component analyses (PCAs) of core traits of (A) 748 global species, (B) 621 arbuscular mycorrhizal (AM) species (blue), and (C) 94 ectomycorrhizal (EM) species (red). NM, nonmycorrhizal. The collaboration gradient (44%) ranges from do-it-yourself roots with high SRL to outsourcing roots with “thick diameters” (D). The conservation gradient (33%) ranges from fast (N) to slow (RTD). For each corner of the root economics space, in (A) we highlight two representative plant species: QV, Quercus virginiana Mill.; CH, Carex humilis Leyss.; CO, Cornus officinalis Siebold & Zucc.; ZM, Zea mays L.; LP, Lathyrus pratensis L.; GB, Ginkgo biloba L.; BL, Betula lenta L.; CP, Cardamine pratensis L. (D) Woody (ocher) and nonwoody (green) species show no distinct pattern within the root economics space (see fig. S4 and table S4). (E) PCA based on bivariate trait relationships. The percentage mycorrhizal colonization (%M) and the CF are positively correlated with D along the collaboration gradient, while root life span is negatively correlated with N along the conservation gradient. See table S1 for PCA scores.

Species associated with AMF were the largest group in the database and were distributed over the entire trait space (Fig. 2A) but differed significantly from both nonmycorrhizal (NM) and EM species (table S4). NM plants aggregated on the do-it-yourself side of the collaboration gradient and on the slow side of the conservation gradient. EM plants showed less variation along the collaboration gradient than AM plants with a tendency toward do-it-yourself and slow as well. A high RTD, indicative of a slow strategy, might partly originate from the fact that EM species are often woody (28) or it might reflect a general slow nutrient cycling in ecosystems dominated by EM species (29, 30). The tendency for EM plants toward do-it-yourself roots with high SRL likely results not only from the nature of the EM symbiosis that is less dependent on cortex area but also from its more recent evolution, as evolutionarily younger species tend to have thinner roots (14, 23, 27, 31). Even so, PCAs that solely represent the root traits of either AM or EM plant species (Fig. 2, B and C, and table S1) show the same dimensions of variation as in the global dataset, highlighting the existence of the same trade-offs within each mycorrhizal type.

Plant species associated with N2-fixing bacteria differed from other species (table S4) by being located on the fast side of the conservation gradient as their roots are rich in N (fig. S2A). Nevertheless, we could still confirm the collaboration gradient as the first PCA axis within this species set (fig. S2, B and C, and table S1). The importance of the collaboration gradient within N2-fixing species might derive from the large P demands of plants associated with N2-fixing bacteria, which leads to either high mycorrhizal dependency or alternative do-it-yourself strategies like cluster roots (24). Investigating different plant growth forms, we found that woody plants span a wider range of variation than nonwoody plants within the global trait space (Fig. 2 and table S4). Still, the two gradients of the root economics space exist within both woody and nonwoody plants (fig S4), indicating that there is wide variation and very similar trade-offs operating irrespective of growth form. Last, the two dimensions of the root economics space are present irrespective of biome (Fig. 3 and table S1), which did not differ from each other in their location within the global trait space (table S4).

Fig. 3 The root economics space is present in different biomes.

Root traits and trait relations are known to vary across biomes (14). We found no respective between group variation within the root economics space (table S4). Still, to test whether the concept is broadly generalizable, we present separate PCAs for biomes spanning arid to tropical. We found that the root economics space was apparent in all of the biomes represented by our species (panels A, B, C, and D). In continental systems, the conservation gradient was represented by principal component 3 (D) instead of principal component 2 (E). See table S1 for principal component analyses. pc, principal component.

To confirm our ecological interpretation of the proposed gradients, we added traits to the PCA that act as proxies for ecological functions (Fig. 2E and table S2). We used percent root length colonized by AMF (%M) as a proxy for the strength of the mycorrhizal symbiosis (32) and CF as a general proxy for the ability of a species to host mycorrhizal fungi (19, 33, 34). We found both %M and CF to be associated with the outsourcing side of the collaboration gradient. To test whether the proposed conservation gradient aligns with the classical fast-slow economics spectrum, we used root life span as a proxy for short- or long-term investment of plant carbon (1, 3537). We found that longer life span was indeed associated with the slow side of the conservation gradient, which is consistent with reports of negative relationships between root life span and N (1, 35, 37).

The decrease in D over evolutionary time (14, 31) suggests a reduced dependence of plants on mycorrhizal fungi. We found that the collaboration gradient was indeed phylogenetically conserved, showing an evolutionary transition from outsourcing to do-it-yourself (Fig. 4 and tables S3 and S5). In contrast, the fast-slow trade-off of the conservation gradient was less pronounced across all plant families in our database (Fig. 4) and also less phylogenetically conserved (table S3). Terrestrial plants coevolved with AMF, causing the mycorrhizal symbiosis to be evolutionarily stable (38, 39). This might explain the finding that the consequences for root morphology and anatomy—being associated with the collaboration gradient—are phylogenetically conserved at high levels. Different explanations have been proposed as to why D gradually decreased with evolutionary time, including a decline in atmospheric CO2 (13), leading to higher water demands and a reduction in the dependence on mycorrhizal fungi (14). Remarkably, this trend very rarely resulted in a complete loss of the mycorrhizal symbiosis but instead led to varying degrees of outsourcing. Following this line of reasoning, evolutionary history might be the reason why the collaboration gradient is the main source of variation in root traits. As the ability to outsource is a major difference between above- and belowground economics, the importance of the collaboration gradient might be the key to explaining decoupling of root and leaf traits, leading to inconsistencies within the plant economics spectrum found in the past (9, 11, 27, 40).

Fig. 4 The collaboration gradient is phylogenetically conserved.

On the left, we display the phylogenetic tree of 1810 species aggregated at a family level with the standardized family mean trait values of the four core traits (center) ranging from low (yellow) to medium (green) to high (blue). The collaboration gradient shows a strong phylogenetic pattern (λ = 0.8, P < 0.001) with a transition from families with thick D to those with high SRL. The phylogenetic signal in the conservation gradient is less pronounced (λ = 0.5, P < 0.001), although still significant (see also table S3). For detailed information about specific clades, see table S5, and for family distribution across clades, see table S6. Pie charts (right) depict the fraction of different mycorrhizal association types within the broader plant phylogenetic clades (indicated by corresponding background colors).


Together, we provide a conceptual framework explaining the mechanistic basis behind root trait covariation and show its ubiquity across biomes, growth forms, and symbiotic partnerships. The root economics space synthesizes recent evidence to illustrate why root trait variation cannot be adequately explained by a one-dimensional spectrum (8, 9, 11, 14, 19, 41). Plant outsourcing of belowground resource acquisition through collaboration with mycorrhizal fungal partners is not just an “extra” dimension in the root economics space but rather is the main dimension of root trait variation, which is fundamentally different from the aboveground economy. This collaboration gradient from do-it-yourself to outsourcing represents an investment in soil exploration by either the root itself or its mycorrhizal fungal partners. It is independent from the conservation gradient, which represents the well-known concept of fast versus slow return on investment. Thus, both gradients depict different facets of root economics and, rather than a single one-dimensional spectrum, encompass a whole root economics space of plant strategies for belowground resource acquisition.



All analyses presented here are based on the GRooT database (15). The GRooT database combines root trait observations from the Fine-Root Ecology Database (FRED) (42) and TRY (43) with additional datasets providing data measured on individual plants for which taxonomical information is available. It includes data on both coarse and fine roots. For the objective of this study, we selected fine roots only, as coarse roots are usually not absorptive and therefore less relevant in the context of root economics (42, 44). We treated roots as fine roots if they met at least one of the following criteria: (i) they were of root orders 1 to 3, (ii) they were classified as “fine roots” by the initial authors, or (iii) their diameter was smaller than 2 mm. Data measured on dead roots were excluded from the analyses. Furthermore, we excluded ferns (Polypodiopsida) because of their very special root morphology that is hardly comparable with vascular plants (45, 46). We only selected data where species-level information was available. During the past decade, a set of root traits was found to be highly informative of root economics: SRL, D, RTD, and root N (811, 14, 20, 27). Hence, we focused our main analyses on those four traits. In addition, we analyzed the percentage of root length colonized by mycorrhizal fungi (%M) – while > 99 % of these data refer to arbuscular mycorrhizal colonization - and the root area occupied by the cortex, i.e., CF as proxies for the strength of mycorrhizal symbiosis, as well as mean root life span. We checked the values of these traits for outliers and excluded values of RTD exceeding 1.0 in further analyses.

Categorical data from GRooT such as main biome type (tropical, temperate, continental, arid, or polar) following the Köppen-Geiger classification, plant woodiness (woody, nonwoody, or facultatively woody), mycorrhizal association (NM, arbuscular mycorrhiza, ectomycorrhiza, or other (e.g., ericoid mycorrhiza), and nitrogen-fixing ability (fixers or nonfixers) were used in the downstream analysis and testing of our conceptual framework. GRooT includes mycorrhizal association data from FungalRoot (47), which did not cover our entire species set. To achieve full data cover, we filled the gaps and did minor annotations based on the following general rules:

(i)Mycorrhizal association is constant within species, hence excluding a facultative mycorrhizal type. In cases where a lack of mycorrhizal colonization has been reported only under specific environmental conditions less suitable for the mycorrhizal symbiosis, we assigned the species as mycorrhizal, while in cases with intraradical hyphae but no evidence for a symbiotic interface, we assigned species to be NM.

(ii)Almost all plants have one type of mycorrhizal association as the dominant one. Therefore, dual mycorrhizal association was only assigned if species show no clear dominance toward one type.

(iii)The mycorrhizal association type is usually constant within a monophyletic genus and often within a family (24, 47, 48). Therefore, we filled remaining gaps with the respective mycorrhizal association type of sister species.

Data processing

All data processing and analyses were done using R 3.6.1 (49). In this study, we analyzed how different root traits are related to each other at the level of plant species; hence, a first step was to calculate species mean values. As root traits were measured in different studies varying in design (e.g., on in situ grown plants versus plants growing in pots), and because most traits varied several orders of magnitude, several steps of data processing were required before calculating species mean trait values. First, to obtain normal distributions, we log-transformed each trait, except for %M and CF, which were scaled to the range of 0 to 1 and arcsine square root transformed. We then Z transformed each trait to a mean of 0 and an SD of 1 to assure variance homogeneity. Furthermore, we corrected for main study design (measurements on plants in situ, in pots, or hydroponics) and the publication in which the trait measurements were first reported (as a proxy for other study specific factors, e.g., plant age, soil conditions, or sample handling). This was done by building a linear mixed model for each trait, where the trait was treated as the response variable, study design as a fixed factor, and publication as random factor. We used residuals of these models in further analyses.

Within some species, the categorical traits woodiness and biome had different data entries (e.g., because the species occurs both in temperate and continental biomes). In those cases, we categorized the species in the biome in which it had most observations, and we categorized its woodiness by its most commonly observed entry in further analyses.

In total, we analyzed information of 1547, 1662, 1361, and 1158 species for D, SRL, RTD, and N, respectively. Scientific names in GRooT are standardized among datasets and brought up to date by querying species names using the Taxonomic Name Resolution Service v4.0 ( (50). We constructed a phylogenetic tree including all species using the backbone phylogeny from Zanne et al. (51) and adding additional missing species with the function “” from the package “phangorn” (52). We calculated Pagel’s λ using the package “picante” and evaluated the strength of the phylogenetic signal for each trait; a large (close to the upper bound of 1.0) Pagel’s λ value indicates higher phylogenetic conservatism (53), whereas a low (close to 0.0) value indicates a lack of phylogenetic conservatism.


As all traits exerted strong phylogenetic signal (table S3), we used phylogenetically informed methods for all analyses. We first assessed bivariate relationships between the four core traits (D, SRL, RTD, and N) and CF to build our conceptual framework (Table 1 and Fig. 1), and we also tested for relationships of these traits with %M and root life span (fig. S1). Sample sizes varied for these bivariate correlations, depending on the number of species with complete information for both involved traits, and ranged from 19 (for the correlation between %M and root life span) to 1402 (for the correlation between D and SRL) (fig. S1). In total, we used 1810 species for these bivariate correlations. We fit phylogenetic generalized least square models using the “pgls” function in the R package “caper” (54, 55) to each pair of traits to conduct phylogenetically corrected regression analyses. Phylogenetically corrected correlation coefficients (r values) were then calculated by taking the square root of the adjusted model r2 and by multiplying this with −1 if the regression coefficient was negative. In cases where the adjusted model r2 was negative, we assigned an r coefficient of 0.

We used phylogenetically informed PCA to identify main dimensions of variation among root economic traits. A phylogenetic PCA was performed for the four core traits D, SRL, RTD, and N using the “phyl.pca” function of the “phytools” package (56). There were 748 species that had complete data for these four core traits. The eigenanalysis uses the correlation structure of the phylogeny to inform its estimates of eigenvalues and eigenvectors (56, 57). To assess whether the PCA results and hence the dimensions of the root economics space were sensitive to biome type, mycorrhizal association, woodiness, or nitrogen-fixing ability, we repeated the above analysis for subsets of different biomes [tropics, temperate, continental, and arid—the polar biome was represented by too few species (n = 5) to perform a reliable analysis], mycorrhizal association type (arbuscular mycorrhiza versus ectomycorrhiza), woodiness (woody versus nonwoody), and nitrogen-fixing ability (present or absent).

We assessed whether roots from species with different mycorrhizal associations (arbuscular mycorrhiza, ectomycorrhiza, arbuscular mycorrhiza, and ectomycorrhiza, i.e., intraspecific variation in mycorrhizal association type, ericoid mycorrhiza, and nonmycorrhiza associated), species from different biomes (temperate, tropical, arid, and continental), woody or nonwoody species, and species that either did or did not associate with bacteria able to fix nitrogen, differed significantly from each other in the global multidimensional PCA space, i.e., in their first two PCA axes (which jointly explained 77% of all trait variation). This was done using a permutational multiple analysis of variance, in which the first two PCA axes were treated as the response variables and mycorrhizal association type, biome, woodiness, or ability to fix nitrogen as the fixed factor. We used Euclidean pairwise distances in PCA space among species and calculated 999 permutations using the “pairwise.adonis” function in the “pairwiseAdonis” package (58). To test for the significance of differences between different categories of mycorrhizal associations and biomes, we used false discovery rates (59) to reduce the likelihood of type I errors due to multiple testing.

Furthermore, we investigated multivariate trait space for seven traits, i.e., the four core traits from the above-described PCAs (D, SRL, RTD, and N) supplemented by three additional traits: CF, %M, and root life span. Observation-based PCA requires each replicate (species) to have complete data for all traits, but there were few species with a complete set of all seven traits. Therefore, we performed an alternative dimensionality reduction analysis based on pairwise correlations between traits. For this analysis, we computed phylogenetically informed pairwise correlations between each of the 21 trait combinations, where each trait combination included a slightly different set of species for which traits were available (fig. S1). We then performed a standard eigenanalysis on this positive definite matrix of phylogenetic correlation coefficients (60).

Visual examination of the distribution of traits across the phylogeny was obtained using the function “phylo.heatmap” in the package “phytools” (56). We further examined the phylogenetic trends observed across broader phylogenetic clades of seed plants using a randomization test to quantitatively compare individual clade trait values to the rest of the phylogeny. The test determines whether the mean trait value observed in a clade deviates significantly from the population mean under the null hypothesis that the trait has a random phylogenetic distribution. To do so, we created an algorithm in R that selected clades sequentially at each node. Because of the large number of species, we selected particular nodes that enveloped important phylogenetic clades with at least 30 species (tree tips) included. For each clade, we calculated the observed mean and kurtosis values as measures of central tendency and dispersion values within clades, respectively. Then, we generated a series of 999 random values shuffling trait values among the tips of the original tree. Significance was calculated after estimating if the observed clade mean or kurtosis were outside the 95% confidence intervals of the clade estimations using the randomized datasets. In the case of the kurtosis, values higher than the randomized mean were interpreted as evidence of underdispersion in the clade (leptokurtic distribution), whereas lower values were considered sign of overdispersion (platykurtic distribution).


Supplementary material for this article is available at



Acknowledgments: We would like to thank I. Mansour for text editing. Funding: This paper is a joint effort of the working group sROOT supported by sDiv, the Synthesis Centre of the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, funded by the German Research Foundation (FZT 118). The sROOT workshops and L.M. were also supported by NWO-Vidi grant 864.14.006. J.B. was supported by DFG grants RI-1815/20-1 and RI 1815/22-1. C.M.I., M.L.M., and FRED were supported by the U.S. Department of Energy’s Office of Science, Biological and Environmental Research Program. Author contributions: J.B., A.W., and L.M. conceived the idea for the project. All authors were involved in collecting datasets, developing the conceptual framework, and interpreting the results. J.B., F.v.d.P., D.C.L., N.G.-R., O.J.V.-B., and L.M.Y. performed the statistical analyses. T.W.K. annotated the mycorrhizal associations. J.B., A.W., C.M.I., and L.M. wrote the first draft of the manuscript. All authors commented on and agreed with the final version of the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. All data analyzed in the study originate from the GRooT database (15), which is publicly available on GitHub ( The R script including all analyses and figure preparations is available from the corresponding author upon reasonable request.

Stay Connected to Science Advances

Navigate This Article