Functional trait heritability and local climatic adaptation among grasses: a meta-analysis

Variation in climate has been demonstrated to be a powerful driver of selection and local adaptation among plant populations. Variation in functional traits among populations can also be indicative of the drivers of local adaptation. However, it is not clear to what extent species exhibit consistent patterns of local adaptation as revealed by common, heritable trait–environment relationships among populations. To address this, we conducted a meta-analysis of common garden studies of grass populations to quantify the degree of heritability of several commonly measured functional traits, and whether demonstrated heritability was driven by climate. We found that leaf size, specific leaf area (SLA) and total biomass all displayed strong broad-sense of heritability. Both leaf area and SLA decreased significantly with increasing temperature seasonality among populations within species, while total biomass increased with increasing annual and dry season precipitation, and decreased with increasing precipitation seasonality. These results indicate similar, consistent drivers of local adaptation among species of grasses. Further information on trait–environment relationships within species could greatly improve our ability to predict broad scale patterns in functional diversity across multiple levels of ecological organization. Expanding the range of traits and regions incorporated in common garden research, in the present case by incorporating root traits and Southern Hemisphere taxa, will provide even greater benefits to the fields of restoration, conservation, and global change ecology.     

Cumulative weather effects can impact across the whole life cycle

Predicting how species will be affected by future climatic change requires the underlying environmental drivers to be identified. As vital rates vary over the lifecycle, structured population models derived from statistical environment–demography relationships are often used to inform such predictions. Environmental drivers are typically identified independently for different vital rates and demographic classes. However, these rates often exhibit positive temporal covariance, suggesting that vital rates respond to common environmental drivers. Additionally, models often only incorporate average weather conditions during a single, a priori chosen time window (e.g. monthly means). Mismatches between these windows and the period when the vital rates are sensitive to variation in climate decrease the predictive performance of such approaches. We used a demographic structural equation model (SEM) to demonstrate that a single axis of environmental variation drives the majority of the (co)variation in survival, reproduction, and twinning across six age–sex classes in a Soay sheep population. This axis provides a simple target for the complex task of identifying the drivers of vital rate variation. We used functional linear models (FLMs) to determine the critical windows of three local climatic drivers, allowing the magnitude and direction of the climate effects to differ over time. Previously unidentified lagged climatic effects were detected in this well‐studied population. The FLMs had a better predictive performance than selecting a critical window a priori, but not than a large‐scale climate index. Positive covariance amongst vital rates and temporal variation in the effects of environmental drivers are common, suggesting our SEM–FLM approach is a widely applicable tool for exploring the joint responses of vital rates to environmental change.     

Does environmental heterogeneity drive functional trait variation? A test in montane and alpine meadows

While community‐weighted means of plant traits have been linked to mean environmental conditions at large scales, the drivers of trait variation within communities are not well understood. Local environmental heterogeneity (such as microclimate variability), in addition to mean environmental conditions, may decrease the strength of environmental filtering and explain why communities support different amounts of trait variation. Here, we assess two hypotheses: first, that more heterogeneous local environments and second, that less extreme environments, should support a broader range of plant strategies and thus higher trait variation. We quantified drivers of trait variation across a range of environmental conditions and spatial scales ranging from sub‐meter to tens of kilometers in montane and alpine plant communities. We found that, within communities, both environmental heterogeneity and environmental means are drivers of trait variation. However, the importance of each environmental factor varied depending on the trait. Our results indicate that larger‐scale trait–climate linkages that hold across communities also apply at small spatial scales, suggesting that microclimate variation within communities is a key driver of community functional diversity. Microclimatic variation provides a potential mechanism for helping to maintain diversity in local communities and also suggests that small‐scale environmental heterogeneity should be measured as a better predictor of functional diversity.     

Correlated evolution between climate and suites of traits along a fast–slow continuum in the radiation of Protea

Evolutionary radiations are responsible for much of Earth's diversity, yet the causes of these radiations are often elusive. Determining the relative roles of adaptation and geographic isolation in diversification is vital to understanding the causes of any radiation, and whether a radiation may be labeled as “adaptive” or not. Across many groups of plants, trait–climate relationships suggest that traits are an important indicator of how plants adapt to different climates. In particular, analyses of plant functional traits in global databases suggest that there is an “economics spectrum” along which combinations of functional traits covary along a fast–slow continuum. We examine evolutionary associations among traits and between trait and climate variables on a strongly supported phylogeny in the iconic plant genus Protea to identify correlated evolution of functional traits and the climatic‐niches that species occupy. Results indicate that trait diversification in Protea has climate associations along two axes of variation: correlated evolution of plant size with temperature and leaf investment with rainfall. Evidence suggests that traits and climatic‐niches evolve in similar ways, although some of these associations are inconsistent with global patterns on a broader phylogenetic scale. When combined with previous experimental work suggesting that trait–climate associations are adaptive in Protea, the results presented here suggest that trait diversification in this radiation is adaptive.     

Plasticity of functional traits varies clinally along a rainfall gradient in Eucalyptus tricarpa

Widespread species often occur across a range of climatic conditions, through a combination of local genetic adaptations and phenotypic plasticity. Species with greater phenotypic plasticity are likely to be better positioned to cope with rapid anthropogenic climate changes, while those displaying strong local adaptations might benefit from translocations to assist the movement of adaptive genes as the climate changes. Eucalyptus tricarpa occurs across a climatic gradient in south‐eastern Australia, a region of increasing aridity, and we hypothesized that this species would display local adaptation to climate. We measured morphological and physiological traits reflecting climate responses in nine provenances from sites of 460 to 1040 mm annual rainfall, in their natural habitat and in common gardens near each end of the gradient. Local adaptation was evident in functional traits and differential growth rates in the common gardens. Some traits displayed complex combinations of plasticity and genetic divergence among provenances, including clinal variation in plasticity itself. Provenances from drier locations were more plastic in leaf thickness, whereas leaf size was more plastic in provenances from higher rainfall locations. Leaf density and stomatal physiology (as indicated by δ¹³C and δ¹⁸O) were highly and uniformly plastic. In addition to variation in mean trait values, genetic variation in trait plasticity may play a role in climate adaptation.     

MULTIVARIATE ADAPTATION BUT NO INCREASE IN COMPETITIVE ABILITY IN INVASIVE GERANIUM CAROLINIANUM L. (GERANIACEAE)

Adaptive evolution can affect the successful establishment of invasive species, but changes in selective pressures, loss of genetic variation in relevant traits, and/or altered trait correlations can make adaptation difficult to predict. We used a common‐garden experiment to assess trait correlations and patterns of adaptation in the invasive plant, Geranium carolinianum, sampled across 20 populations in its native (United States) and invasive (China) ranges. We used multivariate Q_ST – F_ST tests to determine if phenotypic differences between countries are attributable to adaptation. We also compared population‐level variation within each country to assess whether local adaptation resulted in similar multivariate phenotypes in the United States and China. Between countries, most phenotypic differences are indistinguishable from genetic drift, although we detected a signature of adaptation to the colder, drier winters in China. There was no evidence for increases in invasive traits in China. Within countries, strong multivariate adaptation appears to be driven by latitudinal climatic variation in the United States, but not in China. Additionally, adaptive trait combinations as well as their underlying correlations differ between the two countries, indicating that adaptation in invasive populations does not parallel patterns in native populations due to differences in selection pressures, genetic constraints, or both.     

青藏高原东缘常见阔叶木本植物叶片性状对环境因子的响应

叶片性状-环境关系对于预测气候变化对植物的影响至关重要。该研究以青藏高原东缘常见阔叶木本植物为研究对象, 从47个样点采集了332个物种共666个种群的叶片, 测量了15个叶片性状, 调查了该区域木本植物叶片性状的变异程度, 并从种内和种间水平探讨了叶片性状对环境的响应及适应策略。结果表明, 反眏叶片大小的性状均具有较高的变异, 其中, 叶片面积是变异程度最大的性状。除气孔密度外, 大多数叶片性状与海拔显著相关。气候是叶片性状变异的重要驱动因素, 3.3%-29.5%的叶片性状变异由气候因子组合解释。其中, 气温对叶片性状变异解释度最高, 日照时间能解释大部分叶片性状的变异, 而降水量对叶片性状变异的解释度相对较小。与环境(海拔和气候因子)显著相关的叶片性状在种内明显少于种间水平, 可能是植物性状之间的协同变化与权衡使种内性状变异比较小, 从而减弱了种内叶片性状与环境因子的相关性。研究结果总体表明,叶片性状与木本植物对环境的适应策略密切相关, 植物通过选择小而厚的叶片和较短的叶柄以适应高海拔的 环境。     

Trait variation and integration across scales: is the leaf economic spectrum present at local scales?

Trait‐based approaches have taken an increasingly dominant role in community ecology. Although trait‐based strategy dimensions such as the leaf economic spectrum (LES) have been identified primarily at global‐scales, trait variation at the community scale is often interpreted in this context. Here we argue from several lines of evidence that a research priority should be to determine whether global‐scale trait relationships hold at more local scales. We review recent literature assessing trait variation at smaller scales, and then present a case study exploring the relationship between the correlation strength of leaf traits and their similarity in variation structure across ecological scales. We find that the correlation strength between pairs of leaf traits does not predict whether the traits respond similarly to different drivers of variation. Instead, correlation strength only sets an upper bound to the dissimilarity in trait variation structure. With moderate correlation strengths, LES traits largely retain the ability to respond independently to different drivers of phenotypic variation at different scales. Recent literature and our results suggest that LES relationships may not hold at local scales. Clarifying under what conditions and at which scales the LES is consistently expressed is necessary for us to make the most of the emerging trait toolbox.     

The interplay among intraspecific leaf trait variation,niche breadth and species abundance along light and soil nutrient gradients

It is assumed that widespread, generalist species have high phenotypic variation, but we know little about how intraspecific trait variation (ITV) relates to species abundance and niche breadth. In the temperate rainforest of southern Chile, we hypothesized that species with wide niche breadth would exhibit 1) high among‐plot ITV, 2) a strong relationship between trait values and the environment, and 3) a close fit between traits and local environment trait optima. We measured leaf functional traits (leaf area, LMA, leaf N and P concentrations) of saplings in woody species, and compared the relative abundance of each species with its niche breadth, measured as the range of light, soil N and P availability. We used the slope of the linear regression of species’ trait–environment relationships to assess the strength and direction of these relationships, and measured the degree to which species’ trait values track the environmental optimum across plots. In some cases, species having wide niche breadth had high ITV in leaf N and also matched traits (LMA and leaf P) to local optima along the light gradient; they also had high ITV in general and matched leaf P to local optima along the soil P gradient. The relationship between species with wide niche breadth and the strength of intraspecific trait–environment relationships was generally weak and varied depending on the niche dimension and trait in question. Species varied considerably in the strength of trait–environment relationships and total magnitude of ITV, and this variation was not generally strongly related to species abundances or niche breadth patterns. In conclusion, trait variation at the community level is not driven by a few abundant, widely distributed species, but depends on the aggregate trait responses of both abundant and rare species. This makes it difficult to scale individual species trait responses up to the community level.     

Variation in plant functional traits across and within four species of Western Australian Banksia (Proteaceae) along a natural climate gradient

Plant traits are fundamental components of the ecological strategies of plants, relating to how plants acquire and use resources. Their study provides insight into the dynamics of species geographical ranges in changing environments. Here, we assessed the variation in trait values at contrasting points along an environmental gradient to provide insight into the flexibility of species response to environmental heterogeneity. Firstly, we identified how commonly measured functional traits of four congeneric species (Banksia baxteri, B. coccinea, B. media and B. quercifolia) varied along a longitudinal gradient in the South Western Australian Floristic Region. This regional gradient provides significant variation in moisture, temperature and soil nutrients: soil nitrogen content decreases with declining rainfall and increasing temperature. We hypothesized that (i) the regional pattern in trait–environment associations across the species would match those observed on a global scale and (ii) that the direction and slopes of the within‐species relationships would be similar to those across species for each of the measured traits. Along the regional gradient we observed strong shifts in trait values, and cross‐species relationships followed the expected trend: specific leaf area was significantly lower, and leaf N_area and seed dry mass significantly higher, at the drier end of the rainfall gradient. However, traits within species were generally not well correlated with habitat factors: we found weak patterns among populations, either due to the small ecological gradient or perhaps because fine‐scale structuring among populations (at a micro‐evolutionary scale) was low due to high gene flow within species. Understanding how species traits shift as a result of climatic influences, both at the regional (across species) and local (within species) scale, provides insight into plant adaptation to the environment. Such studies have important applications for conservation biology and population management in the face of global change.     

Intraspecific phenotypic variability of plant functional traits in contrasting mountain grasslands habitats

Empirical studies that link plants intraspecific variation to environmental conditions are almost lacking, despite their relevance in understanding mechanisms of plant adaptation, in predicting the outcome of environmental change and in conservation. Here, we investigate intraspecific trait variation of four grassland species along with abiotic environmental variation at high spatial resolution (n = 30 samples per species trait and environmental factor per site) in two contrasting grassland habitats in Central Apennines (Italy). We test for phenotypic adaptation between habitats, intraspecific trait-environment relationships within habitats, and the extent of trait and environmental variation. We considered whole plant, clonal, leaf, and seed traits. Differences between habitats were tested using ANOVA and ANCOVA. Trait-environment relationships were assessed using multiple regression models and hierarchical variance partitioning. The extent of variation was calculated using the coefficient of variation. Significant intraspecific differences in trait attributes between the contrasting habitats indicate phenotypic adaptation to in situ environmental conditions. Within habitats, light, soil temperature, and the availability of nitrate, ammonium, magnesium and potassium were the most important factors driving intraspecific trait-environment relationships. Leaf traits and height growth show lower variability than environment being probably more regulated by plants than clonal traits which show much higher variability. We show the adaptive significance of key plant traits leading to intraspecific adaptation of strategies providing insights for conservation of extant grassland communities. We argue that protecting habitats with considerable medium- and small-scale environmental heterogeneity is important to maintain large intraspecific variability within local populations that finally can buffer against uncertainty of future climate and land use scenarios.     

Adaptation to climate through flowering phenology: a case study in Medicago truncatula

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness‐related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome‐wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.     

Linking trait variation to the environment: critical issues with community‐weighted mean correlation resolved by the fourth‐corner approach

Establishing trait–environment relationships has become routine in community ecology. Here, we demonstrate that the community weighted means correlation (CWM) and its parallel approach in linking trait variation to the environment, the species niche centroid correlation (SNC), have important shortcomings, arguing against their continuing application. Using mathematical derivations and simulations, we show that the two major issues are inconsistent parameter estimation and unacceptable significance rates when only the environment or only traits are structuring species distributions, but they themselves are not linked. We show how both CWM and SNC are related to the fourth‐corner correlation and propose to replace all by the Chessel fourth‐corner correlation, which is the fourth‐corner correlation divided by its maximum attainable value. We propose an appropriate hypothesis testing procedure that is not only unbiased but also has much greater statistical power in detecting trait–environmental relationships. We derive an additive framework in which trait variation is partitioned among and within communities, which can be then modeled against the environment. We finish by presenting a contrast between methods and an application of our proposed framework across 85 lake‐fish metacommunities.     

Ecological genomics meets community‐level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation

Local adaptation is a central feature of most species occupying spatially heterogeneous environments, and may factor critically in responses to environmental change. However, most efforts to model the response of species to climate change ignore intraspecific variation due to local adaptation. Here, we present a new perspective on spatial modelling of organism–environment relationships that combines genomic data and community‐level modelling to develop scenarios regarding the geographic distribution of genomic variation in response to environmental change. Rather than modelling species within communities, we use these techniques to model large numbers of loci across genomes. Using balsam poplar (Populus balsamifera) as a case study, we demonstrate how our framework can accommodate nonlinear responses of loci to environmental gradients. We identify a threshold response to temperature in the circadian clock gene GIGANTEA‐5 (GI5), suggesting that this gene has experienced strong local adaptation to temperature. We also demonstrate how these methods can map ecological adaptation from genomic data, including the identification of predicted differences in the genetic composition of populations under current and future climates. Community‐level modelling of genomic variation represents an important advance in landscape genomics and spatial modelling of biodiversity that moves beyond species‐level assessments of climate change vulnerability.     

Characterizing the contribution of plasticity and genetic differentiation to community‐level trait responses to environmental change

The match between functional trait variation in communities and environmental gradients is maintained by three processes: phenotypic plasticity and genetic differentiation (intraspecific processes), and species turnover (interspecific). Recently, evidence has emerged suggesting that intraspecific variation might have a potentially large role in driving functional community composition and response to environmental change. However, empirical evidence quantifying the respective importance of phenotypic plasticity and genetic differentiation relative to species turnover is still lacking. We performed a reciprocal transplant experiment using a common herbaceous plant species (Oxalis montana) among low‐, mid‐, and high‐elevation sites to first quantify the contributions of plasticity and genetic differentiation in driving intraspecific variation in three traits: height, specific leaf area, and leaf area. We next compared the contributions of these intraspecific drivers of community trait–environment matching to that of species turnover, which had been previously assessed along the same elevational gradient. Plasticity was the dominant driver of intraspecific trait variation across elevation in all traits, with only a small contribution of genetic differentiation among populations. Local adaptation was not detected to a major extent along the gradient. Fitness components were greatest in O. montana plants with trait values closest to the local community‐weighted means, thus supporting the common assumption that community‐weighted mean trait values represent selective optima. Our results suggest that community‐level trait responses to ongoing climate change should be mostly mediated by species turnover, even at the small spatial scale of our study, with an especially small contribution of evolutionary adaptation within species.     

Defining a spectrum of integrative trait‐based vegetation canopy structural types

Vegetation canopy structure is a fundamental characteristic of terrestrial ecosystems that defines vegetation types and drives ecosystem functioning. We use the multivariate structural trait composition of vegetation canopies to classify ecosystems within a global canopy structure spectrum. Across the temperate forest sub‐set of this spectrum, we assess gradients in canopy structural traits, characterise canopy structural types (CST) and evaluate drivers and functional consequences of canopy structural variation. We derive CSTs from multivariate canopy structure data, illustrating variation along three primary structural axes and resolution into six largely distinct and functionally relevant CSTs. Our results illustrate that within‐ecosystem successional processes and disturbance legacies can produce variation in canopy structure similar to that associated with sub‐continental variation in forest types and eco‐climatic zones. The potential to classify ecosystems into CSTs based on suites of structural traits represents an important advance in understanding and modelling structure–function relationships in vegetated ecosystems.     

Native drivers of fish life history traits are lost during the invasion process

Rapid adaptation to global change can counter vulnerability of species to population declines and extinction. Theoretically, under such circumstances both genetic variation and phenotypic plasticity can maintain population fitness, but empirical support for this is currently limited. Here, we aim to characterize the role of environmental and genetic diversity, and their prior evolutionary history (via haplogroup profiles) in shaping patterns of life history traits during biological invasion. Data were derived from both genetic and life history traits including a morphological analysis of 29 native and invasive populations of topmouth gudgeon Pseudorasbora parva coupled with climatic variables from each location. General additive models were constructed to explain distribution of somatic growth rate (SGR) data across native and invasive ranges, with model selection performed using Akaike's information criteria. Genetic and environmental drivers that structured the life history of populations in their native range were less influential in their invasive populations. For some vertebrates at least, fitness‐related trait shifts do not seem to be dependent on the level of genetic diversity or haplogroup makeup of the initial introduced propagule, nor of the availability of local environmental conditions being similar to those experienced in their native range. As long as local conditions are not beyond the species physiological threshold, its local establishment and invasive potential are likely to be determined by local drivers, such as density‐dependent effects linked to resource availability or to local biotic resistance.     

Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species

The wild currant tomato Solanum pimpinellifolium inhabits a wide range of abiotic habitats across its native range of Ecuador and Peru. Although it has served as a key genetic resource for the improvement of domestic cultivars, little is known about the genetic basis of traits underlying local adaptation in this species, nor what abiotic variables are most important for driving differentiation. Here we use redundancy analysis (RDA) and other multivariate statistical methods (structural equation modelling [SEM] and generalized dissimilarity modelling [GDM]) to quantify the relationship of genomic variation (6,830 single nucleotide polymorphisms [SNPs]) with climate and geography, among 140 wild accessions. RDA, SEM and GDM each identified environment as explaining more genomic variation than geography, suggesting that local adaptation to heterogeneous abiotic habitats may be an important source of genetic diversity in this species. Environmental factors describing temporal variation in precipitation and evaporative demand explained the most SNP variation among accessions, indicating that these forces may represent key selective agents. Lastly, by studying how SNP–environment associations vary throughout the genome (44,064 SNPs), we mapped the location and investigated the functions of loci putatively contributing to climatic adaptations. Together, our findings indicate an important role for selection imposed by the abiotic environment in driving genomic differentiation between populations.     

LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.