Skip to main content

Transcriptional deciphering of the metabolic pathways associated with the bioactive ingredients of wolfberry species with different quality characteristics



Wolfberry is rich in carotenoids, flavonoids, vitamins, alkaloids, betaines and other bioactive ingredients. For over 2,000 years, wolfberry has been used in China as a medicinal and edible plant resource. Nevertheless, the content of bioactive ingredients varies by cultivars, resulting in uneven quality across wolfberry cultivars and species. To date, research has revealed little about the underlying molecular mechanism of the metabolism of flavonoids, carotenoids, and other bioactive ingredients in wolfberry.


In this context, the transcriptomes of the Lycium barbarum L. cultivar ‘Ningqi No. 1’ and Lycium chinense Miller were compared during the fruit maturity stage using the Illumina NovaSeq 6000 sequencing platform, and subsequently, the changes of the gene expression profiles in two types of wolfberries were analysed. In total, 256,228,924 clean reads were obtained, and 8817 differentially expressed genes (DEGs) were identified, then assembled by Basic Local Alignment Search Tool (BLAST) similarity searches and annotated using Gene Ontology (GO), Clusters of Orthologous Groups of proteins (KOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). By combining these transcriptome data with data from the PubMed database, 36 DEGs related to the metabolism of bioactive ingredients and implicated in the metabolic pathway of carotenoids, flavonoids, terpenoids, alkaloids, vitamins, etc., were identified. In addition, among the 9 differentially expressed transcription factors, LbAPL, LbPHL11 and LbKAN4 have raised concerns. The protein physicochemical properties, structure prediction and phylogenetic analysis indicated that LbAPL and LbPHL11 may be good candidate genes involved in regulating the flavonoid metabolism pathway in wolfberry.


This study provides preliminary evidence for the differences in bioactive ingredient content at the transcription level among different wolfberry species, as well as a research and theoretical basis for the screening, cloning and functional analysis of key genes involved in the metabolism of bioactive ingredients in wolfberry.

Peer Review reports


Wolfberry, a solanaceous perennial shrub, is widely used as a functional food and herbal medicine for its biological and pharmacological activity [1, 2]. It has been reported to play a significant role in lowering blood sugar and serum lipids, in addition to having antiaging, immunomodulatory, and antitumour properties, together with other functions [3, 4], mainly due to its various bioactive ingredients, such as Lycium barbarum polysaccharides (LBPs), flavonoids, carotenoids, betaines, and alkaloids [5]. LBPs, for example, possess antiaging, antiapoptotic, and anti-inflammatory properties [6, 7]. Carotenoids (β-carotene and zeaxanthin) have the potential to induce cancer cell apoptosis and prevent heart disease and stroke [8]. Consequently, there has recently been considerable interest in the bioactive ingredients found in wolfberry.

It is well known that differences in wolfberry species lead to variations in their bioactive ingredient content. The Lycium genus includes approximately 100 species worldwide, the majority of which are found in North and South America, with 7 species and 2 varieties found in Northwest and North China [9]. Lycium barbarum L. and Lycium chinense Miller, two closely related species, are currently the best-selling wolfberries, with nearly 90% of all commercially available wolfberries belonging to the former, which is the only wolfberry species listed in the Chinese Pharmacopoeia. As a L. barbarum cultivar, ‘Ningqi No. 1’ is the most widely cultivated variety in China because of its superior quality, high yield, strong adaptability and other characteristics. In addition, the fruits of ‘Ningqi No. 1’ possess the unique morphological features of having fruits that are fusiform or oblong in shape (Fig. 1) and this cultivar is commonly grown in northwestern China’s Ningxia Hui Autonomous Region and Xinjiang Uyghur Autonomous Region. Unlike the fruits of ‘Ningqi No. 1’, L. chinense fruits are elliptical in shape (Fig. 1) and can be found in Northwest China and other warm and subtropical countries, such as Korea, Japan, and some European countries [3]. Additionally, there is a difference between ‘Ningqi No. 1’ and L. chinense in the flavour of the fruit, with the fruit of the former tasting only sweet, whereas the fruit of the latter tastes sweet with some bitterness. Remarkably, current findings indicate that the polysaccharide content of L. chinense fruit is approximately 12.36 mg/g, which is three times that of ‘Ningqi No. 1’ fruit. However, the content of betaine and vitamin C in the fruit of L. chinense are 4.89 mg/g and 1.38 mg/g, respectively, which are lower by approximately 4.68% and 25.41%, respectively, than the content of these bioactive ingredients in the fruit of ‘Ningqi No. 1’ [10]. Clearly, the fruits of ‘Ningqi No. 1’ and L. chinense differ in appearance and internal quality. The essence of such uneven fruit quality is attributed to differences in gene expression and to genetic and environmental factors.

Fig. 1
figure 1

Fresh fruit appearance of ‘Ningqi No.1’ (A) and L. chinense (B)

Transcriptomics is a technique for studying gene expression and transcriptional regulation at a large scale, and it is an effective research tool for studying the relationship between cell phenotype and gene function. RNA-Seq has been utilized to discover and identify functional candidate genes involved in plant growth, development and secondary metabolite biosynthesis, as it is an affordable and trustworthy method for transcriptome analysis. For example, Wen et al. [11] found some unigenes that are closely related to chlorogenic acid content, revealing the potential mechanism of chlorogenic acid biosynthesis through RNA-seq along with de novo transcriptome analysis in Korla fragrant pear grown in the Xinjiang region of China. Similarly, a comparative transcriptome analysis of Strombosiopsis tetrandra roots and leaves was performed, revealing potential genes that regulate the biosynthesis of isoquinoline alkaloids [12]. In L. barbarum, numerous enzyme genes associated with carotenoid biosynthesis [3] and phenylpropanoid biosynthesis [13] were identified through transcriptome analysis. Previous research on L. chinense has primarily concentrated on basic biological research such as nutritional component extraction and separation [6], pharmacology and medicinal function [14, 15], and investigation of candidate genes participated in the biosynthesis of vital secondary metabolites [3, 13]. However, there are few studies to reveal the variations in bioactive ingredient content at the molecular level among multiple species of wolfberry. More importantly, this information about gene expression is of great significance for studying the molecular mechanisms of genetic evolution as well as further improving wolfberry species. Therefore, this study identified several differentially expressed genes associated with the metabolism of bioactive ingredients within the fruits of ‘Ningqi No. 1’ and L. chinense through RNA-Seq. The ultimate goal was to elucidate the molecular basis for differences in fruit bioactive components between these two types of wolfberries and to provide basic data and candidate genes for improving wolfberry fruit quality.


RNA-seq analysis

To provide a detailed summary of the transcriptome and gene expression differences between ‘Ningqi No. 1’ and L. chinense fruits, 6 cDNA samples from each of the two types (with each sample having three biological replicates) were prepared and sequenced using Illumina NovaSeq 6000. Following a thorough quality assessment and data screening, 256,228,924 clean reads were obtained. The Q20 (nucleotide ratio with a quality value greater than 20 in the reads) and Q30 (nucleotide ratio with a quality value greater than 30 in the reads) percentages of each sample were greater than 97.28% and 92.30%, respectively. The GC content of each sample (percentage of GC nucleotide content in high-quality reads) was 41.86% ~ 42.30%. The clean reads were then aligned with the designated Lycium barbarum reference genome to obtain mapped data for subsequent transcript assembly, expression computation, and so on. Consequently, the mapping rate ranged from 84.80% to 95.45% (Table 1). Considering that the correlation of biological duplication is very important for analysing transcriptome sequencing data, Pearson correlation analysis was performed on the fruit samples of ‘Ningqi No. 1’ and L. chinense during the maturity stage. The correlation coefficient was 0.7 < R2 < 1 (Fig. 2A), manifesting that the 6 samples presented a relatively high degree of homogeneity of genes within each sample. Overall, these findings showed that RNA-seq produced the high-quality data and reflected a high assembly integrity, which demonstrated that the experimental data could be subsequently analysed.

Table 1 Statistics of sequencing data (N1, N2 and N3 are three replicates of the ‘Ningqi No.1’ fruits. C1, C2 and C3 are three replicates of the L. chinense fruits.)
Fig. 2
figure 2

A Individual correlation analyses between 6 transcriptome samples of ‘Ningqi No. 1’ and L. chinense. B Gene annotations of ‘Ningqi No. 1’ and L. chinense

Functional annotation

A total of 23,414 expressed genes had the most remarkable BLAST matches to known proteins in the six public databases, KEGG, GO, NR, Swiss-Prot, Pfam, and KOG (Fig. 2B). Among the six databases, NR had the most annotated genes with 88.52% (20,726), while the Swiss-Prot database had the lowest rate, at 66.46% (15,561).

DEG analysis of ‘Ningqi No. 1’ and L. chinense

To identify DEGs, a comparative transcriptome analysis of ‘Ningqi No. 1’ and L. chinense fruits was performed using ‘Ningqi No. 1’ fruits as controls. A total of 8817 DEGs were found, with 4036 genes upregulated and 4181 genes downregulated (Fig. 3), indicating that the gene expression patterns of ‘Ningqi No. 1’ and L. chinense differed greatly.

Fig. 3
figure 3

Volcano map (A) and hierarchical cluster analysis (B) of DEGs

As the world’s largest source of classification information and gene function, the GO database is grouped into three categories: Cellular components (CC), Biological processes (BP) and Molecular function (MF). To determine the functional distribution of the genes, GO annotation was used to classify the DEGs from ‘Ningqi No. 1’ and L. chinense, with the results revealing that these genes were divided into three predominant functional categories, CC, BP and MF, with 54 subcategories (Fig. 4). Specifically, the CC category was divided into 16 subcategories, with ‘cell’ and ‘cell part’ containing the most DEGs, at 3699 (41.95%), followed by ‘organelle’ with 2788 (31.62%) DEGs. The BP category was further subdivided into 27 subcategories, with ‘cellular process’ and ‘metabolic process’ being the two main subcategories, including 3040 (34.48%) and 2619 (29.70%) DEGs. The MF category was mapped into 11 GO terms, with the most plentiful subcategories being ‘binding’ and ‘catalytic activity’, which had 3127 (35.47%) and 2965 (33.63%) DEGs, respectively. The functional classification of all the DEGs could be crudely understood after the GO annotation. Furthermore, a GO enrichment analysis was performed to investigate the distribution of DEGs to clarify differences in gene function among samples. The 50 GO terms with a significant enrichment of DEGs showed that those belonging to BP were the most numerous, with a total of 22. The highest proportion of the 50 GO terms was accounted for by the sulfur compound metabolic process (2.78%), transferase activity, transferring hexosyl groups (2.52%) and UDP-glycosyltransferase activity (2.41%) (Fig. 5).

Fig. 4
figure 4

GO classification of DEGs

Fig. 5
figure 5

GO enrichment of DEGs

The KEGG tool can organically combine genomic with functional information and systematically analyse the metabolic pathways of each gene product in the cell, so as to more comprehensively understand the functions of these gene products. Therefore, KEGG annotation and enrichment analysis of the DEGs in ‘Ningqi No. 1’ and L. chinense was performed. The results demonstrated that a total of 139 pathways were annotated, with the most significant enrichment in ‘Metabolic pathways’, ‘Biosynthesis of secondary metabolites’, ‘Brassinosteroid biosynthesis’, ‘Pyruvate metabolism’, ‘Ascorbate and aldarate metabolism’, ‘Amino sugar and nucleotide sugar metabolism’ and ‘Alanine, aspartate and glutamate metabolism’. Notably, the 20 pathways that were significantly enriched included ‘sesquiterpenoid and triterpenoid biosynthesis’, ‘folate biosynthesis’, and ‘ascorbate and aldehyde metabolism’. The DEGs enriched in these metabolic pathways were related to the metabolism of bioactive ingredients, implying that there might be many DEGs correlated with the difference in bioactive component accumulation in fruits (Fig. 6).

Fig. 6
figure 6

KEGG analysis of DEGs

Analysis of DEGs related to the metabolism of bioactive ingredients

The DEGs can each be located in 139 specific metabolic pathway branches, including 16 pathways related to bioactive ingredient metabolism, in the KEGG database [16,17,18] (Fig. 7). The metabolic pathway for ascorbate and aldarate had the most annotated DEGs, with 76 DEGs. Additionally, the flavonoid metabolic pathway, which included ‘flavone and flavonol biosynthesis’, ‘flavonoid biosynthesis’, ‘anthocyanin biosynthesis’ and ‘isoflavonoid biosynthesis’, was annotated with a total of 60 DEGs. The terpenoid metabolic pathway, which includes ‘sesquiterpenoid and triterpenoid biosynthesis’, ‘monoterpenoid biosynthesis’ and ‘diterpenoid biosynthesis’, was also highly annotated, with 51 DEGs, indicating that there may be many good candidate genes associated with the metabolism of bioactive ingredients found in fruits, such as flavonoids and terpenoids.

Fig. 7
figure 7

The KEGG pathways related to bioactive ingredients metabolism

According to the KEGG pathway annotation, combined with the information from the NR and PubMed databases, 36 DEGs associated with the bioactive ingredient metabolism were identified (Table 2), and then whose expression characteristics in ‘Ningqi No. 1’ and L. chinense were determined by the fragments per kilobase of transcript per million mapped reads (FPKM) values of the transcriptome data (Fig. 8). L. chinense had 21 genes with high expression and 15 genes with low expression when compared to ‘Ningqi No. 1’. These DEGs included five key carotenoid biosynthesis enzyme genes that separately encode phytoene synthase 2, phytoene desaturase, lycopene epsilon-cyclase, carotenoid cleavage dioxygenase 4, and 15-cis-zeta-carotene isomerase. It is also worth noting that these DEGs were all upregulated in L. chinense. In addition, there were 17 genes encode 8 enzymes associated with the flavonoid biosynthesis pathway, namely, phenylalanine ammonia-lyase, 4-coumarate–CoA ligase, anthocyanidin 3-O-glucosyltransferase 5-like, chalcone synthase, anthocyanidin synthase, anthocyanidin reductase-like, dihydvroflavonol-4-reductase and flavanone 3-hydroxylase-like. Moreover, four enzyme genes involved in vitamin metabolism, four enzyme genes related to terpenoid metabolism, 5 enzyme genes associated with alkaloid metabolism and 1 enzyme gene involved in phytosterol synthesis were identified. Notably, the vast majority of enzyme genes involved in the metabolism of vitamins, terpenoids, and alkaloids were found to be downregulated in L. chinense.

Table 2 DEGs related to bioactive ingredients metabolism of ‘Ningqi No.1’ and L. chinense
Fig. 8
figure 8

Heatmap of 36 DEGs related to metabolism of active ingredients

Analysis of transcription factors related to bioactive ingredient metabolism

Transcription factors can ensure the expression of target genes at a specific time and in a particular space, and they are important elements of networks that regulate metabolism, development, response to biotic and abiotic stresses, and so on. By analysing the transcriptome data of ‘Ningqi No. 1’ and L. chinense, this study hopes to find the transcription factors that regulate the metabolism of active components in wolfberry and provide fundamental data for understanding their mechanism of action. As a result, among the 9 transcription factors that were differentially expressed, the MYB family transcription factor APL (LbAPL), PHL11 isoform X2 (LbPHL11) and transcription factor KAN4 (LbKAN4) were identified. They were downregulated in L. chinense compared to ‘Ningqi No. 1’.

Prediction of the physicochemical properties and structures of the transcription factors LbAPL, LbPHL11 and LbKAN4

The amino acid sequences of LbAPL, LbPHL11 and LbKAN4 were analysed by the NCBI ( BLAST tool. The results showed that the three proteins all shared the SANT superfamily conserved domain, which belongs to the MYB protein superfamily.

ProtParam tool ( was applied to analyse the physicochemical properties of the LbAPL, LbPHL11 and LbKAN4 proteins. The findings revealed that the amino acid numbers of these transcription factors were 238, 280 and 401, the molecular weights were 26,558.29 Da, 31,415.25 Da and 45,287.20 Da, and the theoretical isoelectric points (pIs) were 9.37, 6.11 and 9.16, respectively. The instability indexes were all above 40, manifesting that they were unstable proteins. These transcription factors were hydrophilic proteins because of their grand average of hydropathicity (GRAVY) less than 0 (Table 3). The subcellular localization prediction results showed that the transcription factors LbAPL, LbPHL11 and LbKAN4 were in the nucleus.

Table 3 Physical and chemical characteristics of Lycium barbarum MYB transcription factor family proteins

Based on their amino acid sequences, the secondary and tertiary structures of the LbAPL, LbPHL11 and LbKAN4 proteins were each predicted using Self-Optimized Prediction Method with Alignment (SOPMA) ( and SWISS-MODEL ( software in the ExPASy database. The results indicated that the alpha helix and random coil were the main components of these proteins’ secondary structure and that these components formed the spatial structure of these MYB proteins. Among the proteins, the LbAPL and LbPHL11 proteins showed high similarity in their tertiary structures (Fig. 9).

Fig. 9
figure 9

Tertiary structure prediction of Lycium barbarum MYB transcription factor family proteins

Phylogenetic analysis of LbAPL, LbPHL11 and LbKAN4 proteins

Several MYB transcription factors participate in the regulation of flavonoid biosynthesis [19, 20]. MEGA 5.05 software was used in this study to construct a phylogenetic tree between the LbAPL, LbPHL11, and LbKAN4 proteins and the amino acid sequences of the Arabidopsis thaliana MYB family proteins (Fig. 10). The results showed that LbAPL had close ortholog with Arabidopsis thaliana MYB106 (at3g01140.1) and MYB48 (at3g46130.1). LbPHL11 had close ortholog with MYB27 (at3g53200.1) and MYB125 (at3g60460.1). LbKAN4 had close ortholog with MYB98 (at4g18770.1) and MYB115 (at5g40360.1). Therefore, the transcription factors LbAPL, LbPHL11 and LbKAN4 may have the same function as the Arabidopsis thaliana MYB protein from the same branch. In addition, it has been reported that MYB27 is the repressor of the anthocyanin biosynthesis pathway, acting as a member of the MBW complex and repressing transcription via its C-terminal EAR motif [21]. MYB48 regulates flavonol biosynthesis mainly in cotyledons [22]. Obviously, MYB27 and MYB48 were closely connected with the synthesis of flavonoids. Therefore, it could be hypothesized that LbAPL and LbPHL11 might perform similar function in the regulation of flavonoid biosynthesis in wolfberry.

Fig. 10
figure 10

The phylogenetic tree of LbAPL, LbPHL11 and LbKAN4 transcription factors of Lycium barbarum and other MYB proteins involved in flavonoid metabolism from different species

Analysis of LbAPL, LbPHL11 and LbKAN4 gene expression patterns

The expression of the LbAPL, LbPHL11 and LbKAN4 genes in the roots, stems, leaves and mature fruits of ‘Ningqi No. 1’ and L. chinense was analysed through RT‒qPCR. The results indicated that the level of LbAPL expression gradually decreased in the roots, stems, leaves and mature fruits of ‘Ningqi No. 1’ and L. chinense (Fig. 11A). In ‘Ningqi No. 1’, LbPHL11 expression was abundant in stems, roots and mature fruit but low in leaves. However, the level of LbPHL11 expression in mature fruits was higher than that in other tissues in L. chinense (Fig. 11B). LbKAN4 expression was highest in the stems of both ‘Ningqi No. 1’ and L. chinense, followed by the roots, and was lowest in the mature fruits (Fig. 11C). Except for LbPHL11 in the leaves of ‘Ningqi No. 1’, all tissues of ‘Ningqi No. 1’ expressed higher levels of LbAPL, LbKAN4 and LbPHL11 than L. chinense.

Fig. 11
figure 11

Relative expression levels of three transcription factors in different tissues of ‘Ningqi No.1’ and L. chinense. A LbAPLB LbPHL11. C LbKAN4. Data are means (± SD) of three biological replicates. Asterisks ‘*’ represent statistical differences in the same index between different species, with p < 0.05 being a significant difference (*p < 0.05, **p < 0.01, ***p < 0.001)

Analysis of carotenoid, flavonoid, and isoflavone content in the fruit maturity stage of the two wolfberry types

We further quantitatively profiled the carotenoids, flavonoids and isoflavones between the two types of wolfberries based on the higher number of DEGs involved in carotenoid and flavonoid biosynthesis. The results suggested that the total carotenoid content detected in ‘Ningqi No. 1’ (0.278 mg/g) was lower than that in L. chinense (0.360 mg/g). The total flavonoid content was found to be at least 3.2 times higher in L. chinense than that in ‘Ningqi No. 1’. On the contrary, a higher amount of isoflavone was observed in ‘Ningqi No. 1’, which increased by approximately 86% compared to L. chinense (Fig. 12).

Fig. 12
figure 12

A Total carotenoid content, B Total flavonoid content, C Total isoflavone content. The bars represent mean (n = 3), and the error bars indicate standard deviation (SD). The asterisks (***) reveal the high significant difference at p < 0.001

Discussion and conclusions

Wolfberry, a traditional medicinal herb and food supplement, contains a high concentration of bioactive constituents such as LBPs, carotenoids, betaines, flavonoids and vitamins [23]. Wolfberries have been reported to have significant antiaging effects, to enrich the liver and kidneys, and to participate in immune regulation and antitumour processes [1, 5]. In recent years, RNA-Seq has become an effective tool for studying gene expression in depth and detecting novel transcripts. In this study, 8817 DEGs were obtained during the fruit maturity stage of ‘Ningqi No. 1’ and L. chinense via Illumina sequencing technology. Numerous transcripts encoding putative genes involved in bioactive ingredient metabolism were explored using annotation. To gain a better understanding of these genes’ specific functions, in combination with the annotation information from the NR database and the PubMed database, 36 DEGs related to bioactive ingredient metabolism were identified that participated in carotenoid, flavonoid, terpenoid, alkaloid, vitamin metabolic pathways, etc. Interestingly, three transcription factors, LbAPL, LbPHL11 and LbKAN4, were differentially expressed in this study.

Carotenoids, which are fat-soluble compounds, are the second largest group of metabolites of wolfberry. They are primarily composed of zeaxanthin (83%), β-cryptoxanthin (7%), mutatoxanthin (1.4%), β-carotene (0.9%), and other nutrients [8], which are essential for human health. In the present research, the transcriptome dataset contained 5 DEGs encoding related enzymes in the carotenoid biosynthesis pathway, including phytoene synthase 2 (PSY2), phytoene desaturase (PDS), 15-cis-zeta-carotene isomerase (Z-ISO), lycopene ε-cyclase (LCYE) and carotenoid cleavage dioxygenase 4 (CCD4). The important functions of these enzymes in carotenoid biosynthesis are as follows: the first step in the plant carotenoid biosynthesis pathway is to form phytoene, a two-step condensation reaction catalysed by PSY involving two GGPP molecules. Then PDS performs a two-step desaturation reaction of phytoene, converting it into phytofluene and ζ-carotene. The role of Z-ISO in transforming 9,15,9'-tri-cis-ζ-carotene to 9,9'-di-cis-ζ-carotene is critical to carotenogenesis in the darkness. In carotenoid biosynthesis, the third step consists of lycopene cyclization, in which LCYE attaches an ε ring to lycopene and forms δ-carotene. Most notably, these five genes were upregulated in L. chinense compared with ‘Ningqi No. 1’ (Fig. 13). Presumably, the increased expression of PSY2, PDS, Z-ISO and LCYE resulted in higher levels of α-carotene, α-cryptoxanthin, lutein, β-carotene, β-cryptoxanthin and zeaxanthin in L. chinense than in ‘Ningqi No. 1’, which led to the significantly higher total carotenoid content in L. chinense than that of ‘Ningqi No. 1’. However, downregulation of the carotenoid cleavage dioxygenase 4 (StCCD4) transcript level in transgenic potato has been shown to result in higher accumulation of violaxanthin and lutein than in the wild type, thereby increasing carotenoid levels [24]. In sweet Osmanthus, carotenoid cleavage dioxygenase 4 (OfCCD4) can cleave carotenoids such as β-carotene and zeaxanthin [25]. In addition, CitCCD4, which is found in citrus, can cleave β-cryptoxanthin and zeaxanthin [26, 27]. According to these reports, CCD4 is a key gene that negatively regulates carotenoid accumulation. Intriguingly, CCD4 was significantly upregulated in L. chinense in this study, which could enable the cleavages of β-carotene, β-cryptoxanthin, zeaxanthin and other carotenoids. Consequently, the content of total carotenoids would be reduced in L. chinense. However, since the identified five genes, including four genes that positively regulate carotenoid biosynthesis and only one negatively, all upregulated in L. chinense, the total carotenoid content in L. chinense could be still higher than that in Ningqi No.1. This hypothesis was consistent with the findings of the present study that the amounts of total carotenoids (0.360 mg/g) in L. chinense was approximately1.29 times higher in ‘Ningqi No. 1’ (0.278 mg/g). Obviously, PSY2, PDS, Z-ISO, LCYE and CCD4 were suggested as good candidate genes for carotenoid biosynthesis in wolfberry, which needs to be further verified. In general, the identification of these crucial genes in the carotenoid biosynthesis pathway may help to elucidate the mechanism of carotenoid biosynthesis in ‘Ningqi No. 1’ and L. chinense. It could be of great importance to cultivate elite wolfberry species with high carotenoid content for reinforcing nutritional performances.

Fig. 13
figure 13

Schematic representation of the carotenoid biosynthesis pathway coupled with expression levels determined by RNA-seq and RT-qPCR analyses for each DEG. The enzyme names are abbreviated as follows: PSY, phytoene synthase; PDS, phytoene desaturase; ZISO, ζ-carotene isomerase; ZDS, ζ-carotene desaturase; CRTISO, carotenoid isomerase; LCYE, lycopene ε-cyclase; LCYB, lycopene β-cyclase; BCH, β-carotene hydroxylase; CCDs, carotenoid cleavage dioxygenases

The major component of goji berries are flavonoids, which have beneficial biological and pharmacological properties. Purified flavonoids extracted from goji berries, for example, dramatically improve the expression of six antiaging genes and play a pivotal role in ageing-related Alzheimer’s disease [28]. The present study revealed that a total of seventeen DEGs encoding seven enzymes, including one PAL, five 4CL, two CHS, two F3H, two DFR, two ANS and three 3GT, were involved in flavonoid biosynthesis (Fig. 14). Six downregulated and eleven upregulated unigenes were found among the seventeen DEGs. Recently, an abundance of experimental data has shown that transcriptional regulation affects active ingredient biosynthesis [29]. In particular, MYB transcription factors were closely relevant to control phenylpropanoid metabolism. For instance, 13 full-length cDNA clones of R2R3-MYB transcription factor (TFs) from E. sagittatum (EsMYB) have been isolated and characterized and were found to regulate the flavonoid biosynthetic pathway [30]. Moreover, it has been demonstrated that GmMYB176 was in connection with flavonoid biosynthesis in soybean [31]. In this study, we found 3 MYB family transcription factors, APL (LbAPL), PHL11 isoform X2 (LbPHL11) and transcription factor KAN4 (LbKAN4), all downregulated in L. chinense. Significantly, MYB27, which is closely homologous to LbPHL11, is an anthocyanin repressor [21]. LbPHL11 could thus play a critical role in anthocyanin biosynthesis in wolfberry. As is known, red wolfberry contains a small amount of anthocyanins. Precisely these identified DEGs can preliminarily explain the lower anthocyanin accumulation in ‘Ningqi No. 1’ and L. chinense. In this study, these six downregulated genes, including one PAL, two CHS, one 4CL, one F3H and one 3GT, were identified as potentially important in inhibiting anthocyanin accumulation in L. chinense. Additionally, the other eleven upregulated genes and the downregulated transcription factor LbPHL11 may play a role in inhibiting anthocyanin biosynthesis in ‘Ningqi No. 1’. Of course, the flavonoid pathway is a complex process in which multiple enzyme metabolites are produced, including anthocyanin. Since they shared the same upstream pathway, these seventeen DEGs were also involved in other metabolic pathways, such as those of flavones, isoflavones and flavonols. In this study, the isoflavone contents in L. chinense were lower compared with ‘Ningqi No. 1’. Nevertheless, the total flavonoid levels detected in L. chinense were significantly higher than that of ‘Ningqi No. 1’, which might have been caused by the differential expression of pivotal enzyme genes in the flavonoid pathway, such as PAL, 4CL, CHS. Notably, since PAL encoded by Lba08g01407 was the first enzyme in the phenylpropanoid biosynthesis pathway, the expression level of PAL was thus the basis of flavonoid biosynthesis. Furthermore, PAL expression levels were considerably higher in ‘Ningqi No. 1’ than in L. chinense, making it an important candidate gene for studying the difference in flavonoid and isoflavone content in the two types of wolfberries. Interestingly, the LbAPL identified in this study shared high ortholog with MYB48, and MYB48 was found to regulate the process of flavonol biosynthesis in previous studies [22], which suggests that LbAPL could be a key transcription factor regulating flavonol content in wolfberry. Of course, future research is needed to confirm the role of these transcription factor genes in flavonoid biosynthesis.

Fig. 14
figure 14

Schematic representation of the flavonoid biosynthesis pathway coupled with heatmaps of relevant genes involved in the pathway. The enzyme names are abbreviated as follows: PAL, phenylalanine ammonia lyase; C4H, cinnamic acid 4-hydroxylase; 4CL, 4-coumarate: CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; FNS, flavone synthase; IFS, isoflavone synthase; F3H, flavanone 3-hydroxylase; F3′5’H, flavanone 3’,5’-hydroxylase; F3’H, flavonoid 3’-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol 4-reductase; ANS, anthocyanidin synthase; 3GT, anthocyanidin 3-O-glucosyltransferase

In addition to carotenoids and flavonoids, wolfberry fruits are also rich in various vitamins, terpenoids and alkaloids. Among them, vitamins have numerous functions, such as scavenging free radicals, promoting hepatocyte proliferation, enhancing nonspecific immunity, and improving disease resistance [8, 32]. In addition, it is well known that many terpenoids are effective components of Chinese herbal medicine; they have been identified along with terpenoid glycosides in wolfberry. As one of the primary bioactive ingredients of wolfberry, alkaloids possess antioxidant properties and liver protection, neuroprotection, and antitumour properties [33,34,35]. In our research, four, four and five DEGs involved in vitamin, terpenoid and alkaloid metabolism pathways were identified, respectively. Three of the five DEGs involved in the alkaloid metabolic pathway participate in indole alkaloid biosynthesis, while the other two are involved in quinoline alkaloid biosynthesis. These results can provide basic data to elucidate the biosynthetic pathway and to conduct a dynamic analysis of the vitamin, terpenoid and alkaloid content in wolfberry.

To our knowledge, traditional breeding goals primarily focus on cultivating properties such as resistance to insect attack, early or late maturation, and increased yield. Enhancing functional and nutritional values is becoming a new trend in breeding objectives as people’s living standards improve. David Chagné et al. [36], for instance, demonstrated that the candidate genes Leucoanthocyanidin reductase (LAR1) and hydroxy cinnamate/quinate transferase (HCT/HQT) are likely to increase the content of polyphenolic compounds in apples, potentially facilitating the development of new apple cultivars containing fruits with higher concentrations of a variety of polyphenolic compounds with human health benefits. Furthermore, using de novo transcriptome assembly, researchers investigated the pivotal enzyme genes in the biosynthesis of chlorogenic acid in Korla fragrant pear, providing a theoretical basis for pear species selection and breeding in Xinjiang [11]. Therefore, this study focused on explaining the molecular basis of the difference between the bioactive components in ‘Ningqi No. 1’ and L. chinense fruits, identifying key genes and providing a wealth of transcriptome data as guidance for functional component research, future breeding of excellent strains and quality improvement in wolfberry.

Materials and methods

Plant materials and growth conditions

The major L. barbarum cultivar ‘Ningqi No. 1’ [37] and L. chinense [38] were cultivated in the standardization management field of the Qixin Wolfberry Seedling Professional Cooperatives in Zhongning, Ningxia Hui Autonomous Region, China (37°53′N, 105°72′E). Using the diagonal sampling method, the fruits of ‘Ningqi No. 1’ and L. chinense were harvested at the maturity stage, when the fruits were 1–2 times larger and bright red in colour compared with those at the green fruit stage (Fig. 1). Three biological replicates were collected, each of which included wolfberry samples from at least three diverse plants. The collected samples were grouped into two parts: one part was used for biochemical validation, and the other part was frozen at -80 °C for transcriptome sequencing analysis.

RNA extraction and Illumina sequencing

Total RNA and mRNA extraction for sequencing was performed using Metwell Biotechnology Co., Ltd. (Wuhan, China). To ensure that the RNA met the quality requirements for constructing a sequencing library, the RNA purity, concentration and integrity were measured using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA), the Qubit® RNA Assay Kit in Qubit®2.0 Fluorometer (Life Technologies, CA, USA) and the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 System (Agilent Technologies, CA, USA). The cDNA library was constructed with the cDNA library construction kit (Beijing Genomics Institute, Shenzhen, China) and then sequenced using the Illumina NovaSeq 6000 platform at Wuhan Metwell Biotechnology Co., Ltd.

Transcriptome data processing and de novo assembly

The raw sequencing image data was converted into raw reads by CASAVA base recognition. To get high-quality data, adapter sequences, empty reads, and low-quality reads with unknown base pairs ‘N’ greater than 10% and reads Q <  = 20 were eliminated from the raw paired-end reads through fastp [39], and then calculated the GC content of the clean reads. The Q20 and Q30 values were also generated by FastQC to evaluate the base quality. Then, HISAT2 2.1.0 software was utilized to align the clean read sequences from each sample to the Lycium barbarum reference genome [40], obtaining information about the location of the reference genome or genes as well as the unique sequence features of the sequencing sample [41]. StringTie version 1.3.4d software was used to de novo assemble the RNA-seq data into a transcriptome [42].

Sequence annotation and classification

New transcript information was extracted from the comparison results of spliced transcripts and genome annotation, and then the sequence of new genes was collected from the genome. Using BlastX, the sequences were explored for annotation against the NCBI nonredundant (NR, protein database with a cut-off E-value of 10–5. Blast2GO (version: 2.5.0, parameters: default) was used to retrieve Gene Ontology (GO) terms from the annotation of high scoring BLAST matches against the NCBI NR protein database (E-value ≤ 1.0 × 10–5, Then, the GO categories were sorted with in-house Perl scripts. The Kyoto Encyclopedia of Genes and Genomes pathways (KEGG, were annotated against the KEGG database using Blastall software (version: 2.2.23, parameters: default). In addition, the sequences were annotated by aligning them in the Cluster of Orthologous Groups of proteins (KOG,, and the manual annotation and reviewing of protein sequences were conducted with the Swiss-Prot ( and Protein family databases (Pfam, The Plant Transcription Factor Database (PlantTFDB) and Pln TFDB were used to annotate and classify the transcription factors (TFs) using iTAK (v1.7 a) [43,44,45].

Analysis of differentially expressed genes

Using ‘Ningqi No. 1’ fruits as a control, the transcriptome data from L. chinense fruits were analysed using high-throughput sequencing technology. The number of reads per gene was counted based on alignment results and information about the location of the gene on the Lycium barbarum reference genome. Fragments per kilobase of transcript per million fragments mapped (FPKM) of each gene were calculated as an index of gene expression level according to gene length and read count mapped to the gene. DESeq2 1.22.1 software was used to analyse the DEGs between the two groups [46, 47]. The false discovery rate (FDR) was calculated using p values that had been adjusted via the Benjamini‒Hochberg method. The DEGs were screened using |log2Fold Change|≥ 1.5 and FDR < 0.05. Then, GO function enrichment, KOG annotation and KEGG pathway enrichment of the DEGs were analysed.

Real-time quantitative RT‒PCR (RT‒qPCR) analysis

Real-time quantitative PCR (RT‒qPCR) was utilized to verify the reliability of the RNA-seq results. Five key candidate genes were chosen for their important roles in carotenoid metabolism. The same RNA samples used in the transcriptome analysis were reverse transcribed with a HiScript® IIQ RT SuperMix for qPCR (+ gDNA wiper) kit (Vazyme Biotech Co. Ltd, Nanjing, China). Primer Premier 5.0 software was used to design the primer that was synthesized by Shanghai Sangong Biotechnology Co., LTD. The primers are listed in Table 4. The constitutively expressed gene actin (GenBank: HQ415754.1) from wolfberry was applied to normalize gene expression as an internal control. Real-time PCR was preformed in a 20 μL reaction mixture containing 50 ng of template cDNA, 10 μL of 2 × Universal SYBR Green Fast qPCR Mix (ABclonal Technology Co., Ltd, Wuhan, China), 0.4 μL of each primer and ddH2O. qPCR amplification was carried out on the qTOWER Real-time PCR Instrument (Analytik Jena AG, Germany). The following thermal cycling conditions were used: 95 °C for 3 min and 40 cycles of 95 °C for 5 s and 60 °C for 30 s. The 2−ΔΔCt method was used to calculate the relative expression of DEGs [48]. We performed 3 independent bioreplicates per sample and 3 technical replicates per bioreplicate to ensure reproducibility and reliability of the experiment.

Table 4 Primers for real-time fluorescence quantitative PCR

Determination of carotenoid, flavonoid, and isoflavone content in the fruit maturity stage of the two wolfberry types

The determination of total carotenoid content was executed in accordance with a previously described method [49], with slight adjustments. The ground samples of 0.1 g wolfberry (triplicate) were extracted via ultrasound in 4 mL extraction solution (petroleum ether: acetone = 4:1) under dark conditions for 15 min at ambient temperature. Afterwards, the extracting solution was filtered and standardized to 4 mL through extraction solution, subsequently, extract was added into a 96—well plate. Finally, the absorbance values were recorded at 450 nm wavelength by a microplate reader (SpectraMax ABS plus, Molecular Devices, CA, USA), and carotenoid content was computed.

The NaNO2-AlCl3-NaOH approach based on the biochemical kit (NMKD0120, Norminkoda Biotechnology Co., Ltd. Wuhan, China) was utilized to determine the total flavonoid [50, 51]. The tissue homogenate of approximately 0.1 g wolfberry was added into 1.5 ml extracting solution, and then ultrasonically extracted in a water bath at 60℃ for 30 min. The suspension was centrifuged at 10,000 × g for 10 min at 25℃ to obtain an aliquot of extract. Subsequently, 1 mL of supernatant was placed in a 10 mL flask and combined with 4 mL deionized water, 0.3 mL NaNO2 (5%), and 0.3 mL AlCl3·6H2O (10%). The mixture was incubated at ambient temperature for 5–6 min followed by the utilization of 1 N NaOH. The solution was then brought up to 10 mL with 80% methanol, and its absorbance was measured at 510 nm wavelength using a spectrophotometer microplate reader. The total flavonoid content was expressed in mg/g FW.

An assay kit (NMKD0121, Norminkoda Biotechnology Co., Ltd. Wuhan, China) was utilized to measure the total isoflavone content at a 260 nm wavelength using a microplate reader (SpectraMax ABS plus, Molecular Devices, CA, USA) [52].

Availability of data and materials

The transcriptome raw reads have been deposited as a BioProject under accession: PRJNA962102 ( The materials are available from the corresponding author on reasonable request after the publication of the work.


  1. Lal M, Chandraker SK, Shukla R. Antimicrobial properties of selected plants used in traditional Chinese medicine. In: Functional and preservative properties of phytochemicals. 2020. p. 119–43.

    Chapter  Google Scholar 

  2. Jiang Y, Fang Z, Leonard W, Zhang P. Phenolic compounds in Lycium berry: composition, health benefits and industrial applications. J Funct Foods. 2021;77:104340.

    Article  CAS  Google Scholar 

  3. Wang G, Du X, Ji J, Guan C, Li Z, Josine TL. De novo characterization of the Lyciumchinense Mill. leaf transcriptome and analysis of candidate genes involved in carotenoid biosynthesis. Gene. 2015;555(2):458–63.

    Article  CAS  PubMed  Google Scholar 

  4. Lopatriello A, Previtera R, Pace S, Werner M, Rubino L, Werz O, Taglialatela-Scafati O, Forino M. NMR-based identification of the major bioactive molecules from an Italian cultivar of Lyciumbarbarum. Phytochemistry. 2017;144:52–7.

    Article  CAS  PubMed  Google Scholar 

  5. Donno D, Mellano MG, Cerutti AK, Beccaro GL. Nutraceuticals in alternative and underutilized fruits as functional food ingredients: ancient species for new health needs. In: Alternative and replacement foods. 2018. p. 261–82.

    Chapter  Google Scholar 

  6. Tian X, Liang T, Liu Y, Ding G, Zhang F, Ma Z. Extraction, structural characterization, and biological functions of Lycium barbarum polysaccharides: a review. Biomolecules. 2019;9(9):389.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Fu YW, Peng YF, Huang XD, Yang Y, Huang L, Xi Y, Hu ZF, Lin S, So KF, Ren CR. Lycium barbarum polysaccharide-glycoprotein preventative treatment ameliorates aversive. Neural Regen Res. 2021;16(3):543–9.

    Article  CAS  PubMed  Google Scholar 

  8. Gao Y, Wei Y, Wang Y, Gao F, Chen Z. Lycium barbarum: a traditional Chinese herb and a promising anti-aging agent. Aging Dis. 2017;8(6):778–91.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Yao R, Heinrich M, Weckerle CS. The genus Lycium as food and medicine: a botanical, ethnobotanical and historical review. J Ethnopharmacol. 2018;212:50–66.

    Article  CAS  PubMed  Google Scholar 

  10. Shu, X. Study on nutrients and antioxidant activity in different germplasm wolfberry fruit (Master’s degree thesis, Northwest A&F University), 2016.

  11. Wen H, Wang W, Jiang X, Wu M, Bai H, Wu C, Shen L. Transcriptome analysis to identify candidate genes related to chlorogenic acid biosynthesis during development of Korla fragrant pear in Xinjiang. Food Sci Human Wellness. 2022;11(4):854–64.

    Article  CAS  Google Scholar 

  12. Zhang Y, Kang Y, Xie H, Wang Y, Li Y, Huang J. Comparative transcriptome analysis reveals candidate genes involved in isoquinoline alkaloid biosynthesis in Stephania tetrandra. Planta Med. 2020;86(17):1258–68.

    Article  CAS  PubMed  Google Scholar 

  13. Zhao S, Tuan PA, Li X, Kim YB, Kim H, Park CG, Yang J, Li CH, Park SU. Identification of phenylpropanoid biosynthetic genes and phenylpropanoid accumulation by transcriptome analysis of Lyciumchinense. BMC Genomics. 2013;14(1):802.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Zhang M, Tang X, Wang F, Zhang Q, Zhang Z. Characterization of Lycium barbarum polysaccharide and its effect on human hepatoma cells. Int J Biol Macromol. 2013;61:270–5.

    Article  CAS  PubMed  Google Scholar 

  15. Luo Q, Li J, Cui X, Yan J, Zhao Q, Xiang C. The effect of Lycium barbarum polysaccharides on the male rats’ reproductive system and spermatogenic cell apoptosis exposed to low-dose ionizing irradiation. J Ethnopharmacol. 2014;154(1):249–58.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.

    Article  CAS  PubMed  Google Scholar 

  19. Wada T, Kunihiro A, Tominaga-Wada R. Arabidopsis CAPRICE (MYB) and GLABRA3 (bHLH) control tomato (Solanum lycopersicum) anthocyanin biosynthesis. PLoS ONE. 2014;9(9):e109093.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Wang L, Lu W, Ran L, Dou L, Yao S, Hu J, Fan D, Li C, Luo K. R2R3-MYB transcription factor MYB6 promotes anthocyanin and proanthocyanidin biosynthesis but inhibits secondary cell wall formation in Populus tomentosa. Plant J. 2019;99(4):733–51.

    Article  CAS  PubMed  Google Scholar 

  21. Albert NW, Davies KM, Lewis DH, Zhang H, Montefiori M, Brendolise C, Boase MR, Ngo H, Jameson PE, Schwinn KE. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell. 2014;26(3):962–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Chen Y, Fu M, Liguo W, Li H, Renzhong L, Zhanji L. Identification of co-expressed modules of cotton genes responding to Verticillium dahliae infection by WGCNA. Acta Agron Sin. 2020;46:668–79.

    Article  Google Scholar 

  23. Lu F, Zhai R, Ruan S, Yang X, Alenyorege EA, Wang Y, Ding Y, Ma H. Enhancement of ultrasound on the dynamic decolorization of wolfberry (Lyciumbarbarum) polysaccharides. LWT Food Sci Technol. 2021;145:111384.

    Article  CAS  Google Scholar 

  24. Bruno M, Beyer P, Al-Babili S. The potato carotenoid cleavage dioxygenase 4 catalyzes a single cleavage of beta-ionone ring-containing carotenes and non-epoxidated xanthophylls. Arch Biochem Biophys. 2015;572:126–33.

    Article  CAS  PubMed  Google Scholar 

  25. Han Y, Wu M, Cao L, Yuan W, Dong M, Wang X, Chen W, Shang F. Characterization of OfWRKY3, a transcription factor that positively regulates the carotenoid cleavage dioxygenase gene OfCCD4 in Osmanthus fragrans. Plant Mol Biol. 2016;91(4–5):485–96.

    Article  CAS  Google Scholar 

  26. Ma G, Zhang L, Matsuta A, Matsutani K, Yamawaki K, Yahata M, Wahyudi A, Motohashi R, Kato M. Enzymatic formation of beta-citraurin from beta-cryptoxanthin and Zeaxanthin by carotenoid cleavage dioxygenase4 in the flavedo of citrus fruit. Plant Physiol. 2013;163(2):682–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Rodrigo MJ, Alquezar B, Alos E, Medina V, Carmona L, Bruno M, Al-Babili S, Zacarias L. A novel carotenoid cleavage activity involved in the biosynthesis of Citrus fruit-specific apocarotenoid pigments. J Exp Bot. 2013;64(14):4461–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Liu J, Meng J, Du J, Liu X, Pu Q, Di D, Chen C. Preparative separation of flavonoids from Goji berries by mixed-mode macroporous adsorption resins and effect on Abeta-expressing and anti-aging genes. Molecules. 2020;25(15):3511.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zheng H, Fu X, Shao J, Tang Y, Yu M, Li L, Huang L, Tang K. Transcriptional regulatory network of high-value active ingredients in medicinal plants. Trends Plant Sci. 2023;28(4):429–46.

    Article  PubMed  Google Scholar 

  30. Huang W, Sun W, Lv H, Xiao G, Zeng S, Wang Y. Isolation and molecular characterization of thirteen R2R3-MYB transcription factors from Epimedium sagittatum. Int J Mol Sci. 2012;14(1):594–610.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Yi J, Derynck MR, Li X, Telmer P, Marsolais F, Dhaubhadel S. A single-repeat MYB transcription factor, GmMYB176, regulates CHS8 gene expression and affects isoflavonoid biosynthesis in soybean. Plant J. 2010;62(6):1019–34.

    CAS  PubMed  Google Scholar 

  32. Zhang Y, Guo Y, Luo Y, Du M, Yin X, Xu X, Zhang G. Integrated metabolomics and transcriptome revealed the effect of fermented Lycium barbarum residue promoting Ovis aries immunity. Front Immunol. 2022;13:889436.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhou ZQ, Fan HX, He RR, Xiao J, Tsoi B, Lan KH, Kurihara H, So KF, Yao XS, Gao H. Lycibarbarspermidines A-O, new dicaffeoylspermidine derivatives from wolfberry, with activities against Alzheimer’s disease and oxidation. J Agric Food Chem. 2016;64(11):2223–37.

    Article  CAS  PubMed  Google Scholar 

  34. Zhu PF, Dai Z, Wang B, Wei X, Yu HF, Yan ZR, Zhao XD, Liu YP, Luo XD. The anticancer activities phenolic amides from the stem of Lyciumbarbarum. Nat Prod Bioprospect. 2017;7(6):421–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Yang Y, Gao L, Niu Y, Li X, Liu W, Jiang X, Liu Y, Zhao Q. Kukoamine A protects against NMDA-induced neurotoxicity accompanied with down-regulation of GluN2B-containing NMDA receptors and phosphorylation of PI3K/Akt/GSK-3beta signaling pathway in cultured primary cortical neurons. Neurochem Res. 2020;45(11):2703–11.

    Article  CAS  PubMed  Google Scholar 

  36. Chagné D, Krieger C, Rassam M, Sullivan M, Fraser J, André C, Pindo M, Troggio M, Gardiner SE, Henry RA, Allan AC, McGhie TK, Laing WA. QTL and candidate gene mapping for polyphenolic composition in apple fruit. BMC Plant Biol. 2012;23(12):12.

    Article  Google Scholar 

  37. Shi J, Chen L, Zheng R, Guan C, Wang Y, Liang W, Yang S, Wang L, Gong L, Zheng G, Huang B. Comparative phenotype and microRNAome in developing anthers of wild-type and male-sterile Lyciumbarbarum L. Plant Sci. 2018;274:349–59.

    Article  CAS  PubMed  Google Scholar 

  38. Oh MR, Jung SJ, Chae SW, Park BH, Lee SO. Lyciumchinense Miller fruit extract lowers liver enzyme levels in subjects with mild hepatic dysfunction: a randomized, double-blind, placebo-controlled clinical trial. Eur J Nutr. 2023;62(3):1415–25.

    Article  CAS  PubMed  Google Scholar 

  39. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Cao YL, Li YL, Fan YF, Li Z, Yoshida K, Wang JY, Ma XK, Wang N, Mitsuda N, Kotake T, Ishimizu T, Tsai KC, Niu SC, Zhang D, Sun WH, Luo Q, Zhao JH, Yin Y, Zhang B, Wang JY, Qin K, An W, He J, Dai GL, Wang YJ, Shi ZG, Jiao EN, Wu PJ, Liu X, Liu B, Liao XY, Jiang YT, Yu X, Hao Y, Xu XY, Zou SQ, Li MH, Hsiao YY, Lin YF, Liang CK, Chen YY, Wu WL, Lu HC, Lan SR, Wang ZW, Zhao X, Zhong WY, Yeh CM, Tsai WC, Van de Peer Y, Liu ZJ. Wolfberry genomes and the evolution of Lycium (Solanaceae). Commun Biol. 2021;4(1):671.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Perez-Rodriguez P, Riano-Pachon DM, Correa LG, Rensing SA, Kersten B, Mueller-Roeber B. PlnTFDB: updated content and new features of the plant transcription factor database. Nucleic Acids Res. 2010;38(Database issue):D822-827.

    Article  CAS  PubMed  Google Scholar 

  44. Jin J, Zhang H, Kong L, Gao G, Luo J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 2014;42(1):D1182-1187.

    Article  CAS  PubMed  Google Scholar 

  45. Zheng Y, Jiao C, Sun H, Rosli HG, Pombo MA, Zhang P, Banf M, Dai X, Martin GB, Giovannoni JJ, Zhao PX, Rhee SY, Fei Z. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol Plant. 2016;9(12):1667–70.

    Article  CAS  PubMed  Google Scholar 

  46. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Varet H, Brillet-Gueguen L, Coppee JY, Dillies MA. SARTools: a DESeq2- and EdgeR-based R pipeline for comprehensive differential analysis of RNA-Seq data. PLoS ONE. 2016;11(6):e0157022.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

  49. Zhao X. Optimize the methods of leaching and exacting the carotenoids in Fructus Lyciumbarbarum L. with ultrasound. Ningxia Med J. 2010;32(05):424–5.

    Article  Google Scholar 

  50. Luan A, Zhang W, Yang M, Zhong Z, Wu J, He Y, He J. Unveiling the molecular mechanism involving anthocyanins in pineapple peel discoloration during fruit maturation. Food Chem. 2023;412:135482.

    Article  CAS  PubMed  Google Scholar 

  51. Sun X, Wang Z, Li X, Du S, Lin D, Shao Y. Effects of Yucca schidigera extract on serum biochemical parameters, humoral immune response, and intestinal health in young pigeons. Front Vet Sci. 2023;9:1077555.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Zhang Y, Guo W, Chen L, Shen X, Yang H, Fang Y, Ouyang W, Mai S, Chen H, Chen S, Hao Q, Yuan S, Zhang C, Huang Y, Shan Z, Yang Z, Qiu D, Zhou X, Cao D, Li X, Jiao Y. CRISPR/Cas9-mediated targeted mutagenesis of GmUGT enhanced soybean resistance against leaf-chewing insects through flavonoids biosynthesis. Front Plant Sci. 2022;13:802716.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors would like to thank Wolfberry Industry Development Center. Also, we would like to thank the reviewers for valuable comments.


This work was supported by project ‘differentiation analysis and functional gene mining of Lycium barbarum core germplasm resources’ from Ningxia Hui Autonomous Region Wolfberry Industry Development Center of China; National Natural Science Foundation of China (31560418); Key R&D plan of Ningxia Hui Autonomous Region (2022BBF02010); The Chinese Central Government Guiding Local Funds for Science and Technology Development (2023FRD05032).

Author information

Authors and Affiliations



RZ, SJY, JZZ and JNT: conceived and designed the experiments. XXL and WQF: performed the experiments. WQF and HG: analyzed the data. XXL: wrote the paper. RZ, YR and JS: revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Corresponding authors

Correspondence to Rui Zheng, Sijun Yue or Jianning Tang.

Ethics declarations

Ethics approval and consent to participate

L. barbarum cultivar ‘Ningqi No.1’ and L. chinense were collected from the Qixin Wolfberry Seedling Professional Cooperatives in Zhongning, Ningxia Hui Autonomous Region, China (37°53′N, 105°72′E), with permission from whose researcher, Mr. Jinzhong Zhu. The plants are freely accessible to Jinzhong Zhu with non-commercial resource purpose. The authors comply with relevant institutional, national, and international guidelines and legislation for plant study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, X., Zheng, R., Radani, Y. et al. Transcriptional deciphering of the metabolic pathways associated with the bioactive ingredients of wolfberry species with different quality characteristics. BMC Genomics 24, 658 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: