- 1College of Geographical Sciences, Fujian Normal University, Fuzhou, China
- 2Fujian Provincial Key Laboratory of Plant Ecophysiology, Fujian Normal University, Fuzhou, China
- 3Key Laboratory of Humid Subtropical Eco-Geographical Process, Ministry of Education, Fuzhou, China
- 4School of Resources and Environment, Anqing Normal University, Anqing, China
Foliage leaves are the primary photosynthetic organ of the majority of vascular plants, and their area vs. biomass scaling relationships provide valuable insights into the capacity and investment in light interception, which is critical to plant growth and performance. The “diminishing returns” hypothesis (DRH), which is based primarily on data from gymnosperms and angiosperms, posits that leaf (lamina) area scales with leaf dry mass. on average with a scaling exponent less than 1.0. However, it remains uncertain whether DRH applies to ferns or whether ecological factors affect the scaling exponents governing fern leaf morphometrics. To address this issue, 182 individuals of 28 subtropical ferns species were studied at low, medium, and high elevations (i.e., 600 m, 900 m, and 1200 m, respectively) in Mount Wuyi National Park, Jiangxi Province, China. The scaling relationships between leaf area and leaf biomass for individual and total leaf of ferns at different elevations were examined by using standardized major axis regression protocols. Analyses of the 28 fern species (using Blomberg K-value protocols) indicated no phylogenetic biases among the species compositions of the three different elevations. In addition, at the individual plant level, individual leaf area (ILA) did not differ significantly among the three different elevations (P > 0.05). However, individual leaf mass (ILM) was significantly higher at 900m than at 1200m (P < 0.05), resulting in a significantly higher leaf mass per area (LMA) at the 900m elevation than at the 600m and 1200m elevations (P < 0.05). The ILA and ILM at the 900m elevation were significantly higher than at the 600m elevation (P < 0.05). At the species level, ILA and ILM did not differ significantly among the three elevations (P > 0.05). The total leaf area per individual (TLA) did not differ significantly across the different elevations (P > 0.05). However, total leaf mass per individual (TLM) did differ significantly (P < 0.05). At the individual plant level, the scaling exponents for ILA vs. ILM and TLA vs. TLM at the three different elevations were all significantly less than 1.0 (P < 0.05), which was consistent with the DRH. At the species level, the scaling exponents for the ILA vs. ILM were significantly less than 1.0 at the middle and high elevations, but not at the low elevation. The scaling exponents of the TLA and TLM were numerically highest in the middle elevation, and all were less than 1.0 for the three elevations. These results indicate that the scaling relationships of leaf area versus mass of subtropical ferns at different elevations support the DRH hypothesis. The study further informs our understanding of the resource allocation strategies of an ancient and diverse plant lineage.
1 Introduction
The relationship between leaf (lamina) biomass and area not only reflects the light-harvesting capability of most land plants, but also their environmental adaptability (Kleiman and Aarssen, 2007; Huang et al., 2020; Guo et al., 2022). Therefore, this scaling relationship has important implications for plant growth and community productivity (Liu et al., 2020).
Typically, dry leaf biomass (ML) characterizes the cost of leaf construction, whereas leaf area (AL) characterizes the potential for photosynthetic gain. Prior work has shown that AL is statistically significantly related to ML and takes the form of a power function, i.e., AL =β MLα, where α is the scaling exponent and β is the scaling (normalization) constant. When α < 1.0, leaf area increases more slowly relative to increasing leaf biomass, a phenomenon called “diminishing returns”. When α > 1.0, there is a gain in leaf area relative to increasing leaf biomass. An isometric (“break even”) relationship occurs when α = 1.0 (Niklas et al., 2007). By determining the leaf area vs. leaf biomass scaling relationships of six functional species groups and 19 individual species, Niklas et al. (2007) confirmed the “diminishing returns” hypothesis (DRH) wherein α < 1.0. Other studies have explored the DRH for plants at different geographic scales and different regions (Shipley, 1995; Pearcy et al., 2005). For example, Shi et al. (2019) showed that leaves area vs. dry mass scaling relationship were in accordance with the DRH in all examined vines leaves. Sun et al. (2017) showed that all five species of bamboo leaves are consistent with the DRH. More importantly, the total leaf area-biomass of whole plant should follow the DRH as the individual leaf area- biomass. Indeed, Zhu et al. (2019) found that DRH applied to the scaling relationship between total leaf area and total leaf biomass in current twigs of 64 subtropical woody species.
However, plant growth and development are affected by the environment (Zhang et al., 2020), and plant biomass exhibits the same or different scaling relationship at different elevations likely in response to hydrothermal conditions (Korner et al., 2001). Thus, it is reasonable to surmise that the scaling of leaf area with respect to leaf biomass may differ across gradients of temperature, water availability, and elevation. Yang et al. (2012) found that leaf size is smaller in species adapted to high elevations compared to species adapted to lower elevations, and report that an elevation gradient can significantly affect the trade-off between leaf area and leaf number. However, Sun et al. (2017) showed that the scaling exponents for the leaf area vs. dry mass scaling relationship are less than 1.0 within and across five species of bamboo and are insensitive to light conditions and elevation. Geng et al. (2019) found that the slope of leaf area and leaf biomass of Disanthus cercidifolius was significantly different from 1.0 at low altitudes, whereas at mid and high altitudes the slope was not significantly different from 1.0, i.e., an isotropic relationship. These different findings indicate that differences in elevation can have an effect on the scaling relationship of plants under different circumstances.
Importantly, the DRH is based primarily on data drawn from gymnosperms and angiosperms. Consequently, little is known about the scaling of non-seed plant species, such as ferns, particularly those growing in tropical forest ecosystems. Ferns represent a unique and extremely ancient evolutionary group of vascular plants (Bierhorst, 1971; Page, 2002) and have an important position and role in the structure, function, and dynamics of forest ecosystems (Li et al., 2000). The response of ferns to external natural conditions is highly sensitive and selective (Schneider et al., 2004; Jin et al., 2021), significantly influencing the dynamic processes of forest communities, and can be important for the occurrence and succession of communities, nutrient cycling, and energy flow (Yan et al., 2004). However, as noted, most current studies on plant functional traits have been focused on gymnosperm and angiosperms species with self-supporting stems (Westoby et al., 2002; Shipley et al., 2006). Far fewer studies have dwelt on the relationship between leaf area and leaf biomass in rhizomatous species. Fern fronds are often the only above-ground part of the plant body and are therefore the primary photosynthetic organs (Mao et al., 2012). Therefore, studying the relationship between frond area and biomass and its adaptation to different environments can inform our understanding of the adaptation strategies of an ancient and diverse lineage growing under different conditions.
The subtropical region of China has high forest cover, soil fertility, and biodiversity. It therefore provides a rich natural environmental context to study fern biology. The typical middle subtropical forest ecosystem of Wuyi Mountain was selected for study because of its huge variety of ferns that are widely distributed. In this study, ferns at different elevations were selected to determine the leaf area vs. biomass scaling relationship and to determine whether there are differences in the numerical values of the scaling exponents of the leaf area vs. leaf biomass scaling relationship at different elevations.
2 Materials and methods
2.1 Study area
The study area is located in the Wuyi Mountain National Park in southeastern of China. The area has a subtropical monsoon climate, and is the largest and most completely preserved meso-subtropical forest ecosystem in the world. The average annual precipitation is 1486-2150 mm and is concentrated in April-June. The average annual temperature is 8.5-18°C. The highest peak in the reserve is Huanggang Mountain. As the elevation decreases, the soil vertical zones are mountain meadow soil zone, yellow soil zone, yellow-red soil zone, and red soil zone (Chen, 2000).
2.2 Sample collection and index determination
Twenty-eight plants with healthy and fully mature leaves were selected for each species at elevations of 600 m, 900 m, and 1200 m, and we verified on website (http://www.sp2000.org.cn/) that all fern species included in the sample were native. The individual plants were excavated and placed into coded sealed bags and then brought back to the laboratory. Leaves were removed from rhizomes and their lamina were scanned using a scanner (EPSON V39, Japan). The scanned images were processed using ImageJ software to obtain leaf area data. The scanned leaf samples were then dried in an oven at 75°C for 48 hours to constant weight and weighed using an electronic scale with an accuracy of 0.001 g. LMA was calculated as the leaf dry weight per leaf area (mg/mm²). The leaf (lamina) area and dry weight were measured and calculated the leaf (lamina) area and dry weight for each fern samples collected. The ILM at individual plant level was calculated as the total leaf mass of each sample divided by the number of leaves. Similarly, the total leaf mass per individual (TLM) was calculated as the sum of the leaf biomass of each fern plant harvested. At the species level, the ILM was calculated as the average value of the TLM of the sample divided by the number of leaves. The TLM was calculated as the average of the sum of the total leaf biomass of each fern harvested at the species level. Similarly, the individual leaf area (ILA) and total leaf area per individual (TLA) were calculated in the same way.
3 Analysis
The data were processed using Excel 2016 software. The mean and standard errors were calculated for the leaf area and dry weight of individual plants using SPSS 19.0. The differences between leaf area, biomass, and mass per area at different elevations were analyzed using one-way ANOVA.
The K-value method proposed by Blomberg et al. (2003) was used to measure the strength of the phylogenetic signal for continuous functional traits. It detects correlations between functional traits and species evolutionary history, and allows comparisons between traits and between phylogenetic trees. If the K-value is greater than 1, the functional traits under investigation show a stronger phylogenetic signal than in the Brownian motion model. If the K-value is less than 1, the functional traits show a weaker phylogenetic signal than the Brownian motion model. Under Brownian motion, evolutionary changes are simply added to the values that occurred in the previous generation or at the last node on the phylogenetic tree. Consequently, lineages that have recently diverged will exhibit a higher degree of similarity (on average), compared to the lineages that are more distantly related (Blomberg et al., 2003). We built a phylogenetic tree using the “V. PhyloMaker2” package in the R(4.0.3)(Jin and Qian, 2022), and the K values in this study were analyzed in R using the phylosignal() function in the package “picante” (Kembel et al., 2010). The data were log10-transformed and leaf area and biomass at different elevations were analyzed in R (4.0.3) using the “smatr” package to determine the scaling relationships among the variables of interest, i.e., AL =β MLα, where β is the scaling (normalization) constant, α is the scaling exponent, and AL and ML are leaf area and leaf biomass, respectively. In addition, we assessed the scaling relationships of shared species across three different elevations. To compare the exponent associated with each of the three elevations, we used one-way ANOVA in SPSS 19.0. The analysis was performed at the species level using the mean of each variable of the species. Origin 2019 was used to plot the data.
4 Results
4.1 Characteristics of fern leaf area and leaf biomass across different elevations
The data indicated that the phylogenetic differences of individual leaf mass (ILM), individual leaf area (ILA), and leaf mass per area (LMA) of the plants from the three elevations of ferns were not significant (S1; Table 1).
In addition, the numerical values of the scaling exponents of two fern species (Phegopteris decursive-pinnata and Dicranopteris pedata) growing at each of the three elevations were used to assess whether exponents at the species level were affected by elevation (Table 2). The scaling exponents of ILA and ILM for P. decursive-pinnata were 0.82, 0.60, 0.78 at 600m, 900m, and 1200m, respectively (Table 2), and the respective scaling exponents for D. dichotoma were 1.60, 0.90, 0.88 at the three elevations (Table 2). Based on their 95% Cis, these results indicated that the values of the scaling exponents of these two species were, by in large, insensitive to elevation, although one of the scaling relationships was statistically insignificant (e.g., D. dichotoma at 600 m; Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-t002.jpg)
Table 2 (S)MATR reduced major axis regression slopes and y-intercepts (α and log β, respectively) for log10-tranformed data for individual leaf area and leaf biomass for Phegopteris decursive-pinnata and Dicranopteris pedata across different elevations.
The leaf biomass of individual species at different elevations was species-specific and differed but not significantly as a function of elevation (Table 3). The mean ILM was 149.75 ± 335.53 mg at the 900m elevation, which was significantly higher than the ILM of 57.41 ± 66.01 mg at the 1200m elevation (Figure 1A). However, no statistically significant difference in ILA was detected across elevations (P > 0.05, Figure 1B), although ILM varied significantly within each elevation. Specifically, the difference between the maximum and minimum values of ILM was 37, 269, and 22 times in three elevations, respectively, whereas LMA varied significantly within each elevation. Specifically, LMA was 15, 11 and 2 times between the maximum and minimum in 600m, 900m and 1200m, respectively (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-t003.jpg)
Table 3 Mean (± SE) number of leaves per plant (N), average biomass per leaf, average area per leaf, and average specific leaf weight for each of 38 species of subtropical ferns.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-g001.jpg)
Figure 1 Leaf traits of 28 species of subtropical ferns at different elevations including: (A) individual leaf mass, (B) individual leaf area, (C) total leaf mass per individual, (D) total leaf area per individual, (E) leaf mass per area. Letters above each plot indicate statistical differences (P < 0.05). No letters marked mean no statistical differences [P > 0.05, see (B)].
The average value of total leaf mass per individual (TLM) at 900m elevation was 428.46 ± 372.10 mg, which was significantly higher than that at 600m elevation (244.10 ± 211.40 mg) and at 1200 elevation (222.71 ± 142.42 mg, P < 0.05, Figure 1C). Total leaf area per individual (TLA) at 600m elevation (21664.52 ± 15145.01 mm²) was significantly lower than that of 900 elevation (29735.54 ± 17287.02 mm², P < 0.05, Figure 1D). The LMA of ferns growing at 900m elevation (0.015 ± 0.013 mg/mm²) was significantly higher than that of ferns at 600m elevation (0.010 ± 0.003 mg/mm²) and at 1200m elevation (0.008 ± 0.002 mg/mm²; P < 0.05, Figure 1E).
4.2 The scaling relationship between fern leaf area and leaf biomass at the individual plant level
At the individual plant level, species manifested a significant ILA vs. ILM scaling relationship at each of the three elevations (Figure 2A). The scaling exponents of ILA vs. ILM at 600 m,900 m, and 1200 m were 0.90 (95% CI = 0.85-0.94), 0.76 (95% CI = 0.73-0.80), and 0.83 (95% CI = 0.71-0.97), respectively. All were significantly less than 1.0 (P < 0.01), thereby validating the DRH (Table 4). TLA vs. TLM scaling relationship was also statistically significant for each of the three elevations. The scaling exponents of TLA vs. TLM at 600m, 900, and 1200m were 0.84 (95% CI = 0.78-0.89), 0.61 (95% CI = 0.55-0.68), and 0.71 (95% CI = 0.65-0.77), respectively. The scaling exponent at the different elevations were significantly less than 1.0 (P < 0.01, Figure 2B). At the different elevations, the scaling (normalization) constant of individual leaves increased with increasing elevation, while the scaling (normalization) constant at the whole plant level showed a tendency to increase and then decrease with increasing elevation (Table 4). These results indicate that the DRH at the level of individual leaves and whole plants leaves remained valid.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-g002.jpg)
Figure 2 The scaling relationship between leaf area and biomass at the individual level. (A) Individual leaf area (ILA) versus individual leaf mass (ILM). (B) Total leaf area per individual (TLA) versus total leaf mass per individual (TLM).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-t004.jpg)
Table 4 (S)MATR reduced major axis regression slopes and y-intercepts (α and logβ, respectively) for log10-tranformed data for the individual leaf area vs. leaf biomass scaling relationships of individual plants and at the species level for different elevations.
4.3 The scaling relationship between leaf area and leaf biomass at the species level
At the species level, the ferns investigated in this study had an isometric scaling relationship at 600m elevation, but had an allometric ILA vs. ILM scaling relationship at 900m and 1200m elevation Figure 3A. The isometric exponent of ILA vs. ILM at 600m elevation was 0.92 (95% CI = 0.83-1.02), which was not significantly different from 1.0 (P > 0.05; Table 4). The scaling exponent of ILA vs. ILM at 900m and 1200m elevations were 0.77 (95% CI = 0.68-0.87) and 0.81 (95% CI = 0.70-0.93), respectively, both of which were significantly less than 1.0 (P < 0.05). In addition, the ILA vs. ILM scaling relationship did not differ significantly among the three different elevations (P > 0.05).
At the species level, the scaling exponents of TLA vs. TLM were 0.86 (95% CI = 0.74-0.99), 0.59 (95% CI = 0.45-0.77), and 0.67 (95% CI = 0.53-0.84) for the 600 m, 900 m, and 1200 m, respectively. Each was significantly less than 1.0 (P < 0.01, Figure 3B; Table 4). However, each exponent was statistically different at the three different elevations. The TLM was not significantly different among the three elevations (P > 0.05), and LMA were different among elevations (P < 0.05), the TLM was significantly different between middle and high elevations (P < 0.05). The TLM at 900 elevation (442.66 ± 88.33 mg) was significantly higher than 1200m elevation (117.14 ± 35.32 mg, P < 0.05). The ILA (75.23 ± 29.40mm²) was significantly higher in the 900m than that at 600m elevation (71.30 ± 17.76 mm², P < 0.05, Figure 3B). Likewise, the scaling (normalization) constant of individual leaves increased with increasing elevation, whereas the scaling (normalization) constant at the whole plant level increased and then decreased with increasing elevation (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1187704/fpls-14-1187704-HTML/image_m/fpls-14-1187704-g003.jpg)
Figure 3 The scaling relationship between leaf area and biomass at the species level. (A) Individual leaf area (ILA) versus individual leaf mass (ILM). (B) Total leaf area per individual (TLA) versus total leaf mass per individual (TLM).
5 Discussion
The scaling relationship between leaf area and leaf biomass is important for understanding the investment trade-off strategy of plants because it reflects a ‘profit’-’investment’ strategy for photosynthesis and thus plant growth (Niklas et al., 2007). Compared with seed plants, there is a paucity of research on the functional characteristics of fern leaves despite the antiquity of this pteridophyte lineage and its ecological importance in tropical biomass. Therefore, studying the relationship between fern leaf area and leaf biomass at different elevations can further inform our understanding of the trade-off between photosynthetic gain versus biomass investment.
Our study shows that the scaling exponents of leaf area versus leaf biomass are consistently less than 1, this result is consistent with previous studies that there is a significant allometric relationship between the lamina area and mass (Sun et al., 2006; Huang et al., 2019; Jiao et al., 2022). However, our results differ at different elevations at both the individual leaf level and at the species level (Table 4). These results support the hypothesis called diminishing returns which predicts that increases in leaf area fail to keep pace with increases in dry biomass investments (Niklas et al., 2007). Additionally, the results indicate that differences in elevation can alter the scaling relationship between individual leaf area (ILA) and biomass (ILM).
5.1 Leaf area vs. leaf biomass relationships at the individual level at different elevations
At the whole plant level, the scaling exponent of the leaf area vs. leaf biomass relationship varied numerically at the three different elevations in this study, indicating that ferns have different strategies for leaf biomass investment at different elevations. This finding is consistent with previous research by Sun et al. (2017) on bamboo, which examined the relationship between leaf area and leaf biomass of bamboo at two different elevations in the Wuyi Mountain, southeast China, and found that the scaling exponent of both elevations was less than 1.0 but numerically different. Additionally, our findings are consistent with those of Geng et al. (2019), who discovered that the leaf area and leaf dry biomass scaling relationship of long-stalked bilberry differed with elevation although they confirmed the diminishing returns hypothesis.
However, prior investigation have reported data that contradict the concept of diminishing returns. For example, Pan et al. (2013) analyzed the relationship between leaf area and leaf biomass of 121 vascular plant species and found that the scaling exponents undergo numerical changes from values less than one to greater than one with increasing elevation. These divergent findings suggest that the scaling relationship between leaf area and leaf biomass is influenced by external environment condition and species composition.
5.2 Relationship between leaf area vs. leaf biomass at the species level at different elevations
At the species level, the relationship between ILA and ILM is isometric at the low elevation and follows an allometric scaling relationship at the middle and high elevations. We speculate that ferns may invest more biomass at higher elevation to increase the support of the lamina, which in turn may improve the light capture capacity and photosynthetic efficiency of leaves (Ackerly et al., 2002; Ishida et al., 2005; Niinemets et al., 2007). The plasticity of plant morphology in response to environmental variations is an important feature of plant that allows plants to allocate water and nutrients and other resources to the most needed tissues and organs by changing its structure and function, thus maximizing plant competitiveness (Li et al., 2013). It is important to bear in mind that changes in the numerical value of scaling exponents reflect adjustments in the proportional allocations in response to external environmental factors (Shi et al., 2015). The Changes in elevation gradient in mountain areas are influenced by many factors such as changes in light conditions, air pressure, temperature, and humidity, and soils (Chen and Wang, 2021). In this study, Wuyi National Park were characterized by aquatic hydrophytic vegetation at the low elevation, a predominantly bamboo forests at the middle elevation, and tall trees forested community are densely distributed at the high elevation. In the context of this study, it is noteworthy that previous studies suggest that the leaf dry mass per unit leaf (LMA) of plants increased with the increasing light availability (Ackerly et al., 2002). The changes in the numerical values of the leaf area vs. dry mass scaling relationship appear to reflect a greater investment in leaf area at the lowest and highest elevations, and a decreased investment at the mid-elevation. At the lowest and highest elevations, light may be the main primary limiting factor affecting plant access to resources. Therefore, the photosynthetic capacity per unit leaf area decreases under these conditions, making it important to increase the size of leaf area lamina. This adaptive response is significant as larger leaves tend to get more light resources, while small leaves can only intercept light through other parts (petiole, leaf angle) and features (leaf shape) an adaptive response (Ninemets et al., 2004). At the mid-elevation, on the contrary as the elevation increases, the temperature gradually decreases, strong radiation and winds have an impact on the productive performance of plants. In this case the plant leaves then need to allocate more biomass to support structures, the size of the leaf area directly affects the area under stress, so there is a significant difference in the proportion of biomass obtained by large and small leaves, while small leaves have a smaller LMA and thus a greater input to the support tissue.
In this study, leaf area, leaf biomass, and LMA were significantly higher at the mid-altitude than at low and high altitudes. Leaf biomass distribution depended on leaf size that larger leaves invest greater proportion of petiole biomass than smaller leaves, which collectively indicates a proportional reduction in lamina area compared to the investment in leaf dry biomass. The LMA increases in resource-poor, light-intense, high-altitude or arid habitats (Hanley et al., 2007). However, the increase in LMA caused by different habitats are extremely different (Wright et al., 2002). Because LMA is a composite parameter, the effect of light intensity on LMA is a synergistic effect of leaf density and thickness (Onoda et al., 2008). LMA increases with increasing light radiation and this phenotypic plasticity show differences between species (Niinemets et al., 1998; Valladares et al., 2000). This is consistent with Liu et al. (2010), which they showed that LMA were significant differences between the two light levels even for leaves of the same individual plant. While the maximum photosynthetic rate of Symplocos setchuensis differed significantly between the two light levels, there was no significant difference in LMA between the two, suggesting that species specificity also affects the expression of plasticity in LMA. This variation is ecologically significant, as plants in low light habitats capture light energy better by increasing leaf spread per unit mass, while in high light habitats they make better use of light energy by increasing dry matter content per unit area (Poorter et al., 2009). We speculate that the foliage leaf area of bamboos distributed at middle altitude is smaller and less shading to light, which is more favorable to the growth of ferns. Conversely, tall trees with dense vegetation and large canopy cover are mainly distributed at high altitude. Due to the shading of trees, the fronds of ferns are exposed to less sunlight, leading to a reduction in the LMA. This finding is consistent with the results of Zhu et al. (2011) on the functional traits of fern fronds under 16 different light habitat conditions. Ferns confined to two habitat conditions, the natural forest shade understory and disturbed open habitats adopted different cost effective strategies for fern fronds due to differences in light conditions. This phenomenonlogy reflects a response to an increase in the availability of light owing to the low foliage density of the bamboo community typical at the mid-elevation of the study site.
6 Conclusion
The results of this study show that, at the individual plant and species level, the scaling relationship for fern leaf area vs. dry biomass follows the diminishing returns hypothesis, but that the numerical values of the scaling exponents differ among three elevations in the study region. We speculate that the differences in the numerical values of the scaling exponents reflect differences in the ambient light variability. Collectively, our data indicate that the leaf area vs. biomass scaling relationship of ferns is consistent with prior studies focusing on seed plant species (angiosperms and gymnosperms) showing that leaf biomass investments are responsive to environmental factors, such as the availability of light.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
SC, DC, QZ and DH conceived the ideas and designed methodology; SC, JL, and JS collected the data. SC and JL analyzed the data; SC, JS, DH and DC led the writing of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (32071555, 32001294, 31971643) and the Key Public Welfare Project of Fujian Provincial Department of Science and Technology (2022R1002002).
Acknowledgments
The authors thank L. Cheng and Y. Fang, for their hospitality during the implementation of our study in the National Nature Reserve of Wuyi Mountain. We are very grateful to Karl J. Niklas for his many constructive suggestions that have helped to improve this paper.
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.1187704/full#supplementary-material
References
Ackerly, D. D., Knight, C. A., Weiss, S. B., Starmer, K. P. (2002). Leaf size, specific leaf area and microhabitat distribution of chaparral plants: contrasting patterns in species level and community level analyses. Oecologia 130, 449–457. doi: 10.1007/s004420100805
Blomberg, S. P., Garland, T., Ives, A. R. (2003). Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57, 717–745. doi: 10.1111/j.0014-3820.2003.tb00285.x
Chen, J. F. (2000). Genetic characteristics and taxonomic classification of soils in the mountain wuyi. Chin. J. Soil Sci. 03, 97–101+145. doi: 10.19336/j.cnki.trtb.2000.03.001
Chen, J. R., Wang, X. L. (2021). Altitude differences of allometric growth elationship of Pinus densata in southeast Tibet. HuNan Agric. Sci. 2, 29–32. doi: 10.16498/j.cnki.hnnykx.2021.002.008
Geng, M. Y., Chen, F. Q., Lü, K., Wang, Y. B., Guan, S. P., Liu, C. B. (2019). Changes of allometric relationships among leaf traits of disanthus cercidifolius var. longipes in different ontogenetic stages and altitude gradients. Guihaia 39, 1387–1397. doi: 10.11931/guihaia.gxzw201810038
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. doi: 10.3389/fpls.2022.832300
Hanley, M. E., Lamont, B. B., Fairbanks, M. M., Rafferty, C. M. (2007). Plant structural traits and their role in anti-herbivore defence. Perspective Plant Ecol. Evol. System 8, 157–178. doi: 10.1016/j.ppees.2007.01.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
Huang, W., Reddy, G. V. P., Li, Y., Larsen, J. B., Shi, P. (2020). Increase in absolute leaf water content tends to keep pace with that of leaf dry mass–evidence from bamboo plants. Symmetry 12, 1345. doi: 10.3390/sym12081345
Ishida, A., Yazaki, K. ,. A. L., Hoe, A. L. (2005). Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, macaranga gigantean. Tree Physiol. 25, 513–522. doi: 10.1093/treephys/25.5.513
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. doi: 10.3389/fpls.2022.860206
Jin, Y., Qian, H. (2022). V.PhyloMaker2:An updated and enlarged r package that can generate very large phylogenies for vascular plants. Plant Diversity 04, 335–339. doi: 10.1016/j.pld.2022.05.005
Jin, D. M., Zhou, X. L., Schneider, H., Wei, H. J., Wei, H. Y., Yan, Y. H. (2021). Functional traits: adaption of ferns in forest. J. Systematics Evolution 59, 1040–1050. doi: 10.1111/jse.12669
Kembel, S. W., Cowan, P. D., Helmus, M. R., Cornwell, W. K., Morlon, H., Ackerly, D. D., et al. (2010). Picante: r tools for integrating phylogenies and ecology. Bioinformat-ics 26, 1463–1464. doi: 10.1093/bioinformatics/btq166
Kleiman, D., Aarssen, L. W. (2007). The leaf size/number trade-off in trees. J. Ecol. 95, 376–382. doi: 10.1111/j.1365-2745.2006.01205.x
Korner, C. (2001). Alpine Plant Life: Functional Plant Ecology of High Mountain Ecosystems. 2nd ed. New York: Springer Berlin.
Li, Y., Zhao, C. Z., Dong, X. G., Hou, Z. J., Ma, X. L., Zhang, Q. (2013). Twig and leaf trait differences in Stellera chamaejasme with slope in alpine grassland. Chin. J. Plant Ecol. 37, 709–717. doi: 10.3724/SP.J.1258.2013.00074
Li, B. G., Zhu, H., Wang, H. (2000). The changing of fern diversity from the fragment rain forest on a holy hill in sipsongpanna. Wuhan Botanical Res. 06, 479–486. doi: 10.3969/j.issn.2095-0837.2000.06.008
Liu, M., Niklas, K. J., Niinemets, Ü, Hölscher, D., Chen, L., Shi, P. (2020). Comparison of the scaling relationships of leaf biomass versus leaf surface area between spring and summer for two deciduous tree species. Forests 11, 1010. doi: 10.3390/f11091010
Liu, F. J., Xiang, S., Yang, X. C., Sun, S. C. (2010). Comparison of leaf mass per area,photosynthetic capacity and chemical defense traits of four evergreen broad-leaved tree species under contrasting light conditions. Chin. J. Appl. Environ. Biol. 16, 462–467. doi: 10.3724/sp.j.1145.2010.00462
Mao, W., Li, Y. L., Zhang, T. H., Zhao, X. Y., Hang, X. X., Song, L. L. (2012). Research advances of plant leaf traits at different ecology scales. J. Desert Res. 32, 33–41.
Niinemets, Ü., Kull, O., Tenhunen, J. D. (1998). An analysis of light effects on foliar morphology, physiology, and light interception in temperate deciduous woody species of contrasting shade tolerance. Tree Physiol. 18, 681–696. doi: 10.1093/treephys/18.10.681
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. Botany 100, 283–303. doi: 10.1093/aob/mcm107
Niklas, K. J., Cobb, E. D., Ninemets, Ü., 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. Sciences USA 104, 8891–8896. doi: 10.1073/pnas.0701135104
Ninemets, Ü., Najwa, A. A., 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
Onoda, Y., Schieving, F., Anten, N. P. R. (2008). Effects of light and nutrientavailability on leaf mechanical properties of plantago major: a conceptual approach. Ann. Botany 101 (5), 727–736. doi: 10.1093/aob/mcn013
Page, C. N. (2002). Ecological strategies in fern evolution: a neopteridological overview. Rev. Palaeobotany Palynol. 119, 1–33. doi: 10.1016/S0034-6667(01)00127-0
Pan, S., Liu, C., Zhang, W. P., Xu, S. 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
Pearcy, R. W., Muraoka, H. F., Valladares, F. (2005). Crown architecture in sun and shade environments: assessing function and trade-offs with a three dimensional model. New Phytologist 166, 79–800. doi: 10.1111/j.1469-8137.2005.01328.x
Poorter, H., Niinemets, Ü., Poorer, L., Wright, I. J., Villar, R. (2009). Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis. New Phytologist 182, 565–588. doi: 10.1111/j.1469-8137.2009.02830.x
Schneider, H., Schuettpelz, E., Pryer, K. M., Cranfill, R., Magallón, S., Lupia, R. (2004). Ferns diversified in the shadow of angiosperms. Nature 428, 553–557. doi: 10.1038/nature02361
Shi, P., Li, Y., Hui, C., Ratkowsky, D. A., Yu, X., Niinemets, Ü. (2019). Does the law of diminishing returns in leaf scaling apply to vines? evidence from 12 species of climbing plants. Global Ecol. Conserv. 21, e00830. doi: 10.1016/j.gecco.2019.e00830
Shi, Y. C., Zhao, C. Z., Song, Q. H., Du, J., Wang, J. W. (2015). Allometric relationship between height and crown width or diameter of platycladus orientalis on different slope aspects of lanzhou northern mountains. Chin. J. Ecol. 34, 1879–1885. doi: 10.13292/j.1000-4890.20150615.009
Shipley, B. (1995). Structured interspecific determinants of specific leaf area in 34 species of herbaceous angiosperms. Funct. Ecol. 9, 312–319. doi: 10.2307/2390579
Shipley, B., Lechowicz, M. J., Reich, W. P. B. (2006). Fundamental trade-offs generating the worldwide leaf economics spectrum. Ecology 87, 535–541. doi: 10.1890/05-1051
Sun, J., Fan, R. R., Niklas, K. J., Zhong, Q. L., Yang, F., Li, M., et al. (2017). “Diminishing returns” in the scaling of leaf area vs. dry mass in wuyi mountain bamboos, southeast China. Am. J. Botany 104, 993–998. doi: 10.3732/ajb.1700068
Sun, S., Jin, D., Shi, P. (2006). The leaf size-twig size spectrum of temperate woody species along an altitudinal gradient: an invariant allometric scaling relationship. Ann. Botany 97, 97–107. doi: 10.1093/aob/mcj004
Valladares, F., Wright, S. J., Lasso, E., Kitajima, K., Pearcy, R. W. (2000). Plastic phenotypic response to light of 16 congeneric shrubs from a Panamanian rainforest. Ecology 81, 1925–1936. doi: 10.1890/0012-9658(2000)081[1925:PPRTLO]2.0.CO;2
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. Systematics 33, 125–159. doi: 10.1146/annurev.ecolsys.33.010802.150452
Wright, I. J., Westoby, M., Reich, P. B. (2002). Covergence towards higher leaf mass per area in dry and nutrient-poor habitats has different consequenses for leaf life span. J. Ecol. 90, 534–543. doi: 10.1046/j.1365-2745.2002.00689.x
Yan, Y. H., Yi, Q. F., Huang, Z. L., Xing, F. W. (2004). The ecological response of fern diversity to vegetation succession in gudoushan nature reserve, guangdong. Biodivs. Sci. 12, 339–347. doi: 10.3321/j.issn:1005-0094.2004.03.006
Yang, D. M., Zhang, J. J., Zhou, D., Qian, M. J., Zheng, Y., Jin, L. M., et al (2012). Leaf and twig functional traits of woody plants and their relationships with environmental change: A review. Chin. J. Ecol. 31, 702–713. doi: CNKI:SUN:STXZ.0.2012-03-032
Zhang, M., Zhang, S., Ye, M., Jiang, L., Vallejos, C. E., Wu, R. (2020). The genetic control of leaf allometry in the common bean, phaseolus vulgaris. BMC Genet. 21, 1. doi: 10.1186/s12863-020-00838-2
Zhu, J. D., Meng, T. T., Ni, J., Su, H. X., Xie, Z. Q., Zhang, S. R., et al. (2011). Within-leaf allometric relationships of mature forests in different bioclimatic zones vary with plant functional types. Chin. J. Plant Ecol. 35, 687–698. doi: 10.3724/SP.J.1258.2011.00687
Zhu, S. D., Li, R. H., Song, J., He, P. C., Liu, H., Berninger, F., et al. (2015). Different leaf cost-benefit strategies of ferns distributed in contrasting light habitats of sub-tropical forests. Ann. Botany 117, 497–506. doi: 10.1093/aob/mcv179
Keywords: leaf area-biomass relationship, diminishing returns hypothesis, elevation, subtropical ecosystem, ferns
Citation: Chen S, Li J, Sun J, Zhong Q, Hu D and Cheng D (2023) “Diminishing returns” and leaf area-biomass scaling of ferns in subtropical ecosystems. Front. Plant Sci. 14:1187704. doi: 10.3389/fpls.2023.1187704
Received: 24 March 2023; Accepted: 06 June 2023;
Published: 27 June 2023.
Edited by:
Alejandro Ordonez, Aarhus University, DenmarkReviewed by:
Peijian Shi, Nanjing Forestry University, ChinaDong-Mei Jin, Shanghai Chenshan Plant Science Research Center (CAS), China
Copyright © 2023 Chen, Li, Sun, Zhong, Hu and Cheng. 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: Dandan Hu, hudanworld@126.com; Dongliang Cheng, chengdl02@aliyun.com