Leaf size in conifers: global associations with climate and evolutionary history
Data files
Jun 17, 2026 version files 4.20 MB
-
climate_data.csv
3.16 MB
-
conifer_leaf.csv
50.85 KB
-
distribution_map.R
2.78 KB
-
functions_Bandow.R
24.69 KB
-
gheatmaps.R
11.74 KB
-
Leslie_et_al_AJB_2018.nex
883.04 KB
-
mrpmm.R
23.44 KB
-
node_labels.csv
31.56 KB
-
README.md
3.03 KB
-
regressions.R
7.65 KB
Abstract
Global analyses of leaf size suggest that large leaves predominate in favourably warm and wet climates, and small leaves occur at climatic extremes. However, these patterns are dominated by data from angiosperms, and may obscure drivers of leaf size variation in older, less speciose groups, such as conifers. Here we employ a novel modelling framework, multi-response phylogenetic mixed models (MRPMM), that identifies trait correlations at both phylogenetic and phylogenetically independent levels to investigate how climate influences leaf size evolution across conifers. We show overall patterns for all conifers, and focus on three groups with distinctive leaf architectures: scale-leaved Cupressaceae, needle-leaved Pinus, and broad-leaved Podocarpus.
We found moderate to strong phylogenetic signal in conifer leaf and climate niche traits. Phylogenetic relationships explained most of the association between leaf size and climate, with temperature revealed as a stronger driver than dry season precipitation, indicating deep-time co-evolutionary associations. These patterns reflect trade-offs associated with leaf hydraulic architecture under contrasting selection pressures. In single-veined leaves, lateral water transport and thus leaf width is constrained, yet this narrow leaf form can be advantageous for survival under climatic extremes such as drought and freezing. In contrast, anatomical innovations (like accessory transfusion tissue in Podocarpus) allow for broader leaves which are favourable in competitive environments, while also potentially making leaves more vulnerable to climate extremes.
Our results support previous evidence for phylogenetic niche conservatism in conifers, where species tend to track their ancestral climatic preferences rather than adapting to new environments. This conservatism, likely controlled by leaf hydraulic architecture, results in strong evolutionary constraints on current bioclimatic distributions and potential responses to changing climates in conifers. This study also highlights the importance of considering phylogenetic impacts on functional trait evolution, especially in evolutionarily conservative groups like conifers.
Dataset DOI: 10.5061/dryad.3tx95x6vw
Corresponding Author Information:
Name: Katya Isabelle Bandow
Institution: University of Tasmania
Email: katya.bandow@utas.edu.au
Data:
Phylogenetic tree (from Leslie et al. 2018) : Leslie_et_al_AJB_2018.nex and node_labels.csv
Leaf size data (various sources, see manuscript): conifer_leaf.csv
Species climate data (extracted from Brown et al. 2021): climate_data.csv
Analysis Code:
All analysis was performed in R Version 4.3.1
Multi-Response Phylogenetic Mixed Models (MRPMM) using the R package MCMCglmm.
Functions to call upon for MRPMM results: functions_Bandow.R
Running models, extracting model outputs, model figures and model diagnositics: mrpmm.R
Non-MRPMM analysis:
Species distribution map: distribution_map.R
Non-phylogenetic leaf size-climatic correlations: regressions.R
Phylogenetic heatmap: gheatmaps.R
Abbreviations:
| Entry | Class | Explanation |
|---|---|---|
| name | character | Species binomial name |
| family | character | Taxonomic family |
| genus | character | Taxonomic genus |
| country | character | Country of record |
| decimalLatitude | numeric | Latitude of observation |
| decimalLongitude | numeric | Longitude of observation |
| leaf_length_min | numeric | Minimum leaf length (mm) |
| leaf_length_max | numeric | Maximum leaf length (mm) |
| length_mean | numeric | Mean leaf length (mm) |
| leaf_width_min | numeric | Minimum leaf width (mm) |
| leaf_width_max | numeric | Maximum leaf width (mm) |
| width_mean | numeric | Mean leaf width (mm) |
| leaf1 | character | Juvenile leaf form |
| leaf2 | character | Adult leaf form |
For detailed observation metadata see Brown et al. 2021.
Climatic variables sourced from WorldClim2, 30":
| Entry | Class | Explanation |
|---|---|---|
| MinTCM | numeric | Bioclimatic variable BIO6: Minimum Temperature of the Coldest Month (°C) |
| PDQ | numeric | Bioclimatic variable BIO17: Precipitation of the Driest Quarter (mm) |
| MaxT | numeric | Bioclimatic variable BIO5: Maximum Temperature of the Warmest Month (°C) |
