Assessing adaptive and plastic responses in growth and functional traits in a 10-year-old common garden experiment with pedunculate oak (Quercus robur L.) suggests that directional selection can drive climatic adaptation

Understanding how tree species will respond to a future climate requires reliable and quantitative estimates of intra-specific variation under current climate conditions. We studied three 10-year-old common garden experiments established across a rainfall and drought gradient planted with nearly 10,000 pedunculate oak (Quercus robur L.) trees from ten provenances with known family structure. We aimed at disentangling adaptive and plastic responses for growth (height and diameter at breast height) as well as for leaf and wood functional traits related to adaptation to dry environments. We used restricted maximum likelihood approaches to assess additive genetic variation expressed as narrow-sense heritability (h2), quantitative trait differentiation among provenances (QST), and genotype-by-environment interactions (GxE). We found strong and significant patterns of local adaptation in growth in all three common gardens, suggesting that transfer of seed material should not exceed a climatic distance of approximately 1°C under current climatic conditions, while transfer along precipitation gradients seems to be less stringent. Moreover, heritability reached 0.64 for tree height and 0.67 for dbh at the dry margin of the testing spectrum, suggesting significant additive genetic variation of potential use for future selection and tree breeding. GxE interactions in growth were significant and explained less phenotypic variation than origin of seed source (4% versus 10%). Functional trait variation among provenances was partly related to drought regimes at provenances origins but had moderate explanatory power for growth. We conclude that directional selection, either naturally or through breeding, is the most likely and feasible outcome for pedunculate oak to adapt to warmer and drier climate conditions in the future.

Intraspecific responses to climate reveal nonintuitive warming impacts on a widespread thermophilic conifer

Many ecologically important forest trees from dry areas have been insufficiently investigated for their ability to adapt to the challenges posed by climate change, which hampers the implementation of mitigation policies. We analyzed 14 common-garden experiments across the Mediterranean which studied the widespread thermophilic conifer Pinus halepensis and involved 157 populations categorized into five ecotypes. Ecotype-specific tree height responses to climate were applied to projected climate change (2071-2100 ad), to project potential growth patterns both locally and across the species' range. We found contrasting ecotypic sensitivities to annual precipitation but comparatively uniform responses to mean temperature, while evidence of local adaptation for tree height was limited to mesic ecotypes. We projected intriguing patterns of response range-wide, implying either height inhibition or stimulation of up to 75%, and deduced that the ecotype currently experiencing more favorable (wetter) conditions will show the largest inhibition. Extensive height reductions can be expected for coastal areas of France, Greece, Spain and northern Africa. Our findings underline the fact that intraspecific variations in sensitivity to precipitation must be considered when projecting tree height responses of dry forests to future climate. The ecotype-specific projected performances call for management activities to ensure forest resilience in the Mediterranean through, for example, tailored deployment strategies.

Adaptation to drought is coupled with slow growth, but independent from phenology in marginal silver fir (Abies alba Mill.) populations

Drought is one of the most important selection pressures for forest trees in the context of climate change. Yet, the different evolutionary mechanisms, and their environmental drivers, by which certain populations become more drought tolerant than others is still little understood. We studied adaptation to drought in 16 silver fir (Abies alba Mill.) populations from the French Mediterranean Alps by combining observations on seedlings from a greenhouse experiment (N = 8,199) and on adult tress in situ (N = 315). In the greenhouse, we followed half-sib families for four growing seasons for growth and phenology traits, and tested their water stress response in a "drought until death" experiment. Adult trees in the field were assessed for δ 13C, a proxy for water use efficiency, and genotyped at 357 SNP loci. SNP data was used to generate a null expectation for seedling trait divergence between populations in order to detect the signature of selection, and 31 environmental variables were used to identify the selective environment. We found that seedlings originating from populations with low soil water capacity grew more slowly, attained a smaller stature, and resisted water stress for a longer period of time in the greenhouse. Additionally, adult trees of these populations exhibited a higher water use efficiency as evidenced by their δ 13C. These results suggest a correlated evolution of the growth-drought tolerance trait complex. Population divergence in bud break phenology was adaptive only in the second growing season, and evolved independently from the growth-drought tolerance trait complex. Adaptive divergence in bud break phenology was principally driven by the inter- and intra-annual variation in temperature at the geographic origin of the population. Our results illustrate the different evolutionary strategies used by populations to cope with drought stress at the range limits across a highly heterogeneous landscape, and can be used to inform assisted migration programs.

Contrasted patterns of local adaptation to climate change across the range of an evergreen oak, Quercus aquifolioides

Long-lived tree species are genetically differentiated and locally adapted with respect to fitness-related traits, but the genetic basis of local adaptation remains largely unresolved. Recent advances in population genetics and landscape genomic analyses enable identification of putative adaptive loci and specific selective pressures acting on local adaptation. Here, we sampled 60 evergreen oak (Quercus aquifolioides) populations throughout the species' range and pool-sequenced 587 individuals at drought-stress candidate genes. We analyzed patterns of genetic diversity and differentiation for 381 single nucleotide polymorphisms (SNPs) from 65 candidate genes and eight microsatellites. Outlier loci were identified by genetic differentiation analysis and genome-environment associations. The response pattern of genetic variation to environmental gradient was assessed by linear isolation-by-distance/environment tests, redundancy analysis, and nonlinear methods. SNPs and microsatellites revealed two genetic lineages: Tibet and Hengduan Mountains-Western Sichuan Plateau (HDM-WSP), with reduced genetic diversity in Tibet lineage. More outlier loci were detected in HDM-WSP lineage than Tibet lineage. Among these, three SNPs in two genes responded to dry season precipitation in the HDM-WSP lineage but not in Tibet. By contrast, genetic variation in the Tibet lineage was related to geographic distance instead of the environment. Furthermore, risk of nonadaptedness (RONA) analyses suggested HDM-WSP lineage will have a better capacity to adapt in the predicted future climate compared with the Tibet lineage. We detected genetic imprints consistent with natural selection and molecular adaptation to drought on the Qinghai-Tibet Plateau (QTP) over a range of long-lived and widely distributed oak species in a changing environment. Our results suggest that different within-species adaptation processes occur in species occurring in heterogeneous environments.

Genomic insights into adaptation to heterogeneous environments for the ancient relictual Circaeaster agrestis (Circaeasteraceae, Ranunculales)

Investigating the interaction between environmental heterogeneity and local adaptation is critical for understanding the evolutionary history of a species, providing the premise for studying the response of organisms to rapid climate change. However, for most species how exactly the spatial heterogeneity promotes population divergence and how genomic variations contribute to adaptive evolution remain poorly understood. We examine the contributions of geographical and environmental variables to population divergence of the relictual, alpine herb Circaeaster agrestis, as well as the genetic basis of local adaptation using RAD-seq and plastome data. We detected significant genetic structure with an extraordinary disequilibrium of genetic diversity among regions, and signals of isolation-by-distance along with isolation-by-resistance. The populations were estimated to begin diverging in the late Miocene, along with a possible ancestral distribution of the Hengduan Mountains and adjacent regions. Both environmental gradient and redundancy analyses revealed significant association between genetic variation and temperature variables. Genome-environment association analyses identified 16 putatively adaptive loci related mainly to biotic and abiotic stress resistance. Our genome-wide data provide new insights into the important role of environmental heterogeneity in shaping genetic structure, and access the footprints of local adaptation in an ancient relictual species, informing future conservation efforts.

Ten Years of Provenance Trials and Application of Multivariate Random Forests Predicted the Most Preferable Seed Source for Silviculture of Abies sachalinensis in Hokkaido, Japan

Research highlights: Using 10-year tree height data obtained after planting from the range-wide provenance trials of Abies sachalinensis, we constructed multivariate random forests (MRF), a machine learning algorithm, with climatic variables. The constructed MRF enabled prediction of the optimum seed source to achieve good performance in terms of height growth at every planting site on a fine scale. Background and objectives: Because forest tree species are adapted to the local environment, local seeds are empirically considered as the best sources for planting. However, in some cases, local seed sources show lower performance in height growth than that showed by non-local seed sources. Tree improvement programs aim to identify seed sources for obtaining high-quality timber products by performing provenance trials. Materials and methods: Range-wide provenance trials for one of the most important silvicultural species, Abies sachalinensis, were established in 1980 at nine transplanting experimental sites. We constructed an MRF to estimate the responses of tree height at 10 years after planting at eight climatic variables at 1 km × 1 km resolution. The model was applied for prediction of tree height throughout Hokkaido Island. Results: Our model showed that four environmental variables were major factors affecting height growth—winter solar radiation, warmth index, maximum snow depth, and spring solar radiation. A tree height prediction map revealed that local seeds showed the best performance except in the southernmost region and several parts of northern regions. Moreover, the map of optimum seed provenance suggested that deployment of distant seed sources can outperform local sources in the southernmost and northern regions. Conclusions: We predicted that local seeds showed optimum growth, whereas non-local seeds had the potential to outperform local seeds in some regions. Several deployment options were proposed to improve tree growth.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Genomic and phenotypic divergence unveil microgeographic adaptation in the Amazonian hyperdominant tree Eperua falcata Aubl. (Fabaceae)

Plant populations can undergo very localized adaptation, allowing widely distributed populations to adapt to divergent habitats in spite of recurrent gene flow. Neotropical trees-whose large and undisturbed populations often span a variety of environmental conditions and local habitats-are particularly good models to study this process. Here, we explore patterns of adaptive divergence from large (i.e., regional) to small (i.e., microgeographic) spatial scales in the hyperdominant Amazonian tree Eperua falcata Aubl. (Fabaceae) under a replicated design involving two microhabitats (~300 m apart) in two study sites (~300 km apart). A three-year reciprocal transplant illustrates that, beyond strong maternal effects and phenotypic plasticity, genetically driven divergence in seedling growth and leaf traits was detected both between seedlings originating from different regions, and between seedlings from different microhabitats. In parallel, a complementary genome scan for selection was carried out through whole-genome sequencing of tree population pools. A set of 290 divergence outlier SNPs was detected at the regional scale (between study sites), while 185 SNPs located in the vicinity of 106 protein-coding genes were detected as replicated outliers between microhabitats within regions. Outlier-surrounding genomic regions are involved in a variety of physiological processes, including plant responses to stress (e.g., oxidative stress, hypoxia and metal toxicity) and biotic interactions. Together with evidence of microgeographic divergence in functional traits, the discovery of genomic candidates for microgeographic adaptive divergence represents a promising advance in our understanding of local adaptation, which probably operates across multiple spatial scales and underpins divergence and diversification in Neotropical trees.

Variation in the phenology of photosynthesis among eastern white pine provenances in response to warming

In higher-latitude trees, temperature and photoperiod control the beginning and end of the photosynthetically active season. Elevated temperature (ET) has advanced spring warming and delayed autumn cooling while photoperiod remains unchanged. We assessed the effects of warming on the length of the photosynthetically active season of three provenances of Pinus strobus L. seedlings from different latitudes, and evaluated the accuracy of the photochemical reflectance index (PRI) and the chlorophyll/carotenoid index (CCI) for tracking the predicted variation in spring and autumn phenology of photosynthesis among provenances. Seedlings from northern, local and southern P. strobus provenances were planted in a temperature-free-air-controlled enhancement (T-FACE) experiment and exposed to ET (+1.5/3°C; day/night). Over 18 months, we assessed photosynthetic phenology by measuring chlorophyll fluorescence, gas exchange, leaf spectral reflectance and pigment content. During autumn, all seedlings regardless of provenance followed the same sequence of phenological events with the initial downregulation of photosynthesis, followed by the modulation of non-photochemical quenching and associated adjustments of zeaxanthin pool sizes. However, the timing of autumn downregulation differed between provenances, with delayed onset in the southern provenance (SP) and earlier onset in the northern relative to the local provenance, indicating that photoperiod at the provenance origin is a dominant factor controlling autumn phenology. Experimental warming further delayed the downregulation of photosynthesis during autumn in the SP. A provenance effect during spring was also observed but was generally not significant. The vegetation indices PRI and CCI were both effective at tracking the seasonal variations of energy partitioning in needles and the differences of carotenoid pigments indicative of the stress status of needles. These results demonstrate that PRI and CCI can be useful tools for monitoring conifer phenology and for the remote monitoring of the length of the photosynthetically active season of conifers in a changing climate.

Plant adaptation to climate change - Where are we?

Climate change poses critical challenges for population persistence in natural communities, agriculture and environmental sustainability, and food security. In this review, we discuss recent progress in climatic adaptation in plants. We evaluate whether climate change exerts novel selection and disrupts local adaptation, whether gene flow can facilitate adaptive responses to climate change, and if adaptive phenotypic plasticity could sustain populations in the short term. Furthermore, we discuss how climate change influences species interactions. Through a more in-depth understanding of these eco-evolutionary dynamics, we will increase our capacity to predict the adaptive potential of plants under climate change. In addition, we review studies that dissect the genetic basis of plant adaptation to climate change. Finally, we highlight key research gaps, ranging from validating gene function, to elucidating molecular mechanisms, expanding research systems from model species to other natural species, testing the fitness consequences of alleles in natural environments, and designing multifactorial studies that more closely reflect the complex and interactive effects of multiple climate change factors. By leveraging interdisciplinary tools (e.g., cutting-edge omics toolkits, novel ecological strategies, newly-developed genome editing technology), researchers can more accurately predict the probability that species can persist through this rapid and intense period of environmental change, as well as cultivate crops to withstand climate change, and conserve biodiversity in natural systems.

Association of genetic and climatic variability in giant sequoia, Sequoiadendron giganteum, reveals signatures of local adaptation along moisture-related gradients

Uncovering the genetic basis of local adaptation is a major goal of evolutionary biology and conservation science alike. In an era of climate change, an understanding of how environmental factors shape adaptive diversity is crucial to predicting species response and directing management. Here, we investigate patterns of genomic variation in giant sequoia, an iconic and ecologically important tree species, using 1,364 bi-allelic single nucleotide polymorphisms (SNPs). We use an F _ST outlier test and two genotype-environment association methods, latent factor mixed models (LFMMs) and redundancy analysis (RDA), to detect complex signatures of local adaptation. Results indicate 79 genomic regions of potential adaptive importance, with limited overlap between the detection methods. Of the 58 loci detected by LFMM, 51 showed strong correlations to a precipitation-driven composite variable and seven to a temperature-related variable. RDA revealed 24 outlier loci with association to climate variables, all of which showed strongest relationship to summer precipitation. Nine candidate loci were indicated by two methods. After correcting for geographic distance, RDA models using climate predictors accounted for 49% of the explained variance and showed significant correlations between SNPs and climatic factors. Here, we present evidence of local adaptation in giant sequoia along gradients of precipitation and provide a first step toward identifying genomic regions of adaptive significance. The results of this study will provide information to guide management strategies that seek to maximize adaptive potential in the face of climate change.

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype-by-environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.

Genomics of Clinal Local Adaptation in Pinus sylvestris Under Continuous Environmental and Spatial Genetic Setting

Understanding the consequences of local adaptation at the genomic diversity is a central goal in evolutionary genetics of natural populations. In species with large continuous geographical distributions the phenotypic signal of local adaptation is frequently clear, but the genetic basis often remains elusive. We examined the patterns of genetic diversity in Pinus sylvestris, a keystone species in many Eurasian ecosystems with a huge distribution range and decades of forestry research showing that it is locally adapted to the vast range of environmental conditions. Making P. sylvestris an even more attractive subject of local adaptation study, population structure has been shown to be weak previously and in this study. However, little is known about the molecular genetic basis of adaptation, as the massive size of gymnosperm genomes has prevented large scale genomic surveys. We generated a both geographically and genomically extensive dataset using a targeted sequencing approach. By applying divergence-based and landscape genomics methods we identified several loci contributing to local adaptation, but only few with large allele frequency changes across latitude. We also discovered a very large (ca. 300 Mbp) putative inversion potentially under selection, which to our knowledge is the first such discovery in conifers. Our results call for more detailed analysis of structural variation in relation to genomic basis of local adaptation, emphasize the lack of large effect loci contributing to local adaptation in the coding regions and thus point out the need for more attention toward multi-locus analysis of polygenic adaptation.

Intraspecific Variation in Drought Response of Three Populations of Cryptocarya alba and Persea lingue, Two Native Species From Mediterranean Central Chile

An increase in the severity of drought events on Mediterranean climates highlights the need of using plant material adapted to drought during restoration efforts. Thus, we investigated between-population morpho-physiological differences in Cryptocarya alba and Persea lingue, two native species from Mediterranean central Chile, for traits that could effectively discriminate population performance in response to water restriction (WR) testing. Three populations from each species were subjected to WR treatment and physiological, morphological, and growth parameters were assessed at the beginning and at the end of the experiment. In C. alba, the most xeric population displayed smaller plants with mesophyllous leaves and lower photosynthetic rates indicating a resource saving strategy. Moreover, the xeric population performed better during WR than the most mesic populations, exhibiting higher water use efficiency (iWUE) and maintenance of growth rates. All C. alba populations responded equally to WR in terms of morphology and biomass partitioning. In contrast, differences among P. lingue populations were subtle at the morpho-physiological level with no apparent relation to provenance environmental conditions, and no morphological traits were affected by WR. However, in response to WR application, the most mesic population was, as observed through reduction in relative growth rates, more affected than xeric populations. We attribute such discrete differences between P. lingue provenances to the lower distributional range of selected populations. Our results show that relative growth rates in both species, and iWUE only in C. alba, exhibited population specific responses upon WR imposition; these results correspond with the environmental conditions found at the origin of each populations. Both traits could further assist in the selection of populations for restoration according to their response to water stress.

Assisted migration across fixed seed zones detects adaptation lags in two major North American tree species

Boreal forests are experiencing dramatic climate change, having warmed 1.0°-1.9°C over the last century. Yet forest regeneration practices are often still dictated by a fixed seed zone framework, in which seeds are both harvested from and planted into predefined areas. Our goal was to determine whether seedlings sourced from southern seed zones in Minnesota USA are already better adapted to northerly seed zones because of climate change. Bur oak (Quercus macrocarpa) and northern red oak (Quercus rubra) seedlings from two seed zones (i.e., tree ecotypes) were planted into 16 sites in two northern seed zones and measured for 3 yr. Our hypotheses were threefold: (1) tree species with more southern geographic distributions would thrive in northern forests where climate has already warmed substantially, (2) southern ecotypes of these species would have higher survival and growth than the northern ecotype in northern environments, and (3) natural selection would favor seedlings that expressed phenotypic and phenological traits characteristic of trees sourced from the more southern seed zone. For both species, survival was high (>93%), and southern ecotypes expressed traits consistent with our climate adaptation hypotheses. Ecotypic differences were especially evident for red oak; the southern ecotype had had higher survival, lower specific leaf area (SLA), faster height and diameter growth, and extended leaf phenology relative to the northern ecotype. Bur oak results were weaker, but the southern ecotype also had earlier budburst and lower SLA than the northern ecotype. Models based on the fixed seed zones failed to explain seedling performance as well as those with continuous predictors (e.g., climate and geographical position), suggesting that plant adaptations within current seed zone delineations do align with changing climate conditions. Adding support for this conclusion, natural selection favored traits expressed by the more southern tree ecotypes. Collectively, these results suggest that state seed sourcing guidelines should be reexamined to permit plantings across seed zones, a form of assisted migration. More extensive experiments (i.e., provenance trails) are necessary to make species-specific seed transfer guidelines that account for climate trends while also considering the precise geographic origin of seed sources.

Microgeographic adaptation and the effect of pollen flow on the adaptive potential of a temperate tree species

In species with long-distance dispersal capacities and inhabiting a large ecological niche, local selection and gene flow are expected to be major evolutionary forces affecting the genetic adaptation of natural populations. Yet, in species such as trees, evidence of microgeographic adaptation and the quantitative assessment of the impact of gene flow on adaptive genetic variation are still limited. Here, we used extensive genetic and phenotypic data from European beech seedlings collected along an elevation gradient, and grown in a common garden, to study the signature of selection on the divergence of eleven potentially adaptive traits, and to assess the role of gene flow in resupplying adaptive genetic variation. We found a significant signal of adaptive differentiation among plots separated by < 1 km, with selection acting on growth and phenological traits. Consistent with theoretical expectations, our results suggest that pollen dispersal contributes to increase genetic diversity for these locally differentiated traits. Our results thus highlight that local selection is an important evolutionary force in natural tree populations and suggest that management interventions to facilitate movement of gametes along short ecological gradients would boost genetic diversity of individual tree populations, and enhance their adaptive potential to rapidly changing environments.

Candidate gene SNP variation in floodplain populations of pedunculate oak (Quercus robur L.) near the species' southern range margin: Weak differentiation yet distinct associations with water availability

Populations residing near species' low-latitude range margins (LLMs) often occur in warmer and drier environments than those in the core range. Thus, their genetic composition could be shaped by climatic drivers that differ from those occurring at higher latitudes, resulting in potentially adaptive variants of conservation value. Such variants could facilitate the adaptation of populations from other portions of the geographical range to similar future conditions anticipated under ongoing climate change. However, very few studies have assessed standing genetic variation at potentially adaptive loci in natural LLM populations. We investigated standing genetic variation at single nucleotide polymorphisms (SNPs) located within 117 candidate genes and its links to putative climatic selection pressures across 19 pedunculate oak (Quercus robur L.) populations distributed along a regional climatic gradient near the species' southern range margin in southeastern Europe. These populations are restricted to floodplain forests along large lowland rivers, whose hydric regime is undergoing significant shifts under modern rapid climate change. The populations showed very weak geographical structure, suggesting extensive genetic connectivity and gene flow or shared ancestry. We identified eight (6.2%) positive F_ST -outlier loci, and genotype-environment association analyses revealed consistent associations between SNP allele frequencies and several climatic variables linked to water availability. A total of 61 associations involving 37 SNPs (28.5%) from 35 annotated genes provided important insights into putative functional mechanisms in our system. Our findings provide empirical support for the role of LLM populations as sources of potentially adaptive variation that could enhance species' resilience to climate change-related pressures.

Synchronization of senescence and desynchronization of flowering in Arabidopsis thaliana

In a recent publication, we proposed that adjusting lifespan in order to synchronize senescence is important for timing of reproduction, and we quantified the synchrony of reproductive timing relative to germination timing. Here, in a second sequential seeding experiment (SSE), the germination timing of Arabidopsis thaliana accessions was manipulated and plants were then grown under two different temperature regimes. Life stage traits of plants in each temperature regime were analysed and it was evaluated whether the cohorts were grouped according to age and/or environmental conditions. While flowering-related traits showed desynchrony among cohorts, striking synchrony in the timing of senescence among cohorts for each group was found. A quantitative trait locus (QTL) analysis using a genotyped population of 'Cvi/Ler' recombinant inbred lines (RILs) was then conducted. Novel and known loci were assigned to flowering and senescence timing. However, senescence synchrony resulted in low variation in senescence time and weak QTL detection for flowering termination. Overlapping flowering and senescence genes with loci affecting either of those traits were found and suggest a potential interdependency of reproductive traits.

Genetic variation in early fitness traits across European populations of silver birch (Betula pendula)

Given that the ecological niche of tree species is typically narrower for earlier life stages, intraspecific genetic variation at early fitness traits may greatly influence the adaptive response of tree populations to changing environmental conditions. In this study, we evaluated genetic variation in early fitness traits among 12 populations of Betula pendula from a wide latitudinal range in Europe (41-55°N). We first conducted a chamber experiment to test for population differences in germination and the effect of pre-chilling treatment on seed dormancy release. We then established three common gardens spread across the species latitudinal range in order to evaluate levels of quantitative genetic variation and genotype-by-environment interaction at different early life traits. Our results showed significant variation in chamber germination rates among populations (0-60 %), with southern populations exhibiting lower germination. Pre-chilling treatments did not generally improve germination success. Population seedling emergence rates in the field were correlated with chamber germination rates, though being an order of magnitude lower, with an average ranging from 0 to 1.3 % across gardens. Highly significant variation was found in field emergence rates among populations, and between seed-crop years within populations, but not among families within populations. Populations differed in seedling height, diameter, slenderness and budburst date, with significant among-family variation. Population latitude was positively associated with chamber germination rate and with seedling emergence rate in one of the central field sites. Overall, genetic, environmental and demographic factors seem to influence the observed high levels of variation in early fitness traits among B. pendula populations. Our results suggest limited regeneration capacity for the study species under drier conditions, but further field trials with sufficient replication over environments and seed crops will improve our understanding of its vulnerability to climate change.

Meta-analysis of the differential effects of habitat fragmentation and degradation on plant genetic diversity

Genetic diversity is a key factor for population survival and evolution. However, anthropogenic habitat disturbance can erode it, making populations more prone to extinction. Aiming to assess the global effects of habitat disturbance on plant genetic variation, we conducted a meta-analysis based on 92 case studies obtained from published literature. We compared the effects of habitat fragmentation and degradation on plant allelic richness and gene diversity (equivalent to expected heterozygosity) and tested whether such changes are sensitive to different life-forms, life spans, mating systems, and commonness. Anthropogenic disturbance had a negative effect on allelic richness, but not on gene diversity. Habitat fragmentation had a negative effect on genetic variation, whereas habitat degradation had no effect. When we examined the individual effects in fragmented habitats, allelic richness and gene diversity decreased, but this decrease was strongly dependent on certain plant traits. Specifically, common long-lived trees and self-incompatible species were more susceptible to allelic richness loss. Conversely, gene diversity decreased in common short-lived species (herbs) with self-compatible reproduction. In a wider geographical context, tropical plant communities were more sensitive to allelic richness loss, whereas temperate plant communities were more sensitive to gene diversity loss. Our synthesis showed complex responses to habitat disturbance among plant species. In many cases, the absence of effects could be the result of the time elapsed since the disturbance event or reproductive systems favoring self-pollination, but attention must be paid to those plant species that are more susceptible to losing genetic diversity, and appropriate conservation should be actions taken.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Intraspecific trait variation across elevation predicts a widespread tree species' climate niche and range limits

Global change is widely altering environmental conditions which makes accurately predicting species range limits across natural landscapes critical for conservation and management decisions. If climate pressures along elevation gradients influence the distribution of phenotypic and genetic variation of plant functional traits, then such trait variation may be informative of the selective mechanisms and adaptations that help define climatic niche limits. Using extensive field surveys along 16 elevation transects and a large common garden experiment, we tested whether functional trait variation could predict the climatic niche of a widespread tree species (Populus angustifolia) with a double quantile regression approach. We show that intraspecific variation in plant size, growth, and leaf morphology corresponds with the species' total climate range and certain climatic limits related to temperature and moisture extremes. Moreover, we find evidence of genetic clines and phenotypic plasticity at environmental boundaries, which we use to create geographic predictions of trait variation and maximum values due to climatic constraints across the western US. Overall, our findings show the utility of double quantile regressions for connecting species distributions and climate gradients through trait-based mechanisms. We highlight how new approaches like ours that incorporate genetic variation in functional traits and their response to climate gradients will lead to a better understanding of plant distributions as well as identifying populations anticipated to be maladapted to future environments.

Adaptive introgression as a driver of local adaptation to climate in European white oaks

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

How well do the spring indices predict phenological activity across plant species?

The spring indices, models that represent the onset of spring season biological activity, were developed using a long-term observational record from the mid-to-late twentieth century of three species of lilacs and honeysuckles contributed by volunteer observers across the nation. The USA National Phenology Network (USA-NPN) produces and freely delivers maps of spring index onset dates at fine spatial scale for the USA. These maps are used widely in natural resource planning and management applications. The extent to which the models represent activity in a broad suite of plant species is not well documented. In this study, we used a rich record of observational plant phenology data (37,819 onset records) collected in recent years (1981-2017) to evaluate how well gridded maps of the spring index models predict leaf and flowering onset dates in (a) 19 species of ecologically important, broadly distributed deciduous trees and shrubs, and (b) the lilac and honeysuckle species used to construct the models. The extent to which the spring indices predicted vegetative and reproductive phenology varied by species and with latitude, with stronger relationships revealed for shrubs than trees and with the Bloom Index compared to the Leaf Index, and reduced concordance between the indices at higher latitudes. These results allow us to use the indices as indicators of when to expect activity across widely distributed species and can serve as a yardstick to assess how future changes in the timing of spring will impact a broad array of trees and shrubs across the USA.

Analyses of Genetic Diversity, Differentiation and Geographic Origin of Natural Provenances and Land Races of Casuarina equisetifolia Based on EST-SSR Markers

Research Highlights: High variation of genetic diversity and differentiation among 27 seed sources within 14 natural provenances and 13 land race samples of Casuarina equisetifolia were found. High proportions of monoecious individuals may be present in some populations, as indicated by severe heterozytote deficiency and inbreeding found in many provenances and land races. The most probable origins of the land races were inferred according to the values of pairwise provenance differentiation and Nei’s genetic distances. Targeted introductions and testing of unrelated new accessions of C. equisetifolia from the Pacific and Philippines was proposed to identify Ralstonia-resistant genotypes. Background and Objectives: Casuarina equisetifolia was introduced to China a hundred years ago and has become a critically important tree species in coastal protection since the 1950s. Despite its importance, patterns of genetic variation, genetic relationships among natural provenances and probable origins of the land races remain unresolved. This has become a concern in China where Ralstonia solanacearum bacterial wilt has devastated plantations that are known to be from a narrow genetic base that urgently needs to be broadened. Materials and Methods: Fourteen natural provenances from Australia, Pacific islands and Southeast Asia, and 13 land race samples from parts of Asia and Africa outside the natural range were genotyped using 13 SSR (Simple Sequence Repeats) markers to characterize their allelic variation and genetic relationship. Results: Significant genetic diversity and differentiation among 27 seed sources within 14 provenances and 13 land race samples of C. equisetifolia was indicated. Significant heterozygote deficiency and inbreeding was indicated for a number of provenances, perhaps indicating a high proportion of monoecious parents in these populations. The most probable origins of the land races of the introduced countries were suggested according to the values of pairwise provenance differentiation (FST) and Nei’s genetic distances. Conclusions: We found significant genetic diversity and genetic differentiation among seed sources of C. equisetifolia. While individual land races do not appear to lack diversity, we were able to infer the origins of some, allowing targeted introductions of unrelated material to be made. In the case of the Chinese land race, targeting and testing new accessions from the Pacific and the Philippines may be a good strategy to identify Ralstonia-resistant genotypes.

"Green pointillism": detecting the within-population variability of budburst in temperate deciduous trees with phenological cameras

Phenological cameras have been used over a decade for identifying plant phenological markers (budburst, leaf senescence) and more generally the greenness dynamics of forest canopies. The analysis is usually carried out over the full camera field of view, with no particular analysis of the variability of phenological markers among trees. Here we show that images produced by phenological cameras can be used to quantify the within-population variability of budburst (WPVbb) in temperate deciduous forests. Using seven site-years of image analyses, we report a strong correlation (r² = 0.97) between the WPVbb determined with a phenological camera and its quantification through ground observation. We show that WPVbb varies strongly (by a factor of 4) from year to year in a given population and that those variations are linked with temperature conditions during the budburst period, with colder springs associated to a higher differentiation of budburst (higher WPVbb) among trees. Deploying our approach at the continental scale, i.e., throughout phenological cameras networks, would improve the understanding of the spatial (across populations) and temporal (across years) variations of WPVbb, which have strong implications on forest functioning, tree fitness and phenological modelling.

Environmental Genome-Wide Association Reveals Climate Adaptation Is Shaped by Subtle to Moderate Allele Frequency Shifts in Loblolly Pine

Understanding the genomic basis of local adaptation is crucial to determine the potential of long-lived woody species to withstand changes in their natural environment. In the past, efforts to dissect the genomic architecture in gymnosperms species have been limited due to the absence of reference genomes. Recently, the genomes of some commercially important conifers, such as loblolly pine, have become available, allowing whole-genome studies of these species. In this study, we test for associations between 87k SNPs, obtained from whole-genome resequencing of loblolly pine individuals, and 270 environmental variables and combinations of them. We determine the geographic location of significant loci and identify their genomic location using our newly constructed ultradense 26k SNP linkage map. We found that water availability is the main climatic variable shaping local adaptation of the species, and found 821 SNPs showing significant associations with climatic variables or combinations of them based on the consistent results of three different genotype–environment association methods. Our results suggest that adaptation to climate in the species might have occurred by many changes in the frequency of alleles with moderate to small effect sizes, and by the smaller contribution of large effect alleles in genes related to moisture deficit, temperature and precipitation. Genomic regions of low recombination and high population differentiation harbored SNPs associated with groups of environmental variables, suggesting climate adaptation might have evolved as a result of different selection pressures acting on groups of genes associated with an aspect of climate rather than on individual environmental variables.

A New Perspective on Ecological Prediction Reveals Limits to Climate Adaptation in a Temperate Tree Species

Forests absorb a large fraction of anthropogenic CO2 emission, but their ability to continue to act as a sink under climate change depends in part on plant species undergoing rapid adaptation. Yet models of forest response to climate change currently ignore local adaptation as a response mechanism. Thus, considering the evolution of intraspecific trait variation is necessary for reliable, long-term species and climate projections. Here, we combine ecophysiology and predictive climate modeling with analyses of genomic variation to determine whether sugar and starch storage, energy reserves for trees under extreme conditions, have the heritable variation and genetic diversity necessary to evolve in response to climate change within populations of black cottonwood (Populus trichocarpa). Despite current patterns of local adaptation and extensive range-wide heritable variation in storage, we demonstrate that adaptive evolution in response to climate change will be limited by a lack of heritable variation within northern populations and by a need for extreme genetic changes in southern populations. Our method can help design more targeted species management interventions and highlights the power of using genomic tools in ecological prediction to scale from molecular to regional processes to determine the ability of a species to respond to future climates.

Population Genomic Analysis of Two Endemic Schizothoracins Reveals Their Genetic Differences and Underlying Selection Associated with Altitude and Temperature

Simple Summary

Diptychus maculates and Gymnodiptychus dybowskii are two rare aboriginal fishes in the Xinjiang Uygur Autonomous Region. In recent years, due to overfishing and habitat fragmentation caused by construction of water conservancy and hydropower projects, the fishery resources have decreased sharply. Understanding the genetic background is of great significance for resource protection. In this study, we revealed the similar trends of population genetic diversities in these two species collected from the Tarim River and the Yili River. In addition, outlier SNPs associated with temperature and altitude were detected in both of them, indicating that Schizothoracinae fishes represented by D. maculates and G. dybowskii were still under the selection pressure of plateau environments.

Abstract

Schizothoracins are a group of cyprinid fishes distributed throughout the Qinghai–Tibet Plateau, which can be classified in three grades: primitive, specialised and highly specialised according to adaptation ability to plateau environments. As the only specialised schizothoracins in Xinjiang, China, Diptychus maculates and Gymnodiptychus dybowskii are ideal materials for adaptive evolution research. Based on single-nucleotide polymorphism (SNP) loci detected by specific-locus amplified fragment (SLAF) technology, the genome-wide genetic diversities of these two species from nine sites in Xinjiang were evaluated. D.maculates in the Muzat River (BM) and G. dybowskii in the Kaidu River (LKG) presented the lowest genetic diversity levels, whereas D. maculates in the Kumarik River (BK) and G.dybowskii in the Kashi River (LK) were just the opposite. Cluster and principal component analysis demonstrated a distant genetic affinity between D. maculates in the Tashkurgan River (BT) and other populations. Outlier SNP loci were discovered both in D. maculates and G. dybowskii. The coalescent Bayenv and latent factor mixed model (LFMM) methods showed that a total of thirteen and eighteen SNPs in D. maculates were associated with altitude and temperature gradient, respectively. No intersection was revealed in G. dybowskii. The results indicated that D. maculates was subject to much greater divergent selection pressure. A strong signal of isolation-by-distance (IBD) was detected across D. maculates (Mantel test, rs = 0.65; p = 0.05), indicating an evident geographical isolation in the Tarim River. Isolation-by-environment (IBE) analysis implied that temperature and altitude selections were more intensive in D. maculates, with greater environmental variation resulting in weak gene flow.

Pleiotropy facilitates local adaptation to distant optima in common ragweed (Ambrosia artemisiifolia)

Pleiotropy, the control of multiple phenotypes by a single locus, is expected to slow the rate of adaptation by increasing the chance that beneficial alleles also have deleterious effects. However, a prediction arising from classical theory of quantitative trait evolution states that pleiotropic alleles may have a selective advantage when phenotypes are distant from their selective optima. We examine the role of pleiotropy in regulating adaptive differentiation among populations of common ragweed (Ambrosia artemisiifolia); a species that has recently expanded its North American range due to human-mediated habitat change. We employ a phenotype-free approach by using connectivity in gene networks as a proxy for pleiotropy. First, we identify loci bearing footprints of local adaptation, and then use genotype-expression mapping and co-expression networks to infer the connectivity of the genes. Our results indicate that the putatively adaptive loci are highly pleiotropic, as they are more likely than expected to affect the expression of other genes, and they reside in central positions within the gene networks. We propose that the conditionally advantageous alleles at these loci avoid the cost of pleiotropy by having large phenotypic effects that are beneficial when populations are far from their selective optima. We further use evolutionary simulations to show that these patterns are in agreement with a model where populations face novel selective pressures, as expected during a range expansion. Overall, our results suggest that highly connected genes may be targets of positive selection during environmental change, even though they likely experience strong purifying selection in stable selective environments.

Gene Expression Modularity Reveals Footprints of Polygenic Adaptation in Theobroma cacao

Separating footprints of adaptation from demography is challenging. When selection has acted on a single locus with major effect, this issue can be alleviated through signatures left by selective sweeps. However, as adaptation is often driven by small allele frequency shifts at many loci, studies focusing on single genes are able to identify only a small portion of genomic variants responsible for adaptation. In face of this challenge, we utilize coexpression information to search for signals of polygenetic adaptation in Theobroma cacao, a tropical tree species that is the source of chocolate. Using transcriptomics and a weighted correlation network analysis, we group genes with similar expression patterns into functional modules. We then ask whether modules enriched for specific biological processes exhibit cumulative effects of differential selection in the form of high FST and dXY between populations. Indeed, modules putatively involved in protein modification, flowering, and water transport show signs of polygenic adaptation even though individual genes that are members of those groups do not bear strong signatures of selection. Modeling of demography, background selection, and the effects of genomic features reveal that these patterns are unlikely to arise by chance. We also find that specific modules are enriched for signals of strong or relaxed purifying selection, with one module bearing signs of adaptive differentiation and an excess of deleterious mutations. Our results provide insight into polygenic adaptation and contribute to understanding of population structure, demographic history, and genome evolution in T. cacao.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

Range margin populations show high climate adaptation lags in European trees

How populations of long-living species respond to climate change depends on phenotypic plasticity and local adaptation processes. Marginal populations are expected to have lags in adaptation (i.e. differences between the climatic optimum that maximizes population fitness and the local climate) because they receive pre-adapted alleles from core populations preventing them from reaching a local optimum in their climatically marginal habitat. Yet, whether adaptation lags in marginal populations are a common feature across phylogenetically and ecologically different species and how lags can change with climate change remain unexplored. To test for range-wide patterns of phenotypic variation and adaptation lags of populations to climate, we (a) built model ensembles of tree height accounting for the climate of population origin and the climate of the site for 706 populations monitored in 97 common garden experiments covering the range of six European forest tree species; (b) estimated populations' adaptation lags as the differences between the climatic optimum that maximizes tree height and the climate of the origin of each population; (c) identified adaptation lag patterns for populations coming from the warm/dry and cold/wet margins and from the distribution core of each species range. We found that (a) phenotypic variation is driven by either temperature or precipitation; (b) adaptation lags are consistently higher in climatic margin populations (cold/warm, dry/wet) than in core populations; (c) predictions for future warmer climates suggest adaptation lags would decrease in cold margin populations, slightly increasing tree height, while adaptation lags would increase in core and warm margin populations, sharply decreasing tree height. Our results suggest that warm margin populations are the most vulnerable to climate change, but understanding how these populations can cope with future climates depend on whether other fitness-related traits could show similar adaptation lag patterns.

Phenotypic integration and life history strategies among populations of Pinus halepensis: an insight through structural equation modelling

BACKGROUND AND AIMS

Understanding inter-population variation in the allocation of resources to specific anatomical compartments and physiological processes is crucial to disentangle adaptive patterns in forest species. This work aims to evaluate phenotypic integration and trade-offs among functional traits as determinants of life history strategies in populations of a circum-Mediterranean pine that dwells in environments where water and other resources are in limited supply.

METHODS

Adult individuals of 51 populations of Pinus halepensis grown in a common garden were characterized for 11 phenotypic traits, including direct and indirect measures of water uptake at different depths, leaf area, stomatal conductance, chlorophyll content, non-structural carbohydrates, stem diameter and tree height, age at first reproduction and cone production. The population differentiation in these traits was tested through analysis of variance (ANOVA). The resulting populations' means were carried forward to a structural equation model evaluating phenotypic integration between six latent variables (summer water uptake depth, summer transpiration, spring photosynthetic capacity, growth, reserve accumulation and reproduction).

KEY RESULTS

Water uptake depth and transpiration covaried negatively among populations, as the likely result of a common selective pressure for drought resistance, while spring photosynthetic capacity was lower in populations originating from dry areas. Transpiration positively influenced growth, while growth was negatively related to reproduction and reserves among populations. Water uptake depth negatively influenced reproduction.

CONCLUSIONS

The observed patterns indicate a differentiation in life cycle features between fast-growing and slow-growing populations, with the latter investing significantly more in reproduction and reserves. We speculate that such contrasting strategies result from different arrays of life history traits underlying the very different ecological conditions that the Aleppo pine must face across its distribution range. These comprise, principally, drought as the main stressor and fire as the main ecological disturbance of the Mediterranean basin.

A multiscale approach to detect selection in nonmodel tree species: Widespread adaptation despite population decline in Taxus baccata L

Detecting the molecular basis of local adaptation and identifying selective drivers is still challenging in nonmodel species. The use of purely population genetic approaches is limited by some characteristics of genetic systems, such as pleiotropy and polygenic control, and parallel evidence from phenotypic-based experimental comparisons is required. In long-lived organisms, the detection of selective pressures might also be precluded by evolutionary lag times in response to the environment. Here, we used the English yew to showcase an example of a multiscale integrative approach in a nonmodel species with limited plant and genomic resources. We combined information from two independent sources, phenotypes in a common environment and genomic data in natural populations, to investigate the signature of selection. Growth differences among populations in a common environment, and phenological patterns of both shoot elongation and male strobili maturation, were associated with climate clines, providing evidence for local adaptation and guiding us in the selection of populations for genomic analyses. We used information on over 25,000 SNPs from c. 1,200 genes to infer the demographic history and to test for molecular signatures of selection at different levels: SNP, gene, and biological pathway. Our results confirmed an overall demographic history of population decline, but we also found evidence for putative local adaptation at the molecular level. We identified or confirmed several candidate genes for positive and negative selection in forest trees, including the pseudo-response regulator 7 (PRR7), an essential component of the circadian clock in plants. In addition, we successfully tested an approach to detect polygenic adaptation in biological pathways, allowing us to identify the flavonoid biosynthesis pathway as a candidate stress-response pathway that deserves further attention in other plants. Finally, our study contributes to the emerging view that explaining contemporary standing genetic variation requires considering adaptation to past climates, especially for long-lived trees.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Forest genomics: Advancing climate adaptation, forest health, productivity, and conservation

Forest ecosystems provide important ecological services and resources, from habitat for biodiversity to the production of environmentally friendly products, and play a key role in the global carbon cycle. Humanity is counting on forests to sequester and store a substantial portion of the anthropogenic carbon dioxide produced globally. However, the unprecedented rate of climate change, deforestation, and accidental importation of invasive insects and diseases are threatening the health and productivity of forests, and their capacity to provide these services. Knowledge of genetic diversity, local adaptation, and genetic control of key traits is required to predict the adaptive capacity of tree populations, inform forest management and conservation decisions, and improve breeding for productive trees that will withstand the challenges of the 21st century. Genomic approaches have well accelerated the generation of knowledge of the genetic and evolutionary underpinnings of nonmodel tree species, and advanced their applications to address these challenges. This special issue of Evolutionary Applications features 14 papers that demonstrate the value of a wide range of genomic approaches that can be used to better understand the biology of forest trees, including species that are widespread and managed for timber production, and others that are threatened or endangered, or serve important ecological roles. We highlight some of the major advances, ranging from understanding the evolution of genomes since the period when gymnosperms separated from angiosperms 300 million years ago to using genomic selection to accelerate breeding for tree health and productivity. We also discuss some of the challenges and future directions for applying genomic tools to address long-standing questions about forest trees.

Evaluating genomic data for management of local adaptation in a changing climate: A lodgepole pine case study

We evaluate genomic data, relative to phenotypic and climatic data, as a basis for assisted gene flow and genetic conservation. Using a seedling common garden trial of 281 lodgepole pine (Pinus contorta) populations from across western Canada, we compare genomic data to phenotypic and climatic data to assess their effectiveness in characterizing the climatic drivers and spatial scale of local adaptation in this species. We find that phenotype-associated loci are equivalent or slightly superior to climate data for describing local adaptation in seedling traits, but that climate data are superior to genomic data that have not been selected for phenotypic associations. We also find agreement between the climate variables associated with genomic variation and with 20-year heights from a long-term provenance trial, suggesting that genomic data may be a viable option for identifying climatic drivers of local adaptation where phenotypic data are unavailable. Genetic clines associated with the experimental traits occur at broad spatial scales, suggesting that standing variation of adaptive alleles for this and similar species does not require management at scales finer than those indicated by phenotypic data. This study demonstrates that genomic data are most useful when paired with phenotypic data, but can also fill some of the traditional roles of phenotypic data in management of species for which phenotypic trials are not feasible.

Experimental migration upward in elevation is associated with strong selection on life history traits

One of the strongest biological impacts of climate change has been the movement of species poleward and upward in elevation. Yet, what is not clear is the extent to which the spatial distribution of locally adapted lineages and ecologically important traits may also shift with continued climate change. Here, we take advantage of a transplant experiment mimicking up-slope seed dispersal for a suite of ecologically diverse populations of yellow monkeyflower (Mimulus guttatus sensu lato) into a high-elevation common garden during an extreme drought period in the Sierra Nevada mountains, California, USA. We use a demographic approach to quantify fitness and test for selection on life history traits in local versus lower-elevation populations and in normal versus drought years to test the potential for up-slope migration and phenotypic selection to alter the distribution of key life history traits in montane environments. We find that lower-elevation populations tend to outperform local populations, confirming the potential for up-slope migration. Although selection generally favored some local montane traits, including larger flowers and larger stem size at flowering, drought conditions tended to select for earlier flowering typical of lower-elevation genotypes. Taken together, this suggests that monkeyflower lineages moving upward in elevation could experience selection for novel trait combinations, particularly under warmer and drier conditions that are predicted to occur with continued climate change.

Genotypic variation in Norway spruce correlates to fungal communities in vegetative buds

The taxonomically diverse phyllosphere fungi inhabit leaves of plants. Thus, apart from the fungi's dispersal capacities and environmental factors, the assembly of the phyllosphere community associated with a given host plant depends on factors encoded by the host's genome. The host genetic factors and their influence on the assembly of phyllosphere communities under natural conditions are poorly understood, especially in trees. Recent work indicates that Norway spruce (Picea abies) vegetative buds harbour active fungal communities, but these are hitherto largely uncharacterized. This study combines internal transcribed spacer sequencing of the fungal communities associated with dormant vegetative buds with a genome-wide association study (GWAS) in 478 unrelated Norway spruce trees. The aim was to detect host loci associated with variation in the fungal communities across the population, and to identify loci correlating with the presence of specific, latent, pathogens. The fungal communities were dominated by known Norway spruce phyllosphere endophytes and pathogens. We identified six quantitative trait loci (QTLs) associated with the relative abundance of the dominating taxa (i.e., top 1% most abundant taxa). Three additional QTLs associated with colonization by the spruce needle cast pathogen Lirula macrospora or the cherry spruce rust (Thekopsora areolata) in asymptomatic tissues were detected. The identification of the nine QTLs shows that the genetic variation in Norway spruce influences the fungal community in dormant buds and that mechanisms underlying the assembly of the communities and the colonization of latent pathogens in trees may be uncovered by combining molecular identification of fungi with GWAS.

Local adaptation of a dominant coastal tree to freshwater availability and solar radiation suggested by genomic and ecophysiological approaches

Local adaptation is often a product of environmental variations in geographical space and has implications for biodiversity conservation. We investigated the role of latitudinal heterogeneity in climate on the organization of genetic and phenotypic variation in the dominant coastal tree Avicennia schaueriana. In a common garden experiment, samples from an equatorial region, with pronounced seasonality in precipitation, accumulated less biomass, and showed lower stomatal conductance and transpiration, narrower xylem vessels, smaller leaves and higher reflectance of long wavelengths by the stem epidermis than samples from a subtropical region, with seasonality in temperature and no dry season. Transcriptomic differences identified between trees sampled under field conditions at equatorial and subtropical sites, were enriched in functional categories such as responses to temperature, solar radiation, water deficit, photosynthesis and cell wall biosynthesis. Remarkably, the diversity based on genome-wide SNPs revealed a north-south genetic structure and signatures of selection were identified for loci associated with photosynthesis, anthocyanin accumulation and the responses to osmotic and hypoxia stresses. Our results suggest the existence of divergence in key resource-use characteristics, likely driven by seasonality in water deficit and solar radiation. These findings provide a basis for conservation plans and for predicting coastal plants responses to climate change.

Towards a New Generation of Trait-Flexible Vegetation Models

Plant trait variability, emerging from eco-evolutionary dynamics that range from alleles to macroecological scales, is one of the most elusive, but possibly most consequential, aspects of biodiversity. Plasticity, epigenetics, and genetic diversity are major determinants of how plants will respond to climate change, yet these processes are rarely represented in current vegetation models. Here, we provide an overview of the challenges associated with understanding the causes and consequences of plant trait variability, and review current developments to include plasticity and evolutionary mechanisms in vegetation models. We also present a roadmap of research priorities to develop a next generation of vegetation models with flexible traits. Including trait variability in vegetation models is necessary to better represent biosphere responses to global change.

An analysis of mating biases in trees

Assortative mating is a deviation from random mating based on phenotypic similarity. As it is much better studied in animals than in plants, we investigate for trees whether kinship of realized mating pairs deviates from what is expected from the set of potential mates and use this information to infer mating biases that may result from kin recognition and/or assortative mating. Our analysis covers 20 species of trees for which microsatellite data is available for adult populations (potential mates) as well as seed arrays. We test whether mean relatedness of observed mating pairs deviates from null expectations that only take pollen dispersal distances into account (estimated from the same data set). This allows the identification of elevated as well as reduced kinship among realized mating pairs, indicative of positive and negative assortative mating, respectively. The test is also able to distinguish elevated biparental inbreeding that occurs solely as a result of related pairs growing closer to each other from further assortativeness. Assortative mating in trees appears potentially common but not ubiquitous: nine data sets show mating bias with elevated inbreeding, nine do not deviate significantly from the null expectation, and two show mating bias with reduced inbreeding. While our data sets lack direct information on phenology, our investigation of the phenological literature for each species identifies flowering phenology as a potential driver of positive assortative mating (leading to elevated inbreeding) in trees. Since active kin recognition provides an alternative hypothesis for these patterns, we encourage further investigations on the processes and traits that influence mating patterns in trees.

Variation in Susceptibility of Tanoak to the NA1 and EU1 Lineages of Phytophthora ramorum, the Cause of Sudden Oak Death

Phytophthora ramorum, the cause of sudden oak death (SOD), kills tanoak (Notholithocarpus densiflorus) trees in southwestern Oregon and California. Two lineages of P. ramorum are now found in wildland forests of Oregon (NA1 and EU1). In addition to the management of SOD in forest ecosystems, disease resistance could be used as a way to mitigate the impact of P. ramorum. The objectives of this study were to (i) characterize the variability in resistance of N. densiflorus among families using lesion length; (ii) determine whether lineage, isolate, family, or their interactions significantly affect variation in lesion length; and (iii) determine whether there are differences among isolates and among families in terms of lesion length. The parameters isolate nested within lineage (isolate[lineage]) and family × isolate(lineage) interaction explained the majority of the variation in lesion length. There was no significant difference between the NA1 and EU1 lineages in terms of mean lesion length; however, there were differences among the six isolates. Lesions on seedlings collected from surviving trees at infested sites were smaller, on average, than lesions of seedlings collected from trees at noninfested sites (P = 0.0064). The results indicate that there is potential to establish a breeding program for tanoak resistance to SOD and that several isolates of P. ramorum should be used in an artificial inoculation assay.

Downhill seed dispersal by temperate mammals: a potential threat to plant escape from global warming

Vertical seed dispersal, i.e. seed dispersal towards a higher or lower altitude, is considered a critical process for plant escape from climate change. However, studies exploring vertical seed dispersal are scarce, and thus, its direction, frequency, and mechanisms are little known. In the temperate zone, evaluating vertical seed dispersal of animal-dispersed plants fruiting in autumn and/or winter is essential considering the dominance of such plants in temperate forests. We hypothesized that their seeds are dispersed towards lower altitudes because of the downhill movement of frugivorous animals following the autumn-to-winter phenology of their food plants which proceeds from the mountain tops to the foot in the temperate zone. We evaluated the vertical seed dispersal of the autumn-fruiting wild kiwi, Actinidia arguta, which is dispersed by temperate mammals. We collected dispersed seeds from mammal faeces in the Kanto Mountains of central Japan and estimated the distance of vertical seed dispersal using the oxygen isotope ratios of the dispersed seeds. We found the intensive downhill seed dispersal of wild kiwi by all seed dispersers, except the raccoon dog (bear: mean −393.1 m; marten: −245.3 m; macaque: −98.5 m; and raccoon dog: +4.5 m). Mammals with larger home ranges dispersed seeds longer towards the foot of the mountains. Furthermore, we found that seeds produced at higher altitudes were dispersed a greater distance towards the foot of the mountains. Altitudinal gradients in autumn-to-winter plant phenology and other mountain characteristics, i.e. larger surface areas and more attractive human crops at lower altitudes compared to higher altitudes, were considered drivers of downhill seed dispersal via animal movement. Strong downhill seed dispersal by mammals suggests that populations of autumn-to-winter fruiting plants dispersed by animals may not be able to sufficiently escape from current global warming in the temperate zone.

Looking for Local Adaptation: Convergent Microevolution in Aleppo Pine (Pinus halepensis)

Finding outlier loci underlying local adaptation is challenging and is best approached by suitable sampling design and rigorous method selection. In this study, we aimed to detect outlier loci (single nucleotide polymorphisms, SNPs) at the local scale by using Aleppo pine (Pinus halepensis), a drought resistant conifer that has colonized many habitats in the Mediterranean Basin, as the model species. We used a nested sampling approach that considered replicated altitudinal gradients for three contrasting sites. We genotyped samples at 294 SNPs located in genomic regions selected to maximize outlier detection. We then applied three different statistical methodologies-Two Bayesian outlier methods and one latent factor principal component method-To identify outlier loci. No SNP was an outlier for all three methods, while eight SNPs were detected by at least two methods and 17 were detected only by one method. From the intersection of outlier SNPs, only one presented an allelic frequency pattern associated with the elevational gradient across the three sites. In a context of multiple populations under similar selective pressures, our results underline the need for careful examination of outliers detected in genomic scans before considering them as candidates for convergent adaptation.

Ecotypic differentiation reveals seed colour-related alkaloid content in a crop wild relative

Crop wild relatives can be a useful source of genotypes that maximise crop survival and yield in specific habitats. Lupinus angustifolius is an annual forb with crop varieties derived from a narrow genetic basis but that are cultivated worldwide. Its seeds have high nutritional value, but they naturally contain alkaloids with anti-nutritive factors. The study of its wild populations can be useful to find genotypes that contribute to higher climate resilience and greater yield under stressing environmental conditions. Using a common garden, we evaluated ecotypic differentiation in four natural populations from two contrasting latitudes in terms of plant biomass, seed mass and number, alkaloid content in seeds for the three main alkaloids present in the plant and seed colour, including its possible influence on post-dispersal predation. Correlations among traits were also assessed. We found differences among populations for all traits except final biomass. Northern populations had lighter seeds and a tendency to yield more seeds when they produced white seeds, compared to southern populations and variegated seeds. Regardless of latitude, populations showed differences in alkaloid concentration, with all three alkaloids found generally in high or low concentrations in each population. Proportion of white seed morphs varied in each population. Seed colour did not influence predator preference. In addition, white seed colour was related to a low alkaloid content. Our results evidence the existence of natural ecotypic differentiation in L. angustifolius not only due to latitudinal range, but also to local environmental factors. White seed coat colour could be used as a visual clue for identification of low-alkaloid genotypes, a priority trait in L. angustifolius breeding programmes.