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.

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla. The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana. In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Background: Pinus koraiensis is an evergreen tree species with strong cold resistance. However, the transcriptomic patterns in response to cold stress are poorly understood for P. koraiensis. In this study, global transcriptome profiles were generated for P. koraiensis under cold stress (− 20 °C) over time by high-throughput sequencing. Results: More than 763 million clean reads were produced, which assembled into a nonredundant data set of 123,445 unigenes. Among them, 38,905 unigenes had homology with known genes, 18,239 were assigned to 54 gene ontology (GO) categories and 18,909 were assigned to 25 clusters of orthologous groups (COG) categories. Comparison of transcriptomes of P. koraiensis seedlings grown at room temperature (20 °C) and low temperature (− 20 °C) revealed 9842 differential expressed genes (DEGs) in the 6 h sample, 9250 in the 24 h sample, and 9697 in the 48 h sample. The number of DEGs in the pairwise comparisons of 6 h, 24 h and 48 h was relatively small. The accuracy of the RNA-seq was validated by analyzing the expression patterns of 12 DEGs by quantitative real-time PCR (qRT-PCR). In this study, 34 DEGs (22 upregulated and 12 downregulated) were involved in the perception and transmission of cold signals, 96 DEGs (41 upregulated and 55 downregulated) encoding 8 transcription factors that regulated cold-related genes expression, and 27 DEGs (17 upregulated and 10 downregulated) were involved in antioxidant mechanisms in response to cold stress. Among them, the expression levels of c63631_g1 (annexin D1), c65620_g1 (alpha-amylase isozyme 3C), c61970_g1 (calcium-binding protein KIC), c51736_g1 (ABA), c58408_g1 (DREB3), c66599_g1 (DREB3), c67548_g2 (SOD), c55044_g1 (CAT), c71938_g2 (CAT) and c11358_g1 (GPX) first increased significantly and then decreased significantly with the extension of stress time. Conclusions: A large number of DEGs were identified in P. koraiensis under cold stress, especially the DEGs involved in the perception and transmission of cold signals, the DEGs encoding TFs related to cold regulation and the DEGs removing ROS in antioxidation mechanisms. The transcriptome and digital expression profiling of P. koraiensis could facilitate the understanding of the molecular control mechanism related to cold responses and provide the basis for the molecular breeding of conifers.

Background: Catalpa bungei is an important tree species used for timber in China and widely cultivated for economic and ornamental purposes. A high-density linkage map of C. bungei would be an efficient tool not only for identifying key quantitative trait loci (QTLs) that affect important traits, such as plant growth and leaf traits, but also for other genetic studies. Results: Restriction site-associated DNA sequencing (RAD-seq) was used to identify molecular markers and construct a genetic map. Approximately 280.77 Gb of clean data were obtained after sequencing, and in total, 25,614,295 single nucleotide polymorphisms (SNPs) and 2,871,647 insertions-deletions (InDels) were initially identified in the genomes of 200 individuals of a C. bungei (7080) × Catalpa duclouxii (16-PJ-3) F₁ population and their parents. Finally, 9072 SNP and 521 InDel markers that satisfied the requirements for constructing a genetic map were obtained. The integrated genetic map contained 9593 pleomorphic markers in 20 linkage groups and spanned 3151.63 cM, with an average distance between adjacent markers of 0.32 cM. Twenty QTLs for seven leaf traits and 13 QTLs for plant height at five successive time points were identified using our genetic map by inclusive composite interval mapping (ICIM). Q16–60 was identified as a QTL for five leaf traits, and three significant QTLs (Q9–1, Q18–66 and Q18–73) associated with plant growth were detected at least twice. Genome annotation suggested that a cyclin gene participates in leaf trait development, while the growth of C. bungei may be influenced by CDC48C and genes associated with phytohormone synthesis. Conclusions: This is the first genetic map constructed in C. bungei and will be a useful tool for further genetic study, molecular marker-assisted breeding and genome assembly.

The C-Repeat Binding Factor (CBF) gene family has been identified and characterized in multiple plant species, and it plays a crucial role in responding to low temperatures. Presently, only a few studies on tree species demonstrate the mechanisms and potential functions of CBFs associated with cold resistance, while our study is a novel report on the multi-aspect differences of CBFs among three tree species, compared to previous studies. In this study, genome-wide identification and analysis of the CBF gene family in Acer truncatum, Acer pseudosieboldianum, and Acer yangbiense were performed. The results revealed that 16 CBF genes (five ApseCBFs, four AcyanCBFs, and seven AtruCBFs) were unevenly distributed across the chromosomes, and most CBF genes were mapped on chromosome 2 (Chr2) and chromosome 11 (Chr11). The analysis of phylogenetic relationships, gene structure, and conserved motif showed that 16 CBF genes could be clustered into three subgroups; they all contained Motif 1 and Motif 5, and most of them only spanned one exon. The cis-acting elements analysis showed that some CBF genes might be involved in hormone and abiotic stress responsiveness. In addition, CBF genes exhibited tissue expression specificity. High expressions of ApseCBF1, ApseCBF3, AtruCBF1, AtruCBF4, AtruCBF6, AtruCBF7, and ApseCBF3, ApseCBF4, ApseCBF5 were detected on exposure to low temperature for 3 h and 24 h. Low expressions of AtruCBF2, AtruCBF6, AtruCBF7 were detected under cold stress for 24 h, and AtruCBF3 and AtruCBF5 were always down-regulated under cold conditions. Taken together, comprehensive analysis will enhance our understanding of the potential functions of the CBF genes on cold resistance, thereby providing a reference for the introduction of Acer species in our country.

Populus euphratica is mainly distributed in desert environments with dry and hot climate in summer and cold in winter. Compared with other poplars, P. euphratica is more resistant to salt stress. It is critical to investigate the transcriptome and molecular basis of salt tolerance in order to uncover stress-related genes. In this study, salt-tolerant treatment of P. euphratica resulted in an increase in osmo-regulatory substances and recovery of antioxidant enzymes. To improve the mining efficiency of candidate genes, the analysis combining both the transcriptome WGCNA and the former GWAS results was selected, and a range of key regulatory factors with salt resistance were found. The PeERF1 gene was highly connected in the turquoise modules with significant differences in salt stress traits, and the expression levels were significantly different in each treatment. For further functional verification of PeERF1, we obtained stable overexpression and dominant suppression transgenic lines by transforming into Populus alba × Populusglandulosa. The growth and physiological characteristics of the PeERF1 overexpressed plants were better than that of the wild type under salt stress. Transcriptome analysis of leaves of transgenic lines and WT revealed that highly enriched GO terms in DEGs were associated with stress responses, including abiotic stimuli responses, chemical responses, and oxidative stress responses. The result is helpful for in-depth analysis of the salt tolerance mechanism of poplar. This work provides important genes for poplar breeding with salt tolerance.