- 1Co-Innovation Centre for Sustainable Forestry in Southern China, Bamboo Research Institute, College of Biology and Environment, Nanjing Forestry University, Nanjing, China
- 2School of Natural Sciences, Macquarie University, Sydney, NSW, Australia
- 3Institute of Botany, Jiangsu Province and Chinese Academy Sciences, Nanjing, China
- 4School of Integrative Plant Science, Cornell University, Ithaca, NY, United States
Foliage leaves are essential for plant survival and growth, and how plants allocate biomass to their leaves reveals their economic and ecological strategies. Prior studies have shown that leaf-age significantly influences leaf biomass allocation patterns. However, unravelling the effects of ontogeny on partitioning biomass remains a challenge because it is confounded by the effects of environmental factors. Here, we aim to elucidate whether leaf-age affects the allocation to the lamina and petiole by examining leaves of known age growing in the same general environmental context. We sampled 2698 Photinia serratifolia leaves developing in the same environment from April to November 2021, representing eight leaf-ages (n > 300 for each leaf-age). Petiole and lamina biomass, and lamina area were measured to evaluate the scaling relationships using reduced major axis regression protocols. The bootstrap percentile method was used to determine the differences in scaling exponents among the different leaf-ages. ANOVA with Tukey’s HSD was used to compare the ratios of petiole and lamina biomass to lamina area across the leaf-ages. Correlation tests were used to determine if exponents, intercepts, and ratios differed significantly across the different leaf-ages. The data indicated that (i) the ratio of petiole and lamina biomass to lamina area and the scaling exponent of lamina biomass versus lamina area correlate positively with leaf-age, and (ii) the scaling exponent of petiole biomass versus lamina area correlates negatively with leaf-age. Leaf maturation process involves an inverse proportional allocation between lamina and petiole biomass for expanding photosynthetic area. This phenomenon underscores the effect of leaf-age on biomass allocation and the importance of adopting an ontogenetic perspective when entertaining plant scaling theories and unravelling the principles governing shifts in biomass allocation throughout the leaf lifespan.
Introduction
Foliage leaves are the primary photosynthetic organs of terrestrial plants and therefore are critical to plant survival and growth (Rascher and Nedbal, 2006; Díaz et al., 2016; Adams III and Terashima, 2018). The photosynthates produced by leaves are utilized to sustain metabolism, support the functions of other organs, and regenerate damaged body parts (Nakamura and Hashimoto, 1988; Lockhart et al., 2003; Wright et al., 2004). However, the availability of raw materials for photosynthesis constrains the production of photosynthates, prompting plants to balance the allocation of resources to leaves, stems, and roots according to the physiological demands of each organ-type (Enquist and Niklas, 2002; Cheng et al., 2015). Prior physiological studies using different measurement methods provide insights into the utilization and transportation of photosynthates (Farrar and Jones, 2000; Minchin and Thorpe, 2003; Jahnke et al., 2009), and the final metabolic outcome of these processes can be assessed by measuring the biomass allocated to each organ- or tissue-type (Poorter et al., 2012).
Eudicot leaves are composed of two basic structural units, the lamina and the petiole, and the distribution of biomass between these two different units reflects a trade-off between their respective functions (Niklas, 1992; Niklas, 1999; Niklas and Enquist, 2002; Pasini and Mirjalili, 2006). This distribution, referred to as “leaf biomass allocation”, provides valuable information about plant economics and ecological strategies, and its variations reveal the ecological and evolutionary processes that shape plant development (Niklas, 1994; Westoby et al., 2002; Niklas et al., 2007).
Previous studies have applied two approaches (i.e., the ratio and scaling approach) to characterize leaf biomass allocation (Poorter et al., 2012; Poorter et al., 2015). The first approach involves determining the quotient of biomass at a specific time to describe the allocation of biomass to a particular structural or functional trait, such as lamina dry mass per unit area (LMA), which represents the leaf-level cost for light-harvesting (Poorter et al., 2009). Wright et al. (2005) quantified climatic influences on key leaf traits globally and reported that plants tend to have higher mean LMA at sites experiencing greater environmental stress (e.g., hotter, drier, and higher irradiance sites).
In contrast, the scaling approach employs a power-law equation to describe the relationship between or among different body parts (Niklas, 1994; Enquist et al., 2007; Niinemets et al., 2007; Niklas et al., 2009; Sack et al., 2012). For example, the lamina biomass (M) and area (A) scaling relationship exhibits a log-log linear relationship with a scaling exponent (i.e., the log-log slope denoted by α) whose numerical value reflects different leaf biomass allocation patterns affecting light interception and mechanical support. An isometric scaling relationship (i.e., α = 1.0) indicates that increases in leaf biomass obtain a one-to-one proportional increase in area (i.e., ). However, Niklas et al. (2007) found that the numerical value of the interspecific scaling exponent often exceeds unity (i.e., α > 1.0), indicating that increases in lamina biomass do not obtain proportional increases in lamina area. This phenomenon has been called “diminishing returns” (also see Milla and Reich, 2007).
However, a limited number of studies have tested leaf-age’s effects on biomass allocation patterns. Jiao et al. (2022) reported that the numerical value of the scaling exponent α increases as Photinia × fraseri “Red Robin” leaves grow older, resulting in a higher M vs. A scaling exponent for leaves sampled in the summer compared to leaves sampled in the spring. However, the proximate cause for the phenomenon was not resolved because it can be explained by seasonal changes in environmental factors (i.e., temperature and precipitation significantly increase in the summer) as well as developmental changes in leaf functional traits (Reich et al., 2014) and leaf growth and structural allocation patterns (Hudson et al., 2011). Therefore, the effect of leaf-age on leaf biomass allocation patterns remains uncertain due to a paucity of direct observation and phenotypic plasticity in response to environmental factors (Jiao et al., 2022; Westoby et al., 2022).
To address this area of uncertainty, we collected the leaves of the evergreen species Photinia serratifolia (Desfontaines) Kalkman in eight leaf-age groups from the spring to the winter of 2021. P. serratifolia was selected due to the new leaves of this species only emerged on the top of branches once a year, there is no need to worry about the investigated leaves being shadowed by another round of new leaves. We employed both the ratio and scaling approaches to investigate the relationships between lamina biomass and area. In addition, the scaling relationships between petiole biomass and other variables of interest (e.g., lamina biomass and lamina area) of the eight different age groups were examined because of their important mechanical and hydraulic functionalities (Sack et al., 2003; Filartiga et al., 2022; Li et al., 2022). The goal of this study was to answer one key question: does leaf-age affect leaf biomass allocation patterns?
Materials and methods
Leaf collection
Newly emerging leaves on three P. serratifolia trees, which all grew in Nanjing Forestry University campus (118°48′35″E, 32°4′67″N), Nanjing, Jiangsu Province, China, were tagged during late February and early March in 2021, and a total of 2698 leaves were collected from April to November, with over 300 leaves sampled for each month. Table 1 provides the sampling information for each leaf-age group, and Figure 1 presents representative examples of the investigated leaf-age groups. Supplementary Figure S1 presents examples of the investigated tagged leaves. The climatic information pertinent to Table 1 was collected from https://www.ncei.noaa.gov (Station Name: NANJING; ID: CHM00058238; 118.9°E, 31.93°N; 15 m a.s.l.).
Leaf measurements
Leaves were scanned to bitmap images at a 600-dpi resolution using a photo scanner (V550, Epson Indonesia, Batam, Indonesia). Adobe Photoshop CS6 (version: 13.0; Adobe, San Jose, CA, USA) was used to obtain lamina black-white images. The planar coordinates of each lamina boundary were extracted using an M-file based on MATLAB (version ≥ 2009a; MathWorks, Natick, MA, USA) developed by Shi et al. (2018). The lamina area of each leaf was calculated using the “bilat” function in the “biogeom” package (version 1.3.5; Shi et al., 2022) using R software (version 4.2.0; R Core Team, 2022). Lamina fresh mass (LFM) and dry mass (LDM), and petiole fresh mass (PFM) and dry mass (PDM) were measured using an electronic balance (Type: ML 204; Mettler Toledo Company, Greifensee, Switzerland). LFM and PFM were determined after drying laminae and petioles to constant weight in an oven (Type: XMTD8222; Jinghong Experimental Equipment Co., Ltd., Shanghai, China) at 80°C for 48 h before measurement.
Statistical methods
Preliminary regression analyses of the untransformed and log-transformed data indicated that scaling relationships statistically complied with power-law functions taking the general form of:
where and represent two interdependent variables of interest (e.g., lamina dry mass and petiole dry mass), and and are the slope and y-intercept of the log-log regression curve (i.e., the scaling exponent and the normalization constant, respectively), respectively. The log-transformed form of Equation 1 takes the linear form (Niklas, 1994; Niklas et al., 2007).
where , , and . Parameters and in Equation 2 were estimated using reduced major axis regression protocols (Niklas, 1994; Smith, 2009). The bootstrap percentile method (using 3000 bootstrapping replicates) was used to test the significance of the difference in the estimated scaling exponents of vs. between any two of the leaf-age groups (Efron and Tibshirani, 1993; Sandhu et al., 2011). Analysis of variance followed by Turkey’s honestly significant difference test with a 0.05 significance level were used to test for the significance of the differences in the quotients of LDM and A, PDM and A, LFM and A, and PFM and A among eight leaf-age groups (Hsu, 1996). Correlation tests were used to test whether there were significant correlations between the slopes (i.e., the scaling exponents) and leaf-age, between the intercepts and leaf-age, and between the ratios of leaf traits and leaf-age. All statistical analyses were performed using R (version 4.2.0; R Core Team, 2022).
Results
The numerical values of the scaling exponents of LDM vs. A were significantly greater than unity (i.e., α > 1.0) for each of the seven leaf-age groups, i.e., the lower bounds of the corresponding 95% confidence intervals (CIs) exceeded unity, with the exception of the first age group whose 95% CIs did include unity. Similarly, the scaling exponents of PDM vs. A were significantly greater than unity for the first and second age groups, whereas the 95% CIs of the scaling exponents for the remaining six age-groups included unity or the upper bounds of the 95% CIs were significant smaller than unity. The goodness of fit for LDM vs. A was significantly better than that of PDM vs. A. The r2 values for LDM vs. A were typically greater than 0.8, with two exceptions (i.e., 0.735 and 0.793), whereas the r2 for PDM vs. A was consistently smaller than 0.8 (Table 2).
Table 2 Fitted results for lamina dry mass (LDM) vs. lamina area (A) and petiole dry mass (PDM) vs. A in eight leaf-age groups.
The scaling exponents of LFM vs. A were significantly greater than unity for all leaf-age groups, with the lower bounds of the corresponding 95% CIs exceeding unity in each age group. The 95% CIs of the scaling exponents of PFM vs. A included unity for the third, seventh, and eighth leaf-age groups; for the remaining age groups, the scaling exponents of PFM vs. A were significantly greater than unity, i.e., the lower bounds of the 95% CIs were greater than unity. The goodness of fit for LFM vs. A was significantly greater than that of PFM vs. A. The r2 values for LFM vs. A were consistently above 0.9, whereas the r2 values for PFM vs. A did not exceed 0.8 (Table 3). The data indicated that the scaling exponents of LDM vs. A tended to be slightly greater than those of LFM vs. A, and that the exponents of PDM vs. A tended to be slightly smaller than those of PFM vs. A. Additionally, the goodness of fit for both lamina fresh mass and petiole fresh mass was better than for their corresponding dry mass counterparts (Tables 2, 3).
Table 3 Fitted results for lamina fresh mass (LFM) vs. lamina area (A) and petiole fresh mass (PFM) vs. A in eight leaf-age groups.
A comparison of the scaling exponents of LDM vs. A, PDM vs. A, LFM vs. A, and PFM vs. A indicated a significant correlation between the scaling exponents and leaf-age (|r| > 0.5, P < 0.05). Additionally, the temporal variation in the lamina mass vs. A scaling exponent and that in the petiole mass vs. A scaling exponent were significantly opposite (Figure 2). The scaling exponents of LDM vs. A, and LFM vs. A numerically increased with increasing leaf-age, and the positive correlation between the LDM vs. A scaling exponent and leaf-age was more robust than that of the LFM vs. A scaling exponent, as indicated by a larger correlation coefficient (Figures 2A, B). Conversely, the scaling exponents of PDM vs. A, and PFM vs. A decreased with increasing leaf-age, and the negative correlation between PFM vs. A scaling exponent and leaf-age was more robust than that between the PDM vs. A scaling exponent and leaf-age, as reflected by a lower correlation coefficient (Figures 2C, D). Although there were slight differences in the scaling exponents of LDM vs. A and LFM vs. A, the influence of leaf-age on those exponents was similar. The same held true for petiole scaling exponents that vary in the opposite direction (Figure 2).
Figure 2 Comparisons of the scaling exponents (α-values) of LDM vs. A (A), of LFM vs. A (B), of PDM vs. A (C), and of PFM vs. A (D). Each boxplot was obtained from 3000 bootstrap replications. The lowercase letters a–c on the top of each box denote the significance of the difference in the scaling exponents between any two leaf-ages at a 0.05 significance level. r is the correlation coefficient for the scaling exponents and leaf-age groups, and P is the significance test parameter. Leaf-age codes correspond to those in Table 1.
A comparison of the intercepts of PDM vs. A, LFM vs. A, and PFM vs. A also reveals a significant correlation between the intercept and leaf-age (|r| > 0.5, P < 0.05), with the exception of the intercept of LDM vs. A whose correlation coefficient with leaf-age is −0.481 (Figure 3). The trend in the variation of the intercept was opposite to the trend in the numerical value of the scaling exponent (Figures 3, 4). This inverse relationship has been reported previously for other scaling relationship but is not always biologically meaningful because it emerges as a mathematical artifact whenever multiple scaling relationships share a common point (, ) for which (Niklas and Hammond, 2019). Nevertheless, the observed variation in the numerical values of intercepts is meaningful because the values influence the quantities of biomass allocated to the construction of laminas given that . In addition, the correlations between the intercept of LDM vs. A and leaf-age, and between PDM vs. A and leaf-age (Figures 3A, C) were less robust than their fresh mass counterparts (Figures 3B, D), as shown by the numerically larger or smaller correlation coefficients.
Figure 3 Comparisons of the intercepts of LDM vs. A (A), of LFM vs. A (B), of PDM vs. A (C), and of PFM vs. A (D). Each boxplot was obtained from 3000 bootstrap replications. The lowercase letters a–c on the top of each box denote the significance of the difference in the scaling exponents between any two leaf-ages at a 0.05 significance level. r is the correlation coefficient for the scaling exponents and leaf-age groups, and P is the significance test parameter. Leaf-age codes correspond to those in Table 1.
Figure 4 Comparisons of the ratios of LDM to A (A), LFM to A (B), PDM to A (C), and PFM to A (D). The lowercase letters a–e on the top of each box denote the significance of the difference in the means between any two leaf-ages based on Tukey’s HSD test at a 0.05 significance level. The numbers above the whiskers represent the coefficients of variation (%). The horizontal solid lines represent the medians, and the asterisks within boxes represent the means. r is the correlation coefficient for the scaling exponents and leaf-age groups, and P is the significance test parameter. Leaf-age codes correspond to those in Table 1.
The ratios of LDM to A, LFM to A, PDM to A, and PFM to A were correlated positively with leaf-age (P < 0.05), and the correlation between the ratios of LFM to A and leaf-age was less statistically robust than the other correlations (Figure 4). Additionally, the positive correlation between ratios of dry mass (including lamina dry mass and petiole dry mass) to lamina area (Figures 4A, C) and leaf-age were more robust than their fresh mass counterparts (Figures 4B, D), as shown by the numerically larger correlation coefficients.
Discussion
This study has documented that there is a difference between the ratio and scaling approaches. As determined by the ratio approach, both the ratio of lamina biomass to lamina area and the ratio of petiole biomass to lamina area increased with increasing leaf-age. In contrast, the scaling approach reveals a statistically significant opposite trend in the scaling exponents dictating lamina biomass vs. lamina area compared to that of petiole biomass vs. lamina area. In this context, it is worth noting that the same ratio can be achieved by either a decrease in the denominator or an increase in the numerator, which makes the interpretation of a ratio potentially ambiguous when taken in isolation. Consequently, it is wise to consider the results of the ratio and scaling approach in tandem to eliminate any potential ambiguities in interpreting results such as those presented here.
Based on these two approaches, this study confirms that leaf-age is an important factor in determining the scaling relationships of leaf biomass allocation patterns. However, previous studies have attributed the variation of leaf biomass allocation patterns to environmental factors. For example, Pan et al. (2013) showed that the numerical value of α increases with altitude across 121 vascular plant species. Likewise, Thakur et al. (2019) observed that the numerical value of α increases with the degree of environmental stress (i.e., higher, drier, open habitats). These studies indicate that leaf development is responsive to local ambient conditions in addition to differing across species. In contrast, our data indicate that environmental differences during the eight months over which leaf were sampled had little or no observable effect. The correlation test of the relationships between the numerical values of the scaling exponents and normalization constants for any of the biomass allocation patterns for lamina and petiole (listed in Tables 2, 3) and the differences in monthly temperature and precipitation among the eight leaf samples (listed on Table 1) failed to reveal any statistically significant relationship (i.e., P >0.05). Although it is important to note in this context that the absence of evidence for an “environmental effect” is based on a small sample size of leaf-age (n = 8) drawn from only one taxon (see Supplementary Figures S2 and S3). In addition, the differences in the monthly temperature or precipitation recorded over the duration of sampling are arguably not extreme.
The foregoing caveats are important because environmental factors are known to affect leaf growth and interact with ontogeny (Hudson et al., 2011; Reich et al., 2014), such as the noticeable shift that is observed in the data collected between the second and fourth months, which might result from increasing temperature or precipitation (Figures 2, 3). Environmental factors can also result in temporal changes in leaf biomass allocation patterns over the lifespans of leaves. For example, older leaves can be shaded by new cohorts of leaves, which can alter light-interception, leaf temperatures, and rates of evapotranspiration. Notably, overshadowing was not significant in our P. serratifolia leaf samples, as the new leaves emerged in the spring and were sampled from the upper canopy. Nevertheless, given the limited research on the effects of leaf-age on leaf biomass allocation patterns, additional studies are required across diverse species and different environmental contexts before any definitive conclusions can be drawn (Figure 5). In addition, the ratio and scaling approaches for describing leaf biomass allocation are based on different assumptions and theories, making the choice of whether to use fresh or dry biomass as a representation of organ biomass uncertain. To address these uncertainties, we discuss these issues separately in the following three sections.
Figure 5 Summary of the effects of leaf-age on leaf anatomical traits and scaling exponents. The illustrated maturation processes are based on prior findings.
Ratio and scaling approaches, and the choice of fresh or dry as a representation of biomass
The difference in the results emerging from the use of the ratio and scaling approaches should be a key consideration when designing a research program (Cheng et al., 2015; Li et al., 2022). The ratio-based approach has the advantage that it reveals plant biomass allocation as a simple proportion or percentage between any two variables of interest (Wright et al., 2005; Poorter et al., 2012), whereas the scaling approach has the advantage of capturing any non-linear or linear relationship between two variables while adjusting for size or age, i.e., provides a simple way of summarizing proportional and size-dependent phenomena (Niklas, 1994; Milla and Reich, 2007; Niklas et al., 2007; Guo et al., 2022).
This study shows that the goodness of fit of the lamina fresh mass vs. lamina area scaling relationship is consistently more statistically robust than that of lamina dry mass vs. lamina area among each of the eight different leaf-age groups. The same is true for the petiole fresh mass vs. lamina area scaling relationship compared to that of petiole dry mass vs. lamina area. This phenomenon highlights the superiority of fresh mass over dry mass in describing leaf biomass allocation with the scaling approach, as previously demonstrated (see Huang et al., 2019; Liu et al., 2020; Guo et al., 2022).
Nevertheless, the foregoing results do not necessarily indicate that fresh mass should be substituted by dry mass under all circumstances. The choice should depend on the specific research question. For example, studies of nutrient cycling and ecosystem dynamics may benefit from the use of leaf dry biomass, as it reflects the long-term accumulation of plant biomass measured primarily in terms of carbon investments (Shipley et al., 2006). Dry biomass is also commonly used due to its practicality and ease of measurement, which makes it more suitable for large-scale studies and comparisons across different ecosystems (Wright et al., 2004; Poorter et al., 2012). On the other hand, studies of plant growth and development may require fresh biomass to accurately capture the plant’s physiological state (Shipley et al., 2005; Niinemets et al., 2007). Fresh biomass is also a more reasonable biomechanical trait because it reflects the loads that must be supported and because it correlates with turgor and therefore the stiffness of most primary plant tissues, particularly the hydrostatic tissues found in leaves (Niklas, 1991a; Niklas, 1991b; Niklas, 1992).
Influence of leaf-age on leaf biomass allocation
Leaf-age significantly influences leaf morphology and chemical composition, such as leaf thickness and cellulose, as leaves grow older (Mediavilla et al., 2011), and the photosynthetic capacity of leaves is reported to decline gradually with age, presumably due to changes in leaf chemical composition, leaf nitrogen content, and CO2 diffusion limitation (Kitajima et al., 1997; Day et al., 2001; Han et al., 2008; Zhang et al., 2008). Consequently, the influence of leaf-age on biomass allocation patterns appears to persist throughout leaf ontogeny until maturation (Figure 5), as demonstrated by the observed variation in the numerical values of scaling exponents for both lamina and petiole biomass allocation with respect to lamina area (Pantin et al., 2012; Jiao et al., 2022), although the variation of leaf biomass allocation in the senescence process is uncertain (Lim et al., 2007; Koyama, 2018). Specifically, a significant opposite trend in the scaling exponents between the two allocation patterns during ontogeny is evident, which highlights the age-related effects on leaf biomass allocation.
The scaling exponent of leaf dry mass vs. lamina area tends to be close to unity in the first leaf-age group (Table 2), probably because juvenile leaves tend to allocate a larger proportion of their biomass to support area expansion and expedite photosynthesis-related processes, such as the development of chloroplasts and other pigments, thereby maximizing their ability to capture light and nutrients (Westoby et al., 2000). Changes in the numerical values of scaling exponents likely reflect the accumulation of secondary cell wall materials and lignification (e.g., xylem and phloem fibers) during maturation. For example, the maturation of Laurus nobilis leaves can be described as a “hardening process” that entails the development of the vascular tissues, and the thickening and lignification of the cell walls in the bundle sheath extensions and the epidermis (Fasseas and Akoumianaki-Ioannidou, 2010; Schuetz et al., 2014). This phenomenology, which also occurs in petioles as a consequence of the vascular tissues within them, is partly responsible for the increase in lamina dry mass per unit area (LMA) and lamina thickness with age, which is reported to correlate with increases in the concentrations of Ca and lower concentrations of N, P, K, and Mg on a dry mass basis (Mediavilla et al., 2011). Typically, a higher LMA also corresponds with thicker palisade layers, which helps to maximize overall light absorption at greater depths within the mesophyll (Coble and Cavaleri, 2017). Therefore, the scaling exponents governing the relationship between lamina mass and area numerically increase (Niinemets, 2001; Westoby et al., 2002).
Despite similar developmental patterns, the vascular development within the lamina differs from that within the petiole, especially in larger leaves with higher hydraulic resistance (Sack et al., 2004; Pantin et al., 2012). The extra leaf lamina construction demands together with biomass constraints during growth (Niinemets et al., 2004; Niinemets et al., 2006) may necessitate the prioritization biomass allocation to the petiole for increased transport capacity (Filartiga et al., 2022). Regardless of the cause, the scaling exponents of petiole biomass versus lamina biomass numerically decrease with leaf-age, indicating a shift in biomass allocation from the lamina to the petiole as leaves reach their full maturity (see Supplementary Figure S4). This observation is consistent with prior studies, e.g., petiole mechanical stiffness increases with leaf-age, as demonstrated in a comparison between young (May) and mature (August) leaves of Populus tremuloides within a single growing season (Niklas, 1991c).
Taking leaf-age into account when assessing environmental effects
Temperature and precipitation are often reported to be the most important environmental factors affecting plant growth, thereby significantly affecting nutrient cycling, productivity, ecosystem fluxes, and other key plant and ecosystem processes (Kulmatiski and Beard, 2013; Hatfield and Prueger, 2015). Therefore, it is not surprising that leaf biomass allocation patterns to different organs or tissues are also significantly affected by temperature and other environmental conditions (Mediavilla et al., 2014; Reich et al., 2014). Two common leaf functional traits that respond to temperature are LDM and A, both of which are reported to increase with increasing temperature. LMA has been shown to respond plastically to abiotic variables such as light, temperature, water and nutrient availability, and atmospheric composition (Niinemets, 2001; Wright et al., 2005; Poorter et al., 2009). Thus, dynamic but predictable environmental conditions exert an influence on leaf biomass allocation patterns (Poorter et al., 2012).
However, as noted, the data reported here indicate that the “environmental effect” on the numerical value of biomass scaling exponents is statistically negligible, although the data presented here are limited both in leaf-age sample size (n = 8) and taxa (n = 1). The limited data also indicate that an “environmental effect” may be overestimated if leaf-age is not considered. At the site examined in this study, ambient temperature and precipitation gradually increase during the rainy season, especially in July and August. Likewise, the scaling exponents of lamina biomass vs. A numerically increased and the petiole biomass vs. A decreased in this period. However, after that period, the numerical values of scaling exponents are largely invariant and insensitive to changes in temperature and precipitation, even though the temperature and precipitation varied oppositely. This behavior is interpreted to indicate that scaling exponents are not responsive to changes in environmental conditions, which contradicts some previous studies (Pan et al., 2013; Thakur et al., 2019). However, additional research is clearly required to track in considerably greater detail the extent to which biomass vs. area scaling relationships respond to changes in environmental conditions, particularly during the early expansion and maturation of laminae (Westoby et al., 2022).
This caveat is particularly important because an increase in leaf biomass must be viewed as an ontogenetic process that is tied to lamina area expansion and structural maturation of both lamina and petiolar tissues. With increasing leaf-age, mechanical and hydraulic tissues are ontogenetically modified altering the biochemical and mechanical properties of cell primary and secondary walls (Wu et al., 2021). Research has also shown that the leaf circadian clock exerts control over biomass allocation patterns and that there is an asynchrony between the young leaves and older leaves (Pantin et al., 2012). Therefore, it is essential to accurately determine the age of leaves, when using either ratio or scaling methods to describe biomass allocation patterns. In passing, it is worth noting that leaf-age also affects other leaf traits. For example, maximum photosynthetic rates and nitrogen content are reported to decrease with leaf-age (Ishida et al., 1999; Han et al., 2008; Oikawa et al., 2008; Wang et al., 2015). Hence, a variety of important processes other than biomass allocation patterns are affected by leaf-age.
Conclusions
The data emerging from this study highlight an inverse proportional allocation between leaf petiole and lamina, and the significant influences of leaf-age on the relation of both lamina and petiole mass to the lamina area during leaf maturation. The data also yield different insights when subjected to two different approaches to quantity biomass allocation patterns (i.e., the ratio approach and the scaling approach). Based on these analyses, we recommend that (i) both the ratio and scaling approaches can be applied when analysing the relevant data, and (ii) that leaf-age should be included as an important variable of interest when evaluating biomass allocation patterns. Given the significant changes in the numerical values of petiole and lamina scaling exponents as a function of leaf maturation, this work demonstrates the necessity of an ontogenetic perspective when exploring biomass allocation patterns.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
XG: Formal analysis, Writing – original draft. JS: Writing – review & editing. PS: Formal analysis, Methodology, Writing – review & editing. YJ: Formal analysis, Investigation, Writing – original draft. QM: Writing – review & editing. JX: Writing – review & editing. KJN: Formal analysis, Methodology, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. JX was supported by the Evaluation on the biodiversity and carbon sequestration of typical wetland ecosystems in Jiangsu Province (No: LYKJ(2022)02).
Acknowledgments
We are thankful to Juan M. Losada and two expert reviewers for their invaluable comments on this manuscript. We also thank Alex Fajardo, David A. Ratkowsky, Kexin Yu, and Xiao Zheng for their valuable help in the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1322245/full#supplementary-material
References
Adams, W. W., III, Terashima, I. (2018). The Leaf: A Platform for Performing Photosynthesis. (Cham: Springer International Publishing). doi: 10.1007/978-3-319-93594-2
Cheng, D., Zhong, Q., Niklas, K. J., Ma, Y., Yang, Y., Zhang, J. (2015). Isometric scaling of above- and below-ground biomass at the individual and community levels in the understorey of a sub-tropical forest. Ann. Bot. 115, 303–313. doi: 10.1093/aob/mcu238
Coble, A. P., Cavaleri, M. A. (2017). Vertical leaf mass per area gradient of mature sugar maple reflects both height-driven increases in vascular tissue and light-driven increases in palisade layer thickness. Tree Physiol. 37, 1337–1351. doi: 10.1093/treephys/tpx016
Day, M. E., Greenwood, M. S., White, A. S. (2001). Age-related changes in foliar morphology and physiology in red spruce and their influence on declining photosynthetic rates and productivity with tree age. Tree Physiol. 21, 1195–1204. doi: 10.1093/treephys/21.16.1195
Díaz, S., Kattge, J., Cornelissen, J. H. C., Wright, I. J., Lavorel, S., Dray, S., et al. (2016). The global spectrum of plant form and function. Nature 529, 167–171. doi: 10.1038/nature16489
Efron, B., Tibshirani, R. J. (1993). An introduction to the bootstrap. (New York: Chapman and Hall/CRC).
Enquist, B. J., Kerkhoff, A. J., Stark, S. C., Swenson, N. G., McCarthy, M. C., Price, C. A. (2007). A general integrative model for scaling plant growth, carbon flux, and functional trait spectra. Nature 449, 218–222. doi: 10.1038/nature06061
Enquist, B. J., Niklas, K. J. (2002). Global allocation rules for patterns of biomass partitioning in seed plants. Science 295, 1517–1520. doi: 10.1126/science.1066360
Farrar, J. F., Jones, D. L. (2000). The control of carbon acquisition by roots. New Phytol. 147, 43–53. doi: 10.1046/j.1469-8137.2000.00688.x
Fasseas, C., Akoumianaki-Ioannidou, A. (2010). Some anatomical features of Laurus Nobilis leaves contributing to its suitability as ornamental plant in temperate climates. Acta Hortic. 885, 99–108. doi: 10.17660/ActaHortic.2010.885.13
Filartiga, A. L., Klimeš, A., Altman, J., Nobis, M. P., Crivellaro, A., Schweingruber, F., et al. (2022). Comparative anatomy of leaf petioles in temperate trees and shrubs: the role of plant size, environment and phylogeny. Ann. Bot. 129, 567–582. doi: 10.1093/aob/mcac014
Guo, X., Niklas, K. J., Li, Y., Xue, J., Shi, P., Schrader, J. (2022). Diminishing returns: A comparison between fresh mass vs. area and dry mass vs. area in deciduous species. Front. Plant Sci. 13, 832300. doi: 10.3389/fpls.2022.832300
Han, Q., Kawasaki, T., Nakano, T., Chiba, Y. (2008). Leaf-age effects on seasonal variability in photosynthetic parameters and its relationships with leaf mass per area and leaf nitrogen concentration within a Pinus densiflora crown. Tree Physiol. 28, 551–558. doi: 10.1093/treephys/28.4.551
Hatfield, J. L., Prueger, J. H. (2015). Temperature extremes: Effect on plant growth and development. Weather Clim. Extrem. 10, 4–10. doi: 10.1016/j.wace.2015.08.001
Huang, W., Ratkowsky, D. A., Hui, C., Wang, P., Su, J., Shi, P. (2019). Leaf fresh weight versus dry weight: Which is better for describing the scaling relationship between leaf biomass and leaf area for broad-leaved plants? Forests 10, 256. doi: 10.3390/f10030256
Hudson, J. M. G., Henry, G. H. R., Cornwell, W. K. (2011). Taller and larger: shifts in Arctic tundra leaf traits after 16 years of experimental warming: Arctic leaf trait responses to warming. Glob. Change Biol. 17, 1013–1021. doi: 10.1111/j.1365-2486.2010.02294.x
Ishida, A., Uemura, A., Koike, N., Matsumoto, Y., Hoe, A. L. (1999). Interactive effects of leaf age and self-shading on leaf structure, photosynthetic capacity and chlorophyll fluorescence in the rain forest tree, Dryobalanops aromatica. Tree Physiol. 19, 741–747. doi: 10.1093/treephys/19.11.741
Jahnke, S., Menzel, M. I., Van Dusschoten, D., Roeb, G. W., Bühler, J., Minwuyelet, S., et al. (2009). Combined MRI-PET dissects dynamic changes in plant structures and functions. Plant J. 59, 634–644. doi: 10.1111/j.1365-313X.2009.03888.x
Jiao, Y., Niklas, K. J., Wang, L., Yu, K., Li, Y., Shi, P. (2022). Influence of leaf age on the scaling relationships of lamina mass vs. area. Front. Plant Sci. 13, 860206. doi: 10.3389/fpls.2022.860206
Kitajima, K., Mulkey, S. S., Wright, S. J. (1997). Decline of photosynthetic capacity with leaf age in relation to leaf longevities for five tropical canopy tree species. Am. J. Bot. 84, 702–708. doi: 10.2307/2445906
Koyama, T. (2018). A hidden link between leaf development and senescence. Plant Sci. 276, 105–110. doi: 10.1016/j.plantsci.2018.08.006
Kulmatiski, A., Beard, K. H. (2013). Woody plant encroachment facilitated by increased precipitation intensity. Nat. Clim. Change 3, 833–837. doi: 10.1038/nclimate1904
Li, Y., Shi, P., Niinemets, Ü., Song, Y., Yu, K., Schrader, J., et al. (2022). Diminishing returns among lamina fresh and dry mass, surface area, and petiole fresh mass among nine Lauraceae species. Am. J. Bot. 109, 377–392. doi: 10.1002/ajb2.1812
Lim, P. O., Kim, H. J., Gil Nam, H. (2007). Leaf senescence. Annu. Rev. Plant Biol. 58, 115–136. doi: 10.1146/annurev.arplant.57.032905.105316
Liu, M., Niklas, K. J., Niinemets, Ü., Höelscher, D., Chen, L., Shi, P. (2020). Comparison of the scaling relationships of leaf biomass versus surface area between spring and summer for two deciduous tree species. Forests 11, 1010. doi: 10.3390/f11091010
Lockhart, B. R., Hodges, J. D., Gardiner, E. S., Ezell, A. W. (2003). Photosynthate distribution patterns in cherrybark oak seedling sprouts. Tree Physiol. 23, 1137–1146. doi: 10.1093/treephys/23.16.1137
Mediavilla, S., González-Zurdo, P., García-Ciudad, A., Escudero, A. (2011). Morphological and chemical leaf composition of Mediterranean evergreen tree species according to leaf age. Trees 25, 669–677. doi: 10.1007/s00468-011-0544-z
Mediavilla, S., Herranz, M., Gonzalez-Zurdo, P., Escudero, A. (2014). Ontogenetic transition in leaf traits: a new cost associated with the increase in leaf longevity. J. Plant Ecol. 7, 567–575. doi: 10.1093/jpe/rtt059
Milla, R., Reich, P. B. (2007). The scaling of leaf area and mass: the cost of light interception increases with leaf size. Proc. R. Soc B Biol. Sci. 274, 2109–2115. doi: 10.1098/rspb.2007.0417
Minchin, P. E. H., Thorpe, M. R. (2003). Using the short-lived isotope 11C in mechanistic studies of photosynthate transport. Funct. Plant Biol. 30, 831–841. doi: 10.1071/FP03008
Nakamura, Y., Hashimoto, H. (1988). Characteristics of photosynthate partitioning during chloroplast development in Avena Leaves. Plant Physiol. 87, 458–462. doi: 10.1104/pp.87.2.458
Niinemets, Ü. (2001). Global-scale climatic control of leaf dry mass per area, density, and thickness in trees and shrubs. Ecology 82, 453–469. doi: 10.1890/0012-9658(2001)082[0453:GSCCOL]2.0.CO;2
Niinemets, Ü., Al Afas, N., Cescatti, A., Pellis, A., Ceulemans, R. (2004). Petiole length and biomass investment in support modify light interception efficiency in dense poplar plantations. Tree Physiol. 24, 141–154. doi: 10.1093/treephys/24.2.141
Niinemets, Ü., Portsmuth, A., Tena, D., Tobias, M., Matesanz, S., Valladares, F. (2007). Do we underestimate the importance of leaf size in plant economics? Disproportional scaling of support costs within the spectrum of leaf physiognomy. Ann. Bot. 100, 283–303. doi: 10.1093/aob/mcm107
Niinemets, Ü., Portsmuth, A., Tobias, M. (2006). Leaf size modifies support biomass distribution among stems, petioles and mid-ribs in temperate plants. New Phytol. 171, 91–104. doi: 10.1111/j.1469-8137.2006.01741.x
Niklas, K. J. (1991a). Biomechanical responses of Chamaedorea and Spathiphyllum petioles to tissue dehydration. Ann. Bot. 67, 67–76. doi: 10.1093/oxfordjournals.aob.a088101
Niklas, K. J. (1991b). Effects of tissue volume and location on the mechanical consequences of dehydration of petioles. Am. J. Bot. 78, 361–369. doi: 10.1002/j.1537-2197.1991.tb15198.x
Niklas, K. J. (1991c). The elastic moduli and mechanics of Populus tremuloides (Salicaceae) petioles in bending and torsion. Am. J. Bot. 78, 989–996. doi: 10.1002/j.1537-2197.1991.tb14503.x
Niklas, K. J. (1992). Plant biomechanics: an engineering approach to plant form and function (Chicago: The University of Chicago press).
Niklas, K. J. (1994). Plant allometry: the scaling of form and process (Chicago: The University of Chicago Press).
Niklas, K. J. (1999). A mechanical perspective on foliage leaf form and function. New Phytol. 143, 19–31. doi: 10.1046/j.1469-8137.1999.00441.x
Niklas, K. J., Cobb, E. D., Niinemets, Ü., Reich, P. B., Sellin, A., Shipley, B., et al. (2007). “Diminishing returns” in the scaling of functional leaf traits across and within species groups. Proc. Natl. Acad. Sci. 104, 8891–8896. doi: 10.1073/pnas.0701135104
Niklas, K. J., Cobb, E. D., Spatz, H.-C. (2009). Predicting the allometry of leaf surface area and dry mass. Am. J. Bot. 96, 531–536. doi: 10.3732/ajb.0800250
Niklas, K. J., Enquist, B. J. (2002). Canonical rules for plant organ biomass partitioning and annual allocation. Am. J. Bot. 89, 812–819. doi: 10.3732/ajb.89.5.812
Niklas, K. J., Hammond, S. T. (2019). On the interpretation of the normalization constant in the scaling equation. Front. Ecol. Evol. 6. doi: 10.3389/fevo.2018.00212
Oikawa, S., Hikosaka, K., Hirose, T. (2008). Does leaf shedding increase the whole-plant carbon gain despite some nitrogen being lost with shedding? New Phytol. 178, 617–624. doi: 10.1111/j.1469-8137.2008.02415.x
Pan, S., Liu, C., Zhang, W., Xu, S., Wang, N., Li, Y., et al. (2013). The scaling relationships between leaf mass and leaf area of vascular plant species change with altitude. PloS One 8, e76872. doi: 10.1371/journal.pone.0076872
Pantin, F., Simonneau, T., Muller, B. (2012). Coming of leaf age: control of growth by hydraulics and metabolics during leaf ontogeny. New Phytol. 196, 349–366. doi: 10.1111/j.1469-8137.2012.04273.x
Pasini, D., Mirjalili, V. (2006). “The optimized shape of a leaf petiole” in Design and Nature III: Comparing Design in Nature with Science and Engineering (The New Forest, UK: WIT Press) 87, 35–45. doi: 10.2495/DN060041
Poorter, H., Jagodzinski, A. M., Ruiz-Peinado, R., Kuyah, S., Luo, Y., Oleksyn, J., et al. (2015). How does biomass distribution change with size and differ among species? An analysis for 1200 plant species from five continents. New Phytol. 208, 736–749. doi: 10.1111/nph.13571
Poorter, H., Niinemets, Ü., Poorter, L., Wright, I. J., Villar, R. (2009). Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis. New Phytol. 182, 565–588. doi: 10.1111/j.1469-8137.2009.02830.x
Poorter, H., Niklas, K. J., Reich, P. B., Oleksyn, J., Poot, P., Mommer, L. (2012). Biomass allocation to leaves, stems and roots: meta-analyses of interspecific variation and environmental control. New Phytol. 193, 30–50. doi: 10.1111/j.1469-8137.2011.03952.x
R Core Team (2022) R: A language and environment for statistical computing. Available at: https://www.r-project.org/.
Rascher, U., Nedbal, L. (2006). Dynamics of photosynthesis in fluctuating light. Curr. Opin. Plant Biol. 9, 671–678. doi: 10.1016/j.pbi.2006.09.012
Reich, P. B., Luo, Y., Bradford, J. B., Poorter, H., Perry, C. H., Oleksyn, J. (2014). Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc. Natl. Acad. Sci. 111, 13721–13726. doi: 10.1073/pnas.1216053111
Sack, L., Cowan, P. D., Jaikumar, N., Holbrook, N. M. (2003). The ‘hydrology’ of leaves: co-ordination of structure and function in temperate woody species. Plant Cell Environ. 26, 1343–1356. doi: 10.1046/j.0016-8025.2003.01058.x
Sack, L., Scoffoni, C., McKown, A. D., Frole, K., Rawls, M., Havran, J. C., et al. (2012). Developmentally based scaling of leaf venation architecture explains global ecological patterns. Nat. Commun. 3, 837. doi: 10.1038/ncomms1835
Sack, L., Streeter, C. M., Holbrook, N. M. (2004). Hydraulic analysis of water flow through leaves of sugar maple and red oak. Plant Physiol. 134, 1824–1833. doi: 10.1104/pp.103.031203
Sandhu, H. S., Shi, P., Kuang, X., Xue, F., Ge, F. (2011). Applications of the bootstrap to insect physiology. Fla. Entomol. 94, 1036–1041. doi: 10.1653/024.094.0442
Schuetz, M., Benske, A., Smith, R. A., Watanabe, Y., Tobimatsu, Y., Ralph, J., et al. (2014). Laccases direct lignification in the discrete secondary cell wall domains of protoxylem. Plant Physiol. 166, 798–807. doi: 10.1104/pp.114.245597
Shi, P., Gielis, J., Quinn, B. K., Niklas, K. J., Ratkowsky, D. A., Schrader, J., et al. (2022). ‘biogeom’: An R package for simulating and fitting natural shapes. Ann. N. Y. Acad. Sci. 1516, 123–134. doi: 10.1111/nyas.14862
Shi, P., Ratkowsky, D., Li, Y., Zhang, L., Lin, S., Gielis, J. (2018). A general leaf area geometric formula exists for plants—evidence from the simplified Gielis equation. Forests 9, 714. doi: 10.3390/f9110714
Shipley, B., Lechowicz, M. J., Wright, I., Reich, P. B. (2006). Fundamental trade-offs generating the worldwide leaf economics spectrum. Ecology 87, 535–541. doi: 10.1890/05-1051
Shipley, B., Vile, D., Garnier, E., Wright, I. J., Poorter, H. (2005). Functional linkages between leaf traits and net photosynthetic rate: reconciling empirical and mechanistic models. Funct. Ecol. 19, 602–615. doi: 10.1111/j.1365-2435.2005.01008.x
Smith, R. J. (2009). Use and misuse of the reduced major axis for line-fitting. Am. J. Phys. Anthropol. 140, 476–486. doi: 10.1002/ajpa.21090
Thakur, D., Rathore, N., Chawla, A. (2019). Increase in light interception cost and metabolic mass component of leaves are coupled for efficient resource use in the high altitude vegetation. Oikos 128, 254–263. doi: 10.1111/oik.05538
Wang, Z., Lu, J., Yang, M., Yang, H., Zhang, Q. (2015). Stoichiometric characteristics of carbon, nitrogen, and phosphorus in leaves of differently aged Lucerne (Medicago sativa) stands. Front. Plant Sci. 6, 01062. doi: 10.3389/fpls.2015.01062
Westoby, M., Falster, D. S., Moles, A. T., Vesk, P. A., Wright, I. J. (2002). Plant ecological strategies: some leading dimensions of variation between species. Annu. Rev. Ecol. Syst. 33, 125–159. doi: 10.1146/annurev.ecolsys.33.010802.150452
Westoby, M., Schrader, J., Falster, D. (2022). Trait ecology of startup plants. New Phytol. 235, 842–847. doi: 10.1111/nph.18193
Westoby, M., Warton, D., Reich, P. B. (2000). The time value of leaf area. Am. Nat. 155, 649–656. doi: 10.1086/303346
Wright, I. J., Reich, P. B., Cornelissen, J. H. C., Falster, D. S., Groom, P. K., Hikosaka, K., et al. (2005). Modulation of leaf economic traits and trait relationships by climate. Glob. Ecol. Biogeogr. 14, 411–421. doi: 10.1111/j.1466-822x.2005.00172.x
Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. Nature 428, 821–827. doi: 10.1038/nature02403
Wu, W., Du, K., Kang, X., Wei, H. (2021). The diverse roles of cytokinins in regulating leaf development. Hortic. Res. 8, 118. doi: 10.1038/s41438-021-00558-3
Keywords: diminishing returns, leaf mass per unit area, ontogeny, scaling growth, trade-offs
Citation: Guo X, Schrader J, Shi P, Jiao Y, Miao Q, Xue J and Niklas KJ (2023) Leaf-age and petiole biomass play significant roles in leaf scaling theory. Front. Plant Sci. 14:1322245. doi: 10.3389/fpls.2023.1322245
Received: 16 October 2023; Accepted: 06 December 2023;
Published: 21 December 2023.
Edited by:
Juan M. Losada, Spanish National Research Council (CSIC), SpainReviewed by:
Mauro Maver, Free University of Bozen-Bolzano, ItalyUlises Rosas, National Autonomous University of Mexico, Mexico
Copyright © 2023 Guo, Schrader, Shi, Jiao, Miao, Xue and Niklas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peijian Shi, cGpzaGlAbmpmdS5lZHUuY24=; Karl J. Niklas, a2puMkBjb3JuZWxsLmVkdQ==