Several studies have focused on cold tolerance in multiple regulated levels. However, a genome-scale molecular analysis of the regulated network under the control of transcription factors (TFs) is still lacking, especially for trees. To comprehensively identify the TFs that regulate cold stress response in the paper mulberry and understand their regulatory interactions, transcriptomic data was used to assess changes in gene expression induced by exposure to cold. Results indicated that 794 TFs, belonging to 47 families and comprising more than 59% of the total TFs of this plant, were involved in the cold stress response. They were clustered into three groups, namely early, intermediate and late responsive groups which contained 95, 550 and 149 TFs, respectively. Among of these differentially expressed TFs, one bHLH, two ERFs and three CAMTAs were considered to be the key TFs functioning in the primary signal transduction. After that, at the intermediate stage of cold stress, there were mainly two biological processes that were regulated by TFs, namely cold stress resistance (including 5 bHLH, 14 ERFs, one HSF, 4 MYBs, 3 NACs, 11 WRKYs and so on) and growth and development of lateral organ or apical meristem (including ARR-B, B3, 5 bHLHs, 2 C2H2, 4 CO-like, 2 ERF, 3 HD-ZIP, 3 YABBYs, G2-like, GATA, GRAS and TCP). In late responsive group, 3 ARR-B, C3H, 6 CO-like, 2 G2-like, 2 HSFs, 2 NACs and TCP. Most of them presented the up-regulated expression at 12 or 24 hours after cold stress implied their important roles for the new growth homeostasis under cold stress. Our study identified the key TFs that function in the regulatory cascades mediating the activation of downstream genes during cold tress tolerance in the paper mulberry. Based on the analysis, we found that the AP2/ERF, bHLH, MYB, NAC and WRKY families might play the central and significant roles during cold stress response in the paper mulberry just as in other species. Meanwhile, many other TF families previously reported as involving in regulation of growth and development, including ARF, DBB, G2-like, GRF, GRAS, LBD, WOX and YAABY exhibited their important potential function in growth regulation under cold stress.

The WOX (WUSCHEL-related homeobox) is a plant-specific transcription factor involved in plant development and stress response. However, few studies have been reported on the WOX gene in woody plants. In this study, 10 BpWOX genes were isolated from paper mulberry by RACE-PCR and categorized into three clades through phylogenetic analysis, ancient, intermediate and WUS clade. Among them, five members had the transcriptional activity detected by yeast one-hybrid and seven were uniquely localized to the nucleus through green fluorescent protein (GFP) observation. The expression patterns of BpWOX genes in different tissues and under diverse treatments were quantified by the qRT-PCR method. Results showed that BpWUS was expressed in the apical bud, stem and root, BpWOX5 and BpWOX7 functioned only in the root tip, and three BpWOXs regulated leaf development redundantly. BpWOX9 and BpWOX10 were induced by indole-3-acetic acid (IAA) or jasmonic acid (JA), while BpWOX2 was repressed by five phytohormones. Interestingly, most BpWOX genes were responsive to the abiotic stress stimuli of drought, salt, cold, and cadmium (CdCl₂). Together, our study revealed that BpWOXs were functionally divergent during paper mulberry development and environmental adaptation, which might be related to their evolutionary relationships. Our work will benefit the systematic understanding of the precise function of WOX in plant development and environmental stress responses.

Heterophylly is regard as an important adaptive mechanism in response to different environments within plants. However, the genetic mechanisms responsible for heterophylly in woody plants are still poorly understood. Herein, the divergence of heterophyllous leaves was investigated at morphogenesis and using microdissection and physiological indexes in paper mulberry, and the genetic basis of heterophylly was further revealed combined with genome-wide association study (GWAS), transcriptome analysis and weighted gene coexpression network analysis (WGCNA). Our results revealed that the flavonoid content and antioxidant activity increased gradually from the entire leaf to the palmatisect leaf, while the hormone content and net photosynthetic rate decreased. Through GWAS and transcriptome analysis, a total of 98 candidate genes and 2338 differentially expressed genes associated with heterophylly were identified. Importantly, we uncovered critical variations in the candidate genes Bp07g0981 (WOX) and Bp07g0920 (HHO), along with significant differences in haplotypes and expression levels among heterophyllous leaves. Our results also suggested that the genes involved in hormone signaling pathways, antioxidant activity, and flavonoid metabolism might be closely related to the heterophylly of paper mulberry, which could account for the physiological data. Indeed, CR-wox mutant lines showed significant changes in leaf phenotypes, and differential expression profile analysis also highlighted the expression of genes related to phytohormones and transcription factors. Together, the genetic variations and candidate genes detected in this study provide novel insights into the genetic mechanism of heterophylly, and would improve the understanding of eco-adaptability in heterophyllous woody plants.

YABs play an important role in the leaf development of the paper mulberry (Broussonetia papyrifera) and of the heterophylly. Thus, we investigated the function of BpYABs. Gene cloning, phylogenetic analysis, motif identification, subcellular localization, transactivation activity assay, qRT-PCR, in situ hybridization, and ectopic expression were used in our study. Six BpYABs were isolated, and four of them had transcriptional activity. BpYAB1, BpYAB3, BpYAB4, and BpYAB5 were localized to the nucleus. BpYAB1 was only expressed in the flower, while BpYAB6 was not expressed in any detected tissues; the four remaining BpYABs were expressed in the bud, leaf and flower, and their expression level decreased with leaf development. Further in situ hybridization showed that BpYAB3 and BpYAB5 were expressed in the vascular tissues and lamina, but neither showed the adaxial–abaxial polarity distribution pattern in the mature leaf lamina. Ectopic expression of BpYAB2, BpYAB3, BpYAB4 and BpYAB5 induced increased expression of AtWOX1 and caused the leaf of Arabidopsis to become smaller and curl downwards. Ectopic expression also led to shorter siliques and smaller seeds, but not for BpYAB5. These results suggest that BpYABs have functional divergency and redundancy in regulating leaf and silique development.

Leaves are essential for photosynthesis and transpiration, directly influencing plant growth and development. Leaf morphology, such as length, width, and area, affects photosynthetic efficiency and transpiration rates. In this study, we investigated the role of PtoWOX1 in leaf morphogenesis by generating both overexpression and CRISPR/Cas9 knockout lines in P. tomentosa. The results showed that PtoWOX1A and PtoWOX1B encode nuclear-localized transcription factors highly expressed in young leaves, particularly in palisade and epidermal cells. Knockout of PtoWOX1 resulted in reduced leaf width and area, enlarged upper epidermal cells, and lower stomatal density. Overexpression led to wrinkled leaf surfaces and reduced margin serration. Anatomical analysis revealed altered palisade cell arrangement and increased leaf thickness in knockout lines, accompanied by higher chlorophyll content and enhanced photosynthetic rates. Additionally, PtoWOX1A interacts with PtoYAB3B, suggesting a complex that regulates leaf margin development. These findings clarify the function of PtoWOX1 in regulating mid-lateral axis development and leaf margin morphology and provide new insights for the molecular breeding of poplar.