Multiple phosphorus acquisition strategies adopted by fine roots in low-fertility soils in Central Amazonia

Ancient Amazon soils are characterised by low concentrations of soil phosphorus (P). Therefore, it is hypothesised that plants may invest a substantial proportion of their resources belowground to adjust their P-uptake strategies, including root morphological, physiological (phosphatase enzyme activities) and biotic (arbuscular mycorrhizal (AM) associations) adaptations. Since these strategies are energy demanding, we hypothesise that trade-offs between morphological traits and root phosphatase exudation and symbiotic associations would occur. Specifically, we expected that plants which invest in finer roots, and therefore have greater ability to explore large soil volumes, would have a high investment in physiological adaptations such as enhanced phosphatase production. In contrast, we expected that plants with predominantly thicker roots would invest more in symbiotic associations, in which carbon is traded for P acquired from AM fungal communities. We collected absorptive roots (<2 mm diameter) from a lowland Central Amazon forest near Manaus, Brazil. We measured fine root diameter, specific root length (SRL), specific root area (SRA), root tissue density (RTD), root phosphatase activity (APase) and arbuscular mycorrhizal (AM) fungi colonisation. Root morphological traits were related to APase activity, with higher APase activity in roots with higher SRL and SRA but lower RTD. However, the degree of AM colonisation was not related to any measured root morphological trait. Fine absorptive roots likely benefit from having low RTD, high SRL, SRA and APase exudation to acquire P efficiently. However, because AM colonisation was not related to root morphology, we suggest that investment in multiple P-uptake strategies is required for maintaining productivity in Central Amazon forests.


Introduction
Amazon forests are highly productive ecosystems, crucial for global carbon (C), nutrient and water cycling (Malhi et al. 2004) and C storage (Pan et al. 2011). Yet, about 60% of the soils in the Amazon Basin are characterised by concentrations of phosphorus (P) and cations that are considered limiting to plant growth (Quesada et al. 2011). Phosphorus is used in an array of plant processes, including photosynthesis and respiration (Raghothama and Karthikeyan 2005;Vance et al. 2003), and, as a result, P availability in soils is usually considered to play an important role controlling growth rates in Amazon forests (Quesada et al. 2010(Quesada et al. , 2012. In young ecosystems, phosphorus is supplied to soil almost exclusively by the weathering of the parent material (Walker and Syers 1976), with small contributions from atmospheric deposition, varying greatly in space and time (Yu et al. 2015). Generally, during soil development, the importance of parent material as a source of nutrients declines considerably over time; total soil P concentration decreases due to weathering and leaching, and the remaining P pools become gradually more recalcitrant, adsorbed or occluded within mineral matrices and organic matter (Quesada et al. 2010;Tiessen et al. 1984;Walker and Syers 1976). Hence, over geological time scales, tropical forests become essentially dependent on the recycling of organic P from litter decomposition.
The concept of P limitation in tropical forests is widely accepted, but still, the mechanisms that allow Amazon forests to thrive in such P-poor environments are not clear. Plants can respond in two ways to overcome P limitation in soils: 1) increasing P-use efficiency aboveground, and/or 2) by improving P-uptake belowground (Kitayama 2013;Vitousek et al. 2010). Plants can enhance aboveground P-use efficiency by increasing the residence time of P within the canopy (Kitayama 2013;Kitayama and Aiba 2002;Kitayama et al. 2000;Vitousek 2004), maximising P resorption in leaves before senescence (Aerts 1996;Hayes et al. 2014;Raghothama and Karthikeyan 2005) or reducing or replacing P in essential compounds with other elements White and Hammond 2008). In addition to aboveground phosphorus efficient use, plants under low-fertility soils may be able to alter their root morphology and physiology, and/or vary their reliance on mycorrhizal fungi to optimise nutrient acquisition (Raghothama 1999;Richardson and Simpson 2011). Although belowground plant nutrient uptake mechanisms could be important contributors to tropical forest productivity (Bardgett et al. 2014), there is limited understanding of the role of different root traits in Puptake in tropical forests, especially in the Amazon.
Phosphorus is very immobile in soils, and for this reason roots preferentially move towards P (Aerts 1999). Root morphological traits, such as total root length, surface area and diameter directly control the efficiency by which P is assimilated by plants (Eissenstat 1992;Eissenstat and Yanai 1997), with small diameter roots being usually efficient at exploring large soil volumes (Bates and Lynch 2001;Hodge 2004;Liu et al. 2015). However, the only form of P assimilated by plants and microorganisms is inorganic phosphate (mainly PO 4 −3 ), which is largely bound in organic compounds within tropical soils (Walker and Syers 1976). Consequently, these compounds cannot be directly taken up by plants and microorganisms and need to be degraded before being assimilated by roots (Lambers et al. 2006). The exudation of phosphatase enzymes (Hinsinger 2001;Olander and Vitousek 2000) and carboxylates (Jones and Oburger 2011;Lambers et al. 2015;Turner et al. 2012) are some of the mechanisms by which roots mobilise organic and inorganic occluded forms of P in soils. The main process through which organic P is hydrolysed in soils is via the activity of phosphatase enzymes released by microbes and plant roots (Hinsinger 2001;Olander and Vitousek 2000;Treseder and Vitousek 2001;Vance et al. 2003). Therefore, the coordination between small root diameter and large SRL and SRA together with phosphatase exudation may give fine roots a competitive advantage over thick roots by enabling them to explore the soil volume and mobilise organic P more efficiently (Kong et al. 2016;Liu et al. 2015).
Species with thick absorptive root systems, on the other hand, are not very efficient at nutrient foraging and are expected to preferentially invest in mycorrhizal associations to meet nutrient demands (Comas et al. 2014;Eissenstat et al. 2015;Kong et al. 2016). The very fine hyphal network typical of arbuscular mycorrhizas (AM) allow the fungi to acquire P away from P-depleted zones around roots, resulting in high inorganic P-uptake in exchange for photosynthetically fixed C from the host plant Hodge 2004;Smith and Read 2010). In addition to their role in P uptake, AM fungi may also benefit plants by altering the microbial community in the rhizosphere, changing soil aggregation and/or decreasing plant susceptibility to pathogens (Herre et al. 2007;Koide 1991;Laliberte et al. 2015). Because of the multiple benefits that AM fungi confer to plants, this symbiotic association is extremely common and it is estimated that about 74% of all plants species are colonised by these fungi (Brundrett 2009;Smith and Read 2010).
Root morphological adaptations, as well as investment in phosphatase enzyme production and association with mycorrhizal fungi are, however, resource-costly (Jakobsen and Rosendahl 1990;McCormack et al. 2015;Treseder and Vitousek 2001). For instance, the construction of fine acquisitive roots might require less C per unit root length but these roots also display higher rates of respiration and faster turnover when compared to coarser roots McCormack et al. 2015). Mycorrhizal fungi can receive up to 25% of C assimilated by plants (Jakobsen and Rosendahl 1990) and root exudates can represent up to half of belowground C allocation (Lynch and Ho 2005). As a result, it has been suggested that many plants do not invest energy in multiple nutrient-acquisition mechanisms, but rather show a trade-off between one strategy at the expense of another (Nasto et al. 2017;Ryan et al. 2012;Ushio et al. 2015). One such trade-off in nutrient foraging strategies between thin and thick absorptive roots has been shown for a few plant species (Chen et al. 2016;Ushio et al. 2015) but the generality of such a mechanism at the ecosystem-level has yet to be tested.
Since the role of fine roots in forest functioning in Amazonia is poorly understood, this study aims to describe root strategies used by plants to overcome Plimitation in Central Amazonia, testing for possible trade-offs between alternative strategies. To address this we investigated root morphological traits represented by root diameter (mm), specific root length (cm mg −1 ), specific root area (cm 2 mg −1 ) and root tissue density (mg cm −3 ), as well as root phosphatase activity and association with arbuscular mycorrhizal fungi. Our hypothesis was that under similar environmental conditions, a trade-off between resource investment in root morphological traits and P-uptake strategies would occur, with different predominant P-uptake strategies being associated with roots of different morphologies. We would therefore expect that root phosphatase activity (APase) should be negatively correlated with root diameter and RTD, but positively correlated with SRL and SRA and that AM colonisation should correlate positively with root diameter and RTD and negatively with SRL and SRA. Under an alternative scenario, we hypothesised the lack of such a trade-off because of the very low P availability in soils and/or the many potential roles that AM fungi could play, resulting in roots investing in multiple P-uptake mechanisms independently of their morphology.

Study site description
This study was carried out in the Cuieiras Reserve at ZF2, ca. 60 km north of Manaus, Amazonas, Brazil, maintained by the National Institute of Amazonian Research (INPA). The vegetation is an old growth, terra firme lowland forest with mean air temperature of 26°C and mean annual precipitation around 2400 mm (Araújo 2002). The plots at the study site are characterised by high species diversity and the soils are classified as Ferralsols (World Reference Base for Soil Resources, 2006), with soil pH (in H 2 O) ranging from 4 to 4.7. The soils in our study area show particularly low P concentrations in relation to other soils across the Amazon basin (Quesada et al. 2010(Quesada et al. , 2011, being on the low range of fertility for Amazon soils. For instance, total P in the top 10 cm of soil along the basin varies from 29 to 968 mg kg −1 (Quesada et al. 2010), with total P in our study site varying from 118 to 217 mg kg −1 in the same soil horizon (Table 1 and Table S1, online resources; methodology as described in the online resources).
We sampled 32 plots, measuring 40 m × 40 m, established in areas with similar soil, vegetation and terrain (all plateaus), with plots at least 50 m away from each other. The similar environmental conditions and consistently low-P concentrations in soils across the plots allowed us to investigate the extent to which, under similar soil conditions, different P-uptake strategies are used by plants at the ecosystem level. Due to the small variation in soil fertility among plots in our study site, relationships between soil P concentrations and root traits are only briefly discussed and are shown in Fig.  S1 (online resources).

Root sampling
Freshly produced fine-roots were sampled in May, 2016 using the ingrowth core technique (Metcalfe et al. 2007). In each plot (n = 32), five 12 cm-diameter, 30 cm-deep, root-free ingrowth cores (2 cm plastic mesh bag) were installed in January-February 2016 in the central 20 m × 20 m plot area. Ingrowth cores were collected after 3-4 months of installation and all fine roots were manually extracted and root-free soil reinserted into the existing holes (Metcalfe et al. 2007). Root samples were divided by depth (0-5, 5-10, 10-20 and 20-30 cm) and cleaned by removing soil particles. Fine roots are usually defined as ≤2 mm in diameter, but because we sampled relatively young roots (< 4 months old), our samples represented only first to third order roots, as classified by Pregitzer et al. (2002). We therefore defined these roots by their function as absorptive fine roots, directly involved in nutrient uptake (McCormack et al. 2015). We restricted our analyses for subsamples from the 0-10 cm depth interval, since the majority of roots were found in the top 10 cm soil layer. Five replicate cores were combined per plot (n = 32) and subsamples were analysed for mycorrhizal colonisation and fine root acid phosphatase activity (APase). Another subsample was scanned for morphology analysis and then dried at 60 o C to determine dry root mass. Root productivity was calculated as dry mass of roots (mg) produced per day for the top 10 cm soil layer.

Root morphology
Root samples from the 0-10 cm layer were scanned twice to determine root morphological traits. First an image was obtained with all roots from the sample, subsequently, some root fragments were picked to analyse for mycorrhizal colonisation and the sample was scanned again to generate a second image on which root fragments used for AM analysis were absent. The morphological traits for the fragments used for AM analysis (methodology as described in the online resources) were calculated using the difference between the two scanned images. In addition, another small subsample of root tips Values represent mean (± SD) and ranges for each P fraction (mg kg −1 ) extracted from all plots (n = 32). Pi = inorganic P; Po = organic P. P fractions are shown in decreasing order of availability, with resin P being the most available fraction and residual P being the least available. Total P is the sum of all fractions; readily available P is calculated as the sum of resin P and bicarbonate Pi and Po (fractions that are very to moderately available for plants); total extractable P is calculated as the sum of all fractions except residual P. SD = standard deviation were used for the enzyme assay and scanned separately. All images were scanned at 600 dpi and analysed using WinRHIZO (WinRHIZO Regular, Regent Instruments, Canada) to determine specific root length (SRL), specific root area (SRA), root tissue density (RTD) and mean root diameter (Metcalfe et al. 2008). SRL (cm mg −1 ) was calculated as the length per unit root dry mass; SRA (cm 2 mg −1 ) was calculated as root superficial area per unit dry mass; RTD (mg cm −3 ) was calculated as root dry mass per unit root volume. With these three sets of root morphological measurements we were able to discuss general morphological traits at plot and ecosystem level but also look at correlations between the root morphological traits, AM colonisation and APase activity.

Root phosphatase activity
Potential root acid phosphomonoesterase activity (APase) was measured for roots from 0 to 10 cm soil depth (n = 32) using a fluorimetric microplate assay adapted following protocols by German et al. (2011) and Turner and Romero (2010). All clean root samples were kept refrigerated (4°C) to avoid tissue degradation and enzyme activities were measured within 2 weeks of root sampling using triplicate subsamples per plot. About 10 mg of each root sample (washed, fresh weight basis) were placed into a sterile 2 ml Eppendorf snapcap vial with 1 ml of buffer and 0.25 ml of Methylumbelliferyl-phosphate (MUF-phosphate; 2 mM). A further identical root sample was prepared as a control by adding 1.25 ml of buffer. In addition, buffer blanks (n = 3) and substrate blanks (n = 3) were prepared as 1.25 ml of buffer (no roots, no substrate) and 1 ml of buffer + 0.25 ml of MUF-phosphate (no roots), respectively. Samples were incubated for 30 min at 25°C while gently shaking, then 50 μL of 1 M NaOH were added to all samples and standard vials to terminate the reactions. Aliquots of the sample solution were pipetted into a black 96-well microplate and 20 min post termination, fluorescence was read on a fluorometer (Tecan Infinite® 200 PRO, Grödig, Austria), at 365 nm excitation and 450 nm emission. Roots were then removed from vials, rinsed with Milli-Q water and dried at 60°C for 72 h. Root APase activity per plot is represented in nmol MUF mg −1 dry root mass h −1 . Based on the images from the scanned root subsamples, root APase activity was also calculated per cm root length and cm 2 root area (online resources, Table S2).

Mycorrhizal colonisation
To determine AM colonisation, root samples from the 0-10 cm soil depth were cleaned and segments of fresh absorptive roots from the first three orders were cut and stored in 50% ethanol. Samples were then transferred to small cassettes and if needed, subsamples were separated by colour/appearance: dark roots, red/brown and clear roots. The clearing and staining processes, designed to highlight only the mycorrhizal structures, were adapted to tropical roots based on Brundrett et al. (1984) and Wurzburger and Wright (2015). The cassettes were placed in a 2.5% KOH solution and autoclaved at1 20°C for~10 min, with processing time depending on the darkness of the roots. If roots were still darkly pigmented after the clearing process, they were placed in alkaline H 2 O 2 solution for further bleaching for3 0 min, checking them constantly so no material was lost due to over clearing or bleaching. Before staining, roots were acidified by placing cassettes in 2% HCl solution for 30 min and then were added to a beaker with Trypan Blue 0.05% for 30 min. When consistently blue, roots were rinsed well and stored in distilled water in the fridge (4°C) for slide preparation. The time between root staining and AM colonisation analysis never exceeded 2 weeks. Roots were placed in rows across the length of the slide and fixed using polyvinyl alcohol (PVA). To quantify mycorrhizal colonisation, 50-100 intersections per sample were read using the microscope micrometer (vertical cross-hair) to analyse the sections (McGonigle et al. 1990). At each intersection the colonisation was scored according to the following categories: no mycorrhizal structures, hyphae only, hyphae + arbuscules, hyphae + vesicles, hyphae + vesicles + arbuscules. Mycorrhizal colonisation was assessed as the percentage of the total root intersections along the root length that had mycorrhizal fungi. Only total AM colonisation (sum of colonisation by hyphae, arbuscules and vesicles) is discussed in detail in the results, and colonisation by each AM fungal structure is shown in Table S2 (online resources).

Statistical analyses
Linear models were used to test for the influence of soil P fractions on root traits. Bivariate relationships between root morphology properties, root phosphatase activity and AM colonisation were described using Standard Major Axis (SMA) line fits using the package smatr (Warton et al. 2012;Warton et al. 2006). Standard Major Axis analysis is generally used in plant allometry studies when there is no clear causation among the variables tested. By using the SMA test, we aimed to estimate the line of best fit between traits, taking into account inherent error associated with both axes. All statistical analyses were conducted in R version 3.3.3. (R Core Team 2017).

Root traits variation
Despite the common soil conditions in our study site (Table 1 and Table S1, online resources), we found marked variation in root morphology (Table 2 and  Table S2, online resources), however, no significant relationships between soil P fractions and root traits were observed, with the exception of a positive relationship between bicarbonate Pi (one of the most readily available forms of P) and RTD (Fig. S1, online resources). Correlations between root morphological traits are not discussed here but are shown in Figs. S2-4 in the online resources for each set of morphological data. Fine root diameter varied from 0.39 to 1.1 mm, with RTD displaying a three-fold variation, ranging from 141.78-419.22 mg cm −3 (Table 2 and Table S2, online resources). Specific root area and root APase activity showed a four-fold variation, varying from 0.14-0.56 cm 2 mg −1 and 15-66 nmol mg −1 h −1 respectively. Among all traits measured, SRL was the most variable, ranging from 0.59 to 4.15 cm mg −1 . All sampled roots showed some degree of AM colonisation, ranging between 10 to 80% of root length colonised by AM, with an average of 44.31% (Table 2 and Table S2, online resources).
Relationship between root morphological properties and P-uptake strategies Root morphological traits were found to be related to the levels of root APase activity (Fig. 1). Root APase activity decreased with increasing RTD (r = 0.14, P = 0.035), and increased with increasing SRL (r = 0.16, P = 0.02) and SRA (r = 0.25, P = 0.003), with no significant relationship between root APase and root diameter (Fig. 1a). In contrast, no relationship between total AM colonisation levels and root morphological traits was found, as shown in Fig. 2 (correlations between root morphological traits and AM colonisation by different fungi structure are not discussed here but are shown in Fig. S5).

Trade-offs between root P-uptake strategies
Since root morphological traits were found to influence root APase activity but not AM colonisation, the existence of a trade-off between both root P-uptake strategies was not supported by our results. This is further demonstrated by the lack of a clear relationship between APase activity and AM colonisation ( Fig. S6 and S7, online resources) although it should be noted that it was Values represent means and ranges for each root trait. SRL = specific root length (cm mg −1 ); SRA = specific root area (cm 2 mg −1 ); RTD = root tissue density (mg cm −3 ), APase = acid phosphatase activity (nmol mg −1 h −1 ) and AM colonisation = arbuscular mycorrhizal fungi colonisation (%). SD = standard deviation. Mean root morphological data (diameter, SRL, SRA and RTD) and root productivity are from the whole ingrowth core sample. APase activity and AM colonisation were determined in roots subsampled from the total ingrowth core sample. All measurements, with the exception of AM colonisation and diameter, are shown in root dry weight basis not possible to carry out the enzyme and mycorrhizal assays on the same root subsamples.

Discussion
Fine root traits are often reported to reflect plant efficiency in acquiring the most limiting soil nutrients (Lambers et al. 2006;Raghothama 1999). Based on a presumed trade-off in energy investments between root trait adaptations to improve P-uptake, we aimed to test if in Central Amazon forests plants would also preferentially display one dominant root P-uptake strategy. Despite the similar topography and the relatively narrow variation in P availability in our study site, root traits did not converge towards a single P-uptake strategy but rather displayed multiple complementary strategies to overcome apparent P limitation. We therefore suggest that the variations in root traits found in this Central Amazon forest reflect the dominant features of a much larger diversity of adaptations displayed by different plants in low fertility soils. APase = acid phosphatase activity (nmol mg −1 dry root h −1 ); SRL = specific root length (cm mg −1 ); SRA = specific root area (cm 2 mg −1 ); RTD = root tissue density (mg cm −3 ). Solid lines indicate significant relationships (P < 0.05)

Variation in root morphology
Although we found SRL and SRA varying across plots, our mean values were similar to the ones found by Metcalfe et al. (2008) in a clay-rich Ferralsol in Para, Brazil. Also using the ingrowth core technique, these authors found SRL varying from 1 to 3.5 cm mg −1 and SRA varying from 0.25 to 1.5 cm 2 mg −1 during a yearround sampling campaign, while our mean values for one sampling interval (3-4 months) were 2.09 cm mg −1 and 0.33 cm 2 mg −1 for SRL and SRA, respectively. (2015) reported mean values of 1.5 cm mg −1 for SRL and 150 mg cm −3 for RTD at a site in Panama, but this was based on measurements made on standing stocks rather than in ingrowth cores. Additionally, Kong et al. (2014) also using standing stock cores, found mean SRL and RTD of 5.6 cm mg −1 and 358.5 mg cm −3 respectively for trees at a subtropical site in China, contrasting with 2.09 cm mg −1 and 239.47 mg cm −3 found in this study. The data reported here and from other studies in tropical forests are, however, not directly comparable because of Fig. 2 Standardised major axis (SMA) relationships between total root AM colonisation (% of root length) and root morphology properties. Each dot represent a mean of root traits per plot (0-10 cm soil depth) extracted from ingrowth cores (ecosystem level measurement). SRL = specific root length (cm mg −1 ); SRA = specific root area (cm 2 mg −1 ); RTD = root tissue density (mg cm −3 ) the difference in methods used (ingrowth core versus standing stock biomass). We aimed to investigate young, absorptive roots, making ingrowth cores the most appropriate method for our study, but further research on established root systems in future studies would also be valuable.

Relationship between root morphology and phosphatase activity
We hypothesised that root phosphatase activity would be related to root morphological traits, with higher APase activity in roots with smaller diameter and lower tissue density and therefore higher SRL and SRA. Although our samples refer to plot means at the ecosystem level, we found partial support for the first hypothesis, since three of the four root morphological traits analysed here were correlated as expected with APase activity, with the exception of root diameter. Our results confirm the findings of Ushio et al. (2015) who found lower APase activity with increasing RTD, and higher APase activity with increasing SRA at the species level in tropical montane forests in Borneo. However, Ushio et al. (2015), studying fine roots up to 2 mm diameter, reported that APase activity was negatively correlated with root diameter, whilst no such relationship was found in this study. The lack of a significant relationship between root diameter and APase found here could be due to the small variation in root diameter from our root subsamples (0.29-0.57 mm; Fig. 1a), since only first order acquisitive roots were analysed. A decline in RTD with increasing APase could be seen as a strategy to maximise belowground biomass investment per root volume (Eissenstat 1992) thereby reserving energy for phosphatase exudation and consequently P-uptake. Despite the low C-costs of construction, fine roots have higher costs of maintenance when compared to coarse roots, displaying high rates of respiration and turnover (Eissenstat and Yanai 1997;McCormack et al. 2015). For instance, because of their lower mechanical strength, fine roots are more susceptible than coarse roots to the attack of soil-borne pathogens and herbivores Ushio et al. 2015). The high investment in producing roots with high SRA and SRL associated with high risks of carbon loss could be compensated by high-P acquisition resulting from higher APase activity (Ushio et al. 2015). Because samples in this study represent many different plant species, correlation coefficient values obtained for APase activity and morphological traits were low, but despite the species-specific noise, our analysis still captures the dominant signal on the relationships tested at the ecosystem level. Nevertheless, our results suggest that investing in fine long roots and APase exudation seem to be an efficient P-uptake mechanism, possibly overcoming the potential costs of root maintenance.
The efficiency of plants in acquiring P generally increases under P limitation and therefore, higher root APase activity is usually found in soils with lower P concentrations (Kitayama 2013;Raghothama and Karthikeyan 2005). Root APase activity at our site varied from 15 to 65 nmol mg −1 h −1 , which was lower than the values found in a montane forest in Borneo (100-200 μmol pNPP g −1 h −1 ; note that nmol mg −1 h −1 and μmol g −1 h −1 are equivalent units) (Kitayama 2013;Ushio et al. 2015) but higher than root APase activity found in plants associated with AM in Panama (10 μmol pNPP g −1 h −1 ) (Steidinger et al. 2015) and also higher than in forests in Costa Rica (5-12 μmol 4-MUF g −1 h −1 ) (Nasto et al. 2014(Nasto et al. , 2017. Total soil P concentrations in forests in Borneo were very low, ranging from 10 to 60 mg kg −1 , which could explain the higher investment of plants in root APase activity reported by Kitayama (2013) and Ushio et al. (2015). In addition, litter decay rates and consequently organic P mineralisation is slower in montane forests because of the lower temperatures (Grubb 1977;Vitousek and Sanford 1986), which could also translate into higher phosphatase exudation in Borneo. On the other hand, total P concentrations in soils in Costa Rica (665-1600 mg kg −1 ) (Nasto et al. 2014), were much higher than those found in this study (118-227 mg kg −1 ). Thus, the low root APase activity in Costa Rica and Panama could therefore suggest that (1) Central American forests are not as P-limited as the Central Amazon forests studied here or (2) the greater P availability allows plants in these Central American forests to invest in Puse/acquisition strategies other than APase activity.
We present some of the first root phosphatase data from Central Amazon forests, with assays undertaken on roots that were stored for 2 weeks prior to analysis. While 2 weeks of storage has been used in studies of enzyme activities in tropical soils (Turner and Romero 2010), it would be useful to carry out further investigations into whether storage time or storage conditions affect root phosphatase measurements on tropical roots. Such tests would help us to evaluate and compare our potential activity values further and to develop methodological recommendations to be employed in future research. Nonetheless, our results still represent real differences among plots, once all the samples received the same treatment.

Relationship between root morphology and total mycorrhizal colonisation
We also hypothesised that thick, dense and short absorptive roots would have higher AM colonisation than small diameter roots. Despite substantial variation in root morphological traits, no significant relationships were observed between AM colonisation and the other root traits, irrespective of whether relationships with total mycorrhizal colonisation rates (Fig. 2) or rates of colonisation by different AM structures (arbuscules, vesicles, hyphae; see supplementary material) were investigated. Total AM root colonisation in our site ranged from 10 to 80%, with a mean of 44%, whilst earlier studies have shown that AM colonisation in Brazilian native woody species could range from absent to about 80% in different levels of soil-P (Siqueira and Saggin-Júnior 2001). Apart from soil conditions, root development (length, diameter, root hair density and nutrientuptake capacity), environmental characteristics (soil nutrient availability, plant competition, soil microbial biota, soil depth, seasonality) and fungal hyphae properties (root colonisation rate, growth rate and nutrient-uptake capacity) have all been reported to influence AM colonization (Hoeksema et al. 2010;Johnson 2010;Li et al. 2017;Smith et al. 2011). Evidence for the positive correlation between root diameter and AM colonisation exists for some species in Central Amazon (St John 1980), but it is possible that relationships at the ecosystem level may be different or harder to detect.
Various studies have suggested that fine and thick absorptive roots display different strategies for nutrient acquisition, with thin roots depending on directly acquiring nutrients and thicker roots depending more on mycorrhizal fungi (Chen et al. 2016;Eissenstat et al. 2015;Kong et al. 2016;Liu et al. 2015). The absence of a significant relationship between AM colonisation and root morphological traits found in this study was also reported for Brazilian native woody species studied under greenhouse conditions (Siqueira and Saggin-Júnior 2001) and in a meta-analyses of studies with plants from different biomes (Maherali 2014). Siqueira andSaggin-Júnior (2001) and Maherali (2014) reported certain coarse-root tree species being completely AM-independent, whereas other species with small root diameter showed high AM colonisation. Arbuscular mycorrhizal fungi form associations within root cortical cells only, and for this reason, root diameter does not necessarily represent an appropriate proxy for mycorrhizal colonisation (Brundrett 2002;Comas et al. 2014;Guo et al. 2008). For roots with similar diameters, those with a greater cortex to stele ratio, are usually expected to support higher AM colonisation (Kong et al. 2017;Valverde-Barrantes et al. 2016). Therefore, differences in key root properties between species may obscure relationships between diameter and colonisation rates. Other relevant properties of roots, such as root architecture (e.g. branching frequency, root hair density and length), could also change in response to AM colonisation (Maherali 2014), but were beyond the scope of this study. As an alternative, there is evidence suggesting that plants lack efficient mechanisms to control AM colonization (Johnson 2010;Valverde-Barrantes et al. 2016), meaning that all root cortical tissues created could be potentially colonized independently of P availability. Furthermore, the relative importance of internal AM structures versus extraradical hyphae for nutrient foraging is still not completely understood (Smith and Read 2010), and since we did not sample soil hyphae we assumed that the percentage of internal colonisation was a proxy for P uptake.
Trade-off between root P-uptake strategies and the consequences for ecosystem functioning under P-limitation Investment in nutrient acquisition by plants demands the allocation of other resources, such as C and nitrogen (N) (Nasto et al. 2017). For instance, it has been suggested that mycorrhizal fungi can receive up to 25% of C assimilated by plants (Jakobsen and Rosendahl 1990) and root exudates can represent up to half of belowground C allocation (Lynch and Ho 2005). In addition, the N costs of P-uptake via APase activity could range between 1 to 16 g N per g of P (Treseder and Vitousek 2001). Based on the possible resource costs of investing in multiple P-uptake strategies simultaneously, we hypothesised the existence of a tradeoff between APase activity and AM colonisation (Nasto et al. 2017). However, because only APase activity was related to root morphological traits with no influence on AM colonisation, we can infer that these two strategies were not negatively correlated as we hypothesised, suggesting that at the ecosystem level there is not a trade-off between these P-uptake mechanisms at our study site.
Since P in tropical soils is largely bound in organic compounds or occluded in secondary minerals (Cross and Schlesinger 1995) that need to be hydrolysed prior to plant uptake, we would expect root APase activity to play a major role in P cycling in this Central Amazon forest if compared to the role of AM fungi. The average level of AM root colonisation in our site was 44%, lower than compared to the 70% AM colonisation observed in nutrient-rich soils in Panama (Wurzburger and Wright 2015). However, independent of root morphology, all plots in our study site displayed moderate levels of AM colonisation, suggesting that besides improving inorganic P uptake, AM association could play other roles. Arbuscular mycorrhizal fungi could benefit plants by increasing the uptake of micronutrients (e.g. zinc, copper), change the microbial community in the rhizosphere, increase plant chemical defence against pathogens and herbivores (Herre et al. 2007;Koide 1991;Laliberte et al. 2015) and also possibly by influencing phosphatase production by plants (Joner et al. 2000;Nasto et al. 2014). Fine roots are more susceptible to soil pathogen and herbivore attack  and the potential role of AM increasing plant defence could therefore explain why small diameter roots were found to be colonised by AM in our study site. The multiple functions of AM fungi associating with roots could also obscure the relationship between AM colonisation rates and P uptake.
Additionally, the lack of correlation (neither negative nor positive) between AM colonisation and root APase could suggest the possible complementary role of both strategies (Steidinger et al. 2015;Turner 2008). In this scenario, roots would benefit by investing in both strategies at the same time, since AM fungi and APase are known to explore different forms of P in soils: roots release phosphatase enzymes to mineralise organic P whilst AM fungi could benefit plants by exploiting inorganic P pools that do not need to be hydrolysed by enzymes or, more importantly, that have already been released by phosphatase enzymes (Nasto et al. 2014). This mechanism may be of special importance for fine absorptive roots in our study site, which displayed intermediate levels of AM colonisation and phosphatase exudation. This could also suggest that fine roots with high SRL and SRA may still not ensure adequate nutrient supply for some species in low fertility soils in the tropics, making them also dependent on AM (Siqueira and Saggin-Júnior 2001).
Investing in both enzyme production and AM association at the same time could also suggest that plants have access to enough resources, such as C and N, to support both strategies. In addition, mechanisms such as the exudation of organic acid/anions by roots could be an important plant reaction under P deficiency (Aoki et al. 2012;Lambers et al. 2006), but were not analysed in our study. Since our measurements were made at the ecosystem level, we were not able to analyse the individual relationships between root morphology, phosphatase activity and AM colonisation at the species level and therefore our interpretations could be different if we could control for host species. For instance, not all tree species have the same P requirements and P-use-efficiencies and therefore the investment in strategies to scavenge P could differ. Furthermore, plant benefits due to association with AM fungi would vary according to the tree and fungus identity, with similar levels of AM colonisation theoretically resulting in different P-uptake rates for different species (Herre et al. 2007;Maherali 2014;Siqueira and Saggin-Júnior 2001). Given the multiple and increasingly recalcitrant forms of P in highly weathered soils, enzyme production, association with AM fungi and a high diversity of root morphological traits, combined could maximize P uptake (Kong et al. 2014;Li et al. 2017;Zemunik et al. 2015), indicating that in soils with very low P availability, such those in Central Amazon, more than one P-uptake strategy is needed. Ultimately, understanding how root traits vary at the ecosystem level and how they relate to each other is critical to expand our knowledge about P cycling mechanisms, especially in the context of Plimitation in tropical forests with very high tree diversity.

Conclusions
In a P-poor tropical forest in Central Amazonia, we show that morphological traits of absorptive fine roots are related to root phosphatase activity. However, these roots are also associated with AM fungi, independent of their morphology. We suggest that: 1) the substantial investment in multiple root P-uptake strategies further emphasises the importance of P-limitation to the function of these forests; 2) mycorrhizal association could have other benefits for plants other than only P-uptake; 3) the relative abundance of other resources such as water, light, carbon and nitrogen may explain how plants are able to invest resources in multiple P-uptake mechanisms; and 4) the high diversity of tree species and the resulting diversity and combinations of P-uptake strategies makes detecting possible trade-offs between individuals or species challenging, which highlights the importance of ecosystem level estimates (i.e. inability of add up all possible strategies in a meaningful way). Overall, we conclude that because of the many forms that P is found in Amazon soils and the multiple steps needed to make P available for plants, roots also display a range of adaptations to enhance P-uptake under limitation. The hydrolysis of organic P by root phosphatase exudation and inorganic P-uptake by AM fungi could be seen as complementary mechanisms contributing to the functioning of this Central Amazon forest, at even very low P availability. mycorrhizal colonisation assessment and Dr. M. Liptak on behalf of EnviroScience Co. (https://www.enviroscienceinc.com/) for the donation that funded the purchase of the scanner. We are also thankful to R. Filho and the Laboratorio Tematico de Plantas e Solos (LTSP) at INPA for the soil P fractionation and the Laboratorio de Manejo Florestal (LMF -INPA) for helping with the logistics during the field campaign. This research was supported by the Brazilian National Council for Scientific and Technological Development (CNPq) and UK's Natural Environment Research Council (NERC), grant number NE/L007223/1; PM also acknowledges support from ARC, DP170104091.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.