Skip to main content

Lipidomic and transcriptomic profiles of glycerophospholipid metabolism during Hemerocallis citrina Baroni flowering



Hemerocallis citrina Baroni (daylily) is a horticultural ornamental plant and vegetable with various applications as a raw material in traditional Chinese medicine and as a flavouring agent. Daylily contains many functional substances and is rich in lecithin, which is mostly composed of glycerophospholipids. To study the comprehensive dynamic changes in glycerophospholipid during daylily flowering and the underlying signalling mechanisms, we performed comprehensive, time-resolved lipidomic and transcriptomic analyses of ‘Datong Huanghua 6’ daylily.


Labelling with PKH67 fluorescent antibodies clearly and effectively helped visualise lipid changes in daylily, while relative conductivity and malonaldehyde content detection revealed that the early stages of flowering were controllable processes; however, differences became non-significant after 18 h, indicating cellular damage. In addition, phospholipase D (PLD) and lipoxygenase (LOX) activities increased throughout the flowering process, suggesting that lipid hydrolysis and oxidation had intensified. Lipidomics identified 558 lipids that changed during flowering, with the most different lipids found 12 h before and 12 h after flowering. Transcriptome analysis identified 13 key functional genes and enzymes in the glycerophospholipid metabolic pathway. The two-way orthogonal partial least squares analysis showed that diacylglycerol diphosphate phosphatase correlated strongly and positively with phosphatidic acid (PA)(22:0/18:2), PA(34:2), PA(34:4), and diacylglycerol(18:2/21:0) but negatively with phospholipase C. In addition, ethanolamine phosphotransferase gene and phospholipid-N-methyltransferase gene correlated positively with phosphatidylethanolamine (PE)(16:0/18:2), PE(16:0/18:3), PE(33:2), and lysophosphatidylcholine (16:0) but negatively with PE(34:1).


Overall, this study elucidated changes in the glycerophospholipid metabolism pathway during the daylily flowering process, as well as characteristic genes, thus providing a basis for future studies of glycerophospholipids and signal transduction in daylilies.

Peer Review reports


Hemerocallis citrina Baroni, also known as daylily, is a perennial monocotyledon that is used as an ornamental flower in China [1], edible vegetable, ingredient in traditional Chinese medicine, and a flavouring agent [2]. Furthermore, daylilies have considerable economic value and potential applications [3] due to the presence of various functional components, including amides, polyphenols, anthraquinone, flavonols, naphthol [4], alkaloids [5], rutin, and zeaxanthin [3]. Daylily is also rich in lecithin, which is mostly composed of glycerophospholipids [6]. Pharmacological studies have revealed that daylilies can exert antidepressant, anti-inflammatory, and antitumour effects and have explored the main functional factors and signalling pathways underlying these effects [7,8,9]. For instance, phenolic extracts of daylily have been shown to increase the concentration of neurotransmitters (5-hydroxy tryptamine, dopamine, and norepinephrine) in the hippocampus and frontal cortex via the nuclear factor kappa B signalling pathway, which is inhibited by antidepressants [7, 8]. As daylilies bloom for just 24 h from bud opening to flower tissue wilting, they are often used as a model to study the circadian rhythm of flowers [10, 11], programmed cell death, and flower senescence [12, 13]. Yang et al. [14] studied the pectin features and cell wall changes in daylily tepals during flower opening and senescence and revealed that remodelling of the cell wall pectin is necessary for the occurrence of senescence. Qing et al. [1] provided the first chromosome-level genome of Hemerocallis citrina Borani that provides new insights into rutin biosynthesis and lack of colchicine. Twenty-three key circadian clock genes were identified, which related to sensitivity to light signal input and gating, and these genes might closely relate to Flower opening time (FOT) in Hemerocallis [10].

During plant senescence, membrane deterioration is commonly associated with the progressive decrease in the membrane phospholipid content. Changes in phospholipase D (PLD) activity have been observed in many physiological processes, including senescence [15]. For instance, the loss of membrane integrity during petal senescence in Tradescantia sp. has been associated with a significant decrease in phospholipid content [16]. In addition, plant senescence is characterised by the degradation of cell membranes, and one of the most characteristic features of membrane degradation is the massive decline in phospholipid content [17]. Biochemical changes in cell membranes during senescence include a simultaneous decrease in the contents of all classes of phospholipids and an increase in that of neutral lipids, and the senescence of carnation flowers is accompanied by increased lipid peroxidation and fatty acid de-esterification [18]. Studies on the senescence of roses revealed that the decline in phosphatidylcholine (PC) content with age was due to its lower biosynthetic capacity in older petals and not due to its increased degradation. The decrease in PC content with age paralleled the overall decline in the content of phospholipids, and cytidine-diphosphocholine (CDP-choline) diacylglycerol phosphorylcholine phosphotransferase (CDP-choline phosphotransferase), which is involved in PC synthesis, might be related to such changes [19]. Triacylglycerol content, along with the content of 60 other lipases, has been reported to increase during natural senescence of Arabidopsis leaves [20]. Bax inhibitor-1 is regarded as a cell death suppressor that can interact with other molecules to alter lipid dynamics [21]. In Arabidopsis, leaf senescence-related 1 may be involved in the transport of specific lipid signalling molecules to regulate senescence in leaves [22].

Omics-based technologies, such as transcriptomics, proteomics, and metabolomics, have recently been used to reveal biological mechanisms at the molecular level. In particular, differentially expressed proteins at different stages of daylily flower senescence have been detected using proteomics [23], while transcriptomics allowed the identification of 23 key genes related to flowering rhythm in Hemerocallis sp. [10]. Lipidomics based on analytical chemistry and statistical analysis have been used to comprehensively study the genes involved in lipid pathways and explain their physiological significance in lipid metabolism through large-scale analysis and quantitative studies [24]. In addition, lipidomics using liquid chromatography combined with tandem mass spectrometry (LC–MS/MS) has enabled the large-scale determination of lipid levels and their associated lipid molecular species [25]. Transcriptomics and lipidomics have also been combined to analyse the characteristics of glycerides during flower development in Arabidopsis [26]. However, few studies have combined transcriptomics and lipidomics to better understand glycerophospholipid biosynthesis and metabolism in daylilies.

Traditionally, research on daylily has only been focused on one or a few substances and genes at any given time point. Moreover, glycerophospholipid-associated signals and molecular events related to senescence during the flowering of daylilies remain unclear. In the present study, we have provided a comprehensive and dynamic overview of lipid metabolism in daylilies by combining lipidomics and transcriptomics to analyse the changes in the glycerophospholipid metabolism pathway during the flowering process. In addition, we identified characteristic genes and lipid categories to provide a basis for future studies on daylilies and their glycerophospholipid contents.


Fluorescence localisation of membrane lipids in petals

First, we visualised the membrane lipids of daylily petals using PKH67 immunofluorescent antibodies (1:500; Sigma Aldrich, St. Louis, MI, USA). The tissue morphology and fluorescence levels of membrane lipids in the cells of the petals changed significantly during development from 24 h before flowering to 30 h after flowering (Fig. 1). The membrane lipids in the cells of the flower tissues were evenly distributed 24 h before flowering; however, the fluorescence intensity increased at 0 h when the tepals were slightly opened and was the highest around the vascular bundles in petals. When the petals opened, membrane lipid fluorescence decreased by varying degrees and became weaker around the vascular bundle, suggesting that membrane lipid levels were reduced to a certain extent. The gap between petal tissues was enlarged, and the intercellular space was evident at 24 h after flowering, indicating that the tissue had been damaged during senescence. As expected, almost no fluorescence was observed when the tissue was observed without immunofluorescence labelling 12 h before flowering (Fig. 1). However, the combination of bright field imaging and PKH67 immunofluorescence labelling allowed effective visualisation of petal membrane lipids. In summary, the bilayer membrane lipids were arranged orderly to perform effective functions before flower opening. In the later stage of flower opening, the molecular structure of membrane lipids was destroyed, tissue cavities appeared, and the membrane function was impaired.

Fig. 1
figure 1

PKH67 fluorescent antibody labelling of daylily petal membrane lipids at different stages of flowering. (− 24 h to + 30 h) Fluorescence images of daylily petals 24 h before flowering to 30 h after flowering. The last two pictures are a bright-field image and a fluorescence image of a daylily petal 12 h before flowering without fluorescence labelling, respectively. White arrows indicate vascular bundles. Red arrows indicate petal tissue cavities. Scale bar = 2 cm for flower images. Scale bar = 100 μm for fluorescence images

Membrane permeability, malondialdehyde (MDA) content, and membrane lipid-related enzyme activity

Next, we measured the relative conductivity and MDA content of daylily samples during flowering. Electrical conductivity, which indicates cell membrane permeability, increased at 0 h when the tepals were slightly open, peaked at 12 h after flowering (17.86%, p < 0.05), and then showed a non-significant decrease during wilting and senescence (Fig. 2A, B). Conversely, MDA content first decreased to 5.24 nmol g−1 fresh weight (FW) 6 h after flowering and then gradually increased. These results suggested that cell membrane lipids were intact and functioned normally from early bud development to tepal opening. Conversely, membrane lipids were damaged, cell permeability increased, and MDA content increased 6 h after flowering due to membrane lipid oxidation; however, the cells could still function normally and resist senescence. At 18 h after flowering, most cell functions were lost, and the increased permeability and membrane lipid peroxidation could not be reversed.

Fig. 2
figure 2

Relative conductivity (A), malondialdehyde (MDA) content (B), and membrane lipid-related enzyme activities (C) during daylily flowering. Different letters indicate significant differences (p < 0.05)

In addition, the activities of lipoxygenase (LOX), PLD, and phospholipid-N-methyltransferase (PLMT) were closely associated with lipid metabolism in the cell membrane (Fig. 2C). Notably, LOX and PLD activities tended to increase during daylily development and senescence, whereas PLMT activity first increased and then decreased. LOX activity was 413.3 U/L 6 h after flowering and then decreased from 18 to 24 h before peaking at 30 h (434.6 U/L). Conversely, PLD activity increased continuously until 6 h before flowering, then decreased slightly at 18 h after flowering, and peaked (1,603.9 U/L) at 30 h after flowering. The activity of both the enzymes decreased approximately 18 h after flowering, possibly due to cellular resistance to lipid degradation and increased enzyme activity during this period. The increase in enzyme activity 24 h after flowering indicated that cells were damaged, reducing their ability to resist senescence, while lipid oxidation increased enzyme activity. Since PLMT plays a key role in phospholipid synthesis during lipid metabolism, the peak in its activity during bud development at 6 h before flowering (186.7 U/L) suggests that phospholipid synthesis is active at this stage, generating large amounts of phospholipids to support membrane integrity and cell function. Collectively, these results demonstrated that various indicators were normal in the early stage of daylily flowering and the cell function was maintained. However, the membrane lipid function was impaired in the later stage of flowering and the cell function was inhibited to varying degrees.

Lipidomic changes during flowering

Generally, daylilies are used for food production 12 h before flowering (− 12 h). Considering the differences observed in lipid metabolites before and after flowering, we detected the lipid groups present in daylily samples at − 24, − 12, 0, and + 12 h of flowering. Lipidomic analysis revealed 558 lipids (Table S1), including 199 glycerolipids, 269 phospholipids, 39 sphingolipids, 50 glycolipids, and one sterol (Table 1). Notably, the greatest number of differences in lipids (68) were observed between − 12 h and + 12 h samples, followed by 62 lipid differences between − 24 h and + 12 h samples, and 21 lipid differences − 24 h and − 12 h samples (Fig. 3A, B). The most common lipid difference between all groups was triacylglycerol(18:3/18:2/18:3). The volcanic map showed that 37 differential lipids were significantly up-regulated and 31 were significantly down-regulated in the + 12 h/-12 h group (P < 0.05); however, 21 species were significantly up-regulated and 20 species were significantly down-regulated in the 0 h/-24 h group (Fig. S1). Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analysis revealed that the different metabolites in each group were enriched in the glycerolipid, glycerophospholipid, and phosphatidylinositol (PI) metabolic pathways (Fig. 3C). These results indicated that these three lipid metabolic pathways play an important role and that the lipid metabolism is highly dynamic and varies during the flowering of daylilies.

Table 1 Lipidomic analytes in daylilies
Fig. 3
figure 3

Lipidome differences and KEGG enrichment analysis. A Statistical analysis of lipid differences between each group. B Heatmap of the lipids that differed between the − 12 h and + 12 h groups. C KEGG pathway enrichment analysis of the different lipids

Transcriptomic changes during flowering

A total of 80.51 Gb clean reads were detected in the daylily transcriptome at − 24, − 12, 0, and + 12 h of flowering, with an average of 6.05–7.04 Gb per sample, 93.9–94.43% Q30 bases, and 46.49% GC content. The reads from each sample were compared to the reference genome to evaluate the dynamic changes during daylily flowering, with an alignment rate of 91.77–92.13%. Subsequent protein-coding gene expression analysis revealed six groups containing 1591, 7238, 10,108, 7092, 9307, and 7795 differentially expressed genes (DEGs) (Fig. 4A). The number of common differential genes in the comparison group which combined by -24, -12, 0 and 12 h were 107 (− 12/ − 24), 390 (+ 12/ − 12), 739 (+ 12/ − 24), 579 (+ 12/0), 574 (0/ − 12), 473 (0/ − 24) respectively, 116 of which were common among the six groups (Fig. 4B). The comparison between the − 24 h and + 12 h groups revealed the greatest number of DEGs, with 5252 upregulated and 4856 downregulated genes (Fig. 4A); the volcanic map is shown in Fig. 4C. KEGG enrichment analysis revealed that most DEGs were related to carbohydrate metabolism, followed by genetic information processing-translation, and lipid metabolism (Fig. 4D). These results suggested strong transcriptional changes during daylily flowering, with lipid metabolism playing an important role.

Fig. 4
figure 4

Transcriptome analysis of differentially expressed genes (DEGs) during daylily flowering. A DEGs in each comparison group. B Venn diagram of DEGs in each group. C Up (red)- and down (green)-regulated DEGs between the − 12 h/ + 12 h groups. Gray indicates no difference. D DEG distribution in metabolic pathways

Next, we screened all DEGs involved in glycerophospholipid metabolism to analyse their significance in each group. Twenty-two common DEGs were found between the two 24 h groups (Table S2) that were divided into 13 categories: glycerol-3-phosphate dehydrogenase 1, PLMT, diacylglycerol kinase (DGK), phospholipase A (PLA), phosphoethanolamine N-methyltransferase (PEAMT), PLD, lysophospholipase A, glycerol-3-phosphate acyltransferase, lysophospholipid acyltransferase, choline kinase (CK), triacylglycerol lipase SDP1-like, diacylglycerol diphosphate phosphatase (DPP), and glycerophosphodiester phosphodiesterase (GDE). The enzymes encoded by these genes play important roles in glycerophospholipid metabolism: PLMT, PEAMT, and CK are rate-limiting enzymes involved in lecithin synthesis, while PLD and DGK can degrade phospholipids to produce signal molecules such as phosphatidic acid (PA) and diacylglycerol (DAG). Six of these screened genes were detected using quantitative real-time PCR (qRT-PCR), and their relative expression was consistent with the results of transcriptome analysis, confirming the reliability of the transcriptome data (Fig. S2). These data indicated that glycerophospholipid metabolism is regulated by many genes at different levels during daylily flowering.

Combined analysis of the glycerophospholipid metabolism pathway

As the transcriptome analysis revealed that the DEGs screened during daylily flowering played important roles in the glycerophospholipid metabolic pathway, we next identified DEGs of the same type in the transcriptome with a correlation coefficient > 0.8 and plotted the glycerophospholipid metabolic pathway (Fig. 5A). Simultaneously, we screened different lipids in the glycerophospholipid metabolic pathway using lipidomics. The combined analysis of these two sets of data revealed that ethanolamine kinase genes (EKI) in the phosphatidylethanolamine (PE) synthesis pathway were initially upregulated during daylily flowering; however, the key genes coding for ethanolamine phosphotransferase (EPT) and ethanolamine-phosphate cytidylyltransferase (PCYT) were downregulated, leading to a significant decrease in PE(16:0/18:2), PE(32:1), PE(16:0/18:3), PE(34:3), PE(32:2), PE(33:2), PE(34:4), and PE(35:5) levels. In the PC pathway, the initial synthesis genes were upregulated; however, in the consecutive methylation pathway, both PEAMT and its key synthesis genes diacylglycerol cholinephosphotransferase (CPT) and PLMT were downregulated, resulting in no significant difference in PC(16:0/18:3), PC(34:2), PC(34:3), PC(36:6), PC(36:2), PC(36:4), PC(34:1), and PC(36:5) levels (Table S1). Furthermore, the expression of lysophosphatidylcholine (LPC)(19:3) and LPC(17:0), a hydrolysate of phosphatidylcholine, increased due to the degradation of the latter following PLA upregulation, while phosphatidylserine (PS) levels were increased by upregulated phosphatidylserine synthase (Fig. 5B). Thus, highly dynamic changes in glycerophospholipids, which are the main components of the cell membrane, were observed during flowering.

Fig. 5
figure 5

Genes and lipids in glycerophospholipid metabolism. A DEG expression in the glycerophospholipid metabolism pathway (Refer to KEGG pathway database). B Differential lipid expression in the glycerophospholipid metabolism pathway

Due to the upregulation of PLD and phospholipase C (PLC) in the glycerophospholipid pathway, the glycerophospholipid hydrolysates PA and DAG, which are important signalling molecules in plants, were temporarily upregulated. However, transcriptome and lipidome analyses based on the complex network of synthesis and degradation pathways showed that the signalling molecules displayed a downward trend during the flowering process. In addition, phosphatidylinositol kinase gene, phosphatidylinositol phosphate kinase gene, and PLC/D were sequentially downregulated in the dual messenger signalling pathway, which could cause the decrease in DAG(36:3), DAG(32:1), DAG(32:2), and DAG(39:8) levels observed in the lipidome analysis. Collectively, these results indicated that signal transduction is gradually downregulated during the later stages of daylily flowering alongside the weakening of cell structure and tissue function. Consistently, digalactosyl diglyceride (DGDG) and sulfoquinovosyl diacylglycerol (SQDG), which play important roles in the chloroplast membrane, were upregulated due to the increased DGDG synthase (DGD) and SQDG synthase (SQD) expression.

To comprehensively explore glycerophospholipid metabolism during daylily flowering, we used two-way orthogonal partial least squares (O2PLS) for combined transcriptome and lipidomic analysis (Fig. 6). A total of 52 transcripts and 69 (Table S4) metabolites in the glycerophospholipid pathway were used to construct the model (R2X = 0.944, R2Y = 0.957). The 15 most closely related genes and lipids are shown in Fig. 6. DPP (HciG00002298 and HciG00063704) was strongly and positively correlated with PA(22:0/18:2), PA(34:2), PA(34:4), and DAG(18:2/21:0), whereas these lipid signalling molecules were strongly and negatively correlated with PLC (HciG00067428). In addition, EPT (HciG00043670), ethanolamine-phosphate cytidylyltransferase (HciG00072905), and PLMT (HciG0003099, HciG0003101) were positively correlated with PE(16:0/18:2), PE(16:0/18:3), PE(33:2), and LPC 16:0) but negatively correlated with PE(34:1). Collectively, these results indicated that different genes regulate different substances during the flowering process and that their synthesis and catabolic pathways form a highly complex network.

Fig. 6
figure 6

O2PLS loadings plot of metabolites and transcripts involved in daylily flowering. Circles and triangles represent individual transcript and metabolite loading values, respectively


Lipids are important structural and functional membrane components that play key roles in maintaining cell homeostasis, membrane transport, and signal transduction [27]. To date, few studies have been reported on glycerophospholipid metabolism and phosphatidylinositol signalling pathway in plants such as Arabidopsis thaliana [26], soybean [28], blueberry [29], and bell peppers [30], tobacco [31]. However, the lipids in daylilies, particularly the dynamic changes in glycerophospholipid and lipid signalling molecules during the flowering process, have not yet been reported. Glycerophospholipids constitute one of the eight lipid groups and include PE, PC, phosphatidylglycerol (PG), PI, PS, PA, and lysophospholipids [32]. The structure, function, composition, and synthesis pathways of glycerophospholipids vary dramatically during cell division, differentiation, and senescence. PE synthesis proceeds via two pathways: the PS decarboxylation pathway, which accounts for a small proportion of PE synthesis in plants, and the CDP-ethanolamine pathway, which is the main source of PE in plants and whose rate-limiting enzyme is EPT [33]. The Arabidopsis thaliana genome encodes two EPT (or AAPT) genes with dual specificity. Plants with mutations in the aapt1 or aapt2 genes show no-growth phenotypes, while double knockout mutants that are hemizygous for either aapt1 or aapt2 display impaired pollen and seed development, leading to embryonic lethality of the double knockout plants [34]. PC synthesis involves two primary pathways: the CDP-choline pathway and the triple successive methylation pathway [35]. The rate-limiting enzyme of the latter is PEAMT [36]. CDP-diacylglycerol is a hydrolysate of PE and PC that can be used to synthesise PG and PI, which are also substrates for PS synthesis [33]. As lipid synthesis and degradation are regulated through a large and complex cellular network, further studies are required to fully elucidate the underlying regulatory mechanisms.

Senescence and adverse damage in plant cells are closely related to membrane lipid peroxidation induced by the accumulation of reactive oxygen species. As MDA is an important product of membrane lipid peroxidation, it can be a useful indicator of the degree of membrane lipid peroxidation and membrane damage [37]. Cell membranes generally exhibit selective permeability; however, when plants undergo aging or are affected by the environment, the electrolytes in the cells experience osmosis, causing an increase in electrical conductivity and ion leakage due to changes in membrane structure and lipid composition [29, 38]. Four lipid-degrading enzymes, PLD, PA phosphatase, lipolysis acylhydrolase, and LOX, are closely associated with plant senescence [39]. Therefore, in the present study, we measured the relative conductivity, MDA content, and PLD and LOX activities, which are closely related to cell membrane characteristics, to detect changes in the cell membrane during daylily flowering. Interestingly, a normal physiological function was detected in daylily cells during the early stages of flowering; however, the cell membrane appeared damaged and cell function was affected during the later stages. Previous studies have demonstrated that PLD activity is positively correlated with the rate of senescence in castor bean leaves and could trigger the calcium signalling pathway to regulate additional cell processes [15]. LOX catalyses lipid oxidation, which produces superoxide radicals capable of inducing lipid peroxidation [39]. In addition, the metabolism of membrane lipids, especially phospholipids, can reduce bilayer volume lipid saturation, thereby affecting membrane conformation, boundary defects, content leakage, and ultimately reducing function [39]. Our findings, together with those of previous studies, indicate that the membrane lipid-related properties of daylilies change dynamically and substantially before and after flowering. Therefore, it is necessary to study the specific lipids and genes that underlie these changes in daylilies.

Unfortunately, studying the changes in membrane lipids in detail [40] is challenging due to the complexity of the membrane lipid molecules and limitations of the available analytical methods. Here, we used both transcriptomics and lipidomics to identify changes in the genes and lipids during daylily flowering. The findings are consistent with the notion that gene expression precedes material accumulation.

Cell membrane glycerophospholipids are a potent source of synthetic lipid messengers, which are characterised by transient accumulation and rapid turnover [31]. Phospholipases (PLs) are key enzymes that catalyse the initiation of lipid hydrolysis and produce signalling molecules [41]. The three main PLs in the glycerophospholipid pathway are PLD, PLC, and PLA. The results of our study suggested that glycerophospholipid hydrolysis was increased during the later stages of flowering.

PA synthesis and decomposition occur via multiple pathways and are therefore highly dynamic and complex [42]. PA promotes the negative curvature of the membrane bilayer and increases the permeability of membrane vesicles in vitro [43]. Furthermore, high PA levels can disrupt membrane stability, thereby causing the loss of membrane integrity and membrane protein function during injury, senescence, freezing, and drought [42]. In this study, the lipidome analysis showed that PA levels decreased during flowering, whereas DAG, monogalactosyl diglyceride (MGDG), and DGDG were upregulated, likely due to PA dephosphorylation by PA phosphatase to generate DAG, which was converted into MGDG by MGDG synthase [26]. MGDG and DGDGs are unique and abundant structural glycolipids [40] in photosynthetic membranes with the delicate framework of the plasma membrane [44]. Here, transcriptome and lipidome analyses revealed that DGD and SQD were upregulated during flowering while the MGDG, DGDG, and SQDG content was increased, which stabilised the photosynthetic system II complex [29]. Unlike that observed for cytoplasmic membranes, senescence does not significantly change the phase transformation or lipid composition of the thylakoid membranes of chloroplasts [39], which may explain why the MGDG, DGDG, and SQDG levels increased in the present study.

Paul Karl Horan (PKH) lipophilic, long-chain carbocyanine dyes are used to label artificial membranes and biofilms. The aliphatic tails of PKH dyes, such as PKH26 and PKH67, are rapidly inserted into the exposed lipid bilayer, resulting in the formation of strong non-covalent interactions that promote long-term retention and stable fluorescence while generally not affecting cell growth, viability, or proliferation. Consequently, PKH dyes have been widely used to label cells for tracking both in vivo and in vitro [45, 46]; however, PKH dyes have rarely been applied to plant materials. In the present study, we used PKH67 to determine the fluorescence localisation of lipids in daylily petals for the first time. The images of PKH fluorescence staining are very clear, indicating that this method is feasible in the study of daylily membrane lipids.


To determine the fluorescence localisation of lipids, we used PKH67 in daylily petals for the first time. The resulting images indicated that membrane lipids were extensively damaged in the later stages of flowering and affected the membrane, internal structure, and functions of cells. Furthermore, the results showed that MDA content and relative conductivity, which are closely related to membrane lipids, decreased after 18 h of flowering and that the cellular changes were difficult to reverse. Conversely, the activity of PLD and LOX increased during daylily flowering, indicating that hydrolysis and lipid oxidation intensified, which may lead to changes in membrane lipid conformation and membrane damage. Thirteen DEGs related to glycerophospholipids were identified in the daylily samples and played important roles in glycerophospholipid pathways while significantly regulating lipids. Daylily flowering is mainly regulated via the DAG and PA signalling pathways through PC and PE hydrolysis as well as the PI signalling pathway. DPP expression was significantly and positively correlated with the contents of PA(22:0/18:2), PA(34:2), PA(34:4), and DAG (18:2/21:0); however, PLC expression was negatively correlated with their content. EPT and PLMT expression were significantly and positively correlated with the contents of PE(16:0/18:2), PE(16:0/18:3), PE(33:2), and LPC(16:0), but negatively correlated with that of PE(34:1). This study provides novel insights into the glycerophospholipid metabolism of daylilies.

Materials and methods


‘Datong Huanghua 6’ plants were cultivated in the Hemerocallis germplasm nursery at Shanxi Agricultural University (Taigu, China). Daylily buds were sampled 24, 18, 12, and 6 h before flowering ( −), at tepal bloom (0 h), and 6, 12, 18, 24, and 30 h after flowering ( +). A total of 10 samples were collected from − 24 to + 30 h. Samples for transcriptome and lipidome analyses were collected at − 24, − 12, 0, and + 12 h.

Immunofluorescence labelling of membrane lipids

Daylily petal samples collected at different time points were cut into 0.5 cm pieces and fixed in 38% formaldehyde: glacial acetic acid: 70% ethanol fixative (1:1:18, v/v) for over 24 h before embedding them in paraffin. The tissue sections were then dried in an oven at 65 ℃ for 2 h, placed in xylene for 10 min, dewaxed twice, and hydrated in an ethanol series (100%, 95%, and 80%) with purified water for 5 min each [47, 48]. After incubation in 0.5% Triton X-100 [prepared in phosphate-buffered saline (PBS)] for 20 min at 25 ℃, the slides were immersed in PBS (pH 7.4) three times for 5 min each and then sealed with 5% bovine serum albumin at 37 ℃ for 30 min. Once the sealing fluid around the tissues had been aspirated, the slides were incubated in diluted membrane lipid fluorescent antibody PKH67 (1:500, Sigma Aldrich) in a humid box at 37 °C for 1.5 h and then sealed with a plate sealing solution (YZB, BaSO) containing an anti-fluorescence quenching agent [49, 50]. The slides were imaged using a fluorescence microscope (CKX53; Olympus, Tokyo, Japan; 4 × 100).

Determination of membrane permeability

The daylily samples were washed separately, rinsed with deionised water, and punched into 30 discs per group after the surface moisture had dried. The discs were divided into three clean test tubes (n = 10 per tube) containing 15 mL of deionised water. The tubes were sealed with plastic wrap, which was punctured with an anatomical needle, and then air pumped for 10 min to create vacuum. The intake valve was then opened slowly to allow the petals to sink. After 30 min, the initial conductivity was measured using a conductometer (5430R; Eppendorf, Hamburg, Germany) with deionised water as a blank control group. After the tubes had been incubated in boiling water for 10 min and cooled to 25 ℃, the final conductivity was measured. Membrane permeability was calculated as the relative conductivity (%) [29].

Determination of MDA content

The daylily samples (1 g) were added to 5 mL of 10% (w/v) trichloroacetic acid in an ice bath, ground into a homogenate, and centrifuged (10,000 rpm) at 4 ℃ for 20 min. Next, 2 mL of supernatant was added to 2 mL of 0.67% (w/v) thiobarbituric acid, boiled for 20 min, and cooled to 25 ℃. After further centrifugation (10,000 rpm), the absorbance of the supernatant was detected at 450, 532, and 600 nm (UV5200, Shanghai Metash Instruments Co., Ltd., China). Distilled water was used as the control. MDA content was estimated using the following formula: MDA (nmol g−1 FW) = [6.45*(A532 – A600) – 0.56*A450]*5 [29].

Determination of membrane lipid-related enzyme activity

Enzyme-linked immunosorbent assays (ELISAs) were used to detect the activities of LOX, PLD, and PLMT, all of which are related to membrane lipid metabolism. Briefly, 0.1 g of daylily was ground in liquid nitrogen, extracted using 1 mL 80% methanol, and incubated at − 20 °C for 12 h. After centrifugation at 8000 rpm and 4 °C for 1 h, the supernatant was collected and purified using a C-18 solid-phase extraction column balanced with 80% methanol (1 mL). The eluent was successively eluted with 100% methanol (5 mL), 100% ether (5 mL), and 100% methanol (5 mL). After the samples had been vacuum-dried, they were incubated in PBS buffer (1 mL total volume; pH 7.4) at 25 ℃ for 30 min, followed by centrifugation at 8000 rpm and 4 ℃ for 15 min. The supernatant was collected to perform ELISAs according to the manufacturer’s instructions (Shanghai Enzyme-linked Biotechnology, Shanghai, China) [30].

Lipidome analysis

First, 60 mg of each daylily bud was incubated at − 20 ℃ for 2 min with 20 μL 0.01 mg/mL internal standard Lyso PC17:0 (methanol solvent) and 300 μL methanol–water (1:1, v/v), vortexed for 30 s, and crushed using a high-throughput tissue crusher for 2 min at 60 Hz. The samples were then vortexed with 300 μL chloroform for 30 s, extracted using sonication in an ice-water bath for 10 min, and incubated at − 20 °C for 20 min. After centrifugation at 12,000 rpm and 4 °C for 10 min, the lower chloroform layer (200 μL) was pipetted into an LC–MS injection vial. Chloroform–methanol (300 μL; 2:1, v/v) was added to the remaining solution; the mixture was vortexed for 30 s and extracted using ultrasonication for 10 min in an ice-water bath. The extracted samples were left to stand for 20 min at − 20 °C and then centrifuged (12,000 rpm); thereafter, 300 μL of the lower chloroform layer was added to the previous LC–MS vial. The pooled sample was then dried in a vacuum desiccator. Thereafter, the samples were reconstituted in 200 μL isopropanol-methanol (1:1, v/v), vortexed for 30 s, sonicated for 3 min, and transferred to a 1.5 mL centrifuge tube which was left at − 20 °C for 2 h. After centrifugation (12,000 rpm), the supernatant (150 μL) was collected for LC–MS/MS analysis [51]. Quality control samples were prepared by mixing all sample extracts in equal volumes to the same total volume for each sample.

LC–MS/MS analysis was performed using a Dionex U3000 ultra-high performance liquid chromatography (UHPLC) system (Thermo Fisher Scientific, Waltham, MA, USA), equipped with an ACQUITY UPLC BEH C18 column (100*2.1 mm, 1.7 μm; Waters, Milford, MA, USA) at 60 ℃. Mobile phase A was acetonitrile: water (60:40, v/v), mobile phase B was acetonitrile: isopropanol (10:90, v/v), and both mobile phases contained 10 mmol/L ammonium formate. The injection volume was 5 μL, and the flow rate was 0.4 mL/min. MS was performed using positive and negative heated electrospray ionisation (HESI) modes in a Q Exactive instrument (Thermo Fisher Scientific). The HESI source conditions for positive ionisation mode were as follows: heater temperature, 350 °C; sheath gas flow rate, 50 arbitrary units (au); auxiliary gas flow rate, 15 au; sweep gas flow rate, 1 au; spray voltage, 3.8 kV; capillary temperature, 320 °C; S-Lens RF level 75%; MS1 scanning range, m/z 135–2000. For negative ionisation mode, all parameters except spray voltage (3.0 kV) were the same as in the positive ionisation mode. The mass charge ratio (m/z) of lipid molecules and lipid fragments was obtained using the following method: ten fragment profiles (MS2 Scan, HCD) were obtained after each full scan. MS1 had a resolution of 70,000 at m/z 200, and MS2 had a resolution of 17,500 at m/z 200. The raw data exported by the Q Exactive LC–MS/MS was read using LipidSearch software (Thermo Fisher Scientific), and the exact mass numbers of MSn and parent ions were read. Based on the parent ion and multistage MS data of each independent sample, the lipid molecular structure, and the addition mode of positive and negative ions were identified. The search results of each independent sample were aligned according to a certain retention time range, and the results were combined into a single report to sort out the original data matrix [52, 53].

Transcriptome analysis

Total RNA was extracted from the − 24, − 12, 0, and + 12 h daylily samples using the mirVana miRNA Isolation Kit (Ambion, Austin, TX, USA) and stored at − 70 ℃. RNA concentration and quality were detected using a NanoDrop ND 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and an Agilent Bioanalyzer 2100 System (Agilent Technologies, Palo Alto, CA, USA), respectively. A cDNA library was constructed using a TruSeq Stranded mRNA LTSample Prep Kit (Illumina, San Diego, CA, USA). Briefly, mRNA was used as a template to synthesise single-stranded cDNA, which was converted into purified double-stranded cDNA using a two-strand synthesis reaction system. After end repair, 3' A tail addition, and sequencing adapter addition, the templates were amplified using PCR. Following quality inspection, 150 bp paired-end reads were generated using an Illumina HiSeq™ 2500 system. The raw reads generated by high-throughput sequencing were subjected to quality control and linker removal. The remaining high-quality clean reads were compared to the daylily reference genome and annotated for proteins. The fragments per kilobase of transcript per million mapped reads of each gene was calculated and used to quantify the expression levels of the annotated genes. The nbinomTest function was used to calculate the p and fold-change values for the different comparisons. DEGs with p < 0.05 and 2 < fold change < 0.5 were subjected to gene ontology and KEGG enrichment analyses to determine their main functions and metabolic pathways, respectively [54, 55]. Lipidomic and transcriptomic profiles of glycerophospholipid metabolism pathway referenced the KEGG pathway database [56].

qRT-PCR validation

Total RNA was extracted from each sample using a TaKaRa MiniBEST Plant RNA Extraction Kit (Takara, Dalian, China). RNA concentration and purity were detected using a NanoDrop ND 2000 (NanoDrop Technologies) spectrophotometer. cDNA was synthesised in an ice box using a Reverse Transcriptase Kit (PrimeScript™ RT, Takara, Dalian, China) according to the manufacturer’s instructions. All primers were designed using Primer Premier 5.0. Relative gene expression levels were calculated from the melting curve fluorescence signals using the 2−ΔΔCt method, with AP4 as the reference gene [57].

Statistical analyses

All experiments were repeated three times. Statistical analyses were performed using Excel 2019. SAS 8E was used to analyse significance, and Origin2019 and TBtools were used for mapping. The O2PLS analysis was performed to integrate metabolome and transcriptome data online ( One-way analysis of variance with Duncan's test was conducted to evaluate significant differences between the means of multiple groups (p < 0.05).

Availability of data and materials

The lipidomics datasets generated and analysed during the current study are available in the metabolights repository, with the website:





Choline kinase




Differentially expressed genes


DGDG synthase


Digalactosyl diglyceride


Diacylglycerol kinase


Diacylglycerol diphosphate phosphatase


Enzyme-linked immunosorbent assays


Ethanolamine phosphotransferase


Fresh weight


Heated electrospray ionisation


Kyoto Encyclopaedia of Genes and Genomes


Liquid chromatography combined with tandem mass spectrometry








Monogalactosyl diglyceride


Two-way orthogonal partial least squares


Phosphatidic acid


Phosphate-buffered saline






Phosphoethanolamine N-methyltransferase






Paul Karl Horan


Phospholipase A


Phospholipase C


Phospholipase D








Quantitative real-time PCR


SQDG synthase


Sulfoquinovosyl diacylglycerol


Ultra-high performance liquid chromatography


  1. Qing Z, Liu J, Yi X, Liu X, Hu G, Lao J, et al. The chromosome-level Hemerocallis citrina Borani genome provides new insights into the rutin biosynthesis and the lack of colchicine. Hortic Res. 2021;8:89.

    Article  CAS  Google Scholar 

  2. Cui H, Zhang Y, Shi X, Gong F, Xiong X, Kang X, et al. The numerical classification and grading standards of daylily (Hemerocallis) flower color. PLoS ONE. 2019;14:e0216460.

    Article  CAS  Google Scholar 

  3. Hsu YW, Tsai CF, Chen WK, Ho YC, Lu FJ. Determination of lutein and zeaxanthin and antioxidant capacity of supercritical carbon dioxide extract from daylily (Hemerocallis disticha). Food Chem. 2011;129:1813–8.

    Article  CAS  Google Scholar 

  4. Liu J, Zhong X, Jiang Y, Yu L, Huang X, Dong Z, et al. Systematic identification metabolites of Hemerocallis citrina Borani by high-performance liquid chromatography/quadrupole-time-of-flight mass spectrometry combined with a screening method. J Pharm Biomed Anal. 2020;186:113314.

    Article  CAS  Google Scholar 

  5. Matsumoto T, Nakamura S, Nakashima S, Ohta T, Yano M, Tsujihata J, et al. γ-Lactam alkaloids from the flower buds of daylily. J Nat Med. 2016;70:376–83.

    Article  CAS  Google Scholar 

  6. Hirashima T, Toyoshima M, Moriyama T, Sato N. Evolution of the phosphatidylcholine biosynthesis pathways in green algae: combinatorial diversity of methyltransferases. J Mol Evol. 2018;86(1):68–76.

    Article  CAS  Google Scholar 

  7. Xu P, Wang KZ, Lu C, Dong LM, Le Zhai J, Liao YH, et al. Antidepressant-like effects and cognitive enhancement of the total phenols extract of Hemerocallis citrina Baroni in chronic unpredictable mild stress rats and its related mechanism. J Ethnopharmacol. 2016;194:819–26.

    Article  CAS  Google Scholar 

  8. Li CF, Chen XQ, Chen SM, Chen XM, Geng D, Liu Q, et al. Evaluation of the toxicological properties and anti-inflammatory mechanism of Hemerocallis citrina in LPS-induced depressive-like mice. Biomed Pharmacother. 2017;91:167–73.

    Article  Google Scholar 

  9. Cichewicz RH, Zhang Y, Seeram NP, Nair MG. Inhibition of human tumor cell proliferation by novel anthraquinones from daylilies. Life sci. 2004;74(14):1791–9.

    Article  CAS  Google Scholar 

  10. Ren Y, Gao Y, Zhang Q. Morning and evening alarm of the circadian clock for flower opening times in Hemerocallis. Plant Sci. 2021;311:110992.

    Article  CAS  Google Scholar 

  11. Nitta K, Yasumoto AA, Yahara T. Variation of flower opening and closing times in F1 and F2 hybrids of daylily (Hemerocallis fulva; Hemerocallidaceae) and nightlily (H. citrina). Am J Bot. 2010;97:261–7.

    Article  Google Scholar 

  12. Guerrero C, de la Calle M, Reid MS, Valpuesta V. Analysis of the expression of two thiolprotease genes from daylily (Hemerocallis spp.) during flower senescence. Plant Mol Biol. 1998;36:565–71.

    Article  CAS  Google Scholar 

  13. Tremblay RR, Bourassa S, Nehmé B, Calvo EL. Daylily protein constituents of the pollen and stigma a proteomics approach. J Plant Physiol. 2017;212:1–12.

    Article  CAS  Google Scholar 

  14. Yang Y, Qin N, Huang J, et al. Dynamic changes of pectin epitopes and daylily tepals during flower opening and senescence of Hemerocallis citrina. Sci Hortic. 2021;288(1):110367.

    Article  CAS  Google Scholar 

  15. Ryu SB, Wang X. Expression of phospholipase D during castor bean leaf senescence. Plant Physiol. 1995;108:713–9.

    Article  CAS  Google Scholar 

  16. Suttle JC, Kende H. Ethylene action and loss of membrane integrity during petal senescence in Tradescantia. Plant Physiol. 1980;65:1067–72.

    Article  CAS  Google Scholar 

  17. Fan L, Zheng S, Wang X. Antisense suppression of phospholipase D alpha retards abscisic acid- and ethylene-promoted senescence of postharvest Arabidopsis leaves. Plant Cell. 1997;9:2183–96.

    Article  CAS  Google Scholar 

  18. Shahri W, Tahir I. Flower senescence-strategies and some associated events. Bot Rev. 2011;77:152–84.

    Article  Google Scholar 

  19. Itzhaki H, Mayak S, Borochov A. Phosphatidylcholine turnover during senescence of rose petals. Plant Physiol Biochem. 1998;36:457–62.

    Article  CAS  Google Scholar 

  20. Troncoso-Ponce MA, Cao X, Yang Z, Ohlrogge JB. Lipid turnover during senescence. Plant Sci. 2013;205–206:13–9.

    Article  CAS  Google Scholar 

  21. Shahri W, Tahir I. Flower senescence: Some molecular aspects. Planta. 2014;239:277–97.

    Article  CAS  Google Scholar 

  22. Feng G, Zhong Y, Zou W. Lipid transporter LSR1 positively regulates leaf senescence in Arabidopsis. Plant Signal Behav. 2022;17:2007328.

    Article  Google Scholar 

  23. Ma G, Shi X, Zou Q, Tian D, An X, Zhu K. iTRAQ-based quantitative proteomic analysis reveals dynamic changes during daylily flower senescence. Planta. 2018;248:859–73.

    Article  CAS  Google Scholar 

  24. Hu T, Zhang JL. Mass-spectrometry-based lipidomics. J Sep Sci. 2018;41:351–72.

    Article  CAS  Google Scholar 

  25. Liang J, Wen F, Liu J. Transcriptomic and lipidomic analysis of an EPA-containing Nannochloropsis sp. PJ12 in response to nitrogen deprivation. Sci Rep. 2019;9:4540.

    Article  CAS  Google Scholar 

  26. Nakamura Y, Teo NZ, Shui G, Chua CH, Cheong WF, Parameswaran S, et al. Transcriptomic and lipidomic profiles of glycerolipids during Arabidopsis flower development. New Phytol. 2014;203:310–22.

    Article  CAS  Google Scholar 

  27. Hou Q, Ufer G, Bartels D. Lipid signalling in plant responses to abiotic stress. Plant Cell Environ. 2016;39:1029–48.

    Article  CAS  Google Scholar 

  28. Zhang G, Ahmad MZ, Chen B, Manan S, Zhang Y, Jin H, Wang X, Zhao J. Lipidomic and transcriptomic profiling of developing nodules reveals the essential roles of active glycolysis and fatty acid and membrane lipid biosynthesis in soybean nodulation. Plant J. 2020;103(4):1351–71.

    Article  CAS  Google Scholar 

  29. Wang Y, Ji S, Dai H, Kong X, Hao J, Wang S, et al. Changes in membrane lipid metabolism accompany pitting in blueberry during refrigeration and subsequent storage at room temperature. Front Plant Sci. 2019;10:829.

    Article  Google Scholar 

  30. Kong X, Wei B, Gao Z, Zhou Y, Shi F, Zhou X, et al. Changes in membrane lipid composition and function accompanying chilling injury in bell peppers. Plant Cell Physiol. 2018;59:167–78.

    Article  CAS  Google Scholar 

  31. Sosa Alderete LG, Flor S, Lucangioli S, Agostini E. Impact of phenol on the glycerophospholipid turnover and potential role of circadian clock in the plant response against this pollutant in tobacco hairy roots. Plant Physiol Biochem. 2020;151:411–20.

    Article  CAS  Google Scholar 

  32. Fahy E, Subramaniam S, Murphy RC, Nishijima M, Raetz CR, Shimizu T, et al. Update of the LIPID MAPS comprehensive classification system for lipids. J Lipid Res. 2009;50(Suppl):S9-14.

    Article  CAS  Google Scholar 

  33. Michaud M, Prinz WA, Jouhet J. Glycerolipid synthesis and lipid trafficking in plant mitochondria. FEBS J. 2017;284:376–90.

    Article  CAS  Google Scholar 

  34. Liu Y, Wang G, Wang X. Role of aminoalcoholphosphotransferases 1 and 2 in phospholipid homeostasis in Arabidopsis. Plant Cell. 2015;27(5):1512–28.

    Article  CAS  Google Scholar 

  35. Tasseva G, Richard L, Zachowski A. Regulation of phosphatidylcholine biosynthesis under salt stress involves choline kinases in Arabidopsis thaliana. FEBS Lett. 2004;566:115–20.

    Article  CAS  Google Scholar 

  36. Hirashima T, Toyoshima M, Moriyama T, Nakamura Y, Sato N. Characterization of phosphoethanolamine-N-methyltransferases in green algae. Biochem Biophys Res Commun. 2017;488:141–6.

    Article  CAS  Google Scholar 

  37. Su LQ, Li JG, Xue H, Wang XF. Super absorbent polymer seed coatings promote seed germination and seedling growth of Caragana korshinskii in drought. J Zhejiang Univ Sci B. 2017;18:696–706.

    Article  CAS  Google Scholar 

  38. Elsherbiny EA, Dawood DH, Safwat NA. Antifungal action and induction of resistance by β-aminobutyric acid against Penicillium digitatum to control green mold in orange fruit. Pestic Biochem Physiol. 2021;171:104721.

    Article  CAS  Google Scholar 

  39. Thompson JE, Froese CD, Madey E, Smith MD, Hong Y. Lipid metabolism during plant senescence. Prog Lipid Res. 1998;37:119–41.

    Article  CAS  Google Scholar 

  40. Welti R, Li W, Li M, Sang Y, Biesiada H, Zhou HE, et al. Profiling membrane lipids in plant stress responses. Role of phospholipase D alpha in freezing-induced lipid changes in Arabidopsis. J Biol Chem. 2002;277:31994–2002.

    Article  CAS  Google Scholar 

  41. Wang C, Wang X. A novel phospholipase D of Arabidopsis that is activated by oleic acid and associated with the plasma membrane. Plant Physiol. 2001;127:1102–12.

    Article  CAS  Google Scholar 

  42. Testerink C, Munnik T. Molecular, cellular, and physiological responses to phosphatidic acid formation in plants. J Exp Bot. 2011;62:2349–61.

    Article  CAS  Google Scholar 

  43. Chen M, Liu R, Huang X, Du Z, Heng S, Zeng W. Characterization of low temperature-induced plasma membrane lipidome remodeling combined with gene expression analysis reveals mechanisms that regulate membrane lipid desaturation in Carica papaya. Sci Hortic. 2020;272:109505.

    Article  CAS  Google Scholar 

  44. Wan H, Liu H, Zhang J, Lyu Y, Li Z, He Y, et al. Lipidomic and transcriptomic analysis reveals reallocation of carbon flux from cuticular wax into plastid membrane lipids in a glossy “Newhall” navel orange mutant. Hortic Res. 2020;7:41.

    Article  CAS  Google Scholar 

  45. PužarDominkuš P, Stenovec M, Sitar S, Lasič E, Zorec R, Plemenitaš A, et al. PKH26 labeling of extracellular vesicles: characterization and cellular internalization of contaminating PKH26 nanoparticles. Biochim Biophys Acta Biomembr. 2018;1860:1350–61.

    Article  CAS  Google Scholar 

  46. Nagyova M, Slovinska L, Blasko J, Grulova I, Kuricova M, Cigankova V, et al. A comparative study of PKH67, DiI, and BrdU labeling techniques for tracing rat mesenchymal stem cells. In Vitro Cell Dev Biol Anim. 2014;50:656–63.

    Article  CAS  Google Scholar 

  47. Moreno-Sanz P, D’Amato E, Nebish A, Costantini L, Grando MS. An optimized histological proceeding to study the female gametophyte development in grapevine. Plant Methods. 2020;16:61.

    Article  CAS  Google Scholar 

  48. Cui Y, Peng Y, Zhang Q, Xia S, Ruan B, Xu Q, et al. Disruption of EARLY LESION LEAF 1, encoding a cytochrome P450 monooxygenase, induces ROS accumulation and cell death in rice. Plant J. 2021;105:942–56.

    Article  CAS  Google Scholar 

  49. Barbier M, Morjani H, Muirhead KA, Ronot X, Boutonnat J. Optimization of PKH67 labeling conditions for proliferation monitoring in daunorubicin-treated leukemic cells. Oncol Rep. 2005;14:553–9.

    Article  CAS  Google Scholar 

  50. Hu W, Ru Z, Zhou Y, Xiao W, Sun R, Zhang S, et al. Lung cancer-derived extracellular vesicles induced myotube atrophy and adipocyte lipolysis via the extracellular IL-6-mediated STAT3 pathway. Biochim Biophys Acta Mol Cell Biol Lipids. 2019;1864:1091–102.

    Article  CAS  Google Scholar 

  51. Kehelpannala C, Rupasinghe TWT, Hennessy T, Bradley D, Ebert B, Roessner U. A comprehensive comparison of four methods for extracting lipids from Arabidopsis tissues. Plant Methods. 2020;16:155.

    Article  CAS  Google Scholar 

  52. Liu J, Li Q, Chen J, Jiang Y. Revealing further insights on chilling injury of postharvest bananas by untargeted lipidomics. Foods. 2020;9:894.

    Article  CAS  Google Scholar 

  53. Wang YS, Yao HY, Xue HW. Lipidomic profiling analysis reveals the dynamics of phospholipid molecules in Arabidopsis thaliana seedling growth. J Integr Plant Biol. 2016;58:890–902.

    Article  CAS  Google Scholar 

  54. Sun S, Xiong XP, Zhu Q, Li YJ, Sun J. Transcriptome sequencing and metabolome analysis reveal genes involved in pigmentation of green-colored cotton fibers. Int J Mol Sci. 2019;20:4838.

    Article  CAS  Google Scholar 

  55. Wang M, Chen L, Liang Z, He X, Liu W, Jiang B, et al. Metabolome and transcriptome analyses reveal chlorophyll and anthocyanin metabolism pathway associated with cucumber fruit skin color. BMC Plant Biol. 2020;20:386.

    Article  CAS  Google Scholar 

  56. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  Google Scholar 

  57. Yu Y, Kou M, Gao Z, Liu Y, Xuan Y, Liu Y, et al. Involvement of phosphatidylserine and triacylglycerol in the response of sweet potato leaves to salt stress. Front Plant Sci. 2019;10:1086.

    Article  Google Scholar 

Download references


We are grateful to Prof.Jin Xu for his help in revising this manuscript. We also thank the Shanghai Luming biological technology co., Ltd. (Shanghai, China) for providing lipidomics services.


This work was supported by the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi [grant number 2020L0693], the Key Research and Development Project of Shanxi Province [grant number 201903D211011], the biological breeding engineering of Shanxi Agricultural University [grant number YZGC122], and Shanxi Province Basic Research Plan [grant number 202103021223384].

Author information

Authors and Affiliations



AG designed and performed the research. YY and NQ collected the samples. FH and KL analysed the transcriptome data. JW and YG analysed the lipidome data. SL and GX advised on the project and participated in the revision of the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Guoming Xing or Sen Li.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with relevant guidelines and regulations. The collection of Hemerocallis citrina Baroni has been permitted by the Hemerocallis germplasm nursery of Shanxi Agricultural University.

Consent for publication

Not applicable.

Competing interests

There have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

 The lipids of lipidomic.

Additional file 2: Figure S1.

Lipidomic differences in daylily.

Additional file 3: Table S2.

Summary of key genes involved in the glycerophospholipid metabolism pathways. 

Additional file 4: Figure S2.

Verification of key genes using qRT-PCR.

Additional file 5:

Table S3. Primers of qRT-PCR.

Additional file 6: Table S4.

 The classes of lipids in lipidomic.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Guo, A., Yang, Y., Wu, J. et al. Lipidomic and transcriptomic profiles of glycerophospholipid metabolism during Hemerocallis citrina Baroni flowering. BMC Plant Biol 23, 50 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Hemerocallis citrina Baroni
  • Lipidomics
  • Transcriptomics
  • Glycerophospholipids