Forests are not only the most predominant of the Earth's terrestrial ecosystems, but are also the core supply for essential products for human use. However, global climate change and ongoing population explosion severely threatens the health of the forest ecosystem and aggravtes the deforestation and forest degradation. Forest genomics has great potential of increasing forest productivity and adaptation to the changing climate. In the last two decades, the field of forest genomics has advanced quickly owing to the advent of multiple high-throughput sequencing technologies, single cell RNA-seq, clustered regularly interspaced short palindromic repeats (CRISPR)-mediated genome editing, and spatial transcriptomes, as well as bioinformatics analysis technologies, which have led to the generation of multidimensional, multilayered, and spatiotemporal gene expression data. These technologies, together with basic technologies routinely used in plant biotechnology, enable us to tackle many important or unique issues in forest biology, and provide a panoramic view and an integrative elucidation of molecular regulatory mechanisms underlying phenotypic changes and variations. In this review, we recapitulated the advancement and current status of 12 research branches of forest genomics, and then provided future research directions and focuses for each area. Evidently, a shift from simple biotechnology-based research to advanced and integrative genomics research, and a setup for investigation and interpretation of many spatiotemporal development and differentiation issues in forest genomics have just begun to emerge.

Summary: Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome‐level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome‐wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic‐based breeding of important traits in perennial woody plants.

Drought stress limits woody species productivity and influences tree distribution. However, dissecting the molecular mechanisms that underpin drought responses in forest trees can be challenging due to trait complexity. Here, using a panel of 300 Chinese white poplar (Populus tomentosa) accessions collected from different geographical climatic regions in China, we performed a genome-wide association study (GWAS) on seven drought-related traits and identified PtoWRKY68 as a candidate gene involved in the response to drought stress. A 12-bp insertion and/or deletion and three nonsynonymous variants in the PtoWRKY68 coding sequence categorized natural populations of P. tomentosa into two haplotype groups, PtoWRKY68^hap1 and PtoWRKY68^hap2. The allelic variation in these two PtoWRKY68 haplotypes conferred differential transcriptional regulatory activities and binding to the promoters of downstream abscisic acid (ABA) efflux and signaling genes. Overexpression of PtoWRKY68^hap1 and PtoWRKY68^hap2 in Arabidopsis (Arabidopsis thaliana) ameliorated the drought tolerance of two transgenic lines and increased ABA content by 42.7% and 14.3% compared to wild-type plants, respectively. Notably, PtoWRKY68^hap1 (associated with drought tolerance) is ubiquitous in accessions in water-deficient environments, whereas the drought-sensitive allele PtoWRKY68^hap2 is widely distributed in well-watered regions, consistent with the trends in local precipitation, suggesting that these alleles correspond to geographical adaptation in Populus. Moreover, quantitative trait loci analysis and an electrophoretic mobility shift assay showed that SHORT VEGETATIVE PHASE (PtoSVP.3) positively regulates the expression of PtoWRKY68 under drought stress. We propose a drought tolerance regulatory module in which PtoWRKY68 modulates ABA signaling and accumulation, providing insight into the genetic basis of drought tolerance in trees. Our findings will facilitate molecular breeding to improve the drought tolerance of forest trees.

Summary: Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein‐encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue‐specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis‐related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker‐assisted breeding efforts to genetically improve woody plants.