The scaling relationship between leaf nitrogen and phosphorus concentrations in vascular epiphytes

Introduction: The scaling relationship between leaf nitrogen (N) and phosphorus (P) concentrations reflects plant adaptation strategies and evolutionary dynamics. While extensively studied in terrestrial plants, vascular epiphytes—a key yet understudied component of global biodiversity—remain poorly understood.

Methods: We compiled leaf N and P data from 38 epiphyte species across tropical seasonal rainforests and subtropical montane forests in China, supplemented by a global literature synthesis. Standardized major axis (SMA) regression analyzed N-P scaling exponents (β) across forest types, functional groups, and habitats.

Results: Epiphytes exhibited a distinct global N-P scaling exponent (β=0.78), significantly higher than terrestrial plants. Facultative epiphytes showed lower β in epiphytic (0.69) versus terrestrial habitats (0.91). No significant variation occurred among functional groups or forest types, suggesting conserved nutrient allocation strategies.

Discussion: The elevated β underscores epiphytes’ reliance on atmospheric nutrient inputs and adaptive P retention. Habitat-driven differences highlight niche specialization, while functional group uniformity reflects stoichiometric constraints of canopy living. These findings redefine epiphyte nutrient ecology within broader plant stoichiometry theory.

Introduction

The scaling relationship between leaf nitrogen (N) and phosphorus (P) concentrations provides an essential framework for understanding plant nutrient allocation and physiological adaptation. Because N and P are key components of proteins and rRNA, their stoichiometric balance reflects fundamental constraints on plant growth and metabolism (Niklas et al., 2005; Rivas-Ubach et al., 2012). This relationship is often described by the power function [N] = α[P]^β, where β represents the scaling exponent (Niklas et al., 2005). Unlike a static N:P ratio, which has been widely used to indicate nutrient limitation but lacks sensitivity to variation across scales (Güsewell, 2004), β quantifies the proportional accumulation of N relative to P across species and environments. Variation in β reveals differences in nutrient investment and adaptive strategies among plant groups (Tian et al., 2018; Guo et al., 2020; Ouyang et al., 2024; Rao et al., 2025), thereby offering a mechanistic perspective on nutrient coupling and stoichiometric evolution in plant.

Generally, β is less than 1, indicating that plants require greater P-rich rRNA to sustain rapid protein synthesis (Sterner and Elser, 2002; Elser et al., 2010). Reich et al. (2010) proposed a general β ≈ 2/3 for terrestrial leaves, yet subsequent studies have revealed substantial variability, driven by both biotic and abiotic factors. Among biological drivers, plant functional types and growth strategies exert a major influence: fast-growing or deciduous species often exhibit lower β values due to their higher P demand and shorter nutrient turnover time, whereas evergreen or stress-tolerant species display higher β values reflecting P-use conservation (Tian et al., 2018; Guo et al., 2020; Rao et al., 2025). Environmental factors such as temperature, precipitation, light intensity, and soil nutrient status also play decisive roles. Warm, moist, and nutrient-rich environments typically favor lower β values, consistent with accelerated P uptake and high growth rates (Guo et al., 2020; 2021), whereas nutrient-poor or P-limited soils promote higher β values, indicating enhanced P investment efficiency and resistance to P limitation (Guo et al., 2021). Additionally, soil N–P allometry can constrain leaf nutrient coupling, implying that leaf-level β partly mirrors the stoichiometric structure of its edaphic environment (Guo et al., 2021). Together, these findings suggest that β is shaped by the interplay of biological traits, resource availability, and environmental filtering. As such, β serves not only as an indicator of nutrient allocation flexibility but also as a functional trait reflecting plant adaptation and evolutionary responses to ecological niches (Tian et al., 2018; Guo et al., 2020; Rao et al., 2025).

Most studies on leaf N-P scaling relationships have focused on terrestrial plants (Reich et al., 2010; Tian et al., 2018; Guo et al., 2021), while vascular epiphytes—plants that depend on other plants for structural support and account for approximately 10% of global plant biodiversity (Zotz et al., 2021b) —remain underexplored in this context. Unlike terrestrial plants, epiphytes lack a root system connected to the soil and do not directly extract nutrients from the soil (Zotz and Hietz, 2001; Zotz, 2016). As a result, their nutrient uptake and allocation strategies may differ significantly. Several studies have highlighted differences in leaf nutrient concentrations between epiphytes and terrestrial plants (Cardelús and Mack, 2010; Wu et al., 2020; Hietz et al., 2022; Sulaiman et al., 2023). For example, a recent study found that LNC and N:P ratios are significantly lower in epiphytes compared to terrestrial plants, while LPC does not exhibit such difference (Hietz et al., 2022). Additionally, epiphytes are expected to have lower relative growth rates than terrestrial plants (Zotz, 2016; Hietz et al., 2022). These observations suggest that β of epiphytes may differ markedly from that of terrestrial plants.

Variation in β has been documented for terrestrial plants (Tian et al., 2018; Guo et al., 2020, 2021), with its dependence on species functional groups, ecoregions and sampling sites (Tian et al., 2018). For epiphytes, significant differences in nutrient status have been observed between epiphytic ferns and seed plants. Specifically, epiphytic ferns generally exhibit higher LNC and LPC compared to seed plants (Cardelús and Mack, 2010; Hietz et al., 2022), potentially leading to differences in β between these two groups. Beyond nutrient status and canopy microclimate, additional physiological mechanisms may also influence the N–P scaling exponent (β). For example, variations in photosynthetic pathways (e.g., C₃, C₄, and CAM (Crassulacean Acid Metabolism) metabolism) can alter N allocation to photosynthetic enzymes and phosphorus investment in energy metabolism, thereby modulating the proportional accumulation of N and P (Elser et al., 2010; de Paula Oliveira et al., 2021). Similarly, hydraulic traits such as xylem conductivity, leaf water potential, and water-use efficiency affect nutrient transport and metabolic rates under different water availabilities, potentially influencing the N–P coupling relationship (Yuan and Chen, 2015; Li et al., 2025). In epiphytes, where water and nutrient acquisition largely depend on atmospheric inputs (Hietz et al., 2002), these physiological traits may play an even greater role in shaping β by mediating trade-offs between photosynthetic efficiency, nutrient conservation, and drought tolerance (Querejeta et al., 2018; Gotsch et al., 2022). Vascular epiphytes communities also experience high species turnover at small spatial scales (Zhao et al., 2015; Amici et al., 2020), and the canopy microclimate plays a significantly role in modulating the physiological processes of individual epiphytes (Petter et al., 2016). These ecological and environmental variations suggest that β values in epiphytes may differ across groups, sites and habitats.

In this study, we aim to address these gaps in knowledge by presenting a comprehensive analysis of leaf N-P scaling relationships in vascular epiphytes. Our objectives are twofold: First, we compare the β values of epiphytes to those of terrestrial plants, testing whether the observed differences in nutrient status and growth rates translate into significant variation in β. Second, we examine β variation across different functional groups and spatial scales to determine if it correlates with functional groups, forest types, and habitats. By integrating global data on leaf N and P concentrations with detailed local measurements from epiphytes in a subtropical moist forest and a tropical seasonal rainforest, this study offers a novel perspective on leaf N-P scaling relationships in epiphytes and contributes to understanding their adaptive strategies in diverse ecological contexts.

Methods

Study site

This study was conducted in the subtropical montane moist forest of Xujiaba, located in the core area of the Ailao Mountains National Natural Reserve (101°01′E, 24°32′N, 2000~2650m altitude) in Yunnan, Southwest (SW) China. The region has a subtropical monsoon climate influenced by both southwest and southeast monsoons (Liu et al., 2003). Meteorological data (2011-2019) from the subtropical Forest Ecosystem Research Station at Xujiaba (2450m) indicate an annual mean precipitation of 1730 mm, with 85% occurring between May and October, and an annual mean temperature of 11.3°C (Gnanamoorthy et al., 2021). Primary forests dominate the area (77.9% of 667km²) is dominated by primary forests, primarily composed of subtropical evergreen broadleaved species such as Lithocarpus hancei, Castanopsis wattii, and Lithocarpus xylocarpus (Li et al., 2013; Wen et al., 2018). These primary forests (trees>3.5cm DBH) have a basal area of 77.31 m²·ha^-1, canopy openness of 5.8%, mean DBH of 19.23 cm, and tree density of 1656 trees·ha^-1 (Li et al., 2013).Secondary forests (~16% of the total area) have regenerated for 50–100 years after disturbance and show smaller basal area (55.17 m²·ha^-¹), lower mean DBH (9.84 cm), higher canopy openness (31%), and greater tree density (5903.33 trees·ha^-1; Li et al., 2013). A 28.5 km² reservoir adjacent to primary forests creates a humid microclimate that supports high epiphyte diversity. Facultative epiphytes capable of both terrestrial and epiphytic growth are common and provide an ideal system for N-P scaling relationships studies (Zotz, 2016; Zhang et al., 2020).

For comparison, samples were also collected from a tropical seasonal rainforest in Bubeng, within the Mengla National Natural Reserve (101°35′E, 21°37′N) in Yunnan, SW China. This region has a typical monsoon climate with a six-month dry season (November-April) and six-month rainy season (May-October), an annual mean precipitation of 1493 mm (84% in the rainy season), and an annual mean temperature of 21.8°C (Cao et al., 2006). The forest is dominated by tree species such as Parashorea chinensis, Pittosporopsis kerrii, Garcinia cowa and Castanopsis echidnocarpa.

Dataset description

To examine variation in N and P concentrations and their scaling relationships, two datasets were compiled (Appendices A and B). The individual-level dataset (Appendix A) includes 1,803 data points for 274 species, combining: (1) 1390 individuals from 236 species (de Paula Oliveira et al., 2021; Hietz et al., 2022), and (2) original measurements from 343 individuals from 20 vascular epiphytes in Ailao Mountains and 70 individuals from 18 vascular epiphytes in Bubeng.

The species-level dataset (Appendix B) was obtained by averaging per species to reduce sampling bias and incorporatingdata from Zotz (2004) and Querejeta et al. (2018), comprising 323 species across 41 families (103 ferns, 223 seed plants).

Vascular epiphyte sampling

In August 2018, vascular epiphytes were collected from primary forests (both interior and reservoir-adjacent) and secondary forests in Xujiaba using ladders, twig shears, and free-hand climbing. One or more conspcific individuals of small-sized epiphytes collected from the same host tree were treated as a single sample.

A total of 343 individuals from 20 species were collected (Supplementary Table S1): 107 individuals from reservoir-adjacent forests, 161 individuals from primary forests, and 75 individuals from secondary forests. Host tree leaves (51 individuals, 16 species) were also sampled, focusing on outer canopy layers. Epiphytic and terrestrial individuals of five facultative epiphyte species (Briggsia longifolia, Cautleya gracilis, Chirita macrophylla, Elatostema monandrum, and Impatiens siculifer) were collected from nearby trees and ground sites. All samples were washed to remove dust and mud, oven-dried at 80 °C for 48 hours, and stored in envelopes.

In July 2018, vascular epiphyte samples of the tropical seasonal rainforest were collected in Bubeng. Five 20×20 m plots were established at the lowest (709 m) and highest (869 m) altitude points of the 20-ha tropical forest dynamics plot (101°34′26″–47″E, 21°36′42″–58″N), as well as at the periphery of the 1.44-ha square plot (101°34′59″E, 21°37′2″N), with plots separated by more than 20 m. A total of 70 epiphyte individuals from 18 species and 126 tree individuals from 18 species were collected. The same processing protocols were applied.

Chemical analysis

For each sample, leaf nitrogen concentration (LNC, in g·kg^-¹) and leaf phosphorus concentration (LPC, in g·kg^-¹) were measured. Samples were ground manually with a mortar to ensure complete homogenization of small leaves and then dried at 65 °C for 24 hours prior to chemical analysis. The leaf chemistry analyses were performed at the Biogeochemical Laboratory of the Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences. LNC was determined using an elemental analyzer (Vario MAX CN, Elemental Analysensysteme GmbH, Germany). For LPC determination, 0.2–0.3 g of each sample was weighed and digested with 10 mL of HNO₃ and 2 mL of HClO₄. The digestion process was carried out at room temperature for 5 hours, followed by heating at 185 °C for 2 hours. After digestion, the samples were diluted to a final volume of 50 mL with deionized water. Mass-based P concentrations were quantified using an inductively coupled plasma-atomic emission spectrometer (iCAP6300, Thermo Fisher Scientific, U.S.A.).

Data analysis

All leaf stoichiometry data were log-transformed prior to analysis due to their strongly left-skewed distributions. The N-P scaling relationships for vascular epiphytes, as well as for epiphytic ferns and seed plants across different scales, were quantified using standardized major axis (SMA) regression (Wright et al., 2005). This was implemented with the ‘sma’ function in the ‘smart’ package (Warton et al., 2012) in R. Likelihood ratio tests were employed to compare β among vascular epiphytes.

Results

Global vascular epiphytes exhibited significant correlations between LNC and LPC at both individual and species levels (Figure 1). At the individual level, the leaf N-P scaling exponent (β) for global epiphytes was 0.78 (95% CI: 0.75–0.81), with R² = 0.62 and p < 0.001 (Figure 1A). The likelihood test indicated that the exponent was significantly different from the theoretical value of 2/3 (H₀: β = 0.667; r = 0.1931, df = 1801, p < 0.001). At the species level, the β was 0.69 (95% CI: 0.65–0.73), R² = 0.55, p < 0.001 (Figure 1B).

Figure 1

Figure 1. Leaf N-P scaling relationships of vascular epiphytes at global scale. (A) Individual level, (B) Species level. Scaling exponents (β) were calculated from the SMA regression between leaf N and P concentrations. The data in parentheses after the β values in the table represent the 95% confidence interval.

In comparisons with trees, epiphytes in the tropical seasonal rainforest displayed significantly lower LNC (13.43 g kg^-¹) and N:P ratio (11.05) than trees (LNC = 19.15 g kg^-¹; N:P = 15.69), while those in the subtropical montane moist forest showed significantly higher LNC (21.47 g kg^-¹) and leaf phosphorus concentration (LPC = 2.07 g kg^-¹), and a lower N:P ratio (10.92) compared with trees (LNC = 16.88 g kg^-¹; LPC = 1.14 g kg^-¹; N:P = 15.17) (Supplementary Table S1). In the subtropical montane moist forest, β for vascular epiphytes was 0.84 (95% CI: 0.78–0.90; R² = 0.57; p < 0.001), compared to 1.01 (95% CI: 0.95–1.07; R² = 0.64; p < 0.001) for trees (Figure 2A). In the tropical seasonal rainforest, β values were 0.67 (95% CI: 0.61–0.74; R² = 0.52; p < 0.001) for epiphytes and 0.60 (95% CI: 0.54–0.66; R² = 0.48; p < 0.001) for trees (Figure 2B). Likelihood ratio tests confirmed that these β differences between epiphytes and trees were not statistically significant (p = 0.25). Facultative epiphytes exhibited significantly lower β values in epiphytic habitats (0.69 ± 0.04, R² = 0.49, p < 0.001) compared to terrestrial habitats (0.91 ± 0.05, R² = 0.63, p < 0.001; Figure 3B). This difference was significant (p = 0.012), suggesting the role of epiphytic habitats in shaping N–P scaling.

Figure 2

Figure 2. Leaf N-P scaling relationships of vascular epiphytes and trees in a subtropical montane forest (a) and a tropical seasonal rainforest (b). The leaf N-P scaling exponents (β) were calculated from the SMA regression between leaf N and P concentrations. The data in parentheses after the β values in the table represent the 95% confidence interval. The likelihood ratio test comparing leaf N–P scaling exponents (β) among forest types for epiphytes is p = 0.315, and for trees is p = 0.061.

Figure 3

Figure 3. Leaf N-P scaling relationships of vascular epiphytes under different habitats at individual level (A), and the epiphytic and terrestrial individuals of facultative epiphytes (B). PFE, primary forest edge near the reservoir; PF, primary forest; SF, secondary forest. The N-P scaling exponents (β) were calculated from the SMA regression between leaf N and leaf P concentrations. The data in parentheses after the β values in the figure represent the 95% confidence interval.

At the global scale, the mean LNC and LPC were 16.09 mg g^-1 and 1.49 mg g^-1 in epiphytic ferns, respectively, compared to 13.81 mg g^-1 and 1.47 mg g^-1 in seed plants (Supplementary Table S1). The corresponding N:P ratios were 13.10 and 11.75 for the two groups, respectively (Supplementary Table S1). Surprisingly, no significant differences in β were observed between global epiphytic ferns and seed plants (Figure 4A). Similarly, in the subtropical montane moist forest, β did not differ significantly between epiphytic ferns and seed plants (Figure 4B), despite the higher nutrient concentrations observed in ferns.

Figure 4

Figure 4. Leaf N-P scaling relationships of epiphytic ferns and seed plants at different scales. (A) Global scale, (B) Local scale. Scaling exponents (β) were calculated from the SMA regression between leaf N and P concentrations. The data in parentheses after the β values in the table represent the 95% confidence interval.

Furthermore, β showed no significant variation across forest types (p = 0.315, Figure 2) or local habitats (p = 0.421, Figure 3A). In the subtropical montane moist forest, β were 0.72 in primary forests adjacent to the reservoir, 0.80 in primary forests, and 0.76 in secondary forests (Figure 3A). Notably, the strength of the leaf N-P correlation declined from 0.62 (R² of the SMA fit) in primary forest forests adjacent to the reservoir to 0.34 in primary forests and to 0.22 in secondary forests (Figure 3A), indicating a tighter coupling of LNC and LPC in habitats with more favorable growing conditions.

Discussion

Uncovering patterns in leaf N-P scaling relationships is crucial for advancing our understanding of epiphyte adaptation and evolution. Our finding reveal distinct and consistent patterns in these relationships across spatial scales, underscoring the strong influence of the epiphytic habitat on the relative accumulation of N and P in leaves (Figure 5). This study provides new mechanistic insight by demonstrating that the epiphytic lifestyle fundamentally alters nutrient scaling relationships compared with terrestrial plants.

Figure 5

Figure 5. Leaf N and P concentrations and their N-P scaling exponents (β) of vascular epiphytes from global to local scales. The Leaf N-P scaling exponents (β) were calculated from the SMA regression between LNC and LPC. The global leaf β values of vascular epiphytes were calculated based on the leaf N and P contents from the global epiphyte functional traits dataset (Hietz et al., 2022), other relevant studies (Zotz, 2004; Querejeta et al., 2018; de Paula Oliveira et al., 2021) and this study. The leaf β values of vascular epiphytes at the regional and local scales were calculated based on the leaf N and P concentrations of vascular epiphytes in the subtropical montane moist forest in Ailao Mountain and the tropical seasonal rainforest in Bubeng in this study. ALM, the subtropical montane moist forest in Ailao Mountains; BBV, the tropical seasonal rainforest in Bubeng. PF (primary forest), PFE (primary forest edge near the reservoir) and SF (secondary forest) are different habitats in ALM. LNC, leaf nitrogen concentrations; LPC, leaf phosphorus concentrations.

Do epiphytes and terrestrial plants differ in β?

Vascular epiphytes exhibit distinct β compared with terrestrial plants. Global analysis reveals a higher β (0.78) for epiphytes than for terrestrial plants (Reich et al., 2010; Tian et al., 2018), indicating a faster relative accumulation of N with respect to P. Mechanistically, in SMA estimation β is proportional to the ratio of the dispersion of log-transformed N to that of log-transformed P (β ≈ SD[log N]/SD[log P]); therefore, greater variability in leaf N combined with relatively stable LPC mathematically yields a higher β (Warton et al., 2006; Guo et al., 2020). In epiphytes, the pronounced sensitively of leaf N to fluctuating N supply under arid or nutrient-poor conditions (Cheaib et al., 2025), together with high spatial–temporal heterogeneity in canopy nutrient inputs (Zotz, 2016), increases variability in leaf N. Meanwhile, epiphytes often exhibit “luxury consumption” of P, maintaining stable LPC despite environmental fluctuations (Winkler and Zotz, 2009; Hadi Alkarawi and Zotz, 2014; Zotz, 2016). Together, these factors explain the higher β observed in epiphytes. Lower N:P further support this interpretation. Our pooled data show that global epiphytes have a lower mean N:P ratio (11.09) than terrestrial plants (15.8, Tian et al., 2018), consistent with dominant N limitation in epiphytic habitats (Güsewell, 2004; Hietz et al., 2022; Liu et al., 2025). The lower N availability in epiphytic habitats may enhance the relative responsiveness of leaf N to nutrient inputs, producing faster proportional N accumulation with respect to P and thereby contributing to the higher β observed in epiphytes. This interpretation aligns with Guo et al. (2020), who proposed that N limitation influences N–P scaling primarily by altering relative investment patterns rather than by directly constraining leaf N concentrations. Supporting this, facultative epiphytes exhibit a lower β in epiphytic habitats (0.69) than in terrestrial habitats (0.91), indicating that when N becomes more available and stable—as in terrestrial environments—absolute N accumulation increases, but relative N variability decreases. Thus, nutrient availability regulates both the magnitude and proportionality of N and P accumulation, reconciling the apparent contrast between low- and high-N conditions.

Beyond nutrient limitation, several physiological and ecological characteristics of epiphytes may shape their leaf N-P scaling relationships. Because these plants rely primarily on intermittent atmospheric inputs rather than soil-derived nutrients (Zotz, 2016), their nutrient acquisition is closely linked to water balance (Querejeta et al., 2018). Morphological adaptations such as trichomes, velamen, and succulent tissues enhance the interception and retention of throughfall and mist, thereby influencing both total nutrient availability and within-tissue allocation (Zotz, 2016). Many epiphytes exhibit CAM or intermediate photosynthetic pathways that enhance water-use efficiency but reduce nutrient turnover (Gilman and Edwards, 2020; Zotz et al., 2021a). This constraint slows P cycling and results in a relative accumulation of N in structural and enzymatic compounds (Schleuss et al., 2020). Such physiological adjustments may explain the elevated β values observed, reflecting N’s dominant role in maintaining metabolic capacity under prolonged drought or fluctuating nutrient supply. In additions, the pronounced microclimatic heterogeneity of epiphytic habitats—characterized by large diurnal fluctuations in temperature and vapor pressure deficit—likely reinforce conservative nutrient-use strategies (Yu et al., 2023). Epiphytes exposed to higher irradiance and desiccation may allocate more N to photoprotective pigments and osmotic-regulating metabolites rather than growth-related functions (Adibah and Ainuddin, 2011; He et al., 2024; Liu et al., 2025). Collectively, these adaptive features suggest that the leaf N–P scaling relationship in epiphytes emerge from an integration of nutrient limitation, water balance, and photosynthetic physiology.

Notably, our study found no significant difference in β between vascular epiphytes and trees in the subtropical montane moist forest or the tropical seasonal rainforest. This similarity may result from two main factors. First, the leaf N–P scaling exponent of plant is likely constrained by macroscale climatic conditions (Tian et al., 2018; Guo et al., 2021). Consequently, both epiphytes and trees within the same forest type may experience similar climatic constraints, leading to convergent β values. Second, most vascular epiphytes in our dataset are herbaceous or fern species, whereas the terrestrial plants are mainly woody trees. In terrestrial ecosystems, herbaceous plants generally exhibit lower β values than woody species (Tian et al., 2018); however, the relatively higher β of epiphytes offsets this difference, resulting in no significant contrast between the epiphytic (herbaceous) and terrestrial (woody) groups. Future research should quantitatively assess the influence of climatic factors on β in epiphytes and compare the β of epiphytic herbaceous species with that of terrestrial herbaceous plants to better understand the combined effects of plant functional type and environment on nutrient-scaling relationships.

Leaf N-P scaling exponents of epiphytes across different scales

Inconsistent with Tian et al. (2018), no significant differences in β values between epiphytic ferns and seed plants at both global and subtropical montane moist forest scales in our study. However, this supports the idea that β is uniform across major functional groups (Reich et al., 2010).

Additionally, no significant difference in β were observed among epiphytes across forest types. However, trees in the tropical seasonal rainforest exhibited a slightly lower β compared to those in the subtropical montane moist forest. This can be attributed to the tropical seasonal rainforest being more P-limited and warmer, which promotes faster P accumulation relative to N in trees, leading to a lower β (Tian et al., 2018a; Guo et al., 2020, 2021). In contrast, for epiphytes, the relative accumulation of N and P in leaves appears to be affected by N limitation in epiphytic habitats (Figures 1, 2, 4). Based on our analysis, epiphytes in these habitats tend to accumulate of N faster relative to P. Consequently, the differences in β between epiphytes in the two forest types are less pronounced than those observed in trees.

At the local scale, β values showed no significant variation across habitats. A previous study in the subtropical montane moist forest of the Ailao Mountains, Southwest China, found that epiphytes in primary forests exhibited higher annual rhizome growth than those in secondary forests (Chen et al., 2019), suggesting faster growth rates and greater P accumulation in primary forests (Guo et al., 2020). Yet, epiphytes in primary forests did not display a lower β compared to those in secondary forests. Combined with our findings, this suggests that the unique characteristics of epiphytic habitats—particularly their reliance on atmospheric nutrient inputs (e.g., through rain and canopy deposition) and detrital canopy soils, coupled with limited access to groundwater (Zotz, 2016)—constrain leaf N-P scaling relationship, favoring a trade-off in nutrient uptake that leans more towards N than P (Guo et al., 2020).

Interestingly, we observed a decreasing trend in the correlation between leaf N and P concentrations from primary forests adjacent to the reservoir to primary forests and secondary forests. This aligns with the view that changes in global water availability can decouple plant N:P (Yuan and Chen, 2015). To cope with limited water supply, a significant fraction of N in epiphytes is allocated to functions other than growth, such as water retention and stress tolerance (Sardans et al., 2012; Zotz, 2016). This allocation likely contributes to the weaker correlation between LNC and LPC in secondary forests.

While this study offers new insights into the mechanisms governing N–P scaling relationships in vascular epiphytes, several limitations should be noted. First, field sampling was confined to two forest types in southwest China, which may not capture spatial variation across broader biogeographic gradients. Second, the taxonomic coverage of epiphytes was uneven, with some functional groups and lineages underrepresented. Moreover, we did not directly measure microclimatic or physiological variables—such as canopy humidity, nutrient uptake rates, or photosynthetic efficiency—that could further elucidate the drivers of N–P coupling. Future work integrating trait-based experiments, and nutrient manipulation studies will help refine these conclusions and strengthen causal inference.

Conclusions

In this study, we investigated the leaf N-P scaling exponent of vascular epiphytes in comparison to terrestrial plants, examined its variation across different scales (Figure 5). Based on our findings, we conclude that (1) global vascular epiphytes exhibit distinct leaf N-P scaling exponents compared to terrestrial plants; (2) there are no significant differences in β across functional groups, forest types, or habitats, indicating the constraints imposed by the epiphytic lifestyle. Overall, this study uncovers key patterns in the leaf N-P scaling relationships of vascular epiphytes, and emphasizes the critical role of the epiphytic habitat in shaping these relationships, thereby advancing our understanding of epiphyte adaptation and evolution.

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 author.

Author contributions

TH: Conceptualization, Data curation, Funding acquisition, Writing – original draft. TTZ: Data curation, Writing – review & editing. DT: Data curation, Writing – review & editing. SLiu: Data curation, Writing – review & editing. SLi: Writing – review & editing. XH: Data curation, Writing – review & editing. YM: Writing – review & editing. WYL: Conceptualization, Funding acquisition, Writing – review & editing.

Funding

The author(s) declared financial support was received for this work and/or its publication. This study was supported by the National Natural Science Foundation of China (42071071, 31770496), the 14th Five-Year Plan of the Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences (E3ZKFF7B), and Jiujiang Natural Science Foundation (2025_001561, 2025_001334, 2025_001377).

Acknowledgments

We thank the members of the Public Technology Service Center of the XTBG for their assistance with the chemical analysis of plant samples. We would also like to thank the National Field Scientific Observation and Research Station of Subtropical Forest Ecosystem in Ailao Mountain and the National Field Scientific Observation and Research Station of Tropical Forest Ecosystem in Xishuangbanna for providing access to the study sites and background information.

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.

The handling editor MZ declared a shared parent affiliation with the authors TH, SLi, YM, WYL at the time of review.

The reviewer MX declared a shared parent affiliation with the authors TH, SLi, YM, WYL to the handling editor at the time of review.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.2025.1712082/full#supplementary-material

References

Adibah, M. S. and Ainuddin, A. N. (2011). Epiphytic plants responses to light and water stress. Asian J. Plant Sci. 10, 97–107. doi: 10.3923/ajps.2011.97.107

Crossref Full Text | Google Scholar

Amici, A. A., Nadkarni, N. M., Williams, C. B., and Gotsch, S. G. (2020). Differences in epiphyte biomass and community composition along landscape and within-crown spatial scales. Biotropica 52, 46–58. doi: 10.1111/btp.12725

Crossref Full Text | Google Scholar

Cao, M., Zou, X., Warren, M., and Zhu, H. (2006). Tropical forests of xishuangbanna, China. Biotropica 38, 306–309. doi: 10.1111/j.1744-7429.2006.00146.x

Crossref Full Text | Google Scholar

Cardelús, C. L. and Mack, M. C. (2010). The nutrient status of epiphytes and their host trees along an elevational gradient in Costa Rica. Plant Ecol. 207, 25–37. doi: 10.1007/s11258-009-9651-y

Crossref Full Text | Google Scholar

Cheaib, A., Waring, E. F., McNellis, R., Perkowski, E. A., Martina, J. P., Seabloom, E. W., et al. (2025). Soil nitrogen supply exerts largest influence on leaf nitrogen in environments with the greatest leaf nitrogen demand. Ecol. Lett. 28, e70015. doi: 10.1111/ele.70015

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Q., Lu, H.-Z., Liu, W.-Y., Wu, Y., Song, L., and Li, S. (2019). Obligate to facultative shift of two epiphytic Lepisorus species during subtropical forest degradation: Insights from functional traits. For. Ecol. Manage. 435, 66–76. doi: 10.1016/j.foreco.2018.12.037

Crossref Full Text | Google Scholar

de Paula Oliveira, R., Zotz, G., Wanek, W., and Franco, A. C. (2021). Leaf trait co-variation and trade-offs in gallery forest C ₃ and CAM epiphytes. Biotropica 53, 520–535. doi: 10.1111/btp.12895

Crossref Full Text | Google Scholar

Elser, J. J., Fagan, W. F., Kerkhoff, A. J., Swenson, N. G., and Enquist, B. J. (2010). Biological stoichiometry of plant production: metabolism, scaling and ecological response to global change. New Phytol. 186, 593–608. doi: 10.1111/j.1469-8137.2010.03214.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gilman, I. S. and Edwards, E. J. (2020). Crassulacean acid metabolism. Curr. Biol. 30, R57–R62. doi: 10.1016/j.cub.2019.11.073

PubMed Abstract | Crossref Full Text | Google Scholar

Gnanamoorthy, P., Song, Q., Zhao, J., Zhang, Y., Liu, Y., Zhou, W., et al. (2021). Altered albedo dominates the radiative forcing changes in a subtropical forest following an extreme snow event. Global Change Biol. 27, 6192–6205. doi: 10.1111/gcb.15885

PubMed Abstract | Crossref Full Text | Google Scholar

Gotsch, S. G., Williams, C. B., Bicaba, R., Cruz-de Hoyos, R., Darby, A., Davidson, K., et al. (2022). Trade-offs between succulent and non-succulent epiphytes underlie variation in drought tolerance and avoidance. Oecologia 198, 645–661. doi: 10.1007/s00442-022-05140-9

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Y., Yan, Z., Gheyret, G., Zhou, G., Xie, Z., and Tang, Z. (2020). The community-level scaling relationship between leaf nitrogen and phosphorus changes with plant growth, climate and nutrient limitation. J. Ecol. 108, 1276–1286. doi: 10.1111/1365-2745.13369

Crossref Full Text | Google Scholar

Guo, Y., Yan, Z., Zhang, Y.-W., Zhou, G., Xie, Z., and Tang, Z. (2021). Environmental constraints on the inter-genus variation in the scaling relationship between leaf nitrogen and phosphorus concentrations. J. Plant Ecol. 14, 616–627. doi: 10.1093/jpe/rtab016

Crossref Full Text | Google Scholar

Güsewell, S. (2004). N: P ratios in terrestrial plants: variation and functional significance. New Phytol. 164, 243–266. doi: 10.1111/j.1469-8137.2004.01192.x

PubMed Abstract | Crossref Full Text | Google Scholar

Hadi Alkarawi, H. and Zotz, G. (2014). Phytic acid in green leaves. Plant Biol. 16, 697–701. doi: 10.1111/plb.12136

PubMed Abstract | Crossref Full Text | Google Scholar

He, R., Liu, Y., Song, C., Feng, G., and Song, J. (2024). Osmotic regulation beyond nitrate nutrients in plant resistance to stress: a review. Plant Growth Regul. 103, 1–8. doi: 10.1007/s10725-023-01093-y

Crossref Full Text | Google Scholar

Hietz, P., Wagner, K., Nunes Ramos, F., Cabral, J. S., Agudelo, C., Benavides, A. M., et al. (2022). Putting vascular epiphytes on the traits map. J. Ecol. 110, 340–358. doi: 10.1111/1365-2745.13802

Crossref Full Text | Google Scholar

Hietz, P., Wanek, W., Wania, R., and Nadkarni, N. M. (2002). Nitrogen-15 natural abundance in a montane cloud forest canopy as an indicator of nitrogen cycling and epiphyte nutrition. Oecologia 131, 350–355. doi: 10.1007/s00442-002-0896-6

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Liu, W.-Y., and Li, D.-W. (2013). Epiphytic lichens in subtropical forest ecosystems in southwest China: Species diversity and implications for conservation. Biol. Conserv. 159, 88–95. doi: 10.1016/j.biocon.2012.12.027

Crossref Full Text | Google Scholar

Li, J., Yan, X., Zhang, P., Zhuo, Z., Wang, X., Huang, K., et al. (2025). Multiple global change factors alter the scaling of nitrogen to phosphorus in alpine plants. Funct. Ecol. 39, 2044–2055. doi: 10.1111/1365-2435.70087

Crossref Full Text | Google Scholar

Liu, W., Fox, J. E. D., and Xu, Z. (2003). Litterfall and nutrient dynamics in a montane moist evergreen broad-leaved forest in Ailao Mountains, SW China. Plant Ecol. 164, 157–170. doi: 10.1023/A:1021201012950

Crossref Full Text | Google Scholar

Liu, Y., Sun, W., Jia, T., Su, T.-H., Wu, S.-S., Zhou, C.-Y., et al. (2025). Leaf nutrient resorption of vascular epiphytes is regulated by stoichiometry and nutrient limitation control strategies. Plant Cell Environ. 48, 4550–4563. doi: 10.1111/pce.15455

PubMed Abstract | Crossref Full Text | Google Scholar

Niklas, K. J., Owens, T., Reich, P. B., and Cobb, E. D. (2005). Nitrogen/phosphorus leaf stoichiometry and the scaling of plant growth. Ecology Letters 8, 636–642. doi: 10.1111/j.1461-0248.2005.00759.x

Crossref Full Text | Google Scholar

Ouyang, M., Tian, D., Niklas, K. J., Yan, Z., Han, W., Yu, Q., et al. (2024). The scaling of elemental stoichiometry and growth rate over the course of bamboo ontogeny. New Phytol. 241, 1088–1099. doi: 10.1111/nph.19408

PubMed Abstract | Crossref Full Text | Google Scholar

Petter, G., Wagner, K., Wanek, W., Sánchez Delgado, E. J., Zotz, G., Cabral, J. S., et al. (2016). Functional leaf traits of vascular epiphytes: vertical trends within the forest, intra- and interspecific trait variability, and taxonomic signals. Funct. Ecol. 30, 188–198. doi: 10.1111/1365-2435.12490

Crossref Full Text | Google Scholar

Querejeta, J. I., Prieto, I., Torres, P., Campoy, M., Alguacil, M. M., and Roldán, A. (2018). Water-spender strategy is linked to higher leaf nutrient concentrations across plant species colonizing a dry and nutrient-poor epiphytic habitat. Environ. Exp. Bot. 153, 302–310. doi: 10.1016/j.envexpbot.2018.06.007

Crossref Full Text | Google Scholar

Rao, M., Cao, J., Yuan, Q., Fan, H., Ye, D., Liu, L., et al. (2025). Divergent scaling relationship between nitrogen and phosphorus of different shrub plant functional groups under different light and soil nutrient conditions. J. Plant Ecol., rtaf182. doi: 10.1093/jpe/rtaf182

Crossref Full Text | Google Scholar

Reich, P. B., Oleksyn, J., Wright, I. J., Niklas, K. J., Hedin, L., and Elser, J. J. (2010). Evidence of a general 2/3-power law of scaling leaf nitrogen to phosphorus among major plant groups and biomes. Proc. R. Soc. B: Biol. Sci. 277, 877–883. doi: 10.1098/rspb.2009.1818

PubMed Abstract | Crossref Full Text | Google Scholar

Rivas-Ubach, A., Sardans, J., Perez-Trujillo, M., Estiarte, M., and Penuelas, J. (2012). Strong relationship between elemental stoichiometry and metabolome in plants. Proc. Natl. Acad. Sci. 109, 4181–4186. doi: 10.1073/pnas.1116092109

PubMed Abstract | Crossref Full Text | Google Scholar

Sardans, J., Rivas-Ubach, A., and Peñuelas, J. (2012). The elemental stoichiometry of aquatic and terrestrial ecosystems and its relationships with organismic lifestyle and ecosystem structure and function: a review and perspectives. Biogeochemistry 111, 1–39. doi: 10.1007/s10533-011-9640-9

Crossref Full Text | Google Scholar

Schleuss, P. M., Widdig, M., Heintz-Buschart, A., Kirkman, K., and Spohn, M. (2020). Interactions of nitrogen and phosphorus cycling promote P acquisition and explain synergistic plant-growth responses. Ecology 101, e03003. doi: 10.1002/ecy.3003

PubMed Abstract | Crossref Full Text | Google Scholar

Sterner, R. W. and Elser, J. J. (2002). Ecological stoichiometry: the biology of elements from molecules to the biosphere (Princeton University Press). doi: 10.1515/9781400885695

Crossref Full Text | Google Scholar

Sulaiman, S. P., Metali, F., and Cicuzza, D. (2023). Tropical terrestrial and epiphytic ferns have different leaf stoichiometry with ecological implications. Plant Biosyst. - Int. J. Dealing all Aspects Plant Biol. 157, 899–906. doi: 10.1080/11263504.2023.2229852

Crossref Full Text | Google Scholar

Tian, D., Yan, Z., Niklas, K. J., Han, W., Kattge, J., Reich, P. B., et al. (2018). Global leaf nitrogen and phosphorus stoichiometry and their scaling exponent. Natl. Sci. Rev. 5, 728–739. doi: 10.1093/nsr/nwx142

Crossref Full Text | Google Scholar

Warton, D. I., Duursma, R. A., Falster, D. S., and Taskinen, S. (2012). smatr 3– an R package for estimation and inference about allometric lines. Methods Ecol. Evol. 3, 257–259. doi: 10.1111/j.2041-210X.2011.00153.x

Crossref Full Text | Google Scholar

Warton, D. I., Wright, I. J., Falster, D. S., and Westoby, M. (2006). Bivariate line-fitting methods for allometry. Biol. Rev. 81, 259–291. doi: 10.1017/S1464793106007007

PubMed Abstract | Crossref Full Text | Google Scholar

Wen, H.-D., Lin, L.-X., Yang, J., Hu, Y.-H., Cao, M., Liu, Y.-H., et al. (2018). Species composition and community structure of a 20 hm2 plot of mid-mountain moist evergreen broad-leaved forest on the Mts. AiLaoshan, Yunnan Province, China. Chin. J. Plant Ecol. 42, 419–429. doi: 10.17521/cjpe.2017.0272

Crossref Full Text | Google Scholar

Winkler, U. and Zotz, G. (2009). Highly efficient uptake of phosphorus in epiphytic bromeliads. Ann. Bot. 103, 477–484. doi: 10.1093/aob/mcn231

PubMed Abstract | Crossref Full Text | Google Scholar

Wright, I. J., Reich, P. B., Cornelissen, J. H. C., Falster, D. S., Garnier, E., Hikosaka, K., et al. (2005). Assessing the generality of global leaf trait relationships. New Phytol. 166, 485–496. doi: 10.1111/j.1469-8137.2005.01349.x

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, Y., Liu, W.-Y., Lu, H.-Z., Li, S., Shen, Y.-X., Liu, W.-G., et al. (2020). Stoichiometric and isotopic flexibility: facultative epiphytes exploit rock and bark interchangeably. Environ. Exp. Bot. 179, 104208. doi: 10.1016/j.envexpbot.2020.104208

Crossref Full Text | Google Scholar

Yu, X., Zhang, Y., Zhao, X., Li, J., Yu, X., Zhang, Y., et al. (2023). Systemic effects of the vapor pressure deficit on the physiology and productivity of protected vegetables. VR 3. doi: 10.48130/VR-2023-0020

Crossref Full Text | Google Scholar

Yuan, Z. Y. and Chen, H. Y. H. (2015). Decoupling of nitrogen and phosphorus in terrestrial plants associated with global changes. Nat. Clim Change 5, 465–469. doi: 10.1038/nclimate2549

Crossref Full Text | Google Scholar

Zhang, T., Liu, W., Hu, T., Tang, D., Mo, Y., and Wu, Y. (2020). Divergent adaptation strategies of vascular facultative epiphytes to bark and soil habitats: insights from stoichiometry. Forests 12, 16. doi: 10.3390/f12010016

Crossref Full Text | Google Scholar

Zhao, M., Geekiyanage, N., Xu, J., Khin, M. M., Nurdiana, D. R., Paudel, E., et al. (2015). Structure of the epiphyte community in a tropical montane forest in SW China. PloS One 10, e0122210. doi: 10.1371/journal.pone.0122210

PubMed Abstract | Crossref Full Text | Google Scholar

Zotz, G. (2004). The resorption of phosphorus is greater than that of nitrogen in senescing leaves of vascular epiphytes from lowland Panama. J. Trop. Ecol. 20, 693–696. doi: 10.1017/S0266467404001889

Crossref Full Text | Google Scholar

Zotz, G. (2016). Plants on plants – the biology of vascular epiphytes (Cham: Springer International Publishing). doi: 10.1007/978-3-319-39237-0

Crossref Full Text | Google Scholar

Zotz, G. and Hietz, P. (2001). The physiological ecology of vascular epiphytes: current knowledge, open questions. J. Exp. Bot. 52, 2067–2078. doi: 10.1093/jexbot/52.364.2067

PubMed Abstract | Crossref Full Text | Google Scholar

Zotz, G., Hietz, P., and Einzmann, H. J. R. (2021a). “Functional ecology of vascular epiphytes,” in Annual plant reviews online (John Wiley & Sons, Ltd), 869–906. doi: 10.1002/9781119312994.apr0777

Crossref Full Text | Google Scholar

Zotz, G., Weigelt, P., Kessler, M., Kreft, H., and Taylor, A. (2021b). EpiList 1.0: a global checklist of vascular epiphytes. Ecology 102, e03326. doi: 10.1002/ecy.3326

PubMed Abstract | Crossref Full Text | Google Scholar