Assembly of forest communities across East Asia--insights from phylogenetic community structure and species pool scaling

Global patterns and climatic determinants of phylogenetic structure of regional fern floras

Knowledge of relationships between phylogenetic structure of a biological assemblage and ecological factors that drive the variation of phylogenetic structure among regions is crucial for understanding the causes of variation in taxonomic composition and richness among regions, but this knowledge is lacking for the global flora of ferns. Here, we fill this critical knowledge gap. We divided the globe into 392 geographic units on land, collated species lists of ferns for each geographic unit, and used different phylogenetic metrics (tip- vs basal-weighted) reflecting different evolutionary depths to quantify phylogenetic structure. We then related taxonomic and phylogenetic structure metrics to six climatic variables for ferns as a whole and for two groups of ferns (old clades vs polypods) reflecting different evolutionary histories across the globe and within each continental region. We found that when old clades and polypods were considered separately, temperature-related variables explained more variation in these metrics than did precipitation-related variables in both groups. When analyses were conducted for continental regions separately, this pattern holds in most cases. Climate extremes have a stronger relationship with phylogenetic structure of ferns than does climate seasonality. Climatic variables explained more variation in phylogenetic structure at deeper evolutionary depths.

Plant community phylogeny responses to protections and its main drivers in boreal forests, China: General pattern and implications

Patterns of the phylogenetic structure have been broadly applied to predict community assembly processes. However, the distribution pattern of evolutionary diversity and its drivers under nature conservation are still poorly understood in boreal forests. Here, we investigated 1738 sampling plots and subplots from distinct protection intensities (PIs) zones in five representative National Nature Reserves (NNRs). Multiple comparisons, redundancy analysis, and linear mixed model were performed to identify the changes in community phylogeny across different PIs and NNRs and the drivers for these variations. Our results showed considerable plant community phylogeny variations in different NNRs. As indicated by SesMPD (standardized mean pairwise distance) and SesMNTD (standardized the mean nearest taxon distance), trees, shrubs, and herbs presented overdispersed, clustered, and random distribution patterns, respectively, in different PIs. Protection resulted in the phylogenetic structure between the nearest species of trees showing a more overdispersed pattern (p < 0.05). Protection decreased the phylogenetically clustered degree between the nearest species of shrubs (p > 0.05), while the herbs still maintained a random pattern. Community traits explained the most to phylogeny variation of different communities (24 %-71 %, p < 0.01), followed by geoclimatic factors (2 %-24 %) and conservation processes (1 %-21 %). The higher mean annual precipitation and under branch height at the lower latitude area accompanied the higher SesMPD and SesMNTD. The higher PIs attended with higher tree SesMPD, and the longer protection time resulted in higher shrub PSR (phylogenetic species richness) and PSV (phylogenetic species variability). Including the location of NNRs, community traits, and years of protection, rather than only emphasizing PI itself, could optimize community phylogenetic structure and preserve the evolutionary potential of biodiversity.

A signal of competitive dominance in mid-latitude herbaceous plant communities

Understanding the main determinants of species coexistence across space and time is a central question in ecology. However, ecologists still know little about the scales and conditions at which biotic interactions matter and how these interact with the environment to structure species assemblages. Here we use recent theoretical developments to analyse plant distribution and trait data across Europe and find that plant height clustering is related to both evapotranspiration (ET) and gross primary productivity. This clustering is a signal of interspecies competition between plants, which is most evident in mid-latitude ecoregions, where conditions for growth (reflected in actual ET rates and gross primary productivities) are optimal. Away from this optimum, climate severity probably overrides the effect of competition, or other interactions become increasingly important. Our approach bridges the gap between species-rich competition theories and large-scale species distribution data analysis.

Relationships between climate and phylogenetic community structure of fossil pollen assemblages are not constant during the last deglaciation

Disentangling the influence of environmental drivers on community assembly is important to understand how multiple processes influence biodiversity patterns and can inform understanding of ecological responses to climate change. Phylogenetic Community Structure (PCS) is increasingly used in community assembly studies to incorporate evolutionary perspectives and as a proxy for trait (dis)similarity within communities. Studies often assume a stationary relationship between PCS and climate, though few studies have tested this assumption over long time periods with concurrent community data. We estimated two PCS metrics-Nearest Taxon Index (NTI) and Net Relatedness index (NRI)-of fossil pollen assemblages of Angiosperms in eastern North America over the last 21 ka BP at 1 ka intervals. We analyzed spatiotemporal relationships between PCS and seven climate variables, evaluated the potential impact of deglaciation on PCS, and tested for the stability of climate-PCS relationships through time. The broad scale geographic patterns of PCS remained largely stable across time, with overdispersion tending to be most prominent in the central and southern portion of the study area and clustering dominating at the longitudinal extremes. Most importantly, we found that significant relationships between climate variables and PCS (slope) were not constant as climate changed during the last deglaciation and new ice-free regions were colonized. We also found weak, but significant relationships between both PCS metrics (i.e., NTI and NRI) and climate and time-since-deglaciation that also varied through time. Overall, our results suggest that (1) PCS of fossil Angiosperm assemblages during the last 21ka BP have had largely constant spatial patterns, but (2) temporal variability in the relationships between PCS and climate brings into question their usefulness in predictive modeling of community assembly.

Using metacommunity ecology to understand environmental metabolomes

Environmental metabolomes are fundamentally coupled to microbially-linked biogeochemical processes within ecosystems. However, significant gaps exist in our understanding of their spatiotemporal organization, limiting our ability to uncover transferrable principles and predict ecosystem function. We propose that a theoretical paradigm, which integrates concepts from metacommunity ecology, is necessary to reveal underlying mechanisms governing metabolomes. We call this synthesis between ecology and metabolomics ‘meta-metabolome ecology’ and demonstrate its utility using a mass spectrometry dataset. We developed three relational metabolite dendrograms using molecular properties and putative biochemical transformations and performed ecological null modeling. Based upon null modeling results, we show that stochastic processes drove molecular properties while biochemical transformations were structured deterministically. We further suggest that potentially biochemically active metabolites were more deterministically assembled than less active metabolites. Understanding variation in the influences of stochasticity and determinism provides a way to focus attention on which meta-metabolomes and which parts of meta-metabolomes are most likely to be important to consider in mechanistic models. We propose that this paradigm will allow researchers to study the connections between ecological systems and their molecular processes in previously inaccessible detail.

Despite growing interest in environmental metabolomics, we lack conceptual frameworks for considering how metabolites vary across space and time in ecological systems. Here, the authors apply (species) community assembly concepts to metabolomics data, offering a way forward in understanding the assembly of metabolite assemblages.

Latitudinal and elevational patterns of phylogenetic structure in forest communities in China's mountains

The phylogenetic structure incorporates both ecological and evolutionary processes to explain assembly of a local community. The "phylogenetic niche conservatism" (PNC) hypothesis suggests that distributions of species along environmental gradients reflect both ancestral traits and ecological fitness of individual species The temperature is generally regarded to change in similar ways along both latitudinal and elevational gradients but with different historical contingence. Therefore, comparing the latitudinal and elevational patterns of phylogenetic structure of communities is of help to depict the effects of ecological and evolutionary processes in shaping the community assembly. In this study, we explored the latitudinal, elevational and climatic patterns of phylogenetic structure of 569 angiosperm tree communities from 38 mountains across China. We found a larger mean abundance-weighted net relatedness index (NRI) than the presence/absence-based NRI; and the NRI decreased when the species pool downscaled from the full pool to county-level pool. The mean family age and phylogenetic species evenness decreased with latitude, and increased with temperature of the coldest month and precipitation; whilst NRI increased with latitude, and decreased with mean temperature of the coldest month. In most mountains, NRI, mean family age and phylogenetic species evenness showed non-significant trends along the elevational gradient. Our results support the main predictions of PNC for the latitudinal gradient, i.e., species tend to be more phylogenetically related to each other and clades are younger in temperate environments, compared to those in tropical environments. We suggested that independent species pools and abundance should be incorporated in analysis to fully represent the phylogenetic structure of communities.

Phylogenetic structure and formation mechanism of shrub communities in arid and semiarid areas of the Mongolian Plateau

The mechanisms of species coexistence within a community have always been the focus in ecological research. Community phylogenetic structure reflects the relationship of historical processes, regional environments, and interactions between species, and studying it is imperative to understand the formation and maintenance mechanisms of community composition and biodiversity. We studied the phylogenetic structure of the shrub communities in arid and semiarid areas of the Mongolian Plateau. First, the phylogenetic signals of four plant traits (height, canopy, leaf length, and leaf width) of shrubs and subshrubs were measured to determine the phylogenetic conservation of these traits. Then, the net relatedness index (NRI) of shrub communities was calculated to characterize their phylogenetic structure. Finally, the relationship between the NRI and current climate and paleoclimate (since the Last Glacial Maximum, LGM) factors was analyzed to understand the formation and maintenance mechanisms of these plant communities. We found that desert shrub communities showed a trend toward phylogenetic overdispersion; that is, limiting similarity was predominant in arid and semiarid areas of the Mongolian Plateau despite the phylogenetic structure and formation mechanisms differing across habitats. The typical desert and sandy shrub communities showed a significant phylogenetic overdispersion, while the steppified desert shrub communities showed a weak phylogenetic clustering. It was found that mean winter temperature (i.e., in the driest quarter) was the major factor limiting steppified desert shrub phylogeny distribution. Both cold and drought (despite having opposite consequences) differentiated the typical desert to steppified desert shrub communities. The increase in temperature since the LGM is conducive to the invasion of shrub plants into steppe grassland, and this process may be intensified by global warming.

Phylogenetic study documents different speciation mechanisms within the Russula globispora lineage in boreal and arctic environments of the Northern Hemisphere

The Russula globispora lineage is a morphologically and phylogenetically well-defined group of ectomycorrhizal fungi occurring in various climatic areas. In this study we performed a multi-locus phylogenetic study based on collections from boreal, alpine and arctic habitats of Europe and Western North America, subalpine collections from the southeast Himalayas and collections from subtropical coniferous forests of Pakistan. European and North American collections are nearly identical and probably represent a single species named R. dryadicola distributed from the Alps to the Rocky Mountains. Collections from the southeast Himalayas belong to two distinct species: R. abbottabadensis sp. nov. from subtropical monodominant forests of Pinus roxburghii and R. tengii sp. nov. from subalpine mixed forests of Abies and Betula. The results suggest that speciation in this group is driven by a climate disjunction and adaptation rather than a host switch and geographical distance.

Estimating species pools for a single ecological assemblage

BACKGROUND

The species pool concept was formulated over the past several decades and has since played an important role in explaining multi-scale ecological patterns. Previous statistical methods were developed to identify species pools based on broad-scale species range maps or community similarity computed from data collected from many areas. No statistical method is available for estimating species pools for a single local community (sampling area size may be very small as ≤ 1 km²). In this study, based on limited local abundance information, we developed a simple method to estimate the area size and richness of a species pool for a local ecological community. The method involves two steps. In the first step, parameters from a truncated negative trinomial model characterizing the distributional aggregation of all species (i.e., non-random species distribution) in the local community were estimated. In the second step, we assume that the unseen species in the local community are most likely the rare species, only found in the remaining part of the species pool, and vice versa, if the remaining portion of the pool was surveyed and was contrasted with the sampled area. Therefore, we can estimate the area size of the pool, as long as an abundance threshold for defining rare species is given. Since the size of the pool is dependent on the rarity threshold, to unanimously determine the pool size, we developed an optimal method to delineate the rarity threshold based on the balance of the changing rates of species absence probabilities in the sampled and unsampled areas of the pool.

RESULTS

For a 50 ha (0.5 km²) forest plot in the Barro Colorado Island of central Panama, our model predicted that the local, if not regional, species pool for the 0.5 km² forest plot was nearly the entire island. Accordingly, tree species richness in this pool was estimated as around 360. When the sampling size was smaller, the upper bound of the 95% confidence interval could reach 418, which was very close to the flora record of tree richness for the island. A numerical test further demonstrated the power and reliability of the proposed method, as the true values of area size and species richness for the hypothetical species pool have been well covered by the 95% confidence intervals of the true values.

CONCLUSIONS

Our method fills the knowledge gap on estimating species pools for a single local ecological assemblage with little information. The method is statistically robust and independent of sampling size, as proved by both empirical and numerical tests.

Patterns of taxonomic, phylogenetic diversity during a long-term succession of forest on the Loess Plateau, China: insights into assembly process

Quantifying the drivers underlying the distribution of biodiversity during succession is a critical issue in ecology and conservation, and also can provide insights into the mechanisms of community assembly. Ninety plots were established in the Loess Plateau region of northern Shaanxi in China. The taxonomic and phylogenetic (alpha and beta) diversity were quantified within six succession stages. Null models were used to test whether phylogenetic distance observed differed from random expectations. Taxonomic beta diversity did not show a regular pattern, while phylogenetic beta diversity decreased throughout succession. The shrub stage occurred as a transition from phylogenetic overdispersion to clustering either for NRI (Net Relatedness Index) or betaNRI. The betaNTI (Nearest Taxon Index) values for early stages were on average phylogenetically random, but for the betaNRI analyses, these stages were phylogenetically overdispersed. Assembly of woody plants differed from that of herbaceous plants during late community succession. We suggest that deterministic and stochastic processes respectively play a role in different aspects of community phylogenetic structure for early succession stage, and that community composition of late succession stage is governed by a deterministic process. In conclusion, the long-lasting evolutionary imprints on the present-day composition of communities arrayed along the succession gradient.

Fast Computations for Measures of Phylogenetic Beta Diversity

For many applications in ecology, it is important to examine the phylogenetic relations between two communities of species. More formally, let [Formula: see text] be a phylogenetic tree and let A and B be two samples of its tips, representing the examined communities. We want to compute a value that expresses the phylogenetic diversity between A and B in [Formula: see text]. There exist several measures that can do this; these are the so-called phylogenetic beta diversity (β-diversity) measures. Two popular measures of this kind are the Community Distance (CD) and the Common Branch Length (CBL). In most applications, it is not sufficient to compute the value of a beta diversity measure for two communities A and B; we also want to know if this value is relatively large or small compared to all possible pairs of communities in [Formula: see text] that have the same size. To decide this, the ideal approach is to compute a standardised index that involves the mean and the standard deviation of this measure among all pairs of species samples that have the same number of elements as A and B. However, no method exists for computing exactly and efficiently this index for CD and CBL. We present analytical expressions for computing the expectation and the standard deviation of CD and CBL. Based on these expressions, we describe efficient algorithms for computing the standardised indices of the two measures. Using standard algorithmic analysis, we provide guarantees on the theoretical efficiency of our algorithms. We implemented our algorithms and measured their efficiency in practice. Our implementations compute the standardised indices of CD and CBL in less than twenty seconds for a hundred pairs of samples on trees with 7 ⋅ 10(4) tips. Our implementations are available through the R package PhyloMeasures.