Summary Large intraspecific functional trait variation strongly impacts many aspects of communities and ecosystems, and is the medium upon which evolution works. Yet intraspecific trait variation is inconsistent and hard to predict across traits, species and locations.We measured within‐species variation in leaf mass per area (LMA), leaf dry matter content (LDMC), branch wood density (WD), and allocation to stem area vs leaf area in branches (branch Huber value (HV)) across the aridity range of seven Australian eucalypts and a co‐occurringAcaciaspecies to explore how traits and their variances change with aridity.Within species, we found consistent increases in LMA, LDMC and WD and HV with increasing aridity, resulting in consistent trait coordination across leaves and branches. However, this coordination only emerged across sites with large climate differences. Unlike trait means, patterns of trait variance with aridity were mixed across populations and species. Only LDMC showed constrained trait variation in more xeric species and drier populations that could indicate limits to plasticity or heritable trait variation.Our results highlight that climate can drive consistent within‐species trait patterns, but that patterns might often be obscured by the complex nature of morphological traits, sampling incomplete species ranges or sampling confounded stress gradients.
more »
« less
This content will become publicly available on December 1, 2025
Extending phylogenetic regression models for comparing within‐species patterns across the tree of life
Abstract Evolutionary biologists characterize macroevolutionary trends of phenotypic change across the tree of life using phylogenetic comparative methods. However, within‐species variation can complicate such investigations. For this reason, procedures for incorporating nonstructured (random) intraspecific variation have been developed.Likewise, evolutionary biologists seek to understand microevolutionary patterns of phenotypic variation within species, such as sex‐specific differences or allometric trends. Additionally, there is a desire to compare such within‐species patterns across taxa, but current analytical approaches cannot be used to interrogate within‐species patterns while simultaneously accounting for phylogenetic non‐independence. Consequently, deciphering how intraspecific trends evolve remains a challenge.Here we introduce an extended phylogenetic generalized least squares (E‐PGLS) procedure which facilitates comparisons of within‐species patterns across species while simultaneously accounting for phylogenetic non‐independence.Our method uses an expanded phylogenetic covariance matrix, a hierarchical linear model, and permutation methods to obtain empirical sampling distributions and effect sizes for model effects that can evaluate differences in intraspecific trends across species for both univariate and multivariate data, while conditioning them on the phylogeny.The method has appropriate statistical properties for both balanced and imbalanced data. Additionally, the procedure obtains evolutionary covariance estimates that reflect those from existing approaches for nonstructured intraspecific variation. Importantly, E‐PGLS can detect differences in structured (i.e. microevolutionary) intraspecific patterns across species when such trends are present. Thus, E‐PGLS extends the reach of phylogenetic comparative methods into the intraspecific comparative realm, by providing the ability to compare within‐species trends across species while simultaneously accounting for shared evolutionary history.
more »
« less
- PAR ID:
- 10588764
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Methods in Ecology and Evolution
- Volume:
- 15
- Issue:
- 12
- ISSN:
- 2041-210X
- Page Range / eLocation ID:
- 2234 to 2246
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Summary Vulnerability to embolism varies between con‐generic species distributed along aridity gradients, yet little is known about intraspecific variation and its drivers. Even less is known about intraspecific variation in tissues other than stems, despite results suggesting that roots, stems and leaves can differ in vulnerability. We hypothesized that intraspecific variation in vulnerability in leaves and stems is adaptive and driven by aridity.We quantified leaf and stem vulnerability ofQuercus douglasiiusing the optical technique. To assess contributions of genetic variation and phenotypic plasticity to within‐species variation, we quantified the vulnerability of individuals growing in a common garden, but originating from populations along an aridity gradient, as well as individuals from the same wild populations.Intraspecific variation in water potential at which 50% of total embolism in a tissue is observed (P50) was explained mostly by differences between individuals (>66% of total variance) and tissues (16%). There was little between‐population variation in leaf/stem P50in the garden, which was not related to site of origin aridity. Unexpectedly, we observed a positive relationship between wild individual stem P50and aridity.Although there is no local adaptation and only minor phenotypic plasticity in leaf/stem vulnerability inQ. douglasii, high levels of potentially heritable variation within populations or strong environmental selection could contribute to adaptive responses under future climate change.more » « less
-
Abstract Gymnosperms encompass a diverse group of mostly woody plants with high ecological and economic value, yet little is known about the scope and organization of fine‐root trait diversity among gymnosperms due to the undersampling of most gymnosperm families and the dominance of angiosperm groups in recent syntheses.New and existing data were compiled for morphological traits (root diameter, length, tissue density, specific root length [SRL] and specific root area [SRA]), the architectural trait branching ratio, root nitrogen content [N] and mycorrhizal colonization. We used phylogenetic least squares regression and principal component analysis to determine trait–trait relationships and coordination across 66 species, representing 11 of the 12 extant gymnosperm families from boreal, temperate, subtropical and tropical biomes. Finally, we compared the relationship between family divergence time and mean trait values to determine whether evolutionary history structured variation in fine‐root traits within the gymnosperm phylogeny.Wide variation in gymnosperm root traits could be largely captured by two primary axes of variation defined by SRL and diameter, and root tissue density and root nitrogen, respectively. However, individual root length and SRA each had significant correlations with traits defining both main axes of variation. Neither mycorrhizal colonization nor root branching ratio were closely related to other traits. We did not observe a directional evolution of mean trait values from older to more recently diverged gymnosperm families.Synthesis. Despite their unique evolutionary history, gymnosperms display a root economic space similar to that identified in angiosperms, likely reflecting common constraints on plants adapting to diverse environments in both groups. These findings provide greater confidence that patterns observed in broad syntheses justly capture patterns of trait diversity among multiple, distinct lineages. Additionally, independence between root architecture and other traits may support greater diversity in below‐ground resource acquisition strategies. Unlike angiosperms, there were no clear trends towards increasingly thin roots over evolutionary time, possibly because of lower diversification rates or unique biogeographic history among gymnosperms, though additional observations are needed to more richly test evolutionary trends among gymnosperms.more » « less
-
Abstract It has become common in evolutionary biology to characterize phenotypes multivariately. However, visualizing macroevolutionary trends in multivariate datasets requires appropriate ordination methods.In this paper we describe phylogenetically aligned component analysis (PACA): a new ordination approach that aligns phenotypic data with phylogenetic signal. Unlike phylogenetic principal component analysis (Phy‐PCA), which finds an alignment of a principal eigenvector that is independent of phylogenetic signal, PACA maximizes variation in directions that describe phylogenetic signal, while simultaneously preserving the Euclidean distances among observations in the data space.We demonstrate with simulated and empirical examples that with PACA, it is possible to visualize the trend in phylogenetic signal in multivariate data spaces, irrespective of other signals in the data. In conjunction with Phy‐PCA, one can visualize both phylogenetic signal and trends in data independent of phylogenetic signal.Phylogenetically aligned component analysis can distinguish between weak phylogenetic signals and strong signals concentrated in only a portion of all data dimensions. We provide empirical examples that emphasize the difference. Use of PACA in studies focused on phylogenetic signal should enable much more precise description of the phylogenetic signal, as a result.Overall, PACA will return a projection that shows the most phylogenetic signal in the first few components, irrespective of other signals in the data. By comparing Phy‐PCA and PACA results, one may glean the relative importance of phylogenetic and other (ecological) signals in the data.more » « less
-
Abstract Understanding how vegetation responds to drought is fundamental for understanding the broader implications of climate change on foundation tree species that support high biodiversity. Leveraging remote sensing technology provides a unique vantage point to explore these responses across and within species.We investigated interspecific drought responses of twoPopulusspecies (P.fremontii,P.angustifolia) and their naturally occurring hybrids using leaf‐level visible through shortwave infrared (VSWIR; 400–2500 nm) reflectance. AsF1hybrids backcross with either species, resulting in a range of backcross genotypes, we heretofore refer to the two species and their hybrids collectively as ‘cross types’. We additionally explored intraspecific variation inP. fremontiidrought response at the leaf and canopy levels using reflectance data and thermal unmanned aerial vehicle (UAV) imagery. We employed several analyses to assess genotype‐by‐environment (G × E) interactions concerning drought, including principal component analysis, support vector machine and spectral similarity index.Five key findings emerged: (1) Spectra of all three cross types shifted significantly in response to drought. The magnitude of these reaction norms can be ranked from hybrids>P. fremontii>P. angustifolia, suggesting differential variation in response to drought; (2) Spectral space among cross types constricted under drought, indicating spectral—and phenotypic—convergence; (3) Experimentally, populations ofP. fremontiifrom cool regions had different responses to drought than populations from warm regions, with source population mean annual temperature driving the magnitude and direction of change in VSWIR reflectance. (4) UAV thermal imagery revealed that watered, warm‐adapted populations maintained lower leaf temperatures and retained more leaves than cool‐adapted populations, but differences in leaf retention decreased when droughted. (5) These findings are consistent with patterns of local adaptation to drought and temperature stress, demonstrating the ability of leaf spectra to detect ecological and evolutionary responses to drought as a function of adaptation to different environments.Synthesis.Leaf‐level spectroscopy and canopy‐level UAV thermal data captured inter‐ and intraspecific responses to water stress in cottonwoods, which are widely distributed in arid environments. This study demonstrates the potential of remote sensing to monitor and predict the impacts of drought on scales varying from leaves to landscapes.more » « less
An official website of the United States government
