| Site1 | |
|---|---|
| Sp. 1 | 4 |
| Sp. 2 | 300 |
| Sp. 3 | 56 |
| Sp. 4 | 23 |
Diversity and Trait-Based Approaches
🧑🏻💻 Masatoshi Katabuchi @ XTBG, CAS
November 11, 2024 XTBG AFEC
We Learn:
Why we use functional and phylogenetic diversity
How to calculate functional and phylogenetic diversity
Trait-based approaches
Community Ecology
Simple diversity indices
Phylogenetic diversity
Functional traits and diversity
R examples
For over a century, field ecologist have been characterizing patterns in ecological communities and trying to draw theoretical inferences form the resulting data.

Why do species occurs in specific locations?
Why do some species coexist while others do not?

Russo, Sabrina E., Stuart J. Davies, David A. King, and Sylvester Tan. “Soil-Related Performance Variation and Distributions of Tree Species in a Bornean Rain Forest.” Journal of Ecology 93, no. 5 (2005): 879–89.
Species
Species + Site information (1950s ~)
Species + Site information + Species information (2000s ~)
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
A vector of species abundance
Species composition
| Site1 | |
|---|---|
| Sp. 1 | 4 |
| Sp. 2 | 300 |
| Sp. 3 | 56 |
| Sp. 4 | 23 |

[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
Which community is more diverse?
Species richness = 2
What is the chance to get the same species?
A: \(\frac{9}{10} \times \frac{8}{9} + \frac{1}{10} \times \frac{0}{9} = 0.8\)
B: \(\frac{5}{10} \times \frac{4}{9} + \frac{5}{10} \times \frac{4}{9} \simeq 0.44\)
Which community is more diverse?

A: \(\frac{9}{10} \times \frac{8}{9} + \frac{1}{10} \times \frac{0}{9} = 0.8\)
B: \(\frac{5}{10} \times \frac{4}{9} + \frac{5}{10} \times \frac{4}{9} \simeq 0.44\)
We prefer that large values indicate more diverse communities.
Diversity of A: 1 - 0.8 = 0.2
Diversity of B: 1 - 0.44 = 0.56
Simpson’s Index of Diversity: \(D = 1 - \Sigma\frac{n_i(n_i - 1)}{N_i(N_i - 1)}\)
Simpson’s Index of Diversity (ver. 2): \(D = 1 - \Sigma p_i^2\)
n: number of individuals of each species, N: total number of individuals of all species, p: relative species abundance
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
Another simple way to describe diversity?

A: \(p_1\) = 0.9, \(p_2\) = 0.1
B: \(p_1\) = 0.5, \(p_2\) = 0.5
Diversity of A: 0.9 \(\times\) 0.1 = 0.09?
Diversity of B: 0.5 \(\times\) 0.5 = 0.25?
Diversity \(\times\) Diversity? What is the unit?
p: relative species abundance
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Sp. 1 | 4 | 0 | 315 | 23 |
| Sp. 2 | 300 | 250 | 0 | 18 |
| Sp. 3 | 56 | 120 | 74 | 0 |
| Sp. 4 | 23 | 18 | 101 | 0 |
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Site 1 | 0.00 | 0.16 | 0.81 | 0.90 |
| Site 2 | 0.16 | 0.00 | 0.79 | 0.92 |
| Site 3 | 0.81 | 0.79 | 0.00 | 0.91 |
| Site 4 | 0.90 | 0.92 | 0.91 | 0.00 |
e.g., Bray–Curtis dissimilarity
\(BC_{ij}=1-2\frac{\sum min\left(S_{A,i}\mbox{, } S_{B,i}\right)}{\sum S_{A,i}+\sum S_{B,i}}\)
Site 1 vs Site 2: 1 - (2 * (0 + 250 + 56 + 18) / (4 +300 + 56 + 23 + 0 + 250 + 120 + 18)) = 0.16
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Abundance | ||||
| Sp. 1 | 4 | 0 | 315 | 23 |
| Sp. 2 | 300 | 250 | 0 | 18 |
| Sp. 3 | 56 | 120 | 74 | 0 |
| Sp. 4 | 23 | 18 | 101 | 0 |
| Env | ||||
| Env. 1 | 780 | 2500 | 480 | 1200 |
| Env. 2 | 21 | 11 | 24 | 19 |
| Env. 3 | 1500 | 1900 | 700 | 4500 |
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Abundance | ||||
| Sp. 1 | 4 | 0 | 315 | 23 |
| Sp. 2 | 300 | 250 | 0 | 18 |
| Sp. 3 | 56 | 120 | 74 | 0 |
| Sp. 4 | 23 | 18 | 101 | 0 |
| Env | ||||
| Elevation (m) | 780 | 2500 | 480 | 1200 |
| MAT (℃) | 21 | 11 | 24 | 19 |
| MAP (mm) | 1500 | 1900 | 700 | 4500 |
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Abundance | ||||
| Sp. 1 | 4 | 0 | 315 | 23 |
| Sp. 2 | 300 | 250 | 0 | 18 |
| Sp. 3 | 56 | 120 | 74 | 0 |
| Sp. 4 | 23 | 18 | 101 | 0 |
| Env | ||||
| Elevation (m) | 780 | 2500 | 480 | 1200 |
| MAT (℃) | 21 | 11 | 24 | 19 |
| MAP (mm) | 1500 | 1900 | 700 | 4500 |
| Trait 1 | Trait 2 | Trait 3 | Trait 4 | |
|---|---|---|---|---|
| Sp 1 | 1.3 | 4.8 | 1.8 | 30.0 |
| Sp 2 | 1.7 | 12.5 | 2.1 | 22.4 |
| Sp 3 | 7.0 | 5.9 | 5.7 | 11.5 |
| Sp 4 | 2.1 | 2.1 | 3.4 | 119.9 |
[1] Vellend, M. (2016). The Theory of Ecological Communities. Princeton University Press
| Site 1 | Site 2 | Site 3 | Site 4 | |
|---|---|---|---|---|
| Abundance | ||||
| Sp. 1 | 4 | 0 | 315 | 23 |
| Sp. 2 | 300 | 250 | 0 | 18 |
| Sp. 3 | 56 | 120 | 74 | 0 |
| Sp. 4 | 23 | 18 | 101 | 0 |
| Env | ||||
| Elevation (m) | 780 | 2500 | 480 | 1200 |
| MAT (℃) | 21 | 11 | 24 | 19 |
| MAP (mm) | 1500 | 1900 | 700 | 4500 |
| Leaf N | Amax | Rdark | LL | |
|---|---|---|---|---|
| Sp 1 | 1.3 | 4.8 | 1.8 | 30.0 |
| Sp 2 | 1.7 | 12.5 | 2.1 | 22.4 |
| Sp 3 | 7.0 | 5.9 | 5.7 | 11.5 |
| Sp 4 | 2.1 | 2.1 | 3.4 | 119.9 |
Trait diversity and its role in species composition
Trait composition-environment relationships
Photosynthetic rates


Community with 1 genus and 3 species (1:3)

Community with 3 genus and 3 species (3:3)

A low genus:species ratio indicates closely related species coexist.
A high genus:species ratio indicates distantly related species coexist.
Community with 3 genus and 3 species (3:3)

Community with 3 genus and 3 species (3:3)

A low genus:species ratio indicates closely related species coexist.
A high genus:species ratio indicates distantly related species coexist.
[1] Swenson, N. G. The assembly of tropical tree communities - the advances and shortcomings of phylogenetic and functional trait analyses. Ecography 36, 264–276 (2013).
[1] Dimensions of Biodiversity: National Science Foundation

In the 1990’s conservation biologists recognized the biodiversity is not only species diversity
[1] Faith D.P. (1992) Conservation evaluation and phylogenetic diversity. Biological Conservation, 61, 1-10.
[1] Faith D.P. (1992) Conservation evaluation and phylogenetic diversity. Biological Conservation, 61, 1-10.

[1] Faith D.P. (1992) Conservation evaluation and phylogenetic diversity. Biological Conservation, 61, 1-10.

[1] Faith D.P. (1992) Conservation evaluation and phylogenetic diversity. Biological Conservation, 61, 1-10.

[1] Petchey, O. L. & Gaston, K. J. Functional diversity (FD), species richness and community composition. Ecology Letters 5, 402–411 (2002).

[1] Webb, C.O., 2000. Exploring the Phylogenetic Structure of Ecological Communities: An Example for Rain Forest Trees. The American Naturalist 156, 145–155. https://doi.org/10.1086/303378
Solution for genus:species = Use phylogenetic trees to estimate the relatedness of coexisting species

Distance matrix
A B C D E
B 1
C 2 2
D 4 4 3
E 5 5 4 2
F 5 5 4 2 1
[1] \(MPD = \frac{1}{n} \sum^{n}_{i} \sum^n_j \delta_{i,j} \; i \neq j\), where \(\delta_{i, j}\) is the pairwise distance between species i and j
Let’s consider greatest possible mean pairwise node distance (MPD) for a community of 4 taxa

A B C D E
B 1
C 2 2
D 4 4 3
E 5 5 4 2
F 5 5 4 2 1
A B E
B 1
E 5 5
F 5 5 1
Greatest MPD for a community of 4 taxa: 22 / 6 pairs = 3.66 (A, B, E, F)
[1] \(MPD = \frac{1}{n} \sum^{n}_{i} \sum^n_j \delta_{i,j} \; i \neq j\), where \(\delta_{i, j}\) is the pairwise distance between species i and j
Greatest possible mean pairwise node distance for a community of 4 taxa: 22 / 6 pairs = 3.66 (A, B, E, F)
Community 1; A, B, C, D
A B C
B 1
C 2 2
D 4 4 3
MPD = (1 + 2 + 2 + 4 + 4 + 3) / 6 = 2.66
NRI = 1 - (2.66 / 3.66) = 0.273
Community 2; A, B, E, F
A B E
B 1
E 5 5
F 5 5 1
MPD = (1 + 5 + 5 + 5+ 5 + 1) / 6 = 3.66
NRI = 1 - (3.66 / 3.66) = 0
Community 1 is more phylogenetically similar.
[1] \(MNTD = \frac{1}{n} \sum^n_i min \delta_{i,j} \; i \neq j\), where \(min \delta_{i, j}\) is the minimum distance between species i and all other species in the community.
Greatest possible nearest nodal distance for a community of 4 taxa = 2 (A, C, D, F; A to C = 2, D to F = 2)

A B C D E
B 1
C 2 2
D 4 4 3
E 5 5 4 2
F 5 5 4 2 1
[1] \(MNTD = \frac{1}{n} \sum^n_i min \delta_{i,j} \; i \neq j\), where \(min \delta_{i, j}\) is the minimum distance between species i and all other species in the community.
Community 1; A, B, C, D
A B C
B 1
C 2 2
D 4 4 3
MNTD = (1 + 1 + 2 + 3) / 4 = 1.75
NTI = 1 - (1.75 / 2.0) = 0.125
Community 2; A, B, E, F
A B E
B 1
E 5 5
F 5 5 1
MNTD = (1 + 1 + 1 + 1) / 4 = 1
NTI = 1 - (1 / 2.0) = 0.5
Community 2 is more phylogenetically similar (in tips).
[1] https://en.wikipedia.org/wiki/Anolis_ecomorphs

We are assuming that related species are ecologically similar
Related species sometimes have very different traits and ecological niches (e.g., grass-bush, trunk, trunk-crown, trunk-ground and twig ecomorphs)
[1] Losos, J. B., Jackman, T. R., Larson, A., De Queiroz, K. & Rodríguez-Schettino, L. Contingency and determinism in replicated adaptive radiations of island lizards. Science 279, 2115–2118 (1998).

A: Functional dendrogram based on ecomorph
B: Phylogeny indicates frequent evolution of traits
They do not match at all (!!)
Phylogenetically similar = Functional (ecologically) similar??
[1] Blomberg, S. P., T. Garland Jr., A. R. Ives (2003) Testing for phylogenetic signal in comparative data: Behavioral traits are more labile. Evolution, 57, 717-745.


[1] Cavender‐Bares, J., Ackerly, D. D., Baum, D. A. & Bazzaz, F. A. Phylogenetic Overdispersion in Floridian Oak Communities. The American Naturalist 163, 823–843 (2004).
[1] Cavender‐Bares, J., Ackerly, D. D., Baum, D. A. & Bazzaz, F. A. Phylogenetic Overdispersion in Floridian Oak Communities. The American Naturalist 163, 823–843 (2004).
[1] Swenson, N. G. The assembly of tropical tree communities - the advances and shortcomings of phylogenetic and functional trait analyses. Ecography 36, 264–276 (2013).

Phylogeny as a proxy for the functional or ecological similarity of species.
Measuring trait data and arraying it on the phylogenetic tree to demonstrate phylogenetic signal in function so that their phylogenetically-based inferences could be supported.
Compared to simply measuring the trait dispersion, this approach is very indirect.
This approach should be avoided! (phylogeny and traits are useful to make meaningful evolutionary inferences)
Measurable properties of plants that are indicative of ecological strategies
“Hard” traits: e.g., Photosynthetic rates

“Soft” traits: e.g., LMA (leaf mass per area)

[1] Reich, P. B. et al. From tropics to tundra: Global convergence in plant functioning. PNAS 94, 13730–13734 (1997).
[2] Wright, I. J., P. B. Reich, M. Westoby et al. The worldwide leaf economics spectrum. Nature 428, 821–827 (2004).
[3] Osnas, J. L. D. et al. Global leaf trait relationships: Mass, area, and the leaf economics spectrum. Science 340, 741–744 (2013).
[4] Katabuchi, M. et al. Decomposing leaf mass into photosynthetic and structural components explains divergent patterns of trait variation within and among plant species. bioRxiv https://doi.org/10.1101/116855 (2023)

[1] McGill, B. J., Enquist, B. J., Weiher, E. & Westoby, M. Rebuilding community ecology from functional traits. Trends in Ecology and Evolution 21, 178–185 (2006).
Non-trait based statement
Trait-based statement
[1] McGill, B. J., Enquist, B. J., Weiher, E. & Westoby, M. Rebuilding community ecology from functional traits. Trends in Ecology and Evolution 21, 178–185 (2006).
[1] Cornwell, W. K., Schwilk, L. D. W. & Ackerly, D. D. A trait-based test for habitat filtering: convex hull volume. Ecology 87, 1465–71 (2006).

California woody-plant communities (43 plots, 54 species, 3 traits)
Is the trait volume of California woody-plant communities significantly less than expected by chance?
[1] Cornwell, W. K., Schwilk, L. D. W. & Ackerly, D. D. A trait-based test for habitat filtering: convex hull volume. Ecology 87, 1465–71 (2006).

Species in 40 out of 43 plots occupied less trait space than would be expected by chance
Consistent with environmental filtering
[1] Cornwell, W. K. & Ackerly, D. D. Community assembly and shifts in plant trait distributions across an environmental gradient in coastal California. Ecological Monographs 79, 109–126 (2009).

[1] Kraft, N. J. B., Valencia, R. & Ackerly, D. D. Functional Traits and Niche-Based Tree Community Assembly in an Amazonian Forest. Science 322, 580–582 (2008).

Yasuni tropical tree communities, 25ha, 625 20m x 20m quadrats, 1089 species!
Consistent with environmental filtering
A: Ridgetops have lower than expected SLA and valleys have higher
B: Seed mass shows broader distribution than expected - Limiting similarity
C: Range of SLA is smaller than expected - Environmental filtering
[1] Katabuchi, M., Kurokawa, H., Davies, S. J., Tan, S. & Nakashizuka, T. Soil resource availability shapes community trait structure in a species-rich dipterocarp forest. Journal of Ecology 100, 643–651 (2012).

| Process / Facotr | Pattern |
|---|---|
| Biotic factor | |
| Environmental filtering | Clustring |
| Abiotic factor (competition) | |
| Limiting similarity / Compitative exclusion | Overdispersion |
| Competitive hierarchy / directional compation | Clustring |
| Herbivors / Prasites / Pathogens | Overdispersion |
| Pollinator-mediated competition | Overdispersion |
| Abiotic factor (facilitation) | |
| Nurse plants | Overdispersion |
| Pollinator facilitation | Clustering |
| Stochastic process | |
| Neutral theory | Random |
[1] \(t_A\) and \(t_B\) are the functional trait values of species A and B
[2] Kunstler, G. et al. Competitive interactions between forest trees are driven by species’ trait hierarchy, not phylogenetic or functional similarity: Implications for forest community assembly. Ecology Letters 15, 831–840 (2012).

Competitive interaction strengths between species will increase with decreasing niche distance, measured as their absolute traits distance \(|t_A - t_B|\)

Competitive effects of species A on species B will increase with increasing \(t_A - t_B\).
[1] Endara, M.-J. et al. The role of plant secondary metabolites in shaping regional and local plant community assembly. Journal of Ecology 110, 34–45 (2022).

6000+ secondary metabolites from nearly 100 species in a diverse Neotropical plant clade across the whole Amazonia
More differences in their defensive chemistry than expected by chance
Plant–herbivore interactions promote species diversity
[1] López-Angulo, J., Swenson, N. G., Cavieres, L. A. & Escudero, A. Interactions between abiotic gradients determine functional and phylogenetic diversity patterns in Mediterranean-type climate mountains in the Andes. Journal of Vegetation Science 29, 245–254 (2018).

Alpine plants in the Andes
Functional dispersion in harsh environments (higher potential solar radiation)
Facilitation tends to dominate interactions when environmental harshness increases

[1] Katabuchi, M. et al. Contrasting outcomes of species- and community-level analyses of the temporal consistency of functional composition. Ecology 98, 2273–2280 (2017).
[2] Swenson, N. G., Hulshof, C. M., Katabuchi, M. & Enquist, B. J. Long-term shifts in the functional composition and diversity of a tropical dry forest: a 30-yr study. Ecological Monographs e01408 (2020) doi:10.1002/ecm.1408.

50ha Forest Dynamics Plot on Barro Colorado Island, Panama
The changes in community-weighted mean (CWM) of wood density over time seem to suggest a community-level response in climate change (i.e, drier conditions).
No clear species-level pattern: just two of 300 species account for 60% of the temporal shifts in CWM, likely due to species-specific pathogens rather than climate change response.

We may need metrics somewhere between functional traits and invasion growth rate (model parameter) to predict coexistence.
[1] Levine, J. I., An, R., Kraft, N. J. B., Pacala, S. W. & Levine, J. M. Why ecologists struggle to predict coexistence from functional traits. Trends in Ecology & Evolution (2024) doi:10.1016/j.tree.2024.10.002.
Why do we use trait and phylogenetic diversity?
We want to quantify ecological similarities and biodiversity dimensions.
We need to be very careful when we use or develop ways to map multiple dimensions of biodiversity to a few dimensions of diversity.
Trait-based approaches
Moving ecology towards more quantitative and predictive methods.
We only focused on a “snapshot” of biodiversity (predicting future biodiversity is another story)
Developing more effective metrics may be necessary, beyond simply measuring traits.
How to calculate trait and phylogenetic diversity anyway?
The slides are made using Quarto with Emi Tanaka’s CSS design and Danyang Dai’s template.
Source code for the slides can be found at mattocci27/phy-fun-div
