The evolutionary assembly of forest communities along environmental gradients: recent diversification or sorting of pre-adapted clades?

Recent studies have demonstrated that ecological processes that shape community structure and dynamics change along environmental gradients. However, much less is known about how the emergence of the gradients themselves shape the evolution of species that underlie community assembly. In this study, we address how the creation of novel environments leads to community assembly via two nonmutually exclusive processes: immigration and ecological sorting of pre-adapted clades (ISPC), and recent adaptive diversification (RAD). We study these processes in the context of the elevational gradient created by the uplift of the Central Andes. We develop a novel approach and method based on the decomposition of species turnover into within- and among-clade components, where clades correspond to lineages that originated before mountain uplift. Effects of ISPC and RAD can be inferred from how components of turnover change with elevation. We test our approach using data from over 500 Andean forest plots. We found that species turnover between communities at different elevations is dominated by the replacement of clades that originated before the uplift of the Central Andes. Our results suggest that immigration and sorting of clades pre-adapted to montane habitats is the primary mechanism shaping tree communities across elevations.

Historical Assembly of Andean Tree Communities.

Patterns of species diversity have been associated with changes in climate across latitude and elevation. However, the ecological and evolutionary mechanisms underlying these relationships are still actively debated. Here, we present a complementary view of the well-known tropical niche conservatism (TNC) hypothesis, termed the multiple zones of origin (MZO) hypothesis, to explore mechanisms underlying latitudinal and elevational gradients of phylogenetic diversity in tree communities. The TNC hypothesis posits that most lineages originate in warmer, wetter, and less seasonal environments in the tropics and rarely colonize colder, drier, and more seasonal environments outside of the tropical lowlands, leading to higher phylogenetic diversity at lower latitudes and elevations. In contrast, the MZO hypothesis posits that lineages also originate in temperate environments and readily colonize similar environments in the tropical highlands, leading to lower phylogenetic diversity at lower latitudes and elevations. We tested these phylogenetic predictions using a combination of computer simulations and empirical analyses of tree communities in 245 forest plots located in six countries across the tropical and subtropical Andes. We estimated the phylogenetic diversity for each plot and regressed it against elevation and latitude. Our simulated and empirical results provide strong support for the MZO hypothesis. Phylogenetic diversity among co-occurring tree species increased with both latitude and elevation, suggesting an important influence on the historical dispersal of lineages with temperate origins into the tropical highlands. The mixing of different floras was likely favored by the formation of climatically suitable corridors for plant migration due to the Andean uplift. Accounting for the evolutionary history of plant communities helps to advance our knowledge of the drivers of tree community assembly along complex climatic gradients, and thus their likely responses to modern anthropogenic climate change.

FunAndes - A functional trait database of Andean plants.

We introduce the FunAndes database, a compilation of functional trait data for the Andean flora spanning six countries. FunAndes contains data on 24 traits across 2,694 taxa, for a total of 105,466 entries. The database features plant-morphological attributes including growth form, and leaf, stem, and wood traits measured at the species or individual level, together with geographic metadata (i.e., coordinates and elevation). FunAndes follows the field names, trait descriptions and units of measurement of the TRY database. It is currently available in open access in the FIGSHARE data repository, and will be part of TRY's next release. Open access trait data from Andean plants will contribute to ecological research in the region, the most species rich terrestrial biodiversity hotspot.

Mature Andean forests as globally important carbon sinks and future carbon refuges.

It is largely unknown how South America's Andean forests affect the global carbon cycle, and thus regulate climate change. Here, we measure aboveground carbon dynamics over the past two decades in 119 monitoring plots spanning a range of >3000 m elevation across the subtropical and tropical Andes. Our results show that Andean forests act as strong sinks for aboveground carbon (0.67 ± 0.08 Mg C ha-1 y-1) and have a high potential to serve as future carbon refuges. Aboveground carbon dynamics of Andean forests are driven by abiotic and biotic factors, such as climate and size-dependent mortality of trees. The increasing aboveground carbon stocks offset the estimated C emissions due to deforestation between 2003 and 2014, resulting in a net total uptake of 0.027 Pg C y-1. Reducing deforestation will increase Andean aboveground carbon stocks, facilitate upward species migrations, and allow for recovery of biomass losses due to climate change.

Elevation and latitude drives structure and tree species composition in Andean forests: Results from a large-scale plot network.

Our knowledge about the structure and function of Andean forests at regional scales remains limited. Current initiatives to study forests over continental or global scales still have important geographical gaps, particularly in regions such as the tropical and subtropical Andes. In this study, we assessed patterns of structure and tree species diversity along ~ 4000 km of latitude and ~ 4000 m of elevation range in Andean forests. We used the Andean Forest Network (Red de Bosques Andinos, https://redbosques.condesan.org/) database which, at present, includes 491 forest plots (totaling 156.3 ha, ranging from 0.01 to 6 ha) representing a total of 86,964 identified tree stems ≥ 10 cm diameter at breast height belonging to 2341 identified species, 584 genera and 133 botanical families. Tree stem density and basal area increases with elevation while species richness decreases. Stem density and species richness both decrease with latitude. Subtropical forests have distinct tree species composition compared to those in the tropical region. In addition, floristic similarity of subtropical plots is between 13 to 16% while similarity between tropical forest plots is between 3% to 9%. Overall, plots ~ 0.5-ha or larger may be preferred for describing patterns at regional scales in order to avoid plot size effects. We highlight the need to promote collaboration and capacity building among researchers in the Andean region (i.e., South-South cooperation) in order to generate and synthesize information at regional scale.

Elevational gradients in ß-diversity reflect variation in the strength of local community assembly mechanisms across spatial scales.

Despite long-standing interest in elevational-diversity gradients, little is known about the processes that cause changes in the compositional variation of communities (ß-diversity) across elevations. Recent studies have suggested that ß-diversity gradients are driven by variation in species pools, rather than by variation in the strength of local community assembly mechanisms such as dispersal limitation, environmental filtering, or local biotic interactions. However, tests of this hypothesis have been limited to very small spatial scales that limit inferences about how the relative importance of assembly mechanisms may change across spatial scales. Here, we test the hypothesis that scale-dependent community assembly mechanisms shape biogeographic ß-diversity gradients using one of the most well-characterized elevational gradients of tropical plant diversity. Using an extensive dataset on woody plant distributions along a 4,000-m elevational gradient in the Bolivian Andes, we compared observed patterns of ß-diversity to null-model expectations. ß-deviations (standardized differences from null values) were used to measure the relative effects of local community assembly mechanisms after removing sampling effects caused by variation in species pools. To test for scale-dependency, we compared elevational gradients at two contrasting spatial scales that differed in the size of local assemblages and regions by at least an order of magnitude. Elevational gradients in ß-diversity persisted after accounting for regional variation in species pools. Moreover, the elevational gradient in ß-deviations changed with spatial scale. At small scales, local assembly mechanisms were detectable, but variation in species pools accounted for most of the elevational gradient in ß-diversity. At large spatial scales, in contrast, local assembly mechanisms were a dominant force driving changes in ß-diversity. In contrast to the hypothesis that variation in species pools alone drives ß-diversity gradients, we show that local community assembly mechanisms contribute strongly to systematic changes in ß-diversity across elevations. We conclude that scale-dependent variation in community assembly mechanisms underlies these iconic gradients in global biodiversity.