J Plant Ecol ›› Advance articles     DOI:10.1093/jpe/rtag240

    Next Articles

Climatic Seasonality and Topographic Complexity Shape Canary Islands Plant Growth Form Distributions

Lucas S. Jansen1, *, Ryan F.A. Brewer1, 2, *, **, Juli Caujapé-Castells3, Luis Valente1, 2, ∧, FredericLens1, 4, ∧   

  1. 1. Naturalis Biodiversity Center, Darwinweg 2, 2333 CR Leiden, The Netherlands;
    2. Groningen Institute for Evolutionary Life Sciences, University of Groningen, Box 11103, 9700 CC Groningen, The Netherlands;
    3. Department of Molecular Biodiversity & DNA Bank of the Canarian Flora, Jardín Botánico Canario 'Viera y Clavijo'-Unidad Asociada CSIC, Cabildo de Gran Canaria, Las Palmas, Gran Canaria, Spain;
    4. Institute of Biology Leiden, Leiden University, Sylviusweg 72, Leiden 2333 BE, The Netherlands
    * Shared First Authors
    ** Corresponding author: Ryan F. A. Brewer (R.F.Brewer@outlook.com)
    ∧ Shared Senior Authors
  • Received:2025-12-18 Revised:2026-07-27 Accepted:2026-09-06 Published:2026-10-17
  • Supported by:
    This work was supported by the Dutch Research Council [OCENW.KLEIN.498] to F.L.

Abstract: Plant growth forms reflect the environmental conditions under which they evolved and persist. Yet, previous studies linking environmental conditions with growth form have overlooked that different growth form types arose through distinct evolutionary pathways; some woody angiosperm lineages evolved woodiness from a herbaceous ancestor, whereas others have retained their woodiness throughout evolutionary history. The Canary Islands are a global hotspot of insular woodiness (woodiness evolved on islands from herbaceous colonisers), raising the question of whether insular woody species occur in different environmental conditions than other growth forms. Here, we use a custom-built pipeline to extract and filter spatial data to analyse the environmental correlates of growth forms in the Canary Islands non-monocot angiosperm flora, accounting for their distinct evolutionary histories and differences in life span. To assess whether and how five growth forms (annual herbaceous, perennial herbaceous, ancestral woodiness, derived woodiness, and insular woodiness) differ in environmental niche space, we applied phylogenetic PCAs to a dataset of ∼600 native non-monocot Canary Island angiosperm species. Unexpectedly, herbaceous and woody growth forms highly overlap in niche space, with subtle differences explained mostly by climatic seasonality (temperature and precipitation) and topographic complexity. Using phylogenetic ANOVAs, we identified significant pairwise differences between growth forms across five environmental variables, with differentiation between combined woody and herbaceous niches, and between the three woody growth forms. Our phylogenetically inferred macroecological approach identifies terrain rugosity as a key environmental variable influencing niche space, and suggests aridity may be less important for determining insular woodiness distributions than previously hypothesised.

This study investigates how plant growth forms relate to environmental conditions in the Canary Islands, a global hotspot of insular woodiness. Using a dataset of ∼600 native non-monocot angiosperm species classified into five growth forms (annual herbaceous, perennial herbaceous, ancestral woody, derived woody, and insular woody), and accounting for their distinct evolutionary histories and lifespans, the authors applied phylogenetic PCAs and ANOVAs to examine environmental niche differentiation. Contrary to expectations, herbaceous and woody growth forms showed substantial niche overlap, with subtle differences driven mainly by climatic seasonality (temperature and precipitation) and topographic complexity. The analysis identified terrain rugosity as a key environmental variable shaping niche space and suggested that aridity may be less important for insular woodiness distributions than previously hypothesised.

Key words: Canary Islands, environmental niche, growth form, phylogenetic ANOVA, phylogenetic principal component analysis