Plant biological processes, such as growth and metabolism, hormone signal transduction, and stress responses, are affected by gene transcriptional regulation. As gene expression regulators, transcription factors activate or inhibit target gene transcription by directly binding to downstream promoter elements. DOF (DNA binding with One Finger) is a classic transcription factor family exclusive to plants that is characterized by its single zinc finger structure. With breakthroughs in taxonomic studies of different species in recent years, many DOF members have been reported to play vital roles throughout the plant life cycle. They are not only involved in regulating hormone signals and various biotic or abiotic stress responses but are also reported to regulate many plant biological processes, such as dormancy, tissue differentiation, carbon and nitrogen assimilation, and carbohydrate metabolism. Nevertheless, some outstanding issues remain. This article mainly reviews the origin and evolution, protein structure, and functions of DOF members reported in studies published in many fields to clarify the direction for future research on DOF transcription factors.

The formation and development of bulblets are crucial to the Lilium genus since these processes are closely related to carbohydrate metabolism, especially to starch and sucrose metabolism. However, little is known about the transcriptional regulation of both processes. To gain insight into carbohydrate-related genes involved in bulblet formation and development, we conducted comparative transcriptome profiling of Lilium davidii var. unicolor bulblets at 0 d, 15 d (bulblets emerged) and 35 d (bulblets formed a basic shape with three or four scales) after scale propagation. Analysis of the transcriptome revealed that a total of 52,901 unigenes with an average sequence size of 630 bp were generated. Based on Clusters of Orthologous Groups (COG) analysis, 8% of the sequences were attributed to carbohydrate transport and metabolism. The results of KEGG pathway enrichment analysis showed that starch and sucrose metabolism constituted the predominant pathway among the three library pairs. The starch content in mother scales and bulblets decreased and increased, respectively, with almost the same trend as sucrose content. Gene expression analysis of the key enzymes in starch and sucrose metabolism suggested that sucrose synthase (SuSy) and invertase (INV), mainly hydrolyzing sucrose, presented higher gene expression in mother scales and bulblets at stages of bulblet appearance and enlargement, while sucrose phosphate synthase (SPS) showed higher expression in bulblets at morphogenesis. The enzymes involved in the starch synthetic direction such as ADPG pyrophosphorylase (AGPase), soluble starch synthase (SSS), starch branching enzyme (SBE) and granule-bound starch synthase (GBSS) showed a decreasing trend in mother scales and higher gene expression in bulblets at bulblet appearance and enlargement stages while the enzyme in the cleavage direction, starch de-branching enzyme (SDBE), showed higher gene expression in mother scales than in bulblets. An extensive transcriptome analysis of three bulblet development stages contributes considerable novel information to our understanding of carbohydrate metabolism-related genes in Lilium at the transcriptional level, and demonstrates the fundamentality of carbohydrate metabolism in bulblet emergence and development at the molecular level. This could facilitate further investigation into the molecular mechanisms underlying these processes in lily and other related species.

Lilium spp. is a bulb flower with worldwide distribution and unique underground organs. The lack of an efficient genetic transformation system for Lilium has been an international obstacle. Because existing model plants lack bulbs, bulb-related gene function verification studies cannot be carried out in model plants. Here, two stable and efficient genetic transformation systems based on somatic embryogenesis and adventitious bud regeneration were established in two Lilium species. Transgenic plants and T-DNA insertion lines were confirmed by β-glucuronidase (GUS) assay, polymerase chain reaction (PCR) and Southern blot. After condition optimization, transformation efficiencies were increased to 29.17% and 4% in Lilium pumilum DC. Fisch. and the Lilium longiflorum ‘White Heaven’, respectively. To further verify the validity of these transformation systems and apply the CRISPR/Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated protein 9) technology in Lilium, the LpPDS gene in the two Lilium species was knocked out. Completely albino, pale yellow and albino–green chimeric mutants were observed. Sequence analysis in the transgenic lines revealed various mutation patterns, including base insertion, deletion and substitution. These results verified the feasibility and high efficiency of both transformation systems and the successful application of the CRISPR/Cas9 system to gene editing in Lilium for the first time. Overall, this study lays an important foundation for gene function research and germplasm improvement in Lilium spp.

Plant somatic embryos are widely used in the fields of germplasm conservation, breeding for genetic engineering and artificial seed production. MicroRNAs (miRNAs) play pivotal roles in somatic embryogenesis (SE) regulation. However, their regulatory roles during various stages of SE remain unclear. In this study, six types of embryogenic samples of Lilium pumilum DC. Fisch., including organogenic callus, embryogenic callus induced for 4 weeks, embryogenic callus induced for 6 weeks, globular embryos, torpedo embryos and cotyledon embryos, were prepared for small RNA sequencing. The results revealed a total of 2,378,760 small RNA reads, among which the most common size was 24 nt. Four hundred and fifty-two known miRNAs, belonging to more than 86 families, 57 novel miRNAs and 40 miRNA*s were identified. The 86 known miRNA families were sorted according to an alignment with their homologs across 24 land plants into the following four categories: 23 highly conserved, 4 moderately conserved, 15 less conserved and 44 species-specific miRNAs. Differentially expressed known miRNAs were identified during various stages of SE. Subsequently, the expression levels of 12 differentially expressed miRNAs and 4 targets were validated using qRT-PCR. In addition, six samples were mixed in equal amounts for transcript sequencing, and the sequencing data were used as transcripts for miRNA target prediction. A total of 66,422 unigenes with an average length of 800 bp were assembled from 56,258,974 raw reads. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that 38,004 and 15,497 unigenes were successfully assigned to GO terms and KEGG pathways, respectively. Among the unigenes, 2,182 transcripts were predicted to be targets for 396 known miRNAs. The potential targets of the identified miRNAs were mostly classified into the following GO terms: cell, binding and metabolic process. Enriched KEGG analysis demonstrated that carbohydrate metabolism was the predominant pathway in Lilium SE. Thus, we performed systemic characterization, homology comparisons and profiling of miRNA expression, and we constructed an miRNA-target network during Lilium SE for the first time. Our findings establish a foundation for the further exploration of critical genes and elucidation of SE in Lilium.

Background: Roses are famous ornamental plants worldwide. Floral coloration is one of the most prominent traits in roses and is mainly regulated through the anthocyanin biosynthetic pathway. In this study, we investigated the key genes and metabolites of the anthocyanin biosynthetic pathway involved in color mutation in miniature roses. A comparative metabolome and transcriptome analysis was carried out on the Neptune King rose and its color mutant, Queen rose, at the blooming stage. Neptune King rose has light pink colored petals while Queen rose has deep pink colored petals. Result: A total of 190 flavonoid-related metabolites and 38,551 unique genes were identified. The contents of 45 flavonoid-related metabolites, and the expression of 15 genes participating in the flavonoid pathway, varied significantly between the two cultivars. Seven anthocyanins (cyanidin 3-O-glucosyl-malonylglucoside, cyanidin O-syringic acid, cyanidin 3-O-rutinoside, cyanidin 3-O-galactoside, cyanidin 3-O-glucoside, peonidin 3-O-glucoside chloride, and pelargonidin 3-O-glucoside) were found to be the major metabolites, with higher abundance in the Queen rose. Thirteen anthocyanin biosynthetic related genes showed an upregulation trend in the mutant flower, which may favor the higher levels of anthocyanins in the mutant. Besides, eight TRANSPARENT TESTA 12 genes were found upregulated in Queen rose, probably contributing to a high vacuolar sequestration of anthocyanins. Thirty transcription factors, including two MYB and one bHLH, were differentially expressed between the two cultivars. Conclusions: This study provides important insights into major genes and metabolites of the anthocyanin biosynthetic pathway modulating flower coloration in miniature rose. The results will be conducive for manipulating the anthocyanin pathways in order to engineer novel miniature rose cultivars with specific colors.

Lilium is an important commercial market flower bulb. qRT-PCR is an extremely important technique to track gene expression levels. The requirement of suitable reference genes for normalization has become increasingly significant and exigent. The expression of internal control genes in living organisms varies considerably under different experimental conditions. For economically important Lilium, only a limited number of reference genes applied in qRT-PCR have been reported to date. In this study, the expression stability of 12 candidate genes including α-TUB, β-TUB, ACT, eIF, GAPDH, UBQ, UBC, 18S, 60S, AP4, FP, and RH2, in a diverse set of 29 samples representing different developmental processes, three stress treatments (cold, heat, and salt) and different organs, has been evaluated. For different organs, the combination of ACT, GAPDH, and UBQ is appropriate whereas ACT together with AP4, or ACT along with GAPDH is suitable for normalization of leaves and scales at different developmental stages, respectively. In leaves, scales and roots under stress treatments, FP, ACT and AP4, respectively showed the most stable expression. This study provides a guide for the selection of a reference gene under different experimental conditions, and will benefit future research on more accurate gene expression studies in a wide variety of Lilium genotypes.