The relict woody plant genus Liriodendron contains two endangered species, Liriodendron chinense and Liriodendron tulipifera. Understanding the molecular mechanisms involved in early embryo development is important for horticultural and ecological research, particularly for the development of improved somatic embryogenesis systems. However, the specific molecular processes underlying embryogenesis in these species remain largely unexplored. To address this, we investigated expression of the WOX (WUSCHEL-related homeobox) gene family of transcription factors throughout somatic embryogenesis. We confirmed expression of eight out of 11 novel candidate LcWOX genes in L.chinense using qRT-PCR and examined spatiotemporal expression patterns of the expressed genes using stable reporter lines that had been transformed with different LcWOX promoters driving GUS expression. We observed embryo developmental stages and expression patterns that broadly correlated with those reported for Arabidopsis somatic embryogenesis. LcWUS was weakly expressed during the transition stage and was predominantly restricted to the apical meristem. LcWOX5 was specifically expressed in the root meristem and restricted to the cotyledons thereafter, and LcWOX4 expression was restricted to the vascular tissue of cotyledonary embryos. In contrast, LcWOX9 was expressed in the embryonic callus and the entire embryonic cell mass, then became restricted to the basal cells, indicating a potential role in regulating embryonic maintenance. Our findings provide insights into spatiotemporally specific WOX transcription and shed new light on potential functions of WOX genes during Liriodendron somatic embryogenesis.

Liriodendron chinense is a tree species of the Magnoliaceae family, an ancient relict plant mainly used for landscaping and timber production due to its excellent material properties and ornamental value. The cytokinin oxidase/dehydrogenase (CKX) enzyme regulates cytokinin levels and plays an important role in plant growth, development, and resistance. However, too-high or too-low temperatures or soil drought can limit the growth of L. chinense, representing a key issue for research. Here, we identified the CKX gene family in the L. chinense genome and examined its transcriptional responses to cold, drought, and heat stresses. A total of five LcCKX genes, distributed on four chromosomes and divided into three phylogenetic groups, were identified across the whole L. chinense genome. Further analysis showed that multiple hormone- and stress-responsive cis-acting elements are located in the promoter regions of LcCKXs, indicating a potential role of these LcCKXs in plant growth, development, and response to environmental stresses. Based on existing transcriptome data, LcCKXs, especially LcCKX5, were found to transcriptionally respond to cold, heat, and drought stresses. Furthermore, quantitative reverse-transcription PCR (qRT-PCR) showed that LcCKX5 responds to drought stress in an ABA-dependent manner in stems and leaves and in an ABA-independent manner in roots. These results lay a foundation for functional research on LcCKX genes in the resistance breeding of the rare and endangered tree species of L. chinense.

The clustered, regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) system is the most widely used gene-editing tool to date. However, its application in the genetic improvement of forestry trees has been largely limited. Here, we first established a highly efficient multi-target editing system in the magnoliid woody plant Liriodendron tulipifera. Using phytoene desaturase gene (PDS) as an example, we systematically compared CRISPR/Cas9 and CRSPR/Cpf1 expression systems for loss-of-function analysis and conducted genetic transformations using transient and stable transformation. Ultimately, our findings indicated that the CRISPR/Cas9 system, when applied to transformation based on single-cell-originated somatic embryogenesis, yielded the highest gene-editing efficiency, with mutation rates of nearly 100%. Furthermore, we obtained a total of 137 regeneration plantlets via somatic embryogenesis, of which 82.48% exhibited an albino phenotype. The Illumina sequencing results of albino seedlings and the callus tissue obtained from dedifferentiation of mutant plants revealed that the mutation at the T1 target site was homozygous. These results indicate that CRISPR/Cas9-based multiplex genome-editing technology can not only accelerate the identification of gene function but also be incorporated into the genetic improvement and breeding of tulip trees, supporting the scale propagation of genome-edited plantlets via somatic embryogenesis.

The phytohormone abscisic acid (ABA) plays important roles in response to abiotic and biotic stresses in plants. Pyrabactin resistance 1-like (PYR/PYL) proteins are well-known as ABA receptors, which are responsible for ABA signal transduction. Nevertheless, the characteristics of PYL genes from Liriodendron chinense, an endangered timber tree, remain unclear in coping with various stresses. In this study, five PYLs were identified from the genome of Liriodendron chinense by sequence alignment and conserved motif analysis, which revealed that these LcPYLs contain a conserved gate and latch motif for ABA binding. The LcPYL promoters possess a series of cis-acting elements involved in response to various hormone and abiotic stresses. Moreover, the transcriptome data of Liriodendron hybrid leaves reveal that LcPYL genes specifically transcript under different abiotic stresses; Lchi11622 transcription was induced by drought and cold treatment, and Lchi01385 and Lchi16997 transcription was upregulated under cold and hot stress, respectively. Meanwhile, the LcPYLs with high expression levels shown in the transcriptomes were also found to be upregulated in whole plants treated with the same stresses tested by qPCR. Moreover, under biotic stress caused by scale insect and whitefly, Liriodendron hybrid leaves exhibited a distinct phenotype including disease spots that are dark green in the middle and yellow on the margin; the qPCR results showed that the relative expression levels of Lchi13641 and Lchi11622 in infected leaves were upregulated by 1.76 and 3.75 folds relative to normal leaves, respectively. The subcellular localizations of these stress-responsive LcPYLs were also identified in protoplasts of Liriodendron hybrid. These results provide a foundation to elucidate the function of PYLs from this elite tree species and assist in understanding the molecular mechanism of Liriodendron hybrid in dealing with abiotic and biotic stresses. In future research, the detailed biological function of LcPYLs and the genetic redundancy between LcPYLs can be explored by gene overexpression and knockout based on this study.

The F-box gene family is one of the largest gene families in plants, and it plays a crucial role in regulating plant development, reproduction, cellular protein degradation, and response to biotic and abiotic stresses. Despite their significance, a comprehensive analysis of the F-box gene family in Liriodendron chinense and other magnoliaceae species has not been reported. In this study, we report for the first time the identification of 144 full-length F-box genes in L. chinense. Based on specific domains and phylogenetic analyses, these genes were divided into 10 distinct subfamilies. We further analyzed their gene structure, conserved domain and chromosome distribution, genome-wide replication events, and collinearity. Additionally, based on GO analysis, we found that F-box genes exhibit functional specificity, with a significant proportion of them being involved in protein binding (GO:0005515), suggesting that F-box genes may play an important role in gene regulation in L. chinense. Transcriptome data and q-PCR results also showed that F-box genes are involved in the development of multiple tissues in L. chinense, regulate the somatic embryogenesis of Liriodendron hybrids, and play a pivotal role in abiotic stress. Altogether, these findings provide a foundation for understanding the biological function of F-box genes in L. chinense and other plant species.