The APETALA2/ethylene-responsive transcription factor (AP2/ERF) family has been extensively investigated because of its significant involvement in plant development, growth, fruit ripening, metabolism, and plant stress responses. To date, there has been little investigation into how the AP2/ERF genes influence flower formation and anthocyanin biosynthesis in Lycoris. Herein, 80 putative LrAP2/ERF transcription factors (TFs) with complete open reading frames (ORFs) were retrieved from the Lycoris transcriptome sequence data, which could be divided into five subfamilies dependent on their complete protein sequences. Furthermore, our findings demonstrated that genes belonging to the same subfamily had structural similarities and conserved motifs. LrAP2/ERF genes were analyzed for playing an important role in plant growth, water deprivation, and flower formation by means of gene ontology (GO) enrichment analysis. The expression pattern of the LrAP2/ERF genes differed across tissues and might be important for Lycoris growth and flower development. In response to methyl jasmonate (MeJA) exposure and drought stress, the expression of each LrAP2/ERF gene varied across tissues and time. Moreover, a total of 20 anthocyanin components were characterized using ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis, and pelargonidin-3-O-glucoside-5-O-arabinoside was identified as the major anthocyanin aglycone responsible for the coloration of the red petals in Lycoris. In addition, we mapped the relationships between genes and metabolites and found that LrAP2/ERF16 is strongly linked to pelargonidin accumulation in Lycoris petals. These findings provide the basic conceptual groundwork for future research into the molecular underpinnings and regulation mechanisms of AP2/ERF TFs in anthocyanin accumulation and Lycoris floral development.

Lycoris radiata, belonging to the Amaryllidaceae family, is a well-known Chinese traditional medicinal plant and susceptible to many stresses. WRKY proteins are one of the largest families of transcription factors (TFs) in plants and play significant functions in regulating physiological metabolisms and abiotic stress responses. The WRKY TF family has been identified and investigated in many medicinal plants, but its members and functions are not identified in L. radiata. In this study, a total of 31 L. radiata WRKY (LrWRKY) genes were identified based on the transcriptome-sequencing data. Next, the LrWRKYs were divided into three major clades (Group I–III) based on the WRKY domains. A motif analysis showed the members within same group shared a similar motif component, indicating a conservational function. Furthermore, subcellular localization analysis exhibited that most LrWRKYs were localized in the nucleus. The expression pattern of the LrWRKY genes differed across tissues and might be important for Lycoris growth and flower development. There were large differences among the LrWRKYs based on the transcriptional levels under drought stress and MeJA treatments. Moreover, a total of 18 anthocyanin components were characterized using an ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis and pelargonidin-3-O-glucoside-5-O-arabinoside as well as cyanidin-3-O-sambubioside were identified as the major anthocyanin aglycones responsible for the coloration of the red petals in L. radiata. We further established a gene-to-metabolite correlation network and identified LrWRKY3 and LrWRKY27 significant association with the accumulation of pelargonidin-3-O-glucoside-5-O-arabinoside in the Lycoris red petals. These results provide an important theoretical basis for further exploring the molecular basis and regulatory mechanism of WRKY TFs in anthocyanin biosynthesis and in response to drought stress and MeJA treatment.

The genus Lycoris (about 20 species) includes important medicinal and ornamental plants. Due to the similar morphological features and insufficient genomic resources, germplasm identification and molecular phylogeny analysis are very limited. Here, we sequenced the complete chloroplast genomes of L. chinensis, L. anhuiensis, and L. aurea; they have very similar morphological traits that make it difficult to identify. The full length of their cp genomes was nearly 158k bp with the same guanine-cytosine content of 37.8%. A total of 137 genes were annotated, including 87 protein-coding genes, 42 tRNAs, and eight rRNAs. A comparative analysis revealed the conservation in sequence size, GC content, and gene content. Some variations were observed in repeat structures, gene expansion on the IR-SC (Inverted Repeat-Single-Copy) boundary regions. Together with the cpSSR (chloroplast simple sequence repeats), these genetic variations are useful to develop molecular markers for germplasm identification. Phylogenetic analysis showed that seven Lycoris species were clustered into a monophyletic group, and closed to Narcissus in Amaryllidaceae. L. chinensis, L. anhuiensis, and L. longituba were clustered together, suggesting that they were very likely to be derived from one species, and had the same ancestor with L. squamigera. Our results provided information on the study of genetic diversity, origins or relatedness of native species, and the identification of cultivars.

Taxol is one of the most effective anticancer drugs in the world that is widely used in the treatments of breast, lung and ovarian cancer. The elucidation of the taxol biosynthetic pathway is the key to solve the problem of taxol supply. So far, the taxol biosynthetic pathway has been reported to require an estimated 20 steps of enzymatic reactions, and sixteen enzymes involved in the taxol pathway have been well characterized, including a novel taxane-10β-hydroxylase (T10βOH) and a newly putative β-phenylalanyl-CoA ligase (PCL). Moreover, the source and formation of the taxane core and the details of the downstream synthetic pathway have been basically depicted, while the modification of the core taxane skeleton has not been fully reported, mainly concerning the developments from diol intermediates to 2-debenzoyltaxane. The acylation reaction mediated by specialized Taxus BAHD family acyltransferases (ACTs) is recognized as one of the most important steps in the modification of core taxane skeleton that contribute to the increase of taxol yield. Recently, the influence of acylation on the functional and structural diversity of taxanes has also been continuously revealed. This review summarizes the latest research advances of the taxol biosynthetic pathway and systematically discusses the acylation reactions supported by Taxus ACTs. The underlying mechanism could improve the understanding of taxol biosynthesis, and provide a theoretical basis for the mass production of taxol.

Lycoris is a summer bulbous flower that commonly needs to go through a long period of vegetative growth for 3 to 5 years before flowering. Plant flowering is regulated by a complex genetic network. Compared with most perennial flowers, knowledge on the molecular mechanism responsible for floral transition in bulbous flowers is lacking, and only a few genes that regulate flowering have been identified with few reports on the floral transition in Lycoris. In this study, we identified many differentially expressed genes (DEGs) and transcription factors (TFs) by RNA-Seq in L. chinensis bulbs of different ages, including one- to four-year-old nonflowering bulbs and four-year-old flowering bulbs. Some DEGs were enriched in Gene Ontology (GO) terms between the three- and four-year-old bulbs, and there most genes were enriched in terms of metabolic process and catalytic activity. In the four-year old bulbs, most of the DEGs that may be involved in flowering were classified under the GO term biological process, which was a totally different result from the vegetative bulbs. Some DEGs between flowering and nonflowering bulbs were enriched in plant hormone signal transduction, including the hormones auxin, cytokinin, abscisic acid, and ethylene, but no DEGs were enriched in the gibberellin pathway. Auxin is the main endogenous phytohormone involved in bulb growth and development, but cytokinin, abscisic acid, and ethylene were shown to increase in flowering bulbs. In addition, energy-metabolism-related genes maintain a high expression level in large bulbs, and some positive regulators (SPL, COL, and AP1) and early flowering genes were also shown to be highly expressed in the meristems of flowering bulbs. It suggested that sugar molecules may be the energy source that regulates the signal transduction of flowering by connecting with phytohormone signaling in Lycoris. A total of 1911 TFs were identified and classified into 89 categories, where the top six families with the largest gene numbers were C2H2, NAC, AP2/ERF-ERF, C3H, MYB-related, and WRKY. Most DEGs were in the AP2/ERF-ERF family, and most of them were downregulated in 4-year-old flowering bulbs. A number of families were reported to be involved in plant flowering, including NAC, AP2/ERF, MYB, WRKY, bZIP, MADS, and NF-Y. These results can act as a genetic resource to aid in the explanation of the genetic mechanism responsible for the flowering of Lycoris and other bulbous flowers.

The transition from vegetative to reproductive growth is important for controlling the flowering of Lycoris radiata. However, the genetic control of this complex developmental process remains unclear. In this study, 18 shoot apical meristem (SAM) samples were collected from early-, mid- and late-flowering populations during floral bud differentiation. The histological analysis of paraffin sections showed that the floral bud differentiation could be divided into six stages; the differentiation time of the early group was earlier than that of the middle and late groups, and the late group was the latest. In different populations, some important differential genes affecting the flowering time were identified by transcriptome profiles of floral bud differentiation samples. Weighted gene co-expression network analysis (WGCNA) was performed to enrich the gene co-expression modules of diverse flowering time populations (FT) and floral bud differentiation stages (ST). In the MEyellow module, five core hub genes were identified, including CO14, GI, SPL8, SPL9, and SPL15. The correlation network of hub genes showed that they interact with SPLs, AP2, hormone response factors (auxin, gibberellin, ethylene, and abscisic acid), and several transcription factors (MADS-box transcription factor, bHLH, MYB, and NAC3). It suggests the important role of these genes and the complex molecular mechanism of floral bud differentiation and flowering time in L. radiata. These results can preliminarily explain the molecular mechanism of floral bud differentiation and provide new candidate genes for the flowering regulation of Lycoris.