NAC transcription factors (TFs) participate in multiple biological processes, including biotic and abiotic stress responses, signal transduction and development. Cold stress can adversely impact plant growth and development, thereby limiting agricultural productivity. Prunus mume, an excellent horticultural crop, is widely cultivated in Asian countries. Its flower can tolerate freezing-stress in the early spring. To investigate the putative NAC genes responsible for cold-stress, we identified and analyzed 113 high-confidence PmNAC genes and characterized them by bioinformatics tools and expression profiles. These PmNACs were clustered into 14 sub-families and distributed on eight chromosomes and scaffolds, with the highest number located on chromosome 3. Duplicated events resulted in a large gene family; 15 and 8 pairs of PmNACs were the result of tandem and segmental duplicates, respectively. Moreover, three membrane-bound proteins (PmNAC59/66/73) and three miRNA-targeted genes (PmNAC40/41/83) were identified. Most PmNAC genes presented tissue-specific and time-specific expression patterns. Sixteen PmNACs (PmNAC11/19/20/23/41/48/58/74/75/76/78/79/85/86/103/111) exhibited down-regulation during flower bud opening and are, therefore, putative candidates for dormancy and cold-tolerance. Seventeen genes (PmNAC11/12/17/21/29/42/30/48/59/66/73/75/85/86/93/99/111) were highly expressed in stem during winter and are putative candidates for freezing resistance. The cold-stress response pattern of 15 putative PmNACs was observed under 4 °C at different treatment times. The expression of 10 genes (PmNAC11/20/23/40/42/48/57/60/66/86) was upregulated, while 5 genes (PmNAC59/61/82/85/107) were significantly inhibited. The putative candidates, thus identified, have the potential for breeding the cold-tolerant horticultural plants. This study increases our understanding of functions of the NAC gene family in cold tolerance, thereby potentially intensifying the molecular breeding programs of woody plants.

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean N_A = 5.481, H_E = 0.548) and moderate population differentiation (pairwise F_st = 0.048–0.078, N_m = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.

Silk production is a classic example of a domestication trait, yet the cell-type-specific driver of its enhancement in the silkworm Bombyx mori remains unresolved. To address this, we integrated extensive bulk RNA-seq data with a single-nucleus RNA-seq atlas of silk glands (SGs) from domestic B. mori and wild B. mandarina for deconvolution analysis. This identified phenotype-associated cell subpopulations (Scissor+ and Scissor− cells) that enrich in B. mori and B. mandarina, respectively. Transcriptomic characterization revealed that B. mori SG cells exhibit a pervasive “pro-synthesis” transcriptional state, with concerted upregulation of silk protein genes and metabolic pathways. Conversely, B. mandarina cells maintained a “protective–adaptive” state, enriched for stress response and xenobiotic metabolism genes. Pseudotime analysis further delineated the cell state transitions, pinpointing key dynamic gene expression linked to high silk yield. Our findings demonstrate that domestication reshaped the silk gland cellular landscape, promoting a systemic shift toward a synthesis-optimized cell state. This study offers a new framework at the cellular level to elucidate the evolution of complex traits under selection.