Front. Plant Sci., 19 July 2017
Sec. Functional Plant Ecology

Biogeochemical and Ecomorphological Niche Segregation of Mediterranean Woody Species along a Local Gradient

  • 1Estación Biológica de Doñana, Consejo Superior de Investigaciones Científicas, Seville, Spain
  • 2Área de Ecología, Facultad de Ciencias, Universidad de Córdoba, Córdoba, Spain
  • 3Instituto de Recursos Naturales y Agrobiología de Sevilla, Consejo Superior de Investigaciones Científicas, Seville, Spain
  • 4CEFE UMR 5175, CNRS – Université de Montpellier – Université Paul-Valéry Montpellier – EPHE, Montpellier, France

According with niche theory the species are specialized in different ecological niches, being able to coexist as result of a differential use of resources. In this context, the biogeochemical niche hypothesis proposes that species have an optimal elemental composition which results from the link between the chemical and morphological traits for the optimum plant functioning. Thus, and attending to the limiting similarity concept, different elemental composition and plant structure among co-occurring species may reduce competition, promoting different functional niches. Different functional habits associated with leaf life-span or growth forms are associated with different strategies for resource uptake, which could promote niche partitioning. In the present study, based on the biogeochemical niche concept and the use of resources in different proportions, we have focused on leaf traits (morphological and chemical) associated with resource uptake, and explored the niche partitioning among functional habits: leaf life-span (deciduous, evergreen, and semideciduous) and growth (tree, shrub, and arborescent-shrub). To this end, we have quantified the hypervolume of the leaf functional trait space (both structure and chemical composition) in a sample of 45 Mediterranean woody species from Sierra Morena Mountains (Spain) growing along a local soil resource gradient. Our results show consistent variation in functional space for woody communities distributed along the environmental gradient. Thus, communities dominated by deciduous trees with faster growth and a predominant acquisitive strategy were characteristic of bottom forests and showed highest leaf biogeochemical space. While semideciduous shrubs and evergreen (arborescent, trees) species, characterized by a conservative strategy, dominated ridge forests and showed smaller functional space. In addition, within each topographical zone or environment type, the foliar biogeochemical niche partitioning would underlie the species ability to coexist by diverging on leaf nutrient composition and resource uptake. Lower niche overlap among functional habits were found, which support that different growth forms and leaf life-habits may facilitate the coexistence of the woody species and niche partitioning along and within the gradient.


According to niche theory, species coexistence is promoted by ecological niche differences (MacArthur and Levins, 1967). While this idea is pivotal in community ecology, empirical evidences are contrasted for several reasons. First, theoretically the niche concept is a long and still-debated concept. Second, operationally, quantifying the whole niche of a species remains elusive. As an alternative, it has been proposed to replace niche axes by phenotypic trait axes in community ecology (McGill’s et al., 2006; Violle and Jiang, 2009). In plants, a myriad of studies has focused on plant morphological traits to analyze species niche partitioning along environmental gradients as well as niche partitioning within communities (e.g., Mason et al., 2011; Maire et al., 2012; Lamanna et al., 2014). However, the interpretation of these studies is often ambiguous given that a given morphological trait can be involved in several, sometimes opposed, assembly processes (de la Riva et al., 2017). Community species composition can be viewed as the result of a selection process named filtering, which combines species interactions (biotic filtering) and environmental adversity (abiotic filtering) operating simultaneously and favoring a trait convergence (McGill’s et al., 2006). At the same time, the process of niche partitioning tends to favor coexistence of species with divergent traits and complementarity in the use of resources (Maire et al., 2012). Therefore, stabilizing niche differences and relative fitness differences are determined by abiotic and biotic environment (HilleRisLambers et al., 2012). In this regard, given that plants all use the same types of resources (light, nutrients, CO2, water), the biochemical and ecomorphological characterization of plant niches may be a way to more accurately discriminate plants’ ecological strategies and niche partitioning along environmental gradients

Any plant to survive must occupy an environment with conditions that they can tolerate (Winemiller et al., 2015). Thus, habitat adaptation is determined primarily by abiotic constraints, which shape the species distribution along variable resource scenarios, and therefore, spatial trait segregation provides a starting point for building a niche scheme and performance measurements to predict how species respond to resource gradients (Winemiller et al., 2015). However, despite this spatial niche partitioning, species coexist within the same environment, because niche partitioning might involve some combination of strategies for resource use. In this context, a proposed mechanism of plant coexistence is based on the ‘biogeochemical niche’ hypothesis (Peñuelas et al., 2008). According to this hypothesis, species have an optimal elemental composition resulting from the link between the chemical composition and morphological traits for an optimum plant functioning (Peñuelas et al., 2008). Thus, differences in elemental composition and morpho-structure among co-occurring species may reduce competition, promoting different functional niche (Sardans et al., 2015). Within a determined environment type, the co-occurring species compete in the same space by different resources (not only nutrients, also by water and light) with different intensity, allowing the species coexistence in the same biome.

A great challenge is to identify groups associated with functional traits that allow for successful trait syndromes and ecological predictions. To this end, the assimilation and conservation of resources is a fundamental dimension. Wright et al. (2004) proposed a scheme, known as the ‘leaf economics spectrum,’ that assumes species segregation along environmental resources gradient, according with a plane of trait variation. Hence, this scheme runs from species with a conservative resource-use strategy (i.e., plants of slow growth, with tissues of high-density and low nutrient content) which are more abundant in low-productivity sites, to species from resource-rich sites with the opposite suite of traits, and associated to rapid resource capture (Wright et al., 2001; Villar et al., 2006; Poorter and Garnier, 2007). In this regard, different functional habits associated with leaf life-span (i.e., evergreen or deciduous) or growth forms (i.e., trees or shrubs) are associated with different strategies for resource uptake, which could promote niche partitioning (Mamolos et al., 1995; Niinemets and Kull, 2003; Lusk and Warton, 2007; Sardans et al., 2015; de la Riva et al., 2016a). Nevertheless, there is a large gap in our understanding on the influence of plant traits (leaf morphology and chemistry), leaf habit and growth forms as drivers of species niche partitioning (Sterck et al., 2011).

A variety of multivariate statistical methods have been developed to understand niche partitioning (e.g., Winemiller et al., 2015). However, more powerful tests of niche theories need to move beyond approaches based on species occurrence and instead to focus explicitly on trait-based approaches, recasting the theories in terms of functional space and diversity. The n-dimensional niche space model (Blonder et al., 2014), which is based on the Hutchinson’s multidimensional niche concept (Hutchinson, 1957), allows to quantify niche spaces by assessing the functional trait space that characterize the phenotypic volume occupied by a set of species (Violle and Jiang, 2009; Lamanna et al., 2014). A set of n variables, which represent key and independent biologically axes, are used to create a n-dimensional space, which is defined as “hypervolume.” The main advantage of this method is that the set of points are projected into a hyperspace, defining a high-dimensional shape which may include holes or complex geometrical features. In addition, this method allows quantifying the proportion of the hypervolumes that overlap in different habitats, i.e., the fraction of them sharing the same functional space (Blonder et al., 2014).

In a previous work, carried out at regional scale and at a whole plant level (i.e., quantifying leaf, root, and stem traits), we demonstrated that soil water scarcity led to lower plant functional space and trait segregation (de la Riva et al., 2017). Thus, communities of woody Mediterranean plants from arid and semi-arid shrublands displayed smaller functional space, with a higher degree of overlap between them. That pattern would reflect the species cluster around adaptive peaks, which are defined by sets of trait combinations associated with a given set of environmental attributes (Winemiller et al., 2015). In the present study, based on the biogeochemical niche concept and the use of resources in different proportions, we have focused on leaf traits (morphological and chemical) associated with resource uptake, and explored the niche partitioning among functional habits (leaf life-span and growth form). We have quantified the functional space based on leaf functional traits in a set of 45 woody species from Sierra Morena Mountains (South Spain) growing along a broad gradient of soil resources. Assuming that differences in leaf nutrient composition and morpho-structure reflect different ecological strategies, associated with niche partitioning, and so permitting coexistence (McGroddy et al., 2004; Wright et al., 2004; Chen et al., 2011), we hypothesized that: (i) As a consequence of different plant adaptations to maximize their fitness under determined environmental conditions; we expect a segregation in the foliar chemical composition and morphology among the species along the explored soil resource gradient (mainly soil water availability); (ii) the major differences in leaf chemical and morphological traits are found among species with contrasting growth forms and leaf life-spans. This implies that coexisting species would share their niches by using different ranges and proportions of resources, assuming trade-offs in resource allocation.

Materials and Methods

Study Area and Sampling Design

The Mediterranean forests studied are located in Sierra Morena Mountains, in the south of Spain (Córdoba province). The area is characterized by a continental-Mediterranean climate with cold, wet winters and dry, warm summers. Mean annual temperature is 17.6°C (with maximum values in summer reaching 40°C) and mean annual precipitation is 536 mm (with a 3-month period in summer without rainfall; data from AEMET for the years 1971–20001). Several shrub and arborescent species, such as Cistus albidus and Quercus coccifera, are abundant in drier soils, while broad-leaf deciduous trees, such as Alnus glutinosa and Fraxinus angustifolia, are dominant in moister soils (see Supplementary Table 1 for details of the 45 studied species). One cultivate species (Cydonia oblonga) has been included, because is naturalized in the study area and could potentially alter the biogeochemical niche of the coexistent species. Twelve sampling sites distributed over four different south-facing slopes were selected along a topographic gradient (from ridges to valley bottoms; see Supplementary Figure 1) with the aim of spanning a broad range of variations in soil resource availability, mainly of soil water (Supplementary Figures 2, 3 and de la Riva et al., 2016b). The species composition was recorded by measuring the cover of each woody species intercepted by four 20-m transects in each sampling site.

Trait Measurements

For trait measurements, we selected all the species appearing in the sampling transects, excluding only those with a relative abundance below 1% (in these cases it was difficult to find at least six individuals per species in the sampling site). This gave a total of 45 selected species, many of them appearing in more than one sampling site (Supplementary Table 1). We measured different leaf traits associated with resource uptake and conservation: two key functional traits related to morphology and associated with light capture and growth rate [specific leaf area (SLA)] and stress tolerance [leaf dry matter content (LDMC)], all the macronutrients (N, C, P, K, S, Ca, and Mg) and the isotopic ratio of carbon, which is related with gas exchange and water-use efficiency (Pérez-Harguindeguy et al., 2013).

Six samples of leaves from different individuals per species and site were collected in spring. Specific leaf area (leaf area per unit dry leaf mass; m2 kg-1) and LDMC (dry mass per unit of water-saturated fresh mass; g g-1) were measured according to the methods recommended by Pérez-Harguindeguy et al. (2013). Chemical composition was determined for a mixture of leaves, from six different individuals per species and sampling site (except N and C which were measured in each individual). N and C concentrations were measured using an elemental analyser (Eurovector EA, 3000; EuroVector SpA, Milan, Italy). The macronutrients P, K, S, Ca, and Mg were extracted by wet oxidation with concentrated HNO3 under pressure in a microwave digester, and analyzed by ICP-OES (Domínguez et al., 2012). Carbon isotope ratio (δ13C; ‰) was measured by combustion at 1020°C using a continuous flow isotope-ratio mass spectrometry system by means of Flash HT Plus elemental analyzer coupled to a Delta-V Advantage isotope ratio mass spectrometer via a CONFLO IV interface (Thermo Fisher Scientific, Bremen, Germany). The analytical measurement errors were ±0.1‰ for δ13C.

Data Analysis

For all the statistical analyses, plant species were sorted into different groups according to (I) their topographical position (environment): Ridge Forest (hereafter RF), Middle-slope Forest (hereafter MF) and Riparian Forest (hereafter RiF); (II) leaf habit: Evergreen (hereafter Ev), Deciduous (hereafter De) and Semideciduous (hereafter Sd); and (III) growth form: Shrubs (hereafter Sr), Trees (hereafter T), Arborescent-Shrubs (hereafter ST), and Climber (C) sensu lato. Semideciduous species are considered those that have the capacity to drop their leaves under severe drought conditions (Zunzunegui et al., 2005; Ciccarelli et al., 2016), and therefore they have functional differences with respect evergreens.

A general principal components analysis (PCA) was performed with the whole set of leaf structural and chemical traits (10 leaf variables) for the 100 observations (measurements) made of woody plants (of 45 species, as some species appear in more than one sampling site), to study the degree of leaf trait variation among them. We used linear mixed models (considering species as a random variable) to assess whether the PCA scores of the first, second and third components (see details in Supplementary Table 2) differed among all possible combinations of the following fixed explanatory factors: environment, leaf habit, and growth form.

The total niche space of the community was calculated by the estimation of the n-dimensional hypervolume (Blonder et al., 2014), from the trait space occupied by the total species group in the different categories (environment, leaf habit, and growth form; the number of observations for climbers was not enough to calculate the hypervolume). In order to reduce the number of dimensions (which is recommended for this analysis), we used the first three PCA axes to calculate the hypervolume for each group, using a multidimensional kernel density estimation (KDE) procedure (see Blonder et al., 2014 for mathematical details). The units of the hypervolumes are reported as the standard deviations of centered and scaled log-transformed trait values, raised to the power of the number of trait dimensions used (sdnumber of dimensions). We also calculated the overlap between the hypervolumes of each group with the correlation analysis of the “hypervolume” package, which compares the similarity between different hypervolumes using the Sørensen index (see Blonder et al., 2015). A rarefaction analysis was performed to control for the effects of species richness on the hypervolume. Thus, for each environment type, leaf habit and growth form, we built 100 randomized communities composed of species drawn (nine observations) from the species pool of that group. Then, we calculated the hypervolume of each sample and performed a one-way ANOVA to compare the hypervolumes of the groups, independently of species richness. In addition, to assess the functional trait overlap within and among groups (leaf habit and growth form) within each environment type, we calculated the hypervolumes and the Sørensen index between each pair of the most representative groups along the gradient.

All these analyses were conducted in the R 2⋅10⋅0 statistical platform (R Development Core Team, 2011), using the packages “vegan” (Oksanen et al., 2013), “nlme” (Pinheiro et al., 2015), and “hypervolume” (Blonder et al., 2014).


Species Variation

The overall PCA showed a clear separation among the 45 woody species in the volume defined by PCA1 (accounting for 48.2% variance), PCA2 (12.6% variance), and PCA3 (8.5% variance) (Figure 1). The PCA1 mostly reflected a gradient of increasing LDMC and C, and decreasing SLA, K and Mg; the PCA2 was mostly related in one extreme (negative values) with high N, P, and S and at the opposite extreme (positive values) with high values of Ca and Mg; the PCA3 was represented mainly by variations in N, S, P, and K (see variables scores in Supplementary Table 2).


FIGURE 1. Principal components analysis (PCA) for the 10 leaf traits – morphological (LDMC and SLA), macronutrients (C, N, P, K, Mg, Ca, S) and carbon isotope ratio (C13) – and 45 woody species.

The woody species inhabiting the ridge forest (RF) had different trait values (well-separated according the first and the third PCA axes) from those living at the Riparian forest (RiF) (Figure 2). There was also a separation in the trait PCA between species grouped by leaf habit: Deciduous (De) were different from Sd and Ev. With regards to growth form, there were differences among the arborescent shrubs (for PCA1), and trees (for PCA3) with arborescent shrubs and climbers (Figure 2); while no significant differences were found for PCA2.


FIGURE 2. Differences among the scores of the first (leaf economic spectrum) and third (macronutrients) PCA axes, for species grouped according with environment (RF-Ridge forest, MF-Medium-slope, and RiF-Riparian forest), leaf life-span (Sd, semideciduous; Ev, evergreen; and De, deciduous) and growth forms (Sr, shrubs; ST, arborescent-shrubs; T, trees; and C, climbers). All the results were significant; Linear mixed models, P < 0.05 (Different letters means significant differences).

Functional Space

The results from the n-dimensional hypervolume approach showed that the functional space was greatest for the riparian communities from riparian forest (Figure 3A), the group of deciduous species (Figure 3B), and the growth form “trees” (Figure 3C), in the three plant dimensions (PCA axes 1, 2, and 3). In addition, after standardizing for species richness, the functional space showed significant variation along the topographic gradient: the hypervolume was significantly greater at the riparian forest (P < 0.001; Figure 4). Semideciduous species showed the lowest hypervolume, while deciduous species had the highest. Among the growth form types, trees had the highest hypervolume (Figure 4).


FIGURE 3. Estimated three-dimensional hypervolumes for: (A) Environment (RF, Ridge forest; MF, Medium-slope; and RiF, Riparian forest), (B) Leaf life-span (Sd, semideciduous; Ev, evergreen; and De, deciduous), (C) Growth forms (Sr, shrubs; ST, arborescent-shrubs; T, trees). (D) Sørensen similarity index among hypervolumes. Each plant dimension was based in each of the first three PCA axes (Figure 1).


FIGURE 4. Size of the hypervolume variations for each environment, leaf life-span and growth form. Different letters indicate significant differences (P < 0.05), according to post hoc multiple pairwise comparisons (Tukey’s test). The ANOVA results were significant in all cases (P < 0.001).

Overlapping Niches

The degree of overlapping among the hypervolumes of the different communities was variable, ranging from 28 to 90% (Figure 3D). Attending to the environment, the overlap degree of the ridge forest RF was greater with middle MF than with Riparian forest RiF. Semideciduous type showed the highest overlap with evergreen type, and the lowest overlap with the deciduous group. Among the growth form types, the lowest degree of overlapping was found among the hypervolumes of trees, and arborescent-shrubs. In general, the degree of hypervolume-overlapping among different growth forms was relatively low, which indicate a trait space occupation more variable for growth forms than for environment or leaf habit.

Differences within the same groups were also remarkable (Figure 5). Only evergreen type was represented at the three environment types; it showed low functional niche space overlap (<0.3) with semideciduous (at the ridge) and (<0.2) with deciduous (at the riparian). The functional space was also very different among growth forms within the same environments (overlap < 0.4); with exception of trees and arborescent-shrubs in ridge forests, which had the highest overlap (Figure 5).


FIGURE 5. Sørensen similarity index for the functional niche space among species, grouped by leaf life-span and growth form, within the same habitat. Values are indicated by color according to the rank in the bottom of the figure.


Environment and Niche Partitioning

We have documented a consistent variation in functional space (based on leaf morphology and chemical composition) for woody communities distributed along the topographical gradient. Thus, communities dominated by deciduous trees with faster growth and a predominant acquisitive strategy (i.e., with high values of SLA and leaf nutrient concentration) were characteristic of riparian forests, while semideciduous shrubs and evergreen (arborescent, trees) species, characterized by a conservative strategy (i.e., with high-density leaves and a predominance of carbon-based structural defenses), dominated ridge forests. Hence, the first strong mechanism of species segregation and niche partitioning was related to habitat suitability. Along that local gradient of stress-productivity, plant establishment might be limited physiologically at one end (by stress tolerance) and by competition at the other end. Overall, leaf trait variation followed the patterns commonly associated to leaf economics spectrum (Grime, 1979; Wright et al., 2004).

This general niche differentiation began with the evolutionary accentuation of trait differences during the expansion of deciduous species on Cretaceous (Axelrod, 1966). The evolutionary and environmental drivers have determined this leaf chemical and morphological composition to improve plant functioning (Reich et al., 2003). However, we detected another trait dimension related with leaf concentration in some macronutrients (third PCA axis) that was not clearly aligned with the leaf economics spectrum. Other factors, like site differences or evolutionary history, might explain the lower values of leaf P and K associated to deciduous trees (Watanabe et al., 2007; Auger and Shipley, 2013). In a previous study on leaf chemical traits in 98 Mediterranean woody species, we also observed that leaf P and K were highly conditioned by the environment and phylogeny (de la Riva et al., accepted). These results suggest that leaf trait variation is not aligned along a single acquisition-conservation axis, and highlights the necessity of considering plant traits from other potentially independent leading dimensions such as those related with leaf chemical composition.

Biogeochemical Niches

Within each topographical zone or environment type, the leaf biogeochemical niche partitioning would underlie the species ability to coexist by diverging on leaf nutrient composition and resource uptake. Thus, riparian communities from valley bottom showed the highest leaf biogeochemical space, which supports the existence of a greater number of functionally different species, as far as the acquisition of resources is concerned. In a regional study, the higher functional diversity (at whole plant level) was also found in wetter environments along the aridity gradient (de la Riva et al., 2017). On the contrary, the lowest diversity could be due to resource scarcity limiting the establishment of species that are not physiologically able to tolerate such abiotic constraints; as a consequence, the range of functional traits is reduced within the range of viable traits that allow plants to persist in that arid environment, in detriment of the functional space (de la Riva et al., 2016c). The more productive environments, however, are able to increase resource heterogeneity, and promote the coexistence of the species with different demand of them (de la Riva et al., 2017).

Niche partitioning within the same growth form or leaf life-span may also facilitate the coexistence among woody species. Riparian forests were dominated by deciduous tree species, which in fact was the life-span type with higher biogeochemical and morphological niche functional space. The dense shade created by the trees in this type of environment likely may induce a high among-species competition for light. However, woody species may also exhibit some degree of plasticity in nutrient uptake to respond to this competition. Plant species adapted to these productive environments, where nutrients have usually intermittent availability, might respond better to these temporal changes showing higher capacities for taking up resources and higher nutrient flexibility (Aerts, 1999; Sardans et al., 2015). Thus, our results suggest that differences in niche functional space (implying different strategies of resource capturing) among coexisting species would potentially allow to buffer competition in these resource-rich environments.

On the contrary, we detected smaller functional space in more stressful environments (i.e., ridge forests). This result could be explained because species growing in poor environments comprise traits which lead to nutrient retention and higher nutrient-use efficiency, which seem to confer a lower capacity to change their functioning in response to environmental changes (Aerts, 1999). In addition, plant size seems to be related with root expansion and nutrient uptake. Thus smaller species (more frequent in resource-poor environments) tend to show lower values of leaf nutrients than trees as a result of their lower capacity to maintain larger root systems and explore larger soil volumes to uptake nutrients (Niinemets and Kull, 2003 and references therein). As a consequence, shrub and arborescent-shrub species had lower niche space than trees. The combined effects of the slow population dynamics in poor environments (Wright, 2002), and their lower capacity to alter their elemental composition independently of the nutrient pulses (Sardans et al., 2015) would facilitate the coexistence of ecologically equivalent species.

Growth Forms and Niche Partitioning

Differences in growth form and leaf life-span may facilitate the coexistence of woody species and their niche partitioning, as showed by the lower functional overlap among them. In spite of the general patterns found along the topographical gradient (e.g., dominance of deciduous trees at riparian forests, and evergreen arborescent-shrubs at ridge forests), strong divergence were also found between their functional spaces within the same habitat type. Thus, the low degree of overlap between the functional space of deciduous trees and evergreen arborescent-shrubs at the riparian forest reflects a niche partitioning, probably promoted by light competition and resource partitioning. The dense-shaded forest understory is a relatively stressful habitat for plants (in terms of light availability), which promotes a resource gradient also related with leaf economics spectrum (Lusk and Warton, 2007). On the one hand, deciduous trees (such as Alnus glutinosa and Fraxinus angustifolia) are effective competitors, since they are able to rapidly acquire nutrients and water, that allow them to grow faster and compete more efficiently for light interception (Lambers and Poorter, 1992; Wright et al., 2001); on the other hand, evergreen shrub species (such as Rubus ulmifolius) had higher LMA, lower nutrient concentration, and are associated with shade tolerance (Williams et al., 1989; Lusk, 2002).

In the opposite extreme of the environmental gradient – at the ridge forests – woody species with semideciduous life-span seem to display more acquisitive patterns than evergreens. Shorter leaf life-span is related to higher relative nutrient requirement and lower resistance to physical hazards (Ryser, 1996; Pérez-Harguindeguy et al., 2013). In this sense, the habit of semideciduous leaf could be considered a facultative strategy, typical of shrubs in dry Mediterranean conditions (Axelrod, 1966; Zunzunegui et al., 2005). That strategy allows them to be effective competitors during favorable conditions, while they shed partly or completely their leaves during summer, so reducing water loss by transpiration (Zunzunegui et al., 2005; Ciccarelli et al., 2016; Bongers et al., 2017). Therefore, the shrub species might be segregated along the resource gradient through a trade-off between growth rate and survival. In summary, changes in plant structure such as growth form or leaf life-span, might reduce competition by differences in the ability to capture and use resources, supporting that limiting similarity was operating within the communities (Stubbs and Wilson, 2004).

Author Contributions

EdlR designed the study, harvested samples, collected the data, analyzed the data, wrote the first draft, and prepared the manuscript. TM, CV, RV, and IP-R designed the study and prepared the manuscript.


This study was funded by the Spanish MEC projects DIVERBOS (CGL2011-30285-C02-01 and C02-02), ECO-MEDIT (CGL2014-53236-R), and RESTECO (CGL2014-52858-R), and European FEDER funds. CV was supported by the European Research Council (ERC) Starting Grant Project “Ecophysiological and biophysical constraints on domestication in crop plants” (Grant ERC-StG-2014-639706-CONSTRAINTS).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


Thanks are due to C. Navarro and M. Olmo for help during field sampling and trait measurements. Leaf multielemental analysis was carried out at the IRNAS, CSIC, isotopic analysis at the LIE, EBD-CSIC, and elemental C and N in the SCAI of the University of Córdoba.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.01242/full#supplementary-material


  1. ^ http://www.aemet.es/es/serviciosclimaticos/datosclimatologicos


Aerts, R. (1999). Interspecific competition in natural plant communities: mechanisms, trade-offs and plant-soil feedbacks. J. Exp. Bot. 50, 29–37. doi: 10.1093/jxb/50.330.29

CrossRef Full Text | Google Scholar

Auger, S., and Shipley, B. (2013). Inter-specific and intra-specific trait variation along short environmental gradients in an old-growth temperate forest. J. Veg. Sci. 24, 419–428. doi: 10.1111/j.1654-1103.2012.01473.x

CrossRef Full Text | Google Scholar

Axelrod, D. I. (1966). Origin of deciduous and evergreen habits in temperate forests. Evolution 20, 1–15. doi: 10.1111/j.1558-5646.1966.tb03339.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Blonder, B., Lamanna, C., Violle, C., and Enquist, B. J. (2014). The n-dimensional hypervolume. Glob. Ecol. Biogeogr. 23, 595–609. doi: 10.1111/geb.12146

CrossRef Full Text | Google Scholar

Blonder, B., Nogués-Bravo, D., Borregaard, M. K., Donoghue, I. I., John, C., Jørgensen, P. M., et al. (2015). Linking environmental filtering and disequilibrium to biogeography with a community climate framework. Ecology 96, 972–985. doi: 10.1890/14-0589.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Bongers, F. J., Olmo, M., Lopez-Iglesias, B., Anten, N., and Villar, R. (2017). Drought responses, phenotypic plasticity and survival of Mediterranean species in two different microclimatic sites. Plant Biol. 19, 386–395. doi: 10.1111/plb.12544

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, F. S., Niklas, K. J., and Zeng, D. H. (2011). Important foliar traits depend on species-grouping: analysis of a remnant temperate forest at the Keerqin Sandy Lands, China. Plant Soil 340, 337–345. doi: 10.1007/s11104-010-0606-9

CrossRef Full Text | Google Scholar

Ciccarelli, D., Picciarelli, P., Bedini, G., and Sorce, C. (2016). Mediterranean Sea cliff plants: morphological and physiological responses to environmental conditions. J. Plant Ecol. 9, 153–164. doi: 10.1093/jpe/rtv042

CrossRef Full Text | Google Scholar

de la Riva, E. G., Olmo, M., Poorter, H., Ubera, J. L., and Villar, R. (2016a). Leaf mass per area (LMA) and its relationship with leaf structure and anatomy in 34 Mediterranean woody species along a water availability gradient. PLoS ONE 11:e0148788. doi: 10.1371/journal.pone.0148788

PubMed Abstract | CrossRef Full Text

de la Riva, E. G., Pérez-Ramos, I. M., Tosto, A., Navarro-Fernández, C. M., Olmo, M., Marañón, T., et al. (2016b). Disentangling the relative importance of species occurrence, abundance and intraspecific variability in community assembly: a trait-based approach at the whole-plant level in Mediterranean forests. Oikos 125, 354–363. doi: 10.1111/oik.01875

CrossRef Full Text | Google Scholar

de la Riva, E. G., Tosto, A., Pérez-Ramos, I. M., Navarro-Fernández, C. M., Olmo, M., Anten, N. P., et al. (2016c). A plant economics spectrum in Mediterranean forests along environmental gradients: Is there coordination among leaf, stem and root traits? J. Veg. Sci. 27, 187–199. doi: 10.1111/jvs.12341

CrossRef Full Text | Google Scholar

de la Riva, E. G., Violle, C., Pérez-Ramos, I. M., Marañón, T., Navarro-Fernández, C. M., Olmo, M., et al. (2017). A multidimensional functional trait approach reveals the imprint of environmental stress in Mediterranean woody communities. Ecosystems doi: 10.1007/s10021-017-0147-7 [Epub ahead of print].

CrossRef Full Text | Google Scholar

Domínguez, M. T., Aponte, C., Pérez-Ramos, I. M., García, L. V., Villar, R., and Marañon, T. (2012). Relationships between leaf morphological traits, nutrient concentrations and isotopic signatures for Mediterranean woody plant species and communities. Plant Soil 357, 407–424. doi: 10.1007/s11104-012-1214-7

CrossRef Full Text | Google Scholar

Grime, J. P. (1979). Plant Strategies and Vegetation Processes. Chichester: Wiley.

Google Scholar

HilleRisLambers, J., Adler, P. B., Harpole, W. S., Levine, J. M., and Mayfield, M. M. (2012). Rethinking community assembly through the lens of coexistence theory. Annu. Rev. Ecol. Evol. Syst. 43, 227–248. doi: 10.1146/annurev-ecolsys-110411-160411

CrossRef Full Text | Google Scholar

Hutchinson, G. E. (1957). Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 22, 415–427. doi: 10.1101/SQB.1957.022.01.039

CrossRef Full Text | Google Scholar

Lamanna, C., Blonder, B., Violle, C., Kraft, N. J. B., Sandel, B., Šímová, I., et al. (2014). Functional trait space and the latitudinal diversity gradient. Proc. Natl. Acad. Sci. U.S.A. 111, 13745–13750. doi: 10.1073/pnas.1317722111

PubMed Abstract | CrossRef Full Text | Google Scholar

Lambers, H., and Poorter, H. (1992). Inherent variation in growth rate between higher plants: a search for physiological causes and ecological consequences. Adv. Ecol. Res. 23, 187–261. doi: 10.1016/S0065-2504(08)60148-8

CrossRef Full Text | Google Scholar

Lusk, C. H. (2002). Leaf area accumulation helps juvenile evergreen trees tolerate shade in a temperate rainforest. Oecologia 132, 188–196. doi: 10.1007/s00442-002-0974-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Lusk, C. H., and Warton, D. I. (2007). Global meta-analysis shows that relationships of leaf mass per area with species shade tolerance depend on leaf habit and ontogeny. New Phytol. 176, 764–774. doi: 10.1111/j.1469-8137.2007.02264.x

PubMed Abstract | CrossRef Full Text | Google Scholar

MacArthur, R., and Levins, R. (1967). The limiting similarity, convergence, and divergence of coexisting species. Am. Nat. 101, 377–385. doi: 10.1086/282505

CrossRef Full Text | Google Scholar

Maire, V., Gross, N., Börger, L., Proulx, R., Wirth, C., and da Silveira Pontes, L. (2012). Habitat filtering and niche differentiation jointly explain species relative abundance within grassland communities along fertility and disturbance gradients. New Phytol. 196, 497–509. doi: 10.1111/j.1469-8137.2012.04287.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mamolos, A. P., Elisseou, G. K., and Veresoglou, D. S. (1995). Depth of root activity of coexisting grassland species in relation to N and P additions, measured using nonradioactive tracers. J. Ecol. 83, 643–652. doi: 10.2307/2261632

CrossRef Full Text | Google Scholar

Mason, N. W., de Bello, F., Doležal, J., and Lepš, J. (2011). Niche overlap reveals the effects of competition, disturbance and contrasting assembly processes in experimental grassland communities. J. Ecol. 99, 788–796. doi: 10.1111/j.1365-2745.2011.01801.x

CrossRef Full Text | Google Scholar

McGill’s, B. J., Enquist, B. J., Weiher, E., and Westoby, M. (2006). Rebuilding community ecology from functional traits. Trends Ecol. Evol. 21, 178–185. doi: 10.1016/j.tree.2006.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

McGroddy, M. E., Daufresne, T., and Hedin, L. O. (2004). Scaling of C: N: P stoichiometry in forests worldwide: implications of terrestrial Redfield-type ratios. Ecology 85, 2390–2401. doi: 10.1890/03-0351

CrossRef Full Text | Google Scholar

Niinemets,Ü, and Kull, K. (2003). Leaf structure vs. nutrient relationships vary with soil conditions in temperate shrubs and trees. Acta Oecol. 24, 209–219. doi: 10.1016/S1146-609X(03)00094-8

CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., et al. (2013). Package ‘Vegan’. Community Ecology Package, Version, 2.

Google Scholar

Peñuelas, J., Sardans, J., Ogaya, R., and Estiarte, M. (2008). Nutrient stoichiometric relations and biogeochemical niche in coexisting plant species: effect of simulated climate change. Pol. J. Ecol. 56, 613–622.

Google Scholar

Pérez-Harguindeguy, N., Díaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234. doi: 10.1071/BT12225

CrossRef Full Text | Google Scholar

Pinheiro, J., Bates, D., DebRoy, S., and Sarkar, D. (2015). nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 3.1-121. Available at: http://CRAN.R-project.org/package=nlme

Google Scholar

Poorter, H., and Garnier, E. (2007). “Ecological significance of inherent variation in relative growth rate and its components,” in Functional Plant Ecology, 2nd Edn, eds F. I. Pugnaire and F. Valladares (Boca Raton, FL: CRC Press), 67–100.

Google Scholar

Reich, P. B., Wright, I. J., Cavender-Bares, J., Craine, J. M., Oleksyn, J., Westoby, M., et al. (2003). The evolution of plant functional variation: traits, spectra, and strategies. Int. J. Plant Sci. 164, 143–164. doi: 10.1086/374368

CrossRef Full Text | Google Scholar

R Development Core Team (2011) R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing

Ryser, P. (1996). The importance of tissue density for growth and life span of leaves and roots: a comparison of five ecologically contrasting grasses. Funct. Ecol. 10, 717–723. doi: 10.2307/2390506

CrossRef Full Text | Google Scholar

Sardans, J., Janssens, I. A., Alonso, R., Veresoglou, S. D., Rillig, M. C., Sanders, T. G. M., et al. (2015). Foliar elemental composition of European forest tree species associated with evolutionary traits and present environmental and competitive conditions. Glob. Ecol. Biogeogr. 24, 240–255. doi: 10.1111/geb.12253

CrossRef Full Text | Google Scholar

Sterck, F., Markesteijn, L., Schieving, F., and Poorter, L. (2011). Functional traits determine trade-offs and niches in a tropical forest community. Proc. Natl. Acad. Sci. U.S.A. 108, 20627–20632. doi: 10.1073/pnas.1106950108

PubMed Abstract | CrossRef Full Text | Google Scholar

Stubbs, W. J., and Wilson, J. B. (2004). Evidence for limiting similarity in a sand dune community. J. Ecol. 92, 557–567. doi: 10.1111/j.0022-0477.2004.00898.x

CrossRef Full Text | Google Scholar

Villar, R., Robleto, J. R., De Jong, Y., and Poorter, H. (2006). Differences in construction costs and chemical composition between deciduous and evergreen woody species are small as compared to differences among families. Plant Cell Environ. 29, 1629–1643. doi: 10.1111/j.1365-3040.2006.01540.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Violle, C., and Jiang, L. (2009). Towards a trait-based quantification of species niche. J. Plant Ecol. 2, 87–93. doi: 10.1093/jpe/rtp007

CrossRef Full Text | Google Scholar

Watanabe, T., Broadley, M. R., Jansen, S., White, P. J., Takada, J., Satake, K., et al. (2007). Evolutionary control of leaf element composition in plants. New Phytol. 174, 516–523. doi: 10.1111/j.1469-8137.2007.02078.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, K., Field, C. B., and Mooney, H. A. (1989). Relationships among leaf construction cost, leaf longevity, and light environments in rain-forest plants of the genus Piper. Am. Nat. 133, 199–211. doi: 10.1086/284910

CrossRef Full Text | Google Scholar

Winemiller, K. O., Fitzgerald, D. B., Bower, L. M., and Pianka, E. R. (2015). Functional traits, convergent evolution, and periodic tables of niches. Ecol. Lett. 18, 737–751. doi: 10.1111/ele.12462

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, I. J., Reich, P. B., and Westoby, M. (2001). Strategy shifts in leaf physiology, structure and nutrient content between species of high- and low-rainfall and high- and low-nutrient habitats. Funct. Ecol. 15, 423–434. doi: 10.1046/j.0269-8463.2001.00542.x

CrossRef Full Text | Google Scholar

Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. Nature 428, 821–827. doi: 10.1038/nature02403

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. S. (2002). Plant diversity in tropical forests: a review of mechanisms of species coexistence. Oecologia 130, 1–14. doi: 10.1007/s004420100809

PubMed Abstract | CrossRef Full Text | Google Scholar

Zunzunegui, M., Barradas, M. D., Ain-Lhout, F., Clavijo, A., and Novo, F. G. (2005). To live or to survive in Doñana dunes: adaptive responses of woody species under a Mediterranean climate. Plant Soil 273, 77–89. doi: 10.1007/s11104-004-6806-4

CrossRef Full Text | Google Scholar

Keywords: functional ecology, growth forms, hypervolume, carbon isotope ratio, LDMC, leaf habit, nutrients, SLA

Citation: de la Riva EG, Marañón T, Violle C, Villar R and Pérez-Ramos IM (2017) Biogeochemical and Ecomorphological Niche Segregation of Mediterranean Woody Species along a Local Gradient. Front. Plant Sci. 8:1242. doi: 10.3389/fpls.2017.01242

Received: 29 April 2017; Accepted: 30 June 2017;
Published: 19 July 2017.

Edited by:

Raquel Esteban, University of the Basque Country, Spain

Reviewed by:

Manuel J. Macía, Universidad Autonoma de Madrid, Spain
Georgina Conti, Instituto Multidisciplinario de Biologia Vegetal, Argentina

Copyright © 2017 de la Riva, Marañón, Violle, Villar and Pérez-Ramos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Enrique G. de la Riva, enga70@gmail.com