# Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient

Cristina C. Bastias, Claire Fortunel, Fernando Valladares, Christopher Baraloto, Raquel Benavides, William Cornwell, Lars Markesteijn, Alexandre A. de Oliveira, Jeronimo B. B. Sansevero, Marcel C. Vaz, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Cristina C. Bastias, Claire Fortunel, Fernando Valladares, Christopher Baraloto, Raquel Benavides, et al. Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient. PLoS ONE, 2017, 12 (2), pp.1-16. 10.1371/journal.pone.0172495. hal-02625462

HAL Id: hal-02625462

https://hal.inrae.fr/hal-02625462

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient

Cristina C. Bastias¹*, Claire Fortunel²,³, Fernando Valladares¹,⁴, Christopher Baraloto⁵, Raquel Benavides¹,⁶, William Cornwell⁷, Lars Markesteijn⁸,⁹,¹⁰, Alexandre A. de Oliveira¹¹, Jeronimo B. B. Sansevero¹²,¹³, Marcel C. Vaz²,³, Nathan J. B. Kraft²,³

¹ Departamento de Biogeografía y Cambio Global. Museo Nacional de Ciencias Naturales- CSIC, Madrid, Spain, ² Department of Biology. University of Maryland, College Park, MD, United States of America,

³ Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA, United States of America, ⁴ Área de Biodiversidad y Conservación, Universidad Rey Juan Carlos, Móstoles, Madrid, Spain,

⁵ INRA, UMR "Ecologie des Forets de Guyane", Kourou Cedex 97387, French Guiana, ⁶ Institut für Biologie, Geobotanik, Albert-Ludwigs-Universität Freiburg, Freiburg, Germany, ⁷ Evolution & Ecology Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Sydney, New South Wales 2052, Australia, ⁸ Smithsonian Tropical Research Institute, Balboa, Ancón, Panamá,

⁹ Department of Zoology, University of Oxford, South Parks Road, Oxford, United Kingdom, ¹⁰ School of Environment, Natural Resources and Geography, Bangor University, Bangor LL57 2DG, United Kingdom,

¹¹ Departamento de Ecologia, Universidade de São Paulo, São Paulo, SP, Brasil, ¹² Instituto de Pesquisas Jardim Botânico do Rio de Janeiro, Jardim Botânico, Rio de Janeiro, Brazil, ¹³ Universidade Federal Rural do Rio de Janeiro–UFRRJ. Departamento de Ciencias Ambientais–DCA. Instituto de Florestas–, Seropédica, Rio de Janeiro, Brazil

* crbasc@gmail.com

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

## OPEN ACCESS

**Citation:** Bastias CC, Fortunel C, Valladares F, Baraloto C, Benavides R, Cornwell W, et al. (2017) Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient. PLoS ONE 12(2): e0172495. doi:10.1371/journal.pone.0172495

**Editor:** Sylvain Delzon, INRA—University of Bordeaux, FRANCE

**Received:** June 29, 2016

**Accepted:** February 6, 2017

**Published:** February 27, 2017

**Copyright:** This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** Leaf data used to come to this manuscript resulted from many different studies involving different funding sources: European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 265171; Spanish-funded project REMEDINAL3eCM (S2013/MAE-2719); Ecometas excellence network (CGL2014-53840-REDT); a FPU grant by the

## Abstract

Disentangling the mechanisms that shape community assembly across diversity gradients is a central matter in ecology. While many studies have explored community assembly through species average trait values, there is a growing understanding that intraspecific trait variation (ITV) can also play a critical role in species coexistence. Classic biodiversity theory hypothesizes that higher diversity at species-rich sites can arise from narrower niches relative to species-poor sites, which would be reflected in reduced ITV as species richness increases. To explore how ITV in woody plant communities changes with species richness, we compiled leaf trait data (leaf size and specific leaf area) in a total of 521 woody plant species from 21 forest communities that differed dramatically in species richness, ranging from boreal to tropical rainforests. At each forest, we assessed ITV as an estimate of species niche breadth and we quantified the degree of trait overlap among co-occurring species as a measure of species functional similarity. We found ITV was relatively invariant across the species richness gradient. In addition, we found that species functional similarity increased with diversity. Contrary to the expectation from classic biodiversity theory, our results rather suggest that neutral processes or equalizing mechanisms can be acting as potential drivers shaping community assembly in hyperdiverse forests.
Spanish Government (AP2010-5600); a Marie
Curie IEF fellowship (FP7-PEOPLE-2011-IEF. no.
302445); a Investissement d'Avenir' grant
managed by Agence Nationale de la Recherche
(CEBA, ref. ANR-10-LABX-0025); a collaborative
NSF DEB-0743103/0743800 and Institut National
de la Recherche Agronomique (INRA); Instituto
Interamericano para Pesquisas em Mudanças
Globais (IAI); CAPES; Wageningen graduate
school, Production Ecology and Resource
Conservation (PE&RC); Center for Tropical Forest
Science and the University of California, Berkeley,
Department of Integrative, Biology; a NSF graduate
fellowship and a DDIG.

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

**Introduction**

The relative importance of ecological factors in shaping plant communities across species diversity gradients is the subject of longstanding debate in ecology [1–4] that has been recently invigorated by the lens of functional trait diversity [5,6]. The use of traits in a community ecology context hinges on the hypothesis that there is a link between traits and the breadth and position of species’ realized niches [7,8]. Trait-based studies have often used a trait mean approach (i.e. assigning all conspecific individuals a species average trait value) to examine community assembly mechanisms [9,10]. The implicit assumption in many of these studies is that interspecific trait differences are much larger than intraspecific trait differences [11,12]. However, there is increasing evidence that community assembly at local scales depends critically on the extent of intraspecific trait variation (ITV) [13–16]. Recently, the scientific community has reconsidered the importance of ITV [17–19] and its non-negligible contribution to the total trait variability, being sometimes as important as interspecific trait variation [20,21]. Even when interspecific trait differences are larger, incorporating ITV can improve the answer to key questions about the assembly and functioning of plant communities [15,22].

The extent of ITV among species in a community is expected to vary depending on community attributes such as the number of co-occurring species or the community trait diversity [13,19]. Previous studies have suggested that ITV should be greater in species-poor than species-rich communities [23,24]. If the biotic pressure via competitive interactions is lower in species-poor than in species-rich communities, conspecific individuals in species-poor communities could occupy a greater extension of available trait space (i.e. substantial extent of ITV) without increasing interspecific interactions (Fig 1A and 1A'). As the number of co-occurring species increases under the assumption of all co-occurring species with equal fitness (i.e. flat fitness landscape), species' trait breadths are expected to be reduced (i.e. decline of ITV) to accommodate more species without increasing the potential for interspecific competition by resource use, consistent with classical niche theory [19,25] (Fig 1B). However, if the extent of ITV does not change from species-poor to species-rich communities or even increases (for example, see 'individual variation' theory by Clark [26]) and the trait range within a community in turn does not increase when species richness increases (Fig 1C), species' trait overlaps would be expected to increase in more diverse communities (Fig 1C'). If traits map to resource use (and stabilizing niche differences, sensu Chesson [27]), species can coexist more readily by being functionally distinct, thereby promoting trait dissimilarity among species for coexistence [25]. However, some have argued that the lack of interspecific dissimilarity could lead to neutral dynamics (i.e. all individuals are considered functionally equivalent), reducing or removing the role of niche differences in shaping community assembly outcomes [28,29], or alternatively, communities structured primarily by the acting of equalizing mechanisms in a non-neutral model [30]. On the other hand, under an alternative niche differentiation model, specifically if species are differentiating along a landscape represented by combinations of peaks (high-fitness) and deep valleys (low-fitness) (i.e. multi-peak fitness landscape), a decrease of ITV is not predicted with species richness [31] despite a limiting similarity principle playing a role to select species on each peak.

While the relationships between species richness, ITV and trait overlap have important implications for community assembly and the maintenance of species richness [33], they have been poorly studied due to the substantial effort required to measure ITV, especially in species-rich systems [34]. Moreover, the few studies testing these relationships so far have shown contrasting results. For instance, some studies have found a decrease in the extent of ITV in specific leaf area (SLA) accompanied by declining SLA overlap with increased species richness [35,36]; whereas others have found no change in ITV in SLA or an increase of ITV and overlap
<figure><img src="image_5.png" /><figcaption>**Fig 1.** A schematic of possible changes in the extent of ITV and the degree of trait overlap (measured as the similarity between species; [32]) when species richness increases assuming a flat fitness landscape. (A) A substantial extent of ITV is expected in species-poor communities without increasing interspecific interactions since the biotic pressure is low. Therefore, low values of trait overlap are expected in species-poor communities leading to a right-skewed distribution (A'). (B) The extent of ITV is reduced, remaining trait overlap fairly constant to accommodate more species. It translates into a right-skewed distribution to low trait overlap values in species-rich forests (B'). (C) Trait overlap is increased without a change in the extent of ITV, resulting in greater degree of functional similarity among individuals. A left-skewed distribution of trait overlap values would be expected in species-rich forests (C').</figcaption></figure>

doi:10.1371/journal.pone.0172495.g001

using multiple approach (included SLA) with species richness [37]. These contradictory results show evidence that there is a need to improve our understanding of these relationships between species richness, ITV and trait overlap, exploring all possible approaches for a given trait in combination with the attributes of the study system in order to infer community assembly mechanisms across diversity gradients.

Here, we explore how the extent of ITV and trait overlap changes across a broad species richness gradient. We compile leaf trait data for 521 woody plant species found in 21 forest communities that varied in species richness from 1 to 284 species per hectare (from boreal and temperate to tropical forests). In accordance with the classic niche theory under a flat fitness landscape scenario, we hypothesize that the extent of ITV will decrease when more species are added to the community (Fig 1B) [25,27,38]. By extension, this hypothesis would suggest a reduction of the contribution of local intraspecific variation to the total trait diversity of the community with increasing species richness. Alternatively, if ITV does not decline and the overall length of the trait gradient does not increase in diverse communities, we expect a higher overlap in trait distributions among species with increasing species richness (Fig 1C). This would translate into a left-skewed distribution of trait overlap values, i.e. higher median values of trait overlap as well as a higher proportion of pairs of species with high trait overlap in species-rich than species-poor forests (Fig 1C'). By extension, ITV could be expected to contribute to the total trait diversity similarly or even more than interspecific trait differences in rich communities.
# Materials and methods

## Data collection

We used two leaf traits to quantify ITV: leaf size refers to individual leaf area (cm²) and specific leaf area (SLA), defined as the ratio of leaf area to dry mass (cm²·g⁻¹) [39]. We focused on these traits, firstly, because different empirical studies of woody plants from a wide range of environments have shown that SLA and leaf size are weakly or not correlated across species or at species level [40,41]. Secondly, these key leaf traits have been widely used in species distribution studies across gradients to predict future shifts in individual species distributions or even in species-realized niches due to their strong response to abiotic and biotic changes [41,42]. Moreover, SLA and leaf size are traits with an important ecological significance in relation to plant economics and plant resource acquisition: plants investing in greater SLA values increase light-capture efficiency, but are more vulnerable to high temperatures, drought, nutrient-limitation and herbivory [43]. On the other hand, leaf size is more related to the thermal conductance of the leaf boundary layer. Smaller leaf size helps to keep optimal leaf temperature and a higher water balance efficiency, especially under high solar radiation and low water availability conditions [44]. Finally, these traits can easily be measured on a large number of individuals.

We measured leaf traits in 3712 individuals of 521 woody plant species in 21 forest communities with contrasting species richness, ranging from the species-poor boreal and temperate European forests (1–4 species/ha) to hyperdiverse rainforests in Ecuador and Brazil (>200 species/ha) (Table 1). Species richness was calculated for each forest as the number of species with a diameter at breast height (DBH) ≥ 10 cm as it is customary [45–48]. We also obtained climatic variables (the mean annual temperature and annual precipitation) per each forest community from the Worldclim global climate models [49]. We included leaf trait data for species from a given community using the following criteria: (i) at least five individuals per species were measured in each forest in order to estimate ITV (see [50]), but it was higher whenever was possible (Table 1). (ii) To minimize the influence of ontogenetic variation on ITV, we restricted our sampling to understory individuals (saplings and small-stature trees), whose DBH is smaller than 20 cm and height less than 10 m. To restrict the sampling to understory individuals may imply either an underestimation or overestimation in the number of sampled species compared to the species richness of the community (calculated taking into account individuals with a DBH >10cm). An underestimation may occur since we are sub-sampling woody species from a forest layer, but also an overestimation since the arrival of propagules from neighbor canopies can lead to new established individuals in the understory of the canopy of the focal community. (iii) Due to the considerable effort required in sampling at ITV level in hyperdiverse forests, we only conducted trait measurements on a subset of the common species in tropical forests (those species that we found at least 5 individuals) (see Table 1 for the number of sampled species in each forest). This restriction in tropical forests may result in an underestimation of the extent of ITV and trait overlap among species in these forests since we do not have trait data available for the whole range of naturally occurring trait variation (i.e. we do not consider those rare species that occupy unique trait space compared to common species [51,52]). Besides, common species could show less trait variation [44], which further contributes to underestimating actual ITV. (iv) Also, it is important to mention that trait data used here were previously collected for other specific goals (see Table 1 for original references), but in general, individuals and leaves were chosen under standardized abiotic conditions (i.e. recently matured and fully expanded leaves) [39]. As a consequence, we may underestimate the actual ITV since ITV measured here is controlled for two main sources of variation: phenotypic plasticity in response to local abiotic conditions (i.e. we biased ITV towards natural standardized conditions) and ontogenetic variations (i.e. we biased ITV at a
<table><caption>Table 1. Description of each forest community included in the study.</caption><thead><tr><th>Type of forest</th><th>Location</th><th>Latitude</th><th>Longitude</th><th>Mean annual T<sup>ra</sup> (°C)</th><th>Pp (mm)</th><th>SR</th><th>Representative families of sampled species</th><th>No. species measured with ≥ 5 individuals</th><th>Individuals measured per species (Min, Max)</th><th>Original references</th></tr></thead><tbody><tr><td>Boreal forest</td><td>Joensuu, North Karelia (Finland)</td><td>62.616</td><td>29.89</td><td>2.1</td><td>628</td><td>1</td><td>Pinaceae, Betulaceae</td><td>3</td><td>[19, 55]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Boreal forest</td><td>Joensuu, North Karelia (Finland)</td><td>62.504</td><td>29.76</td><td>2.1</td><td>628</td><td>2</td><td>Pinaceae, Betulaceae</td><td>2</td><td>[5, 64]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Boreal forest</td><td>Joensuu, North Karelia (Finland)</td><td>62.558</td><td>30.16</td><td>2.1</td><td>628</td><td>3</td><td>Pinaceae, Betulaceae</td><td>2</td><td>[7, 8]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mountainous beech forest</td><td>Carpathian mountains (Romania)</td><td>47.295</td><td>26.05</td><td>5.6</td><td>689</td><td>1</td><td>Pinaceae, Fagaceae, Sapindaceae</td><td>4</td><td>[11, 54]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mountainous beech forest</td><td>Carpathian mountains (Romania)</td><td>47.294</td><td>26.05</td><td>5.6</td><td>689</td><td>2</td><td>Pinaceae, Fagaceae, Sapindaceae</td><td>4</td><td>[8, 60]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mountainous beech forest</td><td>Carpathian mountains (Romania)</td><td>47.292</td><td>26.05</td><td>5.6</td><td>689</td><td>3</td><td>Pinaceae, Fagaceae, Sapindaceae</td><td>3</td><td>[10, 37]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mountainous beech forest</td><td>Carpathian mountains (Romania)</td><td>47.291</td><td>26.05</td><td>5.6</td><td>689</td><td>4</td><td>Pinaceae, Fagaceae, Sapindaceae</td><td>2</td><td>[5, 20]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mediterranean mixed forest</td><td>Alto Tajo Natural Park (Spain)</td><td>40.731</td><td>-2.25</td><td>9.9</td><td>533</td><td>1</td><td>Pinaceae, Fagaceae</td><td>4</td><td>[6, 57]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mediterranean mixed forest</td><td>Alto Tajo Natural Park (Spain)</td><td>40.713</td><td>-2.19</td><td>9.9</td><td>533</td><td>2</td><td>Pinaceae, Fagaceae</td><td>4</td><td>[21, 70]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Mediterranean mixed forest</td><td>Alto Tajo Natural Park (Spain)</td><td>40.698</td><td>-2.13</td><td>9.9</td><td>533</td><td>3</td><td>Pinaceae, Fagaceae</td><td>2</td><td>[5, 10]</td><td>(Bastias et al. unpublished data)</td></tr><tr><td>Tropical lowland dry deciduous forest</td><td>Inpa, Concepcion, Santa Cruz (Bolivia)</td><td>-16.117</td><td>-61.72</td><td>23.5</td><td>1124</td><td>34</td><td>Fabaceae, Flacourtiaceae, Euphorbiaceae</td><td>52</td><td>[5, 10]</td><td>[48]</td></tr><tr><td>Riparian, chaparral, broadleaf evergreen forest</td><td>Jasper Ridge Biological Preserve (California, USA)</td><td>37.4</td><td>-122.25</td><td>13.8</td><td>598</td><td>54</td><td>Fagaceae, Rosaceae, Rhamnaceae</td><td>43</td><td>[5, 42]</td><td>[42]</td></tr><tr><td>Tropical lowland semi-deciduous seasonal moist forest</td><td>Soberania National Park (Panama)</td><td>9.162</td><td>-79.75</td><td>26</td><td>2553</td><td>131</td><td>Fabaceae, Piperaceae, Rubiaceae</td><td>16</td><td>[5, 6]</td><td>(Markesteijn, unpublished data)</td></tr><tr><td>Lowland tropical rainforest</td><td>Acarouany (French Guiana)</td><td>5.544</td><td>-53.81</td><td>26.5</td><td>2237</td><td>148</td><td>Annonaceae, Burseraceae, Lecythidaceae</td><td>11</td><td>[3, 22]</td><td>[10,53]</td></tr><tr><td>Lowland tropical rainforest</td><td>Paracou (French Guiana)</td><td>5.272</td><td>-52.93</td><td>25.8</td><td>2821</td><td>150</td><td>Euphorbiaceae, Fabaceae, Lecythidaceae</td><td>35</td><td>[3, 39]</td><td>[10,53]</td></tr><tr><td>Lowland tropical rainforest</td><td>BAFOG (French Guiana)</td><td>5.494</td><td>-53.99</td><td>26.4</td><td>2460</td><td>156</td><td>Annonaceae, Burseraceae, Lecythidaceae</td><td>11</td><td>[3, 25]</td><td>[10,53]</td></tr><tr><td>Lowland tropical rainforest</td><td>Nouragues (French Guiana)</td><td>4.087</td><td>-52.67</td><td>24.8</td><td>3337</td><td>197</td><td>Lecythidaceae, Malvaceae, Sapotaceae</td><td>24</td><td>[3, 25]</td><td>[10,53]</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Type of forest</th><th>Location</th><th>Latitude</th><th>Longitude</th><th>Mean annual T<sup>ra</sup> (°C)</th><th>Pp (mm)</th><th>SR</th><th>Representative families of sampled species</th><th>No. species measured with ≥ 5 individuals</th><th>Individuals measured per species (Min, Max)</th><th>Original references</th></tr></thead><tbody><tr><td>Lowland tropical rainforest</td><td>Montagne Tortue (French Guiana)</td><td>4.219</td><td>-52.41</td><td>24.6</td><td>3591</td><td>213</td><td>Sapotaceae</td><td>14</td><td>[3, 11]</td><td>[10,53]</td></tr><tr><td>Lowland tropical rainforest</td><td>Saut Lavilette (French Guiana)</td><td>4.151</td><td>-52.2</td><td>25.7</td><td>3590</td><td>224</td><td>Annonaceae, Sapotaceae</td><td>10</td><td>[3, 28]</td><td>[10,53]</td></tr><tr><td>Evergreen lowland tropical rainforest</td><td>Yasuní National Park (Ecuador)</td><td>0.683</td><td>-76.4</td><td>25</td><td>3129</td><td>251</td><td>Euphorbiaceae, Annonaceae, Fabaceae, Myrtaceae</td><td>59</td><td>[3, 21]</td><td>[9]</td></tr><tr><td>Lowland tropical rainforest</td><td>Biological Dynamics of Forest Fragments Project Reserve (BDFFP). Florestal and Cabo Frio (Central Amazon, Brazil)</td><td>-2.433</td><td>-59.83</td><td>27</td><td>2410</td><td>284</td><td>Fabaceae, Lecythidaceae, Sapotaceae</td><td>16</td><td>[3, 7]</td><td>[54] Oliveira, unpublished data</td></tr></tbody></table>

Pp: annual precipitation. SR: number of species with a DBH ≥10cm/ha.

doi:10.1371/journal.pone.0172495.t001

single ontogenetic stage: understory individuals). Overall, all criteria were consistently in direction of the underestimation of actual ITV.

Sampling permissions were granted by Delegación Provincial de la Conserjería de Agricultura y Medio Ambiente (Guadalajara-Castilla La Mancha) for the Mediterranean mixed forest in Alto Tajo Natural Park (Spain), by Instituto Boliviano de Investigación Forestal (IBIF) and the logging companies INPA Parket Ltd. and Planet La Chonta Investment Ltda for the tropical lowland dry deciduous forest in Bolivia, by Ministerio de Ambiente de Panama (MiAmbiente) for the tropical lowland semi-deciduous seasonal moist forest in Panama, by ICMBio/SISBIO- license number 23191-1 for the lowland tropical Atlantic Forest (Poço das Antas Biological Reserve-Southeastern Brazil), by Ministerio del Ambiente of Ecuador for the evergreen lowland tropical rainforest in Yasuní National Park (Ecuador) and by ICMBio/SISBIO- license number 18757-1 for the lowland tropical rainforest in Central Amazon (Brasil). For the rest of forests were not required specific permissions. The authorities responsible of these areas were informed and they expressed their consent to this sampling. Moreover, sampling did not involve endangered or protected species.

## Statistical analyses

We used the coefficient of variation (CV; 100 * standard deviation / mean) as an estimate of ITV. Because ITV may be influenced by the number of individuals sampled, we performed a rarefaction analysis in order to account for differences in sample size among species within and among forest communities [55]. This rarefaction analysis generated an expected trait value for each species in each forest by randomly drawing five individuals from the total pool of individuals of each species. We repeated this re-sampling process 1000 times for each species in each forest community. We then calculated the CV for each species from the average of the expected trait values generated by 1000 randomizations. To be sure of unbiased statistics estimated from rarefaction analysis, we checked both the community' rarefaction curves did not cross (S1 Fig) and also, species ranks in ITV were the same across sample sizes (S2 and S3
Figs)[55]. Moreover, we confirm that species with smaller sample sizes did not have systemati-
cally lower ITV values (S4 Fig for leaf size and S5 Fig for SLA). We performed a generalized linear
mixed model (GLMM, [56]) using ITV as response variable, species richness as an explanatory
variable together with the mean annual temperature and annual precipitation as explanatory
covariates in order to account for climate differences among forests. Type of forest was included
as a random factor to control for other intrinsic characteristics of each community.

We also assessed the degree of trait overlap among co-occurring species and its relationship
with species richness. Trait overlap is defined as the overlapping area between two trait distri-
bution curves [32] and calculated by (1) assuming that trait values of a species are normally
distributed around the mean [25,57] or (2) using kernel density estimators, which do not
assume any particular shape of the trait distribution [58]. Using a normal distribution rather
than kernel density tends to overestimate trait overlap, but at the same time it is considered
more robust to small sample size (i.e. in our case, species with 5 individuals) than kernel distri-
bution [58]. Because of these concerns, we estimated trait overlap with both normal and kernel
density approaches using the R function Trova [32]. For analyses with normal distributions, a
mean and standard deviation of the traits is required for each species in each forest commu-
nity. Given the differences among species in sampling intensity, we first ran a rarefaction anal-
ysis for each species by randomly re-sampling 5 individuals per species per site, repeated 1000
times. We then calculated the mean and standard deviation from the average of expected trait
values from 1000 randomizations for each species and forest community. Trait overlap figures
for both methods range from 0 to 1, where values close to 1 indicate a high overlap between
species or a high trait similarity. We calculated the median and the proportion of values
obtained with low (less than 0.25 out of 1) and high trait overlap (higher than 0.75 out of 1) as
a categorical description of the distribution of trait overlap values for each community from
both methods. Finally, we applied linear regression models using categorical parameter
description of the distribution of trait overlap values as response variables and species richness
as explanatory variable. Categorical parameters were square root transformed to improve
normality.

All statistical analyses were carried out in R v. 3.2.1 [59] using the packages lme4 [60] and MuMln [61].

**Results**

Intraspecific trait variability and species richness

We found considerable ITV for both leaf traits among species co-occurring in all forest com-
munities (Fig 2A and 2B). Accordingly, ITV for leaf size and SLA did not vary consistently
with species richness (Fig 2A and 2B; Table 2). None significant effects were also observed for
climatic covariates (mean annual temperature and annual precipitation) on the extent of ITV
for both studied traits (Table 2).

Trait overlap and species richness

For both leaf traits, median trait overlaps between species for each forest were significantly lower when assuming normal trait distribution than using kernel density estimators (Wilcoxon-Signed Rank test: n = 21; Z = 3.7; p-value < 0.001 for leaf size; n = 21, Z = 3.8, p-value < 0.001 for SLA). We found increasing trait overlap with species richness for both leaf traits using normal trait distribution (Fig 3A and 3B). We obtained similar results using kernel density estimators (r² = 0.64, p-value < 0.001 for leaf size, S6A Fig; and r² = 0.44, p-value < 0.001 for SLA, S6B Fig). We found that the proportion of species pairs with very low trait overlap (< 0.25) decreased significantly with species richness for both traits assuming normal trait distribution (Fig 3C and 3D) and
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Effect of species richness on the extent of ITV (estimated as coefficient of variation, CV) for (A) leaf size and (B) SLA. n = number of species measured for each forest community with ≥ 5 individuals measured. Empty circles indicate the average of ITV values for each forest community and richness level. The number of sampled species (n) may be both lower than the species richness of the community since we sub-sampled a forest layer, but also, potentially larger than the species richness of that community since the arrival of seeds from neighbor canopies may result in the establishment of new individuals that previously were not part of the main canopy of the focal community.</figcaption></figure>

doi:10.1371/journal.pone.0172495.g002
<table><caption>**Table 2.** Summary table of the GLMM testing the change in the extent of ITV with species richness (SR), mean annual temperature (Mean. Ann. T<sup>ra</sup>) and annual precipitation (Pp) for both leaf size and SLA.</caption><thead><tr><th>Response</th><th>Exp. variable</th><th>Estimate</th><th>SD. Error</th><th>t-value</th><th>p-value</th><th>R<sup>2</sup><sub>m</sub></th><th>R<sup>2</sup><sub>c</sub></th></tr></thead><tbody><tr><td rowspan="4">Leaf size n = 321</td><td>Intercept</td><td>5.098</td><td>0.379</td><td>13.436</td><td>0.000***</td><td rowspan="4">0.010</td><td rowspan="4">0.213</td></tr><tr><td>SR</td><td>0.027</td><td>0.004</td><td>0.706</td><td>0.480</td></tr><tr><td>Mean. Ann. T<sup>ra</sup></td><td>0.003</td><td>0.031</td><td>0.119</td><td>0.905</td></tr><tr><td>Pp</td><td>-0.000</td><td>0.000</td><td>-0.577</td><td>0.564</td></tr><tr><td rowspan="4">SLA n = 315</td><td>Intercept</td><td>3.726</td><td>0.302</td><td>12.338</td><td>0.000***</td><td rowspan="4">0.009</td><td rowspan="4">0.243</td></tr><tr><td>SR</td><td>0.001</td><td>0.003</td><td>0.469</td><td>0.639</td></tr><tr><td>Mean. Ann. T<sup>ra</sup></td><td>0.013</td><td>0.025</td><td>0.512</td><td>0.459</td></tr><tr><td>Pp</td><td>-0.000</td><td>0.000</td><td>-0.740</td><td>0.609</td></tr></tbody></table>

n = number of species included in the analyses; R<sup>2</sup><sub>m</sub>: marginal-R<sup>2</sup>; R<sup>2</sup><sub>c</sub>: conditional-R<sup>2</sup>.

(**): p-value <0.001

(**): p-value <0.01

(*): p-value < 0.05

(.): p-value <0.1.

doi:10.1371/journal.pone.0172495.t002

kernel density approximation (r<sup>2</sup> = 0.72, p-value < 0.001 for leaf size, S6C Fig and r<sup>2</sup> = 0.27, p-value < 0.001 for SLA, S6D Fig). Moreover, we found that the proportion of species pairs with high trait overlap (> 0.75) increased significantly with species richness for both traits assuming normal distribution (Fig 3E and 3F). However, under kernel approach, the relationship between the proportion of species pairs with high trait overlap with species richness was only marginally significant for leaf size (r<sup>2</sup> = 0.14, p-value = 0.07 for leaf size, S6E Fig) and none relationship was observed for SLA (S6F Fig).

Discussion

A largely unanswered question in biodiversity theory is whether ITV actually varies with species richness as predicted by classical niche theory [19,25] assuming all co-occurring species with equal fitness, that postulates that species would show narrower trait breadths (i.e. decrease of ITV) with increased species richness in order to avoid competition. Contrary to niche theory, our results showed a lack of relationship between ITV and species richness for leaf size and SLA, suggesting that the species' niches did not exhibit tight packing of the trait space in species-rich forests. In addition, we found greater trait overlap in species-rich communities for both traits, reflected here by an increase of the median values in trait overlap, with decreasing proportion of species pairs with low trait overlap and thus increasing proportion of species pairs with high trait overlap as species richness increased. Our results did not support the predictions from the principle of limiting similarity, which predicts a higher spread of trait values (i.e. trait dissimilarity among species) at the community level [42,62] since co-occurring individuals with high similarity in ecological requirements are more likely to face competitive exclusion [25].

The few studies testing the species richness-ITV-trait overlap relationships have so far focused on a single study system, plant growth form or functional trait, and found inconsistent results [35–37]. Some studies have shown declines in ITV and reduction of trait overlap with increasing species richness. For example, Hulshof *et al.* [35] found that the ratio between intraspecific and interspecific variability (a good proxy of trait overlap among species [19,32]) in SLA in woody plant communities decreased with increasing species richness. A similar pattern was also found by Kumordzi *et al.* [36] when studying variation in SLA of understory vegetation across different
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Linear regression models of the median values of trait overlap (panel A, B) and the proportion of low (less than 0.25; panel C, D) and high (> 0.75; panel E, F) values of degree of trait overlap between species for each forest community against species richness for both leaf size (left) and SLA (right). Trait similarity was calculated by assuming normal trait distribution of species with ≥ 5 individuals measured. (***): p-value <0.001; (**): p-value <0.01; (*): p-value < 0.05; (.): p-value <0.1.</figcaption></figure>

doi:10.1371/journal.pone.0172495.g003

boreal forest communities differing in species diversity. Felten *et al.* [63] found a decrease of niche breadth and niche overlap in temperate grasslands with increasing species richness, indicating complementarity of soil N use from different soil depth. In contrast, others have shown increasing ITV and trait overlap with increasing species richness, in concordance with our findings. For example, using a multi-trait approach (including SLA) in limestone grasslands, Le Bagousse-Pinguet *et al.* [37] found increases in both ITV and the ratio between intraspecific and
interspecific variability with increasing species richness, with no effect of environment on ITV in agreement with our results. They suggest asymmetric light competition among competing species as a potential explanation to this pattern. It may merely be due to the increased probability of having individual plants from a shade-intolerant species being slightly more tolerant than an individuals of a supposed shade-tolerant species; this situation would induce ITV and would minimize the differences in plant fitness and competitive ability among co-occurring species [36]. Finally, a recent global meta-analysis Siefert et al. [16] reports that the relative extent of ITV to the total community trait variance decreases with increasing species richness, but this pattern is mainly due to an increase in interspecific variance and consequently, in the total community variance, whereas the absolute extent of ITV remains fairly constant with species richness. While the analysis from Siefert et al. [16] supports partly our results, it includes different plant growth forms (i.e. both herbaceous and woody plants) from multiple community types (from grasslands to forests). Differences in ITV between plant forms can be expected since long-lived woody plants may present limited plasticity (i.e. less ITV) due to higher investment in longer lifespan tissues over their lifetimes compared to short-lived herbaceous species. Unlike, our study is based on a single growth form (freestanding woody plants in forest communities) on which debates about community assembly processes across the latitudinal gradient have been especially focused in recent years (e.g. [4,28,29,33,64,65]). It is important to note that a different way to calculate the ITV is used in our study in comparison with these aforementioned studies. Whereas our study measures ITV as the coefficient of variation at individual species level (i.e. absolute ITV at species level), others measure the mean intraspecific trait variance at community level (including ITV of all coexisting species in relation to total community trait variance). This difference in ITV measurement could explain in somehow differences in patterns found with respect to our results [16,35].

Our findings showing an increased functional similarity in hyperdiverse forests suggest, as Hubbell and Chave have argued (e.g. [28,29,66]), that it is more likely that individuals in high diversity forests are more functionally similar to each other than individuals in lower diversity forests. On the one hand, this can lead to a greater chance for neutral or nearly neutral dynamics in these more diverse communities [29]. On the other hand, if a higher degree of functional similarity between species pairs in diverse forests translates into smaller fitness differences between species (sensu Chesson [27]), only modest stabilizing niche differences between species (e.g. resource partitioning, density-dependent effects or population density fluctuations) would be required to drive community dynamics in a non-neutral fashion [38]. Unfortunately, these questions cannot be resolved without deeper understanding of how trait differences in woody plant communities relate to fitness and stabilizing niche differences. While recent experimental works have made these links for algal [67] and annual plant communities [68], considerable logistical barriers remain in long-lived plants as woody communities.

We suggest that our findings contradict niche theory and the principle of limiting similarity (promoting trait dissimilarity among species), but only based on previous assumptions of a flat fitness landscape (i.e. all species with equal fitness). However, an alternative conclusion could be reached assuming a multi-peak fitness landscapes [69]. Under multi-peak fitness landscapes, species on each peak may have been selected by the interplay of different processes among which competitive limit to similarity. In this scenario, species can reach equal fitness and maintain similar ITV [31] despite species richness increases (i.e. no decrease of ITV with species richness is expected).

In this study we explored shifts in the extent of ITV and trait overlap along a broad species richness gradient, but it is important to mention that our results may be in part limited since we did not carry out trait measurements on all individuals or species that were part of the whole community, particularly in species-rich communities. This likely may translates into an
underestimation of the extent of ITV and trait overlap between species in diverse forests. Nev-
ertheless, the high degree of species rarity in our tropical forests [70] makes it difficult to reach
a complete range for ITV. Moreover, it is important to take into account that we are using two
dimensions of the plant ecological strategy (SLA and leaf size) and they may not be good prox-
ies of plant species' realized niches in species-rich forests. In other words, they may not be cap-
turing niche differences among species, being the competitive ability for limiting resource use
determined by other key traits that we did not take into account. Further analyses incorporat-
ing other traits representative of different plant strategy axes, such as architecture traits and
woody density, and even a multi-trait approach could improve our understanding about the
traits that best relate to fitness and, thus, drive niche differentiation [70]. Furthermore, local
biotic and abiotic factors also have effects on the extent of ITV by selecting a particular subset
of trait values according to the local environment [71]. Future studies integrating other local
environmental factors, such as crowding, light availability or water availability [72–74], as well
as the environmental heterogeneity [75] would improve our understanding of the factors driv-
ing the relationship between ITV and species richness.

To conclude, our study highlights the key role that trait variability within species can play in
understanding community assembly along biodiversity gradients, and emphasizes the value of
estimating intraspecific variability in studies exploring trait diversity at the community level.
We found an increase of functional similarity among co-occurring species in more diverse com-
munities, to the widely recognized classical niche theory. Our study points to neutral processes
or equalizing mechanisms to explain that species with similar ecological requirements can be
present in the same community at the same time.

**Supporting information**

**S1 Dataset.** This file contains data belonging to the article "Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient" by Cristina C. Bastias, Claire Fortunel, Fernando Valladares, Christopher Baraloto, Raquel Benavides, William Cornwell, Lars Markesteijn, Alexandre A. de Oliveira, Jeronimo B.B Sansevero, Marcel C. Vaz, Nathan J. B. Kraft. The first sheet contains leaf trait data for each individual within each species and type of forest. The second sheet contains data from the rarefaction analyses: the rarefied mean, the rarefied coefficient of variation (CV; Fig 2) and the rarefied standard error (sd) for both leaf size (LS) and SLA at species level. Also, it contains the sample size (n) per each species and the species richness for each community (SR). The third sheet presents the overlap data assuming normal trait distribution for both leaf size (LS) and SLA: median, proportion of low (≤ 0.25) and high (> 0.75) values of degree of trait overlap between species per each type of forest (Fig 3). The last sheet presents the overlap data under a kernel approach for both leaf size (LS) and SLA: median, proportion of low (≤ 0.25) and high (> 0.75) values of degree of trait overlap between species per each type of forest (S6 Fig).
(XLSX)

**S1 Fig. Community accumulation curves at sample size of 5 individuals.** Dashed lines are 95% confident.
(TIF)

**S2 Fig. Species ranks at sample size of 10 individuals per species (ITV) vs species ranks at sample size of 5 individuals (rarefied ITV) in order to detect bias in the ITV estimate by a small sample size.** R² close to 1 means no bias (i.e. similar ITV values obtained for a species using 10 individuals and using 5 individuals).
(TIF)
**S3 Fig.** Species ranks at sample size of 20 individuals per species (ITV) vs species ranks at sample size of 5 individuals (rarefied ITV) in order to detect bias in the ITV estimate by a small sample size. R² close to 1 means no bias (i.e. similar ITV values obtained for a species using 20 individuals and using 5 individuals).

(TIF)

**S4 Fig.** Boxplot of the ITV values for leaf size grouped in 7 categories of sample size per species (No. of individuals). n: number of observations in each category.

(TIF)

**S5 Fig.** Boxplot of the ITV values for SLA grouped in 7 categories of sample size per species (No. of individuals). n: number of observations in each category.

(TIF)

**S6 Fig.** Linear regression models of the median values of trait overlap (panel A, B) and the proportion of low (less than 0.25; panel C, D) and high (> 0.75; panel E, F) values of degree of trait overlap between species for each forest community against species richness for both leaf size (left) and SLA (right). Trait similarity was calculated by kernel density approach using species with ≥ 5 individuals measured. (***): p-value <0.001; (**): p-value <0.01; (*): p-value < 0.05; (.): p-value <0.1.

(TIF)

**Acknowledgments**

We would like to thank all colleagues and technicians that assisted with field and laboratory work. This is the study 709 of the Technical Series of the Biological Dynamics of Forest Fragments Project.

**Author Contributions**

**Conceptualization:** CCB NK FV CF.

**Formal analysis:** CCB NK CF.

**Funding acquisition:** CCB CF FV CB RB WC LM AAdO JS MV NJBK.

**Investigation:** CCB CF FV CB RB WC LM AAdO JS MV NJBK.

**Methodology:** CCB CF NJBK.

**Project administration:** CCB CF FV NJBK.

**Supervision:** NJBK FV.

**Visualization:** CCB.

**Writing – original draft:** CCB CF NJBK.

**Writing – review & editing:** CCB CF FV CB RB WC LM AAdO JS MV NJBK.

**References**

1. Stevens GC. The Latitudinal Gradient in Geographical Range: How so Many Species Coexist in the Tropics. Am Nat. 1989; 133: 240–256.

2. Platnick NI. Patterns of biodiversity: tropical vs temperate. J Nat Hist. 1991; 25: 1083–1088.

3. Hill JL, Hill RA. Why are tropical rain forests so species rich? Classifying, reviewing and evaluating theories. Prog Phys Geogr. 2001; 25: 326–354.
4. Myers JA, Chase JM, Jiménez I, Jørgensen PM, Araujo-Murakami A, Paniagua-Zambrana N, et al. Beta-diversity in temperate and tropical forests reflects dissimilar mechanisms of community assembly. Ecol Lett. 2013; 16: 151–157. doi: 10.1111/ele.12021 PMID: 23113954

5. Ackerly DD, Cornwell WK. A trait-based approach to community assembly: partitioning of species trait values into within- and among-community components. Ecol Lett. 2007; 10: 135–145. doi: 10.1111/j. 1461-0248.2006.01006.x PMID: 17257101

6. Adler PB, Fajardo A, Kleinhesselink AR, Kraft NJB. Trait-based tests of coexistence mechanisms. Scherber C, editor. Ecol Lett. 2013; 16: 1294–1306. doi: 10.1111/ele.12157 PMID: 23910482

7. McGill B, Enquist B, Weiher E, Westoby M. Rebuilding community ecology from functional traits. Trends Ecol Evol. 2006; 21: 178–185. doi: 10.1016/j.tree.2006.02.002 PMID: 16701083

8. Weiher E, Freund D, Bunton T, Stefanski A, Lee T, Bentivenga S. Advances, challenges and a developing synthesis of ecological community assembly theory. Philos Trans R Soc Lond B Biol Sci. 2011; 366: 2403–2413. doi: 10.1098/rstb.2011.0056 PMID: 21768155

9. Kraft NJB, Valencia R, Ackerly DD. Functional traits and niche-based tree community assembly in an amazonian forest. Science. 2008; 322: 580–582. doi: 10.1126/science.1160662 PMID: 18948539

10. Fortunel C, Paine CET, Fine PVA, Kraft NJB, Baraloto C. Environmental factors predict community functional composition in Amazonian forests. J Ecol. 2014; 102: 145–155.

11. Garnier E, Laurent G, Bellmann A, Debain S, Berthelier P, Ducout B, et al. Consistency of species ranking based on functional leaf traits. New Phytol. 2001; 152: 69–83.

12. Laughlin DC, Fulé PZ, Huffman DW, Crouse J, Laliberté E. Climatic constraints on trait-based forest assembly. J Ecol. 2011; 99: 1489–1499.

13. Albert CH, Thuiller W, Yoccoz NG, Soudant A, Boucher F, Saccone P, et al. Intraspecific functional variability: extent, structure and sources of variation. J Ecol. 2010; 98: 604–613.

14. Paine CET, Baraloto C, Chave J, Hérault B. Functional traits of individual trees reveal ecological constraints on community assembly in tropical rain forests. Oikos. 2011; 120: 720–727.

15. Kraft NJB, Crutsinger GM, Forrestel EJ, Emery NC. Functional trait differences and the outcome of community assembly: an experimental test with vernal pool annual plants. Oikos. 2014; 123: 1391–1399.

16. Siefert A, Violle C, Chalmandrier L, Albert CH, Taudiere A, Fajardo A, et al. A global meta-analysis of the relative extent of intraspecific trait variation in plant communities. Chase J, editor. Ecol Lett. 2015; 18: 1406–1419. doi: 10.1111/ele.12508 PMID: 26415616

17. Lake JK, Ostling A. Comment on "Functional traits and niche-based tree community assembly in an Amazonian forest." Science. 2009; 324: 1015–1015.

18. Bolnick DI, Amarasekare P, Araújo MS, Bürger R, Levine JM, Novak M, et al. Why intraspecific trait variation matters in community ecology. Trends Ecol Evol. 2011; 26: 183–192. doi: 10.1016/j.tree.2011.01. 009 PMID: 21367482

19. Violle C, Enquist BJ, McGill BJ, Jiang L, Albert CH, Hulshof C, et al. The return of the variance: intraspecific variability in community ecology. Trends Ecol Evol. 2012; 27: 244–252. doi: 10.1016/j.tree.2011. 11.014 PMID: 22244797

20. 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. doi: 10.1111/j.1461-0248.2010.01476.x PMID: 20482582

21. Ainley LB, Vergés A, Bishop MJ. Congruence of intraspecific variability in leaf traits for two co-occurring estuarine angiosperms. Oecologia. 2016;

22. Jung V, Violle C, Mondy C, Hoffmann L, Muller S. Intraspecific variability and trait-based community assembly. J Ecol. 2010; 98: 1134–1140.

23. Whitham TG, Bailey JK, Schweitzer JA, Shuster SM, Bangert RK, LeRoy CJ, et al. A framework for community and ecosystem genetics: from genes to ecosystems. Nat Rev Genet. 2006; 7: 510–523. doi: 10.1038/nrg1877 PMID: 16778835

24. Crutsinger GM, Collins MD, Fordyce JA, Gompert Z, Nice CC, Sanders NJ. Plant Genotypic Diversity Predicts Community Structure and Governs an Ecosystem Process. Science. 2006; 313: 966–968. doi: 10.1126/science.1128326 PMID: 16917062

25. MacArthur R, Levins R. The limiting similarity, convergence, and divergence of coexisting species. Am Nat. 1967; 377–385.

26. Clark JS. Individuals and the Variation Needed for High Species Diversity in Forest Trees. Science. 2010; 327: 1129–1132. doi: 10.1126/science.1183506 PMID: 20185724

27. Chesson P. Mechanisms of Maintenance of Species Diversity. Annu Rev Ecol Syst. 2000; 31: 343–366.

28. Hubbell SP. Neutral theory in community ecology and the hypothesis of functional equivalence. Funct Ecol. 2005; 19: 166–172.
29. Hubbell S.P. The Unified Neutral Theory of Biodiversity and Biogeography. Princeton, NJ: Princeton University Press; 2001.

30. Adler PB, HilleRisLambers J, Levine JM. A niche for neutrality. Ecol Lett. 2007; 10: 95–104. doi: 10.1111/j.1461-0248.2006.00996.x PMID: 17257097

31. Svensson EI, Eroukhmanoff F, Friberg M. Effects of natural and sexual selection on adaptive population divergence and premating isolation in a damselfly. Evolution. 2006; 60: 1242–1253. PMID: 16892974

32. de Bello F, Carmona CP, Mason NWH, Sebastià M-T, Lepš J. Which trait dissimilarity for functional diversity: trait means or trait overlap? Zobel M, editor. J Veg Sci. 2013; 24: 807–819.

33. Lamanna C, Blonder B, Violle C, Kraft NJB, Sandel B, imova I, et al. Functional trait space and the latitudinal diversity gradient. Proc Natl Acad Sci. 2014; 111: 13745–13750. doi: 10.1073/pnas.1317722111 PMID: 25225365

34. Baraloto C, Timothy Paine CE, PatiÂ±o S, Bonal D, HÃault B, Chave J. Functional trait variation and sampling strategies in species-rich plant communities. Funct Ecol. 2010; 24: 208–216.

35. Hulshof CM, Violle C, Spasojevic MJ, McGill B, Damschen E, Harrison S, et al. Intra-specific and inter-specific variation in specific leaf area reveal the importance of abiotic and biotic drivers of species diversity across elevation and latitude. J Veg Sci. 2013; 24: 921–931.

36. Kumordzi BB, de Bello F, Freschet GT, Le Bagousse-Pinguet Y, Lepš J, Wardle DA. Linkage of plant trait space to successional age and species richness in boreal forest understorey vegetation. Gibson D, editor. J Ecol. 2015; 103: 1610–1620.

37. Le Bagousse-Pinguet Y, de Bello F, Vandewalle M, Leps J, Sykes MT. Species richness of limestone grasslands increases with trait overlap: evidence from within- and between-species functional diversity partitioning. J Ecol. 2014; 102: 466–474.

38. Adler PB, HilleRisLambers J, Levine JM. A niche for neutrality. Ecol Lett. 2007; 10: 95–104. doi: 10.1111/j.1461-0248.2006.00996.x PMID: 17257097

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

40. Fonseca CR, Overton JM, Collins B, Westoby M. Shifts in trait-combinations along rainfall and phosphorus gradients. J Ecol. 2000; 88: 964–977.

41. Ackerly D, Knight C, Weiss S, Barton K, Starmer K. Leaf size, specific leaf area and microhabitat distribution of chaparral woody plants: contrasting patterns in species level and community level analyses. Oecologia. 2002; 130: 449–457.

42. Cornwell WK, Ackerly DD. Community assembly and shifts in plant trait distributions across an environmental gradient in coastal California. Ecol Monogr. 2009; 79: 109–126.

43. Westoby M. A leaf-height-seed (LHS) plant ecology strategy scheme. Plant Soil. 1998; 199: 213–227.

44. Miller PC, Stoner WA. Canopy Structure and Environmental Interactions. In: Solbrig OT, Jain S, Johnson GB, Raven PH, editors. Topics in Plant Population Biology. Macmillan Education UK; 1979. pp. 428–458. Available: http://link.springer.com/chapter/10.1007/978-1-349-04627-0_19

45. Avery Thomas Eugene, Burkhart Harold. Forest Measurements. In: McGraw Hill Education [Internet]. 2002 [cited 29 Jun 2016]. Available: http://www.mheducation.com/highered/product.0073661767.html

46. Ricklefs RE, He F. Region effects influence local tree species diversity. Proc Natl Acad Sci. 2016; 113: 674–679. doi: 10.1073/pnas.1523683113 PMID: 26733680

47. Yeboah D, Chen HYH, Kingston S. Tree species richness decreases while species evenness increases with disturbance frequency in a natural boreal forest landscape. Ecol Evol. 2016; 6: 842–850. doi: 10.1002/ece3.1944 PMID: 26865971

48. Markesteijn L, Poorter L, Bongers F. Light-dependent leaf trait variation in 43 tropical dry forest tree species. Am J Bot. 2007; 94: 515–525. doi: 10.3732/ajb.94.4.515 PMID: 21636421

49. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005; 25: 1965–1978.

50. Albert CH, Grassein F, Schurr FM, Vieilledent G, Violle C. When and how should intraspecific variability be considered in trait-based plant ecology? Perspect Plant Ecol Evol Syst. 2011; 13: 217–225.

51. Umaña MN, Zhang C, Cao M, Lin L, Swenson NG. Commonness, rarity, and intraspecific variation in traits and performance in tropical tree seedlings. Suding K, editor. Ecol Lett. 2015; 18: 1329–1337. doi: 10.1111/ele.12527 PMID: 26415689

52. Mouillot D, Bellwood DR, Baraloto C, Chave J, Galzin R, Harmelin-Vivien M, et al. Rare Species Support Vulnerable Functions in High-Diversity Ecosystems. Mace GM, editor. PLoS Biol. 2013; 11: e1001569. doi: 10.1371/journal.pbio.1001569 PMID: 23723735
53. Baraloto C, Hardy OJ, Paine CET, Dexter KG, Cruaud C, Dunning LT, et al. Using functional traits and phylogenetic trees to examine the assembly of tropical tree communities: *Assembly of tropical tree communities*. J Ecol. 2012; 100: 690-701.

54. Vaz, Marcel Carita. Disentangling the mechanisms that shape community assembly across diversity gradients is a fundamental question in ecology. Master's Dissertation. 2011.

55. Gotelli NJ, Colwell RK. Estimating species richness. Biol Divers Front Meas Assess. 2011; 12: 39-54.

56. Breslow NE, Clayton DG. Approximate Inference in Generalized Linear Mixed Models. J Am Stat Assoc. 1993; 88: 9-25.

57. Leps J, De Bello F, Lavorel S, Berman S. Quantifying and interpreting functional diversity of natural communities: practical considerations matter. Preslia. 2006; 78: 481-501.

58. Mouillot D, Stubbs W, Faure M, Dumay O, Tomasini JA, Wilson JB, et al. Niche overlap estimates based on quantitative functional traits: a new family of non-parametric indices. Oecologia. 2005; 145: 345-353. doi: 10.1007/s00442-005-0151-z PMID: 16001221

59. R Development Core Team. R: A language and environment for statistical computing. [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2013. Available: http://www.R-project.org.

60. Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, et al. Package “lme4.” R Found Stat Comput Vienna. 2014; Available: http://cran.r-mirror.de/web/packages/lme4/lme4.pdf

61. Barton K. MuMln: Multi-model inference. R package version 1.15.6. Available on-line at http://CRAN.R-project.org/package=MuMln. 2015;

62. De Bello F, Thuiller W, Lepš J, Choler P, Clément J-C, Macek P, et al. Partitioning of functional diversity reveals the scale and extent of trait convergence and divergence. J Veg Sci. 2009; 20: 475-486.

63. Felten S von, Hector A, Buchmann N, Niklaus PA, Schmid B, Scherer-Lorenzen M. Belowground Nitrogenc Partitioning in Experimental Grassland Plant Communities of Varying Species Richness. Ecology. 2009; 90: 1389-1399. PMID: 19537558

64. HilleRisLambers JHR, Clark JS, Beckage B. Density-dependent mortality and the latitudinal gradient in species diversity. Nature. 2002; 417: 732-735. doi: 10.1038/nature00809 PMID: 12066182

65. Kraft NJB, Comita LS, Chase JM, Sanders NJ, Swenson NG, Crist TO, et al. Disentangling the Drivers of β Diversity Along Latitudinal and Elevational Gradients. Science. 2011; 333: 1755-1758. doi: 10.1126/science.1208584 PMID: 21940897

66. Chave J. Neutral theory and community ecology: Neutral theory and community ecology. Ecol Lett. 2004; 7: 241-253.

67. Narwani A, Alexandrou MA, Oakley TH, Carroll IT, Cardinale BJ. Experimental evidence that evolutionary relatedness does not affect the ecological mechanisms of coexistence in freshwater green algae. De Meester L, editor. Ecol Lett. 2013; 16: 1373-1381. doi: 10.1111/ele.12182 PMID: 24112458

68. Kraft NJB, Godoy O, Levine JM. Plant functional traits and the multidimensional nature of species coexistence. Proc Natl Acad Sci. 2015; 112: 797-802. doi: 10.1073/pnas.1413650112 PMID: 25561561

69. Gavrilets S. Fitness landscapes and the origin of species. Princeton Univ. Press, Princeton, NJ.; 2004.

70. Spasojevic MJ, Suding KN. Inferring community assembly mechanisms from functional diversity patterns: the importance of multiple assembly processes: *Functional diversity along gradients*. J Ecol. 2012; 100: 652-661.

71. Diaz S, Cabido M, Casanoves F. Plant functional traits and environmental filters at a regional scale. J Veg Sci. 1998; 9: 113-122.

72. Whitmore T.C. A review of some aspects of tropical rain forest seedling ecology with suggestions for further enquiry. UNESCO, Paris, France, and also by Parthenon Publishing Group Inc., New York, USA; 1996.

73. Ackerly D. Functional Strategies of Chaparral Shrubs in Relation to Seasonal Water Deficit and Disturbance. Ecol Monogr. 2004; 74: 25-44.

74. Le Bagousse-Pinguet Y, Börger L, Quero J-L, García-Gómez M, Soriano S, Maestre FT, et al. Traits of neighbouring plants and space limitation determine intraspecific trait variability in semi-arid shrublands. J Ecol. 2015; 103: 1647-1657.

75. Galloway L.F. Response to natural environmental heterogeneity: maternal effects and selection on life-history characters and plasticities in Mimulus guttatus. Evolution. 1995 49: 1095-1107.