Size‐dependent and environment‐mediated shifts in leaf traits of a deciduous tree species in a subtropical forest

Abstract Aims Understanding the joint effects of plant development and environment on shifts of intraspecific leaf traits will advance the understandings of the causes of intraspecific trait variation. We address this question by focusing on a widespread species Clausena dunniana in a subtropical broad‐leaved forest. Methods We sampled 262 individuals of C. dunniana at two major topographic habitat types, the slope and hilltop, within the karst forests in Maolan Nature Reserve in southwestern China. We measured individual plant level leaf traits (i.e., specific leaf area (SLA), leaf area, leaf dry‐matter content (LDMC), and leaf thickness) that are associated with plant resource‐use strategies. We adopted a linear mixed‐effects model in which the plant size (i.e., the first principal component of plant basal diameter and plant height) and environmental factors (i.e., topographic habitat, canopy height, and rock‐bareness) were used as independent variables, to estimate their influences on the shifts of leaf traits. Key Results We found that (1) plant size and the environmental factors independently drove the intraspecific leaf trait shifts of C. dunniana, of which plant size explained less variances than environmental factors. (2) With increasing plant size, C. dunniana individuals had increasingly smaller SLA but larger sized leaves. (3) The most influential environmental factor was topographic habitat; it drove the shifts of all the four traits examined. Clausena dunniana individuals on hilltops had leaf traits representing more conservative resource‐use strategies (e.g., smaller SLA, higher LDMC) than individuals on slopes. On top of that, local‐scale environmental factors further modified leaf trait shifts. Conclusions Plant size and environment independently shaped the variations in intraspecific leaf traits of C. dunniana in the subtropical karst forest of Maolan. Compared with plant size, the environment played a more critical role in shaping intraspecific leaf trait variations, and potentially also the underlying individual‐level plant resource‐use strategies.


| INTRODUC TI ON
Intraspecific trait variation is an increasingly important topic in traitbased ecology (Bolnick et al., 2011;Moran et al., 2016;Westerband et al., 2021). Understanding how plant development (hereafter simply development) and environment jointly shape plant leaf traits, will advance the understandings of the causes of intraspecific trait variation (Shipley et al., 2016;Violle et al., 2012). In turn, it would shed light on the understandings of species distribution (Grime, 1965), the assembly of ecological communities (Jung et al., 2010;Siefert et al., 2015), and the modeling of ecological processes (Funk et al., 2017;Moran et al., 2016).
For long-lived woody angiosperms, plant size is commonly used as a proxy of developmental stage (e.g., He & Yan, 2018;Martin & Thomas, 2013). In the forest, as a plant increases in size, it generally faces increasingly higher irradiance and hydraulic resistance at higher forest strata (Rozendaal et al., 2006;Thomas & Bazzaz, 1999), as a result, leaf traits usually show corresponding shifts toward more conservative resource-use strategy (e.g., Cavaleri et al., 2010;Dayrell et al., 2018;He & Yan, 2018;Ishida et al., 2005). Larger conspecific individuals usually have smaller sized leaves (Koch et al., 2004), smaller specific leaf area (SLA) (Kenzo et al., 2015;Thomas & Winner, 2002), as well as higher leaf dry-matter content (LDMC) (Park et al., 2019). These size-dependent trait shifts are interpreted as plant adopting an increasingly more conservative light and wateruse strategies in response to high light and drought stress, through acclimation and plasticity. However, there are exceptional cases (e.g., He & Yan, 2018;Thomas & Ickes, 1995). For example, Thomas and Ickes (1995) found larger conspecific individuals had larger sized leaves for understory treelet species, suggesting a more liberal resource-use strategy.
Evidently, development and environment are two major drivers of intraspecific trait variations. However, it remains unclear how development and environment jointly shape the intraspecific leaf trait shifts (Tredennick et al., 2018). On the one hand, researchers suggest that development and environment interact in mediating intraspecific trait variations, attributing to ontogenetic plasticity and acclimation to environmental heterogeneity (Russo & Kitajima, 2016). On the other hand, empirical evidence was equivocal, some found development and environment interactively (e.g., Dayrell et al., 2018), others found independently (e.g., Fortunel et al., 2020;Liu et al., 2020), in shaping the intraspecific trait variations. Furthermore, previous studies that estimated the relative importance of development and environment in shaping variations in intraspecific traits also found inconsistent results (e.g., Fajardo & Piper, 2011;Fortunel et al., 2020;Liu et al., 2020). For example, Fajardo and Piper (2011) found environment dominates over plant size in shaping leaf mass per area (LMA) of a deciduous angiosperm tree in southern Andes, but Liu et al. (2020) found that the opposite was true for a conifer tree in the forest of northern China.
An extensive zone of karst forests lies in southwestern China (Geekiyanage et al., 2018), but knowledge of karst forests leaf traits is limited (Geekiyanage et al., 2019). Under a subtropical monsoon climate, rugged topography of the karst forest results in fast soil erosion and leaching, accompany with slow soil formation on the carbonate bedrock (Li et al., 2008), together create sharp gradients of light, water, and soil nutrients (Bonacci et al., 2009;Geekiyanage et al., 2018;Guo et al., 2017). A large piece of the karst forest is preserved in Maolan National Nature Reserve (hereafter Maolan; Zhou, 1987). The karst forest of Maolan grows mainly on a geomorphological feature termed as karst peakcluster, which represents a cluster of several steep-sided karst hills sharing a common base (Sweeting, 1995). The peak-cluster karst forest of Maolan can be classified into a series of topographic habitats, starting from the high light, low moisture, and infertile hilltop toward low light, moist, and fertile slope and foothill (Gan & Mu, 1987;Zhang & Zhang, 1987;Zhou, 1987). Apart from this regular pattern, smaller local-scale environmental conditions in the karst forests are largely irregular, such as canopy condition (e.g., canopy height; Long et al., 2005) and microtopography (e.g., rock outcrops; Zhang et al., 2010). The karst forest experiences

K E Y W O R D S
Clausena dunniana, intraspecific trait variation, karst forest, leaf area, Maolan National Nature Reserve, specific leaf area

J E L C L A S S I F I C A T I O N
Functional ecology high frequency of canopy disturbance (i.e., treefall gap; Long et al., 2005), resulting in varied canopy height and light penetration to forest lower strata. Rocky outcrops (or rock-bareness) are a common feature and constrain the availability of soil water and nutrients (Huang et al., 2009;Zhang et al., 2007). Previous studies have confirmed the effects of these environmental factors on plant distribution and specialization (e.g., Geekiyanage et al., 2018;Guo et al., 2017;Li et al., 2019;Long et al., 2005;Zhang et al., 2013).
However, it remains unclear how these environmental factors and plant development jointly influence the variations in intraspecific traits in the karst forests.
Here, we explore how development and environment jointly (including their interactions and relative importance) drive the shifts of intraspecific leaf traits of Clausena dunniana, which is a common species widely distributed in the understory of the subtropical broad-leaved forest on the karst hills of Maolan. Specifically, we test the joint effects of plant size, topographic habitat, and local environmental conditions on the shifts of leaf traits (i.e., SLA, leaf area (LA), LDMC, leaf thickness (L th )) of C. dunniana. We predict that leaf traits of C. dunniana would shift toward more conservative resourceuse strategies (i.e., smaller leaf area and SLA, higher LDMC and leaf thickness) with increasing plant size, and this trend be more pronounced in high light, low moisture, and infertile environments (i.e., hilltop, low canopy, and high rock-bareness locations), compared to dark, moist, and fertile environments (e.g., slope, high canopy, and low rock-bareness locations).

107°52′10″-108°05′40″E) is located in Libo County of Guizhou
Province in southwestern China. Local annual mean temperature is 15.3°C, with January mean of 5.2°C and July mean of 23.5°C. Mean annual precipitation is 1,752.5 mm. The bedrock of the reserve consists of Carboniferous and Permian limestones and dolomites, with elevation varying from 430 to 1,078 m a.s.l. (Zhou, 1987). In the karst peak-cluster of Maolan, the height from the hilltop (i.e., the top of individual hill that forms the peak-cluster) to the bottom ranges from 150 to 300 m (Li & Li, 1987). More than 90% of the reserve is forested (Zhou, 1987). The study species dunniana is a dominant tree species distributed on the whole topographic gradients of the peak-cluster, with more even distribution across topographic habitats from hilltop to bottom, compared with most other tree species (Qin et al., 2018).

| Plant sampling and trait measurement
In late July of 2020, using a strip transect (60 m × 10 m) method (Buckland et al., 2007), we sampled 262 individuals of C. dunniana at two major topographic habitat types, the slope (152 individuals) and the hilltop (110 individuals), within the forests of six peak-clusters (as sites, sample size ranged from 33 to 66 individuals) evenly distributed in two regions (sample size were 144 and 118 individuals respectively) about 10 km apart within the reserve ( Figure 1). The slopes are located between the hilltop and the foothill. The hilltops include the top of the hills as well as ridges that connect nearby hilltops. These hilltops had thinner soil, higher soil calcium, and lower soil phosphorous concentrations compared with the slopes (Jin et al. unpublished data). For each sampled individual plant, we estimated its neighborhood canopy height and rock-bareness rate, and used it to represent local-scale environmental conditions. We estimated canopy height using a measuring pole mounted with a rod level.
For each sampled individual plant, we measured the basal diameter (0.1 m aboveground) of its stem/trunk using an electronic digital vernier caliper (PD-151, Pro'sKit, Shanghai, China) for individuals with basal diameter <5 cm, or using a diameter tape for larger individuals (i.e., basal diameter ≥5 cm). We measured the height of the sampled individual using a pole with a rod level. With respect to leaf sampling, we randomly sampled 2-10 fully expanded intact compound leaves located at the outer layer of the upper portion of the crown of the plant using a tree pruner (LZ5625, VMP, Shandong, China). The sampled leaves were placed in black plastic bags, stored in portable refrigerating box at above 0°C temperature and transported back to the laboratory in the same day. Then, the leaves were placed in fresh water for rehydration overnight. In the next day, the sampled compound leaves were scanned using a flatbed scanner, and leaf area was measured as the sum of one-sided projected fresh lamina surface area of all the leaflets without petiole and rachis, using ImageJ software (Schneider et al., 2012). For each sampled compound leaf, three leaflets were randomly selected. For each leaflet, thickness was measured as the mean of three measurements taken between the tip and base of the lamina, avoiding the major veins, using an electronic digital vernier caliper (PD-151, Pro'sKit, Shanghai, China). Then the average thickness of the three leaflets was calculated to represent the leaf thickness of the compound leaf. Leaf dry-matter content was calculated as the sum of dry leaf mass divided by the sum of fresh weight of the sampled compound leaves for each individual plant.

| Statistical analyses
The plant size-dependent and environmental factors (i.e., topographic habitat, canopy height, and rock-bareness rate) mediated shifts in leaf traits (i.e., SLA, leaf area, LDMC, and leaf thickness) were estimated by linear mixed-effects model using the "lmer" func- A random intercept (μ jk ) was specified for site j nested within region k, as random effect. The model structure was: where β 1 through β 4 are fixed-effect coefficients representing variation in SLA, LA, LDMC, or L th due to variation in plant size, topographic habitat, canopy height, and rock-bareness rate, respectively; γ represents fixed-effect coefficient for second-order interaction of these factors; and α represents the mean trait value of the individual plant averaged over all sites and regions. The random term μ jk was assumed to be normally distributed as jk ∼ N 0, . Before analysis, all continuous predictor variables were standardized by subtracting the mean and dividing by the standard deviation (i.e., mean = 0, SD = 1; Table 1). We tested variance inflation factors (VIF) for the predictor variables as well as their two-way interactions for the four models using the "vif.mer" function in the ieco package (Helmus, 2021) in R; we found all the VIF values were less than 3 for all the four models, which were substantially smaller than 10 (Dormann et al., 2013), suggesting that multicollinearity among the explanatory variables is not an issue in our study.  Furthermore, a variation partitioning analysis was conducted to estimate the relative importance of plant size and the three environmental factors in explaining the variances in the leaf traits, using the "varpart" function in the vegan package. All the analyses were conducted in R 3.6.1 (R Core Team, 2019).

| RE SULTS
The most supported models of the shifts of the four leaf traits explained significant amounts of variance in the data. The fixed part of the models explained a proportion of the variance (R 2 m ) ranged from 0.054 to 0.211, the total proportion of the variance explained by the models (R 2 c ) including the random part ranged from 0.146 to 0.491 (Table 2). However, the most supported models retained no two-way interactions, indicating that there was lack of significant interaction effects between plant size and environment on the leaf trait shifts of C. dunniana (Table 2). In the most supported model, the shifts of SLA and leaf area were driven by both plant size and environmental factors, the shifts of LDMC and leaf thickness were driven only by environmental factors (Table 2). Furthermore, variation partitioning analysis showed that plant size explained a smaller proportion of the variance in three of the four leaf traits than did environmental factors ( Figure 2).
With respect to plant size effects on leaf trait shifts, we found that larger sized C. dunniana individuals have smaller SLA and larger sized leaves (Figure 3a,b and Table 2). The plant size was not retained in the most supported model of LDMC and leaf thickness, suggesting a lack of plant size-related shifts of these two leaf traits (Table 2).
On the other hand, among the three environmental factors, the most important factor was topographic habitat, it propelled the shifts of all the four leaf traits we examined (Table 2). Specifically, comparing between topographic habitats, C. dunniana individuals in the slope habitat had larger SLA, lower LDMC, smaller sized, and thinner leaves than C. dunniana individuals on hilltops (Figure 4 and Table 2). On top of topographic habitat influence, leaf traits were further modified by local-scale environmental conditions. Canopy height was retained in the most supported model of leaf thickness, and rock-bareness rate was retained in the most supported model of SLA (Table 2). Specifically, at locations with increasingly higher canopies, C. dunniana individuals had increasingly thinner leaves ( Figure 3c and Table 2); at locations with increasingly higher rockbareness rates, C. dunniana individuals had increasingly smaller SLA ( Figure 3d and Table 2).
With respect to the independent effects of plant size and environment on the four leaf traits of C. dunniana examined, we found that size-dependent shifts were less frequent and overall, less influential than environment-mediated shifts. Our finding is consistent with the finding of environment dominated over plant size in controlling leaf trait of a widespread tree species Nothofagus pumilio in the southern Andes of Chile (Fajardo & Piper, 2011), but our finding is contrary to a previous study on a conifer tree species Pinus koraiensis in a temperate forest in northern China (Liu et al., 2020).
Since the relative importance of environment in driving intraspecific trait variations likely increase with the spatial extent and environment heterogeneity of a study system (Spasojevic et al., 2016), we suspect that the inconsistent results of these three studies at least partly stemmed from the difference in spatial extent and environmental heterogeneity encompassed. Specifically, Liu et al.'s study was conducted in a 9-ha forest plot lying on a gentle topography (Xu & Jin, 2012), which might have led to the relatively small role of environment in driving intraspecific leaf trait shifts. On the other hand, Fajardo and Piper's study covered large elevation ranges located at two distant regions, and the present study was conducted on rugged karst peak-clusters. These two studies covered much larger spatial extents and probably also larger environmental heterogeneities, hence facilitated more prominent roles of environment in driving intraspecific leaf trait shifts.
Generally, forest plants face increasingly higher light and hydraulic stresses as plant sizes increase and are closer to the canopy (Rozendaal et al., 2006;Thomas & Bazzaz, 1999). During the process of development, these stresses are expected to push the leaf traits toward greater conservative resource-use strategies (Dayrell et al., 2018). Our finding of the pattern of SLA declining F I G U R E 2 Venn diagram of proportion of variance of the leaf traits explained by plant size and the environmental factors SLA, specific leaf area; LDMC, leaf dry-matter content. Portion (a) is explained only by plant size, portion (b) is explained by both plant size and the environment, (c) is explained only by the environment F I G U R E 3 Shifts in leaf traits in response to plant size, canopy height, and rock-bareness rate as estimated by the most supported linear mixed-effects models SLA, specific leaf area; LDMC, leaf dry-matter content. Each empty circle represents an individual plant. The fitted line represents predicted values by the most supported linear mixed-effects model as shown in Table 2. The shaded region indicates 95% confidence interval of the predicted values with increasing plant size of C. dunniana, agrees with this expectation and is in line with the pattern widely reported from forests in other parts of the world (e.g., Kenzo et al., 2015;Martin & Thomas, 2013;Park et al., 2019;Thomas & Winner, 2002). However, the leaf area of C. dunniana enlarged with increasing plant size, which is contrary to expectation (Koch et al., 2004). As a matter of fact, the size-dependent increase in leaf area is not uncommon in forest plants (e.g., He & Yan, 2018;Ishida et al., 2005;Thomas & Ickes, 1995). For example, Thomas and Ickes (1995) found a differential size-dependent shift in leaf area between canopy and understory plant species in a Malaysian rain forest. Specifically, the leaf area of understory treelet species tended to increase with increasing plant size, whereas canopy tree species tended to show the reverse pattern in their study. Thomas and Ickes (1995) suggested that the contrary patterns might be partly due to the difference between canopy and understory species in adult stature and crown exposure. Canopy species are expected to have high crown exposure as adults, and they invest heavily on vertical growth and show leaf traits that are increasingly advantageous in higher irradiance environments (e.g., smaller sized leaves) toward the canopy. On the other hand, the small asymptotic heights of understory species mean low opportunity of high crown exposure as adults, and they invest increasingly more to horizontal growth, such as adopting leaf traits that capture more irradiance in understory environment (e.g., larger sized leaves).
Another possible driver of size-dependent shifts in leaf traits is reproductive onset (Thomas, 2010). Specifically, before reproductive onset, a plant invests resources mainly on vegetative growth, and adopts an increasingly acquisitive resource-use strategy as the plant meets with increasing light availability and achieving higher water-use efficiency during its development. After reproductive onset, the plant allocates more resource to reproduction, and switches to a more conservative leaf resource-use strategy. Under this scenario, one is expected to find a unimodal pattern of leaf trait shifts during development, peaked approximately at the reproductive onset stage (Thomas, 2010). In the present study, our sampling covered a wide plant size range (with a basal diameter range of 0.36-15.7 cm) that probably extends beyond the reproductive onset stage of C. dunniana, yet the monotonic shifts of leaf traits were apparent, therefore, we suspect that reproduction might not be an influential factor of plant size-related leaf trait shifts for C. dunniana.
Previous studies conducted in the karst forests found topography was closely related to plant distribution (Guo et al., 2017;Zhang et al., 2013) and trait specialization (Geekiyanage et al., 2018). In this study, we found topography was critical for leaf trait shifts of C. dunniana, which extended our understanding of the topographic effects on leaf traits in the karst forests from the interspecific level (Geekiyanage et al., 2018) to the intraspecific level. Specifically, Geekiyanage et al. (2018) found that the plant species specialized in karst hilltops had leaf traits representing more conservative resource-use strategy compared with species specialized in slopes and foothills. We found this topography-related leaf trait pattern also holds for the conspecifics, shown as the individuals of C. dunniana on hilltops had leaf trait values representing more conservative resource-use strategy, such as lower SLA and higher LDMC, compared to individuals on slopes.
F I G U R E 4 Differences in leaf traits between topographic habitats as estimated by the most supported linear mixed-effects models SLA, specific leaf area; LDMC, leaf dry-matter content. Boxplot shows the median, the first and third quartiles of the observed leaf traits, with whiskers extending to 1.5 times the interquartile range, and empty circles represent outliers. Significant difference (p < .05 as estimated by the most supported linear mixed-effects models in Table 2) between groups are indicated by different lowercase letters above the boxplots Local-scale rock-bareness rate mediated effects on leaf trait shifts likely stemmed from underlying water and edaphic conditions (Zhang et al., 2007). Since rock-bareness constrains the availability of water and soil nutrients in the karst forest (Huang et al., 2009), the decline of SLA with increasing rock-bareness in the neighborhood was possibly a plant response to greater water and edaphic stresses, and shifted the ecological strategies toward more conservative resource-use strategies. Furthermore, as the optimization theory prediction of leaf thickness be positively related to light (Pérez-Harguindeguy et al., 2016), the decrease in leaf thickness with increasing local-scale canopy height, as was found in this study, is possibly due to a response to the reduction in exposure and light with increasing canopy height (Clark et al., 1996).
In conclusion, our study shows that plant size and environment

CO N FLI C T O F I NTE R E S T
All authors declare no conflicts of interest.

DATA AVA I L A B I L I T Y S TAT E M E N T
The data used in this study are archived in Dryad Digital Repository: https://doi.org/10.5061/dryad.rfj6q 57c4.