Evolution of the leaf economics spectrum in herbs: Evidence from environmental divergences in leaf physiology across Helianthus (Asteraceae)

The leaf economics spectrum (LES) describes a major axis of plant functional trait variation worldwide, defining suites of leaf traits aligned with resource‐acquisitive to resource‐conservative ecological strategies. The LES has been interpreted to arise from leaf‐level trade‐offs among ecophysiological traits common to all plants. However, it has been suggested that the defining leaf‐level trade‐offs of the LES may not hold within specific functional groups (e.g., herbs) nor within many groups of closely related species, which challenges the usefulness of the LES paradigm across evolutionary scales. Here, we examine the evolution of the LES across 28 species of the diverse herbaceous genus Helianthus (the sunflowers), which occupies a wide range of habitats and climate variation across North America. Using a phylogenetic comparative approach, we find repeated evolution of more resource‐acquisitive LES strategies in cooler, drier, and more fertile environments. We also find macroevolutionary correlations among LES traits that recapitulate aspects of the global LES, but with one major difference: leaf mass per area is uncorrelated with leaf lifespan. This indicates that whole‐plant processes likely drive variation in leaf lifespan across Helianthus, rather than leaf‐level trade‐offs. These results suggest that LES patterns do not reflect universal physiological trade‐offs at small evolutionary scales.     

Within‐species patterns challenge our understanding of the leaf economics spectrum

The utility of plant functional traits for predictive ecology relies on our ability to interpret trait variation across multiple taxonomic and ecological scales. Using extensive data sets of trait variation within species, across species and across communities, we analysed whether and at what scales leaf economics spectrum (LES) traits show predicted trait–trait covariation. We found that most variation in LES traits is often, but not universally, at high taxonomic levels (between families or genera in a family). However, we found that trait covariation shows distinct taxonomic scale dependence, with some trait correlations showing opposite signs within vs. across species. LES traits responded independently to environmental gradients within species, with few shared environmental responses across traits or across scales. We conclude that, at small taxonomic scales, plasticity may obscure or reverse the broad evolutionary linkages between leaf traits, meaning that variation in LES traits cannot always be interpreted as differences in resource use strategy.     

Trait variation and integration across scales: is the leaf economic spectrum present at local scales?

Trait‐based approaches have taken an increasingly dominant role in community ecology. Although trait‐based strategy dimensions such as the leaf economic spectrum (LES) have been identified primarily at global‐scales, trait variation at the community scale is often interpreted in this context. Here we argue from several lines of evidence that a research priority should be to determine whether global‐scale trait relationships hold at more local scales. We review recent literature assessing trait variation at smaller scales, and then present a case study exploring the relationship between the correlation strength of leaf traits and their similarity in variation structure across ecological scales. We find that the correlation strength between pairs of leaf traits does not predict whether the traits respond similarly to different drivers of variation. Instead, correlation strength only sets an upper bound to the dissimilarity in trait variation structure. With moderate correlation strengths, LES traits largely retain the ability to respond independently to different drivers of phenotypic variation at different scales. Recent literature and our results suggest that LES relationships may not hold at local scales. Clarifying under what conditions and at which scales the LES is consistently expressed is necessary for us to make the most of the emerging trait toolbox.     

Biogeographic constraints on the world‐wide leaf economics spectrum

Aim The world‐wide leaf economic spectrum (LES) describes tight coordination of leaf traits across global floras, reported to date as being largely independent of phylogeny and biogeography. Here, we present and test an alternative, historical perspective that predicts that biogeography places significant constraints on global trait evolution. These hypothesized constraints could lead to important deviations in leaf trait relationships between isolated floras that were influenced by different magnitudes of genetic constraint and selection. Location Global, including floristic regions of the Northern and Southern Hemispheres, eastern North America, East Asia (EAS), the Hawaiian Islands and tropical mainland floras. Methods We use a large leaf‐trait database (GLOPNET) and species native distribution data to test for variation in leaf trait relationships modulated by floristic region, controlling for climatic differences. Standardized major axis analyses were used to evaluate biogeographic effects on bivariate relationships between LES traits, including relationships of photosynthetic capacity and dark respiration rate (A_mass–R_d‐mass), leaf lifespan and mass per area ratio (LL–LMA), and photosynthetic capacity and nitrogen content (A_mass–N_mass). Results Independent of climate or biome, floras of different evolutionary histories exhibited different leaf trait allometries. Floras of the Northern Hemisphere exhibited greater rates of return on resource investment (steeper slopes for the trait relationships analysed), and the more diverse temperate EAS flora exhibited greater slopes or intercepts in leaf trait relationships, with the exception of the A_mass–N_mass relationship. In contrast to our hypothesis, plants of the floristically isolated Hawaiian Islands exhibited a similar A_mass–N_mass relationship to those of mainland tropical regions. Main conclusions Differences in leaf trait allometries among global floristic regions support a historical perspective in understanding leaf trait relationships and suggest that independent floras can exhibit different tradeoffs in resource capture strategies.     

The leaf economics spectrum’s morning coffee: plant size-dependent changes in leaf traits and reproductive onset in a perennial tree crop

Background and AimsSize-dependent changes in plant traits are an important source of intraspecific trait variation. However, there are few studies that have tested if leaf trait co-variation and/or trade-offs follow a within-genotype leaf economics spectrum (LES) related to plant size and reproductive onset. To our knowledge, there are no studies on any plant species that have tested whether or not the shape of a within-genotype LES that describes how traits covary across whole plant sizes, is the same as the shape of a within-genotype LES that represents environmentally driven trait plasticity.MethodsWe quantified size-dependent variation in eight leaf traits in a single coffee genotype (Coffea arabica var. Caturra) in managed agroecosystems with different environmental conditions (light and fertilization treatments), and evaluated these patterns with respect to reproductive onset. We also evaluated if trait covariation along a within-genotype plant-size LES differed from a within-genotype environmental LES defined with trait data from coffee growing in different environmental conditions.Key ResultsLeaf economics traits related to resource acquisition – maximum photosynthetic rates (A) and mass-based leaf nitrogen (N) concentrations – declined linearly with plant size. Structural traits – leaf mass, leaf thickness, and leaf mass per unit area (LMA) – and leaf area increased with plant size beyond reproductive onset, then declined in larger plants. Three primary LES traits (mass-based A, leaf N and LMA) covaried across a within-genotype plant-size LES, with plants moving towards the ‘resource-conserving’ end of the LES as they grow larger; in coffee these patterns were nearly identical to a within-genotype environmental LES.ConclusionsOur results demonstrate that a plant-size LES exists within a single genotype. Our findings indicate that in managed agroecosystems where resource availability is high the role of reproductive onset in driving within-genotype trait variability, and the strength of covariation and trade-offs among LES traits, are less pronounced compared with plants in natural systems. The consistency in trait covariation in coffee along both plant-size and environmental LES axes indicates strong constraints on leaf form and function that exist within plant genotypes.     

Geographical variation and the role of climate in leaf traits of a relict tree species across its distribution in China

• Intraspecific trait variation and trait–climate relationships are crucial for understanding a species’ response to climate change. However, these phenomena have rarely been studied for tree species. Euptelea pleiospermum is a relict tree species with a wide distribution in China that offers a novel opportunity to examine such relationships.
• Here, we measured 13 leaf traits of E. pleiospermum in 20 sites across its natural distribution in China. We investigated the extent of trait variation at local and regional scales, and developed geographic and climate models to explain trait variation at the regional scale.
• We documented intraspecific trait variation among leaf traits of E. pleiospermum at local and regional scales. Five traits exhibited relatively high trait variation: leaf area, leaf density and three leaf economic traits (leaf dry matter content, specific leaf area [SLA] and leaf phosphorus concentration). Significant trait–geography correlations were mediated by local climate. Most leaf trait variation could be explained (from 24% to 64%) by geographic or climate variables, except leaf width, leaf thickness, leaf dry matter content and leaf length–width ratio. Latitude and temperature were the strongest predictors of trait variation throughout the distribution of E. pleiospermum in China, and temperature explained more leaf trait variation than precipitation. In particular, we showed that leaves had longer petiole lengths, higher SLA and lower densities in northern E. pleiospermum populations. We suggest that northern E. pleiospermum populations are adapting to higher latitudinal environments via high growth rate (higher SLA) and low construction investment strategies (lower leaf densities), benefitting northern migration.
• Overall, we demonstrate that intraspecific trait variation reflects E. pleiospermum response to the local environment. We call for consideration of intraspecific trait variation to examine specific climate response questions. In addition, provenance experiments using widely distributed species are needed to separate trait variation resulting from genetic differentiation and plastic responses to environmental change.

     

Integrating trait and phylogenetic distances to assess scale‐dependent community assembly processes

Biodiversity is structured by multiple mechanisms that are dependent, at least in part, on ecological similarities and differences among species. Integrating traits and phylogenies in diversity metrics may provide deeper insight into community assembly processes across spatial scales. However, different traits are influenced by processes at different spatial scales, and it is not clear how trait‐spatial scale mismatches skew our ability to detect assembly patterns. An additional complexity is how phylogenetic distances, which might capture unmeasured traits, reflect spatially dependent processes. Here we analyze a freshwater zooplankton dataset from 91 ponds and show that different traits are associated with processes at different spatial scales. We first assessed the response of individual traits to processes at both α‐ and β‐scales, and then quantified the power of different combinations of traits and phylogenetic distances to reveal environmental and spatial drivers of α‐ and β‐diversity. We found that explanatory power was maximised when we accounted for environmental and spatial drivers with single, but different traits for α‐ and β‐diversity. Using the most appropriate trait for each spatial scale outperformed phylogenetic information, but phylogenetic information outperformed the same traits when these were used at the wrong spatial scale, and all outperformed taxonomic analyses that ignore trait and phylogenetic information. We demonstrate that accounting for species’ similarities and differences provides important information about dominant assembly mechanisms at different spatial scales, and that phylogeny is especially useful when measured traits are uninformative at a given spatial scale or when there is lack of trait data. Our study also indicates, however, that trait‐scale mismatches among phylogenetically conserved traits may affect the performance of phylogenetic indices compared to indices that account only for the best single trait at each spatial scale.     

Trait biogeography of marine copepods – an analysis across scales

Functional traits, rather than taxonomic identity, determine the fitness of individuals in their environment: traits of marine organisms are therefore expected to vary across the global ocean as a function of the environment. Here, we quantify such spatial and seasonal variations based on extensive empirical data and present the first global biogeography of key traits (body size, feeding mode, relative offspring size and myelination) for pelagic copepods, the major group of marine zooplankton. We identify strong patterns with latitude, season and between ocean basins that are partially (c. 50%) explained by key environmental drivers. Body size, for example decreases with temperature, confirming the temperature‐size rule, but surprisingly also with productivity, possibly driven by food‐chain length and size‐selective predation. Patterns unrelated to environmental predictors may originate from phylogenetic clustering. Our maps can be used as a test‐bed for trait‐based mechanistic models and to inspire next‐generation biogeochemical models.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

Fine root tradeoffs between nitrogen concentration and xylem vessel traits preclude unified whole‐plant resource strategies in Helianthus

Recent work suggests variation in plant growth strategies is governed by a tradeoff in resource acquisition and use, ranging from a rapid resource acquisition strategy to a resource‐conservative strategy. While evidence for this tradeoff has been found in leaves, knowledge of root trait strategies, and whether they reflect adaptive differentiation across environments, is limited. In the greenhouse, we investigated variation in fine root morphology (specific root length and tissue density), chemistry (nitrogen concentration and carbon:nitrogen), and anatomy (root cross‐sectional traits) in populations of 26 Helianthus species and sister Phoebanthus tenuifolius. We also compared root trait variation in this study with leaf trait variation previously reported in a parallel study of these populations. Root traits varied widely and exhibited little phylogenetic signal, suggesting high evolutionary lability. Specific root length and root tissue density were weakly negatively correlated, but neither was associated with root nitrogen, providing little support for a single axis of root trait covariation. Correlations between traits measured in the greenhouse and native site characteristics were generally weak, suggesting a variety of equally viable root trait combinations exist within and across environments. However, high root nitrogen was associated with lower xylem vessel number and cross‐sectional area, suggesting a tradeoff between nutrient investment and water transport capacity. This led to correlations between root and leaf traits that were not always consistent with an acquisition–conservation tradeoff at the whole‐plant level. Given that roots must balance acquisition of water and nutrients with functions like anchorage, exudation, and microbial symbioses, the varied evidence for root trait covariation likely reflects the complexity of interacting selection pressures belowground. Similarly, the lack of evidence for a single acquisition–conservation tradeoff at the whole‐plant level likely reflects the vastly different selection pressures shaping roots and leaves, and the resources they are optimized to obtain.     

Common trade‐offs between xylem resistance to cavitation and other physiological traits do not hold among unrelated Populus deltoides ×Populus nigra hybrids

We examined the relationships between xylem resistance to cavitation and 16 structural and functional traits across eight unrelated Populus deltoides×Populus nigra genotypes grown under two contrasting water regimes. The xylem water potential inducing 50% loss of hydraulic conductance (Ψ₅₀) varied from ?1.60 to ?2.40 MPa. Drought‐acclimated trees displayed a safer xylem, although the extent of the response was largely genotype dependant, with Ψ₅₀ being decreased by as far as 0.60 MPa. At the tissue level, there was no clear relationship between xylem safety and either xylem water transport efficiency or xylem biomechanics; the only structural trait to be strongly associated with Ψ₅₀ was the double vessel wall thickness, genotypes exhibiting a thicker double wall being more resistant. At the leaf level, increased cavitation resistance was associated with decreased stomatal conductance, while no relationship could be identified with traits associated with carbon uptake or bulk leaf carbon isotope discrimination, a surrogate of intrinsic water‐use efficiency. At the whole‐plant level, increased safety was associated with higher shoot growth potential under well‐irrigated regime only. We conclude that common trade‐offs between xylem resistance to cavitation and other physiological traits that are observed across species may not necessarily hold true at narrower scales.     

Australian Tropical and Subtropical Rain Forest Community Assembly: Phylogeny,Functional Biogeography,and Environmental Gradients

To compare community assemblage patterns in tropical northeastern and subtropical central eastern Australia across selected gradients and scales, we tested the relationship of species traits with phylogenetic structure, and niche breadth. We considered phylogenetic relationships across current‐day species in assemblages in relation to rain forest species pool sizes, and trait values along gradients including elevation and latitude. Trait values were quantified across scales for seed size, leaf area, wood density and maximum height at maturity for 1137 species and 596 assemblages using trait gradient analysis (TGA). Local assemblages of subtropical species had narrower trait ranges, and higher niche breadth values than corresponding assemblages of tropical species. Leaf size and seed size increased at low latitudes, and community phylogenetic structure was most strongly correlated with seed traits in the subtropics, reflecting dispersal and re‐colonization processes. Elevation accounted for little of the variance in community phylogenetic structure or trait variation across local and regional scales. Stable moist forest areas retained many species from ancestral rain forest lineages across a range of temporally conserved habitats; species within assemblages were less related; and rain forest assemblages had higher functional diversity, but lower niche breadth. This suggests that on average, assemblages of species in stable areas had greater trait variation and narrower distributions. Historic and recent rain forest contraction and re‐expansion can result in recolonized areas that are dominated by species that are more related (phylogenetically) than by chance, have smaller, widely dispersed seeds, and greater niche breadth (broader distributions).     

Constraints on trait combinations explain climatic drivers of biodiversity: the importance of trait covariance in community assembly

Trade‐offs maintain diversity and structure communities along environmental gradients. Theory indicates that if covariance among functional traits sets a limit on the number of viable trait combinations in a given environment, then communities with strong multidimensional trait constraints should exhibit low species diversity. We tested this prediction in winter annual plant assemblages along an aridity gradient using multilevel structural equation modelling. Univariate and multivariate functional diversity measures were poorly explained by aridity, and were surprisingly poor predictors of community richness. By contrast, the covariance between maximum height and seed mass strengthened along the aridity gradient, and was strongly associated with richness declines. Community richness had a positive effect on local neighbourhood richness, indicating that climate effects on trait covariance indirectly influence diversity at local scales. We present clear empirical evidence that declines in species richness along gradients of environmental stress can be due to increasing constraints on multidimensional phenotypes.     

Common ancestry or environmental trait filters: cross‐continental comparisons of trait–habitat relationships in tropical anuran amphibian assemblages

Aim To investigate whether trait–habitat relations in biological communities converge across three global regions. The goal is to assess the role of habitat templets in shaping trait assemblages when different assembly mechanisms are operating and to test whether trait–habitat relations reflect a common evolutionary history or environmental trait filters. Location Guiana Shield, South America; Upper Guinea Forest Block, West Africa; Borneo rain forests, Southeast Asia. Methods We compared large anuran amphibian data sets at both the regional and cross‐continental scale. We applied a combination of three‐table ordinations (RLQ) and permutation model‐based multivariate fourth‐corner statistics to test for trait–habitat relationships at both scales and used phylogenetic comparative methods to quantify phylogenetic signal in traits that enter these analyses. Results Despite the existence of significant trait–habitat links and congruent trait patterns, we did not find evidence for the existence of a universal trait–habitat relationship at the assemblage level and no clear sign for cross‐continental convergence of trait–habitat relations. Patterns rather varied between continents. Despite the fact that a number of traits were conserved across phylogenies, the phylogenetic signal varied between regions. Trait–habitat relations therefore not only reflect a common evolutionary history, but also more recently operating environmental trait filters that ultimately determine the trait composition in regional assemblages. Main conclusions Integrating trait–habitat links into analyses of biological assemblages can enhance the predictive power and general application of species assembly rules in community and macroecology, particularly when phylogenetic comparative methods are simultaneously applied. However, in order to predict trait composition based on habitat templets, trait–habitat links cannot be assumed to be universal but rather have to be individually established in different regions prior to model building. Only then can direct trait‐based approaches be useful tools for predicting fundamental community patterns.     

Calcium plus magnesium indicates digestibility: the significance of the second major axis of plant chemical variation for ecological processes

Plant variation in nutrient concentrations encompasses two major axes. The first is connected to nitrogen (N) and phosphorus (P), reflects growth rate and has been designated as the leaf economics spectrum (LES) while the second follows the gradient in calcium (Ca) and magnesium (Mg) and mirrors cell structural differences. Here, we tested in grasslands whether the sum Ca + Mg concentrations is a better indicator of digestibility than LES constituents. Structural equation modelling revealed that the total effect size of N (0.30) on digestibility was much lower than that of Ca + Mg (0.58). The N effect originated predominantly from sampling date (biomass ageing), while the Ca + Mg effect largely from phylogenetic composition (proportion of monocots). Thus, plant variation in partially substitutable divalent cations seems to play a significant role in biomass digestion by ruminants. This finding contests, together with litter decomposition studies, the prominent role of the LES for understanding both fundamental ecological processes.     

Leaf trait integration mediates species richness variation in a species-rich neotropical forest domain

Environmental gradients are known to drive changes in mean trait values, but changes in the trait integration strength across local communities are less well understood, particularly with regard to possible links with species richness variation. Here, we tested if climate, soil, and topography gradients drive species richness indirectly via constraints on trait integration in the Atlantic Forest of South America. We evaluated seven traits (from leaf, wood, seed, and plant size) of 1456 species occurring across 84 local communities. Generalized least square models and a path model were applied to test direct and indirect relationships. Correlations were higher between leaf traits (average r?=?0.28) and lower when other traits were included (average r?=?0.16). In line with this result, species richness was related to a multivariate index of interspecific trait integration (ITI) computed for leaf traits, but not to the ITI for all the seven traits. Abiotic gradients influenced species richness both directly and indirectly through the leaf trait integration. A total of 33% and 26% of the variation in species richness and ITI, respectively, were explained by the models, with climatic conditions showing higher contribution than topographic and edaphic factors. These results support a significant but reduced environmental selection role behind the trait-based community assembly and may suggest that other processes are involved in the constrain of trait integration at larger spatial scales. In addition, different directional trends in trait–trait relationships across local communities suggest that global trait relationships may not necessarily hold at local contexts.

     

Greater than the sum of the parts: how the species composition in different forest strata influence ecosystem function

The mechanisms underpinning forest biodiversity‐ecosystem function relationships remain unresolved. Yet, in heterogeneous forests, ecosystem function of different strata could be associated with traits or evolutionary relationships differently. Here, we integrate phylogenies and traits to evaluate the effects of elevational diversity on above‐ground biomass across forest strata and spatial scales. Community‐weighted means of height and leaf phosphorous concentration and functional diversity in specific leaf area exhibited positive correlations with tree biomass, suggesting that both positive selection effects and complementarity occur. However, high shrub biomass is associated with greater dissimilarity in seed mass and multidimensional trait space, while species richness or phylogenetic diversity is the most important predictor for herbaceous biomass, indicating that species complementarity is especially important for understory function. The strength of diversity‐biomass relationships increases at larger spatial scales. We conclude that strata‐ and scale‐ dependent assessments of community structure and function are needed to fully understand how biodiversity influences ecosystem function.     

Highlighting convergent evolution in morphological traits in response to climatic gradient in African tropical tree species: The case of genus Guibourtia Benn.

相似文献   

Global patterns of intraspecific leaf trait responses to elevation

Elevational gradients are often used to quantify how traits of plant species respond to abiotic and biotic environmental variations. Yet, such analyses are frequently restricted spatially and applied along single slopes or mountain ranges. Since we know little on the response of intraspecific leaf traits to elevation across the globe, we here perform a global meta‐analysis of leaf traits in 109 plant species located in 4 continents and reported in 71 studies published between 1983 and 2018. We quantified the intraspecific change in seven morpho‐ecophysiological leaf traits along global elevational gradients: specific leaf area (SLA), leaf mass per area (LMA), leaf area (LA), nitrogen concentration per unit of area (Narea), nitrogen concentration per unit mass (Nmass), phosphorous concentration per unit mass (Pmass) and carbon isotope composition (δ¹³C). We found LMA, Narea, Nmass and δ¹³C to significantly increase and SLA to decrease with increasing elevation. Conversely, LA and Pmass showed no significant pattern with elevation worldwide. We found significantly larger increase in Narea, Nmass, Pmass and δ¹³C with elevation in warmer regions. Larger responses to increasing elevation were apparent for SLA of herbaceous compared to woody species, but not for the other traits. Finally, we also detected evidences of covariation across morphological and physiological traits within the same elevational gradient. In sum, we demonstrate that there are common cross‐species patterns of intraspecific leaf trait variation across elevational gradients worldwide. Irrespective of whether such variation is genetically determined via local adaptation or attributed to phenotypic plasticity, the leaf trait patterns quantified here suggest that plant species are adapted to live on a range of temperature conditions. Since the distribution of mountain biota is predominantly shifting upslope in response to changes in environmental conditions, our results are important to further our understanding of how plants species of mountain ecosystems adapt to global environmental change.     

Expression of anatomical leaf traits in homoploid hybrids between deciduous and evergreen species of Vaccinium

We investigated the anatomical expression of leaf traits in hybrids between evergreen Vaccinium vitis‐idaea and deciduous V. myrtillus. We compared parents from four populations with their respective F1 hybrids and tested whether (i) transgression can be the source of novel anatomical traits in hybrids; (ii) expression of transgressive traits is more probable for traits with similar values in parents and intermediate for more distinct values, as predicted by theory; and (iii) independent origin of hybrids leads to identical trait expression profiles among populations. We found that anatomical leaf traits can be divided into four categories based on their similarity to parents: intermediate, parental‐like, transgressive and non‐significant. Contrary to the common view, parental‐like trait values were equally important in shaping the hybrid profile, as were intermediate traits. Transgression was revealed in 17/144 cases and concerned mainly cell and tissue sizes. As predicted by theory, we observed transgressive segregation more often when there was little phenotypic divergence, but intermediate values when parental traits were differentiated. It is likely that cell and tissue sizes are phylogenetically more conserved due to stabilising selection, whereas traits such as leaf thickness and volume fraction of the intercellular spaces, showing a consistent intermediate pattern across populations, are more susceptible to directional selection. Hybrid populations showed little similarity in expression profile, with only three traits identically expressed across all populations. Thus local adaptation of parental species and specific genetic background may be of importance.