Gene movement and genetic association with regional climate gradients in California valley oak (Quercus lobata Née) in the face of climate change

Plant conservation in the age of genome editing: opportunities and challenges

Numerous plant taxa are threatened by habitat destruction or overexploitation. To overcome these threats, new methods are urgently needed for rescuing threatened and endangered plant species. Here, we review the genetic consequences of threats to species populations. We highlight potential advantages of genome editing for mitigating negative effects caused by new pathogens and pests or climate change where other approaches have failed. We propose solutions to protect threatened plants using genome editing technology unless absolutely necessary. We further discuss the challenges associated with genome editing in plant conservation to mitigate the decline of plant diversity.

Environmental genome-wide association studies across precipitation regimes reveal that the E3 ubiquitin ligase MBR1 regulates plant adaptation to rainy environments

In an era characterized by rapidly changing and less-predictable weather conditions fueled by the climate crisis, understanding the mechanisms underlying local adaptation in plants is of paramount importance for the conservation of species. As the frequency and intensity of extreme precipitation events increase, so are the flooding events resulting from soil water saturation. The subsequent onset of hypoxic stress is one of the leading causes of crop damage and yield loss. By combining genomics and remote sensing data, it is now possible to probe natural plant populations that have evolved in different rainfall regimes and look for molecular adaptation to hypoxia. Here, using an environmental genome-wide association study (eGWAS) of 934 non-redundant georeferenced Arabidopsis ecotypes, we have identified functional variants of the gene MED25 BINDING RING-H2 PROTEIN 1 (MBR1). This gene encodes a ubiquitin-protein ligase that regulates MEDIATOR25 (MED25), part of a multiprotein complex that interacts with transcription factors that act as key drivers of the hypoxic response in Arabidopsis, namely the RELATED TO AP2 proteins RAP2.2 and RAP2.12. Through experimental validation, we show that natural variants of MBR1 have different effects on the stability of MED25 and, in turn, on hypoxia tolerance. This study also highlights the pivotal role of the MBR1/MED25 module in establishing a comprehensive hypoxic response. Our findings show that molecular candidates for plant environmental adaptation can be effectively mined from large datasets. This thus supports the need for integration of forward and reverse genetics with robust molecular physiology validation of outcomes.

History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

Phylogeographic analysis reveals extensive genetic variation of native grass Elymus nutans (Poaceae) on the Qinghai-Tibetan plateau

Introduction

Elymus nutans holds ecological and pastoral significance due to its adaptability and nutritional value, the Qinghai-Tibet Plateau (QTP) is a key hub for its genetic diversity. To conserve and harness its genetic resources in highland ecosystems, a thorough assessment is vital. However, a comprehensive phylogeographic exploration of E. nutans is lacking. The objective of this study was to unravel the genetic diversity, adaptation, and phylogenetics of E. nutans populations.

Methods

Encompassing 361 individuals across 35 populations, the species' genetic landscape and dynamic responses to diverse environments were decoded by using four chloroplast DNA (cpDNA) sequences and nine microsatellite markers derived from the transcriptome.

Results and discussion

This study unveiled a notable degree of genetic diversity in E. nutans populations at nuclear (I = 0.46, He = 0.32) and plastid DNA levels (Hd = 0.805, π = 0.67). Analysis via AMOVA highlighted genetic variation predominantly within populations. Despite limited isolation by distance (IBD), the Mekong-Salween Divide (MSD) emerged as a significant factor influencing genetic differentiation and conserving diversity. Furthermore, correlations were established between external environmental factors and effective alleles of three EST-SSRs (EN5, EN57 and EN80), potentially linked to glutathione S-transferases T1 or hypothetical proteins, affecting adaptation. This study deepens the understanding of the intricate relationship between genetic diversity, adaptation, and environmental factors within E. nutans populations on the QTP. The findings shed light on the species' evolutionary responses to diverse ecological conditions and contribute to a broader comprehension of plant adaptation mechanisms.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

The effect of environmental factors on the genetic differentiation of Cucurbita ficifolia populations based on whole-genome resequencing

BACKGROUND

Cucurbita ficifolia is one of the squash species most resistant to fungal pathogens, and has especially high resistance to melon Fusarium wilt. This species is therefore an important germplasm resource for the breeding of squash and melon cultivars.

RESULTS

Whole-genome resequencing of 223 individuals from 32 populations in Yunnan Province, the main cucurbit production area in China, was performed and 3,855,120 single-nucleotide polymorphisms (SNPs) and 1,361,000 InDels were obtained. SNP analysis suggested that levels of genetic diversity in C. ficifolia were high, but that different populations showed no significant genetic differentiation or geographical structure, and that individual C. ficifolia plants with fruit rinds of a similar color did not form independent clusters. A Mantel test conducted in combination with geographical distance and environmental factors suggested that genetic distance was not correlated with geographical distance, but had a significant correlation with environmental distance. Further associations between the genetic data and five environmental factors were analyzed using whole-genome association analysis. SNPs associated with each environmental factor were investigated and genes 250 kb upstream and downstream from associated SNPs were annotated. Overall, 15 marker-trait-associated SNPs (MTAs) and 293 genes under environmental selection were identified. The identified genes were involved in cell membrane lipid metabolism, macromolecular complexes, catalytic activity and other related aspects. Ecological niche modeling was used to simulate the distribution of C. ficifolia across time, from the present and into the future. We found that the area suitable for C. ficifolia changed with the changing climate in different periods.

CONCLUSIONS

Resequencing of the C. ficifolia accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs identified in this study suggest that environmental factors mediated the formation of the population structure of C. ficifolia in China. These SNPs and Indels might also contribute to the variation in important pathways of genes for important agronomic traits such as yield, disease resistance and stress tolerance. Moreover, the genome resequencing data and the genetic markers identified from 223 accessions provide insight into the genetic variation of the C. ficifolia germplasm and will facilitate a broad range of genetic studies.

Genomic analyses indicate resilience of a commercially and culturally important marine gastropod snail to climate change

Genomic vulnerability analyses are being increasingly used to assess the adaptability of species to climate change and provide an opportunity for proactive management of harvested marine species in changing oceans. Southeastern Australia is a climate change hotspot where many marine species are shifting poleward. The turban snail, Turbo militaris is a commercially and culturally harvested marine gastropod snail from eastern Australia. The species has exhibited a climate-driven poleward range shift over the last two decades presenting an ongoing challenge for sustainable fisheries management. We investigate the impact of future climate change on T. militaris using genotype-by-sequencing to project patterns of gene flow and local adaptation across its range under climate change scenarios. A single admixed, and potentially panmictic, demographic unit was revealed with no evidence of genetic subdivision across the species range. Significant genotype associations with heterogeneous habitat features were observed, including associations with sea surface temperature, ocean currents, and nutrients, indicating possible adaptive genetic differentiation. These findings suggest that standing genetic variation may be available for selection to counter future environmental change, assisted by widespread gene flow, high fecundity and short generation time in this species. We discuss the findings of this study in the content of future fisheries management and conservation.

Genotype-environment associations to reveal the molecular basis of environmental adaptation

A fundamental goal in plant biology is to identify and understand the variation underlying plants' adaptation to their environment. Climate change has given new urgency to this goal, as society aims to accelerate adaptation of ecologically important plant species, endangered plant species, and crops to hotter, less predictable climates. In the pre-genomic era, identifying adaptive alleles was painstaking work, leveraging genetics, molecular biology, physiology, and ecology. Now, the rise of genomics and new computational approaches may facilitate this research. Genotype-environment associations (GEAs) use statistical associations between allele frequency and environment of origin to test the hypothesis that allelic variation at a given gene is adapted to local environments. Researchers may scan the genome for GEAs to generate hypotheses on adaptive genetic variants (environmental genome-wide association studies). Despite the rapid adoption of these methods, many important questions remain about the interpretation of GEA findings, which arise from fundamental unanswered questions on the genetic architecture of adaptation and limitations inherent to association-based analyses. We outline strategies to ground GEAs in the underlying hypotheses of genetic architecture and better test GEA-generated hypotheses using genetics and ecophysiology. We provide recommendations for new users who seek to learn about the molecular basis of adaptation. When combined with a rigorous hypothesis testing framework, GEAs may facilitate our understanding of the molecular basis of climate adaptation for plant improvement.

Transcriptomic and metabolomic analysis of autumn leaf color change in Fraxinus angustifolia

Fraxinus angustifolia is a type of street tree and shade tree with ornamental value. It has a beautiful shape and yellow or reddish purple autumn leaves, but its leaf color formation mechanism and molecular regulation network need to be studied. In this study, we integrated the metabolomes and transcriptomes of stage 1 (green leaf) and stage 2 (red-purple leaf) leaves at two different developmental stages to screen differential candidate genes and metabolites related to leaf color variation. The results of stage 1 and stage 2 transcriptome analysis showed that a total of 5,827 genes were differentially expressed, including 2,249 upregulated genes and 3,578 downregulated genes. Through functional enrichment analysis of differentially expressed genes, we found that they were involved in flavonoid biosynthesis, phenylpropanoid biosynthesis, pigment metabolism, carotene metabolism, terpenoid biosynthesis, secondary metabolite biosynthesis, pigment accumulation, and other biological processes. By measuring the metabolites of Fraxinus angustifolia leaves, we found the metabolites closely related to the differentially expressed genes in two different periods of Fraxinus angustifolia, among which flavonoid compounds were the main differential metabolites. Through transcriptome and metabolomics data association analysis, we screened nine differentially expressed genes related to anthocyanins. Transcriptome and qRT-PCR results showed that these nine genes showed significant expression differences in different stages of the sample, and we speculate that they are likely to be the main regulatory factors in the molecular mechanism of leaf coloration. This is the first time that we have analyzed the transcriptome combination metabolome in the process of leaf coloration of Fraxinus angustifolia, which has important guiding significance for directional breeding of colored-leaf Fraxinus species and will also give new insights for enriching the landscape.

Landscape genetics of a sub-alpine toad: climate change predicted to induce upward range shifts via asymmetrical migration corridors

Climate change is expected to have a major hydrological impact on the core breeding habitat and migration corridors of many amphibians in the twenty-first century. The Yosemite toad (Anaxyrus canorus) is a species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. Despite living entirely on federal lands, it has recently faced severe extirpations, yet our understanding of climatic influences on population connectivity is limited. In this study, we used a previously published double-digest RADseq dataset along with numerous remotely sensed habitat features in a landscape genetics framework to answer two primary questions in Yosemite National Park: (1) Which fine-scale climate, topographic, soil, and vegetation features most facilitate meadow connectivity? (2) How is climate change predicted to influence both the magnitude and net asymmetry of genetic migration? We developed an approach for simultaneously modeling multiple toad migration paths, akin to circuit theory, except raw environmental features can be separately considered. Our workflow identified the most likely migration corridors between meadows and used the unique cubist machine learning approach to fit and forecast environmental models of connectivity. We identified the permuted modeling importance of numerous snowpack-related features, such as runoff and groundwater recharge. Our results highlight the importance of considering phylogeographic structure, and asymmetrical migration in landscape genetics. We predict an upward elevational shift for this already high-elevation species, as measured by the net vector of anticipated genetic movement, and a north-eastward shift in species distribution via the network of genetic migration corridors across the park.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Why Are Species’ Traits Weak Predictors of Range Shifts?

We examine the evidence linking species’ traits to contemporary range shifts and find they are poor predictors of range shifts that have occurred over decades to a century. We then discuss reasons for the poor performance of traits for describing interspecific variation in range shifts from two perspectives: ( a) factors associated with species’ traits that degrade range-shift signals stemming from the measures used for species’ traits, traits that are typically not analyzed, and the influence of phylogeny on range-shift potential and ( b) issues in quantifying range shifts and relating them to species’ traits due to imperfect detection of species, differences in the responses of altitudinal and latitudinal ranges, and emphasis on testing linear relationships between traits and range shifts instead of nonlinear responses. Improving trait-based approaches requires a recognition that traits within individuals interact in unexpected ways and that different combinations of traits may be functionally equivalent.

Waste not, want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single-nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black-capped vireo source-sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.

Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers

BACKGROUND

The malaria vector Anopheles minimus has been influenced by external stresses affecting the survival rate and vectorial capacity of the population. Since An. minimus habitats have continuously undergone ecological changes, this study aimed to determine the population genetic structure and the potential gene flow among the An. minimus populations in Thailand.

METHODS

Anopheles minimus was collected from five malaria transmission areas in Thailand using Centers for Disease Control and Prevention (CDC) light traps. Seventy-nine females from those populations were used as representative samples. The partial mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb) gene sequences were amplified and analyzed to identify species and determine the current population genetic structure. For the past population, we determined the population genetic structure from the 60 deposited COII sequences in GenBank of An. minimus collected from Thailand 20 years ago.

RESULTS

The current populations of An. minimus were genetically divided into two lineages, A and B. Lineage A has high haplotype diversity under gene flow similar to the population in the past. Neutrality tests suggested population expansion of An. minimus, with the detection of abundant rare mutations in all populations, which tend to arise from negative selection.

CONCLUSIONS

This study revealed that the population genetic structure of An. minimus lineage A was similar between the past and present populations, indicating high adaptability of the species. There was substantial gene flow between the eastern and western An. minimus populations without detection of significant gene flow barriers.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L

BACKGROUND

Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst.

RESULTS

We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech.

CONCLUSIONS

Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

Quercus Conservation Genetics and Genomics: Past, Present, and Future

Quercus species (oaks) have been an integral part of the landscape in the northern hemisphere for millions of years. Their ability to adapt and spread across different environments and their contributions to many ecosystem services is well documented. Human activity has placed many oak species in peril by eliminating or adversely modifying habitats through exploitative land usage and by practices that have exacerbated climate change. The goal of this review is to compile a list of oak species of conservation concern, evaluate the genetic data that is available for these species, and to highlight the gaps that exist. We compiled a list of 124 Oaks of Concern based on the Red List of Oaks 2020 and the Conservation Gap Analysis for Native U.S. Oaks and their evaluations of each species. Of these, 57% have been the subject of some genetic analysis, but for most threatened species (72%), the only genetic analysis was done as part of a phylogenetic study. While nearly half (49%) of published genetic studies involved population genetic analysis, only 16 species of concern (13%) have been the subject of these studies. This is a critical gap considering that analysis of intraspecific genetic variability and genetic structure are essential for designing conservation management strategies. We review the published population genetic studies to highlight their application to conservation. Finally, we discuss future directions in Quercus conservation genetics and genomics.

Adaptive divergence despite low effective population size in a peripherally isolated population of the pygmy rabbit, Brachylagus idahoensis

Local adaptation can occur when spatially separated populations are subjected to contrasting environmental conditions. Historically, understanding the genetic basis of adaptation has been difficult, but increased availability of genome-wide markers facilitates studies of local adaptation in non-model organisms of conservation concern. The pygmy rabbit (Brachylagus idahoensis) is an imperiled lagomorph that relies on sagebrush for forage and cover. This reliance has led to widespread population declines following reductions in the distribution of sagebrush, leading to geographic separation between populations. In this study, we used >20,000 single nucleotide polymorphisms, genotype-environment association methods, and demographic modeling to examine neutral genetic variation and local adaptation in the pygmy rabbit in Nevada and California. We identified 308 loci as outliers, many of which had functional annotations related to metabolism of plant secondary compounds. Likewise, patterns of spatial variation in outlier loci were correlated with landscape and climatic variables including proximity to streams, sagebrush cover, and precipitation. We found that populations in the Mono Basin of California probably diverged from other Great Basin populations during late Pleistocene climate oscillations, and that this region is adaptively differentiated from other regions in the southern Great Basin despite limited gene flow and low effective population size. Our results demonstrate that peripherally isolated populations can maintain adaptive divergence.

Genome-Wide Variation in DNA Methylation Predicts Variation in Leaf Traits in an Ecosystem-Foundational Oak Species

Epigenetic modifications such as DNA methylation are a potential mechanism for trees to respond to changing environments. However, it remains controversial the extent to which DNA methylation impacts ecologically important traits that influence fitness. In this study, we used reduced-representation bisulfite sequencing to associate genomic and epigenomic variation with seven phenotypic traits related to growth, leaf function, and disease susceptibility in 160 valley oak (Quercus lobata) saplings planted across two common gardens in California. We found that DNA methylation was associated with a significant fraction of phenotypic variance in plant height, leaf lobedness, powdery mildew infection, and trichome density. Two of the seven traits were significantly associated with DNA methylation in the CG context, three traits were significantly associated with CHG methylation, and two traits were significantly associated with CHH methylation. Notably, controlling for genomic variation in SNPs generally reduced the amount of trait variation explained by DNA methylation. Our results suggest that DNA methylation may serve as a useful biomarker to predict phenotypic variation in trees, though it remains unclear the degree to which DNA methylation is a causal mechanism driving phenotypic variation in forest tree species.

Comparing Genetic Diversity in Three Threatened Oaks

Genetic diversity is a critical resource for species’ survival during times of environmental change. Conserving and sustainably managing genetic diversity requires understanding the distribution and amount of genetic diversity (in situ and ex situ) across multiple species. This paper focuses on three emblematic and IUCN Red List threatened oaks (Quercus, Fagaceae), a highly speciose tree genus that contains numerous rare species and poses challenges for ex situ conservation. We compare the genetic diversity of three rare oak species—Quercus georgiana, Q. oglethorpensis, and Q. boyntonii—to common oaks; investigate the correlation of range size, population size, and the abiotic environment with genetic diversity within and among populations in situ; and test how well genetic diversity preserved in botanic gardens correlates with geographic range size. Our main findings are: (1) these three rare species generally have lower genetic diversity than more abundant oaks; (2) in some cases, small population size and geographic range correlate with genetic diversity and differentiation; and (3) genetic diversity currently protected in botanic gardens is inadequately predicted by geographic range size and number of samples preserved, suggesting non-random sampling of populations for conservation collections. Our results highlight that most populations of these three rare oaks have managed to avoid severe genetic erosion, but their small size will likely necessitate genetic management going forward.

Genetic, Morphological, and Environmental Differentiation of an Arid-Adapted Oak with a Disjunct Distribution

The patterns of genetic and morphological diversity of a widespread species can be influenced by environmental heterogeneity and the degree of connectivity across its geographic distribution. Here, we studied Quercus havardii Rydb., a uniquely adapted desert oak endemic to the Southwest region of the United States, using genetic, morphometric, and environmental datasets over various geographic scales to quantify differentiation and understand forces influencing population divergence. First, we quantified variation by analyzing 10 eastern and 13 western populations from the disjunct distribution of Q. havardii using 11 microsatellite loci, 17 morphological variables, and 19 bioclimatic variables. We then used regressions to examine local and regional correlations of climate with genetic variation. We found strong genetic, morphological and environmental differences corresponding with the large-scale disjunction of populations. Additionally, western populations had higher genetic diversity and lower relatedness than eastern populations. Levels of genetic variation in the eastern populations were found to be primarily associated with precipitation seasonality, while levels of genetic variation in western populations were associated with lower daily temperature fluctuations and higher winter precipitation. Finally, we found little to no observed environmental niche overlap between regions. Our results suggest that eastern and western populations likely represent two distinct taxonomic entities, each associated with a unique set of climatic variables potentially influencing local patterns of diversity.

Genetic data improves niche model discrimination and alters the direction and magnitude of climate change forecasts

Ecological niche models (ENMs) have classically operated under the simplifying assumptions that there are no barriers to gene flow, species are genetically homogeneous (i.e., no population-specific local adaptation), and all individuals share the same niche. Yet, these assumptions are violated for most broadly distributed species. Here, we incorporate genetic data from the widespread riparian tree species narrowleaf cottonwood (Populus angustifolia) to examine whether including intraspecific genetic variation can alter model performance and predictions of climate change impacts. We found that (1) P. angustifolia is differentiated into six genetic groups across its range from México to Canada and (2) different populations occupy distinct climate niches representing unique ecotypes. Comparing model discriminatory power, (3) all genetically informed ecological niche models (gENMs) outperformed the standard species-level ENM (3-14% increase in AUC; 1-23% increase in pROC). Furthermore, (4) gENMs predicted large differences among ecotypes in both the direction and magnitude of responses to climate change and (5) revealed evidence of niche divergence, particularly for the Eastern Rocky Mountain ecotype. (6) Models also predicted progressively increasing fragmentation and decreasing overlap between ecotypes. Contact zones are often hotspots of diversity that are critical for supporting species' capacity to respond to present and future climate change, thus predicted reductions in connectivity among ecotypes is of conservation concern. We further examined the generality of our findings by comparing our model developed for a higher elevation Rocky Mountain species with a related desert riparian cottonwood, P. fremontii. Together our results suggest that incorporating intraspecific genetic information can improve model performance by addressing this important source of variance. gENMs bring an evolutionary perspective to niche modeling and provide a truly "adaptive management" approach to support conservation genetic management of species facing global change.

LEA 3: Factor models in population genetics and ecological genomics with R

A major objective of evolutionary biology is to understand the processes by which organisms have adapted to various environments, and to predict the response of organisms to new or future conditions. The availability of large genomic and environmental data sets provides an opportunity to address those questions, and the R package LEA has been introduced to facilitate population and ecological genomic analyses in this context. By using latent factor models, the program computes ancestry coefficients from population genetic data and performs genotype-environment association analyses with correction for unobserved confounding variables. In this study, we present new functionalities of LEA, which include imputation of missing genotypes, fast algorithms for latent factor mixed models using multivariate predictors for genotype-environment association studies, population differentiation tests for admixed or continuous populations, and estimation of genetic offset based on climate models. The new functionalities are implemented in version 3.1 and higher releases of the package. Using simulated and real data sets, our study provides evaluations and examples of applications, outlining important practical considerations when analysing ecological genomic data in R.

Spatial Genetic Structure and Demographic History of the Dominant Forest Oak Quercus fabri Hance in Subtropical China

Oak trees (Quercus L.) are important models for estimating abiotic impacts on the population structure and demography of long life span tree species. In this study, we generated genetic data for 17 nuclear microsatellite loci in 29 natural populations of Quercus fabri to estimate the population genetic structure. We also integrated approximate Bayesian computation (ABC) and ecological niche analysis to infer the population differentiation processes and demographic history of this oak species. The genetic analyses indicated two genetic clusters across the 29 populations collected, where most approximately corresponded to the intraspecific differentiation among populations from western and eastern China, whereas admixed populations were mainly found in central mountains of China. The best model obtained from hierarchical ABC simulations suggested that the initial intraspecific divergence of Q. fabri potentially occurred during the late Pliocene (ca. 3.99 Ma) to form the two genetic clusters, and the admixed population group might have been generated by genetic admixture of the two differentiated groups at ca. 53.76 ka. Ecological analyses demonstrated clear differentiation among the Q. fabri population structures, and association estimations also indicated significant correlations between geography and climate with the genetic variation in this oak species. Our results suggest abiotic influences, including past climatic changes and ecological factors, might have affected the genetic differentiation and demographic history of Q. fabri in subtropical China.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

Landscape Genetics of Plants: Challenges and Opportunities

Dispersal is one of the most important but least understood processes in plant ecology and evolutionary biology. Dispersal of seeds maintains and establishes populations, and pollen and seed dispersal are responsible for gene flow within and among populations. Traditional views of dispersal and gene flow assume models that are governed solely by geographic distance and do not account for variation in dispersal vector behavior in response to heterogenous landscapes. Landscape genetics integrates population genetics with Geographic Information Systems (GIS) to evaluate the effects of landscape features on gene flow patterns (effective dispersal). Surprisingly, relatively few landscape genetic studies have been conducted on plants. Plants present advantages because their populations are stationary, allowing more reliable estimates of the effects of landscape features on effective dispersal rates. On the other hand, plant dispersal is intrinsically complex because it depends on the habitat preferences of the plant and its pollen and seed dispersal vectors. We discuss strategies to assess the separate contributions of pollen and seed movement to effective dispersal and to delineate the effects of plant habitat quality from those of landscape features that affect vector behavior. Preliminary analyses of seed dispersal for three species indicate that isolation by landscape resistance is a better predictor of the rates and patterns of dispersal than geographic distance. Rates of effective dispersal are lower in areas of high plant habitat quality, which may be due to the effects of the shape of the dispersal kernel or to movement behaviors of biotic vectors. Landscape genetic studies in plants have the potential to provide novel insights into the process of gene flow among populations and to improve our understanding of the behavior of biotic and abiotic dispersal vectors in response to heterogeneous landscapes.

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.

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.

Multiple axes of ecological vulnerability to climate change

Observed ecological responses to climate change are highly individualistic across species and locations, and understanding the drivers of this variability is essential for management and conservation efforts. While it is clear that differences in exposure, sensitivity, and adaptive capacity all contribute to heterogeneity in climate change vulnerability, predicting these features at macroecological scales remains a critical challenge. We explore multiple drivers of heterogeneous vulnerability across the distributions of 96 vegetation types of the ecologically diverse western US, using data on observed climate trends from 1948 to 2014 to highlight emerging patterns of change. We ask three novel questions about factors potentially shaping vulnerability across the region: (a) How does sensitivity to different climate variables vary geographically and across vegetation classes? (b) How do multivariate climate exposure patterns interact with these sensitivities to shape vulnerability patterns? (c) How different are these vulnerability patterns according to three widely implemented vulnerability paradigms-niche novelty (decline in modeled suitability), temporal novelty (standardized anomaly), and spatial novelty (inbound climate velocity)-each of which uses a distinct frame of reference to quantify climate departure? We propose that considering these three novelty paradigms in combination could help improve our understanding and prediction of heterogeneous climate change responses, and we discuss the distinct climate adaptation strategies connected with different combinations of high and low novelty across the three metrics. Our results reveal a diverse mosaic of climate change vulnerability signatures across the region's plant communities. Each of the above factors contributes strongly to this heterogeneity: climate variable sensitivity exhibits clear patterns across vegetation types, multivariate climate change data reveal highly diverse exposure signatures across locations, and the three novelty paradigms diverge widely in their climate change vulnerability predictions. Together, these results shed light on potential drivers of individualistic climate change responses and may help to inform effective management strategies.

Experimental DNA Demethylation Associates with Changes in Growth and Gene Expression of Oak Tree Seedlings

Epigenetic modifications such as DNA methylation, where methyl groups are added to cytosine base pairs, have the potential to impact phenotypic variation and gene expression, and could influence plant response to changing environments. One way to test this impact is through the application of chemical demethylation agents, such as 5-Azacytidine, which inhibit DNA methylation and lead to a partial reduction in DNA methylation across the genome. In this study, we treated 5-month-old seedlings of the tree, Quercus lobata, with foliar application of 5-Azacytidine to test whether a reduction in genome-wide methylation would cause differential gene expression and change phenotypic development. First, we demonstrate that demethylation treatment led to 3–6% absolute reductions and 6.7–43.2% relative reductions in genome-wide methylation across CG, CHG, and CHH sequence contexts, with CHH showing the strongest relative reduction. Seedlings treated with 5-Azacytidine showed a substantial reduction in new growth, which was less than half that of control seedlings. We tested whether this result could be due to impact of the treatment on the soil microbiome and found minimal differences in the soil microbiome between two groups, although with limited sample size. We found no significant differences in leaf fluctuating asymmetry (i.e., deviations from bilateral symmetry), which has been found in other studies. Nonetheless, treated seedlings showed differential expression of a total of 23 genes. Overall, this study provides initial evidence that DNA methylation is involved in gene expression and phenotypic variation in seedlings and suggests that removal of DNA methylation affects plant development.

Adaptational lag to temperature in valley oak (Quercus lobata) can be mitigated by genome-informed assisted gene flow

Climate change over the next century is predicted to cause widespread maladaptation in natural systems. This prediction, as well as many sustainable management and conservation practices, assumes that species are adapted to their current climate. However, this assumption is rarely tested. Using a large-scale common garden experiment combined with genome-wide sequencing, we found that valley oak (Quercus lobata), a foundational tree species in California ecosystems, showed a signature of adaptational lag to temperature, with fastest growth rates occurring at cooler temperatures than populations are currently experiencing. Future warming under realistic emissions scenarios was predicted to lead to further maladaptation to temperature and reduction in growth rates for valley oak. We then identified genotypes predicted to grow relatively fast under warmer temperatures and demonstrated that selecting seed sources based on their genotype has the potential to mitigate predicted negative consequences of future climate warming on growth rates in valley oak. These results illustrate that the belief of local adaptation underlying many management and conservation practices, such as using local seed sources for restoration, may not hold for some species. If contemporary adaptational lag is commonplace, we will need new approaches to help alleviate predicted negative consequences of climate warming on natural systems. We present one such approach, "genome-informed assisted gene flow," which optimally matches individuals to future climates based on genotype-phenotype-environment associations.

Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages

BACKGROUND

Understanding the origin of genetic variation is the key to predict how species will respond to future climate change. The genus Quercus is a species-rich and ecologically diverse woody genus that dominates a wide range of forests and woodland communities of the Northern Hemisphere. Quercus thus offers a unique opportunity to investigate how adaptation to environmental changes has shaped the spatial genetic structure of closely related lineages. Furthermore, Quercus provides a deep insight into how tree species will respond to future climate change. This study investigated whether closely related Quercus lineages have similar spatial genetic structures and moreover, what roles have their geographic distribution, ecological tolerance, and historical environmental changes played in the similar or distinct genetic structures.

RESULTS

Despite their close relationships, the three main oak lineages (Quercus sections Cyclobalanopsis, Ilex, and Quercus) have different spatial genetic patterns and occupy different climatic niches. The lowest level and most homogeneous pattern of genetic diversity was found in section Cyclobalanopsis, which is restricted to warm and humid climates. The highest genetic diversity and strongest geographic genetic structure were found in section Ilex, which is due to their long-term isolation and strong local adaptation. The widespread section Quercus is distributed across the most heterogeneous range of environments; however, it exhibited moderate haplotype diversity. This is likely due to regional extinction during Quaternary climatic fluctuation in Europe and North America.

CONCLUSIONS

Genetic variations of sections Ilex and Quercus were significantly predicted by geographic and climate variations, while those of section Cyclobalanopsis were poorly predictable by geographic or climatic diversity. Apart from the different historical environmental changes experienced by different sections, variation of their ecological or climatic tolerances and physiological traits induced varying responses to similar environment changes, resulting in distinct spatial genetic patterns.

Geographic isolation and climatic variability contribute to genetic differentiation in fragmented populations of the long-lived subalpine conifer Pinus cembra L. in the western Alps

Background

Genetic processes shape the modern-day distribution of genetic variation within and between populations and can provide important insights into the underlying mechanisms of evolution. The resulting genetic variation is often unequally partitioned within species’ distribution range and especially large differences can manifest at the range limit, where population fragmentation and isolation play a crucial role in species survival. Despite several molecular studies investigating the genetic diversity and differentiation of European Alpine mountain forests, the climatic and demographic constrains which influence the genetic processes are often unknown. Here, we apply non-coding microsatellite markers to evaluate the sporadic peripheral and continuous populations of cembra pine (Pinus cembra L.), a long-lived conifer species that inhabits the subalpine treeline ecotone in the western Alps to investigate how the genetic processes contribute to the modern-day spatial distribution. Moreover, we corroborate our findings with paleoecological records, micro and macro-remains, to infer the species’ possible glacial refugia and expansion scenarios.

Results

Four genetically distinct groups were identified, with Bayesian and F_ST based approaches, across the range of the species, situated in the northern, inner and south-western Alps. We found that genetic differentiation is substantially higher in marginal populations than at the center of the range, and marginal stands are characterized by geographic and genetic isolation due to spatial segregation and restricted gene flow. Moreover, multiple matrix regression approaches revealed effects of climatic heterogeneity in species’ spatial genetic pattern. Also, population stability tests indicated that all populations had experienced a severe historical bottleneck, no heterozygosity excess was detected, suggesting that more recently population sizes have remained relatively stable.

Conclusions

Our study demonstrated that cembra pine might have survived in multiple glacial refugia and subsequently recolonized the Alps by different routes. Modern-day marginal populations, at the edge of the species’ range, could maintain stable sizes over long periods without inbreeding depression and preserve high amounts of genetic variation. Moreover, our analyses indicate that climatic variability has played a major role in shaping differentiation, in addition to past historical events such as migration and demographic changes.

Abiotic genetic adaptation in the Fagaceae

Fagaceae can be found in tropical and temperate regions and contain species of major ecological and economic importance. In times of global climate change, tree populations need to adapt to rapidly changing environmental conditions. The predicted warmer and drier conditions will potentially result in locally maladapted populations. There is evidence that major genera of the Fagaceae are already negatively affected by climate change-related factors such as drought and associated biotic stressors. Therefore, knowledge of the mechanisms underlying adaptation is of great interest. In this review, we summarise current literature related to genetic adaptation to abiotic environmental conditions. We begin with an overview of genetic diversity in Fagaceae species and then summarise current knowledge related to drought stress tolerance, bud burst timing and frost tolerance in the Fagaceae. Finally, we discuss the role of hybridisation, epigenetics and phenotypic plasticity in adaptation.

Regional patterns of genetic structure and environmental differentiation in willow populations (Salix humboldtiana Willd.) from Central Mexico

AIM

To infer the geological and climatic factors that have shaped the genetic diversity and structure of a willow species (Salix humboldtiana) in three basins of Central Mexico.

LOCATION

Central Mexico.

METHODS

We collected samples from 11 populations across two hydrological basins (Balsas and Lerma) and one population from another basin (Ameca) within the Mexican Central Plateau (MCP). Individuals were analyzed using sequences of two chloroplast DNA (cpDNA) regions and eight nuclear simple sequence repeats (nSSR). Population genetic diversity and structure were determined from these data. To evaluate whether genetic structure was associated with ecological niche differentiation, we determined whether there is niche equivalence, overlap, or divergence between the Balsas and Lerma basins. Also, we evaluated the relative contributions of geographic distribution and climatic variation on population genetic structuring through redundancy analysis (RDA) and partial RDA.

RESULTS

Both cpDNA and nSSRs data indicated the presence of three highly differentiated genetic groups, mostly geographically congruent with the three main hydrological basins. According to nSSRs, the three genetic groups can be further subdivided into eight subgroups corresponding to different rivers within the main basins. The niche equivalency test showed that the niches of the species in the Balsas and Lerma basins are significantly nonequivalent. The RDA indicated a significant association of genetic variation among populations with climate variables (particularly those related to the precipitation regime), while controlling for geographic distribution.

MAIN CONCLUSIONS

The genetic structure of S. humboldtiana is strongly associated with the historical and current geological configuration of the basins and the rivers within basins. The observed hierarchical genetic differentiation can be due to gene flow limitation resulting from physical barriers to the dispersal of S. humboldtiana, but also to some degree of isolation by environment, as suggested by the significant association between genetic variation among populations and precipitation regime.

Habitat preference differentiates the Holocene range dynamics but not barrier effects on two sympatric, congeneric trees (Tristaniopsis, Myrtaceae)

Niche partitioning can lead to differences in the range dynamics of plant species through its impacts on habitat availability, dispersal, or selection for traits that affect colonization and persistence. We investigated whether niche partitioning into upland and riparian habitats differentiates the range dynamics of two closely related and sympatric eastern Australian trees: the mountain water gum (Tristaniopsis collina) and the water gum (T. laurina). Using genomic data from SNP genotyping of 480 samples, we assessed the impact of biogeographic barriers and tested for signals of range expansion. Circuit theory was used to model isolation-by-resistance across three palaeo-environment scenarios: the Last Glacial Maximum, the Holocene Climate Optimum and present-day (1950-2014). Both trees showed similar genetic structure across historically dry barriers, despite evidence of significant environmental niche differentiation and different post-glacial habitat shifts. Tristaniopsis collina exhibits the signature of serial founder effects consistent with recent or rapid range expansion, whilst T. laurina has genetic patterns consistent with long-term persistence in geographically isolated populations despite occupying a broader bioclimatic niche. We found the minor influence of isolation-by-resistance on both species, though other unknown factors appear to shape genetic variation. We postulate that specialized recruitment traits (adapted to flood-disturbance regimes) rather than habitat availability limited post-glacial range expansion in T. laurina. Our findings indicate that niche breadth does not always facilitate range expansion through colonization and migration across barriers, though it can promote long-term persistence in situ.

Pervasive admixture between eucalypt species has consequences for conservation and assisted migration

Conservation management often uses information on genetic population structure to assess the importance of local provenancing for ecological restoration and reintroduction programs. For species that do not exhibit complete reproductive isolation, the estimation of population genetic parameters may be influenced by the extent of admixture. Therefore, to avoid perverse outcomes for conservation, genetically informed management strategies must determine whether hybridization between species is relevant, and the extent to which observed population genetic patterns are shaped by interspecific versus intraspecific gene flow. We used genotyping by sequencing to identify over 2,400 informative single nucleotide polymorphisms across 18 populations of Eucalyptus regnans F. Muell., a foundation tree species of montane forests in south-eastern Australia. We used these data to determine the extent of hybridization with another species, Eucalyptus obliqua L'Hér., and investigate how admixture influences genetic diversity parameters, by estimating metrics of genetic diversity and examining population genetic structure in datasets with and without admixed individuals. We found hybrid individuals at all sites and two highly introgressed populations. Hybrid individuals were not distributed evenly across environmental gradients, with logistic regression identifying hybrids as being associated with temperature. Removal of hybrids resulted in increases in genetic differentiation (F ST), expected heterozygosity, observed heterozygosity and the inbreeding coefficient, and different patterns of isolation by distance. After removal of hybrids and introgressed populations, mountain ash showed very little population genetic structure, with a small effect of isolation by distance, and very low global F ST(0.03). Our study shows that, in plants, decisions around provenancing of individuals for restoration depend on knowledge of whether hybridization is influencing population genetic structure. For species in which most genetic variation is held within populations, there may be little benefit in planning conservation strategies around environmental adaptation of seed sources. The possibility for adaptive introgression may also be relevant when species regularly hybridize.

Landscape genomics provides evidence of climate-associated genetic variation in Mexican populations of Quercus rugosa

Local adaptation is a critical evolutionary process that allows plants to grow better in their local compared to non-native habitat and results in species-wide geographic patterns of adaptive genetic variation. For forest tree species with a long generation time, this spatial genetic heterogeneity can shape the ability of trees to respond to rapid climate change. Here, we identify genomic variation that may confer local environmental adaptations and then predict the extent of adaptive mismatch under future climate as a tool for forest restoration or management of the widely distributed high-elevation oak species Quercus rugosa in Mexico. Using genotyping by sequencing, we identified 5,354 single nucleotide polymorphisms (SNPs) genotyped from 103 individuals across 17 sites in the Trans-Mexican Volcanic Belt, and, after controlling for neutral genetic structure, we detected 74 F ST outlier SNPs and 97 SNPs associated with climate variation. Then, we deployed a nonlinear multivariate model, Gradient Forests, to map turnover in allele frequencies along environmental gradients and predict areas most sensitive to climate change. We found that spatial patterns of genetic variation were most strongly associated with precipitation seasonality and geographic distance. We identified regions of contemporary genetic and climatic similarities and predicted regions where future populations of Q. rugosa might be at risk due to high expected rate of climate change. Our findings provide preliminary details for future management strategies of Q. rugosa in Mexico and also illustrate how a landscape genomic approach can provide a useful tool for conservation and resource management strategies.

Invoking adaptation to decipher the genetic legacy of past climate change

Persistence of natural populations during periods of climate change is likely to depend on migration (range shifts) or adaptation. These responses were traditionally considered discrete processes and conceptually divided into the realms of ecology and evolution. In a milestone paper, Davis and Shaw (2001) Science 292:673 argued that the interplay of adaptation and migration was central to biotic responses to Quaternary climate, but since then there has been no synthesis of efforts made to set up this research program. Here we review some of the salient findings from molecular genetic studies assessing ecological and evolutionary responses to Quaternary climate change. These studies have revolutionized our understanding of population processes associated with past species migration. However, knowledge remains limited about the role of natural selection for local adaptation of populations to Quaternary environmental fluctuations and associated range shifts, and for the footprints this might have left on extant populations. Next-generation sequencing technologies, high-resolution paleoclimate analyses, and advances in population genetic theory offer an unprecedented opportunity to test hypotheses about adaptation through time. Recent population genomics studies have greatly improved our understanding of the role of contemporary adaptation to local environments in shaping spatial patterns of genetic diversity across modern-day landscapes. Advances in this burgeoning field provide important conceptual and methodological bases to decipher the historical role of natural selection and assess adaptation to past environmental variation. We suggest that a process called "temporal conditional neutrality" has taken place: some alleles favored in glacial environments become selectively neutral in modern-day conditions, whereas some alleles that had been neutral during glacial periods become under selection in modern environments. Building on this view, we present a new integrative framework for addressing the interplay of demographic and adaptive evolutionary responses to Quaternary climate dynamics, the research agenda initially envisioned by Davis and Shaw (2001) Science 292:673.

Genetic diversity of Avena ventricosa populations along an ecogeographical transect in Cyprus is correlated to environmental variables

Avena ventricosa Balansa ex Coss. is considered the C-genome donor of the cultivated hexaploid oat and is a ‘priority’ species for conservation, since it has limited geographic distribution and the only recorded populations in Europe are present in Cyprus. The current study attempts to characterize the genetic structure and fragmentation of the species via the application of genotypic markers. It was revealed that the genetic variety was mainly allocated among the populations collected, since clustering obtained was according to the geographic origin of the samples and the habitat. Species distribution modeling showed that the most important climatic variable defining A. ventricosa distribution is the mean diurnal temperature. Furthermore, significant association of the genetic structure to environmental variables was detected; overall, a negative association to precipitation was confirmed, while significant correlations of genetic structure and the temperature at the time of anthesis and germination were established. The safeguarding of this valuable genetic resource is discussed.

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change.

Drivers of synchrony of acorn production in the valley oak (Quercus lobata) at two spatial scales

We investigated spatial synchrony of acorn production by valley oaks (Quercus lobata) among individual trees at the within-population, local level and at the among-population, statewide level spanning the geographic range of the species. At the local level, the main drivers of spatial synchrony were water availability and flowering phenology of individual trees, while proximity, temperature differences between trees, and genetic similarity failed to explain a significant proportion of variance in spatial synchrony. At the statewide level, annual rainfall was the primary driver, while proximity was significant by itself but not when controlling for rainfall; genetic similarity was again not significant. These results support the hypothesis that environmental factors, the Moran effect, are key drivers of spatial synchrony in acorn production at both small and large geographic scales. The specific environmental factors differed depending on the geographic scale, but were in both cases related to water availability. In addition, flowering phenology, potentially affecting either density-independent pollination failure (the pollination Moran effect) or density-dependent pollination efficiency (pollen coupling), plays a key role in driving spatial synchrony at the local geographic scale.

Vulnerability of eastern US tree species to climate change

Climate change is expected to alter the distribution of tree species because of critical environmental tolerances related to growth, mortality, reproduction, disturbances, and biotic interactions. How this is realized in 21st century remains uncertain, in large part due to limitations on plant migration and the impacts of landscape fragmentation. Understanding these changes is of particular concern for forest management, which requires information at an appropriately fine spatial resolution. Here we provide a framework and application for tree species vulnerability to climate change in the eastern United States that accounts for influential drivers of future distributions. We used species distribution models to project changes in habitat suitability at 800 m for 40 tree species that vary in physiology, range, and environmental niche. We then developed layers of adaptive capacity based on migration potential, forest fragmentation, and propagule pressure. These were combined into metrics of vulnerability, including an overall index and spatially explicit categories designed to inform management. Despite overall favorable changes in suitability, the majority of species and the landscape were considered vulnerable to climate change. Vulnerability was significantly exacerbated by projections of pests and pathogens for some species. Northern and high-elevation species tended to be the most vulnerable. There were, however, some notable areas of particular resilience, including most of West Virginia. Our approach combines some of the most important considerations for species vulnerability in a straightforward framework, and can be used as a tool for managers to prioritize species, areas, and actions.