Testing for functional significance of traits: Effect of the light environment in tropical tree saplings

Abstract Functional traits have been examined to explain the growth rates of forest communities in different sites. However, weak or nonexistent relations are often found, especially due to the following methodological aspects: 1) lack of an environmental context (e.g., light, water, or nutrient supply), 2) use of nonfunctional traits, 3) an approach that does not contemplate phenotypic integration, and 4) neglect of intraspecific variation. Here we measured relative growth rates, crown, and leaf traits in saplings of six tropical tree species growing in two light environments (Gap and Understory) to test whether contrasting light environments modulates trait–trait and trait–growth relationships. Moreover, we tested whether models that integrate traits of different dimensions of the plant (crown and leaf) improve the strength of trait–growth relations. Light availability changed both trait–trait and trait–growth relationships. Overall, in Understory, crown traits (crown length and total leaf area) have a stronger effect on growth rates, while physiological traits related to nutrient acquisition (nitrogen concentration), photochemical efficiency (chlorophyll pigments and chlorophyll a fluorescence), and biochemical efficiency (potassium use efficiency) are strong in Gap. Models including multiple traits explained growth rates better in Gap (up to 62%) and Understory (up to 47%), but just in Gap the best model comprises traits that are representative of different dimensions of the plant. Synthesis. We advanced the knowledge behind the light effects on tree sapling by posit that trait–trait and trait–growth relationships vary across light environments. Therefore, light availability is a key environmental factor to be considered when choosing the set of traits to be measured in functional approach studies using tropical tree saplings. In compliance with the phenotype integration hypothesis, functional traits are better predictors of growth rates when grouped in a set of traits of different dimensions of the plant that represent different functional mechanisms.


| INTRODUC TI ON
Functional traits are by definition any morphological, physiological, or phenological measurements performed at an individual level that indirectly affect fitness due to effects on some aspect of plant performance (i.e., growth, survival, or reproduction) (Violle et al., 2007).
Ecologists have long used traits to explain demography performance affected by abiotic and biotic environmental conditions around the world, especially growth rate (Chave et al., 2009;Paine et al., 2015;Wright et al., 2004). However, studies that found a moderate or strong trait-growth rate relationship are an exception . In this sense, recent studies have concentrated efforts in understanding why weak relationships are frequent and defining standardized protocols for studies of trait-based community ecology (Caruso et al., 2020;Swenson et al., 2020;Worthy & Swenson, 2019;Yang et al., 2018). Some methodological aspects of sampling may affect the ability of traits to predict growth rates, such as 1) lack of environmental context, 2) use of nonfunctional traits, 3) approaches that do not contemplate phenotypic integration, and 4) neglect of intraspecific variation (Yang et al., 2018;Yang et al., 2020;Worthy & Swenson, 2019).
To understand how traits determine plant growth, it is necessary to understand how traits correlate to form plant phenotype and how environmental context modulates these relationships (He et al., 2020;McGill et al., 2006). Globally, trait-trait relationships were ordered in a single axis of variation, where species with high specific leaf area, high photosynthetic rate, high nitrogen content, and low wood density (among other traits) have acquisitive strategies and a high growth, while the opposite trend reflects conservative strategies and low growth (Chave et al., 2009;Reich, 2014;Wright et al., 2004). However, these trends of covariation may change among environments (Nicotra et al., 1997;Poorter et al., 2014;Wright et al., 2005), from global to local scale (Messier et al., 2017) and from interspecific to intraspecific analysis (Laughlin et al., 2017).
For example, a maximum photosynthetic rate (A max ) increases positively according to specific leaf area (SLA) among species (Wright et al., 2004), but across a light gradient the intraspecific correlation is negative, as acclimation of leaves to high light increases A max but reduces SLA (Kenzo et al., 2015;Santos & Ferreira, 2020b).
Furthermore, trait-growth relationships may also change according to light environment (Poorter, 1999). For example, in high-light environments, physiological traits related to photosynthetic nutrient use efficiency are highly related to growth (Guimarães et al., 2018), while in low-light environments, morphological traits related to light interception are highly related to growth (Liu et al., 2016). Studies that examine the effects of environmental contexts on trait-trait and trait-growth relationships are necessary, especially in tropical forests (Worthy & Swenson, 2019).
Beyond the environmental context, a good selection of traits is fundamental. The traits commonly measured and the type of analyses conducted represent little the complexity of the growth process.
Most studies developed to date have focused on measuring soft traits due to an easy applicability and a lower required investment (Gibert et al., 2016;Worthy & Swenson, 2019;Yang et al., 2018).
However, hard traits are more informative of physiological processes that drive plant performance because they highlight, for example, the different steps of photosynthesis (Santos & Ferreira, 2020b;Volaire et al., 2020;Worthy & Swenson, 2019;Yang et al., 2018).
Additionally, a set of traits can better explain demographic rates than a single trait does (Li et al., 2017) because plant growth depends on a complex integration of multiple traits. Therefore, models that allow integration of many traits measured in different dimensions of the plant, such as crown, leaf, roots, and stem (the phenotypic integration hypothesis), that represent mechanisms related to demographically complex processes, such as growth, may have a greater explanation power (Yang et al., 2020).
The great variation in outcomes found in the literature on functional trait-demographic rate relations is one of the most important topics in plant biology (Salguero-Gómez et al., 2018). Considering light, the most limiting factor to plant growth in tropical forests (Goldstein et al., 2016;Graham et al., 2003;Nicotra et al., 1999;Wagner et al., 2017), the role the availability of this resource plays in trait-growth relationship remains unclear. Thus, to fill this gap, an analysis that controls light variation and species composition can better demonstrate how traits determine the performance of plants (Wright et al., 2010). Therefore, based on a two-year experiment of secondary forest enrichment planting with six tropical tree species introduced in artificial gaps, we selected individuals growing in contrasting light environments (Gap and Understory) and then measured thirty traits (from crown and leaf) at an individual level.
We analyzed the relationships among traits and between traits and growth rate and tested the following three questions. First, how does light environment modulate trait-trait relationships in tropical tree saplings? We hypothesize that the plant phenotype is manifested by different trait-to-trait relationships at the different light environment. Second, does the light environment modulate the relationships between traits and growth rate? We hypothesize that light environment changes the plant phenotype in a way (first hypothesis) that affects traits most related to growth. Third, do trait-growth relationships improve by using an integrated set of traits, measured in different dimensions of the plant (crown and leaf)? We hypothesize that a set of traits of different dimensions of plants-the whole phenotype hypothesis-and physiological traits (e.g., photosynthetic process) represents growth better in both environments.

| Study site
This study was conducted in a secondary forest enrichment planting.
However, the area was abandoned, giving way to a naturally regenerating forest. Therefore, silvicultural treatments were applied during the last quarter of 2016 to begin the enrichment planting experiment in a 19-year secondary forest. The site is at the Fazenda Experimental da Universidade Federal do Amazonas (02º38'S, 60º03.5'W), 38 km north of Manaus, Amazonas, Brazil. The annual rainfall in Manaus is 2,350 mm, the monthly average air temperature ranges from 26.4 to 28.5ºC, and the average air humidity ranges from 75% in the dry season to 85% in the wet season (data from 1988-2018INMET, 2019).
Rainfall seasonality is generally moderate, with a dry season (rainfall <100 mm) between August and September (Sombroek, 2001).

| Experimental design
The light manipulation from the secondary forest to the subsequently enrichment planting consisted of understory slashing: removal of all herbaceous plants and small trees (DBH < 5 cm) at the beginning of the experiment and twice a year, and canopy refinement: bringing trees down (DBH ≥ 5 cm) by the canopy and subcanopy with a chainsaw in progressive levels of basal area reduction (0, 20, 40, 60, 80, and 100%). Each level of canopy refinement was applied to a plot of 2,318 m 2 (61 × 38 m), and understory slashing was applied in a subplot half the size of the main canopy refinement plot ( Figure S1). Treatment combinations were replicated five times (blocks), each with 12,768 m 2 (114 × 112 m). Six tropical tree species were selected by the criteria of economic and social importance, life history strategy, and availability of seeds following a minimum number of five matrices. The species are of two ecological groups: long-lived pioneer (Cedrela fissilis, Tabebuia rosea, and Swietenia macrophylla), and partial shade-tolerant (Bertholletia excelsa, Carapa guianensis, and Hymenaea courbaril) (Chazdon, 2014;Finegan, 1992;Swaine & Whitmore, 1988). A preview study within the same experiment showed a wide range of growth responses of these species to a light gradient created by thinning (Santos & Ferreira, 2020a). The seedlings were planted in these plots in March 2017. In each subplot, five seedlings of each species were planted, 3 × 3 m apart, with an edge of 10 m between plots and 11 m between subplots ( Figure S1). During planting, the soil surrounding each seedling was fertilized with P 2 O 5 (46 g), N (11.6 g), KCl (12 g), and micronutrients (10 g of FTE-BR12: 1.8% B; 0.8% Cu; 3.0% Fe; 2.0% Mn; 0.1% Mo). The application methods and doses that ensured the minimal nutrient requirements of tropical trees were obtained from a literature review (Alvarado, 2015;Campoe et al., 2014;Furtini Neto et al., 2000;Resende et al., 2005). Further details are provided in Santos et al. (2020) and Santos and Ferreira (2020a).
For this study, we selected only two subplot combinations with the most contrasting light environments: 1) clear-cut (100% of canopy refinement) and understory slashed subplots-hereinafter named Gap (photosynthetic photon flux density-PPFD = 27.9 mol/m 2 day −1 ), and 2) plots without basal area reduction (0% of canopy refinement) combined with understory slashing-hereinafter named Understory (PPFD = 1.3 mol/m 2 day −1 ) (Santos et al., 2020) (Figure S1). We selected only subplots with understory slashed with the objective to reduce the eventual effects of herbaceous competition on sapling growth.

| Growth rate measurements
We monitored saplings planted in this experiment for two years in bi-monthly campaigns from March 2017 to March 2019. In each monitoring campaign, we measured the root collar diameter (D) (5 cm above the soil) and total height (H) of each sapling. The aboveground biomass was approximated by D 2 H, as Kohyama and Hotta (1990) proposed. We calculated annual relative growth in biomass for each individual according to Hunt (1978).
where RG x = relative growth in biomass, lnX 2 − lnX 1 = increase in logarithmic biomass between two measurements, and T 2 -T 1 = interval between measurements.

| Trait measurements
We measured a set of crown and leaf traits at the end of the second year after planting ( Table 2). We measure all live individuals in each plot of each treatment selected and used individuals values for all analyses. In total, we measured the traits in 113 individuals in Gap, ranging from 16 to 23 per species, and 89 individuals in Understory subplots, ranging from seven to 22 per species. The differences in sampling between light conditions and across species are due the greater mortality of individuals in Understory, especially of shadeintolerant species (long-lived pioneers). All leaf traits were measured on healthy, fully expanded leaves visually without herbivory in the middle third of the crown. We grouped traits in two main groups (crown and leaf) and leaf traits in six sets of traits: leaf display, nutrient acquisition, light absorption efficiency, light use efficiency, gas exchange, and nutrient use efficiency.

| Crown
We measured crown diameter in two orthogonal directions and crown length (CL) from the first leaf insertion (or branch) to the top of the crown. Then, we calculated mean crown diameter (MCD, average of two orthogonal crown diameters), crown projection area (CPA, using the formula for ellipse area), crown ratio (CR, MCD/D), crown length ratio (CLR, CL/H), and relative crown length (RCL, CL/ MCD) (Li et al., 2017).
We estimated total leaf area (TLA) as follows: First, we counted the numbers of leaves in each plant and measured leaf length and width in the field for 30% of leaves. Then, we used leaf-size and species-specific models that relate leaf area to leaf length and width to estimate the area of leaves measured in the field (length and width). The parameters of the models were adjusted using leaves previously measured in laboratory by a leaf area meter (CI-202, CID, from 0.86 to 0.99 (Table S1). Finally, we estimated the total leaf area of each plant using the proportion of leaf area in 30% of leaves. From TLA, we also calculated leaf area index (LAI) as the ratio between TLA and CPA (Poorter, 1999).

| Leaf display
We measured the leaf area of a range of three to six leaves for each individual using a leaf area meter (CI-202, CID, Inc. Beds) and then calculated the average to obtain the individual leaf area of each plant (LA). Ten leaf disks (2.83 cm 2 ) were cut from three leaves per plant, saturated in distillated water, and oven-dried at 65ºC. Specific leaf area (SLA) was calculated as the ratio between fresh leaf area (cm 2 ) and oven-dried mass (g). Leaf dry matter content (LDMC) was calculated as the ratio between oven-dried mass (mg) and water-saturated mass (g) (Pérez-Harguindeguy et al., 2013).

| Nutrient acquisition
We estimated mass-based nitrogen, phosphorus, and potassium concentration in three to six leaves per plant. Leaf nitrogen (N mass ) was determined by the Kjeldahl method using digestion, distillation, and titration (Bremner, 1996). Leaf phosphorus (P mass ) was determined by spectrophotometry (λ = 750 nm) following the molybdate method (Murphy & Riley, 1962). Leaf potassium (K mass ) was obtained by atomic absorption spectrometry (1100B, PerkinElmer).

| Light absorption efficiency
We estimated mass-based chloroplast pigment concentrations in three leaves of each plant. Chlorophyll a (Chl a), chlorophyll b (Chl b), and carotenoids (Car c+x ) concentrations were extracted with acetone (10 ml of 80% acetone with 0.05 g of MgCO 3 per 0.1 g of fresh leaf) followed by filtration (Lichtenthaler & Wellburn, 1983).
The concentrations were calculated by absorbance reading at three wavelengths (663, 645, and 480 nm; Biochrom Libra s50 UV/Vis) following Hendry and Price (1993). Total chlorophyll concentration (Chl a + b) and chlorophyll a and b ratio (Chl a/b) were calculated by sum and ratio of Chl a and Chl b concentrations, respectively.

| Light use efficiency
The fluorescence parameters were measured in twelve leaves per plant between 8 a.m. and 10 a.m. using a portable fluorimeter (PEA, MK2-9600, Hansatech) adjusted to emit a saturating light pulse of 3,000 µmol/m 2 s −1 at a wavelength of 650 nm for 1 s. Before the measurements, the leaves were acclimated to the dark during 30 min for a complete oxidation of the photosynthetic electron transport chain. Maximum quantum yield of PSII (F V /F M ), ABS-based performance index (PI ABS ), and total performance index (PI total ) were calculated by performing JIP test applied to chlorophyll a polyphasic transient stages (Strasser et al., 1995;Strasser et al., 1999;Strasser et al., 2010;Tsimilli-Michael & Strasser, 2008  PI total is an integrated index of the PI ABS and the efficiency at which electrons reduce the end acceptors in the PSI.

| Gas exchange
We measured leaf gas exchange on one leaf per plant using an infrared gas analyzer (IRGA's, LI-6400XT, LI-COR). Prior to recording, tests were performed on at least three leaves per plant to ensure the selection of leaves with the maximum values of A max and g s . The measurements were made between 8 a.m. and 11 a.m. The IRGA's chamber was adjusted to a flow rate of 400 µmol/s, CO 2 concentration of 400 µmol/mol, H 2 O vapor concentration of 21 mmol/mol, leaf temperature of 31ºC, and photosynthetic photon flux density (PPFD) of 2,000 µmol/m 2 s −1 for maximum photosynthetic rate (A max ), stomatal conductance (g s ), and transpiration rate (E), and 0 µmol/m 2 s −1 for respiration rate in the dark (R d ).

| Nutrient use efficiency
We calculated carbon use efficiency (CUE) as the ratio between A max and R d , and nitrogen (PNUE), phosphorus (PPUE), and potassium (PKUE) use efficiency as the ratios between maximum photosynthetic rate per unit mass (A max × SLA × 0.1) and the nutrient concentrations in molar mass units.

| Statistical analyses
We assessed trait and growth rate variation in the 5th and 95th percentile To evaluate whether trait-growth relationships were modulated by light environment, we used an individual-level analysis with traits measured on the same individuals that they were used to predict growth. We adjusted mixed effects models with growth rate as the dependent variable, traits as the independent variable, light environment as the independent categorical variable, and we tested species and plots as random effects. We performed this analysis for each trait separately. We observed interaction effects on relationships with a significance level (p < 0.05). When we observe significant interaction effects on relations, we compare the slopes between light environments and showed uncertainty around the estimates. The best models (random effects definition) were selected based on the Akaike's information criterion (AIC). In all models, we observed the coefficients of determination both marginal (for fixed effects) and conditional (for fixed and random effects). Logarithmic transformations were performed to meet linearity and normality of assumptions whenever necessary.  Figure S2). Then, we used the scores of general PCA axes that jointly explained most of the variation (~60%) as independent variables in the models.
Secondly, we performed PCA for traits of each group (crown, leaf display, nutrients acquisition, light absorption efficiency, light use efficiency, gas exchange, and nutrient use efficiency) separately (categorized PCAs). Next, we used the scores of the first and second axis of each categorized PCA as independent variables in the models. We performed categorized PCAs because they summarize plant mechanisms in a few axes of variation and can improve the explanatory power of models due aggregate the explanation power of many traits, the same idea is applied to general PCA.
Finally, we selected the best model for all traits based on the AIC.
We also observed the best model with a single trait to compare how much the multiple models improve the power of explanation of growth rates. All traits were standardized, subtracting each individual trait value from the trait mean and then dividing by its standard deviation. Coefficient of determination marginal and significance level (P value) were observed to determine the strength of the relations.
The PTNs were fitted using the igraph package (He et al., 2020), and mixed effects models were fitted using the nlme package (Zuur et al., 2009). All analyses were performed using the statistical software R, version 3.5.1 (R Core Team, 2018).

| Effects of light environment on trait-growth relationships
Light environments modulated trait-growth rate relationships for some traits, as the significant interaction effects show ( Figure 2).
Overall, traits related to light interception (CL and TLA) had higher positive effects on growth rates in Understory than in Gap, with slopes of the relations being 52% to 66% higher. However, for CLR and RCL, the slope of the relation in Gap was 122% and 233% greater than times higher. Table S2 shows the values of slopes and uncertainties.
TA B L E 1 Name, abbreviation (AB), units, and 5th and 95th percentiles of relative growth rate and thirty measured traits in Gap and Understory

| General and categorized principal component analyses
Most (crown) of trait variation ( Figure S4).

| Relationships among growth rates and an integrated set of traits
The combinations of traits and the power of explanation of models changed according to light environments and the set of traits (single, multiple, or multiple with PCA scores, Table 2). First, the best models with a single trait explained 50% of growth variation in Gap (CL) and 36% of growth variation in Understory (MCD). Secondly, the best multiple traits model including all traits in their original form explained 62% of growth variation at Gap (MCD + CL + CR + PI total ) and 47% at Understory (MCD + CL + CR). Thirdly, the models including the first four axis of general PCA explained 41% of growth variation in Gap (PC1 + PC2 + PC3 + PC4) and 5% in Understory (PC1). Finally, the models including the first and second axes of categorized PCAs explained 52% of growth variation in Gap (PC1 crown + PC2 crown + PC1 light use efficiency + PC2 light use efficiency ) and 28% in Understory (PC2 crown + PC2 nutrients acquisition ). Therefore, all F I G U R E 1 Plant trait network analyses (PTNs) of crown traits and leaf traits in two light environments. Modules that reflect biological mechanisms of plants are indicated by the colors of circles and shaded regions. Lines represent a significant relationship (connection) between traits. Black lines represent the connection between traits of the same module. Red lines represent the connection of traits of different modules. PTNs with a higher number of connections (lines) have a higher edge density. For abbreviations, see Table 1 in the section Results

F I G U R E 2
Effect of light environment on trait-growth relationships for traits with significant interaction effects. Points displaying values of growth and trait of individuals. Shading displaying 95% confidence intervals. For abbreviations, see Table 1 in the section Results models explained better growth rates in Gap than in Understory and multiple trait models improved growth prediction, however not using a PCA score approach.

| D ISCUSS I ON
The relationships between functional traits and demographic rates consist in the main framework of trait-based plant ecology. They are a major challenge in studies on tropical tree ecology (Chave et al., 2009;Wright et al., 2004Wright et al., , 2010. However, ecologists have reported many traits as poor predictors of tree growth, leading to questioning whether these traits are indeed functional or not (Poorter et al., 2018;Swenson et al., 2020;Worthy & Swenson, 2019;Yang et al., 2018). Three core reasons have been stated for the lack of success of traits in determining demographic outcomes in traitbased functional approach ecology: frequently ignored contextual information (e.g., environmental gradients), functional processes (e.g., photosynthesis) that determine plant performance are not completely explained using only soft traits or a single trait, and focus on species relative to individuals (Worthy & Swenson, 2019;Yang et al., 2018). Here we combined multiple morphophysiological traits and provided evidence on the importance of including the environment context in studies of trait-based community ecology due its role in modulating the trait-growth relationships of tropical tree saplings. Moreover, we reinforce that growth is the result of the expression of many traits, but may be decoupled of an integrated plant phenotypic expression in light-limited environments. Even though our results are for six species, these species show a wide variation in the whole plant economic spectrum (Reich, 2014), as values of traits show (Table 1).

| Relations among traits change according to light environment
Light environment affected relationships among traits and consequently integrated phenotype (Figure 1). Growing in Gap, all plant crown traits were coordinated and formed a single module, while in Understory only crown traits related to crown depth formed a dense module. Crown depth in Understory can reflect a conservative mechanism because we observed that plants with deep crowns had older leaves inserted at the lowest stem height and some since planting. Plants growing in poor environments show a long leaf lifespan and tend to be more successful in surviving because they compensate for the high payback times needed to return carbon investments on leaf construction (Kitajima & Poorter, 2010;. Evidently, photochemical and biochemical efficiency traits were grouped into two distinct modules in Gap. However, these modules are similar and have many connected traits, proving that photochemical and biochemical efficiency are strongly coordinated (Chou et al., 2020;Santos et al., 2019). Traits related to biochemical and photosynthetic nutrient use efficiency were "hub traits" in both light environments. However, LDMC and SLA had a high connectivity (>50%) only in the network of Understory. Therefore, photosynthetic efficiency is key in both Understory and Gap; however, in a light-limited understory, a conservative strategy related Notes: Random effects of the best models are controlled for intercepts of species. Coefficient of determination marginal (R 2 m ) of relationships. P values of traits in the best models are * p < 0.05, ** p < 0.01, *** p < 0.001, not significant (n.s.). p ≥ 0.05. For abbreviations, see Table 1 in the section Results. TA B L E 2 Results of mixed effect models examining the relations between growth rates and an integrated set of traits and between growth rates in both light environments to physically structured and durable leaves for protection against herbivores and long-term carbon investment also play a role (Kitajima, 1994;Kitajima & Poorter, 2010).

| Relationships between growth rates and functional traits are modulated by light environment
Overall, in Understory, light interception traits (LA, CL, and TLA) play a more important role in growth than in the Gap (Figure 2). The disproportional effects of light foraging traits are expected in more limited light environments and have been reported as a strategy to optimize the interception of diffuse light or sun flecks reaching the forest floor (Delagrange et al., 2006;Iida et al., 2014;Niklas, 1989). On the other hand, traits representing mechanisms of nutrients acquisition (N mass ), photochemical efficiency (Chl a, Chl b, Car c+x Chl a + b, PI ABS ), and biochemical efficiency (PKUE) more strongly affected the growth rates in Gap than in Understory. Plants growing in Gap have a greater amount of light energy for use in the physiological process (Chazdon & Fetcher, 1984). In this sense, it is more important for saplings to have a high potential for processing this energy in Gap than in Understory (Figure 2). Our results support the hypothesis that the relationships between functional traits and growth rates of tropical tree saplings depend on the environmental context (Worthy & Swenson, 2019) and specifically light environment.

| Growth rates are better explained by a multiple set of traits
Our results sustain the hypothesis that an integrated set of traits measured in different dimensions of the plant better explains plant growth. Moreover, the set of traits related to growth and the strength of relations changed between contrasting light environments (Table 2). Growth models that include multiple traits of different dimensions of the plant in the original form of traits or in composite traits better explain growth rates because they allow performing a phenotypic integration, that is, the integration of multiple traits that form the phenotype (Yang et al., 2020). The selection of traits in different plant organs (e.g., crown, leaf, stem, branches, and roots), which are related to different mechanisms that affect the growth process, is interesting for phenotypic integration, making it possible to have a broad view of the main factors that determine plant growth (Rosas et al., 2019;Swenson et al., 2020;Wright et al., 2005;Yang et al., 2020). The production ecology model (PEM) theory states that plant productivity depends on the availability of natural resources and on the efficiency of plant in acquiring and using those resources (Binkley et al., 2004;Monteith, 1977). Thus, measuring traits related to the different stages of PEM is important to understand the factors that determine the differential performance of plants between light environments. Light interception traits at a crown level (MCD, CL, and CR) were included in the best models in both Gap and Understory, but one light use trait (PI total ) was also present in the best model for Gap. Therefore, an integrated approach with multiple trait models improves the growth prediction in the sense that more mechanisms are represented in high and low-light environments.
Although composite traits (such as those summarized on a PCA axis) could be potential predictors of growth rates (Liu et al., 2016;Yu et al., 2020), in our study they did not capture the variation of the main processes related to growth, mainly in a light-limited environment ( Table 2). The main leaf traits (hubs) observed in plant network traits did not compose the best models in both environments. In Understory, the best multiple trait model comprised no leaf trait. This result may be interpreted as a decoupling between the integrated phenotype expression of leaf traits and growth performance. In this sense, the mechanisms behind leaf traits could be expressed as a strategy to maximize another component of plant fitness (e.g., sapling survival), mainly in a light-limited environment that demands conservative strategies (Kunstler et al., 2016;Wright et al., 2010). Moreover, the strength of relationships was higher for Gap than for Understory (Table 2). These results corroborate with the hypothesis Poorter et al. (2018) proposed. The authors stated that low-light availability harms the expression of growth and traits, resulting in a convergence of these values and hence weak or nonexistent relationships.

| Implications for functional approach studies
Gap dynamics is one of the main factors that leads to biodiversity in tropical forests. However, to investigate demographic rates in the gap dynamics is expensive and time-consuming (Chazdon, 2014;Connell, 1978). Thereby, the use of functional traits to uncover mechanisms that control the dynamics of forest communities following appropriate methodologies is an option for quick results in ecology studies (Salguero-Gómez et al., 2018;Swenson et al., 2020;Worthy & Swenson, 2019;Yang et al., 2018).
We add knowledge to this functional approach by showing that a complex process such as growth is better represented by multiple trait models, with traits measured in different dimensions of the plant. Thus, future studies measuring traits in different plant organs (i.e., crown, leaf, stem, branches, and roots) that are related to different processes of plant growth can provide a broader view of factors that determine plant growth. Following this approach, our results also highlight the importance in considering the environmental context (e.g., light supply) during the selection of traits to be measured. Traits measured at crown level are more associated with growth rates of saplings in low-light environments but are also good predictors of growth rates in high-light environments.
Nevertheless, PI total , an index of the entire photochemical efficiency accessed by chlorophyll a fluorescence, is a good predictor of growth rates in a high-light environment. Chlorophyll a fluorescence parameters are related to photochemical stage of photosynthesis and have been considered a proxy of carbon assimilation in tropical trees along the vertical profile of Central Amazonian forest (Santos et al., 2019). Therefore, chlorophyll a fluorescence parameters should be incorporated in functional studies because, in addition to its power to explain sapling growth, it also represents a fast and noninvasive method. Finally, due to the capacity of plants to achieve similar performances with alternative designs of traits (Dias et al., 2020;Marks & Lechowicz, 2006;Worthy et al., 2020), future studies should consider the intraspecific traitgrowth relation of each species and/or ecologic group.

| CON CLUS IONS
This study provides evidences on the importance of considering the environmental context (light supply) when deciding which traits to measure in future functional ecological researches. Functional traits seem to be better predictors of growth performance in high-light environments, especially when grouped in multiple traits. Multiple trait measurements in different dimensions of the plant that represent different processes of plant growth can contribute significantly to key ecological questions related to the main mechanisms behind tropical forest dynamics.

ACK N OWLED G M ENTS
We would like to thank the Fazenda Experimental da Universidade

CO N FLI C T O F I NTE R E S T
None declared.