The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Picea koraiensis is major silvicultural and timber species in northeast China, and its distribution area is an important transition zone for genus spruce migration. The degree of intraspecific differentiation of P. koraiensis is high, but population structure and differentiation mechanisms are not clear. In this study, 523,761 single nucleotide polymorphisms (SNPs) were identified in 113 individuals from 9 populations of P. koraiensis by genotyping-by-sequencing (GBS). Population genomic analysis showed that P. koraiensis was divided into three geoclimatic regions: Great Khingan Mountains climatic region, Lesser Khingan Mountains climatic region, and Changbai Mountain climatic region. Mengkeshan (MKS) population on the northern edge of the distribution area and Wuyiling (WYL) population located in the mining area are two highly differentiated groups. Selective sweep analysis showed that MKS and WYL populations had 645 and 1126 selected genes, respectively. Genes selected in the MKS population were associated with flowering and photomorphogenesis, cellular response to water deficit, and glycerophospholipid metabolism; genes selected in the WYL population were associated with metal ion transport, biosynthesis of macromolecules, and DNA repair. Climatic factors and heavy metal stress drives divergence in MKS and WYL populations, respectively. Our findings provide insights into adaptive divergence mechanisms in Picea and will contribute to molecular breeding studies.

To evaluate the effects of donor ages on growth and stress resistance of 6-year-old seedlings propagated from 5-, 2000-, and 3000-year-old Platycladus orientalis donors with grafting, cutting, and seed sowing, growth indicators and physiological and transcriptomic analyses were performed in 6-year-old seedlings in winter. Results showed that basal stem diameters and plant heights of seedlings of the three propagation methods decreased with the age of the donors, and the sown seedlings were the thickest and tallest. The contents of soluble sugar, chlorophyll, and free fatty acid in apical leaves of the three propagation methods were negatively correlated with donor ages in winter, while the opposite was true for flavonoid and total phenolic. The contents of flavonoid, total phenolic, and free fatty acid in cutting seedlings were highest in the seedlings propagated in the three methods in winter. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of differentially expressed genes showed phenylpropanoid biosynthesis and fatty acid metabolism pathways, and their expression levels were up-regulated in apical leaves from 6-year-old seedlings propagated from 3000-year-old P. orientalis donors. In addition, hub genes analysis presented that C4H, OMT1, CCR2, PAL, PRX52, ACP1, AtPDAT2, and FAD3 were up-regulated in cutting seedlings, and the gene expression levels decreased in seedlings propagated from 2000- and 3000-year-old donors. These findings demonstrate the resistance stability of cuttings of P. orientalis and provide insights into the regulatory mechanisms of seedlings of P. orientalis propagated from donors at different ages in different propagation methods against low-temperature stress.

In woody plants, bark is an important protective tissue which can participate in photosynthesis, manage water loss, and transport assimilates. Studying the bark anatomical traits can provide insight into plant environmental adaptation strategies. However, a systematic understanding of the variability in bark anatomical traits and their drivers is lacking in woody plants. In this study, the bark anatomical traits of 23 Picea species were determined in a common garden experiment. We analyzed interspecific differences and interpreted the patterns in bark anatomical traits in relation to phylogenetic relationships and climatic factors of each species according to its global distribution. The results showed that there were interspecific differences in bark anatomical traits of Picea species. Phloem thickness was positively correlated with parenchyma cell size, possibly related to the roles of parenchyma cells in the radial transport of assimilates. Sieve cell size was negatively correlated with the radial diameter of resin ducts, and differences in sieve cells were possibly related to the formation and expansion of resin ducts. There were no significant phylogenetic signals for any bark anatomical trait, except the tangential diameter of resin ducts. Phloem thickness and parenchyma cell size were affected by temperature-related factors of their native range, while sieve cell size was influenced by precipitation-related factors. Bark anatomical traits were not significantly different under wet and dry climates. This study makes an important contribution to our understanding of variability in bark anatomical traits among Picea species and their ecological adaptations.

The effects of tree age on the growth of cutting seedlings propagated from ancient trees have been an important issue in plant breeding and cultivation. In order to understand seedling growth and stress resistance stability, phenotypic measurements, physiological assays, and high-throughput transcriptome sequencing were performed on sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old Platycladus orientalis donors. In this study, the growth of cutting seedlings propagated from ancient trees was significantly slower; the soluble sugar and chlorophyll contents gradually decreased with the increase in the age of donors, and the flavonoid and total phenolic contents of sown seedlings were higher than those of cutting seedlings. Enrichment analysis of differential genes showed that plant hormone signal transduction, the plant–pathogen interaction, and the flavone and flavonol biosynthesis pathways were significantly up-regulated with the increasing age of cutting seedlings propagated from 300- and 700-year-old donors. A total of 104,764 differentially expressed genes were calculated using weighted gene co-expression network analysis, and 8 gene modules were obtained. Further, 10 hub genes in the blue module were identified, which revealed that the expression levels of JAZ, FLS, RPM1/RPS3, CML, and RPS2 increased with the increase in tree age. The results demonstrated that the age of the donors seriously affected the growth of P. orientalis cutting seedlings and that cutting propagation can preserve the resistance of ancient trees. The results of this study provide important insights into the effects of age on asexually propagated seedlings, reveal potential molecular mechanisms, and contribute to an improvement in the level of breeding and conservation of ancient germplasm resources of P. orientalis trees.