Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Intraspecific perspective of phenotypic coordination of functional traits in Scots pine

  • Bárbara Carvalho ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    babic.barbosa@gmail.com

    Affiliations Departamento Biogeografía y Cambio Global, LINCGlobal, Museo Nacional de Ciencias Naturales, (MNCN-CSIC), Madrid, Spain, Área de Biodiversidad y Conservación, Universidad Rey Juan Carlos, Móstoles, Madrid, Spain

    ⨯
  • Cristina C. Bastias,

    Roles Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Departamento Biogeografía y Cambio Global, LINCGlobal, Museo Nacional de Ciencias Naturales, (MNCN-CSIC), Madrid, Spain

    ⨯
  • Adrián Escudero,

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation Área de Biodiversidad y Conservación, Universidad Rey Juan Carlos, Móstoles, Madrid, Spain

    ⨯
  • Fernando Valladares,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

    Affiliations Departamento Biogeografía y Cambio Global, LINCGlobal, Museo Nacional de Ciencias Naturales, (MNCN-CSIC), Madrid, Spain, Área de Biodiversidad y Conservación, Universidad Rey Juan Carlos, Móstoles, Madrid, Spain

    ⨯
  • Raquel Benavides

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Departamento Biogeografía y Cambio Global, LINCGlobal, Museo Nacional de Ciencias Naturales, (MNCN-CSIC), Madrid, Spain

    ⨯

Abstract

Functional traits have emerged as a key to understand species responses to environmental conditions. The concerted expression of multiple traits gives rise to the phenotype of each individual, which is the one interacting with the environment and evolving. However, patterns of trait covariation and how they vary in response to environmental conditions remain poorly understood, particularly at the intraspecific scale. Here, we have measured traits at different scales and in different organs, and analysed their covariation in a large number of conspecifics distributed in two contrasting environments. We expected significant correlations among traits, not only within clusters of traits as found in global, multispecies studies, but also among clusters, with more relationships within clusters, due to genetic constraints, and among clusters due to more coordinated phenotypes than community level, multispecies studies. We surveyed 100 Pinus sylvestris trees in a Mediterranean mountainous area distributed in two contrasting elevations. We measured 13 functional traits, in three clusters (leaf, stem and whole-plant traits), and analysed their variation and coordination. We found significant coordination among traits belonging to different clusters that reveals coordinated phenotypes. However, we found fewer correlations within trait clusters than initially expected. Trait correlation structures (number, intensity and type of correlations among traits) differed among individuals at different elevations. We observed more correlations within trait clusters at low elevation compared to those at high elevation. Moreover, the higher number of correlations among different trait clusters and the lower trait variation at the higher elevation suggests that variability decreases under more stressful conditions. Altogether, our results reveal that traits at intraspecific scale are coordinated in a broad network and not only within clusters of traits but also that this trait covariation is significantly affected by environmental conditions.

Citation: Carvalho B, Bastias CC, Escudero A, Valladares F, Benavides R (2020) Intraspecific perspective of phenotypic coordination of functional traits in Scots pine. PLoS ONE 15(2): e0228539. https://doi.org/10.1371/journal.pone.0228539

Editor: Harald Auge, Helmholtz Centre for Environmental Research - UFZ, GERMANY

Received: June 26, 2019; Accepted: January 17, 2020; Published: February 13, 2020

Copyright: © 2020 Carvalho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Data are available in figshare: (dx.doi.org/10.6084/m9.figshare.8326016).

Funding: BC thanks for a research scholarship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES [DOC-PLENO - Programa de Doutorado Pleno no Exterior. Grant Agreement No 99999.001266/2015-02]. All authors acknowledge support from the European Union Horizon 2020 Research and Innovation Programme Project GenTree [Grant Agreement No. 676876]; REMEDINAL TC-CM [Autonomous Community of Madrid, S2018/EMT-4338]; COMEDIAS FEDER [CGL2017-83170-R, Spanish Ministry of Science, Innovation and Universities] and LINCGlobal [4540-143AP] projects.

Competing interests: The authors have declared that no competing interests exist.

Introduction

One of the major challenges in plant ecology over the last decades has been the development of a general plant classification framework based on adaptive plant strategies [1,2]. The dominant species-based or taxonomic perspective of communities ecology have been gradually replaced by an approach based on functional traits, which can capture general adaptive features in a continuous framework [3]. Plant functional traits are morphological, phenological or physiological features able to confer competiveness and induce niche differences among coexisting plant species [4]. Therefore, they are considered indicators of the ecological role and realised niche of each species within plant communities [3,5,6]. Differences in functional traits among species are then associated with specific ecological strategies, competitiveness and niche breadth and shape whose study are indispensable to elucidate the mechanisms underlying the assembly of plant communities [7,8].

Relationships among functional traits in plant species define coordinated spectra related to common variations in functional strategies of plants. One of the most studied spectrum is the so-called leaf economics spectrum (LES), which predicts how leaf traits should vary across environmental gradients and co-vary among themselves. It describes a universal axis from “slow” to “fast” resource use strategies of plant species [9–13] confronting low-cost, short-lived leaves with rapid return of carbon and nutrients vs. costly long-lived leaves with slow returns. Similarly, the wood economics spectrum (WS) reflects the trade-offs among important wood functions like transport safety, transport efficiency, and mechanical support. Within WS, wood density works as an integrator of wood properties, related to the mechanical support and resistance to embolism in drought periods [14]. As a result of the interaction of the different traits belonging to distinct plant organs, different phenotypes are displayed at the individual level [15,16], which is likely to have a profound adaptive significance in contrasting environments [17,18].

Searching for comprehensive plant strategies, some works addressed trait covariation studies among these spectra or cluster of traits [19–23]. However, results are contrasting depending on the spatial scale of the study. For instance, studies working at wide environmental scales with multiple species showed that leaf and wood economics spectraare not coupled [20,21,24] claiming that trait covariation at wide environmental scales occurs mainly within spectra and independently among trait clusters (Fig 1a). However, Messier et al. [23] evidenced important coordination of traits from different trait clusters when the analysis is conducted at the community level, i.e. when coexisting species have overcome the same environmental filters and limitations (Fig 1b). These findings concur with many others describing relationships between WS and LES traits at the community and intraspecific scales, which showed that denser wood is positively related with more long-lasting leaves [22,25–28]. These studies suggest that traits from different clusters shape a comprehensive plant economics spectrum [29].

thumbnail
Download:
Fig 1. Conceptual figure with general patterns depicted from previous studies of trait correlations (either positive or negative, continuous vs. dashed lines) at different scales (from local to wide environmental spatial scale).

Relationships among traits at (a) large scale where trait clusters (PT: plant size traits, ST: stem traits and LT: leaf traits) of numerous species are mostly independent from each other (elaborated from Baraloto et al., 2010; Fortunel et al., 2012; Díaz et al., 2016). Relationships among traits at (b) the community, local scale with increased level of phenotypic coordination particularly among trait clusters (elaborated from Bucci et al., 2004; Santiago et al., 2004; Ishida et al., 2008; de la Riva et al., 2016; Messier et al., 2017b). Relationships among traits (c) at the intraspecific, local level: we expected an intermediate situation from a) and b) with more relationships within cluster than b) due to genetic constraints, and among clusters than a) due to more integrated phenotypes.

https://doi.org/10.1371/journal.pone.0228539.g001

Intraspecific trait variation is gaining increased recognition [8, 30] and the relationships among traits within species are also expected to change in response to different environmental conditions [31,32]. Plants respond to biotic and abiotic conditions and their phenotypic expressions are bound to coordination rules acting at different spatial scales [33,34]. The number of trait combinations and therefore the number of phenotypes are not unlimited; some are constrained by genetic or biophysical factors, while others are just unsuitable for a given environment and discarded by natural selection [35]. Coordinated changes among plant organs (and traits) can be the result of natural selection that hampers one organ to evolve independently of the other [16,36].

Environmental conditions within the geographical and altitudinal distributions of a plant species vary determining important changes in the phenotypic expressions among populations [37,38]. For instance, at high elevation, plants face changes in physical conditions like temperature decrease, frost, snow accumulation, strong winds or high irradiance [39]. Accordingly, several studies demonstrated trait adjustments aligned to a ‘economic’ strategy searching for resource conservation (persistence), like increasing leaf mass per area (LMA) and carbon isotope discrimination (δ13C), and lower leaf area and photosynthetic pigments due to low temperatures [40,41]. Moreover, nutrients uptake limitations in cold environments affect the growth rates at high elevation, likely decreasing the leaf nitrogen content and favouring a reduction in the specific leaf area [42,43]. Evidence also showed changes in stem properties with increasing elevation, thus wood density tends to be higher [28] and bark thickness higher protecting against frost [44] compared to individuals at low elevations.

Here, we measured traits in different plant systems (hereafter, trait clusters), i.e. stem, leaf and whole-plant traits, in individual tress of Pinus sylvestris (Scots pine) in a Mediterranean mountainous area in Spain at two different elevations. We firstly explored the overall trait network at the intraspecific level, analysing trait correlations within and among trait clusters. Secondly, we compared trait coordination in trees at the two very close sites but with contrasting elevations, to quantify the environmental impact on trait networks. We assumed that the environmental differences between both study sites are mainly driven by elevation. We expected to find a high number of correlations among traits within trait clusters, higher than those found in multispecies studies due to the genetic constraints within a species. We also expected a large number of correlations among the three trait clusters, similarly to multispecies studies at the community level due to the existence of important relationships among clusters also at the individual level (Fig 1c).

Material and methods

Study area

The study was performed at the Sierra de Guadarrama National Park (Madrid, Spain) during June and July 2016. Particularly, we chose two sites very close geographically (ca 7km) but with contrasting elevations, namely Pingarrón at 1900m (40° 48´50”N; 3° 58´12”W) and Ventorrillo Biological Station at 1440m (40° 45´31”N; 4° 00´49”W). Both sites represent natural, dominant and well-conserved stands of P. sylvestris with no signs of recent management and an uneven age structure. They represent the natural highest and lowest elevational distribution of the species in the area. In both sites, the climate is mountainous Mediterranean with wet, cold winters and warm dry summers. Annual mean temperature is 7.9°C in Ventorrillo and 6.5°C in Pingarrón and annual precipitation is 897 mm in Ventorrillo and 1242 mm in Pingarrón. Mean precipitation during the three driest months is 22 mm in Pingarrón and 24 mm in Ventorrillo (Worldclim; [45]). The bedrock in the area is mainly composed of granite and gneiss, and soils are acid and relatively homogeneous, predominantly humic cambisol soils with leptosol at higher-elevation sites [46]. Field work was carried out in the framework of an official research grant from the Autonomous Region of Madrid (REMEDINAL TE-CM (S2018/EMT-4338) and no specific authorization for the activities included in this paper was requested.

Pingarrón site is facing north and the understorey vegetation includes oromediterrenean shrublands, such as Juniperus communis L. subsp alpina (Neilr.) Čelak, Cytisus oromediterraneus (G. López & C.E. Jarvis) Rivas Mart. and Adenocarpus complicatus (L.) Link. Ventorrillo site has a south-west exposure and includes a layer of deciduous Pyrenean oak (Quercus pyrenaica Willd.) under the dominant pine canopy, and the understorey vegetation is mainly composed of Cistus laurifolius L., Cytisus scoparius L. and Genista cinerascens Lange.

Sampling design and functional trait measurements

We conducted an individual-based sampling, with 50 individuals sampled at Pingarrón (1900m) and 50 individuals at Ventorrillo (1440m). We randomly selected individuals among the mature and healthy ones and collected a fully sun-exposed branch from the top of the crown. We also extracted two 5-mm diameter wood cores at 50cm above the ground, using an increment borer.

In each tree, we measured 13 functional traits involved in relevant plant ecological functions. A summary of traits included in our data set and description of their ecological significance are described in Table 1 and related literature (S1 Table). We selected a wide variety of traits with different natures (morphological and chemical) classified in three trait clusters critical for the global plant spectrum: i) Plant size-related traits (whole-plant traits) that reflect the ability to compete for resources: plant height (m), crown depth (m) and diameter at breast height (DBH, cm); ii) Stem traits, related to transport and defence functions: bark thickness (mm) and trunk wood density (WD, g/cm3), and iii) Leaf traits, that balance construction costs vs. growth potential: leaf dry matter content (LDMC; mg/g), specific leaf area (SLA; mm2/mg), leaf nitrogen content (LNC; %), leaf carbon content (LCC; %), leaf carbon isotope (δ13C, ‰), chlorophyll a (μg/g), chlorophyll b (μg/g) and beta-carotene (μg/g). Detailed methods for trait measurements are provided in the Supporting Information (S1 Appendix). Before further analyses, we compared DBH, plant height and age between both populations to discard differences due to different tree size and ontogeny on traits (S2 Appendix). Note that we did not remove correlated traits, because we aimed at evaluating the phenotypic coordination (statistical integration [47]). We assumed that correlated traits exhibited either different patterns of correlation with other traits, or different patterns of change along elevation [48].

thumbnail
Download:
Table 1. Study traits and their functional significance.

https://doi.org/10.1371/journal.pone.0228539.t001

Statistical analyses

We first described the magnitude of intraspecific trait variability (ITV) for each trait by calculating the coefficient of variation (CV) and its 95% confidence limits using bootstrapping with replacement (with 500 replicates). We performed a t-test to analyse differences in trait values between sites. The traits (SLA, LNC, Chlorophyll a, Chlorophyll b, Beta-carotene, stem WD and DBH) were transformed with a Box-Cox transformation as appropriate, to meet normality assumption using ‘AID’ R package and boxcoxnc function.

We ran a correlation analysis to evaluate pairwise relationships among all traits, pooling data of individuals from both sites and separately. Secondly, we compared the three correlation matrices using a chi-square test. This test provides a χ2 value that represents the difference between a pair of correlation matrices, hypothesizing that they do not differ. Then, we also estimated how closely connected are our trait clusters from each other using network analyses. These analyses are quantitative approaches assessing the connectivity and distance of interconnected objects. The objects, our traits, are represented as nodes and their connectivity, their correlations, are represented as edges linking them. We constructed network graphs and carried out network analyses with the ‘igraph’ package using only the statistically significant trait correlations (p < .05). We used the modularity() function to measure the structure of the network for both populations. This metric calculates the fraction of edges within the defined clusters minus the expected fractions if the edges were random [49,50]. High modularity (close to Q = 1, the maximum) reflects dense connections within a trait cluster and weak connections among clusters. If the number of connections is close to zero (Q = 0) the trait network has low modularity.

To assess the influence of elevation on all traits and their relationships, we performed an ordination analysis. Trait data for ordination analysis were standardized using the R function scale. We initially used a Detrended Correspondence Analysis (DCA), to estimate the axes length in units of average standard deviation [51,52]. As the length of the first DCA axis was relatively short (Standard Deviation units; uSD = 0.13), we conducted a Redundancy Analysis (RDA) [53], which assumes linear relationships between elevation and the elements of the trait matrix [54,55]. In our case, the trait matrix was constrained by elevation. Total variation explained (TVE) of data set was the value of canonical extracted axes (Σcons) using the constraining data matrix [56]. Finally, we performed a Monte Carlo permutation test (1000 randomizations) to determine the accuracy of the relationship between the elevation and the trait data sets.

All statistical analyses were conducted in R version 3.5.1 [57]. Comparison of trait correlation matrices between sites was performed by using the cortest.mat function in the 'psych' package [58]. Multivariate analysis was conducted using the 'vegan' package [59].

Results

Traits and their correlations

The magnitude of variation at the intraspecific level depended on the study trait (Tables 1 and 2). Both morphological (LDMC, SLA) and chemical leaf traits (%LNC, %LCC and δ13C) and WD exhibited relatively low intraspecific variability (ITV) with coefficients of variation (CV) below 15% (Table 2). In contrast, plant size traits (plant height, crown depth and DBH); bark thickness, photosynthetic pigments (chlorophyll a, b and beta-carotene) showed substantial variability with CV between 29% and 50% (Table 2).

thumbnail
Download:
Table 2. Descriptive statistics of the study traits measured on Pinus sylvestris individual trees.

https://doi.org/10.1371/journal.pone.0228539.t002

The correlation structure with pooled data had an edge density of 32% of the total pairwise correlations (Fig 2a; S2a Table). We found positive significant correlations within trait clusters, i.e. among the three architectural traits (plant height, DBH and crown depth) and between the two stem traits (WD and bark thickness). Regarding leaf cluster, the results showed some correlations describing the trade-off between acquisition (SLA, %LCN, Chlorophyll a) and conservation (LDMC, %LCC, δ13C) of resources (Fig 2a; S2a Table). The three photosynthetic pigments were positively correlated among them, without any other connection with other leaf trait.

thumbnail
Download:
Fig 2. Correlograms of traits measured in individual Scots pines calculated with pooled data and data from both study sites separately (Pingarrón at 1900m and Ventorrillo at 1440m).

Solid lines represent positive correlations, and dashed lines negative correlations. Line thicknesses are proportional to the correlation strength. Only significant correlations are shown. Plant size traits: plant height, crown depth and DBH: Stem traits: WD and bark thickness; Leaf traits: LDMC, SLA, δ13C, %LCC, %LNC, Chlorophyll a. Chlorophyll b and Beta-carotene. Different colours indicate different trait clusters and their different functions (yellow: stem traits, primarily mechanical support; blue: plant size traits, primarily competitive ability; green: leaf traits, primarily resource acquisition and conservation).

https://doi.org/10.1371/journal.pone.0228539.g002

We found also correlations among clusters. Plant size traits and stem traits were positively correlated (Fig 2a; S2a Table). Among leaf traits, LDMC was positively correlated to plant height, DBH and bark thickness; leaf carbon content (%LCC) increased with increasing DBH and bark thickness, and SLA correlated negatively with WD, crown and plant size (crown depth and plant height) (Fig 2a; S2a Table). Modularity test for pooled data showed that the study traits presented more connections within clusters than across them (Q = 0.10).

Elevation effect on traits and their correlations

The results from RDA ordination showed that elevation accounted for a low and significant fraction of the total variation explained (5%) when all traits were simultaneously considered (S2 Fig; Σcons = 0.64; TVE = 0.049; Monte Carlo test: F1,87 = 4.56; p = 0.001). The first RDA axis represents the traits which differed between the two sites, i.e. chlorophyll a, LDMC and SLA. Despite the small fraction of the total variance explained by altitudinal condition, some leaf trait values were significantly different between both sites (SLA, LDMC and photosynthetic pigments) and bark thickness (Table 2). Particularly, individuals at high elevation presented a significantly greater bark thickness and higher SLA; while individuals at low elevations showed a significantly higher LDMC and photosynthetic pigments content (Table 2; S1 Fig).

The correlation analysis independently run for each site showed a similar number of significant trait correlations, 21 for individuals at high-elevation (26% of the total pairwise correlations tested) and 22 for individuals at low-elevation (28%). We also found a greater number of trade-offs, significant negative correlations, among different clusters at high elevation (13) compared with low elevation (8) (Fig 2b and 2c; S2b Table). Results from the modularity test reflected that traits at high elevation (Q = 0.0034) have denser connections among clusters and sparse connections within clusters than traits at low elevation (Q = 0.23).

Correlations between plant size and stem traits were strong and similar for individuals at different elevations. Results also showed similar patterns with leaf traits for both sites, namely a positive correlation among leaf photosynthetic pigments, and a negative correlation between SLA and LDMC and between SLA and δ13C (Fig 2b and 2c; S2b Table). However, we found some different outcomes between both elevations, resulting in different correlation matrices (χ2 = 238.38, df = 156, p < 0.00001). For instance, at high elevation, we found that SLA was negatively (and LDMC positively) correlated with bark thickness, crown depth and plant height. Moreover, chlorophyll b was positively related to crown depth and beta-carotene showed a positive correlation with plant height (Fig 2b; S2b Table). Among all significant correlations found at low elevation, SLA was positively related with chlorophyll a, and LDMC was positively related with %LCC and plant height. Photosynthetic pigments, chlorophyll a and b, were also correlated negatively with δ13C and positively with bark thickness. Bark thickness, in turn, was positively related with DBH and %LCC. (Fig 2c; S2b Table).

Discussion

Our results showed that the three trait clusters (whole-plant, stem and leaf traits) are highly coordinated among them (Fig 2), reflecting a plant economics spectrum [29]. Although we found fewer and weaker covariations than we had expected, it is clear that traits are connected in a broad network. These results evidence that the existence of trade-offs organizing the diversity of possible plant phenotypes into individual trait clusters found across species and at large spatial scales, it is not maintained at the intraspecific level, at least for the study species.

Studies at the community level showed that trade-offs proposed by LES between acquisition and resources conservation are not always found [12,60–62], reflecting different sensitivity of LES traits to different scales depending on the main environmental, genetic and biophysical drivers affecting them [12]. Natural selection shapes phenotypes with correlated traits based on the interaction of genes and the environmental conditions. It means that certain trait combinations may be favoured or discarded in given environments [63,64]. At the intraspecific scale, evidence analysing this pattern is scarce and highly species-specific [6,65–68]. For example, it is remarkable the lack of relationships in our data between SLA and other leaf traits, like N content, to which has been frequently associated [69]. This lack of correlations can be explained by the coniferous nature of our study species, following the outcomes described by Lusk et al. [70], who showed no relationship between leaf N and SLA in other conifers, in a comprehensive study including many angiosperms and conifers. This result supports the idea that LDMC is a more reliable leaf trait to reflect the leaf trait cluster as suggested by Wilson et al. [71], especially for conifers, than SLA, which depends greatly on the spatial dimensions of leaves that are highly variable [72]. Similarly, we found remarkable the lack of correlation among photosynthetic pigments and other leaf traits, like N content, a pattern described in literature [73–75]. Looking at the variability of leaf traits (Table 2), we suggest that it might be associated with their low variability (all below 15%, excluding the pigments) compared to the other studied traits and findings from other studies addressed at the intraspecific level [30] that encompassed a wider environmental gradient.

The plant size traits were tightly correlated, being taller those trees with thicker trunks and larger crowns. This plant size cluster is complemented with stem traits, as bigger trees had denser wood and thicker bark. Wood density is considered a good indicator of mechanical resistance and hydraulic safety [76–79]. Namely, it is associated to the resistance to embolism, because the greater is the wood density, the narrower are the conduits that confer embolism resistance [80–82]. Cold and drought imperil plants to embolism, and plants in dry or cold areas usually have narrower water-conducting conduits compared to plants thriving in moist, warm areas. Thus, narrow conduits have been considered a key adaptation to cold and drought conditions, together with short plant height [82]. In the line with our results, other studies have found that bark thickness also correlates positively with tree size, bestowing trunk mechanical strength and protection against frost and drought [44,83–86]. Similarly, previous studies have shown the relationship between bark thickness and DBH [85,86] with increasing bark accumulation as trees grow.

Moreover, our results showed interactions among leaf traits and other trait clusters within a species, similarly to studies analysing simultaneously trait clusters in different species [19,21–23]. For instance, increasing LDMC and leaf C content appeared in bigger trees. Greater SLA was also connected with smaller trees with lower WD. Leaf traits were measured in needles collected from the top part of the crown that did not necessary meant fully exposed needles in dominated trees. Therefore, this correlation would reflect the effect of light availability (greater in taller trees) on leaf characteristics [87]. Altogether, this outcome supports our first prediction that all traits measured on leaves, stem and plant size co-vary shaping a broad network of coordinated trait clusters. In other words, they endorsed the existence of a whole-plant economics spectrum.

Our study showed a different pattern in the trait coordination in the two close sites in contrasting elevation, reflecting different sensitivity of traits and their covariation to environmental conditions. An elevation difference of 450 m entailed, amongst others, different rainfall and temperature values (annual mean temperature of 6.5°C vs 7.9 °C and mean annual precipitation of 1242 mm vs. 897 mm), which triggered significant differences in functional traits among pines from both sites. Individuals growing at higher elevation presented greater bark thickness than those at lower elevation (S1 Fig), potential response to harsher conditions in terms of snow, wind at the top of the mountain [84]. Additionally, individuals growing at the lower elevation had leaves with lower SLA and greater LDMC (S1 Fig). These leaf characteristics are frequently found in dry sites [88], indicating a conservative strategy to extend the life-span of expensive organs (needles) in poor or adverse environments. Nevertheless, low elevation trees had also a higher amount of photosynthetic pigments, indicative of higher productivity. The explanation of this apparent inconsistency may lie in the most challenging conditions that plants experience at high elevations, as a combination of high irradiance and low temperatures [39]. Coldness at high elevations inhibits enzymatic reactions reducing carbon absorption, without affecting light capture, triggering a protection mechanism, photoprotection, which eventually implies a reduction of the amount of chlorophylls [89,90].

Two possible factors can underlie these changes in trait variation between both sites (S1 Fig): phenotypic plasticity and local genetic variation. These two factors are not mutually exclusive and they probably act together to help plants to more effectively deal with the different environmental conditions [91]; however, with our data, we cannot discern them. At the site level, we observed that elevation, as a proxy of other environmental variables such temperature, vegetation and stoniness, accounted for a small but significant 5% of total trait variation. Indeed, trait variation (in CV) was greater at low elevation (Table 2) which suggests that trait variability decreases under more stressful conditions. An explanation would be a greater amount of available resources in mild sites that would provide more opportunities (i.e. increase the available niche space) to individuals, that could thrive in different micro-sites [5,7,38]. On the other hand, at high elevation we found a greater number of correlations between clusters revealing that pines up there displayed a more coordinated phenotype compared with those at low elevation. This entails that elevation, as a proxy of others environmental factors, played an important role in the trait coordination, detectable in close individuals of Scots pine. Moreover, despite gene flow is expected between both stands due to the anemophilous pollination of pines and the proximity between them (ca 7km), we acknowledge that genetic differences may partly explain the remaining unexplained variation for phenotypic differences between both sites. An integrated phenotype means that any change in one trait implies changes in the whole phenotypic response of an individual [36]. This scaling relationship can be the result of natural selection acting over individuals at the treeline [92], favouring the coordination and adjustment of the whole phenotype to the environmental conditions at such elevation. Our results agree with those hypothesised by Keddy [93] and Violle et al. [4] which posed that environmental factors act as filters at the individual level. On the contrary, individuals at low elevation did not show the same trait coordination, evidencing that changes in a given trait would not necessarily unleash shifts in traits from a different cluster.

This work represents one of the first attempts to analyse phenotypic coordination within a species under contrasting environmental conditions, in a comprehensive way with a large array of traits. The significant correlations among three trait clusters related to important dimensions of the global plant spectrum, revealed the existence of a more complex and coordinated phenotypes than those found in single-spectrum or one trait cluster approaches. In addition, our results showed differences in the trait correlation networks induced by contrasting environmental conditions that provide evidence for different sensitivity of traits and their covariation to the external conditions.

Supporting information

S1 Appendix. Detailed methods.

https://doi.org/10.1371/journal.pone.0228539.s001

(DOC)

S2 Appendix. Descriptive statistics for the plant size traits measured on Pinus sylvestris individual trees.

Confidence limits surrounding coefficients of variation (CV) were calculated through bootstrapping with replacement (with n = 500 replicates). Comparison of traits values: a) plant height; b) diameter at breast height (DBH) and c) Age at the two different elevations using a t-test are also shown.

https://doi.org/10.1371/journal.pone.0228539.s002

(DOCX)

S1 Table. Study traits, the ecological function they represent and literature describing them.

https://doi.org/10.1371/journal.pone.0228539.s003

(DOCX)

S2 Table. Trait correlations between traits measured in a) pooled data of Pinus sylvestris and b) distinguishing the two study sites (under the diagonal trees from Pingarrón at 1900 m, and above the diagonal trees from Ventorrillo at 1440 m).

It was used the Pearson´s rank correlation analyses. ‘ns’ means ‘not significant’. Brackets indicate the number of individuals at each correlation; number at left side is for Pingarrón correlations and at right side for Ventorrillo correlations.

https://doi.org/10.1371/journal.pone.0228539.s004

(DOCX)

S1 Fig. Bean plot of the traits that were significantly different between both elevations using a t-test (Pingarrón 1900m, Ventorrillo 1440 m).

(SLA: specific leaf area; LDMC: leaf dry matter content).

https://doi.org/10.1371/journal.pone.0228539.s005

(DOC)

S2 Fig. Ordination diagram of traits (RDA) in relation to the elevation of both surveyed sites.

The orientation of line directions indicates the sign of the correlation among the traits, and the length is related to the strength. Triangles indicate the centroid of both populations Pingarrón in pink (1900m) and Ventorrillo in green (1440m). (DBH: diameter at breast height; Stem WD: stem wood density; LDMC: leaf dry matter content; SLA: specific leaf area; LNC%: leaf nitrogen content, LCC%: leaf carbon content; δ13C leaf carbon isotope discrimination).

https://doi.org/10.1371/journal.pone.0228539.s006

(DOC)

Acknowledgments

We are very grateful to David Quiroga and Patricia Alonso for the invaluable technical help. Also we thank to Margot Kaye and Sennaquerib Rodriguez for their help in the laboratory and field work.

References

  1. 1. Westoby M. A leaf-height-seed (LHS) plant ecology strategy scheme. Plant and Soil. 1998; 199: 213–227.
  2. 2. Díaz S, Kattge J, Cornelissen JHC, et al. The global spectrum of plant form and function. Nature. 2016; 529: 167–171. pmid:26700811
  3. 3. McGill BJ, Enquist BJ, Weiher E, Westoby M. Rebuilding community ecology from functional traits. Trends Ecol Evol. 2006; 21:178–185. pmid:16701083
  4. 4. Violle C, Navas ML, Vile D, et al. Let the concept of trait be functional! Oikos. 2007; 116: 882–892.
  5. 5. Violle C, Jiang L. Towards a trait-based quantification of species niche. J Plant Ecol. 2009; 2: 87–93.
  6. 6. Vasseur F, Violle C, Enquist BJ, Granier C, Vile D. A common genetic basis to the origin of the leaf economics spectrum and metabolic scaling allometry. Ecol Lett. 2012; 15: 1149–1157. pmid:22856883
  7. 7. Sides CB, Enquist BJ, Ebersole JJ, Smith MN, Henderson AN, Sloat LL. Revisiting darwins hypothesis: Does greater intraspecific variability increase species ecological breadth? Am J Bot. 2014; 101: 56–62. pmid:24343815
  8. 8. Escudero A, Valladares F. Trait-based plant ecology: moving towards a unifying species coexistence theory: Features of the Special Section. Oecologia. 2016; 180: 919–922. pmid:26897604
  9. 9. Reich PB, Walters MB, Ellsworth DS. From tropics to tundra: Global convergence in plant functioning. Proc Natl Acad Sci U.S.A. 1997; 94: 13730–13734. pmid:9391094
  10. 10. Reich PB, Ellsworth DS, Walters MB, et al. Generality of leaf trait relationships: a test across six biomes. Ecology. 1999; 80: 1955–1969.
  11. 11. Wright IJ, Reich PB, Westoby M, et al. The worldwide leaf economics spectrum. Nature. 2004; 428: 821–827. pmid:15103368
  12. 12. Messier J, McGill BJ, Enquist BJ, Lechowicz MJ. Trait variation and integration across scales: is the leaf economic spectrum present at local scales? Ecography. 2017a; 40: 685–697.
  13. 13. Onoda Y, Wright IJ, Evans JR, et al. Physiological and structural tradeoffs underlying the leaf economics spectrum. New Phytol. 2017; 214: 1447–1463. pmid:28295374
  14. 14. Chave J, Coomes D, Jansen S, Lewis SL, Swenson NG, Zanne AE. Towards a worldwide wood economics spectrum. Ecol Lett. 2009; 12: 351–366. pmid:19243406
  15. 15. Murren C. The integrated phenotype. Integr Comp Biol. 2002; 52: 64–76.
  16. 16. Pigliucci M. Phenotypic integration: studying the ecology and evolution of complex phenotypes. Ecol Lett. 2003; 6: 265–272.
  17. 17. Auger S, Shipley B. Inter-specific and intra-specific trait variation along short environmental gradients in an old-growth temperate forest. J Veg Sci. 2012; 24: 419–428.
  18. 18. Umaña MN, Swenson NG. Does trait variation within broadly distributed species mirror patterns across species? A case study in Puerto Rico. Ecology. 2019 Apr 29:e02745. pmid:31032887
  19. 19. Freschet T, Cornelissen JHC, Van Logtestijn RSP. Evidence of the ‘plant economics spectrum’ in a subarctic flora. J Ecol. 2010; 98: 362–373.
  20. 20. Baraloto C, Paine CET, Poorter L, et al. Decoupled leaf and stem economics in rain forest trees. Ecol Lett. 2010; 13: 1338–1347. pmid:20807232
  21. 21. Fortunel C, Fine PVA, Baraloto C, Ecologie UMR, Cedex K. Leaf, stem and root tissue strategies across 758 Neotropical tree species. Funct Ecol. 2012; 26: 1153–1161.
  22. 22. de La Riva EG, Tosto A, Pérez-Ramos IM, et al. A plant economics spectrum in Mediterranean forests along environmental gradients: is there coordination among leaf, stem and root traits? J Veg Sci. 2016; 27: 187–199.
  23. 23. Messier J, Lechowicz MJ, McGill BJ, et al. Interspecific integration of trait dimensions at local scales: the plant phenotype as an integrated network. J Ecol. 2017b; 105: 1775–1790.
  24. 24. Wright IJ, Ackerly DD, Bongers D, Harms KE, et al. Relationships among ecologically important dimensions of plant trait variation in seven neotropical forests. Ann Bot. 2007; 99:1003–1015. pmid:16595553
  25. 25. Bucci SJ, Goldstein G, Meinzer FC, et al. Functional convergence in hydraulic architecture and water relations of tropical savanna trees: from leaf to whole plant. Tree Physiol. 2004; 24: 891–899. pmid:15172839
  26. 26. Santiago LS, Goldstein G, Meinzer FC, et al. Leaf photosynthetic traits scale with hydraulic conductivity and wood density in Panamanian forest canopy trees. Oecologia. 2004; 140: 543–550. pmid:15232729
  27. 27. Ishida A, Nakano T, Yazaki K, Matsuki S, Koike N, Lauenstein DL, et al. Coordination between leaf and stem traits related to leaf carbon gain and hydraulics across 32 drought-tolerant angiosperms. Oecologia. 2008; 156: 193–202. pmid:18297313
  28. 28. Fajardo A, Piper FI. Intraspecific trait variation and covariation in a widespread tree species (Nothofagus pumilio) in southern Chile. New Phytol. 2011; 189: 259–271. pmid:21039558
  29. 29. Reich PB. The world-wide “fast-slow” plant economics spectrum: a traits manifesto. J Ecol. 2014; 102: 275–301.
  30. 30. Bastias CC, Fortunel C, Valladares F, Baraloto C, Benavides R, et al. Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient. PLoS ONE. 2017; 12(2):e0172495. https://doi.org/10.1371/journal.pone.0172495. pmid:28241033
  31. 31. Gianoli E, Palacio-López K. Phenotypic integration may constrain phenotypic plasticity in plants. Oikos. 2009; 118: 1924–1928.
  32. 32. Salgado-Negret B, Canessa R, Valladares F, Armesto JJ, Pérez F. Functional traits variation explains the distribution of Aextoxicon punctatum (Aextoxicaceae) in pronounced moisture gradients within fog-dependent forest fragments. Front Plant Sci. 2015; 6: 1–11.
  33. 33. Premoli AC, Brewer CA. Environmental v. genetically driven variation in ecophysiological traits of Nothofagus pumilio from contrasting elevations. Aust J Bot. 2007; 55: 585–591.
  34. 34. Messier J, McGill BJ, Lechowicz MJ. How do traits vary across ecological scales? A case for trait-based ecology. Ecol Lett. 2010; 13: 838–848. pmid:20482582
  35. 35. Walsh B, Blows MW. Abundant genetic variation plus strong selection multivariate genetic constraints: a geometric view of adaptation. Annu Rev Ecol Evol Syst. 2009; 40: 41–59.
  36. 36. Chapin S III, Autumm K, Pugnaire F. Evolution of suites of traits in response to environmental stress. Am Nat. 1993; 142 (s1): S78–S92.
  37. 37. de Bello F, Lavorel S, Albert CH, et al. Quantifying the relevance of intraspecific trait variability for functional diversity. Methods Ecol Evol. 2011; 2: 163–174.
  38. 38. Violle C, Enquist BJ, McGill BJ, et al. The return of the variance: Intraspecific variability in community ecology. Trends Ecol Evol. 2012; 27: 244–252. pmid:22244797
  39. 39. Körner C. The use of “altitude” in ecological research. Trends Ecol Evol. 2007; 22: 569–574. pmid:17988759
  40. 40. Hultine R, Marshall D. Altitude trends in conifer leaf morphology and stable carbon isotope composition. Oecologia. 2000; 123: 32–40. pmid:28308741
  41. 41. Qiang WY, Wang XL, Chen T, et al. Variations of stomatal density and carbon isotope values of Picea crassifolia at different altitudes in the Qilian Mountains. Trees. 2003; 17: 258–262.
  42. 42. Birmann K, Körner C. Nitrogen status of conifer needles at the alpine treeline. Plant Ecol Divers. 2009; 2: 233–241.
  43. 43. Fajardo A, McIntire EJB. Reversal of multicentury tree growth improvements and loss of synchrony at mountain tree lines point to changes in key drivers. J Ecol. 2012; 100: 782–794.
  44. 44. Arco Molina JG, Hadad MA, Patón Domínguez D, Roig FA. Tree age and bark thickness as traits linked to frost ring probability on Araucaria araucana trees in northern Patagonia. Dendrochronologia. 2016; 37: 116–125.
  45. 45. Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017; 37: 4302–4315.
  46. 46. Franco-Múgica F, García Antón M, Sainz Ollero H. Vegetation dynamics and human impact in the Sierra de Guadarrama, Central System, Spain. Holocene 1998; 8: 69–82.
  47. 47. Armbruster WS, Pelabon C, Bolstad GH, Hansen TF. Integrated phenotypes: understanding trait covariation in plants and animals. Philos Trans Royal Soc B. 2014; 369: 20130245.
  48. 48. Nicotra AB, Chazdon RL, Schlichting CD. Patterns of genotypic variation and phenotypic plasticity of light response in two tropical Piper (Piperaceae) species. Am J Bot. 1997; 84: 1542–1552. pmid:21708557
  49. 49. Clauset A, Newman MEJ, Moore C. Finding community structure in very large networks. Phys Rev E. 2004; 70: 066111: pmid:15697438
  50. 50. Newman M, Girvan M. Finding and evaluating community structure in networks. Phys Rev E. 2004; 69: 026113. pmid:14995526
  51. 51. Hill MO, Gauch HG. Detrended correspondence analysis, an improved ordination technique. Vegetatio. 1980; 42: 47–58.
  52. 52. Gauch HG. Noise reduction by eigenvector ordinations. Ecology. 1982; 63: 1643–1649.
  53. 53. Legendre P, Anderson MJ. Distance-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecol Monogr. 1999; 69: 1–24.
  54. 54. Boccard D, Legendre P, Drapeau P. Partialling out the spatial component of ecological variation. Ecology. 1992; 73: 1045–1055.
  55. 55. Peres-Neto PR, Legendre P, Dray S, Borcard D. Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology. 2006; 87: 2614–2625. pmid:17089669
  56. 56. Økland RH. On the variation explained by ordination and constrained ordination axes. J Veg Sci. 1999; 10: 131–136.
  57. 57. R Core Team. R: A language and environment for statistical computing. Vienna. R Foundation for Statistical Computing. 2018. http://www.R-project.org/.
    • 58. Revelle WL. psych: Procedures for personality and psychological research. R Package version 1.6.12. 2017. https://personality-project.org/r/psych/.
      • 59. Oksanen J, Blanchet FG, Friendly M, et al. Vegan: Community Ecology Package. R package version 2.4.2. 2017. https://github.com/vegandevs/vegan.
        • 60. Edwards EJ, Chatelet DS, Sack L, Donoghue MJ. Leaf life span and the leaf economic spectrum in the context of whole plant architecture. J Ecol. 2014; 102: 328–336.
        • 61. Mason CM, Goolsby EW, Humphreys DP, Donovan LA. Phylogenetic structural equation modelling reveals no need for an “origin” of the leaf economics spectrum. Ecol. Lett. 2016; 19: 54–61. pmid:26563777
        • 62. Chauvin KM, Asner GP, Martin RE, Kress WJ, Wright SJ, Field CB. Decoupled dimensions of leaf economic and anti-herbivore defense strategies in a tropical canopy tree community. Oecologia. 2018; 186: 765–782. pmid:29302802
        • 63. Berg RL. The ecological significance of correlation pleiades. Evolution. 1960; 14: 171–180.
        • 64. Cheverud J. Quantitative genetics and developmental constraints on evolution by selection. J Theor Biol. 1984; 110:155–171. pmid:6492829
        • 65. Blonder B, Vasseur F, Violle C, Shipley B, Enquist BJ, Vile D. Testing models for the leaf economics spectrum with leaf and whole-plant traits in Arabidopsis thaliana. AoB Plants. 2015, plv049, https://doi.org/10.1093/aobpla/plv049.
        • 66. Niinemets Ü. Is there a species spectrum within the worldwide leaf economics spectrum? Major variations in leaf functional traits in the Mediterranean sclerophyll Quercus ilex. New Phytol. 2015; 205: 79–96. pmid:25580487
        • 67. Martin AR, Rapidel B, Roupsard O, Van den Meersche K, de Melo Virginio Filho E, et al. Intraspecific trait variation across multiple scales: the leaf economics spectrum in coffee. Funct Ecol. 2017; 31: 604–612.
        • 68. Fajardo A, Siefert A. Intraspecific trait variation and the leaf economics spectrum across resource gradients and levels of organization. Ecology. 2018; 99: 1024–1030. pmid:29603183
        • 69. Evans RD. Physiological mechanisms influencing plant nitrogen isotope composition. Trends Plant Sci. 2001; 6: 121–126. pmid:11239611
        • 70. Lusk CH, Wright I, Reich PB. Photosynthetic differences contribute to competitive advantage of evergreen angiosperm trees over evergreen conifers in productive habitats. New Phytol. 2003; 160: 329–336.
        • 71. Wilson PJ, Thompson K, Hodgson JG. Specific leaf area and leaf dry matter content as alternative predictors of plant strategies. New Phytol. 1999; 143: 155–162.
        • 72. Pérez-Harguindeguy N, Diaz S, Garnier E, et al. New Handbook for standardized measurment of plant functional traits worldwide. Aust J Bot. 2013; 61:167–234.
        • 73. Reich P, Walters M, Ellsworth D. Leaf life-span in relation to leaf, plant, and stand characteristics among diverse ecosystems. Ecol Monogr. 1992; 62: 365–392.
        • 74. Lei TT, Tabuchi R, Kitao M, Koike T. Functional relationship between chlorophyll content and leaf reflectance, and light-capturing efficiency of Japanese forest species. Physiol Plant. 1996; 96: 411–418.
        • 75. Niinemets Ü, Portsmuth A, Tena D, et al. Do we underestimate the importance of leaf size in plant economics? Disproportional scaling of support costs within the spectrum of leaf physiognomy. Ann Bot. 2007; 100: 283–303. pmid:17586597
        • 76. Chave J, Muller-Landau HC, Baker TR, Easdale TA, ter Steege H, Webb OC. Regional and phylogenetic variation of wood density across 2456 neotropical tree species. Ecol Appl. 2006; 16: 2356–2367. pmid:17205910
        • 77. Sperry JS, Meinzer FC, McCulloh KA. Safety and efficiency conflicts in hydraulic architecture: Scaling from tissues to trees. Plant Cell Environ. 2008; 31: 632–645. pmid:18088335
        • 78. Martínez-Vilalta J, Cochard H, Mencuccini M, et al. Hydraulic adjustment of Scots pine across Europe. New Phytol. 2009; 184: 353–364. pmid:19674333
        • 79. Olson ME, Soriano D, Rosell JA, Anfodillo T, Donoghue MJ, Edwards EJ. et al. Plant height and hydraulic vulnerability to drought and cold. Proc Natl Acad Sci U.S.A. 2018; 115: 7551–7556. pmid:29967148
        • 80. Bennett AC, Mcdowell NG, Allen CD, Anderson-Teixeira KJ. Larger trees suffer most during drought in forests worldwide. Nat Plants. 2015; 1: 1–5.
        • 81. Greenwood S, Ruiz-Benito P, Martínez-Vilalta J, Lloret F, Kitzberger T, Allen CD, et al. Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area. Ecol Lett. 2017; 20: 539–553. pmid:28220612
        • 82. Fajardo A, McIntire EJB, Olson ME. When short stature is an asset in trees. Trends Ecol Evol. 2019; 34: 193–199. pmid:30447938
        • 83. Vines RG. Heat transfer through bark, and the resistance of trees to fire. Aust J Bot. 1968; 16: 499–514.
        • 84. Niklas KJ. The mechanical role of bark. Am J Bot 1999; 86: 465–469. pmid:10205065
        • 85. Paine CET, Stahl C, Courtois EA, Patiño S, Sarmiento C, Baraloto C. Functional explanations for variation in bark thickness in tropical rain forest trees. Funct Ecol. 2010; 24: 1202–1210.
        • 86. Brando PM, Nepstad DC, Balch JK, et al. Fire-induced tree mortality in a neotropical forest: The roles of bark traits, tree size, wood density and fire behavior. Glob Change Biol. 2012; 18:630–641.
        • 87. Niinemets Ü, Kull O. Effects of light availability and tree size on the architecture of assimilative surface in the canopy of Picea abies—variation in needle morphology. Tree Physiol. 1995; 15: 307–315. pmid:14965954
        • 88. Toivonen JM, Horna V, Kessler M, Ruokolainen K, Hertel D. Interspecific variation in functional traits in relation to species climatic niche optima in Andean Polylepis (Rosaceae) tree species: evidence for climatic adaptations. Funct Plant Biol. 2014; 41: 301–312.
        • 89. Zarter CR, Demmig-Adams B, Ebbert V, Adamska I, Adams WW III. Photosynthetic capacity and light harvesting efficiency during the winter-to-spring transition in subalpine conifers. New Phytol. 2006; 172: 283–292. pmid:16995916
        • 90. Verhoeven A. Sustained energy dissipation in winter evergreens. New Phytol. 2014; 201: 57–65.
        • 91. Valladares F, Matesanz S, Guilhaumon F, et al. The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecol. Lett. 2014; 17: 1351–1364. pmid:25205436
        • 92. Premoli AC. Isozyme polymorphisms provide evidence of clinal variation with elevation in Nothofagus pumilio. J Hered. 2003; 94: 218–226. pmid:12816962
        • 93. Keddy PA. Assembly and response rules: two goals for predictive community ecology. J Veg Sci. 1992; 3: 157–164.