skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: 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
Award ID(s):
2140720 2146220
PAR ID:
10588764
Author(s) / Creator(s):
;
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
  1. 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
  2. Abstract Traits underlie organismal responses to their environment and are essential to predict community responses to environmental conditions under global change. Species differ in life‐history traits, morphometrics, diet type, reproductive characteristics and habitat utilization.Trait associations are widely analysed using phylogenetic comparative methods (PCM) to account for correlations among related species. Similarly, traits are measured for some but not all species, and missing continuous traits (e.g. growth rate) can be imputed using ‘phylogenetic trait imputation’ (PTI), based on evolutionary relatedness and trait covariance. However, PTI has not been available for categorical traits, and estimating covariance among traits without ecological constraints risks inferring implausible evolutionary mechanisms.Here, we extend previous PCM and PTI methods by (1) specifying covariance among traits as a structural equation model (SEM), and (2) incorporating associations among both continuous and categorical traits. Fitting a SEM replaces the covariance among traits with a set of linear path coefficients specifying potential evolutionary mechanisms. Estimated parameters then represent regression slopes (i.e. the average change in trait Y given an exogenous change in trait X) that can be used to calculate both direct effects (X impacts Y) and indirect effects (X impacts Z and Z impacts Y).We demonstrate phylogenetic structural‐equation mixed‐trait imputation using 33 variables representing life history, reproductive, morphological, and behavioural traits for all >32,000 described fishes worldwide. SEM coefficients suggest that one degree Celsius increase in habitat is associated with an average 3.5% increase in natural mortality (including a 1.4% indirect impact that acts via temperature effects on the growth coefficient), and an average 3.0% decrease in fecundity (via indirect impacts on maximum age and length). Cross‐validation indicates that the model explains 54%–89% of variance for withheld measurements of continuous traits and has an area under the receiver‐operator‐characteristics curve of 0.86–0.99 for categorical traits.We use imputed traits to classify all fishes into life‐history types, and confirm a phylogenetic signal in three dominant life‐history strategies in fishes. PTI using phylogenetic SEMs ensures that estimated parameters are interpretable as regression slopes, such that the inferred evolutionary relationships can be compared with long‐term evolutionary and rearing experiments. 
    more » « less
  3. Abstract Plant ecological strategies are shaped by numerous functional traits and their trade‐offs. Trait network analysis enables testing hypotheses for the shifting of trait correlation architecture across communities differing in climate and productivity.We built plant trait networks (PTNs) for 118 species within six communities across an aridity gradient, from forest to semi‐desert across the California Floristic Province, based on 34 leaf and wood functional traits, representing hydraulic and photosynthetic function, structure, economics and size. We developed hypotheses for the association of PTN parameters with climate and ecosystem properties, based on theory for the adaptation of species to low resource/stressful environments versus higher resource availability environments with greater potential niche differentiation. Thus, we hypothesized that across community PTNs, trait network connectivity (i.e., the degree that traits are intercorrelated) and network complexity (i.e., the number of trait modules, and the degree of trait integration among modules) would be lower for communities adapted to arid climates and higher for communities adapted to greater water availability, similarly to trends expected for phylogenetic diversity, functional richness and productivity. Further, within given PTNs, we hypothesized that traits would vary strongly in their network connectivity and that the traits most centrally connected within PTNs would be those with the least across‐species variation.Across communities from more arid to wetter climates, PTN architecture varied from less to more interconnected and complex, in association with functional richness, but PTN architecture was independent of phylogenetic diversity and ecosystem productivity. Within the community PTNs, traits with lower species variation were more interconnected.Synthesis. The responsiveness of PTN architecture to climate highlights how a wide range of traits contributes to physiological and ecological strategies with an architecture that varies among plant communities. Communities in more arid environments show a lower degree of phenotypic integration, consistent with lesser niche differentiation. Our study extends the usefulness of PTNs as an approach to quantify tradeoffs among multiple traits, providing connectivity and complexity parameters as tools that clarify plant environmental adaptation and patterns of trait associations that would influence species distributions, community assembly, and ecosystem resilience in response to climate change. 
    more » « less
  4. Summary Global change has accelerated local species extinctions and colonizations, often resulting in losses and gains of evolutionary lineages with unique features. Do these losses and gains occur randomly across the phylogeny?We quantified: temporal changes in plant phylogenetic diversity (PD); and the phylogenetic relatedness (PR) of lost and gained species in 2672 semi‐permanent vegetation plots in European temperate forest understories resurveyed over an average period of 40 yr.Controlling for differences in species richness, PD increased slightly over time and across plots. Moreover, lost species within plots exhibited a higher degree of PR than gained species. This implies that gained species originated from a more diverse set of evolutionary lineages than lost species. Certain lineages also lost and gained more species than expected by chance, with Ericaceae, Fabaceae, and Orchidaceae experiencing losses and Amaranthaceae, Cyperaceae, and Rosaceae showing gains. Species losses and gains displayed no significant phylogenetic signal in response to changes in macroclimatic conditions and nitrogen deposition.As anthropogenic global change intensifies, temperate forest understories experience losses and gains in specific phylogenetic branches and ecological strategies, while the overall mean PD remains relatively stable. 
    more » « less
  5. Abstract Global change is impacting biodiversity across all habitats on earth. New selection pressures from changing climatic conditions and other anthropogenic activities are creating heterogeneous ecological and evolutionary responses across many species' geographic ranges. Yet we currently lack standardised and reproducible tools to effectively predict the resulting patterns in species vulnerability to declines or range changes.We developed an informatic toolbox that integrates ecological, environmental and genomic data and analyses (environmental dissimilarity, species distribution models, landscape connectivity, neutral and adaptive genetic diversity, genotype‐environment associations and genomic offset) to estimate population vulnerability. In our toolbox, functions and data structures are coded in a standardised way so that it is applicable to any species or geographic region where appropriate data are available, for example individual or population sampling and genomic datasets (e.g. RAD‐seq, ddRAD‐seq, whole genome sequencing data) representing environmental variation across the species geographic range.To demonstrate multi‐species applicability, we apply our toolbox to three georeferenced genomic datasets for co‐occurring East African spiny reed frogs (Afrixalus fornasini, A. delicatusandA. sylvaticus) to predict their population vulnerability, as well as demonstrating that range loss projections based on adaptive variation can be accurately reproduced from a previous study using data for two European bat species (Myotis escaleraiandM. crypticus).Our framework sets the stage for large scale, multi‐species genomic datasets to be leveraged in a novel climate change vulnerability framework to quantify intraspecific differences in genetic diversity, local adaptation, range shifts and population vulnerability based on exposure, sensitivity and landscape barriers. 
    more » « less