Identification of Suitable Reference Genes for RT-qPCR Assays in Liriodendron chinense (Hemsl.) Sarg

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.

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.

Identification of Accurate Reference Genes for qRT-PCR Analysis of Gene Expression in Eremochloa ophiuroides under Multiple Stresses of Phosphorus Deficiency and/or Aluminum Toxicity

Centipedegrass (Eremochloa ophiuroides (Munro.) Hack.) is a species originating in China and is an excellent warm-season turfgrass. As a native species in southern China, it is naturally distributed in the phosphorus-deficient and aluminum-toxic acid soil areas. It is important to research the molecular mechanism of centipedegrass responses to phosphorus-deficiency and/or aluminum-toxicity stress. Quantitative Real-Time PCR (qRT-PCR) is a common method for gene expression analysis, and the accuracy of qRT-PCR results depends heavily on the stability of internal reference genes. However, there are still no reported stable and effective reference genes for qRT-PCR analysis of target genes under the acid-soil-related stresses in different organs of centipedegrass. For scientific rigor, the gene used as a reference for any plant species and/or any stress conditions should be first systematically screened and evaluated. This study is the first to provide a group of reliable reference genes to quantify the expression levels of functional genes of Eremochloa ophiuroides under multiple stresses of P deficiency and/or aluminum toxicity. In this study, centipedegrass seedlings of the acid-soil-resistant strain 'E041' and acid-soil-sensitive strain 'E089' were used for qRT-PCR analysis. A total of 11 candidate reference genes (ACT, TUB, GAPDH, TIP41, CACS, HNR, EP, EF1α, EIF4α, PP2A and actin) were detected by qRT-PCR technology, and the stability of candidate genes was evaluated with the combination of four internal stability analysis software programs. The candidate reference genes exhibited differential stability of expression in roots, stems and leaves under phosphorus-deficiency and/or aluminum-toxicity stress. On the whole, the results showed that GAPDH, TIP41 and HNR were the most stable in the total of samples. In addition, for different tissues under various stresses, the selected reference genes were also different. CACS and PP2A were identified as two stable reference genes in roots through all three stress treatments (phosphate deficiency, aluminum toxicity, and the multiple stress treatment of aluminum toxicity and phosphate deficiency). Moreover, CACS was also stable as a reference gene in roots under each treatment (phosphate deficiency, aluminum toxicity, or multiple stresses of aluminum toxicity and phosphate deficiency). In stems under all three stress treatments, GAPDH and EIF4α were the most stable reference genes; for leaves, PP2A and TIP41 showed the two highest rankings in all three stress treatments. Finally, qRT-PCR analysis of the expression patterns of the target gene ALMT1 was performed to verify the selected reference genes. The application of the reference genes identified as internal controls for qRT-PCR analysis will enable accurate analysis of the target gene expression levels and expression patterns in centipedegrass under acid-soil-related stresses.

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.

Genome-wide identification and characterization of LcCCR13 reveals its potential role in lignin biosynthesis in Liriodendron chinense

INTRODUCTION

Wood formation is closely related to lignin biosynthesis. Cinnamoyl-CoA reductase (CCR) catalyzes the conversion of cinnamoyl-CoA to cinnamaldehydes, which is the initiation of the lignin biosynthesis pathway and a crucial point in the manipulation of associated traits. Liriodendron chinense is an economically significant timber tree. Nevertheless, the underlying mechanism of wood formation in it remains unknown; even the number of LcCCR family members in this species is unclear.

MATERIALS AND RESULTS

This study aimed to perform a genome-wide identification of genes(s) involved in lignin biosynthesis in L. chinense via RT-qPCR assays and functional verification. Altogether, 13 LcCCR genes were identified that were divided into four major groups based on structural and phylogenetic features. The gene structures and motif compositions were strongly conserved between members of the same groups. Subsequently, the expression patterns analysis based on RNA-seq data indicated that LcCCR5/7/10/12/13 had high expression in the developing xylem at the stem (DXS). Furthermore, the RT-qPCR assays showed that LcCCR13 had the highest expression in the stem as compared to other tissues. Moreover, the overexpression of the LcCCR13 in transgenic tobacco plants caused an improvement in the CCR activity and lignin content, indicating that it plays a key role in lignin biosynthesis in the stems.

DISCUSSION

Our research lays a foundation for deeper investigation of the lignin synthesis and uncovers the genetic basis of wood formation in L. chinense.

Liriodendron chinense LcMAX1 regulates primary root growth and shoot branching in Arabidopsis thaliana

Strigolactones (SLs) play prominent roles in regulating shoot branching and root architecture in model plants. However, their roles in non-model (particularly woody) plants remain unclear. Liriodendron chinense is a timber tree species widely planted in southern China. The outturn percentage and wood quality of L. chinense are greatly affected by the branching characteristics of its shoot, and the rooting ability of the cuttings is key for its vegetative propagation. Here, we isolated and analyzed the function of the MORE AXILLARY GROWTH 1 (LcMAX1) gene, which is involved in L. chinense SL biosynthesis. RT-qPCR showed that LcMAX1 was highly expressed in the roots and axillary buds. LcMAX1 was located in the endoplasmic reticulum (ER) and nucleus. LcMAX1 ectopic expression promoted primary root growth, whereas there were no phenotypic differences in shoot branching between transgenic and wild-type (WT) A. thaliana plants. LcMAX1 overexpression in the max1 mutant restored them to the WT A. thaliana phenotypes. Additionally, AtPIN1, AtPIN2, and AtBRC1 expressions were significantly upregulated in transgenic A. thaliana and the max1 mutant. It was therefore speculated that LcMAX1 promotes primary root growth by regulating expression of auxin transport-related genes in A. thaliana, and LcMAX1 inhibits shoot branching by upregulating expression of AtBRC1 in the max1 mutant. Altogether, these results demonstrated that the root development and shoot branching functions of LcMAX1 were similar to those of AtMAX1. Our findings provide a foundation for obtaining further insights into root and branch development in L. chinense.

Genome-wide identification of XTH genes in Liriodendron chinense and functional characterization of LcXTH21

Liriodendron chinense is a relic tree species of the family Magnoliaceae with multiple uses in timber production, landscape decoration, and afforestation. L. chinense often experiences drought stress in arid areas. However, the molecular basis underlying the drought response of L. chinense remains unclear. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in drought stress resistance. Hereby, to explore the drought resistance mechanism of L. chinense, we identify XTH genes on a genome-wide scale in L. chinense. A total of 27 XTH genes were identified in L. chinense, and these genes were classified into three subfamilies. Drought treatment and RT-qPCR analysis revealed that six LcXTH genes significantly responded to drought stress, especially LcXTH21. Hence, we cloned the LcXTH21 gene and overexpressed it in tobacco via gene transfer to analyze its function. The roots of transgenic plants were more developed than those of wild-type plants under different polyethylene glycol (PEG) concentration, and further RT-qPCR analysis showed that LcXTH21 highly expressed in root compared to aboveground organs, indicating that LcXTH21 may play a role in drought resistance through promoting root development. The results of this study provide new insights into the roles of LcXTH genes in the drought stress response. Our findings will also aid future studies of the molecular mechanisms by which LcXTH genes contribute to the drought response.

Overexpression of LtuHB6 from Liriodendron tulipifera causes lobed-leaf formation in Arabidopsis thaliana

Liriodendron tulipifera L. is an ornamental tree species with extraordinarily lobed leaves. However, the mechanisms underlying lobed leaf formation in plants remain unclear. The transcription factor, ARABIDOPSIS THALIANA HOMEBOX 6 (HB6), plays a role in regulating leaf margin development. HB6 is involved in cell division and differentiation of developmental organs and negatively regulates abscisic acid (ABA) signal transmission under external abiotic stress; it is unclear whether HB6 performs a pivotal role in leaf morphogenesis in L. tulipifera. In this study, full-length LtuHB6 from L. tulipifera was heterologously expressed in tobacco and Arabidopsis thaliana; its expression pattern was analyzed to determine its potential role in leaf development. In addition, LtuHB6 is localized in the nucleus and cell membrane of tobacco leaves. The expression of LtuHB6 was highest in mature leaves compared to the other stages of leaf development (bud growth, young leaves, and leaf senescence). Transgenic A. thaliana plants overexpressing LtuHB6 exhibited an abnormal phenotype with lobed leaves. Moreover, LtuHB6 overexpression significantly affected the expression of seven genes related to leaf serration in the initial stage of leaf primordia and altered the expression levels of hormonal genes. Our findings indicate that LtuHB6 is an essential regulatory factor in L. tulipifera lobed-leaf formation and is involved in regulating and responding to hormones.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01254-9.

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 the LcCUC2-like gene in Arabidopsis thaliana alters the cotyledon morphology and increases rosette leaf number

BACKGROUND

The unique 'mandarin jacket' leaf shape is the most famous trait of Liriodendron chinense and this characteristic gives L. chinense aesthetic and landscaping value. However, the underlying regulatory mechanism of genes involved in the leaf development of L. chinense has remained unclear.

METHODS

Based on transcriptome data of leaves at different developmental stages from L. chinense, we identified differentially expression genes (DEGs) functioning in leaf development. A candidate gene named LcCUC2-like (LcCUC2L) had high similarity in sequence with Arabidopsis thaliana CUC2, and used for further research. We isolated the full-length LcCUC2L gene and its promoter from L. chinense. Subsequently, we analyzed the function of the LcCUC2L gene and its promoter activity via transformation into A. thaliana.

RESULTS

In this study, we found that the LcCUC2L and AtCUC2 are homologous in sequence but not homologous in function. Unlike the role of AtCUC2 in leaf serration and SAM formation, the LcCUC2L mainly regulates cotyledon development and rosette leaf number. Histochemical β-glucuronidase (GUS) staining revealed that LcCUC2L was expressed in the cotyledons of A. thaliana seedlings, indicating that the LcCUC2L may play a role in cotyledon development. Ectopic expression of LcCUC2L resulted in long, narrow cotyledons without petioles, abnormal lamina epidermis cells and defective vascular tissue in cotyledons, and these results were consistent with the LcCUC2L expression pattern. Further analysis showed that overexpression of LcCUC2L also induced numerous rosette leaves. Also, LcCUC2L and other related genes showed a severe response in L. chinense by introducing exogenous auxin stimulation, partly revealed that LcCUC2L affects the leaf development by regulating the auxin content.

CONCLUSIONS

These results suggest that LcCUC2L may play a critical role in leaf development and morphogenesis in L. chinense, and our findings provide insight into the molecular mechanisms of leaf development in L. chinense.

The Roles of microRNA-Long Non-coding RNA-mRNA Networks in the Regulation of Leaf and Flower Development in Liriodendron chinense

The leaf and the flower are vital plant organs owing to their roles in photosynthesis and reproduction. Long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and transcription factors (TFs) are very important to the development of these organs. Liriodendron chinense is a common ornamental tree species in southern China with an unusual leaf shape and tulip-like flowers. The genetic mechanisms underlying leaf and flower development in L. chinense and the miRNA-lncRNA-TF regulatory networks are poorly studied. Through the integration and analysis of different types of sequencing data, we identified the miRNA-lncRNA-TF regulatory networks that were related to leaf and flower development. These networks contained 105 miRNAs, 258 lncRNAs, 393 TFs, and 22 endogenous target mimics. Notably, lch-lnc7374-miR156h-SPL3 and lch-lnc7374-miR156j-SPL9 were potential regulators of stamen and pistil development in L. chinense, respectively. miRNA-lncRNA-mRNA regulatory networks were shown to impact anther development, male and female fertility, and petal color by regulating the biosynthesis of phenylpropanoid metabolites. Phenylpropanoid metabolite biosynthesis genes and TFs that were targeted by miRNAs and lncRNAs were differentially expressed in the leaf and flower. Moreover, RT-qPCR analysis confirmed 22 differentially expressed miRNAs, among which most of them showed obvious leaf or flower specificity; miR157a-SPL and miR160a-ARF module were verified by using RLM-RACE, and these two modules were related to leaf and flower development. These findings provide insight into the roles of miRNA-lncRNA-mRNA regulatory networks in organ development and function in L. chinense, and will facilitate further investigation into the regulatory mechanisms of leaf and flower development in L. chinense.

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.

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 survey and identification of AP2/ERF genes involved in shoot and leaf development in Liriodendron chinense

BACKGROUND

Liriodendron chinense is a distinctive ornamental tree species due to its unique leaves and tulip-like flowers. The discovery of genes involved in leaf development and morphogenesis is critical for uncovering the underlying genetic basis of these traits. Genes in the AP2/ERF family are recognized as plant-specific transcription factors that contribute to plant growth, hormone-induced development, ethylene response factors, and stress responses.

RESULTS

In this study, we identified 104 putative AP2/ERF genes in the recently released L. chinense genome and transcriptome database. In addition, all 104 genes were grouped into four subfamilies, the AP2, ERF, RAV, and Soloist subfamilies. This classification was further supported by the results of gene structure and conserved motif analyses. Intriguingly, after application of a series test of cluster analysis, three AP2 genes, LcERF 94, LcERF 96, and LcERF 98, were identified as tissue-specific in buds based on the expression profiles of various tissues. These results were further validated via RT-qPCR assays and were highly consistent with the STC analysis. We further investigated the dynamic changes of immature leaves by dissecting fresh shoots into seven discontinuous periods, which were empirically identified as shoot apical meristem (SAM), leaf primordia and tender leaf developmental stages according to the anatomic structure. Subsequently, these three candidates were highly expressed in SAM and leaf primordia but rarely in tender leaves, indicating that they were mainly involved in early leaf development and morphogenesis. Moreover, these three genes displayed nuclear subcellular localizations through the transient transformation of tobacco epidermal cells.

CONCLUSIONS

Overall, we identified 104 AP2/ERF family members at the genome-wide level and discerned three candidate genes that might participate in the development and morphogenesis of the leaf primordium in L. chinense.

Screening of stable internal reference gene of Quinoa under hormone treatment and abiotic stress

Real-time quantitative polymerase chain reaction is the most commonly used method to accurately detect gene expression patterns. The method requires stable internal reference genes to standardize the data. However, studies have shown that there is no stable expression of internal reference genes in different tissues and under different treatments. Therefore, in order to study the optimal reference genes of quinoa under different hormones and abiotic stress, leaves and stems from quinoa seedlings treated with low temperature (4 °C), salt (200 mmol/L) and abscisic acid (200 mmol/L) were used as experimental materials. Using ACT-1, eIF, EF1α, GAPDH, TUA, TUB-9, TUB-1, H2A and L8-1 as candidate reference genes, the expression stability of these 9 quinoa candidate reference genes under different hormone treatment and abiotic stress was evaluated by using geNorm, NormFinder and BestKeeper software. The results showed that TUB-1 gene under salt stress, L8-1 gene under low temperature stress, EF-1α gene induced by ABA. PLIM2c WLIM1and WLIM2b were selected to verify the candidate internal reference genes, and finally the expression of GAPDH was most unstable under the three treatments, which was not suitable to be the internal reference gene of quinoa under specific conditions, while EF1α showed good stability under the three different treatments and was suitable to be used as the internal reference gene. In conclusion, the results of this study could provide an important reference for quantifying the expression level of reference genes in quinoa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01094-z.

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.

Isolation, expression, and functional analysis of the geranylgeranyl pyrophosphate synthase (GGPPS) gene from Liriodendron tulipifera

Terpenoids are important secondary metabolites in plants and are involved in stress responses and pollinator attraction. Geranylgeranyl pyrophosphate synthase (GGPPS) is a key synthase in the 2C-methyl-D-erythritol-4-phosphate (MEP) pathway of terpenoid synthesis, catalyzing the synthesis of diterpenoids. Liriodendron tulipifera is a nectar plant in North America. Little is known about the key genes involved in the biosynthetic pathways of terpenoids, the precursors of most compounds related to nectar, fragrance and coloring in flowers in L. tulipifera. In this study, the LtuGGPPS2 gene and its promoter (LtuGGPPS2-pro) were cloned from L. tulipifera. The results of sequence alignment showed that the LtuGGPPS2 gene is highly homologous to GGPPS genes of other plants. Subcellular localization analysis showed that the LtuGGPPS2 protein localizes to chloroplasts, suggesting that the LtuGGPPS2 gene is probably related to carotenoid and chlorophyll synthesis. Based on tissue expression profiles revealed by RT-qPCR, the expression level of the LtuGGPPS2 gene was highest in petals. These results were consistent with the changes in volatile and nonvolatile terpenoids in the flowers of L. tulipifera. GUS staining to examine the LtuGGPPS2 promoter indicated that it is responsive to hormones. Overexpression of the LtuGGPPS2 gene increased the carotenoid content and GGPPS enzyme activity in Arabidopsis thaliana, indicating that LtuGGPPS2 is the key terpenoid synthase in the flowers of L. tulipifera. Our findings lay a foundation for further functional analysis of the LtuGGPPS2 gene and deeper investigation of the terpenoid biosynthetic pathway in L. tulipifera.

Morphological, phenological, and transcriptional analyses provide insight into the diverse flowering traits of a mutant of the relic woody plant Liriodendron chinense

Flowering is crucial to plant reproduction and controlled by multiple factors. However, the mechanisms underlying the regulation of flowering in perennial plants are still largely unknown. Here, we first report a super long blooming 1 (slb1) mutant of the relict tree Liriodendron chinense possessing a prolonged blooming period of more than 5 months, in contrast to the 1 month blooming period in the wild type (WT). Phenotypic characterization showed that earlier maturation of lateral shoots was caused by accelerated axillary bud fate, leading to the phenotype of continuous flowering in slb1 mutants. The transcriptional activity of genes related to hormone signaling (auxin, cytokinin, and strigolactone), nutrient availability, and oxidative stress relief further indicated active outgrowth of lateral buds in slb1 mutants. Interestingly, we discovered a unique FT splicing variant with intron retention specific to slb1 mutants, representing a potential causal mutation in the slb1 mutants. Surprisingly, most slb1 inbred offspring flowered precociously with shorter juvenility (~4 months) than that (usually 8-10 years) required in WT plants, indicating heritable variation underlying continuous flowering in slb1 mutants. This study reports an example of a perennial tree mutant that flowers continuously, providing a rare resource for both breeding and genetic research.

Selection and Optimization of Reference Genes for MicroRNA Expression Normalization by qRT-PCR in Chinese Cedar (Cryptomeria fortunei) under Multiple Stresses

MicroRNA (miRNA) expression analysis is very important for investigating its functions. To date, no research on reference genes (RGs) for miRNAs in gymnosperms, including Cryptomeria fortunei, has been reported. Here, ten miRNAs (i.e., pab-miR159a, cln-miR162, cas-miR166d, pab-miR395b, ppt-miR894, cln-miR6725, novel1, novel6, novel14 and novel16) and three common RGs (U6, 5S and 18S) were selected as candidate RGs. qRT-PCR was used to analyse their expressions in C. fortunei under various experimental conditions, including multiple stresses (cold, heat, drought, salt, abscisic acid and gibberellin) and in various tissues (roots, stems, tender needles, cones and seeds). Four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) were employed to assess the stability of candidate RG expression; the geometric mean and RefFinder program were used to comprehensively evaluate RG stability. According to the results, novel16, cln-miR6725, novel1 and U6 were the most stable RGs for studying C. fortunei miRNA expression. In addition, the expression of three target miRNAs (aly-miR164c-5p, aly-miR168a-5p and smo-miR396) was examined to verify that the selected RGs are suitable for miRNA expression normalisation. This study may aid further investigations of miRNA expression/function in the response of C. fortunei to abiotic stress and provides an important basis for the standardisation of miRNA expression in other gymnosperm species.

Selection and Verification of Appropriate Reference Genes for Expression Normalization in Cryptomeria fortunei under Abiotic Stress and Hormone Treatments

Cryptomeria fortunei has become one of the main timber afforestation species in subtropical high-altitude areas of China due to its fast growth, good material quality, and strong adaptability, showing broad application prospects. Quantitative real-time PCR (qRT-PCR) is the most accurate and widely used gene expression evaluation technique, and selecting appropriate reference genes (RGs) is essential for normalizing qRT-PCR results. However, suitable RGs for gene expression normalization in C. fortunei have not been reported. Here, we tested the expression stability for 12 RGs in C. fortunei under various experimental conditions (simulated abiotic stresses (cold, heat, drought, and salinity) and hormone treatments (methyl jasmonate, abscisic acid, salicylic acid, and gibberellin) and in different tissues (stems, tender needles, needles, cones, and seeds) using four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper). Then, geometric mean rankings from these algorithms and the RefFinder program were used to comprehensively evaluate RG stability. The results indicated CYP, actin, UBC, and 18S as good choices for studying C. fortunei gene expression. qRT-PCR analysis of the expression patterns of three target genes (CAT and MAPK1/6) further verified that the selected RGs were suitable for gene expression normalization. This study provides an important basis for C. fortunei gene expression standardization and quantification.

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.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Overexpression of a Novel LcKNOX Transcription Factor from Liriodendron chinense Induces Lobed Leaves in Arabidopsis thaliana

Liriodendron chinense is a common ornamental tree that has attractive leaves, which is a valuable trait for use in landscape architecture. In this work, we aimed to identify the potential genes that control and regulate the development of L. chinense leaf lobes. Based on the transcriptome data for the leaf developmental stages we previously generated, two candidate genes were identified in this study. KNOTTED-LIKE HOMEOBOX(KNOX), encoding homeobox family proteins, play a large role in leaf lobe and leaf complexity regulation. Here, two full length KNOX genes from L. chinense were amplified and named LcKNOX1 and LcKNOX6 according to their sequence similarities with the respective Arabidopsis thaliana KNOX family genes. Overexpression vectors were constructed and subsequently transformed into wild type (WT) A. thaliana. Additionally, LcKNOX6 was expressed in tobacco leaves to examine its subcellular localization, and the 35S::LcKNOX6 transgenic A. thaliana leaf cells were imaged with the use of SEM. The expression of several genes that participate in KNOX gene regulation were validated by quantitative real-time PCR. The results show that LcKNOX1 produces almost the same phenotype as that found in WT A. thaliana. Notably, the LcKNOX6-1 lines presented deep leaf lobes that were similar to L. chinense leaf lobes. Two 35S::LcKNOX6 lines induced an abnormal growth phenotype whose seeds were abortive. In short, these results indicate that the LcKNOX6 gene might affect leaf development in A. thaliana and provide insights into the regulation of L. chinense leaf shaping.

Selection of Suitable Reference Genes in Pinus massoniana Lamb. Under Different Abiotic Stresses for qPCR Normalization

The normalization of data by choosing suitable reference genes is fundamental for obtaining accurate and reliable results in quantitative real-time polymerase chain reaction (qPCR) analyses. In this study, the expression stability of 12 candidate reference genes of Pinus massoniana under different abiotic stresses was evaluated using four statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. The results indicate that the following genes could be used as reference genes under different treatments: Actin 2 (ACT2) and F-box family gene (F-box) for salinity treatment, cyclophilin (CYP) and alpha-tubulin (TUA) for ABA treatment, actin 7 (ACT7) and CYP for drought treatment, actin 1 (ACT1) and ACT7 for cold treatment, ACT1 and CYP for heat treatment, and TUA and ACT2 for the “Total” group. To validate the suitability of the selected reference genes in this study, the Short-Root protein (SHR), Alpha-pinene synthase (APS), and Pyrabactin resistance-like protein (PYL) gene expression patterns were analyzed. The expression patterns had significant biases when the most unstable reference genes were used for normalization, compared with when the optimum reference gene or gene combinations were used for normalization. These results will be beneficial for further studies on gene transcription in early-stage, unlignified seedlings of P. massoniana.