Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

A transcriptome atlas of zygotic and somatic embryogenesis in Norway spruce

Somatic embryogenesis (SE) is a powerful model system for studying embryo development and an important method for scaling up availability of elite and climate-adapted genetic material of Norway spruce (Picea abies L. Karst). However, there are several steps during the development of the somatic embryo (Sem) that are suboptimal compared to zygotic embryo (Zem) development. These differences are poorly understood and result in substantial yield losses during plant production, which limits cost-effective large-scale production of SE plants. This study presents a comprehensive data resource profiling gene expression during zygotic and somatic embryo development to support studies aiming to advance understanding of gene regulatory programmes controlling embryo development. Transcriptome expression patterns were analysed during zygotic embryogenesis (ZE) in Norway spruce, including separated samples of the female gametophytes and Zem, and at multiple stages during SE. Expression data from eight developmental stages of SE, starting with pro-embryogenic masses (PEMs) up until germination, revealed extensive modulation of the transcriptome between the early and mid-stage maturing embryos and at the transition of desiccated embryos to germination. Comparative analysis of gene expression changes during ZE and SE identified differences in the pattern of gene expression changes and functional enrichment of these provided insight into the associated biological processes. Orthologs of transcription factors known to regulate embryo development in angiosperms were differentially regulated during Zem and Sem development and in the different zygotic embryo tissues, providing clues to the differences in development observed between Zem and Sem. This resource represents the most comprehensive dataset available for exploring embryo development in conifers.

The evolutionary adaptation of wood-decay macrofungi to host gymnosperms differs from that to host angiosperms

Wood-decay macrofungi play a vital role in forest ecosystems by promoting nutrient cycling and soil structure, and their evolution is closely related to their host plants. This study investigates the potential evolutionary adaptation of wood-decay macrofungi to their host plants, focusing on whether these relationships differ between gymnosperms and angiosperms. While previous research has suggested non-random associations between specific fungi and plant deadwood, direct evidence of evolutionary adaptation has been lacking. Our study, conducted in a subtropical region, utilized metabarcoding techniques to identify deadwood species and associated fungi. We found significant evidence of evolutionary adaptation when considering all sampled species collectively. However, distinct patterns emerged when comparing angiosperms and gymnosperms: a significant evolutionary adaptation was observed of wood-decay macrofungi to angiosperms, but not to gymnosperms. This variation may be due to the longer evolutionary history and more stable species interactions of gymnosperms, as indicated by a higher modularity coefficient (r = .452), suggesting greater specialization. In contrast, angiosperms, being evolutionarily younger, displayed less stable and more coevolving interactions with fungi, reflected in a lower modularity coefficient (r = .387). Our findings provide the first direct evidence of differential evolutionary adaptation dynamics of these fungi to angiosperms versus gymnosperms, enhancing our understanding of forest ecosystem carbon cycling and resource management.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

The impact of insect egg deposition on Pinus sylvestris transcriptomic and phytohormonal responses to larval herbivory

Plants can improve their resistance to feeding damage by insects if they have perceived insect egg deposition prior to larval feeding. Molecular analyses of these egg-mediated defence mechanisms have until now focused on angiosperm species. It is unknown how the transcriptome of a gymnosperm species responds to insect eggs and subsequent larval feeding. Scots pine (Pinus sylvestris L.) is known to improve its defences against larvae of the herbivorous sawfly Diprion pini L. if it has previously received sawfly eggs. Here, we analysed the transcriptomic and phytohormonal responses of Scots pine needles to D. pini eggs (E-pine), larval feeding (F-pine) and to both eggs and larval feeding (EF-pine). Pine showed strong transcriptomic responses to sawfly eggs and-as expected-to larval feeding. Many egg-responsive genes were also differentially expressed in response to feeding damage, and these genes play an important role in biological processes related to cell wall modification, cell death and jasmonic acid signalling. EF-pine showed fewer transcriptomic changes than F-pine, whereas EF-treated angiosperm species studied so far showed more transcriptional changes to the initial phase of larval feeding than only feeding-damaged F-angiosperms. However, as with responses of EF-angiosperms, EF-pine showed higher salicylic acid concentrations than F-pine. Based on the considerable overlap of the transcriptomes of E- and F-pine, we suggest that the weaker transcriptomic response of EF-pine than F-pine to larval feeding damage is compensated by the strong, egg-induced response, which might result in maintained pine defences against larval feeding.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Effects of drought stress memory on the accumulation of stress-protective compounds in naturally grown pine and spruce

Forest trees are subjected to multiple stressors during their long lifetime and therefore require effective and finely regulated stress-protective systems. Stressors can induce protective systems either directly or with the involvement of stress memory mechanisms. Stress memory has only begun to be uncovered in model plants and is unexplored in coniferous species. Therefore, we studied the possible role of stress memory in the regulation of the accumulation of stress-protective compounds (heat shock proteins, dehydrins, proline) in the needles of naturally grown Scots pine and Norway spruce trees subjected to the subsequent action of long-term (multiyear) and short-term (seasonal) water shortages. Although the water deficit was relatively mild, it significantly influenced the pattern of expression of stress memory-related heat shock factor (HSF) and SWI/SNF genes, indicating the formation of stress memory in both species. In spruce, dehydrin accumulation was increased by water shortage in a manner compatible with Type II stress memory. The accumulation of HSP40 in spruce needles was positively influenced by long-term water shortage, but this increase was unlikely to be of biological importance due to the concomitant decrease in HSP70, HSP90 and HSP101 accumulation. Finally, proline accumulation was negatively influenced by short-term water deficit in spruce. In pine, no one protective compound accumulated in response to water stress. Taken together, the results indicate that the accumulation of stress-protective compounds was generally independent of stress memory effects both in pine and in spruce.

Comparative transcriptome analysis of juniper branches infected by Gymnosporangium spp. highlights their different infection strategies associated with cytokinins

BACKGROUND

Gymnosporangium asiaticum and G. yamadae can share Juniperus chinensis as the telial host, but the symptoms are completely different. The infection of G. yamadae causes the enlargement of the phloem and cortex of young branches as a gall, but not for G. asiaticum, suggesting that different molecular interaction mechanisms exist the two Gymnosporangium species with junipers.

RESULTS

Comparative transcriptome analysis was performed to investigate genes regulation of juniper in responses to the infections of G. asiaticum and G. yamadae at different stages. Functional enrichment analysis showed that genes related to transport, catabolism and transcription pathways were up-regulated, while genes related to energy metabolism and photosynthesis were down-regulated in juniper branch tissues after infection with G. asiaticum and G. yamadae. The transcript profiling of G. yamadae-induced gall tissues revealed that more genes involved in photosynthesis, sugar metabolism, plant hormones and defense-related pathways were up-regulated in the vigorous development stage of gall compared to the initial stage, and were eventually repressed overall. Furthermore, the concentration of cytokinins (CKs) in the galls tissue and the telia of G. yamadae was significantly higher than in healthy branch tissues of juniper. As well, tRNA-isopentenyltransferase (tRNA-IPT) was identified in G. yamadae with highly expression levels during the gall development stages.

CONCLUSIONS

In general, our study provided new insights into the host-specific mechanisms by which G. asiaticum and G. yamadae differentially utilize CKs and specific adaptations on juniper during their co-evolution.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C₁₈ Δ⁹-elongase and a C₂₀ Δ⁵-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ⁵-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Building resiliency in conifer forests: Interior spruce crosses among weevil resistant and susceptible parents produce hybrids appropriate for multi-trait selection

Tree planting programs now need to consider climate change increasingly, therefore, the resistance to pests plays an essential role in enabling tree adaptation to new ranges through tree population movement. The weevil Pissodes strobi (Peck) is a major pest of spruces and substantially reduces lumber quality. We revisited a large Interior spruce provenance/progeny trial (2,964 genotypes, 42 families) of varying susceptibility, established in British Columbia. We employed multivariate mixed linear models to estimate covariances between, and genetic control of, juvenile height growth and resistance traits. We performed linear regressions and ordinal logistic regressions to test for impact of parental origin on growth and susceptibility to the pest, respectively. A significant environmental component affected the correlations between resistance and height, with outcomes dependent on families. Parents sourced from above 950 m a.s.l. elevation negatively influenced host resistance to attacks, probably due to higher P. engelmannii proportion. For the genetic contribution of parents sourced from above 1,200 m a.s.l., however, we found less attack severity, probably due to a marked mismatch in phenologies. This clearly highlights that interspecific hybrid status might be a good predictor for weevil attacks and delineates the boundaries of successful spruce population movement. Families resulting from crossing susceptible parents generally showed fast-growing trees were the most affected by weevil attacks. Such results indicate that interspecific 'hybrids' with a higher P. glauca ancestry might be genetically better equipped with an optimized resource allocation between defence and growth and might provide the solution for concurrent improvement in resistance against weevil attacks, whilst maintaining tree productivity.

Integrative analysis in Pinus revealed long-term heat stress splicing memory

Due to the current climate change, many studies have described main drivers in abiotic stress. Recent findings suggest that alternative splicing (AS) has a critical role in controlling plant responses to high temperature. AS is a mechanism that allows organisms to create an assortment of RNA transcripts and proteins using a single gene. However, the most important roles of AS in stress could not be rigorously addressed because research has been focused on model species, covering only a narrow phylogenetic and lifecycle spectrum. Thus, AS degree of diversification among more dissimilar taxa in heat response is still largely unknown. To fill this gap, the present study employs a systems biology approach to examine how the AS landscape responds to and 'remembers' heat stress in conifers, a group which has received little attention even though their position can solve key evolutionary questions. Contrary to angiosperms, we found that potential intron retention may not be the most prevalent type of AS. Furthermore, our integrative analysis with metabolome and proteome data places splicing as the main source of variation during the response. Finally, we evaluated possible acquired long-term splicing memory in a diverse subset of events, and although this mechanism seems to be conserved in seed plants, AS dynamics are divergent. These discoveries reveal the particular way of remembering past temperature changes in long-lived plants and open the door to include species with unique features to determine the extent of conservation in gene expression regulation.

Conifer Biotechnology: An Overview

The peculiar characteristics of conifers determine the difficulty of their study and their great importance from various points of view. However, their study faces numerous important scientific, methodological, cultural, economic, social, and legal challenges. This paper presents an approach to several of those challenges and proposes a multidisciplinary scientific perspective that leads to a holistic understanding of conifers from the perspective of the latest technical, computer, and scientific advances. This review highlights the deep connection that all scientific contributions to conifers can have in each other as fully interrelated communicating vessels.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome. Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.

Phylogenetic Occurrence of the Phenylpropanoid Pathway and Lignin Biosynthesis in Plants

The phenylpropanoid pathway serves as a rich source of metabolites in plants and provides precursors for lignin biosynthesis. Lignin first appeared in tracheophytes and has been hypothesized to have played pivotal roles in land plant colonization. In this review, we summarize recent progress in defining the lignin biosynthetic pathway in lycophytes, monilophytes, gymnosperms, and angiosperms. In particular, we review the key structural genes involved in p-hydroxyphenyl-, guaiacyl-, and syringyl-lignin biosynthesis across plant taxa and consider and integrate new insights on major transcription factors, such as NACs and MYBs. We also review insight regarding a new transcriptional regulator, 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase, canonically identified as a key enzyme in the shikimate pathway. We use several case studies, including EPSP synthase, to illustrate the evolution processes of gene duplication and neo-functionalization in lignin biosynthesis. This review provides new insights into the genetic engineering of the lignin biosynthetic pathway to overcome biomass recalcitrance in bioenergy crops.

Quantitative analysis of differential dehydrin regulation in pine and spruce seedlings under water deficit

Dehydrins are well-known components of plant responses to different stresses that cause dehydration, including drought, freezing, salinity, etc. In conifers, the dehydrin gene family is very large, implying that the members of this family have important physiological functions in conifer stress tolerance. However, dehydrin gene expression displays a wide range of responses to stress, from thousand-fold increased expression to decreased expression, and it is generally unknown how regulatory systems are connected at the mRNA and protein levels. Therefore, we studied these aspects of dehydrin regulation in Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) H. Karst) seedlings under polyethylene glycol 6000-induced osmotic stress ranging from relatively low (culture medium water potential of -0.15 MPa) to very high (-1.0 MPa) intensities. In pine, the major dehydrin protein was Dhn1 in both the roots and needles, and in spruce, two isoforms of the Dhn4 protein were the major dehydrins; both of these proteins are AESK-type dehydrins. The genes encoding these major proteins were highly expressed even under control conditions; surprisingly, we also observed several highly expressed dehydrin genes that were not abundantly translated. Under osmotic stress, the most prominent expression changes were observed for the dehydrin genes with low basal expression levels, whereas highly expressed genes generally demonstrated rather modest changes in expression. We report proposed constitutive physiological functions of the AESK-type dehydrins in Pinaceae plants.

Clair JB, Crepeau MW, Salzberg SL, Langley CH, Allen B, Neale DB. Dissecting the Polygenic Basis of Cold Adaptation Using Genome-Wide Association of Traits and Environmental Data in Douglas-fir

Understanding the genomic and environmental basis of cold adaptation is key to understand how plants survive and adapt to different environmental conditions across their natural range. Univariate and multivariate genome-wide association (GWAS) and genotype-environment association (GEA) analyses were used to test associations among genome-wide SNPs obtained from whole-genome resequencing, measures of growth, phenology, emergence, cold hardiness, and range-wide environmental variation in coastal Douglas-fir (Pseudotsuga menziesii). Results suggest a complex genomic architecture of cold adaptation, in which traits are either highly polygenic or controlled by both large and small effect genes. Newly discovered associations for cold adaptation in Douglas-fir included 130 genes involved in many important biological functions such as primary and secondary metabolism, growth and reproductive development, transcription regulation, stress and signaling, and DNA processes. These genes were related to growth, phenology and cold hardiness and strongly depend on variation in environmental variables such degree days below 0c, precipitation, elevation and distance from the coast. This study is a step forward in our understanding of the complex interconnection between environment and genomics and their role in cold-associated trait variation in boreal tree species, providing a baseline for the species' predictions under climate change.

Efficient Multi-Sites Genome Editing and Plant Regeneration via Somatic Embryogenesis in Picea glauca

Conifers are the world's major source of timber and pulpwood and have great economic and ecological value. Currently, little research on the application of CRISPR/Cas9, the commonly used genome-editing tool in angiosperms, has been reported in coniferous species. An efficient CRISPR/Cas9 system based on somatic embryogenesis (SEis) suitable for conifers could benefit both fundamental and applied research in these species. In this study, the SpCas9 gene was optimized based on codon bias in white spruce, and a spruce U6 promoter was cloned and function-validated for use in a conifer specific CRISPR/Cas9 toolbox, i.e., PgCas9/PaU6. With this toolbox, a genome-editing vector was constructed to target the DXS1 gene of white spruce. By Agrobacterium-mediated transformation, the genome-editing vector was then transferred into embryogenic tissue of white spruce. Three resistant embryogenic tissues were obtained and used for regenerating plants via SEis. Albino somatic embryo (SE) plants with mutations in DXS1 were obtained in all of the three events, and the ratios of the homozygous and biallelic mutants in the 18 albino mutants detected were 22.2% in both cases. Green plants with mutations in DXS1 were also produced, and the ratios of the DXS1 mutants to the total green plants were 7.9, 28, and 13.5%, respectively, among the three events. Since 22.7% of the total 44 mutants were edited at both of the target sites 1 and 2, the CRISPR/Cas9 toolbox in this research could be used for multi-sites genome editing. More than 2,000 SE plants were regenerated in vitro after genome editing, and part of them showed differences in plant development. Both chimerism and mosaicism were found in the SE plants of white spruce after genome editing with the CRISPR/Cas9 toolbox. The conifer-specific CRISPR/Cas9 system developed in this research could be valuable in gene function research and trait improvement.

Genomic basis of white pine blister rust quantitative disease resistance and its relationship with qualitative resistance

The genomic architecture and molecular mechanisms controlling variation in quantitative disease resistance loci are not well understood in plant species and have been barely studied in long-generation trees. Quantitative trait loci mapping and genome-wide association studies were combined to test a large single nucleotide polymorphism (SNP) set for association with quantitative and qualitative white pine blister rust resistance in sugar pine. In the absence of a chromosome-scale reference genome, a high-density consensus linkage map was generated to obtain locations for associated SNPs. Newly discovered associations for white pine blister rust quantitative disease resistance included 453 SNPs involved in wide biological functions, including genes associated with disease resistance and others involved in morphological and developmental processes. In addition, NBS-LRR pathogen recognition genes were found to be involved in quantitative disease resistance, suggesting these newly reported genes are qualitative genes with partial resistance, they are the result of defeated qualitative resistance due to avirulent races, or they have epistatic effects on qualitative disease resistance genes. This study is a step forward in our understanding of the complex genomic architecture of quantitative disease resistance in long-generation trees, and constitutes the first step towards marker-assisted disease resistance breeding in white pine species.

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.