The precision and reliability of reverse transcription quantitative polymerase chain reaction (RT-qPCR) depend mainly on suitable reference genes; however, reference genes have not yet been identified for Liriodendron chinense (Hemsl.) Sarg. In this study, the expression stability of 15 candidate reference genes, ACT7, ACT97, UBQ1, eIF2, eIF3, HIS, BIG, AGD11, EFG, GAPDH, CYP, RPL25, UBC, RPB1, and TUB, was tested across multiple organs of L. chinense using four algorithms, geNorm, NormFinder, BestKeeper, and RefFinder. To understand the difference between the selected reference genes and the unsuitable candidate reference genes, the expression level of a target gene, LcPAT7, was normalized across various plant samples. ACT97 and eIF3 represented the best combination across all samples tested, while AGD11 and UBQ1 were unsuitable for normalization in this case. In the vegetative organ subset, ACT97, ACT7, and GAPDH showed the highest expression stability. For floral organs, UBC and eIF3 were the most stable reference genes. Unsuitable reference genes underestimated the expression levels of a target gene, LcPAT7. This study identified two reference genes (ACT97 and eIF3) for the precise and reliable normalization of L. chinense RT-qPCR data across various organs. Our work provides an effective framework for quantifying gene expression in L. chinense.

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.

Real-time quantitative PCR (RT-qPCR) has a high sensitivity and strong specificity, and is widely used in the analysis of gene expression. Selecting appropriate internal reference genes is the key to accurately analyzing the expression changes of target genes by RT-qPCR. To find out the most suitable internal reference genes for studying the gene expression in Broussonetia papyrifera under abiotic stresses (including drought, salt, and ZnSO₄ treatments), seven different tissues of B. papyrifera, as well as the roots, stems, and leaves of B. papyrifera under the abiotic stresses were used as test materials, and 15 candidate internal reference genes were screened based on the transcriptome data via RT-qPCR. Then, the expression stability of the candidate genes was comprehensively evaluated through the software geNorm (v3.5), NormFinder (v0.953), BestKeeper (v1.0), and RefFinder. The best internal reference genes and their combinations were screened out according to the analysis results. rRNA and Actin were the best reference genes under drought stress. Under salt stress, DOUB, HSP, NADH, and rRNA were the most stable reference genes. Under heavy metal stress, HSP and NADH were the most suitable reference genes. EIF3 and Actin were the most suitable internal reference genes in the different tissues of B. papyrifera. In addition, HSP, rRNA, NADH, and UBC were the most suitable internal reference genes for the abiotic stresses and the different tissues of B. papyrifera. The expression patterns of DREB and POD were analyzed by using the selected stable and unstable reference genes. This further verified the reliability of the screened internal reference genes. This study lays the foundation for the functional analysis and regulatory mechanism research of genes in B. papyrifera.