Transcriptome analysis and functional validation reveal the novel role of LhCYCL in axillary bud development in hybrid Liriodendron.

Shoot branching significantly influences yield and timber quality in woody plants, with hybrid Liriodendron being particularly valuable due to its rapid growth. However, understanding of the mechanisms governing shoot branching in hybrid Liriodendron remains limited. In this study, we systematically examined axillary bud development using morphological and anatomical approaches and selected four distinct developmental stages for an extensive transcriptome analysis. A total of 9,449 differentially expressed genes have been identified, many of which are involved in plant hormone signal transduction pathways. Additionally, we identified several transcription factors downregulated during early axillary bud development, including a noteworthy gene annotated as CYC-like from the TCP TF family, which emerged as a strong candidate for modulating axillary bud development. Quantitative real-time polymerase chain reaction results confirmed the highest expression levels of LhCYCL in hybrid Liriodendron axillary buds, while histochemical ß-glucuronidase staining suggested its potential role in Arabidopsis thaliana leaf axil development. Ectopic expression of LhCYCL in A. thaliana led to an increase of branches and a decrease of plant height, accompanied by altered expression of genes involved in the plant hormone signaling pathways. This indicates the involvement of LhCYCL in regulating shoot branching through plant hormone signaling pathways. In summary, our results emphasize the pivotal role played by LhCYCL in shoot branching, offering insights into the function of the CYC-like gene and establishing a robust foundation for further investigations into the molecular mechanisms governing axillary bud development in hybrid Liriodendron.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense.

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Physiological and transcriptomic analysis revealed that the accumulation of reactive oxygen species caused the low temperature sensitivity of Liriodendron × sinoamericanum.

Liriodendron × sinoamericanum is widely cultivated in southern China as an excellent wood and garden ornamental trees. However, its intolerance to low temperature limits its application to high latitudes. Understanding the molecular mechanism of low temperature sensitivity of Liriodendron × sinoamericanum is very important for its further application. In this study, combined with physiological and transcriptomic analysis, it was revealed that low temperature stress can lead to water loss and decreased photosynthetic capacity of Liriodendron × sinoamericanum leaves. The accelerated accumulation of reactive oxygen species (ROS) caused by the imbalance of cell REDOX homeostasis is one of the important reasons for the low temperature sensitivity. Further analysis showed that several transcription factors could be involved in regulating the synthesis and degradation of ROS, among which LsNAC72 and LsNAC73a could regulate the accumulation of O2- and H2O2 in leaves by affecting the expression level of LsAPX, LsSOD, LsPAO, and LsPOD.

Genomic survey and expression analysis of LcARFs reveal multiple functions to somatic embryogenesis in Liriodendron.

BACKGROUND: Auxin response factors (ARFs) are critical transcription factors that mediate the auxin signaling pathway and are essential for regulating plant growth. However, there is a lack of understanding regarding the ARF gene family in Liriodendron chinense, a vital species in landscaping and economics. Thus, further research is needed to explore the roles of ARFs in L. chinense and their potential applications in plant development. RESULT: In this study, we have identified 20 LcARF genes that belong to three subfamilies in the genome of L. chinense. The analysis of their conserved domains, gene structure, and phylogeny suggests that LcARFs may be evolutionarily conserved and functionally similar to other plant ARFs. The expression of LcARFs varies in different tissues. Additionally, they are also involved in different developmental stages of somatic embryogenesis. Overexpression of LcARF1, LcARF2a, and LcARF5 led to increased activity within callus. Additionally, our promoter-GFP fusion study indicated that LcARF1 may play a role in embryogenesis. Overall, this study provides insights into the functions of LcARFs in plant development and embryogenesis, which could facilitate the improvement of somatic embryogenesis in L. chinense. CONCLUSION: The research findings presented in this study shed light on the regulatory roles of LcARFs in somatic embryogenesis in L. chinense and may aid in accelerating the breeding process of this tree species. By identifying the specific LcARFs involved in different stages of somatic embryogenesis, this study provides a basis for developing targeted breeding strategies aimed at optimizing somatic embryogenesis in L. chinense, which holds great potential for improving the growth and productivity of this economically important species.

The identification and functional characterization of the LcMCT gene from Liriodendron chinense reveals its potenatial role in carotenoids biosyanthesis.

Terpenoids are not only important component of plant floral scent, but also indispensable elements in the formation of floral color. The petals of Liriodendron chinense are rich in tetraterpene carotenoids and release large amounts of volatile monoterpene and sesquiterpene compounds during full blooming stage. However, the mechanism of terpenoid synthesis is not clear in L. chinense. In this study, we identified a LcMCT gene and characterized its potential function in carotenoids biosynthesis. A total of 2947 up-regulated differentially expressed genes (DEGs) were discerned from the transcriptomic data of L. chinense petals, with a significant enrichment of DEGs related to plant hormone signal transduction and terpenoid backbone biosynthesis. After comprehensive analysis on these DEGs, the LcMCT gene was selected for subsequent function characterization. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) results showed that LcMCT was expressed at the highest level in the petals during full blooming stage, suggesting a possible role in carotenoids biosynthesis and volatile terpenoid biosynthesis. Subcellular localization showed that the LcMCT protein was localized in the chloroplast. Overexpression of LcMCT in Arabidopsis thaliana affected the expression levels of MEP pathway genes. Moreover, the MCT enzyme activity and carotenoids contents in transgenic A. thaliana were increased by 69.27% and 15.57%, respectively. These results suggest that LcMCT promotes the biosynthesis of terpenoid precursors via the MEP pathway. Our work lays a foundation for exploring the mechanism of terpenoid synthesis in L. chinense.

Genome-wide identification of GROWTH-REGULATING FACTORs in Liriodendron chinense and functional characterization of LcGRF2 in leaf size regulation.

GROWTH-REGULATING FACTORs (GRFs) play a pivotal role in the regulation of leaf size in plants and have been widely reported in plants. However, their specific functions in leaf size regulation in Liriodendron chinense remains unclear. Therefore, in this study, we identified GRF genes on a genome-wide scale in L. chinense to characterize the roles of LcGRFs in regulating leaf size. A total of nine LcGRF genes were identified, and these genes exhibited weak expression in mature leaves but strong expression in shoot apex. Notably, LcGRF2 exhibited the highest expression level in the shoot apex of L. chinense. Further RT-qPCR assay revealed that the expression level of LcGRF2 gradually decreased along with the leaf development process, and also displayed a gradient along the leaf proximo-distal and medio-lateral axes. Furthermore, overexpression of LcGRF2 in Arabidopsis thaliana resulted in increased leaf size, and significantly up-regulated the expression of genes involved in cell division like AtCYCD3;1, AtKNOLLE, and AtCYCB1;1, indicating that LcGRF2 may influence leaf size by promoting cell proliferation. This work contributes to a better understanding of the roles and molecular mechanisms of LcGRFs in the regulation of leaf size in L. chinense.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories.

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Genome-Wide Identification and Expression Analysis of TPS Gene Family in Liriodendron chinense.

Terpenoids play a key role in plant growth and development, supporting resistance regulation and terpene synthase (TPS), which is the last link in the synthesis process of terpenoids. Liriodendron chinense, commonly called the Chinese tulip tree, is a rare and endangered tree species of the family Magnoliaceae. However, the genome-wide identification of the TPS gene family and its transcriptional responses to development and abiotic stress are still unclear. In the present study, we identified a total of 58 TPS genes throughout the L. chinense genome. A phylogenetic tree analysis showed that they were clustered into five subfamilies and unevenly distributed across six chromosomes. A cis-acting element analysis indicated that LcTPSs were assumed to be highly responsive to stress hormones, such as methyl jasmonate (MeJA) and abscisic acid (ABA). Consistent with this, transcriptome data showed that most LcTPS genes responded to abiotic stress, such as cold, drought, and hot stress, at the transcriptional level. Further analysis showed that LcTPS genes were expressed in a tissue-dependent manner, especially in buds, leaves, and bark. Quantitative reverse transcription PCR (qRT-PCR) analysis confirmed that LcTPS expression was significantly higher in mature leaves compared to young leaves. These results provide a reference for understanding the function and role of the TPS family, laying a foundation for further study of the regulation of TPS in terpenoid biosynthesis in L. chinense.

Overexpression of the Liriodendron tulipifera BOP2 Gene (LtuBOP2) Affects Leaf Margin Development in Transgenic Arabidopsis thaliana.

BLADE-ON-PETIOLE 2 (BOP2) plays a pivotal role in leaf morphogenesis. Liriodendron tulipifera is a suitable model for exploring the molecular mechanisms underlying leaf serration formation, which are largely unknown. Here, we isolated the full-length LtuBOP2 gene and its promoter from L. tulipifera and characterized its function in leaf morphogenesis through multidimensional approaches. The spatiotemporal expression pattern of LtuBOP2 indicated the high expression of LtuBOP2 in stems and leaf buds. We constructed LtuBOP2 promoter, fused the promoter sequences to the ß-glucuronidase (GUS) gene, and then transformed them into Arabidopsis thaliana. Histochemical GUS staining results indicated that GUS activity was higher in petioles and the main vein. LtuBOP2 overexpression in A. thaliana caused moderate serration in the leaf tip, owing to the increased number of abnormal lamina epidermal cells and defective vascular tissue, thus indicating a novel role of BOP2. The ectopic expression of LtuBOP2 in A. thaliana promoted the expression of the lateral organ boundary gene ASYMMETRIC LEAVES2 (AS2) and inhibited JAGGED (JAG) and CUP-SHAPED COTYLEDON2 (CUC2) expression to establish leaf proximal-distal polarity. Moreover, LtuBOP2 participated in leaf serration formation by promoting the antagonistic relationship between KNOX I and hormones during leaf margin development. Our findings revealed the role of LtuBOP2 in the proximal-distal polarity formation and development of leaf margin morphology, providing new insights into the regulatory mechanisms of the leaf formation development of L. tulipifera.

Genomic Survey of Heat Shock Proteins in Liriodendron chinense Provides Insight into Evolution, Characterization, and Functional Diversities.

Heat shock proteins (HSPs) are conserved molecular chaperones whose main role is to facilitate the regulation of plant growth and stress responses. The HSP gene family has been characterized in most plants and elucidated as generally stress-induced, essential for their cytoprotective roles in cells. However, the HSP gene family has not yet been analyzed in the Liriodendron chinense genome. In current study, 60 HSP genes were identified in the L. chinense genome, including 7 LchiHSP90s, 23 LchiHSP70s, and 30 LchiHSP20s. We investigated the phylogenetic relationships, gene structure and arrangement, gene duplication events, cis-acting elements, 3D-protein structures, protein-protein interaction networks, and temperature stress responses in the identified L. chinense HSP genes. The results of the comparative phylogenetic analysis of HSP families in 32 plant species showed that LchiHSPs are closely related to the Cinnamomum kanehirae HSP gene family. Duplication events analysis showed seven segmental and six tandem duplication events that occurred in the LchiHSP gene family, which we speculated to have played an important role in the LchiHSP gene expansion and evolution. Furthermore, the Ka/Ks analysis indicated that these genes underwent a purifying selection. Analysis in the promoter region evidenced that the promoter region LchiHSPs carry many stress-responsive and hormone-related cis-elements. Investigations in the gene expression patterns of the LchiHSPs using transcriptome data and the qRT-PCR technique indicated that most LchiHSPs were responsive to cold and heat stress. In total, our results provide new insights into understanding the LchiHSP gene family function and their regulatory mechanisms in response to abiotic stresses.

Identification and characteristics of SnRK genes and cold stress-induced expression profiles in Liriodendron chinense.

BACKGROUND: The sucrose non-fermenting 1 (SNF1)-related protein kinases (SnRKs) play a vivid role in regulating plant metabolism and stress response, providing a pathway for regulation between metabolism and stress signals. Conducting identification and stress response studies on SnRKs in plants contributes to the development of strategies for tree species that are more tolerant to stress conditions. RESULTS: In the present study, a total of 30 LcSnRKs were identified in Liriodendron chinense (L. chinense) genome, which was distributed across 15 chromosomes and 4 scaffolds. It could be divided into three subfamilies: SnRK1, SnRK2, and SnRK3 based on phylogenetic analysis and domain types. The LcSnRK of the three subfamilies shared the same Ser/Thr kinase structure in gene structure and motif composition, while the functional domains, except for the kinase domain, showed significant differences. A total of 13 collinear gene pairs were detected in L. chinense and Arabidopsis thaliana (A. thaliana), and 18 pairs were detected in L. chinense and rice, suggesting that the LcSnRK family genes may be evolutionarily more closely related to rice. Cis-regulation element analysis showed that LcSnRKs were LTR and TC-rich, which could respond to different environmental stresses. Furthermore, the expression patterns of LcSnRKs are different at different times under low-temperature stress. LcSnRK1s expression tended to be down-regulated under low-temperature stress. The expression of LcSnRK2s tended to be up-regulated under low-temperature stress. The expression trend of LcSnRK3s under low-temperature stress was mainly up-or down-regulated. CONCLUSION: The results of this study will provide valuable information for the functional identification of the LcSnRK gene in the future.

Genome-Wide Identification and Expression Analysis of SnRK2 Gene Family in Dormant Vegetative Buds of Liriodendron chinense in Response to Abscisic Acid, Chilling, and Photoperiod.

Protein kinases play an essential role in plants' responses to environmental stress signals. SnRK2 (sucrose non-fermenting 1-related protein kinase 2) is a plant-specific protein kinase that plays a crucial role in abscisic acid and abiotic stress responses in some model plant species. In apple, corn, rice, pepper, grapevine, Arabidopsis thaliana, potato, and tomato, a genome-wide study of the SnRK2 protein family was performed earlier. The genome-wide comprehensive investigation was first revealed to categorize the SnRK2 genes in the Liriodendron chinense (L. chinense). The five SnRK2 genes found in the L. chinense genome were highlighted in this study. The structural gene variants, 3D structure, chromosomal distributions, motif analysis, phylogeny, subcellular localization, cis-regulatory elements, expression profiles in dormant buds, and photoperiod and chilling responses were all investigated in this research. The five SnRK2 genes from L. chinense were grouped into groups (I-IV) based on phylogeny analysis, with three being closely related to other species. Five hormones-, six stress-, two growths and biological process-, and two metabolic-related responsive elements were discovered by studying the cis-elements in the promoters. According to the expression analyses, all five genes were up- and down-regulated in response to abscisic acid (ABA), photoperiod, chilling, and chilling, as well as photoperiod treatments. Our findings gave insight into the SnRK2 family genes in L. chinense and opened up new study options.

Identification, Phylogenetic and Expression Analyses of the AAAP Gene Family in Liriodendron chinense Reveal Their Putative Functions in Response to Organ and Multiple Abiotic Stresses.

In this study, 52 AAAP genes were identified in the L. chinense genome and divided into eight subgroups based on phylogenetic relationships, gene structure, and conserved motif. A total of 48 LcAAAP genes were located on the 14 chromosomes, and the remaining four genes were mapped in the contigs. Multispecies phylogenetic tree and codon usage bias analysis show that the LcAAAP gene family is closer to the AAAP of Amborella trichopoda, indicating that the LcAAAP gene family is relatively primitive in angiosperms. Gene duplication events revealed six pairs of segmental duplications and one pair of tandem duplications, in which many paralogous genes diverged in function before monocotyledonous and dicotyledonous plants differentiation and were strongly purification selected. Gene expression pattern analysis showed that the LcAAAP gene plays a certain role in the development of Liriodendron nectary and somatic embryogenesis. Low temperature, drought, and heat stresses may activate some WRKY/MYB transcription factors to positively regulate the expression of LcAAAP genes to achieve long-distance transport of amino acids in plants to resist the unfavorable external environment. In addition, the GAT and PorT subgroups could involve gamma-aminobutyric acid (GABA) transport under aluminum poisoning. These findings could lay a solid foundation for further study of the biological role of LcAAAP and improvement of the stress resistance of Liriodendron.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis.

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

Overexpression of Liriodendron tulipifera JAG Gene (LtuJAG) Changes Leaf Shapes in Transgenic Arabidopsis thaliana.

In Arabidopsis thaliana, JAGGED (JAG) is a transcription inhibitor that controls the development of leaf polarity and regulates the expression of genes controlling lateral organ formation. Liriodendron tulipifera is an ornamental tree with extraordinary tulip-shaped flowers and goose web-like leaves, this is one of the suitable plants for morphological development research. To investigate the potential functions of the LtuJAG gene, we isolated the full-length LtuJAG from L. tulipifera, transferred it into A. thaliana via agrobacterium-mediated transformation, and monitored its expression pattern. Subcellular localization showed that LtuJAG was located in the nucleus. RT-qPCR assays indicated that LtuJAG was expressed mainly in leaf buds and flowers, but not in mature leaves and stems. GUS staining results showed that LtuJAG was expressed in the shoot apical meristem (SAM). Overexpressing LtuJAG changed A. thaliana leaf shapes, causing a moderate serration and a slight asymmetric distribution in the medio-lateral and proximal-distal axes. Ectopic expression of LtuJAG induced the expression of lateral organ boundary suppressors JAGGED LATERAL ORGANS (JLO) and ARABIDOPSIS THALIANA HOMEOBOX1 (ATH1). It also repressed the expression of the apical meristem suppressor class-1 KNOX gene (KNOX I) and altered endogenous hormone levels. Our results suggest that LtuJAG plays a role in negatively regulating leaf polarity formation in L. tulipifera.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress.

BACKGROUND: Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions. RESULTS: In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus. CONCLUSIONS: This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Transcriptomic profiling of the floral fragrance biosynthesis pathway of Liriodendron and functional characterization of the LtuDXR gene.

Floral fragrance, which has the function of attracting pollinators, is a class of volatile secondary metabolites mainly released by the secretory tissue of petals. Terpenoids are key components of floral volatile substances. Previous studies have shown that there are significant differences in the concentration and composition of volatile floral fragrances, especially terpenoids, between Liriodendron chinense and L. tulipifera. At present, the mechanism by which the synthesis of floral fragrance is regulated in Liriodendron remains unexplored. In this study, we analyzed the differentially expressed genes (DEGs) of L. chinense and L. tulipifera, and identified 130 DEGs related to terpenoid synthesis. A KEGG enrichment analysis of DEGs related to terpenoid biosynthesis revealed that the monoterpenoid biosynthesis pathway was the most significant. We cloned the LtuDXR gene from L. tulipifera using RACE technology. RT-qPCR results showed that the expression of the LtuDXR gene was the highest in the early florescence petals, indicating that the LtuDXR gene may play a role in the synthesis of volatile terpenoids. Subcellular localization showed that the LtuDXR protein is mainly localized in the chloroplast. Overexpression of LtuDXR in Arabidopsis thaliana significantly increased the plant height, DXR enzyme activity, and carotenoid content. In this study, we identified and functionally characterized LtuDXR, which is involved in terpenoid synthesis in Liriodendron. Our work lays the foundation for further exploration of the molecular mechanism by which terpenoid biosynthesis is regulated in Liriodendron.

Genome-Wide Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Petal Pigment Synthesis in Liriodendron.

The MYB transcription factor family is one of the largest families in plants, and its members have various biological functions. R2R3-MYB genes are involved in the synthesis of pigments that yield petal colors. Liriodendron plants are widely cultivated as ornamental trees owing to their peculiar leaves, tulip-like flowers, and colorful petals. However, the mechanism underlying petal coloring in this species is unknown, and minimal information about MYB genes in Liriodendron is available. Herein, this study aimed to discern gene(s) involved in petal coloration in Liriodendron via genome-wide identification, HPLC, and RT-qPCR assays. In total, 204 LcMYB superfamily genes were identified in the Liriodendron chinense genome, and 85 R2R3-MYB genes were mapped onto 19 chromosomes. Chromosome 4 contained the most (10) R2R3-MYB genes, and chromosomes 14 and 16 contained the fewest (only one). MEME analysis showed that R2R3-MYB proteins in L. chinense were highly conserved and that their exon-intron structures varied. The HPLC results showed that three major carotenoids were uniformly distributed in the petals of L. chinense, while lycopene and ß-carotene were concentrated in the orange band region in the petals of Liriodendron tulipifera. Furthermore, the expression profiles via RT-qPCR assays revealed that four R2R3-MYB genes were expressed at the highest levels at the S3P/S4P stage in L. tulipifera. This result combined with the HPLC results showed that these four R2R3-MYB genes might participate in carotenoid synthesis in the petals of L. tulipifera. This work laid a cornerstone for further functional characterization of R2R3-MYB genes in Liriodendron plants.

Gibberellin Oxidase Gene Family in L. chinense: Genome-Wide Identification and Gene Expression Analysis.

GAox is a key enzyme for the transformation of gibberellins, and belongs to the 2-ketoglutarate dependent dioxygenase gene family (2ODD). However, a systematic analysis of GAox in the angiosperm L. chinense has not yet been reported. Here, we identified all LcGAox gene family members in L. chinense, which were classified into the three subgroups of GA20ox, C19GA2ox, and C20GA2ox. Comparison of the gene structure, conserve motifs, phylogenetic relationships, and syntenic relationships of gibberellin oxidase gene families in different species indicated that the gene functional differences may be due to the partial deletion of their domains during evolution. Furthermore, evidence for purifying selection was detected between orthologous GAox genes in rice, grape, Arabidopsis, and L. chinense. Analysis of the codon usage patterns showed that mutation pressure and natural selection might have induced codon usage bias in angiosperms; however, the LcGAox genes in mosses, lycophytes, and ambarella plants exhibited no obvious codon usage preference. These results suggested that the gibberellin oxidase genes were more primitive. The gene expression pattern was analyzed in different organs subjected to multiple abiotic stresses, including GA, abscisic acid (ABA), and chlormequat (CCC) treatment, by RNA-seq and qRT-PCR, and the stress- and phytohormone-responsive cis-elements were counted. The results showed that the synthesis and decomposition of GA were regulated by different LcGAox genes in the vegetative and reproductive organs of L. chinense, and only LcGA2ox1,4, and 7 responded to the NaCl, polyethylene glycol, 4 °C, GA, ABA, and CCC treatment in the roots, stems, and leaves of seedlings at different time periods, revealing the potential role of LcGAox in stress resistance.

Genome-wide identification of MIKC-type genes related to stamen and gynoecium development in Liriodendron.

The organogenesis and development of reproductive organs, i.e., stamen and gynoecium, are important floral characteristics that are closely related to pollinators and reproductive fitness. As a genus from Magnoliaceae, Liriodendron has only two relict species: L. chinense and L. tulipifera. Despite the similar flower shapes of these species, their natural seed-setting rates differ significantly, implying interspecies difference in floral organogenesis and development. MADS-box genes, which participate in floral organogenesis and development, remain unexplored in Liriodendron. Here, to explore the interspecies difference in floral organogenesis and development and identify MADS-box genes in Liriodendron, we examined the stamen and gynoecium primordia of the two Liriodendron species by scanning electron microscopy combined with paraffin sectioning, and then collected two types of primordia for RNA-seq. A total of 12 libraries were constructed and 42,268 genes were identified, including 35,269 reference genes and 6,999 new genes. Monoterpenoid biosynthesis was enriched in L. tulipifera. Genome-wide analysis of 32 MADS-box genes was conducted, including phylogenetic trees, exon/intron structures, and conserved motif distributions. Twenty-six genes were anchored on 17 scaffolds, and six new genes had no location information. The expression profiles of MIKC-type genes via RT-qPCR acrossing six stamen and gynoecium developmental stages indicates that the PI-like, AG/STK-like, SEP-like, and SVP-like genes may contribute to the species-specific differentiation of the organogenesis and development of reproductive organs in Liriodendron. Our findings laid the groundwork for the future exploration of the mechanism underlying on the interspecific differences in reproductive organ development and fitness in Liriodendron.