Explore the vast range of titles available.

When confronted with an adaptive challenge, such as extreme temperature, closely related species frequently evolve similar phenotypes using the same genes. Although such repeated evolution is thought to be less likely in highly polygenic traits and distantly related species, this has not been tested at the genome scale. We performed a population genomic study of convergent local adaptation among two distantly related species, lodgepole pine and interior spruce. We identified a suite of 47 genes, enriched for duplicated genes, with variants associated with spatial variation in temperature or cold hardiness in both species, providing evidence of convergent local adaptation despite 140 million years of separate evolution. These results show that adaptation to climate can be genetically constrained, with certain key genes playing nonredundant roles.

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.

Background High-throughput re-sequencing is rapidly becoming the method of choice for studies of neutral and adaptive processes in natural populations across taxa. As re-sequencing the genome of large numbers of samples is still cost-prohibitive in many cases, methods for genome complexity reduction have been developed in attempts to capture most ecologically-relevant genetic variation. One of these approaches is sequence capture, in which oligonucleotide baits specific to genomic regions of interest are synthesized and used to retrieve and sequence those regions. Results We used sequence capture to re-sequence most predicted exons, their upstream regulatory regions, as well as numerous random genomic intervals in a panel of 48 genotypes of the angiosperm tree Populus trichocarpa (black cottonwood, or ‘poplar’). A total of 20.76Mb (5%) of the poplar genome was targeted, corresponding to 173,040 baits. With 12 indexed samples run in each of four lanes on an Illumina HiSeq instrument (2x100 paired-end), 86.8% of the bait regions were on average sequenced at a depth ≥10X. Few off-target regions (>250bp away from any bait) were present in the data, but on average ~80bp on either side of the baits were captured and sequenced to an acceptable depth (≥10X) to call heterozygous SNPs. Nucleotide diversity estimates within and adjacent to protein-coding genes were similar to those previously reported in Populus spp., while intergenic regions had higher values consistent with a relaxation of selection. Conclusions Our results illustrate the efficiency and utility of sequence capture for re-sequencing highly heterozygous tree genomes, and suggest design considerations to optimize the use of baits in future studies.

Genecological studies in widespread tree species have revealed steep genetic clines along environmental gradients for climate-related traits. In a changing climate, the ecological and economic importance of conifers necessitates an appraisal of how molecular genetic variation shapes quantitative trait variation, and one of the most promising approaches to answer this question is association mapping. We phenotyped a wide collection of 410 individuals of the widely distributed conifer Sitka spruce rangewide (Picea sitchensis) for budset timing and autumn cold hardiness, and genotyped these individuals for a panel of 768 single nucleotide polymorphisms (SNPs) representing > 200 expressed nuclear genes. After correcting for population structure, associations were detected in 28 of the candidate genes, which cumulatively explained 28 and 34% of the phenotypic variance in cold hardiness and budset, respectively. Most notable among the associations were five genes putatively involved in light signal transduction, the key pathway regulating autumn growth cessation in perennials. Many SNPs with phenotypic associations were also correlated with at least one climate variable. This study represents a significant step toward the goal of characterizing the genomic basis of adaptation to local climate in conifers, and provides an important resource for breeding and conservation genetics in a changing climate.

Adaptation to climate across latitude and altitude reflects shared climatic constraints, which may lead to parallel adaptation. However, theory predicts that higher gene flow should favor more concentrated genomic architectures, which would lead to fewer locally maladapted recombinants. We used exome capture to resequence the gene space along a latitudinal and two altitudinal transects in the model tree Populus trichocapra. Adaptive trait phenotyping was coupled with FST outlier tests and sliding window analysis to assess the degree of parallel adaptation as well as the genomic distribution of outlier loci. Up to 51% of outlier loci overlapped between transect pairs and up to 15% of these loci overlapped among all three transects. Genomic clustering of adaptive loci was more pronounced for altitudinal than latitudinal transects. In both altitudinal transects, there was a larger number of these ‘islands of divergence’, which were on average longer and included several of exceptional physical length. Our results suggest that recapitulation of genetic clines over latitude and altitude involves extensive parallelism, but that steep altitudinal clines generate islands of divergence. This suggests that physical proximity of genes in coadapted complexes may buffer against the movement of maladapted alleles from geographically proximal but climatically distinct populations.

Species distribution models predict a wholesale redistribution of trees in the next century, yet migratory responses necessary to spatially track climates far exceed maximum post‐glacial rates. The extent to which populations will adapt will depend upon phenotypic variation, strength of selection, fecundity, interspecific competition, and biotic interactions. Populations of temperate and boreal trees show moderate to strong clines in phenology and growth along temperature gradients, indicating substantial local adaptation. Traits involved in local adaptation appear to be the product of small effects of many genes, and the resulting genotypic redundancy combined with high fecundity may facilitate rapid local adaptation despite high gene flow. Gene flow with preadapted alleles from warmer climates may promote adaptation and migration at the leading edge, while populations at the rear will likely face extirpation. Widespread species with large populations and high fecundity are likely to persist and adapt, but will likely suffer adaptational lag for a few generations. As all tree species will be suffering lags, interspecific competition may weaken, facilitating persistence under suboptimal conditions. Species with small populations, fragmented ranges, low fecundity, or suffering declines due to introduced insects or diseases should be candidates for facilitated migration.

American chestnut was once a foundation species of eastern North American forests, but was rendered functionally extinct in the early 20th century by an exotic fungal blight (Cryphonectria parasitica). Over the past 30 years, the American Chestnut Foundation (TACF) has pursued backcross breeding to generate hybrids that combine the timber‐type form of American chestnut with the blight resistance of Chinese chestnut based on a hypothesis of major gene resistance. To accelerate selection within two backcross populations that descended from two Chinese chestnuts, we developed genomic prediction models for five presence/absence blight phenotypes of 1,230 BC3F2 selection candidates and average canker severity of their BC3F3 progeny. We also genotyped pure Chinese and American chestnut reference panels to estimate the proportion of BC3F2 genomes inherited from parent species. We found that genomic prediction from a method that assumes an infinitesimal model of inheritance (HBLUP) has similar accuracy to a method that tends to perform well for traits controlled by major genes (Bayes C). Furthermore, the proportion of BC3F2 trees' genomes inherited from American chestnut was negatively correlated with the blight resistance of these trees and their progeny. On average, selected BC3F2 trees inherited 83% of their genome from American chestnut and have blight resistance that is intermediate between F1 hybrids and American chestnut. Results suggest polygenic inheritance of blight resistance. The blight resistance of restoration populations will be enhanced through recurrent selection, by advancing additional sources of resistance through fewer backcross generations, and by potentially by breeding with transgenic blight‐tolerant trees.

Climate is the primary driver of the distribution of tree species worldwide, and the potential for adaptive evolution will be an important factor determining the response of forests to anthropogenic climate change. Although association mapping has the potential to improve our understanding of the genomic underpinnings of climatically relevant traits, the utility of adaptive polymorphisms uncovered by such studies would be greatly enhanced by the development of integrated models that account for the phenotypic effects of multiple single-nucleotide polymorphisms (SNPs) and their interactions simultaneously. We previously reported the results of association mapping in the widespread conifer Sitka spruce (Picea sitchensis). In the current study we used the recursive partitioning algorithm ‘Random Forest’ to identify optimized combinations of SNPs to predict adaptive phenotypes. After adjusting for population structure, we were able to explain 37% and 30% of the phenotypic variation, respectively, in two locally adaptive traits—autumn budset timing and cold hardiness. For each trait, the leading five SNPs captured much of the phenotypic variation. To determine the role of epistasis in shaping these phenotypes, we also used a novel approach to quantify the strength and direction of pairwise interactions between SNPs and found such interactions to be common. Our results demonstrate the power of Random Forest to identify subsets of markers that are most important to climatic adaptation, and suggest that interactions among these loci may be widespread.

Societal Impact Statement Over four billion American chestnut trees have been killed as a result of an introduced pathogen, the chestnut blight fungus. Recently, transgenic blight‐tolerant American chestnut trees have been produced by inserting a gene from wheat into the American chestnut genome. Pending federal approval to use these transgenic trees for large‐scale forest restoration, this would be the first instance where a transgenic approach has been used to restore a tree species that has been rendered functionally extinct by an introduced pathogen. With the help of citizen scientists, we estimate that large‐scale forest restoration using blight‐tolerant American chestnut trees is possible within the next few decades. Summary Breeding transgenic blight‐tolerant American chestnuts with susceptible wild‐type (WT) trees is potentially an efficient method to rescue the genetic diversity and adaptive capacity of the American chestnut population for large‐scale restoration. To develop a breeding plan to diversify a transgenic blight‐tolerant population, we simulated pedigrees to estimate inbreeding coefficients and effective population size in scenarios involving outcrossing 1–4 transgenic founders to a maximum of 1,500 WT trees over 1–5 generations. We also simulated marker‐assisted introgression scenarios to minimize the extent of the transgenic founder genome, especially on the transgene carrier chromosome. Simulations suggest that the effective population size may be increased to >500, and the average inbreeding coefficient reduced to <0.01, by outcrossing a single transgenic founder over five generations to 2, 25, 50, 150, and 450 (677 total) WT parents. Three generations of marker‐assisted introgression is predicted to decrease the length of the founder genome to between 7% and 13% of the transgene carrier chromosome length as compared to 42% with event selection (ES) only. Transgenic outcross selections may be intercrossed to select progeny homozygous for the transgene for planting in seed orchards. A diversified population of transgenic blight‐tolerant American chestnut is estimated to be available for use in large‐scale forest restoration 20–35 years after federal approval to distribute the trees. In contrast, trees from earlier generations would be available almost immediately after federal approval for personal or horticultural plantings. Methods to accelerate pollen production and outcrossing are discussed. Over four billion American chestnut trees have been killed as a result of an introduced pathogen, the chestnut blight fungus. Recently, transgenic blight‐tolerant American chestnut trees have been produced by inserting a gene from wheat into the American chestnut genome. Pending federal approval to use these transgenic trees for large‐scale forest restoration, this would be the first instance where a transgenic approach has been used to restore a tree species that has been rendered functionally extinct by an introduced pathogen. With the help of citizen scientists, we estimate that large‐scale forest restoration using blight‐tolerant American chestnut trees is possible within the next few decades.

Gene flow and effective population size (Ne) should depend on a population's position within its range: those near the edges are expected to have smaller Ne and lower relative emigration rates, whereas those nearer the centre should have larger Ne and higher relative emigration rates. In species with continuous ranges, this phenomenon may limit the ability of peripheral populations to respond to divergent selection. Here, we employ Sitka spruce as a model to test these predictions. We previously genotyped 339 single nucleotide polymorphisms (SNPs) in 410 individuals from 13 populations, and used these data to identify putative targets of divergent selection, as well as to explore the extent to which central–peripheral structure may impede adaptation. Fourteen SNPs had outlier FST estimates suggestive of divergent selection, of which nine were previously associated with phenotypic variation in adaptive traits (timing of autumn budset and cold hardiness). Using coalescent simulations, we show that populations from near the centre of the range have higher effective populations sizes than those from the edges, and that central populations contribute more migrants to marginal populations than the reverse. Our results suggest that while divergent selection appears to have shaped allele frequencies among populations, asymmetrical movement of alleles from the centre to the edges of the species range may affect the adaptive capacity of peripheral populations. In southern peripheral populations, the movement of cold-adapted alleles from the north represents a significant impediment to adaptation under climate change, while in the north, movement of warm-adapted alleles from the south may enhance adaptation.