Phenology determines leaf functional traits across Rhododendron species in the Sikkim Himalaya

Elevation gradients provide an ideal setting to infer species' functional trait responses to predicted future climate change. In plants, leaf functional traits help assess their capacity to cope with varying resources. Variation in abiotic conditions over short vertical distances can influence plant phenology, particularly leafing and flowering durations, and leaf functional traits at both inter- and intra-specific levels. However, studies examining relationships between leaf functional traits and phenology duration along elevation gradients are limited. We tested the relationship between leaf size, leaf thickness, specific leaf area, and leafing durations in 10 Rhododendron species in the Sikkim Himalaya. All the investigated traits varied significantly across species, but intra-specific variation in functional traits was observed only among a few. Leaf size and thickness showed significant negative relationships with elevation and a comparative phylogenetic method exhibited a strong relationship between leaf traits and leafing duration. We observed higher leaf thickness and size in species with longer leafing durations and less overlap in leafing and flowering durations. In contrast, species with shorter leafing durations and relatively more overlap in their flowering and leafing durations showed lower leaf thickness and leaf size. Leaf traits such as leaf thickness and leaf size also exhibited a strong phylogenetic signal across 10 Rhododendron species. Overall, from our findings, we infer that along an elevation gradient, the magnitude of leaf trait responses to future increases in temperature may vary depending on species phenology durations and phylogeny.