
==== Front
Plant Divers
Plant Divers
Plant Diversity
2096-2703
2468-2659
Kunming Institute of Botany, Chinese Academy of Sciences

S2468-2659(24)00038-6
10.1016/j.pld.2024.03.005
Research Paper
Woody plant taxonomic, functional, and phylogenetic diversity decrease along elevational gradients in Andean tropical montane forests: Environmental filtering and arrival of temperate taxa
Bañares-de-Dios Guillermo gbanaresdd@gmail.com
a⁎
Macía Manuel J. bc
Arellano Gabriel de
Granzow-de la Cerda Íñigo fa
Vega-Álvarez Julia g
Arnelas Itziar hi
Espinosa Carlos I. i
Salinas Norma jk
Cayuela Luis al
a Departamento de Biología y Geología, Física y Química Inorgánica, Universidad Rey Juan Carlos, Calle Tulipán s/n, Móstoles, Madrid, ES-28933, Spain
b Departamento de Biología, Área de Botánica, Universidad Autónoma de Madrid, Calle Darwin 2, Madrid, ES-28049, Spain
c Centro de Investigación en Biodiversidad y Cambio Global (CIBC-UAM), Universidad Autónoma de Madrid, Calle Darwin 2, Madrid, ES-28049, Spain
d Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI 48109, USA
e Oikobit LLC, www.oikobit.com, Albuquerque, NM, 87120, USA
f Real Jardín Botánico – CSIC, Plaza Murillo 2, Madrid, ES-28014, Spain
g Departamento de Biología Animal, Ecología, Parasitología, Edafología y Química Agrícola, Universidad de Salamanca, Calle Donantes de Sangre s/n, Salamanca, ES-37007, Spain
h Departamento de Biodiversidad, Ecología y Evolución, Universidad Complutense de Madrid, Calle José Antonio Nováis 12, Ciudad Universitaria, ES-28040, Spain
i Departamento de Ciencias Biológicas y Agropecuarias, Universidad Técnica Particular de Loja, A. P. 11-01-608 Loja, Ecuador
j Sección Química, Pontificia Universidad Católica del Perú, A. P. Lima 32, Lima, Peru
k ECI, School of Geography and Environment, University of Oxford, OX1 3QY, Oxfordshire, United Kingdom
l Instituto de Investigación en Cambio Global (IICG-URJC), Universidad Rey Juan Carlos, Calle Tulipán s/n, Móstoles, Madrid, ES- 28933, Spain
⁎ Corresponding author. gbanaresdd@gmail.com
25 3 2024
7 2024
25 3 2024
46 4 491501
8 10 2023
15 3 2024
18 3 2024
© 2024 Kunming Institute of Botany, Chinese Academy of Sciences. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co., Ltd.
2024
Kunming Institute of Botany, Chinese Academy of Sciences
https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Mountains are paramount for exploring biodiversity patterns due to the mosaic of topographies and climates encompassed over short distances. Biodiversity research has traditionally focused on taxonomic diversity when investigating changes along elevational gradients, but other facets should be considered. For first time, we simultaneously assessed elevational trends in taxonomic, functional, and phylogenetic diversity of woody plants in Andean tropical montane forests and explored their underlying ecological and evolutionary causes. This investigation covered four transects (traversing ca. 2200 m a.s.l.) encompassing 114 plots of 0.1 ha across a broad latitudinal range (ca. 10°). Using Hill numbers to quantify abundance-based diversity among 37,869 individuals we observed a consistent decrease in taxonomic, functional, and phylogenetic diversity as elevation increased, although the decrease was less pronounced for higher Hill orders. The exception was a slight increase in phylogenetic diversity when dominant species were over-weighted. The decrease in taxonomic and functional diversity might be attributed to an environmental filtering process towards highlands, where the increasingly harsher conditions exclude species and functional strategies. Besides, the differences in steepness decrease between Hill orders suggest that rare species disproportionately contribute to functional diversity. For phylogenetic diversity the shifting elevational trend between Hill orders indicates a greater than previously considered influence in central Andean highlands of tropical lowlands originated species with strong niche conservatism relative to distantly related temperate lineages. This could be explained by a decreasing presence and abundance of temperate, extratropical taxa towards the central Andes relative to northern or southern Andes, where they are more prevalent.

Highlights

• Patterns of taxonomic, functional and phylogenetic diversity in tropical woody plants were assessed for the first time across elevational gradients in Andean tropical montane forests.

• Taxonomic, functional, and phylogenetic diversity decrease with elevation with only the exception of a slight increase in phylogenetic diversity when dominant species are over-weighted.

• The decrease in taxonomic and functional diversity might be the result of environmental filtering process.

• The increase in phylogenetic diversity only when focus is given to dominant species might be explained by a lesser presence of extratropical taxa towards the central Andes relative to northern or southern Andes.

Keywords

Altitudinal gradients
Andes
Environmental filtering
Out of the Tropical Lowlands
Functional rarity
Tropical Niche Conservatism
==== Body
pmc1 Introduction

Biodiversity distribution patterns and their causes have been central to natural history and ecology since they were first conceived (Caldas, 1803; Humboldt and Bonpland, 1805) until the present (Storch et al., 2005; Sanmartín, 2012). Mountains provide invaluable insight for developing a general theory on biological diversity to explain these patterns (Körner, 2004; Sanders and Rahbek, 2012). Compared to latitudinal gradients, elevational gradients encompass more climate regimes, greater topographic complexity, and wider environmental oscillations over shorter geographical distances (Antonelli et al., 2018; Rahbek et al., 2019a). Thus, mountains are exceptionally biodiverse: where they comprise only 25% of Earth's terrestrial surface but harbour ca. 87% of terrestrial biodiversity, with a large proportion of endemicity (Rahbek et al., 2019b; Wambulwa et al., 2021). In particular, Andean tropical montane forests (TMFs) host one-sixth of all plant species, where 44% of which are endemic (Myers et al., 2000; Hughes and Eastwood, 2006). Therefore, Andean TMFs are among the biodiversity hotspots with highest conservation priority (Mittermeier et al., 2005). Lastly, from a biogeographic perspective, Andean TMFs are unique natural laboratories that provide rigorous tests of diversity trends along altitude gradients and comparisons of these trends among sites (Lomolino, 2001; Guo et al., 2013; Graham et al., 2014) because they allow to replicate ample elevational ranges over a large continental north-south corridor that has witnessed historic processes over millions of years (Antonelli et al., 2009; Malhi et al., 2010; Tito et al., 2020).

Biodiversity has traditionally been explored in terms of taxonomic diversity by treating all Linnaean binomials as ecologically equivalent and evolutionary independent. This approach cannot yield much information about the ultimate drivers of biodiversity distribution patterns (Pavoine and Bonsall, 2011). A more accurate depiction of diversity should also consider functional and phylogenetic diversity (McGill et al., 2006; Swenson, 2011). Functional diversity uses species ecological, morphological, and physiological features to understand where species can thrive, how they interact, or how they contribute to ecosystem functioning (Petchey and Gaston, 2006; Violle et al., 2007; Cadotte et al., 2011). Phylogenetic diversity tracks the origins of taxa lineages through deep historical and evolutionary processes that shaped their distribution patterns and thereby provides insights into the key events in biodiversity distribution (Cavender-Bares et al., 2009; Gerhold et al., 2018). Although these three diversity facets are intrinsically related (e.g., closely related taxa tend to share more traits; Moles et al., 2005) their responses to the environment can differ (Devictor et al., 2010; Cisneros et al., 2014). However, few studies have simultaneously considered the three types of diversity, particularly for plants under changing environmental conditions, such as those found along elevational gradients (but see Tanaka and Sato, 2015; Chun and Lee, 2017). Ignoring one diversity facet not only restricts our understanding of the causes of biodiversity distribution but also hinders conservation strategies, as focusing only in one dimension may lead to the neglect of unique species, functional roles and/or evolutionary histories (Brum et al., 2017).

Previous studies on biodiversity patterns along elevational gradients in different biomes and distinct taxonomic groups have arrived at contrasting results regarding taxonomic diversity, but largely agree on functional and phylogenetic trends. Historically, the prevalent pattern for taxonomic diversity is a decrease with elevation, although some studies have described species richness peaks at medium elevations (Rahbek, 1995; McCain, 2005; Dani et al., 2023). The most common pattern for functional diversity also comprises narrowing towards fewer functional strategies to cope with harsher environmental conditions as elevation increases (Swenson et al., 2012; de Bello et al., 2013; Thakur and Chawla, 2019). The elevational pattern for phylogenetic diversity is specific to the key biogeographic historical events of each site (Ricklefs and Schluter, 1993). In the particular case of the Andes, two hypotheses have been proposed to explain the historical formation of biological communities, reflecting distinct patterns of phylogenetic diversity: a) the establishment of Amazonian clades originated under warm conditions at tropical lowlands –typically exhibiting strong conservatism of ancestral narrow thermal tolerances– which, during the Pleistocene, adapted to the comparatively cooler temperatures of the newly raised highlands resulting from the Andean uplift (Latham and Ricklefs, 1993; Kerkhoff et al., 2014); and b) the arrival of cold pre-adapted lineages –originated much earlier in temperate latitudes– via migration along the Andean high-elevation corridor (Jablonski et al., 2006). The former (a) is the cornerstone of the Tropical Niche Conservatism (TNC, Wiens and Donoghue, 2004) hypothesis, which predicts an elevational decrease in phylogenetic diversity (Tolmos et al., 2022), whereas the latter (b) is the core of the Out of the Tropical Lowlands (OTL, Qian and Ricklefs, 2016) hypothesis, which predicts an elevational increase in phylogenetic diversity (González-Caro et al., 2020).

In this study, we investigated plant diversity in Andean TMFs at four transects ranging elevational gradients from ca. 800–3100 m a.s.l. and spanning ca. 1750 km in latitude across Ecuador, Peru, and Bolivia. We used census data from 114 plots (containing ca. 38,000 individuals and 2000 taxa) and analysed them using Hill numbers to quantify different types of abundance-based diversity indices. We aimed to 1) assess the elevational trends in taxonomic, functional, and phylogenetic diversity, 2) examine the relationships between different diversity facets along elevation and 3) explore their underlying ecological and evolutionary causing processes. To the best of our knowledge, this is one of the most ambitious elevational studies of the three facets of plant diversity and a similar study has yet to be conducted along the Andean range.

2 Materials and methods

2.1 Study area

We studied four transects ranging elevational gradients of Andean TMFs located in remote regions: Podocarpus National Park (Ecuador), Río Abiseo National Park (Peru), Madidi National Park (Bolivia), and Pilón-Lajas Biosphere Reserve (Bolivia) (Fig. 1). All four transects contained a similar elevational range (ca. 2200 m), within which we defined three elevational belts: premontane (900–1300 m), montane (1800–2250 m), and upper (2700–3100 m). Within each belt, we established nine to ten plots of 0.1 ha (50 m × 20 m) as described by Arellano et al. (2016), with a total of 114 plots (Tables 1 and S1). Plots were located in mature and apparently undisturbed forests at distances of at least 300 m. In each plot we inventoried all woody stems (including trees, shrubs, lianas, hemiepiphytes, palms, and tree ferns, but excluding woody Poaceae) with a diameter at breast height ≥2.5 cm, with a total of 37,869 individuals belonging to 2021 taxa (71% individuals identified at species level, 18% at genera and 11% at family) (Table 1, Appendix A).Fig. 1 Locations of the four transects ranging elevational gradients in the tropical montane forests (TMFs) studied along the central Andes (blue dots). Census data from n = 114 plots of 0.1 ha located in remote regions: Podocarpus National Park (Ecuador), Río Abiseo National Park (Peru), Madidi National Park (Bolivia), and Pilón-Lajas Biosphere Reserve (Bolivia).

Fig. 1

Table 1 Summary of elevational profiles and composition of 114 Andean TMFs plots located in four transects ranging elevational gradients.n_inds.: number of individual stems; n_taxa: number of taxa. Specific detail for each plot is available at Table S1. All four transects contained a similar elevational range, within which we defined three elevational belts: premontane, montane, and upper.

Table 1Transect (n_plots)	Elevational profile	Composition	
Belt	Range (m a.s.l.)	n_inds.	n_taxa	
Podocarpus
-Ecuador-
(30)	Premontane	1033–1249	3.153	452	
Montane	1806–2217	2.917	389	
Upper	2674–2906	3.840	187	
Gradient extent: 1873 m	1033–2906	9.910	832	
Río Abiseo
-Peru-
(30)	Premontane	745–1122	2.770	355	
Montane	2093–2233	3.339	230	
Upper	2721–2980	3.099	148	
Gradient extent: 2235 m	745–2980	9.208	585	
Madidi
-Bolivia-
(27)	Premontane	1224–1531	3.111	269	
Montane	2019–2303	2.219	205	
Upper	2801–3106	4.832	88	
Gradient extent: 1882 m	1224–3106	10.162	466	
Pilón-Lajas
-Bolivia-
(27)	Premontane	1166–1496	2.892	275	
Montane	2139–2281	2.751	156	
Upper	2789–3061	2.944	119	
Gradient extent: 1895 m	1166–3061	8.587	462	

Environmental characterization of the plots was conducted using bioclimatic variables from the CHELSA climatological dataset (Karger et al., 2017) (Table S2). Mean annual temperature is the most representative variable for studies along elevational gradients because it represents the fundamental environmental change associated with elevation (Körner, 2007). In our study sites, mean annual temperature exhibits high correlation with elevation (r > 0.98) and the thermal ranges of transects vary slightly: 12 °C in Podocarpus, 9 °C in Río Abiseo, and 10 °C in both Madidi and Pilón-Lajas. Alternative climatic variables related to precipitation are not as relevant for our study because 1) moist TMFs, characterized by persistent fog and high precipitation throughout most of the year, are not subjected to relevant seasonal water stress, and 2) at each transect, plots were located along the same river basin to minimize spatial rainfall variability inherent to the complex topography of the Andes, thus leaving temperature as the key source of abiotic variation along our elevational gradients (Laraque et al., 2007; Malhi et al., 2010). Small changes in mean annual temperature are expected to have a major influence on species distributions, especially in the tropics where species have narrow thermal tolerances (Janzen, 1967).

2.2 Taxonomic and functional characterization

We collected voucher specimens in the field, which were later identified and deposited in local herbaria (HUTPL and LPB [Thiers, 2022]). Taxonomic standardization was conducted using the ‘Taxonstand’ R package (Cayuela et al., 2012; v.2.2). A summary of taxonomic composition for elevational belts can be found in Table 1 (specific details for each plot are available at Table S1).

Functional characterization was only conducted for two transects: Podocarpus (Ecuador) and Río Abiseo (Peru). For each taxon, we measured the specific leaf area, leaf thickness, and branch wood density as described by Cornelissen et al. (2003). These functional traits have been widely used for assessing the key woody plant functional strategy axes (Chave et al., 2009; Wright et al., 2004). We calculated mean trait values for each taxon (Báez et al., 2022). On average for each plot, functional data for the three traits were available for 65 taxa (86%) in Podocarpus and 63 taxa (92%) in Río Abiseo.

2.3 Phylogenetic tree generation

We generated a phylogenetic tree (Appendix A) using the ‘V.PhyloMaker2’ R package (Jin and Qian, 2022), which allowed us to prune our species list from the largest currently available dated mega-phylogeny for vascular plants.

For phylogenetic diversity calculations, we only considered taxa identified to the species or genus level, i.e., 703 taxa (88%) in Podocarpus, 416 taxa (71%) in Río Abiseo, 453 taxa (97%) in Madidi, and 457 taxa (98%) in Pilón-Lajas. To include taxa absent from the original mega-phylogeny in the tree, we joined them at the half-way point of the genus/family branch (for species/genera) using the parameter ‘scenarios’ = S3 in the ‘phylo.maker’ function.

2.4 Diversity measurements

Attribute diversity comprises a unified framework for assessing biodiversity by integrating three diversity facets: taxonomic, functional, and phylogenetic diversity (Table S3) (Chao et al., 2014a). Each type of diversity is measured in different entities (species for taxonomic, species pairwise Euclidean distances for functional, and species tree branch length for phylogenetic), but all three can be transformed into directly comparable units by using generalized Hill numbers. Hill numbers (Hill, 1973) quantify abundance-based taxonomic diversity using the parameter, q, which defines the sensitivity to species abundances (Jost, 2006; Chao et al., 2014b). For q = 0, all taxa are treated equally regardless of their abundance (i.e., taxa richness), whereas larger q values give greater importance to the relative abundances of taxa (i.e., more weight to abundant/dominant taxa and less weight to rare taxa). The attribute diversity framework extends the usage of Hill numbers to functional (Chiu and Chao, 2014) and phylogenetic (Chao et al., 2010) diversity by weighting each entity based on either its relative functional distance compared with other taxa (functional diversity) or the relative length of the phylogenetic tree branch for that taxon (phylogenetic diversity). We calculated the attribute diversity indexes for Hill numbers of the orders q = 0, q = 1 and q = 2 at the plot level using the ‘renyi’ function of the ‘vegan’ R package (Oksanen et al., 2019; v.2.5–6) for taxonomic diversity, the code from Chiu and Chao (2014) for functional diversity and the ‘ChaoPD’ function of the ‘entropart’ R package (Chao et al., 2010; v.1.6–10) for phylogenetic diversity.

2.5 Statistical analyses

First, to acknowledge the influence of sample size on diversity metrics we conducted Pearson's correlations between number of trees per plot and the three facets of diversity (taxonomic, functional, and phylogenetic) for different Hill orders.

Next, we explored Pearson's correlations among the three types of diversity for Hill numbers of different orders.

Then, we used generalized linear models (GLMs) to investigate the influence of elevation on the three facets of diversity for these Hill orders. For each diversity facet and Hill order, we fitted five models that included all possible combinations of effects: 1) no effect (i.e., null model), 2) only elevation, 3) only transect, 4) similar effect of elevation across transects (i.e., without interaction between elevation and transect or parallel slopes model), and 5) differential effect of elevation among transects (i.e., with interaction between elevation and transect or different slopes model). We used a negative binomial error distribution (and a log-link function) to account for overdispersion. Alternative models were compared based on the corrected Akaike information criterion (AICc) and we considered the models within two AICc units of the model with the lowest AICc value. Among these models, we selected the most complex and calculated the explained deviance (D2) as a measure of goodness of fit for GLMs. Analyses were conducted using ‘MASS’ (Venables and Ripley, 2002; v.7.3–51.6) and ‘MuMIn’ (Barton, 2013; v.1.43.17) R packages.

All the R Statistical Software analyses were conducted using v.4.0.2 (R Core Team, 2021).

3 Results

Pearson's correlations coefficients between number of trees per plot and any of the three facets of diversity (taxonomic, functional, and phylogenetic) for different Hill orders was very low, ranging between −0.22 and 0.08 and hence there was no need to consider the effect of sample size on diversity metrics.

When all transects were analysed together, taxonomic, and functional diversity were always significantly and positively correlated for q = 0, q = 1, and q = 2 (Fig. 2 and Table S4). However, taxonomic, and phylogenetic diversity, and functional and phylogenetic diversity were significantly and positively correlated only for q = 0 and q = 1. After analysing the transects independently, we found the following.• Taxonomic vs. functional diversity: always significantly and positively correlated for q = 0, q = 1, and q = 2.

• Taxonomic vs. phylogenetic diversity: significantly and positively correlated for q = 0 and q = 1 (except in Podocarpus for q = 1).

• Functional and phylogenetic diversity: significantly and positively correlated for q = 0 and q = 1 (except in Podocarpus for q = 1).

Fig. 2 Pairwise comparisons of the three types of diversity (taxonomic, functional, and phylogenetic) for each of the four transects ranging elevational gradients in Andean tropical montane forests. Predictions of generalized linear models (GLMs) relating a) functional versus taxonomic diversity, b) phylogenetic versus taxonomic diversity, and c) functional versus phylogenetic diversity patterns for Hill orders of q = 0 (left-hand charts), q = 1 (central charts), and q = 2 (right-hand charts). Numbers represent the taxonomic, functional, and phylogenetic entities, a non-dimensional unit result of using the attribute diversity framework. Shaded areas represent 95% confidence intervals based on the model predictions. GLMs were fitted using a negative binomial error distribution and a log-link function. Note that the scale of some axes differ so that charts would be visually informative.

Fig. 2

Overall, results were consistent between transects. The three types of diversity values decreased from lower to higher elevations for q = 0, q = 1, and q = 2. This trend was strongest for q = 0 and weakest for q = 2, as indicated by both the slope of the regression curve (Fig. 3) and the explained deviance (Table 2). The only exception was phylogenetic diversity for q = 2, which appeared to increase slightly with the elevation. The best fitted models always included both elevation and transect (with or without interaction) as predictors, except for phylogenetic diversity for q = 2 (Fig. 3; Table 2, Appendix A).Fig. 3 Patterns of taxonomic, functional, and phylogenetic diversity in Andean tropical montane forests along the four transects ranging elevational gradients. Predictions of generalized linear models (GLMs) relating a) taxonomic, b) functional, and c) phylogenetic diversity patterns with elevation (m a.s.l.) for Hill orders of q = 0 (left-hand charts), q = 1 (central charts), and q = 2 (right-hand charts). Numbers represent the taxonomic, functional, and phylogenetic entities, a non-dimensional unit result of using the attribute diversity framework. Shaded area represents 95% confidence intervals based on the model predictions. GLMs were fitted using a negative binomial error distribution and a log-link function. Note that the scale of some axes differ so that charts would be visually informative.

Fig. 3

Table 2 Corrected Akaike information criterion (AICc) values for generalized linear models (GLMs) of the taxonomic, functional, and phylogenetic diversity values calculated for Hill numbers of different order relative to distinct predictors for 114 plots of Andean TMFs. Values for the best fitted models are highlighted in bold. The explained deviance D2 was calculated only for the best and most complex model (if more than one model was the best, the differences in terms of deviance between the most and least complex would be minimal).

Table 2Taxonomic diversity	AICc	D2	
	q = 0	q = 1	q = 2	q = 0	q = 1	q = 2	
–	1059.49	988.25	908.37	-	-	-	
Elevation	959.23	895.57	830.82	-	-	-	
Transect	1048.93	973.07	894.85	-	-	-	
Elevation + Transect	939.53	868.04	810.33	0.68	-	-	
Elevation: Transect	942.22	864.05	808.73	-	0.68	0.61	
Functional diversity	q = 0	q = 1	q = 2	q = 0	q = 1	q = 2	
–	1297.02	1204.95	1119.08	-	-	-	
Elevation	1262.95	1164.07	1084.03	-	-	-	
Transect	1286.12	1196.45	1111.83	-	-	-	
Elevation + Transect	1251.23	1157.86	1079.89	-	0.55	0.48	
Elevation: Transect	1249.44	1162.64	1083.67	0.58	–	–	
Phylogenetic diversity	q = 0	q = 1	q = 2	q = 0	q = 1	q = 2	
–	1919.31	1489.19	1303.35	-	-	-	
Elevation	1823.22	1475.77	1304.43	-	-	-	
Transect	1917.71	1475.89	1294.10	-	-	-	
Elevation + Transect	1810.94	1460.23	1295.25	-	-	-	
Elevation: Transect	1811.94	1460.02	1295.56	0.65	0.28	0.14	

4 Discussion

4.1 Relationships between the three facets of diversity along elevation

Taxonomic, functional, and phylogenetic diversity responded to elevation in a similar manner (Fig. 3), and they were significantly positively correlated in most cases (Fig. 2 and Table S4). These positive correlations may be expected due to changes in species richness along elevational gradients, where more species are likely to exhibit more functional strategies and belong to more phylogenetic lineages (Losos, 2008). However, this is not always necessarily true and empirical studies often show that two communities with equal number of species can differ greatly in terms of number of strategies and/or lineages (Devictor et al., 2010; Safi et al., 2011) and that a community exhibiting similar strategies may foster distantly related lineages (trait convergence) (Freckleton and Jetz, 2009). Few studies have quantified the three types of diversity simultaneously along elevational gradients for plant communities (e.g., Tanaka and Sato, 2015; Chun and Lee, 2017). In the tropical Andes, this type of study has only been conducted for vertebrates (e.g., Dehling et al., 2014 for birds; Cisneros et al., 2014 for bats; and Dreiss et al., 2015 for rodents). Congruence between the three types of diversity along elevational gradients has not been reported for plants, although pairwise correlations between two facets of diversity have been found: e.g., functional and phylogenetic diversity (Tanaka and Sato, 2015; Chun and Lee, 2017). Regardless of the existence or absence of correlations between the different facets of diversity, all three provide complementary insights into the mechanisms that drive biodiversity patterns –as will be discussed in the following–, which may be obscured if only one is considered (Forest et al., 2007; Swenson, 2013).

4.2 Taxonomic diversity

Our results indicated a monotonic decrease in taxonomic diversity along elevation across all study transects that was consistent for different Hill orders (Fig. 3a), thereby suggesting that this decrease occurred regardless of whether taxa were weighted based on their abundance (q = 1 and q = 2) or not (q = 0).

General elevational decreases in taxonomic diversity of woody plants similar to those found in the present study have been widely reported (Kessler, 2001; Sandoya et al., 2017; Nottingham et al., 2018), being the elevational decrease in temperature the ultimate cause of this trend (Grytnes and McCain, 2007). However, studies have also detected a non-monotonic change in diversity along elevation, specifically a hump-shaped pattern where richness peaks at middle elevations (Gentry, 1995; Rahbek, 1995; Girardin et al., 2014; Salazar et al., 2015). This hump-shaped pattern in taxonomic diversity with elevation can be explained by the relatively stable condensation zone (cloud belt) that characterizes the middle elevations of tropical elevational gradients, which provides favourable conditions for most organisms despite the lower temperatures relative to premontane elevations (Kessler and Kluge, 2008; Huaraca Huasco et al., 2014). If this hump was present along our transects, our sampling design would have been unable to detect it because i) we did not sample below 1000 m a.s.l. (except for a very few plots in the Río Abiseo transect that did show an indication of a hump), and ii) there was a gap in our sampling around 1500 m a.s.l., (where the highest diversity could be expected). These sampling limitations may have concealed a potential hump in species richness (Rahbek, 2005; Nogués-Bravo et al., 2008).

4.3 Functional diversity

Our results indicated a monotonic decrease in functional diversity along elevation across all transects that was consistent for different Hill orders (Fig. 3b). Theory predicts that the range of successful functional strategies narrows along an environmental stress gradient (Weiher and Keddy, 1995; Cornwell and Ackerly, 2009). In the case of elevation, the upslope increasingly harsher and more restrictive conditions –mainly in terms of temperature and resource availability (Stevens, 1992; Sundqvist et al., 2013)– is expected to lead to the selective survival of taxa with conservative strategies –promoting storage and defence– that enable them to cope with the environment, and to the cull of these with acquisitive ones –enhancing photosynthesis and growth– (Poorter et al., 2009; Machac et al., 2011). This pattern has been detected previously in TMFs both for individual functional traits –including leaf size decreases (Schneider et al., 2003; Salinas et al., 2011; Llerena-Zambrano et al., 2021) and wood density increases (Swenson and Enquist, 2007; Slik et al., 2010)–, and for overall community functional diversity (Duivenvoorden and Cuello, 2012; Schellenberger Costa et al., 2017). This evidence has been used in support of environmental filtering as an essential driver of plant community assembly (Swenson et al., 2012; Bañares-de-Dios et al., 2020).

The decrease in functional diversity was less pronounced for Hill numbers of higher order. Our results showed that Hill numbers that gave more weight to the most common entities (q = 2) accounted for two times less diversity than when the entities were strictly weighted by their abundance (q = 1) and four times less than when all entities received the same weight (q = 0) (Fig. 3b). These findings reflect how dominant/abundant species tend to be located at the core of the functional hyperspace whereas rare species tend to be peripheral. Dominant/abundant species exhibit strong traits–environment relationships enabling them a good performance in a given environment, thereby leading to greater demographic success. Conversely, rare species display sub-optimal traits –these species are more reliant on occasional migration events or stochastic environmental fluctuations–, which may explain their poorer demographic success (Umaña et al., 2015, 2017). In highly diverse ecosystems, such as tropical rainforests or coral reefs, rare species exhibit the less common traits (functional rarity), and thus disproportionately contribute to functional diversity, being irreplaceable and their loss not susceptible to being cushioned by the high species richness of the ecosystem as usually assumed. The reason for this is that the most common species have been shown to only support a relatively limited functionality (functional redundancy) (Mouillot et al., 2013, 2014; Leitão et al., 2016). Indeed, functional rarity sustains key ecosystem processes and services and provides an invaluable buffer against climate change and human disturbance (Cadotte et al., 2011; Violle et al., 2017). Therefore, the justification for conserving rare taxa should extend far beyond their taxonomic uniqueness, past, charisma or cultural significance (Tucker et al., 2012).

4.4 Phylogenetic diversity

Our results indicate that phylogenetic diversity decreased along elevation for q = 0 and q = 1 but increased slightly for q = 2 (Fig. 3c). In general, the overall decrease agrees with the predictions of the Tropical Niche Conservatism (TNC) hypothesis, which suggests that biological communities in the highlands are phylogenetically a subset of the communities found in the lowlands, and thus a decrease in phylogenetic diversity is expected with elevation (Wiens and Graham, 2005; Donoghue, 2008). This pattern has been observed for other biota in Andean TMFs (e.g., birds (Graham et al., 2009) and moths (Brehm et al., 2013)), but previous studies of Andean woody plants seem to give more support to the Out of the Tropical Lowlands (OTL) hypothesis (González-Caro et al., 2014, 2020; Tiede et al., 2016; Ramírez et al., 2019). Under OTL, the presence of temperate extratropical lineages within high elevation communities would be responsible for an increased phylogenetic alpha-diversity and turnover (Qian and Ricklefs, 2016), which is the pattern we found when considering phylogenetic diversity for the most abundant species (q = 2).

Therefore, our results are not in full agreement with either hypothesis and endorse neither the support nor the rejection of the 10.13039/100014596 TNC or OTL (see also Tolmos et al., 2022). In this sense an alternative intermediate hypothesis termed Environmental Crossroads (Neves et al., 2020) has been proposed, which posits that floras from different origins –tropical and temperate– co-exist at middle elevations, at the edges –lower and upper, respectively– of their thermal tolerances (Griffiths et al., 2021). In any case, further research is needed, including determining the phylogenetic turnover and nestedness, as well as biogeographical origin of clades (Linan et al., 2021) and changes in their age along the elevational gradient (Qian, 2014; Griffiths et al., 2020). A sampling schema that replicates transects along a broader latitudinal range, such as that used in the present study, would also help to assess the relative importance of the TNC vs. OTL. For example, several extratropical lineages such as Quercus prevail at high elevations in Colombia, which may have partly skewed the results reported by González-Caro et al. (2014, 2020) or Ramírez et al. (2019) based on studies that were latitudinally restricted to the northern Andes and which found evidence in agreement of OTL. Quercus is a Holarctic genus that migrated to South America in the last ca. 0.5 My (Hooghiemstra and Van Der Hammen, 2004) and there may have not been sufficient time for it to reach the central Andes, where niche conservatism and the upslope migration of clades originated in tropical lowlands could have been the predominant mechanisms responsible for population by Andean flora.

In general, extratropical lineages that have migrated from temperate to tropical regions along the Andes, such as Quercus or Alnus in Colombia or Araucaria and Nothofagus in the southern Andes, are better adapted to cold conditions than strictly tropical lineages, and thus they are more likely to thrive in the uppermost TMF communities (Segovia et al., 2020). These taxa were absent from the transects that we studied but we found other floristic elements of Laurasian (e.g., Prunus and Cornus) and Gondwanan (e.g., Podocarpus and Retrophyllum) temperate origin, some of which tended to be dominant locally in the past and nowadays are far less abundant, with the only exception of certain fragments of relict forests (see Vicuña-Miñano, 2005; Yaguana et al., 2012; Huamantupa-Chuquimaco et al., 2017). In this regard, the fact that our transects were located on isolated, protected areas difficult to access may have favoured the preservation of the native flora, allowing us to detect a slight increase in phylogenetic diversity at the highest elevations when giving more weigh to the most abundant taxa (q = 2) (Fig. 3c). Indeed, this slight increase resembles the phylogenetic pattern expected under the OTL. Therefore, as one moves far away from the central towards the northern or southern Andes, it would not be surprising to find an increase in the positive relationship between phylogenetic diversity and elevation, which may be detected by simply considering the relative abundances of taxa (q = 1) or equally considering all taxa (q = 0) given the greater prevalence of extratropical lineages.

Finally, it is important to note that the metrics used in this study to calculate phylogenetic diversity are strongly correlated with the number of lineages in recent evolutionary time but have weaker correlations with the number of lineages deeper in the evolutionary history of an assemblage (Dexter et al., 2019). Thus, our results might differ with the use of different metrics or the inclusion of different phylogenetic depths (Bose et al., 2019).

5 Conclusions

We found an overall decrease in taxonomic, functional, and phylogenetic diversity with elevation in Andean tropical montane forests (TMFs) (the only exception was the phylogenetic diversity when we over-weighted dominant species). Our experimental design is particularly robust because we analysed four elevational gradients over a broad latitudinal range, thus avoiding the potential bias of other studies that focused on the much better known but geographically constrained northern Andean TMFs (Pérez-Escobar et al., 2022). Declines in taxonomic and functional diversity with elevation are consistent with the expected outcomes of an environmental filtering process in which temperature is the main driver of environmental variation. Differences in the steepness of patterns of functional diversity between different Hill orders suggest that rare species play key and unique roles in the ecosystem, thereby supporting the conservation of rare taxa for reasons that extend far beyond their taxonomic or historical singularity. In addition, these findings suggest that niche conservatism and the upslope migration of tropical clades may be the predominant mechanisms responsible for shaping the flora of the central Andean highlands, given the limited arrival of temperate taxa in comparison to northern Andes. Further research across the entire Andean range, by integrating analyses of phylogenetic turnover and nestedness, as well as the biogeographical origin of clades and changes in their age along the elevational gradient will be important for elucidating the biogeographic history of the Andean slopes. The fact that the three types of diversity were positively correlated demonstrates that preserving TMF areas that harbour high taxonomic diversity will also preserve high functional and phylogenetic diversity in one of the most threatened biodiversity hotspots on Earth, which is still poorly understood compared to tropical lowlands (ForestPlots.net 2021) and threatened by anthropogenic intervention and global warming (Mathez-Stiefel et al., 2017; Feeley et al., 2020; Salinas et al., 2021).

CRediT authorship contribution statement

Guillermo Bañares-de-Dios: Writing - review & editing, Writing - original draft, Visualization, Validation, Supervision, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. Manuel J. Macía: Writing - review & editing, Writing - original draft, Supervision, Project administration, Methodology, Investigation, Funding acquisition, Data curation, Conceptualization. Gabriel Arellano: Writing - review & editing, Writing - original draft, Validation, Investigation, Funding acquisition, Data curation, Conceptualization. Íñigo Granzow-de la Cerda: Writing - review & editing, Writing - original draft, Validation, Investigation, Data curation, Conceptualization. Julia Vega-Álvarez: Writing - review & editing, Visualization, Formal analysis. Itziar Arnelas: Writing - review & editing, Validation, Investigation, Data curation. Carlos I. Espinosa: Writing - review & editing, Project administration. Norma Salinas: Writing - review & editing, Project administration. Luis Cayuela: Writing - review & editing, Writing - original draft, Visualization, Validation, Supervision, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.

Declaration of competing interest

Authors declare no conflict of interests.

Appendix A Supplementary data

The following is the supplementary data to this article:Multimedia component 1

Multimedia component 1

Multimedia component 2

Multimedia component 2

Acknowledgements

Guillermo Bañares was funded through grants from the Spanish 10.13039/100009950 Ministry of Education (FPU14/05303), Escuela Internacional de Doctorado - 10.13039/501100007511 Universidad Rey Juan Carlos (Doctor Internacional 2017) and the Education and Research Department of Madrid Autonomous Region Government (REMEDINAL TE; S2018/EMT-4338). The study was supported through three grants from the Spanish Ministries of Economy and Competitiveness and Science and Technology (CGL2013-45634-P, CGL2016-75414-P, and PID2019-105064 GB-I00), and a grant from Centro de Estudios de América Latina (CEAL) at Universidad Autónoma de Madrid and 10.13039/100010784 Banco Santander . We are indebted to those who helped with fieldwork: Stalin Japón, Wilson Remacha, Percy Malqui, “Rosho” Tamayo, Reinerio Ishuiza, Manuel Marca, Carlos Salas, José Sánchez, Anselmo Vergaray, Gonzalo Bañares, Ángel Delso, Julia González de Aledo, Mara Paneghel, Melecio Sullca, Raúl Huasurco, Honorato Pinto, Santiago Mamani, Aníbal Mamani, Armando Torrez, Justo Salas, Juan Mamani, Beto Apaza Coronel, Marcelo Reguerín, Óscar Quiroga, Juan Carlos Mamani, Rogelio Mamani, Vladimir Chura, Martín Chacón, and many others. We are especially thankful to Alex Nina, Jorge Armijos, Pablo Soliz, Leslie Cayola, Alfredo Fuentes, Peter Jørgensen and Tolentino Cueva, and the people of Los Alisos (Pataz, Peru) for their invaluable help. In addition, we are very grateful to Ministerio del Ambiente (MAE) in Ecuador, Servicio Nacional de Áreas Naturales Protegidas (SERNANP) in Peru, particularly Hugo Macedo, Vladimir Ramírez, Octavio Pecho, and Jhonny Ramos, and Ministerio de Medio Ambiente y Agua (MMAyA) in Bolivia. We extend our thanks to all the national park rangers who helped us, especially Rafael Galán, Guillermo Aguilar, Tito Heras, Percy Franco, Berardo Rojas, and Grover Benites. Permits to work in protected areas were granted by national authorities: Ecuador (MAE-DNB-CM-2015-0016; N° 001-2019-IC-FLO-FAU-DPAZCH-UPN-VS/MA), Peru (001-2016-SERNANP-PNRA-10.13039/100007415 JEF ), and Bolivia (MDRAyMA-VBRFMA- DGBAP-UAPVS N°2869/08).

Peer review under responsibility of Editorial Office of Plant Diversity.

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.pld.2024.03.005.
==== Refs
References

Antonelli A. Kissling W.D. Flantua S.G.A. Geological and climatic influences on mountain biodiversity Nat. Geosci. 11 2018 718 725
Antonelli A. Nylander J.A.A. Persson C. Tracing the impact of the Andean uplift on Neotropical plant evolution Proc. Natl. Acad. Sci. U.S.A. 106 2009 9749 9754 19470489
Arellano G. Cala V. Fuentes A. A standard protocol for woody plant inventories and soil characterisation using temporary 0.1-Ha plots in tropical forests J. Trop. For. Sci. 28 2016 508 516
Báez S. Cayuela L. Macía M.J. FunAndes – a functional trait database of Andean plants Sci. Data 9 2022 511 35987763
Barton K. MuMIn: Multi-model inference. R package version 1.43.17 2013
Bañares-de-Dios G. Macía M.J. Granzow-de la Cerda Í. Linking patterns and processes of tree community assembly across spatial scales in tropical montane forests Ecology 101 2020 e03058
de Bello F. Lavorel S. Lavergne S. Hierarchical effects of environmental filters on the functional structure of plant communities: a case study in the French Alps Ecography 36 2013 393 402
Bose R. Ramesh B.R. Pélissier R. Phylogenetic diversity in the Western Ghats biodiversity hotspot reflects environmental filtering and past niche diversification of trees J. Biogeogr. 46 2019 145 157
Brehm G. Strutzenberger P. Fiedler K. Phylogenetic diversity of geometrid moths decreases with elevation in the tropical Andes Ecography 36 2013 1247 1253
Brum F.T. Graham C.H. Costa G.C. Global priorities for conservation across multiple dimensions of mammalian diversity Proc. Natl. Acad. Sci. U.S.A. 114 2017 7461 7466
Cadotte M.W. Carscadden K. Mirotchnick N. Beyond species: functional diversity and the maintenance of ecological processes and services J. Appl. Ecol. 48 2011 1079 1087
Caldas F. Memoria sobre la nivelación de las plantas que se cultivan en la vecindad del Ecuador 1803
Cavender-Bares J. Kozak K.H. Fine P.V.A. The merging of community ecology and phylogenetic biology Ecol. Lett. 12 2009 693 715 19473217
Cayuela L. Granzow-de la Cerda I. Albuquerque F.S. TAXONSTAND: an R package for species names standardisation in vegetation databases Methods Ecol. Evol. 3 2012 1078 1083
Chao A. Chiu C.H. Jost L. Phylogenetic diversity measures based on Hill numbers Phil. Trans. Biol. Sci. 365 2010 3599 3609
Chao A. Chiu C.H. Jost L. Unifying species diversity, phylogenetic diversity, functional diversity, and related similarity and differentiation measures through hill numbers Annu. Rev. Ecol. Evol. Syst. 45 2014 297 324
Chao A. Gotelli N.J. Hsieh T. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies Ecol. Monogr. 84 2014 45 67
Chave J. Coomes D. Jansen S. Towards a worldwide wood economics spectrum Ecol. Lett. 12 2009 351 366 19243406
Chiu C.H. Chao A. Distance-based functional diversity measures and their decomposition: a framework based on Hill numbers PLoS One 9 2014 e100014
Chun J. Lee C. Disentangling the local-scale drivers of taxonomic, phylogenetic and functional diversity in woody plant assemblages along elevational gradients in South Korea PLoS One 12 2017 e0185763
Cisneros L.M. Burgio K.R. Dreiss L.M. Multiple dimensions of bat biodiversity along an extensive tropical elevational gradient J. Anim. Ecol. 83 2014 1124 1136 24428636
Cornelissen J.H.C. Lavorel S. Garnier E. A handbook of protocols for standardised and easy measurement of plant functional traits worldwide Aust. J. Bot. 51 2003 335 380
Cornwell W.K. Ackerly D.D. Community assembly and shifts in the distribution of functional trait values across an environmental gradient in coastal California Ecol. Monogr. 79 2009 109 126
Dani R.S. Divakar P.K. Baniya C.B. Diversity and composition of plant species along elevational gradient: research trends Biodivers. Conserv. 32 2023 2961 2980
Dehling D.M. Fritz S.A. Töpfer T. Functional and phylogenetic diversity and assemblage structure of frugivorous birds along an elevational gradient in the tropical Andes Ecography 37 2014 1047 1055
Devictor V. Mouillot D. Meynard C. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world Ecol. Lett. 13 2010 1030 1040 20545736
Dexter K.G. Segovia R.A. Griffiths A.R. Exploring the concept of lineage diversity across North American forests Forests 10 2019 520
Donoghue M.J. A phylogenetic perspective on the distribution of plant diversity Proc. Natl. Acad. Sci. U.S.A. 105 2008 11549 11555 18695216
Dreiss L.M. Burgio K.R. Cisneros L.M. Taxonomic, functional, and phylogenetic dimensions of rodent biodiversity along an extensive tropical elevational gradient Ecography 38 2015 1 13
Duivenvoorden J.F. Cuello N.L. Functional trait state diversity of Andean forests in Venezuela changes with altitude J. Veg. Sci. 23 2012 1105 1113
Feeley K. Martínez-Villa J. Pérez T. The thermal tolerances, distributions, and performances of tropical montane tree species Front. Forests Glob. Change 3 2020 1 11
Forest F. Grenyer R. Rouget M. Preserving the evolutionary potential of floras in biodiversity hotspots Nature 445 2007 757 760 17301791
ForestPlots.netBlundo C. Carilla J. Taking the pulse of Earth's tropical forests using networks of highly distributed plots Biol. Conserv. 260 2021 108849
Freckleton R.P. Jetz W. Space versus phylogeny: disentangling phylogenetic and spatial signals in comparative data Proc. R. Soc. B-Biol. Sci. 276 2009 21 30
Gentry A.H. Patterns of diversity and floristic composition in Neotropical montane forests Churchill S. Balsev H. Forero E. Biodiversity and Conservation of Neotropical Montane Forests 1995 The New York Botanical Garden New York 103 126
Gerhold P. Carlucci M.B. Procheş Ş. The deep past controls the phylogenetic structure of present, local communities Annu. Rev. Ecol. Evol. Syst. 49 2018 477 497
Girardin C.A.J. Farfan-Rios W. Garcia K. Spatial patterns of above-ground structure, biomass and composition in a network of six Andean elevation transects Plant Ecol. Divers. 7 2014 161 171
González-Caro S. Duque Á. Feeley K.J. The legacy of biogeographic history on the composition and structure of Andean forests Ecology 101 2020 1 11
González-Caro S. Umaña M.N. Álvarez E. Phylogenetic alpha and beta diversity in tropical tree assemblages along regional-scale environmental gradients in northwest South America J. Plant Ecol. 7 2014 145 153
Graham C.H. Carnaval A.C. Cadena C.D. The origin and maintenance of montane diversity: integrating evolutionary and ecological processes Ecography 37 2014 711 719
Graham C.H. Parra J. Rahbek C. Phylogenetic structure in tropical hummingbird communities Proc. Natl. Acad. Sci. U.S.A. 106 2009 19673 19678 19805042
Griffiths A.R. Silman M.R. Farfan-Rios W. Evolutionary diversity peaks at mid-elevations along an Amazon-to-Andes elevation gradient Front. Ecol. Evol. 9 2021 1 10
Griffiths A.R. Silman M.R. Farfan Rios W. Evolutionary heritage shapes tree distributions along an Amazon-to-Andes elevation gradient Biotropica 53 2020 38 50
Grytnes J.A. McCain C.M. Elevational trends in biodiversity Levin S.A. Encyclopedia of Biodiversity 2007 Elsevier 1 8
Guo Q. Kelt D.A. Sun Z. Global variation in elevational diversity patterns Sci. Rep. 3 2013 3007 24157658
Hill M. Diversity and evenness: a unifying notation and its consequences Ecology 54 1973 427 432
Hooghiemstra H. Van Der Hammen T. Quaternary Ice-Age dynamics in the Colombian Andes: developing an understanding of our legacy Phil. Trans. Roy. Soc. Lond. B 359 2004 173 181 15101574
Huamantupa-Chuquimaco I. Luza-Victorio M. Alfaro-Curitumay M. Diversidad y Biomasa Arbórea en los Bosques Andinos del Santuario Nacional del Ampay, Apurímac – perú Q’Euña - Sociedad Botánica del Cusco 8 2017 7 26
Huaraca Huasco W. Girardin C.A.J. Doughty C.E. Seasonal production, allocation and cycling of carbon in two mid-elevation tropical montane forest plots in the Peruvian Andes Plant Ecol. Divers. 7 2014 125 142
Hughes C. Eastwood R. Island radiation on a continental scale: exceptional rates of plant diversification after uplift of the Andes Proc. Natl. Acad. Sci. U.S.A. 103 2006 10334 10339 16801546
Humboldt A.V. Bonpland A. Essai sur la geographie des plantes 1805 Chez Levrault Schoell Paris
Jablonski D. Kaustuv R. Valentine J.W. Out of the Tropics: evolutionary diversity gradient Science 314 2006 102 106 17023653
Janzen D.H. Why mountain passes are higher in the tropics Am. Nat. 101 1967 233 249
Jin Y. Qian H. V.PhyloMaker2: an updated and enlarged R package that can generate very large phylogenies for vascular plants Plant Divers. 44 2022 335 339 35967255
Jost L. Entropy and diversity Oikos 113 2006 363 375
Karger D.N. Conrad O. Böhner J. Climatologies at high resolution for the earth’ s land surface areas Sci. Data 4 2017 170122
Kerkhoff A.J. Moriarty P. Weiser M. The latitudinal species richness gradient in the New World woody angiosperms is consistent with the tropical conservatism hypothesis Proc. Natl. Acad. Sci. U.S.A. 111 2014 8125 8130 24847062
Kessler M. Patterns of diversity and range size of selected plant groups along an elevational transect in the Bolivian Andes Biodivers. Conserv. 10 2001 1897 1921
Kessler M. Kluge J. Tropical mountain forest: patterns and processes in a Biodiversity hotspot Gradstein S.R. Homeier J. Kluge J. Biodiversity and Ecology Series vol. 2 2008 Univ. Göttingen, Göttingen 7 24
Körner C. Mountain biodiversity, its causes and function Ambio 33 2004 11 17
Körner C. The use of ’altitude’ in ecological research Trends Ecol. Evol. 22 2007 569 574 17988759
Laraque A. Ronchail J. Pombosa R. Heterogeneous distribution of rainfall and discharge regimes in the Ecuadorian Amazon basin J. Hydrometeorol. 8 2007 1364 1381
Latham R. Ricklefs R.E. Global patterns of tree species in moist forests: energy-diversity theory does not account for variation in species richness Oikos 67 1993 325 333
Leitão R.P. Zuanon J. Villéger S. Rare species contribute disproportionately to the functional structure of species assemblages Proc. R. Soc. B-Biol. Sci. 283 2016 20160084
Linan A.G. Myers J.A. Edwards C.E. The evolutionary assembly of forest communities along environmental gradients: recent diversification or sorting of pre-adapted clades? New Phytol. 232 2021 2506 2519 34379801
Llerena-Zambrano M. Ordoñez J.C. Llambí L.D. Minimum temperature drives community leaf trait variation in secondary montane forests along a 3000-m elevation gradient in the tropical Andes Plant Ecol. Divers. 14 2021 47 63
Lomolino M.V. Elevation gradients of species-density: historical and prospective views Global Ecol. Biogeogr. 10 2001 3 13
Losos J.B. Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species Ecol. Lett. 11 2008 995 1003 18673385
Machac A. Janda M. Dunn R.R. Elevational gradients in phylogenetic structure of ant communities reveal the interplay of biotic and abiotic constraints on diversity Ecography 34 2011 364 371
Malhi Y. Silman M. Salinas N. Introduction: elevation gradients in the tropics: laboratories for ecosystem ecology and global change research Global Change Biol. 16 2010 3171 3175
Mathez-Stiefel S.L. Peralvo M. Báez S. Research priorities for the conservation and sustainable governance of andean forest landscapes. Mountain Res Dev. 37 2017 323 339
McCain M.C. Elevational gradients in diversity of small mammals Ecology 86 2005 366 372
McGill B.J. Enquist B.J. Weiher E. Rebuilding community ecology from functional traits Trends Ecol. Evol. 21 2006 178 185 16701083
Mittermeier R.A. Myers N. Mittermeier C.G. Hotspots Revisited: Earth's Biologically Richest and Most Endangered Terrestrial Ecoregions 2005 Conservation International Washington DC
Moles A.T. Ackerly D.D. Webb C.O. A brief history of seed size Science 307 2005 576 580 15681384
Mouillot D. Bellwood D.R. Baraloto C. Rare species support vulnerable functions in high-diversity ecosystems PLoS Biology 11 2013 e1001569
Mouillot D. Villéger S. Parravicini V. Functional over-redundancy and high functional vulnerability in global fish faunas on tropical reefs Proc. Natl. Acad. Sci. U.S.A. 111 2014 13757 13762 25225388
Myers N. Mittermeier R.A. Mittermeier C.G. Biodiversity hotspots for conservation priorities Nature 403 2000 853 858 10706275
Neves D.M. Dexter K.G. Baker T.R. Evolutionary diversity in tropical tree communities peaks at intermediate precipitation Sci. Rep. 10 2020 1 7 31913322
Nogués-Bravo D. Araújo M.B. Romdal T. Scale effects and human impact on the elevational species richness gradients Nature 453 2008 216 219 18464741
Nottingham A.T. Fierer N. Turner B.L. Microbes follow Humboldt: temperature drives plant and soil microbial biodiversity patterns from the Amazon to the Andes Ecology 99 2018 2455 2466 30076592
Oksanen J. Blanchet F.G. Friendly M. https://CRAN.R-project.org/package=vegan 2019
Pavoine S. Bonsall M.B. Measuring biodiversity to explain community assembly: a unified approach Biol. Rev. Camb. Phil. Soc. 86 2011 792 812
Pérez-Escobar O.A. Zizka A. Bermúdez M.A. The Andes through time: evolution and distribution of Andean floras Trends Plant Sci. 27 2022 364 378 35000859
Petchey K.J. Gaston O.L. Functional diversity: back to basics and looking forward Ecol. Lett. 9 2006 741 758 16706917
Poorter H. Niinemets Ü. Poorter L. Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis New Phytol. 182 2009 565 588 19434804
Qian H. Contrasting relationships between clade age and temperature along latitudinal versus elevational gradients for woody angiosperms in forests of South America J. Veg. Sci. 25 2014 1208 1215
Qian H. Ricklefs R.E. Out of the tropical lowlands: latitude versus elevation Trends Ecol. Evol. 31 2016 738 741 27523604
R Core Team R: A Language and Environment for Statistical Computing 2021 R Foundation for Statistical Computing https://www.r-project.org/
Rahbek C. The elevational gradient of species richness: a uniform pattern? Ecography 18 1995 200 205
Rahbek C. The role of spatial scale and the perception of large-scale species-richness patterns Ecol. Lett. 8 2005 224 239
Rahbek C. Borregaard M.K. Antonelli A. Building mountain biodiversity: geological and evolutionary processes Science 365 2019 1114 1119 31515384
Rahbek C. Borregaard M.K. Colwell R.K. Humboldt's enigma: what causes global patterns of mountain biodiversity? Science 365 2019 1108 1113 31515383
Ramírez S. González-Caro S. Phillips J. The influence of historical dispersal on the phylogenetic structure of tree communities in the tropical Andes Biotropica 51 2019 500 508
Ricklefs R.E. Schluter D. Ricklefs R.E. Schluter D. Species Diversity in Ecological Communities: Historical and Geographical Perspectives 1993 University of Chicago Press Chicago 523 525
Safi K. Cianciaruso M.V. Loyola R.D. Understanding global patterns of mammalian functional and phylogenetic diversity Phil. Trans. Biol. Sci. 366 2011 2536 2544
Salazar L. Homeier J. Kessler M. Diversity patterns of ferns along elevational gradients in Andean tropical forests Plant Ecol. Divers. 8 2015 13 24
Salinas N. Cosio E.G. Silman M. Editorial: tropical montane forests in a changing environment Front. Plant Sci. 12 2021 1652
Salinas N. Malhi Y. Meir P. The sensitivity of tropical leaf litter decomposition to temperature: results from a large-scale leaf translocation experiment along an elevation gradient in Peruvian forests New Phytol. 189 2011 967 977 21077887
Sanders N.J. Rahbek C. The patterns and causes of elevational diversity gradients Ecography 35 2012 1 3
Sandoya V. Pauchard A. Cavieres L.A. Natives and non-natives plants show different responses to elevation and disturbance on the tropical high Andes of Ecuador Ecol. Evol. 7 2017 7909 7919 29043044
Sanmartín I. Historical biogeography: evolution in time and space Evo. Edu. Outreach 5 2012 555 568
Schellenberger Costa D. Gerschlauer F. Pabst H. Community-weighted means and functional dispersion of plant functional traits along environmental gradients on Mount Kilimanjaro J. Veg. Sci. 28 2017 684 695
Schneider J. Zipp D. Gaviria J. Successional and mature stands in an upper Andean rain forest transect of Venezuela: do leaf characteristics of woody species differ J. Trop. Ecol. 19 2003 251 259
Segovia R.A. Pennington R.T. Baker T.R. Freezing and water availability structure the evolutionary diversity of trees across the Americas Sci. Adv. 6 2020 eaaz5373
Slik J.W.F. Aiba S.I. Brearley F.Q. Environmental correlates of tree biomass, basal area, wood specific gravity and stem density gradients in Borneo's tropical forests Global Ecol. Biogeogr. 19 2010 50 60
Stevens G.C. The elevational gradient in altitudinal range: an extension of Rapoport's latitudinal rule to altitude Am. Nat. 140 1992 893 911 19426029
Storch D. Marquet P.A. Gaston K.J. Untangling an entangled bank Science 307 2005 684 686 15692039
Sundqvist M.K. Sanders N.J. Wardle D.A. Community and ecosystem responses to elevational gradients: processes, mechanisms, and insights for global change Annu. Rev. Ecol. Evol. Syst. 44 2013 261 280
Swenson N.G. The role of evolutionary processes in producing biodiversity patterns, and the interrelationships between taxonomic, functional, and phylogenetic biodiversity Am. J. Bot. 98 2011 472 480 21613140
Swenson N.G. The assembly of tropical tree communities - the advances and shortcomings of phylogenetic and functional trait analyses Ecography 36 2013 264 276
Swenson N.G. Enquist B.J. Ecological and evolutionary determinants of a key plant functional trait, wood density and its community wide variation across latitude and elevation Am. J. Bot. 94 2007 451 459 21636415
Swenson N.G. Enquist B.J. Pither J. The biogeography and filtering of woody plant functional diversity in North and South America Global Ecol. Biogeogr. 21 2012 798 808
Tanaka T. Sato T. Taxonomic, phylogenetic and functional diversities of ferns and lycophytes along an elevational gradient depend on taxonomic scales Plant Ecol. 216 2015 1597 1609
Thakur D. Chawla A. Functional diversity along elevational gradients in the high altitude vegetation of the western Himalaya Biodivers. Conserv. 28 2019 1977 1996
Thiers B Index Herbariorum: A global directory of public herbaria and associated Staff. New York Botanical Garden’s Virtual Herbarium http://sweetgum.nybg.org/ih/
Tiede Y. Homeier J. Cumbicus N. Phylogenetic niche conservatism does not explain elevational patterns of species richness, phylodiversity and family age of tree assemblages in andean rainforest Erdkunde 70 2016 83 106
Tito R. Vasconcelos H.L. Feeley K.J. Mountain ecosystems as natural laboratories for climate change experiments Frontiers Forests Glob. Change 3 2020 1 8
Tolmos M.L. Kreft H. Ramirez J. Water and energy availability mediate biodiversity patterns along an elevational gradient in the tropical Andes J. Biogeogr. 49 2022 712 726
Tucker C.M. Cadotte M.W. Davies T.J. Incorporating geographical and evolutionary rarity into conservation prioritization Conserv. Biol. 26 2012 593 601 22624595
Umaña M.N. Zhang C. Cao M. Commonness, rarity, and intraspecific variation in traits and performance in tropical tree seedlings Ecol. Lett. 18 2015 1329 1337 26415689
Umaña M.N. Zhang C. Cao M. A core-transient framework for trait-based community ecology: an example from a tropical tree seedling community Ecol. Lett. 20 2017 619 628 28371151
Venables W. Ripley B. Modern Applied Statistics with S 2002 Springer New York
Vicuña-Miñano E.E. Las Podocarpáceas de los bosques montanos del noroccidente peruano Rev. Peru. Biol. 12 2005 283 288
Violle C. Navas M.-L. Vile D. Let the concept of trait be functional Oikos 116 2007 882 892
Violle C. Thuiller W. Mouquet N. Functional rarity: the ecology of outliers Trends Ecol. Evol. 32 2017 356 367 28389103
Wambulwa M.C. Milne R. Wu Z.T. Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: current knowledge and future perspectives Ecol. Evol. 11 2021 10794 10812 34429882
Weiher E. Keddy P.A. Assembly rules, null models, and trait dispersion: new questions from old patterns Oikos 74 1995 159 164
Wiens J. Donoghue M.J. Historical biogeography, ecology and species richness Trends Ecol. Evol. 19 2004 639 644 16701326
Wiens J. Graham C.H. Niche conservatism: integrating evolution, ecology and conservation biology Annu. Rev. Ecol. Evol. Syst. 36 2005 519 539
Wright I.J. Reich P.B. Westoby M. The worldwide leaf economics spectrum Nature 428 2004 821 827 15103368
Yaguana C. Lozano D. Neill D. Diversidad florística y estructura del bosque nublado del Río Numbala, Zamora-Chinchipe, Ecuador Revista Amazónica Ciencia y Tecnología 1 2012 226 247
