A genome-guided strategy for climate resilience in American chestnut restoration populations.

American chestnut (Castanea dentata) is a deciduous tree species of eastern North America that was decimated by the introduction of the chestnut blight fungus (Cryphonectria parasitica) in the early 20th century. Although millions of American chestnuts survive as root collar sprouts, these trees rarely reproduce. Thus, the species is considered functionally extinct. American chestnuts with improved blight resistance have been developed through interspecific hybridization followed by conspecific backcrossing, and by genetic engineering. Incorporating adaptive genomic diversity into these backcross families and transgenic lines is important for restoring the species across broad climatic gradients. To develop sampling recommendations for ex situ conservation of wild adaptive genetic variation, we coupled whole-genome resequencing of 384 stump sprouts with genotype-environment association analyses and found that the species range can be subdivided into three seed zones characterized by relatively homogeneous adaptive allele frequencies. We estimated that 21 to 29 trees per seed zone will need to be conserved to capture most extant adaptive diversity. We also resequenced the genomes of 269 backcross trees to understand the extent to which the breeding program has already captured wild adaptive diversity, and to estimate optimal reintroduction sites for specific families on the basis of their adaptive portfolio and future climate projections. Taken together, these results inform the development of an ex situ germplasm conservation and breeding plan to target blight-resistant breeding populations to specific environments and provides a blueprint for developing restoration plans for other imperiled tree species.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage.

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

Validation of an Alternative Small Stem Assay for Blight Resistance in Chestnut Seedlings and Recommendations for Broader Use.

We evaluated an alternative small stem assay (AltSSA) for blight resistance in backcross hybrid chestnut trees (Castanea dentata/mollissima). Whereas standard small stem assays (SSAs) are done by inoculating small incisions in stems, in our AltSSA, 4- to 5-mm stems are cut off, and the exposed (living) stem tips are inoculated with discs of Cryphonectria parasitica inoculum and temporarily covered with plastic sleeves. Intended primarily for forward selection, this method was designed to be easy to implement, to consistently induce cankering, and to better enable seedling recovery via the development of lateral shoots from the lower stem. After 90+ days, cankers are evaluated and removed, and seedlings are prepared for out-planting. Previous results showed that AltSSAs performed at least as well as a common SSA method in distinguishing resistant and susceptible types. In this follow-up analysis of 35 lines of backcross seedlings studied in 2020 and 2021, we showed that mean orange zone canker length (OZCL) and a multifactor principal components analysis-based blight resistance index gave results consistent with predictions derived from two methods of blight resistance phenotyping and percentage of American chestnut ancestry of the parents of each line. As expected, based upon the apparent polygenic inheritance of blight resistance in backcross chestnut trees, mean OZCL of backcross families ranged from intermediate (F1 hybrid-level) to low (wild-type American chestnut-level). Consistent with prior results, canker production was near 100%, survivorship after out-planting was very high, and postinoculation stem dieback was not apparently related to the stem tip inoculations. Altogether, these results suggest that the AltSSA is a viable method for early detection of relative blight resistance in seedlings and may enable a reduction in the numbers of trees out-planted and placed under care for long-term evaluation and breeding. Thus, the AltSSA can prevent time, resources, and orchard space from being used on susceptible trees.

Frozen in time: Rangewide genomic diversity, structure, and demographic history of relict American chestnut populations.

American chestnut (Castanea dentata) was once the most economically and ecologically important hardwood species in the eastern United States. In the first half of the 20th century, an exotic fungal pathogen-Cryphonectria parasitica-decimated the species, killing billions of chestnut trees. Two approaches to developing blight-resistant American chestnut populations show promise, but both will require introduction of adaptive genomic diversity from wild germplasm to produce diverse, locally adapted restoration populations. Here we characterize population structure, demographic history, and genomic diversity in a range-wide sample of 384 wild American chestnuts to inform conservation and breeding with blight-resistant varieties. Population structure analyses suggest that the chestnut range can be roughly divided into northeast, central, and southwest populations. Within-population genomic diversity estimates revealed a clinal pattern with the highest diversity in the southwest, which likely reflects bottleneck events associated with Quaternary glaciation. Finally, we identified genomic regions under positive selection within each population, which suggests that defence against fungal pathogens is a common target of selection across all populations. Taken together, these results show that American chestnut underwent a postglacial expansion from the southern portion of its range leading to three extant genetic populations. These populations will serve as management units for breeding adaptive genetic variation into the blight-resistant tree populations for targeted reintroduction efforts.

Exome resequencing and GWAS for growth, ecophysiology, and chemical and metabolomic composition of wood of Populus trichocarpa.

BACKGROUND: Populus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation. RESULTS: We identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value < 6.1479E-8) in association with diameter, height, leaf carbon and nitrogen contents, and δ15N. The number of significant SNPs ranged from 2 to 220 per trait. Additionally, multiple-marker analyses by sliding-windows tests detected between 6 and 192 significant windows for the analyzed traits. The significant SNPs resided within genes that encode proteins belonging to different functional classes as such protein synthesis, energy/metabolism and DNA/RNA metabolism, among others. CONCLUSIONS: SNP-markers within genes associated with traits of importance for biomass production were detected. They contribute to characterize the genomic architecture of P. trichocarpa biomass required to support the development and application of marker breeding technologies.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding.

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

Evidence for extensive parallelism but divergent genomic architecture of adaptation along altitudinal and latitudinal gradients in Populus trichocarpa.

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.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii).

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

Leveraging whole-genome sequencing to estimate telomere length in plants.

Changes in telomere length are increasingly used to indicate species' response to environmental stress across diverse taxa. Despite this broad use, few studies have explored telomere length in plants. Thus, evaluation of new approaches for measuring telomeres in plants is needed. Rapid advances in sequencing approaches and bioinformatic tools now allow estimation of telomere content from whole-genome sequencing (WGS) data, a proxy for telomere length. While telomere content has been quantified extensively using quantitative polymerase chain reaction (qPCR) and WGS in humans, no study to date has compared the effectiveness of WGS in estimating telomere length in plants relative to qPCR approaches. In this study, we use 100 Populus clones re-sequenced using short-read Illumina sequencing to quantify telomere length comparing three different bioinformatic approaches (Computel, K-seek and TRIP) in addition to qPCR. Overall, telomere length estimates varied across different bioinformatic approaches, but were highly correlated across methods for individual genotypes. A positive correlation was observed between WGS estimates and qPCR, however, Computel estimates exhibited the greatest correlation. Computel incorporates genome coverage into telomere length calculations, suggesting that genome coverage is likely important to telomere length quantification when using WGS data. Overall, telomere estimates from WGS provided greater precision and accuracy of telomere length estimates relative to qPCR. The findings suggest WGS is a promising approach for assessing telomere length and, as the field of telomere ecology evolves, may provide added value to assaying response to biotic and abiotic environments for plants needed to accelerate plant breeding and conservation management.

Los cambios en la longitud de los telómeros se utilizan cada vez más para indicar la respuesta de las especies al estrés ambiental en diversos taxones. A pesar de este amplio uso, pocos estudios han explorado la longitud de los telómeros en las plantas. Por lo tanto, es necesario evaluar nuevos enfoques para medir los telómeros en las plantas. Los rápidos avances en los enfoques de secuenciación y las herramientas bioinformáticas ahora permiten estimar el contenido de los telómeros a partir de datos de secuenciación del genoma completo (WGS), un indicador de la longitud de los telómeros. Si bien el contenido de los telómeros se ha cuantificado ampliamente mediante la reacción en cadena de la polimerasa cuantitativa (qPCR) y WGS en humanos, ningún estudio hasta la fecha ha comparado la efectividad de WGS para estimar la longitud de los telómeros en plantas en relación con los enfoques de qPCR. En este estudio, utilizamos cien clones de álamos (Populus) resecuenciados mediante secuenciación Illumina de lectura corta para cuantificar la longitud de los telómeros comparando tres diferentes enfoques bioinformáticos, Computel, K-seek y TRIP, además de qPCR. En general, las estimaciones de la longitud de los telómeros variaron según los diferentes enfoques bioinformáticos, pero la longitud de los telómeros estuvo altamente correlacionada entre los métodos para genotipos individuales. Se observó una correlación positiva entre las estimaciones de WGS y qPCR; sin embargo, las estimaciones de Computel mostraron la mayor correlación. Computel incorpora la cobertura del genoma en los cálculos de la longitud de los telómeros, lo que sugiere que la cobertura del genoma probablemente es importante para la cuantificación de la longitud de los telómeros cuando se utilizan datos de WGS. En general, las estimaciones de los telómeros de WGS proporcionaron mayor precisión y exactitud de las estimaciones de la longitud de los telómeros en relación con la qPCR. Los hallazgos sugieren que WGS es un enfoque prometedor para evaluar la longitud de los telómeros y, a medida que evoluciona el campo de la ecología de los telómeros, puede proporcionar un valor agregado para analizar la respuesta a ambientes bióticos y abióticos de las plantas necesarias para acelerar los programas de mejoramiento genético y conservación.

Timber DNA release using focused ultrasound extraction (FUSE) for genetic species identification.

The use of genetic data for timber species and population assignment is a powerful tool for combating the illegal timber trade, but the challenges of extracting DNA from timber have prevented the routine use of genetics as a supply chain management tool. To overcome these challenges, we explored the feasibility of focused ultrasound extraction (FUSE) for rapid DNA release from timber. Using high-pressure ultrasound pulses, FUSE generates a cavitation bubble cloud that disintegrates samples into acellular debris, resulting in the mechanical release of DNA. In this work, FUSE was applied to white oak (Quercus alba) timber shavings to test the feasibility of using FUSE for timber DNA extraction for the first time. Results showed that FUSE processing disintegrated the tissue samples and released significant quantities of DNA. After five minutes of tissue processing DNA quantities of 0.21 ± 0.02 ng/mg, 0.99 ± 0.32 ng/mg, and 0.14 ± 0.01 ng/mg, were released from medium, coarse, and combination shaving groups, respectively. Amplification and sequencing of regions within the matK and rbcL chloroplast genes confirmed that the quality of DNA prepared with FUSE was suitable for PCR and short-read sequencing applications. Overall, these results show that FUSE can serve as a DNA sample preparation method capable of releasing high-quality DNA from timber in a fraction of the time required by conventional extraction methods. Based on the improved efficiency of DNA release with FUSE, ongoing work aims to develop this technology into portable systems that can be used to rapidly prepare timber samples for genetic species identification.

The genetic architecture of repeated local adaptation to climate in distantly related plants.

Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the 'cost of complexity' theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study.

Evidence for extensive heterotrophic metabolism, antioxidant action, and associated regulatory events during winter hardening in Sitka spruce.

BACKGROUND: Cold acclimation in woody perennials is a metabolically intensive process, but coincides with environmental conditions that are not conducive to the generation of energy through photosynthesis. While the negative effects of low temperatures on the photosynthetic apparatus during winter have been well studied, less is known about how this is reflected at the level of gene and metabolite expression, nor how the plant generates primary metabolites needed for adaptive processes during autumn. RESULTS: The MapMan tool revealed enrichment of the expression of genes related to mitochondrial function, antioxidant and associated regulatory activity, while changes in metabolite levels over the time course were consistent with the gene expression patterns observed. Genes related to thylakoid function were down-regulated as expected, with the exception of plastid targeted specific antioxidant gene products such as thylakoid-bound ascorbate peroxidase, components of the reactive oxygen species scavenging cycle, and the plastid terminal oxidase. In contrast, the conventional and alternative mitochondrial electron transport chains, the tricarboxylic acid cycle, and redox-associated proteins providing reactive oxygen species scavenging generated by electron transport chains functioning at low temperatures were all active. CONCLUSIONS: A regulatory mechanism linking thylakoid-bound ascorbate peroxidase action with "chloroplast dormancy" is proposed. Most importantly, the energy and substrates required for the substantial metabolic remodeling that is a hallmark of freezing acclimation could be provided by heterotrophic metabolism.

DNA release from plant tissue using focused ultrasound extraction (FUSE).

Premise: Sample preparation in genomics is a critical step that is often overlooked in molecular workflows and impacts the success of downstream genetic applications. This study explores the use of a recently developed focused ultrasound extraction (FUSE) technique to enable the rapid release of DNA from plant tissues for genetic analysis. Methods: FUSE generates a dense acoustic cavitation bubble cloud that pulverizes targeted tissue into acellular debris. This technique was applied to leaf samples of American chestnut (Castanea dentata), tulip poplar (Liriodendron tulipifera), red maple (Acer rubrum), and chestnut oak (Quercus montana). Results: We observed that FUSE can extract high quantities of DNA in 9-15 min, compared to the 30 min required for control DNA extraction methods. FUSE extracted DNA quantities of 24.33 ± 6.51 ng/mg and 35.32 ± 9.21 ng/mg from American chestnut and red maple, respectively, while control methods yielded 6.22 ± 0.87 ng/mg and 11.51 ± 1.95 ng/mg, respectively. The quality of the DNA released by FUSE allowed for successful amplification and next-generation sequencing. Discussion: These results indicate that FUSE can improve DNA extraction efficiency for leaf tissues. Continued development of this technology aims to adapt to field-deployable systems to increase the cataloging of genetic biodiversity, particularly in low-resource biodiversity hotspots.

Targeted enrichment of the black cottonwood (Populus trichocarpa) gene space using sequence capture.

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.

Divergent selection and heterogeneous migration rates across the range of Sitka spruce (Picea sitchensis).

Gene flow and effective population size (N(e)) should depend on a population's position within its range: those near the edges are expected to have smaller N(e) and lower relative emigration rates, whereas those nearer the centre should have larger N(e) 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 F(ST) 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.

Metabolic dynamics during autumn cold acclimation within and among populations of Sitka spruce (Picea sitchensis).

• Autumnal cold acclimation in conifers is a complex process, the timing and extent of which vary widely along latitudinal gradients for many tree species and reflect local adaptation to climate. Although previous studies have detailed some aspects of the metabolic remodelling that accompanies cold acclimation in conifers, little is known about global metabolic dynamics, or how these changes vary among phenotypically divergent populations. • Using untargeted GC-MS metabolite profiling, we monitored metabolic dynamics during autumnal cold acclimation in three populations of Sitka spruce from the southern, central, and northern portions of the species range, which differ in both the timing and extent of cold acclimation. • Latitudinal variation was evident in the nature, intensity, and timing of metabolic events. Early development of strong freezing tolerance in the northern population was associated with a transient accumulation of amino acids. By late autumn, metabolic profiles were highly similar between the northern and central populations, whereas profiles for the southern population were relatively distinct. • Our results provide insight into the metabolic architecture of latitudinal adaptive variation in autumn acclimation and show that different mechanisms are the basis of early October cold hardiness and autumn-acclimated cold hardiness.

Toward genomic selection in Pinus taeda: Integrating resources to support array design in a complex conifer genome.

PREMISE: An informatics approach was used for the construction of an Axiom genotyping array from heterogeneous, high-throughput sequence data to assess the complex genome of loblolly pine (Pinus taeda). METHODS: High-throughput sequence data, sourced from exome capture and whole genome reduced-representation approaches from 2698 trees across five sequence populations, were analyzed with the improved genome assembly and annotation for the loblolly pine. A variant detection, filtering, and probe design pipeline was developed to detect true variants across and within populations. From 8.27 million variants, a total of 642,275 were evaluated and 423,695 of those were screened across a range-wide population. RESULTS: The final informatics and screening approach delivered an Axiom array representing 46,439 high-confidence variants to the forest tree breeding and genetics community. Based on the annotated reference genome, 34% were located in or directly upstream or downstream of genic regions. DISCUSSION: The Pita50K array represents a genome-wide resource developed from sequence data for an economically important conifer, loblolly pine. It uniquely integrates independent projects that assessed trees sampled across the native range. The challenges associated with the large and repetitive genome are addressed in the development of this resource.

Widespread, ecologically relevant genetic markers developed from association mapping of climate-related traits in Sitka spruce (Picea sitchensis).

• 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.

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

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

Optimizing genomic selection for blight resistance in American chestnut backcross populations: A trade-off with American chestnut ancestry implies resistance is polygenic.

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.