Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites

Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites Zahra Nemati, Seyyedeh‌ Sanam Kazemi Shahandashti, Adriana Garibay-Hernández, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-3948859/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 06 May, 2024 Read the published version in BMC Plant Biology → Version 1 posted 8 You are reading this latest preprint version Abstract Background : The increasing demand for saffron metabolites in various commercial industries, including medicine, food, cosmetics, and dyeing, is driven by the discovery of their diverse applications. Saffron, derived from Crocus sativus stigmas, is the most expensive spice, and there is a need to explore additional sources to meet global consumption demands. In this study, we focused on yellow-flowering crocuses and examined their tepals to identify saffron-like compounds. Results : Through metabolomic and transcriptomic approaches, our investigation provides valuable insights into the biosynthesis of compounds in yellow-tepal crocuses that are similar to those found in saffron. The results of our study support the potential use of yellow-tepal crocuses as a source of various crocins (crocetin glycosylated derivatives) and flavonoids. Conclusions : Our findings suggest that yellow-tepal crocuses have the potential to serve as a viable excessive source of some saffron metabolites. The identification of crocins and flavonoids in these crocuses highlights their suitability for meeting the demands of various industries that utilize saffron compounds. Further exploration and utilization of yellow-tepal crocuses could contribute to addressing the growing global demand for saffron-related products. Saffron-like compounds Yellow-tepal crocuses Metabolomics Transcriptomics Crocins Flavonoids Alternative sources Figures Figure 1 Figure 2 Figure 3 Figure 4 Highlights This study uncovers yellow-tepal crocuses as sources of saffron-like compounds, meeting industry demand. Metabolomic and transcriptomic research reveals insights into crocins and flavonoids biosynthesis, offering a promising production avenue. Background The spice saffron consists of vivid crimson threads which are indeed dried stigmas of Crocus sativus . To obtain one kilogram of saffron 150,000-160,000 flowers are required. Literally, the small amount of saffron obtained per plant, along with its manual harvesting, leads to saffron being the most expensive spice in the world. From culinary uses and art crafts to medicine and health, the value of saffron has risen over the course of time. While saffron fulfills different ranges of commercial demands these days, its application is not limited to a spice or colorant [ 1 ]. It has a longstanding reputation in traditional medicine and researchers found biomedical and pharmacological properties of its metabolites. Saffron chemical composition is characterized as being rich in apocarotenoids and comprising kaempferol-derivatives as its sole flavonoid species. The apocarotenoids mainly consist of crocetin and its glucosyl esters, also named crocins (which occur as trans- and cis- isomers), 2,6,6-trimethyl-4-hydroxy-1-carboxyaldehyde-1-cyclohexene (HTCC), picrocrocin (which is glucosylated HTCC), and safranal, which is a product of picrocrocin catabolism [ 2 ]. Crocins and safranal are mainly responsible for saffron’s color and its unique aroma, respectively [ 3 ]. These metabolites also exert pharmacological effects on the prevention and/or improvement of various diseases related to the gastrointestinal, cardiovascular, endocrine, nervous, and immune systems [ 4 ]. Belonging to plant apocarotenoids [ 5 ], crocins are found in crocin-producing plants such as crocuses, Buddleja davidii , and Gardenia jasminoides [ 6 ]. The interaction among general isoprenoid pathways (MEP: methylerythritol phosphate and MVA: mevalonate pathways), carotenoid, and apocarotenoid biosynthetic pathways [ 7 ], results in the presence of a wide variety of crocins in crocin-producing plants. In case of saffron stigmas, these include crocin-I ( t -C-4: trans -crocetin di(β-D-gentibiosyl) ester), with the highest abundance more than 70% of all crocins, crocin-II ( t -C-3: trans -crocetin (β-D-gentibiosyl) (β-D-glucosyl) ester), crocin-III ( t- C-2: trans -crocetin (β-D-gentibiosyl) ester), crocin-IV ( t- C-2: trans -crocetin di(β-D-glucosyl) ester), and crocin-V ( t- C-1: trans -crocetin (β-D-glucosyl) ester) [ 8 , 9 ]. These five forms differ regarding the position and number of glucosyl moieties attached to the crocetin aglycone backbone [ 10 ]. In the general isoprenoid pathways, the precursors of carotenoids are generated via MEP and MVA, followed by the activity of enzymes from the upstream carotenoid biosynthetic pathway such as phytoene synthase (PSY), lycopene β-cyclase (LCYB), and β-carotene hydrolase (BHY). The next stages of saffron apocarotenoid synthesis are accomplished by three key enzymes including the carotenoid cleavage dioxygenase (CCD2), aldehyde dehydrogenase (ALDH), as well as crocetin and HTCC glucosyltransferases (UGTs) [ 11 ]. Being also commercially valuable compounds, kaempferol and its glycosylated derivatives were also identified in crocuses. In fact, kaempferol and its derivatives comprise 90% of the total flavonoid content in tepals of C. sativus [ 12 , 13 ]. Quercetin glycosides are the second most abundant flavonoid class in crocus tepals, comprising 5–10% of the flavonoid compounds, whereas other components like dihydrokaempferol glycosides and naringenin are observed in negligible amounts [ 12 ]. At the initial stage of the flavonoid biosynthetic pathway, phenylalanine is converted into p- coumaryl-CoA by the sequential activity of phenylalanine ammonia lyase, cinnamic acid 4-hydrolase, and p -coumaric acid:CoA ligase. The chalcone synthase catalyzes the first committed step into flavonoid biosynthesis through the condensation of p- coumaroyl-CoA with malonyl-CoA, resulting in the formation of naringenin chalcone, which the chalcone isomerase transforms into naringenin that is further metabolized into multiple flavonoids in a plant lineage-specific manner. Flavonoids are extensively decorated by diverse chemical reactions, including the glycosylation of flavonoid aglycones that is catalyzed by glycosyltransferases [ 14 – 16 ]. Regarding the growing demand for saffron, tepals of C. sativus , which are the major by-products of saffron production, could be an additional source of saffron-like specialized metabolites [ 17 ]. The presence of crocin as well as kaempferol derivatives in tepals of C. sativus was confirmed through multiple analytical approaches such as liquid chromatography coupled to diode array detection (HPLC-DAD) and infrared (IR) spectroscopy [ 18 ]. These studies identified at least six different crocins in C. sativus tepals and the remarkable amount of kaempferol derivatives of 126 mg g -1 dry weight [ 18 ]. Owing to these metabolites’ presence in saffron tepals, they have been suggested as an alternative or supplementary medicine against some diseases [ 19 ] and effectively applied as natural pigments in natural and synthetic substrates [ 20 ]. [ 6 ] C. sativus tepals are not the only extra source of these metabolites. There are about 250 accepted species within the genus Crocus (Iridaceae) around the world [ 21 – 25 ]. They contain a broad spectrum of natural products including flavonoids, coumarins, additional phenolic glycosides, alkaloids, monoterpenoids, sesquiterpenoids, carotenoids, α-hydroxyacids, fatty alcohols, fatty aldehydes, fatty esters, fatty acids, and alkanes [ 26 ]. Within Crocus , compounds similar to those from saffron stigmas may occur in the tepals of different species. Being divided into spring- and autumn-flowering crocuses, the color of their flower’s ranges from white over yellow in spring, to purple and blue-lilac in autumn species [ 6 ]. In general, the presence of carotenoids and flavonoids in flowers leads to petal pigmentation in the range of yellow to orange [ 27 ]. While saffron purple tepals contain crocin and flavonoids, it is possible that the tepals of yellow crocuses also contain saffron-like metabolites due to their yellow color and thus, may be a potential source of these metabolites. Since saffron stigmas and petals are the sources of saffron main metabolites in autumn, yellow-tepal crocuses could increase the yield per year if they contain saffron-like metabolites. In this study, we compared the semi-polar metabolite profiles of the four yellow-tepal Crocus taxa C. chrysanthus , C. flavus subsp. dissectus , C. graveolens , and C. korolkowii via LC-PDA analysis to investigate the presence of compounds similar to those in saffron stigmas. We also provide insights into the biosynthesis of crocins and kaempferol derivatives in yellow-tepal crocuses through their transcriptomic analysis with special focus on the active genes involved in carotenoid and flavonoid metabolism. Results Apocarotenoid and flavonoid profiles of yellow tepals compared to saffron stigmas To investigate the presence of apocarotenoids and flavonoids similar to those of C. sativus stigmas in yellow-tepal crocus plants, the composition of C. sativus stigmas was initially assessed through LC-PDA-MS (Fig. 1 ). We annotated a total of ten major compounds in C. sativus stigmas comprising both, apocarotenoids and flavonoids: (i) safranal, (ii) HTCC (picrocrocin precursor), (iii) picrocrocin, (iv) four trans -crocins including the well-known saffron crocin trans -crocetin-digentibiosyl ester ( t -C-4), (v) one cis -crocin ( c- C-4), and (vi) two major flavonoids, kaempferol 3 -O- sophoroside 7 -O- glucoside (K1) and kaempferol 3 -O- sophoroside (K2). The most abundant crocins in the stigma samples were t- C-3, t- C-4, and c- C-4, while the crocins t- C-2 and t- C-5 were identified in lower abundance (Fig. 1 ). After having established this baseline, we compared via LC-PDA the composition of the stigmas from C. sativus to that of yellow-tepals of C. chrysanthus , C. flavus subsp. dissectus , C. graveolens , and C. korolkowii to determine whether similar compounds were found (Fig. 2 , 3 ). A heatmap shows the metabolite composition of C. sativus stigmas compared to that from tepals of yellow-tepal crocuses (Fig. 3 , Table S1 ). We did not find the monoterpenoids HTCC, picrocrocin, and safranal in any of the yellow-tepal samples. In contrast, the yellow-tepals from C. chrysanthus and C. korolkowii shared the capability of synthesizing trans- crocins that are also found in saffron stigmas. The tepals from these two species displayed the t -C-5 crocin, whereas only C. korolkowii also showed t -C-4 and t -C-2, which are the major crocins found in saffron stigmas (Fig. 1 , [ 34 ]). Interestingly, the yellow-tepals from all the Crocus species displayed additional crocins in their profiles, as confirmed by their UV-Vis spectra (Fig. 2 , S1 , S2). These major crocins detected in the tepal samples did not correspond to the dominant crocins in C. sativus stigmas (Fig. S2). The tepals displayed rich flavonoid profiles that clearly differ from those of saffron stigmas (Fig. 2 , S3). Only C. korolkowii displayed the two kaempferol derivatives that we also found in the saffron stigmas, kaempferol 3-O-sophoroside and kaempferol 3-O-sophoroside 7-O-glucoside (K2 and K1, respectively; Fig. S3). To gain insights into the biosynthetic genes related to saffron-like metabolites in the four yellow-tepal Crocus species, the tepals of C. chrysanthus , C. flavus subsp. dissectus , C. graveolens , C. korolkowii were used for transcriptomic analysis. A total of 177,738,884 clean reads were obtained across all the samples. We compared the assembled sequences to known C. sativus genes involved in apocarotenoid and flavonoid metabolism (Table 1 ). The candidate genes were categorized into three major groups involved saffron apocarotenoid and flavonoid metabolism: carotenoid cleavage dioxygenases (CCDs), uridine diphosphate glucosyltransferases (UGTs), and aldehyde dehydrogenases (ALDHs). A CCD, CCD2, catalyzes the first step of apocarotenoid biosynthesis in saffron, being responsible for the oxidative cleavage of zeaxanthin and the production of HTTC and crocetin dialdehyde [ 2 ]. ALDHs are responsible for the transformation of crocetin dialdehyde into crocetin. The UGTs are responsible for glycosylation of low molecular weight substrates such as crocetin, HTTC, and flavonoid aglycones to produce crocin, picrocrocin, and kaempferol derivatives, respectively. The phylogenetic trees of these gene families are provided, indicating their roles in saffron metabolites pathways (Fig. 4 A). Genes associated with the CCDs family were selected from both C. sativus and other close-related crocus species to encompass a wider range of genetic diversity, as CCD2 is the one characterized/known in C. sativus . Conversely, the selection of genes linked to the UGTs and ALDHs families was limited to those exclusively found in C. sativus , given the substantial abundance of genes from these families within this particular species. Table 1 Genes related to apocarotenoid and kaempferol biosynthesis pathways of C. sativus in yellow-tepal crocuses Name GenBank accession no. C. korolkowii C. chrysantus C. graveolens C. flavus ssp. dissectus E-value/%identitiy E-value/identitiy% E-value/identitiy% E-value/identitiy% C. sativus carotenoid cleavage dioxygenase 2 (CCD2) KJ541749 1e-161 / 70 0.0 / 73 3.00e-116 / 77 1e-46 / 74 C. sativus carotenoid cleavage dioxygenase 2L (CCD2L) KP887110 1e-161 / 66 0.0 / 73 1.00e-115 / 77 2e-46 / 74 C. sativus chromoplast carotenoid cleavage dioxygenase 4b (CCD4b) EU523663.1 No hits found 2e-142/87% 0.0/87% No hits found C. sativus chromoplast carotenoid cleavage dioxygenase 4a (CCD4a) EU523662.1 No hits found 2e-142/87% 0.0/87% No hits found C. sativus carotenoid cleavage dioxygenase 1 OL606625.1 0.0/97% 0.0/96% 0.0/96% 3E-139/97% C. ancyrensis carotenoid cleavage dioxygenase (CCD4c) KP792758.1 0.0/89% 0.0/96% No hits found 3E-115/90% C. ancyrensis carotenoid cleavage dioxygenase (CCD4a/b) KP792757.1 4E-133/98% 0.0/99% 0.0/90% 8E-087/89% C. ancyrensis carotenoid cleavage dioxygenase (CCD1) KP792755.1 0.0/96% 0.0/98% 0.0/98% 9E-146/99% C. sativus carotenoid cleavage dioxygenase 1 (CCD1) MN540633.1 0.0/97% 0.0/96% 0.0/96% 3E-139/97% C. sativus chromoplast carotenoid cleavage dioxygenase 4b (CCD4b) EU523663.1 No hits found 2e-142/87% 0.0/87% No hits found C. sativus UGT709G1 mRNA, complete cds (uridine diphosphate glycosyltransferase) KX385186.1 3e-62 / 67 0.0 / 61 3e-101 / 38 6e-06 / 40 C. sativus glucosyltransferase 2 (GLT2) mRNA AY262037.1 6e-65 / 74 1e-159 / 61 2.00e-49 / 30 7e-46 / 66 C. sativus crocetin glucosyltransferase 74AD1 mRNA MF596166.1 3e-74 / 36 0.0 / 81 2.00e-56 / 31 2e-51 / 71 C. sativus beta-carotene hydroxylase to a 915-base sequence of most likely codons CAC95130 1e-144 / 87 3e-144 / 86 2.00e-13 / 90 0.62 / 50 C. sativus aldehyde dehydrogenase 2B4 (ALDH2B4) mRNA MG672523.1 0.0 / 98 0.0 / 44 1.00e-125 / 96 2.7 / 47 C. sativus aldehyde dehydrogenase ALDH2C4 mRNA, complete MF596160.1 0.0 / 81 0.0 / 97 2.00e-104 / 82 0.024 / 39 C. sativus aldehyde dehydrogenase 3I1 mRNA, complete MF596165.1 0.0 / 96 4e-102 / 58 1.00e-50 / 94 5.8 / 47 C. sativus aldehyde dehydrogenase ALDH5F1 mRNA MF596161.1 0.0 / 96 4e-102 / 38 1.00e-50 / 94 5.8 / 47 C. sativus aldehyde dehydrogenase 6B2 (ALDH6B2) mRNA MG672524.1 0.0 / 97 0.0 / 97 3.00e-93 / 97 1.4 / 35 C. sativus Aldehyde dehydrogenase ALDH7B4 mRNA, complete MF596162.1 0.0 / 99 0.0 / 99 4.00e-171 / 99 2.3 / 41 C. sativus Glucosyltransferase (UGT91P3) MZ190170.1 4e-90 / 62 7e-121 / 86 2.00e-43 / 92 0.68 / 30 C. sativus glucosyltransferase (UGT91P6) MZ190175.1 4e-117 / 79 1e-100 / 45 6.00e-38 / 31 0.010 / 42 C. sativus glucosyltransferase (UGT91K3) MZ190174.1 1e-69 / 72 0.0 / 43 7.00e-40 / 30 1e-07 / 45 C. sativus glucosyltransferase (UGT91K2) MZ190173.1 5e-67 / 73 0.0 / 78 1.00e-36 / 31 2e-08 / 47 C. sativus glucosyltransferase (UGT91P5) MZ190172.1 2e-97 / 64 1e-121 / 82 1.00e-43 / 92 5.0 / 33 glucosyltransferase (UGT91P4) MZ190171.1 2e-134 / 81 1e-106 / 46 3.00e-39 / 31 0.011 / 45 C. sativus glucosyltransferase (Kaempferol) HE793682.1 0.0 / 92 0.0 / 90 5.00e-104 / 91 2e-05 / 29 C. sativus flavonoid glucosyltransferase (GT45) gene FJ194947.1 0.0 / 69 0.0 / 70 2.00e-47 / 28 8e-36 / 81 C. sativus UDP-glucose-dependent flavonoid UGT703B1 KJ381079 7e-92 / 41 0.0 / 89 0.0 / 68 1e-10 / 34 While only a few of these genes were annotated in C. flavus subsp. dissectus with a low identity percentage, meaningful values for all the candidate genes were observed in C. chrysanthus , C. graveolens , and C. korolkowii (Table 1 ). Transcripts of these three species showed a considerable similarity to the candidate saffron genes. When considering the identity percentage of above 70%, the highest number of reads with the highest identity percentage were observed in C. korolkowii. In fact, the transcripts of C. korolkowii were easily aligned to the main genes involved in the apocarotenoid and flavonoid pathways, including CCD2, UGTs and ALDHs. Among the assembled transcripts, we observed potential isoforms of some candidate genes including a flavonoid glucosyltransferase (GT45), UDP-glucose-dependent flavonoid glucosyltransferase (UGT703B1), carotenoid cleavage dioxygenase 2 (CCD2), and beta-carotene hydroxylase in C. chrysanthus , glucosyltransferase 2 (GLT2) in C. chrysanthus, C. graveolens and C. korolkowii , aldehyde dehydrogenase 2B4 (ALDH2B4) and glucosyltransferase (kaempferol) in C. chrysanthus and C. korolkowii (Table S2?). Discussion As more research is conducted on the potential health benefits of saffron value-added metabolites, the demand for these compounds is likely to increase. Finding alternative sources of these metabolites could help ensure a stable supply for the market. While historically the stigmas of some other crocuses like C. cartwrightianus , the known ancestor of saffron [ 35 , 36 ], were used as wild saffron, only C. sativus stigmas contain the most considerable amount of these metabolites. Since previous investigations confirmed the presence of value-added metabolites in saffron tepals, tepals of other crocuses are a potential source worth exploring. In the present work, using the targeted metabolomics and transcriptomics approaches, four yellow-tepal crocus species were investigated: C. flavus subsp. dissectus , C. graveolens , C. korolkowii and C. chrysanthus . By examining the saffron apocarotenoid pathway, we could gain a deeper understanding of the genes involved in producing saffron-like metabolites. The presence of various types of crocins was proved in all four species. However, major crocins of C. sativus ( t -C-2, t- C-3, t- C-4) were only detected in C. chrysanthus and C. korolkowii (Fig. 3 ), which is also supported by our transcriptomic analysis (Table 1 ). We analyzed three key enzymes of the apocarotenoid pathway, CCDs, ALDHs, and UGTs. We observed that CCD2, the enzyme catalyzing the first committed step of crocin biosynthesis, is present in the tepals from all of the four yellow-tepal crocuses. We also identified transcripts for several ALDH homologs including ALDH2B4, ALDH2C4, ALDH3I1, ALDH5F1, ALDH6B2, ALDH7B4 (Fig. 4 A, Table 1 ). The presence of CCD2, ALDH and crocin-related UGTs transcripts was confirmed in tepals of C. flavus subsp. sissectus and C. graveolens , despite no saffron-like crocins being detected in their LC-PDA analysis. The absence of detected stigma-like crocins in LC-PDA analysis prompted us to explore a combination of three different hypotheses to understand the underlying reason. The post-harvest degradation of crocins could be considered as the first explanation for the absence of crocins. Since all transcripts of the crocin pathway genes were observed in C. graveolens , crocins might have become degraded subsequently. The effect of environmental factors on crocins were investigated by HPLC-DAD-MS [ 37 ] where light and temperature were introduced as the main elements that could degrade crocins after only one week. Due to the detection of other types of crocins in all four species, as well as their detection in saffron stigmas, this hypothesis is not likely. We considered the possibility that other compounds and CCD enzymes, rather than those directly involved in the synthesis of saffron-like crocins, may play a more important role in yellow tepal color. Besides crocins, there are various apocarotenoids such as β-cyclocytral and β-ionone which are CCD cleavage products that provide specific aromas and colors ranging from yellow to red in fruits and flowers [ 38 ]. CCD2 belongs to the larger family of CCD enzymes which also includes CCD1, CCD4, CCD7, and CCD8 [ 39 ]. As it is shown in our phylogenetic analysis (Fig. 4 A), candidate genes for CsCCD and CCDs of other crocin-producing plants clustered separately while BdCCD1 clustered in CsCCD. CsCCDs clustered into two groups led by CCD1 and CCD4. The role of CCD1 and CCD4 in providing unique aromas and colors through carotenoid degradation was demonstrated in different plants [ 40 ]. For instance, in Medicago truncatula CCD1 is responsible for yellow color [ 41 , 42 ], while in petals of chrysanthemum, CmCCD4a results in colorless compounds in white flowers [ 43 ]. In contrast, the reddish color in citrus fruits is attributed to CitCCD4 [ 44 ]. Moreover, previous studies validated the presence of CCD1a-, CCD1b- and CCD4a/b-encoding genes in the stigmas of C. sativus and in tepals of C. chrysanthus and C. korolkowii [ 45 ]. CsCCD1 and CsCCD4 convert β-carotene into β-ionone and β-cyclocitral [ 11 ] that play significant roles in the synthesis of aroma and flavor. In our study, CCD1 and CCD4 transcripts were also observed in all the studied yellow-tepal species. The presence of these transcripts could contribute to their yellow color along with the presence of crocins (Fig. 4 ). What is evident is that crocins, although present, are found in low abundance on a dry weight basis compared to saffron stigmas. However, it is essential to consider them as potential commercial sources due to their significantly higher mass compared to saffron stigmas. The third and the most probable hypothesis refers to the glycosylation of plant secondary metabolites that is done by UGTs displaying the well-known PSPG (plant secondary-product-glycosyltransferase) box [ 34 ]. While more than a hundred UGT-encoding genes exist in each plant genome [ 46 ], the saffron-related ones were identified recently [ 2 ] that are shown in phylogenetic tree (Fig. 4 A). While UGT74AD1 (also named GL2) allows the formation of crocins with one and two glucose molecules, UGT91P3 does not use crocetin as a substrate and generate crocins with more than two glucose molecules [ 33 ]. Studies demonstrated the expression patterns of the six identified genes (UGT91K2, UGT91K3, UGT91P3, UGT91P4, UGT91P5 and UGT91P6) of the UGT91 subfamily in several tissues of C. sativus [ 33 ]. While only UGT91P3 showed higher expression levels in the stigma, higher expression levels of UGT91K3 were detected in leaves [ 33 ], suggesting the involvement of the latter in the synthesis of other compounds rather than crocins, which are absent in the leaves. As it is shown in the phylogenetic tree, these genes are grouped together (Fig. 4 A). UGT91K2, UGT91K3 and UGT91P3 had expression patterns similar to UGT74AD1[ 33 ]. Based on the Blast analysis of the crocin-related UGTs, we observed a higher identity percentage among UGT74AD1 (GL2), UGT91K2, UGT91K3, UGT91P3, UGT91P4, UGT91P5 and UGT91P in C. korolkowii (Table 1 ). It could potentially confirm the observed variety of crocins at the metabolome level. In contrast, a lower number of these genes with a high identity percentage were found in C. chrysanthus and the decreasing trend continued in C. graveolens and reached to almost non-significant ones in C. flavus subsp. dissectus . This observation could also suggest that these species are increasingly distantly related to saffron. Our metabolome results provide evidence explaining why we still find crocins in yellow tepals, but these differ significantly from those found in the stigmas. The main reason behind these differences could be attributed to variations in the substrate affinity and catalytic properties of the putative crocetin UGTs. Notably, all the detected crocins in the tepals exhibited higher polarity compared to those present in saffron stigmas, leading to their earlier appearance in the LC chromatogram. This higher polarity suggests a greater degree of glycosylation. Investigating glucose availability in the tissues, as glucose is used for glycosylating the crocetin backbone, could be a compelling aspect to explore further. This might shed light on the discrepancy between tepals and stigmas in terms of crocin composition and abundance. This is in line with another study on some yellow-tepal crocuses showing that the crocins in tepals possess higher glycosylation levels in comparison to the crocins in C. sativus stigmas [ 6 ]. While crocins in saffron stigmas contain up to five and six glucose molecules [ 47 , 48 ], crocins in yellow-tepal crocuses include up to eight sugar molecules [ 6 ]. In fact, these highly glycosylated crocins demonstrate differences in retention times indicating different arrangements of the glucose molecules on the ends of the crocetin [ 49 ]. Moreover, previous studies showed that these crocins are also present in the stigmas of spring crocuses and because of their absence in the stigma of autumn crocuses could be considered a distinguishing factor between spring and autumn crocuses [ 50 ]. The apocarotenoid HTCC is also glycosylated by UGTs to produce picrocrocin, the precursor of safranal. UGT709G1 catalyzes the HTCC glucosyltransferase reaction [ 2 ] Fig. 4 ). Not HTCC, neither picrocrocin nor safranal were detected in yellow-tepals samples (Fig. 3 ). Nevertheless, at the transcriptome level, the UGT709G1 was identified in almost all samples. While the final enzyme in safranal production is still unknown [ 51 , 52 ], we could only check HTCC glucosyltransferase that eventually produces picrocrocin. Considering the unknown safranal pathway, several possibilities could be contemplated to explain the absence of safranal despite the presence of UGT709G1. Furthermore, it has been reported that the presence of picrocrocin and safranal in stigmas is greatly influenced by the developmental stage of the flowers [ 14 ]. Conducting a time-course analysis of tepals could potentially shed light on the ability of crocus tepals to produce these monoterpenoids. Alternatively, the substrates of UGT709G1 might differ in a species- and/or tissue-specific manner, leading to variations in the production of certain compounds. To provide a comprehensive view of the presence of saffron value-added metabolites in yellow-tepal crocuses, the presence of saffron stigma-like kaempferols was also investigated. Our results confirmed the presence in C. korolkowii of two kaempferol derivatives also identified in C. sativus stigmas. Alongside with kaempferol glycosides, which are the sole flavonoids in saffron stigmas, a higher diversity of flavonoids, such as anthocyanins, quercetin- and flavone-related compounds (luteolin, tricin, acacetin, apigenin, and scutellarein) are found in saffron tepals. [ 53 , 54 ]. The detection of multiple compounds with maximum UV-Vis absorption values at 320 nm in all tepal samples, supports the presence of a higher diversity of flavonoids in the yellow-tepal crocuses. Analysis of C. chrysanthus tepals has proved the existence of kaempferol 3-O glucoside, [ 55 ]. Other flavonoids such as luteolin-derivatives in C. corsicus and C. minimus , and tricin in C. heuffelianus and C. korolkowii have been reported [ 56 ]. In our study, we found flavonoid glucosyltransferase transcripts with the highest to the lowest E-value and identity percentage in C. korolkowii , C. chrysanthus , C. graveolens and C. flavus subsp. dissectus , respectively. Ultimately, tepals of C. korolkowii with three proven types of saffron-like crocins and two kaempferol-derivatives showed the highest similarity to saffron stigmas. It is important to highlight that additional crocin species were found in yellow tepals which despite differing from the ones in saffron, support the metabolic capability of these tissues and species to synthesize crocins. The identified crocins differ in the decorations of the crocetin backbone, likely displaying high glycosylation levels. Although the crocins found in the tepals are different from those in saffron stigmas, they could still be a potential source of pigments that meet industrial demands. The identification of additional CCD transcripts in the yellow crocuses may also be involved in the synthesis of non-crocin pigments, which could also be of commercial interest. Moreover, our results support the presence of diverse flavonoid in tepals, which could be exploited as a source of natural products for nutritional and medicinal purposes. Conclusions Previous observations noted that saffron tepals contain metabolites similar to the ones in stigma, leading to recent exploration of tepals as a potential compound source. Our study focused on the metabolomic and transcriptomic analysis of saffron stigma-like apocarotenoids and kaempferol-derivatives in the tepals of four common yellow-tepal crocuses. Among these crocuses, the metabolites in C. korolkowii tepals displayed the highest similarity to those in saffron stigmas, containig three verified stigma-like crocins and two kaempferol-derivatives. However, most crocins in yellow tepals differed structurally from those in saffron, featuring variations in crocetin backbone decorations and increased glycosylation. Despite these distinctions, the identified crocins in tepals remain a potential source for industrial pigments. Moreover, our results indicated a diverse array of flavonoids in tepals, presenting an additional opportunity for exploration. These flavonoids hold promise as a source for the development of natural products with applications in nutrition and medicine. In summary, while crocins in yellow-tepal crocus tepals differ from those in saffron stigmas, they still hold potential for industrial pigments. Additionally, the presence of diverse flavonoids in tepals expands the potential applications of these crocuses, offering opportunities for the development of natural products with nutritional and medicinal benefits. Methods Plant materials Stigmas of the autumn-flowering saffron crocus ( C. sativus ) along with tepal samples of four species each of the yellow, spring-flowering species C. flavus subsp. dissectus , C. graveolens and C. korolkowii were grown and collected in IPK Gatersleben greenhouses while C. chrysanthus was provided by Vladimir Randjelovic from the University of Niš. Samples for metabolite and transcriptome analyses were collected with two different methods. For metabolic analysis, tepals and stigmas were fast-dried in silica gel and stored dry till their analysis. For the transcriptome analyses, tissue samples of tepals and stigmas were stabilized in RNA later (Qiagen, Germany) to preserve RNA integrity during storage. Metabolite extraction and analysis Soluble semi-polar metabolites were extracted from tepal and stigma samples. All extraction steps were performed at room temperature under darkness. Firstly, samples were transferred into 2-ml Eppendorf tubes and weighted (5 to 15 mg). These were added with 100 µl of 50% (v/v) methanol per mg of material. Samples were grinded using a Precellys homogeneizer (Bertin Instruments, France) using 1-1.2 mm-diameter zirconium silicate beads (Mühlmeier GmbH, Germany) with four cycles of 20 s each, at a frequency of 6000 Hz. Samples were vortexed and extracted through continuous agitation at 500 rpm and room temperature for 20 min under dark conditions (Eppendorf ThermoMixer, Germany). Samples were centrifuged at 17,000 g for 5 min and the supernatant was transferred into a new tube. The pellet was re-extracted two more times, each by adding 100 µl of 50% (v/v) methanol per mg of material, followed by vortexing, shaking (20 min, 500 rpm), and centrifugation. In each step, the supernatants were combined with the first one. Extracts were covered with Argon and stored at -20 °C prior to analysis. Semi-polar metabolites were analyzed by reversed phase ultra-performance liquid chromatography coupled with photodiode array detection (RP-UPLC-PDA), using an Acquity UPLC system (Waters, Germany) equipped with an Acquity UPLC PDA eλ detector (Waters, Germany). Extracts were first centrifuged at 17,000 g for 10 min, and 100 µl of supernatant were transferred into glass vials for analysis. Sample injection volumes of 5 µl were employed for UPLC analysis, using a partial-loop with needle-overfill (PLNO) injection mode with a 10 µl loop. Compounds were separated in an Acquity UPLC BEH phenyl column (130 Å, 2.1 x 100 mm, 1.7 μm; Waters, Germany) combined with an Acquity UPLC BEH phenyl VanGuard pre-column (130 Å, 2.1 x 5 mm, 1.7 μm; Waters, Germany) using the following gradient: from 10 to 36% of solvent B for the first 3.9 min, isocratic hold from 3.9 to 5 min, from 36 to 97% B from 5 to 6 min, and an isocratic hold from 6 to 7 min to clean the column; after each run, the column was equilibrated to the starting conditions (10% B); solvent A was LC-MS grade water (CHEMSOLUTE, Th. Geyer, Germany) with 0.5% (v/v) formic acid, and solvent B was LC-MS grade acetonitrile (CHEMSOLUTE, Th. Geyer, Germany) with 0.5% (v/v) formic acid. The column temperature was maintained at 35°C and the flow was set to 500 μl min -1 . PDA-detection was performed in a range between 210 and 800 nm, at a resolution of 1.2 nm and a sampling rate of 20 points s -1 . Within the major Crocus semi-polar metabolites, crocins, picrocrocin, and flavonoids and safranal, were detected at 440, 250, and 320 nm, respectively. The identity of crocin ( trans -crocetin-digentibiosyl ester), picrocrocin, and safranal, was confirmed by retention time and UV absorption spectra with commercial standards. Processing and analysis of the acquired PDA spectra was done with the Empower 3 software (Waters, Germany). To generate a reliable reference composition, the extracts from C. sativus stigmas were further analyzed via ESI-UHR-QTOF-MS (ElectroSpray Ionization-Ultra-High-Resolution-Quadrupole Time of Flight-Mass Spectrometry) by coupling a maXis Impact ESI-QTOF MS (Bruker Daltonik GmbH; Germany) to the RP-UPLC-PDA system. The MS analyses were performed in positive and negative ionization modes. Small molecules ( i.e., safranal, picrocrocin, flavonoids; 50-1000 m/z ) were analyzed using the MS1 and MS/MS settings for barley and sunflower phenylpropanoids described in Garibay-Hernández et al. 2021. Large molecules ( i.e., crocins; 50-1500 m/z ) were analyzed with the MS1 settings in positive ionization mode for anthocyanins as described in Garibay-Hernández et al. 2021. The MS/MS analysis of large molecules in positive mode was performed in auto MS/MS using CID (Collision-Induced Dissociation) with the following settings: absolute area threshold: 5000 counts; exclusion activation: 15 spectra; exclusion release: 30 s; collision energy values (z = 1, 2, 3; isolation mass = 500; width = 8): 35, 25, 20 eV; collision energy values (z = 1, 2, 3; isolation mass = 1000; width = 10): 50, 40, 35 eV. The analysis of large molecules in negative mode was done with the following MS1 settings: 50–1000 m/z ; capillary voltage: 3.5 kV; nebulizer: 3 bar; dry gas: 8 l min − 1 ; dry temperature: 200 °C; hexapole RF (Ratio Frequency) voltage: 150 Vpp (V peak-to-peak); funnel 1 RF: 400 Vpp; funnel 2 RF: 400 Vpp; pre-pulse storage time: 8 μs; transfer time: 60 μs; low mass: 40 m/z ; collision cell RF: 800 Vpp; collision energy: 8 eV. The MS/MS analysis of large molecules in negative mode was done in auto MS/MS using CID as follows: absolute area threshold: 5000 counts; exclusion activation: 2 spectra; exclusion release: 12 s; collision energy values (z = 1; isolation mass = 500; width = 6): 20 eV; collision energy values (z = 1; isolation mass = 1000; width = 8): 20 eV. The Compass HyStar 3.2 SR2 software (Bruker Daltonik GmbH, Germany) was used to operate and coordinate LC-PDA-MS data acquisition. The analysis of MS data was performed using the Compass DataAnalysis 4.4 SR1 package (Bruker Daltonik GmbH). RNA extraction and transcriptome sequencing Total RNA was extracted using the Qiagen RNeasy plant mini kit. Total RNA purity and concentration were determined by the Qubit 1 2.0 Flurometer (Life Technologies, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, USA). A minimum amount of 50 ng μl -1 RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using the llumina TruSeq RNA Sample Prep Kit v2 for Illumina sequencing platforms. Transcriptome assembly and annotation The library preparations were sequenced on an Illumina Novaseq 6000 platform and paired-end reads were generated. In the quality control step, raw reads of fastq format were first processed through trimmomatic [28]. All the downstream analyses were based on clean data with high quality. The transcriptome was assembled using Trinity-v2.9.0 [29] with min_kmer_cov set to 2 by default and all other parameters set default. Genes linked to crocin and kaempferol biosynthesis in C. sativus were selected from [7, 17, 30–33]. These genes are listed in Table 1 and were subsequently used as a database in TBlastN and potential orthologs of these genes in the Trinity transcriptome assemblies were searched using the TBlastN parameters with e-value cut-of 0.01. Genes were then selected as potential full-length orthologs if they contained an ORF that had an overlap with at least 90% of the reference gene. Shorter sequences were listed as partial genes (Table S2). Phylogenetic trees for three main genes (CCDs, ALDHs and UGTs) are calculated using MEGA3 with ClustalW alignment along with maximum likelihood with bootstrap 500 to generate the trees. Abbreviations ALDH: Aldehyde dehydrogenase CHY-β: β-Carotene hydrolase BLAST: Basic Local Alignment Search Tool CCD: Carotenoid Cleavage Dioxygenase c-C: Cis-crocetin ESI-UHR-QTOF-MS: ElectroSpray Ionization-Ultra-High-Resolution-Quadrupole Time of Flight-Mass Spectrometry GT: Glycosyltransferase HPLC-DAD-MS: High-Performance Liquid Chromatography coupled to Diode Array Detection and Mass Spectrometry HTCC: 2,6,6-Trimethyl-4-hydroxy-1-carboxyaldehyde-1-cyclohexene K1: Kaempferol 3-O-sophoroside 7-O-glucoside K2: Kaempferol 3-O-sophoroside LC: Liquid Chromatography LC-PDA: Liquid Chromatography coupled to Photodiode Array LCYB: Lycopene β-cyclase LC-MS: Liquid Chromatography-Mass Spectrometry MEP: Methylerythritol phosphate MS: Mass Spectrometry MVA: Mevalonate pathways ORF: Open Reading Frame PDA: Photodiode Array Detection PSY: Phytoene synthase PSPG: Plant Secondary-Product-Glycosyltransferase QTOF: Quadrupole Time of Flight RNA: Ribonucleic Acid RP-UPLC-PDA: Reversed Phase Ultra-Performance Liquid Chromatography coupled with Photodiode Array Detection t-C: Trans-crocetin TBLASTN: Translated Basic Local Alignment Search Tool with Nucleotides UGTs: Uridine Diphosphate Glucosyltransferases Declarations Data availability The RNA-seq data for this study can be accessed at the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra) under accession number PRJNA926329. Acknowledgments We extend our thanks to Elena Brueckner for her expertise in metabolic profiling, Dr. Reinhard M. Fritsch for providing C . flavus subsp. dissectus samples, and Prof. Vladimir Randjelovic for sharing C . chrysanthus specimens. Funding This research was funded by Deutsche Forschungsgemeinschaft (DFG) for the financial support of part of this work. Author contributions ZN, FRB, HPM, and BU planned the experiment; ZN, SSKS, MHWS analysed of the transcriptome data; AGH performed the metabolite analysis. The first version of the manuscript was prepared by SSKS. All authors participated in reviewing and approving the final version of the paper. Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Competing interests The authors declare no competing interests. References Kazemi-Shahandashti S-S, Mann L, El-Nagish A, Harpke D, Nemati Z, Usadel B, et al. Ancient Artworks and Crocus Genetics Both Support Saffron’s Origin in Early Greece. Front Plant Sci. 2022;13:834416. Diretto G, Ahrazem O, Rubio-Moraga Á, Fiore A, Sevi F, Argandoña J, et al. UGT709G1: a novel uridine diphosphate glycosyltransferase involved in the biosynthesis of picrocrocin, the precursor of safranal in saffron ( Crocus sativus ). New Phytol. 2019;224:725–40. Zeka K, Ruparelia KC, Continenza MA, Stagos D, Vegliò F, Arroo RRJ. Petals of Crocus sativus L. as a potential source of the antioxidants crocin and kaempferol. Fitoterapia. 2015;107:128–34. Bisti S, Maccarone R, Falsini B. Saffron and retina: neuroprotection and pharmacokinetics. Vis Neurosci. 2014;31:355–61. Abdullaev FI. Cancer Chemopreventive and Tumoricidal Properties of Saffron ( Crocus sativus L). Exp Biol Med (Maywood). 2002;227:20–5. Rubio Moraga A, Ahrazem O, Rambla JL, Granell A, Gómez Gómez L. Crocins with high levels of sugar conjugation contribute to the yellow colours of early-spring flowering crocus tepals. PLoS ONE. 2013;8:e71946. Ji A, Jia J, Xu Z, Li Y, Bi W, Ren F, et al. Transcriptome-Guided Mining of Genes Involved in Crocin Biosynthesis. Front Plant Sci. 2017;8:518. Chen Y, Zhang H, Li Y-X, Cai L, Huang J, Zhao C, et al. Crocin and geniposide profiles and radical scavenging activity of gardenia fruits ( Gardenia jasminoides Ellis) from different cultivars and at the various stages of maturation. Fitoterapia. 2010;81:269–73. Ishizuka F, Shimazawa M, Umigai N, Ogishima H, Nakamura S, Tsuruma K, et al. Crocetin, a carotenoid derivative, inhibits retinal ischemic damage in mice. Eur J Pharmacol. 2013;703:1–10. Song Y-N, Wang Y, Zheng Y-H, Liu T-L, Zhang C. Crocins: A comprehensive review of structural characteristics, pharmacokinetics and therapeutic effects. Fitoterapia. 2021;153:104969. Baba SA, Mohiuddin T, Basu S, Swarnkar MK, Malik AH, Wani ZA, et al. Comprehensive transcriptome analysis of Crocus sativus for discovery and expression of genes involved in apocarotenoid biosynthesis. BMC Genomics. 2015;16:698. Nørbæk R, Brandt K, Nielsen JK, Ørgaard M, Jacobsen N. Flower pigment composition of Crocus species and cultivars used for a chemotaxonomic investigation. Biochem Syst Ecol. 2002;30:763–91. Lee G-A, Choi K-C, Hwang K-A. Kaempferol, a phytoestrogen, suppressed triclosan-induced epithelial-mesenchymal transition and metastatic-related behaviors of MCF-7 breast cancer cells. Environ Toxicol Pharmacol. 2017;49:48–57. Moraga AR, Rambla JL, Ahrazem O, Granell A, Gómez-Gómez L. Metabolite and target transcript analyses during Crocus sativus stigma development. Phytochemistry. 2009;70:1009–16. Saito K, Yonekura-Sakakibara K, Nakabayashi R, Higashi Y, Yamazaki M, Tohge T, et al. The flavonoid biosynthetic pathway in Arabidopsis: Structural and genetic diversity. Plant Physiol Biochem. 2013;72:21–34. Duan L, Ding W, Liu X, Cheng X, Cai J, Hua E, et al. Biosynthesis and engineering of kaempferol in Saccharomyces cerevisiae . Microb Cell Fact. 2017;16:165. Winterhalter P, Straubinger M. Saffron-Renewed interest in an ancient spice. Food Reviews Int. 2000;16:39–59. Ravishankar B, Shukla VJ. Indian systems of medicine: a brief profile. Afr J Tradit Complement Altern Med. 2007;4:319–37. Hosseini A, Razavi BM, Hosseinzadeh H. Saffron ( Crocus sativus ) petal as a new pharmacological target: a review. Iran J Basic Med Sci. 2018;21:1091–9. Naidis P, Lykidou S, Mischopoulou M, Vouvoudi E, Nikolaidis NF. Study of the dyeing properties of saffron and ultrafiltrated saffron powders, as colourants for natural and synthetic fibres. Color Technol. 2023;139:565–77. Harpke D, Carta A, Tomović G, Ranđelović V, Ranđelović N, Blattner FR, et al. Phylogeny, karyotype evolution and taxonomy of Crocus series Verni (Iridaceae). Plant Syst Evol. 2015;301:309–25. Harpke D, Kerndorff H, Pasche E, Peruzzi L. Neotypification of the name Crocus biflorus Mill. (Iridaceae) and its consequences in the taxonomy of the genus. Phytotaxa. 2016;260:131. Peruzzi L. Crocus heuffelianus (Iridaceae), a new record for the Italian flora. Phytotaxa. 2016;261:291. Raca I, Harpke D, Shuka L, Ranđelović V. A new species of Crocus ser. Verni (Iridaceae) with 2 n = 12 chromosomes from the Balkans. Plant Biosystems - Int J Dealing all Aspects Plant Biology. 2022;156:36–42. Shuka D, Kerndorff H, Waminal N, Ranđelovič V, Harpke D. Two new species of Crocus series Scardici from the Western Balkans. preprint. Evolutionary Biology; 2022. Mohtashami L, Amiri MS, Ramezani M, Emami SA, Simal-Gandara J. The genus Crocus L.: A review of ethnobotanical uses, phytochemistry and pharmacology. Ind Crops Prod. 2021;171:113923. Zhu C, Bai C, Sanahuja G, Yuan D, Farré G, Naqvi S, et al. The regulation of carotenoid pigmentation in flowers. Arch Biochem Biophys. 2010;504:132–41. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52. Baba SA, Ashraf N. Apocarotenoids of Crocus sativus L: From biosynthesis to pharmacology. Singapore: Springer Singapore; 2016. Ahrazem O, Diretto G, Argandoña J, Rubio-Moraga Á, Julve JM, Orzáez D, et al. Evolutionarily distinct carotenoid cleavage dioxygenases are responsible for crocetin production in Buddleja davidii . J Exp Bot. 2017;68:4663–77. Demurtas OC, Frusciante S, Ferrante P, Diretto G, Azad NH, Pietrella M, et al. Candidate Enzymes for Saffron Crocin Biosynthesis Are Localized in Multiple Cellular Compartments. Plant Physiol. 2018;177:990–1006. López-Jimenez AJ, Frusciante S, Niza E, Ahrazem O, Rubio-Moraga Á, Diretto G, et al. A New Glycosyltransferase Enzyme from Family 91, UGT91P3, Is Responsible for the Final Glucosylation Step of Crocins in Saffron ( Crocus sativus L). Int J Mol Sci. 2021;22:8815. Hughes J, Hughes MA. Multiple secondary plant product UDP-glucose glucosyltransferase genes expressed in cassava (Manihot esculenta Crantz) cotyledons. DNA Seq. 1994;5:41–9. Nemati Z, Blattner FR, Kerndorff H, Erol O, Harpke D. Phylogeny of the saffron-crocus species group, Crocus series Crocus (Iridaceae). Mol Phylogenet Evol. 2018;127:891–7. Nemati Z, Harpke D, Gemicioglu A, Kerndorff H, Blattner FR. Saffron (Crocus sativus) is an autotriploid that evolved in Attica (Greece) from wild Crocus cartwrightianus . Mol Phylogenet Evol. 2019;136:14–20. Suchareau M, Bordes A, Lemée L. Improved quantification method of crocins in saffron extract using HPLC-DAD after qualification by HPLC-DAD-MS. Food Chem. 2021;362:130199. Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54:733–49. Timmins JJB, Kroukamp H, Walker RSK, Pretorius IS, Paulsen IT. Comparative Evaluation of Secreted Plant Carotenoid Cleavage Dioxygenase 1 (CCD1) Enzymes in Saccharomyces cerevisiae . Fermentation. 2022;8:395. Ohmiya A, Kato M, Shimada T, Nashima K, Kishimoto S, Nagata M. Molecular Basis of Carotenoid Accumulation in Horticultural Crops. Hort J. 2019;88:135–49. Klingner A, Hundeshagen B, Kernebeck H, Bothe H. Localization of the yellow pigment formed in roots of gramineous plants colonized by arbuscular fungi. Protoplasma. 1995;185:50–7. Fester T, Schmidt D, Lohse S, Walter MH, Giuliano G, Bramley PM, et al. Stimulation of carotenoid metabolism in arbuscular mycorrhizal roots. Planta. 2002;216:148–54. Ohmiya A, Kishimoto S, Aida R, Yoshioka S, Sumitomo K. Carotenoid cleavage dioxygenase (CmCCD4a) contributes to white color formation in chrysanthemum petals. Plant Physiol. 2006;142:1193–201. Ma G, Zhang L, Matsuta A, Matsutani K, Yamawaki K, Yahata M, et al. Enzymatic formation of β-citraurin from β-cryptoxanthin and Zeaxanthin by carotenoid cleavage dioxygenase4 in the flavedo of citrus fruit. Plant Physiol. 2013;163:682–95. Rubio A, Rambla JL, Santaella M, Gómez MD, Orzaez D, Granell A, et al. Cytosolic and plastoglobule-targeted carotenoid dioxygenases from Crocus sativus are both involved in beta-ionone release. J Biol Chem. 2008;283:24816–25. Le Roy J, Huss B, Creach A, Hawkins S, Neutelings G. Glycosylation Is a Major Regulator of Phenylpropanoid Availability and Biological Activity in Plants. Front Plant Sci. 2016;7:735. Pfister S, Meyer P, Steck A, Pfander H. Isolation and Structure Elucidation of Carotenoid – Glycosyl Esters in Gardenia Fruits ( Gardenia jasminoides Ellis) and Saffron ( Crocus sativus Linne). J Agric Food Chem. 1996;44:2612–5. Dufresne C, Cormier F, Dorion S, Niggli UA, Pfister S, Pfander H. Glycosylation of encapsulated crocetin by a Crocus sativus L. cell culture. Enzym Microb Technol. 1999;24:453–62. Tarantilis PA, Tsoupras G, Polissiou M. Determination of saffron (Crocus sativus L.) components in crude plant extract using high-performance liquid chromatography-UV-visible photodiode-array detection-mass spectrometry. J Chromatogr A. 1995;699:107–18. Castillo R, Fernández J-A, Gómez-Gómez L. Implications of carotenoid biosynthetic genes in apocarotenoid formation during the stigma development of Crocus sativus and its closer relatives. Plant Physiol. 2005;139:674–89. Kanakis CD, Daferera DJ, Tarantilis PA, Polissiou MG. Qualitative determination of volatile compounds and quantitative evaluation of safranal and 4-hydroxy-2,6,6-trimethyl-1-cyclohexene-1-carboxaldehyde (HTCC) in Greek saffron. J Agric Food Chem. 2004;52:4515–21. Moras B, Loffredo L, Rey S. Quality assessment of saffron ( Crocus sativus L.) extracts via UHPLC-DAD-MS analysis and detection of adulteration using gardenia fruit extract ( Gardenia jasminoides Ellis). Food Chem. 2018;257:325–32. Li C-Y, Wu T-S. Constituents of the pollen of Crocus sativus L. and their tyrosinase inhibitory activity. Chem Pharm Bull (Tokyo). 2002;50:1305–9. Mykhailenko O, Kovalyov V, Goryacha O, Ivanauskas L, Georgiyants V. Biologically active compounds and pharmacological activities of species of the genus Crocus : A review. Phytochemistry. 2019;162:56–89. Zengin G, Aumeeruddy MZ, Diuzheva A, Jekő J, Cziáky Z, Yıldıztugay A, et al. A comprehensive appraisal on Crocus chrysanthus (Herb.) Herb. flower extracts with HPLC-MS/MS profiles, antioxidant and enzyme inhibitory properties. J Pharm Biomed Anal. 2019;164:581–9. Harborne JB, Williams CA. 6-Hydroxyflavones and Other Flavonoids of Crocus . Z für Naturforschung C. 1984;39:18–23. Additional Declarations No competing interests reported. Supplementary Files YellowtepalcrocusesSupplementaryFiguresTables.docx Supplementary Information Fig. S1. Representative crocins present in tissues from Crocus species. Representative UV-Vis absorption spectra for major crocins detected in the analyzed samples. Fig. S2. Major crocins in C. sativus stigmas and yellow tepals from multiple Crocus species. Fig. S3. Analysis of flavonoids and safranal in tissues from multiple Crocus species. Fig. S4. Analysis of HTTC and picrocrocin in tissues from multiple Crocus species. Table S1. Table of metabolites quantities detected in C. sativus and yellow-tepal crocuses. Table S2. List of all potential isoform genes resulted by TBlastN. Cite Share Download PDF Status: Published Journal Publication published 06 May, 2024 Read the published version in BMC Plant Biology → Version 1 posted Editorial decision: Revision requested 19 Mar, 2024 Reviews received at journal 08 Mar, 2024 Reviewers agreed at journal 28 Feb, 2024 Reviewers invited by journal 28 Feb, 2024 Editor invited by journal 14 Feb, 2024 Editor assigned by journal 14 Feb, 2024 Submission checks completed at journal 14 Feb, 2024 First submitted to journal 11 Feb, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-3948859","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":272896573,"identity":"857b2e0c-799c-47d3-b688-61a19926d159","order_by":0,"name":"Zahra Nemati","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABHElEQVRIiWNgGAWjYNACAxseBh4Yh52BDUQxNuDTccAgDaolAYiZidLCcJiBeC3m7N2Jnz8UnJcx5znA/OLnD5t8+WbmYw8+7mCQ7cehxbLn7GaJAwa3eSx7G9gsexLSLDccZks3nHmGwXgmDmsMbuRuAGsxOM//zYAn4bCBATOPmTRvG0PihgM4tNx/u/nHAYNzQC0MbIZ/Ev4byDfzf5P+C9SyH5eWG7zbgLYc4DE428D8mCfhgAHDYR42aUaQLTj8YnAmd5vFGYNkHsueA2zMMmnJBgaH2cwNe9skjGfgsuX42c03Kv7Y2ZvzJDB/fGNjZyDf3vzswc82G9l+HN5H6GVgYJNA4kvgVImshfkDYWWjYBSMglEwEgEAA1RcnOCBLp4AAAAASUVORK5CYII=","orcid":"","institution":"Leibniz Institute of Plant Genetics and Crop Plant Research (IPK)","correspondingAuthor":true,"prefix":"","firstName":"Zahra","middleName":"","lastName":"Nemati","suffix":""},{"id":272896574,"identity":"57194125-dba3-4e58-9133-45f16ba92b7a","order_by":1,"name":"Seyyedeh‌ Sanam Kazemi Shahandashti","email":"","orcid":"","institution":"CEPLAS, Forschungszentrum Jülich GmbH","correspondingAuthor":false,"prefix":"","firstName":"Seyyedeh‌","middleName":"Sanam Kazemi","lastName":"Shahandashti","suffix":""},{"id":272896575,"identity":"c95a105d-a9fa-4b7a-b931-3f8415335c44","order_by":2,"name":"Adriana Garibay-Hernández","email":"","orcid":"","institution":"Technische Universität Kaiserslautern","correspondingAuthor":false,"prefix":"","firstName":"Adriana","middleName":"","lastName":"Garibay-Hernández","suffix":""},{"id":272896576,"identity":"d417c3bd-cb27-443b-b16f-db6bb9b86b55","order_by":3,"name":"Hans-Peter Mock","email":"","orcid":"","institution":"Leibniz Institute of Plant Genetics and Crop Plant Research (IPK)","correspondingAuthor":false,"prefix":"","firstName":"Hans-Peter","middleName":"","lastName":"Mock","suffix":""},{"id":272896577,"identity":"d76c04af-7da4-43c8-ab49-2cc30eade319","order_by":4,"name":"Maximilian H-W Schmidt","email":"","orcid":"","institution":"Geisenheim University","correspondingAuthor":false,"prefix":"","firstName":"Maximilian","middleName":"H-W","lastName":"Schmidt","suffix":""},{"id":272896578,"identity":"2bebce5d-cc4e-4d0e-8acc-e5ab407567c2","order_by":5,"name":"Björn Usadel","email":"","orcid":"","institution":"CEPLAS, Forschungszentrum Jülich GmbH","correspondingAuthor":false,"prefix":"","firstName":"Björn","middleName":"","lastName":"Usadel","suffix":""},{"id":272896579,"identity":"1367ae84-98e5-4dba-833e-fbbcf523735a","order_by":6,"name":"Frank.R Blattner","email":"","orcid":"","institution":"Leibniz Institute of Plant Genetics and Crop Plant Research (IPK)","correspondingAuthor":false,"prefix":"","firstName":"Frank.R","middleName":"","lastName":"Blattner","suffix":""}],"badges":[],"createdAt":"2024-02-11 16:14:43","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-3948859/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-3948859/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12870-024-05036-1","type":"published","date":"2024-05-07T03:57:57+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":51209179,"identity":"94fd345c-fc64-4073-bb88-f53038ca23d6","added_by":"auto","created_at":"2024-02-16 04:35:26","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":581558,"visible":true,"origin":"","legend":"\u003cp\u003eMajor soluble semi-polar metabolites in \u003cem\u003eC. sativus\u003c/em\u003e stigmas. (A) Representative LC-PDA isoplot chromatogram of semi-polar compounds detected within the range of 230-550 nm in extracts from \u003cem\u003eC. sativus\u003c/em\u003e stigmas. Three major compound classes displaying maximum UV-Vis absorption values at specific wavelengths are indicated. (B) LC-PDA chromatograms extracted at 440, 320, and 250 nm, showing three major groups of semi-polar metabolites in \u003cem\u003eC. sativus\u003c/em\u003e stigmas. (C) Representative UV-Vis absorption spectra for selected stigma metabolites.\u003c/p\u003e","description":"","filename":"floatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/97a00acba7f39c6e3c00988e.png"},{"id":51209182,"identity":"102a1c30-fc35-4070-9707-0aadd549af2a","added_by":"auto","created_at":"2024-02-16 04:35:27","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":698970,"visible":true,"origin":"","legend":"\u003cp\u003eMajor soluble semi-polar metabolites in yellow-tepal Crocus species compared to \u003cem\u003eC. sativus\u003c/em\u003e stigmas. LC-PDA isoplot chromatograms of semi-polar compounds detected within the range of 230-550 nm in extracts from \u003cem\u003eC. sativus\u003c/em\u003e stigmas and yellow tepals from different Crocus species. The UV-Vis absorption regions where crocins and flavonoids display their maximum absorbance are indicated with red and blue arrows, respectively. Specific crocins are labeled (red color); their UV-Vis absorption spectra are shown in Supplementary Fig. S1.\u003c/p\u003e","description":"","filename":"floatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/abe2b0a4d3a3c4c20c183c68.png"},{"id":51209185,"identity":"39978610-4303-4ec0-ae2b-79ba2bd57aac","added_by":"auto","created_at":"2024-02-16 04:35:30","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":215949,"visible":true,"origin":"","legend":"\u003cp\u003eHeatmap of saffron-like metabolites. Metabolites (relative abundance) are represented on the right-hand side of the plot and the Crocus species are indicated in the bottom of the heatmap.\u003c/p\u003e","description":"","filename":"floatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/7323150f82ebbfff72e51c8d.png"},{"id":51209180,"identity":"91137e39-7be0-4310-9aff-37cbc9d83d1a","added_by":"auto","created_at":"2024-02-16 04:35:27","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":289154,"visible":true,"origin":"","legend":"\u003cp\u003eGeneral legend for Figure. A. Phylogenetic trees of \u003cem\u003eC. sativus\u003c/em\u003e gene families involved in the biosynthesis of value-added compounds. B. Major value-added metabolites of \u003cem\u003eC. sativus\u003c/em\u003e stigma are shown in the circle along with their pathways and the related genes and metabolites. In this study, the metabolites of the tepals of \u003cem\u003eC. korolkowii\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, \u003cem\u003eC. flavus\u003c/em\u003e ssp. \u003cem\u003edissectus\u003c/em\u003eand \u003cem\u003eC. chrysanthus\u003c/em\u003e were investigated to check if they are similar to the ones in saffron stigma. The yellow flowers in the figure indicate the metabolites or genes that were detected in the yellow-tepals crocuses.\u003c/p\u003e","description":"","filename":"floatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/3dd73d9127eac4b114db564e.png"},{"id":56140235,"identity":"f9cb160f-fed9-495a-91ba-91461f0b1606","added_by":"auto","created_at":"2024-05-09 04:09:54","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1631167,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/e28df999-59a9-4ef3-a3d3-3bc1816f752d.pdf"},{"id":51209184,"identity":"eff1073d-ec21-44cf-9ed7-6a1874c987ba","added_by":"auto","created_at":"2024-02-16 04:35:28","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":735746,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupplementary Information\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFig. S1\u003c/strong\u003e.\u003cstrong\u003e \u003c/strong\u003eRepresentative crocins present in tissues from \u003cem\u003eCrocus \u003c/em\u003especies. Representative UV-Vis absorption spectra for major crocins detected in the analyzed samples. \u003cstrong\u003eFig. S2\u003c/strong\u003e. Major crocins in \u003cem\u003eC. sativus\u003c/em\u003e stigmas and yellow tepals from multiple \u003cem\u003eCrocus\u003c/em\u003e species. \u003cstrong\u003eFig. S3\u003c/strong\u003e. Analysis of flavonoids and safranal in tissues from multiple \u003cem\u003eCrocus\u003c/em\u003e species. \u003cstrong\u003eFig. S4\u003c/strong\u003e. Analysis of HTTC and picrocrocin in tissues from multiple \u003cem\u003eCrocus \u003c/em\u003especies. \u003cstrong\u003eTable S1\u003c/strong\u003e. Table of metabolites quantities detected in \u003cem\u003eC. sativus\u003c/em\u003e and yellow-tepal crocuses. \u003cstrong\u003eTable S2\u003c/strong\u003e. List of all potential isoform genes resulted by TBlastN.\u003c/p\u003e","description":"","filename":"YellowtepalcrocusesSupplementaryFiguresTables.docx","url":"https://assets-eu.researchsquare.com/files/rs-3948859/v1/3e2d9aba5fb39a816283897a.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites","fulltext":[{"header":"Highlights","content":"\u003cp\u003eThis study uncovers yellow-tepal crocuses as sources of saffron-like compounds, meeting industry demand. Metabolomic and transcriptomic research reveals insights into crocins and flavonoids biosynthesis, offering a promising production avenue.\u003c/p\u003e"},{"header":"Background","content":"\u003cp\u003eThe spice saffron consists of vivid crimson threads which are indeed dried stigmas of \u003cem\u003eCrocus sativus\u003c/em\u003e. To obtain one kilogram of saffron 150,000-160,000 flowers are required. Literally, the small amount of saffron obtained per plant, along with its manual harvesting, leads to saffron being the most expensive spice in the world. From culinary uses and art crafts to medicine and health, the value of saffron has risen over the course of time. While saffron fulfills different ranges of commercial demands these days, its application is not limited to a spice or colorant [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. It has a longstanding reputation in traditional medicine and researchers found biomedical and pharmacological properties of its metabolites. Saffron chemical composition is characterized as being rich in apocarotenoids and comprising kaempferol-derivatives as its sole flavonoid species. The apocarotenoids mainly consist of crocetin and its glucosyl esters, also named crocins (which occur as \u003cem\u003etrans-\u003c/em\u003e and \u003cem\u003ecis-\u003c/em\u003e isomers), 2,6,6-trimethyl-4-hydroxy-1-carboxyaldehyde-1-cyclohexene (HTCC), picrocrocin (which is glucosylated HTCC), and safranal, which is a product of picrocrocin catabolism [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Crocins and safranal are mainly responsible for saffron\u0026rsquo;s color and its unique aroma, respectively [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. These metabolites also exert pharmacological effects on the prevention and/or improvement of various diseases related to the gastrointestinal, cardiovascular, endocrine, nervous, and immune systems [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eBelonging to plant apocarotenoids [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e], crocins are found in crocin-producing plants such as crocuses, \u003cem\u003eBuddleja davidii\u003c/em\u003e, and \u003cem\u003eGardenia jasminoides\u003c/em\u003e [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. The interaction among general isoprenoid pathways (MEP: methylerythritol phosphate and MVA: mevalonate pathways), carotenoid, and apocarotenoid biosynthetic pathways [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e], results in the presence of a wide variety of crocins in crocin-producing plants. In case of saffron stigmas, these include crocin-I (\u003cem\u003et\u003c/em\u003e-C-4: \u003cem\u003etrans\u003c/em\u003e-crocetin di(β-D-gentibiosyl) ester), with the highest abundance more than 70% of all crocins, crocin-II (\u003cem\u003et\u003c/em\u003e-C-3: \u003cem\u003etrans\u003c/em\u003e-crocetin (β-D-gentibiosyl) (β-D-glucosyl) ester), crocin-III (\u003cem\u003et-\u003c/em\u003eC-2: \u003cem\u003etrans\u003c/em\u003e-crocetin (β-D-gentibiosyl) ester), crocin-IV (\u003cem\u003et-\u003c/em\u003eC-2: \u003cem\u003etrans\u003c/em\u003e-crocetin di(β-D-glucosyl) ester), and crocin-V (\u003cem\u003et-\u003c/em\u003eC-1: \u003cem\u003etrans\u003c/em\u003e-crocetin (β-D-glucosyl) ester) [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e, \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. These five forms differ regarding the position and number of glucosyl moieties attached to the crocetin aglycone backbone [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. In the general isoprenoid pathways, the precursors of carotenoids are generated via MEP and MVA, followed by the activity of enzymes from the upstream carotenoid biosynthetic pathway such as phytoene synthase (PSY), lycopene β-cyclase (LCYB), and β-carotene hydrolase (BHY). The next stages of saffron apocarotenoid synthesis are accomplished by three key enzymes including the carotenoid cleavage dioxygenase (CCD2), aldehyde dehydrogenase (ALDH), as well as crocetin and HTCC glucosyltransferases (UGTs) [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eBeing also commercially valuable compounds, kaempferol and its glycosylated derivatives were also identified in crocuses. In fact, kaempferol and its derivatives comprise 90% of the total flavonoid content in tepals of \u003cem\u003eC. sativus\u003c/em\u003e [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e, \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Quercetin glycosides are the second most abundant flavonoid class in crocus tepals, comprising 5\u0026ndash;10% of the flavonoid compounds, whereas other components like dihydrokaempferol glycosides and naringenin are observed in negligible amounts [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. At the initial stage of the flavonoid biosynthetic pathway, phenylalanine is converted into \u003cem\u003ep-\u003c/em\u003ecoumaryl-CoA by the sequential activity of phenylalanine ammonia lyase, cinnamic acid 4-hydrolase, and \u003cem\u003ep\u003c/em\u003e-coumaric acid:CoA ligase. The chalcone synthase catalyzes the first committed step into flavonoid biosynthesis through the condensation of \u003cem\u003ep-\u003c/em\u003ecoumaroyl-CoA with malonyl-CoA, resulting in the formation of naringenin chalcone, which the chalcone isomerase transforms into naringenin that is further metabolized into multiple flavonoids in a plant lineage-specific manner. Flavonoids are extensively decorated by diverse chemical reactions, including the glycosylation of flavonoid aglycones that is catalyzed by glycosyltransferases [\u003cspan additionalcitationids=\"CR15\" citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eRegarding the growing demand for saffron, tepals of \u003cem\u003eC. sativus\u003c/em\u003e, which are the major by-products of saffron production, could be an additional source of saffron-like specialized metabolites [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. The presence of crocin as well as kaempferol derivatives in tepals of C. sativus was confirmed through multiple analytical approaches such as liquid chromatography coupled to diode array detection (HPLC-DAD) and infrared (IR) spectroscopy [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. These studies identified at least six different crocins in \u003cem\u003eC. sativus\u003c/em\u003e tepals and the remarkable amount of kaempferol derivatives of 126 mg g\u003csup\u003e-1\u003c/sup\u003e dry weight [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Owing to these metabolites\u0026rsquo; presence in saffron tepals, they have been suggested as an alternative or supplementary medicine against some diseases [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e] and effectively applied as natural pigments in natural and synthetic substrates [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e] \u003cem\u003eC. sativus\u003c/em\u003e tepals are not the only extra source of these metabolites. There are about 250 accepted species within the genus \u003cem\u003eCrocus\u003c/em\u003e (Iridaceae) around the world [\u003cspan additionalcitationids=\"CR22 CR23 CR24\" citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. They contain a broad spectrum of natural products including flavonoids, coumarins, additional phenolic glycosides, alkaloids, monoterpenoids, sesquiterpenoids, carotenoids, α-hydroxyacids, fatty alcohols, fatty aldehydes, fatty esters, fatty acids, and alkanes [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Within \u003cem\u003eCrocus\u003c/em\u003e, compounds similar to those from saffron stigmas may occur in the tepals of different species. Being divided into spring- and autumn-flowering crocuses, the color of their flower\u0026rsquo;s ranges from white over yellow in spring, to purple and blue-lilac in autumn species [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. In general, the presence of carotenoids and flavonoids in flowers leads to petal pigmentation in the range of yellow to orange [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. While saffron purple tepals contain crocin and flavonoids, it is possible that the tepals of yellow crocuses also contain saffron-like metabolites due to their yellow color and thus, may be a potential source of these metabolites. Since saffron stigmas and petals are the sources of saffron main metabolites in autumn, yellow-tepal crocuses could increase the yield per year if they contain saffron-like metabolites.\u003c/p\u003e \u003cp\u003eIn this study, we compared the semi-polar metabolite profiles of the four yellow-tepal \u003cem\u003eCrocus\u003c/em\u003e taxa \u003cem\u003eC. chrysanthus\u003c/em\u003e, \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, and \u003cem\u003eC. korolkowii\u003c/em\u003e via LC-PDA analysis to investigate the presence of compounds similar to those in saffron stigmas. We also provide insights into the biosynthesis of crocins and kaempferol derivatives in yellow-tepal crocuses through their transcriptomic analysis with special focus on the active genes involved in carotenoid and flavonoid metabolism.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eApocarotenoid and flavonoid profiles of yellow tepals compared to saffron stigmas\u003c/h2\u003e \u003cp\u003eTo investigate the presence of apocarotenoids and flavonoids similar to those of \u003cem\u003eC. sativus\u003c/em\u003e stigmas in yellow-tepal crocus plants, the composition of \u003cem\u003eC. sativus\u003c/em\u003e stigmas was initially assessed through LC-PDA-MS (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eWe annotated a total of ten major compounds in \u003cem\u003eC. sativus\u003c/em\u003e stigmas comprising both, apocarotenoids and flavonoids: (i) safranal, (ii) HTCC (picrocrocin precursor), (iii) picrocrocin, (iv) four \u003cem\u003etrans\u003c/em\u003e-crocins including the well-known saffron crocin \u003cem\u003etrans\u003c/em\u003e-crocetin-digentibiosyl ester (\u003cem\u003et\u003c/em\u003e-C-4), (v) one \u003cem\u003ecis\u003c/em\u003e-crocin (\u003cem\u003ec-\u003c/em\u003eC-4), and (vi) two major flavonoids, kaempferol 3\u003cem\u003e-O-\u003c/em\u003esophoroside 7\u003cem\u003e-O-\u003c/em\u003eglucoside (K1) and kaempferol 3\u003cem\u003e-O-\u003c/em\u003esophoroside (K2). The most abundant crocins in the stigma samples were \u003cem\u003et-\u003c/em\u003eC-3, \u003cem\u003et-\u003c/em\u003eC-4, and \u003cem\u003ec-\u003c/em\u003eC-4, while the crocins \u003cem\u003et-\u003c/em\u003eC-2 and \u003cem\u003et-\u003c/em\u003eC-5 were identified in lower abundance (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAfter having established this baseline, we compared via LC-PDA the composition of the stigmas from \u003cem\u003eC. sativus\u003c/em\u003e to that of yellow-tepals of \u003cem\u003eC. chrysanthus\u003c/em\u003e, \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, and \u003cem\u003eC. korolkowii\u003c/em\u003e to determine whether similar compounds were found (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e2\u003c/span\u003e,\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e3\u003c/span\u003e). A heatmap shows the metabolite composition of \u003cem\u003eC. sativus\u003c/em\u003e stigmas compared to that from tepals of yellow-tepal crocuses (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). We did not find the monoterpenoids HTCC, picrocrocin, and safranal in any of the yellow-tepal samples. In contrast, the yellow-tepals from \u003cem\u003eC. chrysanthus\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e shared the capability of synthesizing \u003cem\u003etrans-\u003c/em\u003ecrocins that are also found in saffron stigmas. The tepals from these two species displayed the \u003cem\u003et\u003c/em\u003e-C-5 crocin, whereas only \u003cem\u003eC. korolkowii\u003c/em\u003e also showed \u003cem\u003et\u003c/em\u003e-C-4 and \u003cem\u003et\u003c/em\u003e-C-2, which are the major crocins found in saffron stigmas (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e1\u003c/span\u003e, [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]). Interestingly, the yellow-tepals from all the \u003cem\u003eCrocus\u003c/em\u003e species displayed additional crocins in their profiles, as confirmed by their UV-Vis spectra (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e2\u003c/span\u003e, \u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, S2). These major crocins detected in the tepal samples did not correspond to the dominant crocins in \u003cem\u003eC. sativus\u003c/em\u003e stigmas (Fig. S2).\u003c/p\u003e \u003cp\u003eThe tepals displayed rich flavonoid profiles that clearly differ from those of saffron stigmas (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e2\u003c/span\u003e, S3). Only \u003cem\u003eC. korolkowii\u003c/em\u003e displayed the two kaempferol derivatives that we also found in the saffron stigmas, kaempferol 3-O-sophoroside and kaempferol 3-O-sophoroside 7-O-glucoside (K2 and K1, respectively; Fig. S3).\u003c/p\u003e \u003cp\u003eTo gain insights into the biosynthetic genes related to saffron-like metabolites in the four yellow-tepal \u003cem\u003eCrocus\u003c/em\u003e species, the tepals of \u003cem\u003eC. chrysanthus\u003c/em\u003e, \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, \u003cem\u003eC. korolkowii\u003c/em\u003e were used for transcriptomic analysis. A total of 177,738,884 clean reads were obtained across all the samples. We compared the assembled sequences to known \u003cem\u003eC. sativus\u003c/em\u003e genes involved in apocarotenoid and flavonoid metabolism (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The candidate genes were categorized into three major groups involved saffron apocarotenoid and flavonoid metabolism: carotenoid cleavage dioxygenases (CCDs), uridine diphosphate glucosyltransferases (UGTs), and aldehyde dehydrogenases (ALDHs). A CCD, CCD2, catalyzes the first step of apocarotenoid biosynthesis in saffron, being responsible for the oxidative cleavage of zeaxanthin and the production of HTTC and crocetin dialdehyde [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. ALDHs are responsible for the transformation of crocetin dialdehyde into crocetin. The UGTs are responsible for glycosylation of low molecular weight substrates such as crocetin, HTTC, and flavonoid aglycones to produce crocin, picrocrocin, and kaempferol derivatives, respectively. The phylogenetic trees of these gene families are provided, indicating their roles in saffron metabolites pathways (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). Genes associated with the CCDs family were selected from both \u003cem\u003eC. sativus\u003c/em\u003e and other close-related crocus species to encompass a wider range of genetic diversity, as CCD2 is the one characterized/known in \u003cem\u003eC. sativus\u003c/em\u003e. Conversely, the selection of genes linked to the UGTs and ALDHs families was limited to those exclusively found in \u003cem\u003eC. sativus\u003c/em\u003e, given the substantial abundance of genes from these families within this particular species.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eGenes related to apocarotenoid and kaempferol biosynthesis pathways of C. sativus in yellow-tepal crocuses\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eName\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eGenBank accession no.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eC. korolkowii\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eC. chrysantus\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eC. graveolens\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eC. \u003cem\u003eflavus\u003c/em\u003e ssp. dissectus\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eE-value/%identitiy\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eE-value/identitiy%\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eE-value/identitiy%\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eE-value/identitiy%\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e carotenoid cleavage dioxygenase 2 (CCD2)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKJ541749\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1e-161 / 70\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 73\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3.00e-116 / 77\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1e-46 / 74\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e carotenoid cleavage dioxygenase 2L (CCD2L)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKP887110\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1e-161 / 66\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 73\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-115 / 77\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2e-46 / 74\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e chromoplast carotenoid cleavage dioxygenase 4b (CCD4b)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEU523663.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2e-142/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e chromoplast carotenoid cleavage dioxygenase 4a (CCD4a)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEU523662.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2e-142/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e carotenoid cleavage dioxygenase 1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eOL606625.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0/97%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e3E-139/97%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. ancyrensis\u003c/em\u003e carotenoid cleavage dioxygenase (CCD4c)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKP792758.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0/89%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e3E-115/90%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. ancyrensis\u003c/em\u003e carotenoid cleavage dioxygenase (CCD4a/b)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKP792757.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e4E-133/98%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0/99%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/90%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e8E-087/89%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. ancyrensis\u003c/em\u003e carotenoid cleavage dioxygenase (CCD1)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKP792755.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0/98%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/98%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e9E-146/99%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e carotenoid cleavage dioxygenase 1 (CCD1)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMN540633.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0/97%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/96%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e3E-139/97%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e chromoplast carotenoid cleavage dioxygenase 4b (CCD4b)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEU523663.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2e-142/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0/87%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNo hits found\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e UGT709G1 mRNA, complete cds (uridine diphosphate glycosyltransferase)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKX385186.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e3e-62 / 67\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 61\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3e-101 / 38\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e6e-06 / 40\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase 2 (GLT2) mRNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAY262037.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e6e-65 / 74\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1e-159 / 61\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-49 / 30\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e7e-46 / 66\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e crocetin glucosyltransferase 74AD1 mRNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMF596166.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e3e-74 / 36\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 81\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-56 / 31\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2e-51 / 71\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e beta-carotene hydroxylase to a 915-base sequence of most likely codons\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCAC95130\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1e-144 / 87\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3e-144 / 86\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-13 / 90\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.62 / 50\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e aldehyde dehydrogenase 2B4 (ALDH2B4) mRNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMG672523.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 98\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 44\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-125 / 96\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2.7 / 47\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e aldehyde dehydrogenase ALDH2C4 mRNA, complete\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMF596160.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 81\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 97\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-104 / 82\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.024 / 39\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e aldehyde dehydrogenase 3I1 mRNA, complete\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMF596165.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 96\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e4e-102 / 58\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-50 / 94\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e5.8 / 47\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e aldehyde dehydrogenase ALDH5F1 mRNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMF596161.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 96\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e4e-102 / 38\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-50 / 94\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e5.8 / 47\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e aldehyde dehydrogenase 6B2 (ALDH6B2) mRNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMG672524.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 97\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 97\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3.00e-93 / 97\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.4 / 35\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e Aldehyde dehydrogenase ALDH7B4 mRNA, complete\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMF596162.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 99\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 99\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e4.00e-171 / 99\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2.3 / 41\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e Glucosyltransferase (UGT91P3)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190170.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e4e-90 / 62\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e7e-121 / 86\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-43 / 92\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.68 / 30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase (UGT91P6)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190175.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e4e-117 / 79\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1e-100 / 45\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e6.00e-38 / 31\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.010 / 42\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase (UGT91K3)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190174.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1e-69 / 72\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 43\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e7.00e-40 / 30\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1e-07 / 45\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase (UGT91K2)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190173.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e5e-67 / 73\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 78\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-36 / 31\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2e-08 / 47\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase (UGT91P5)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190172.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2e-97 / 64\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1e-121 / 82\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1.00e-43 / 92\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e5.0 / 33\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eglucosyltransferase (UGT91P4)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMZ190171.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2e-134 / 81\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1e-106 / 46\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3.00e-39 / 31\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.011 / 45\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e glucosyltransferase (Kaempferol)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eHE793682.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 92\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 90\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e5.00e-104 / 91\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2e-05 / 29\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e flavonoid glucosyltransferase (GT45) gene\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eFJ194947.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.0 / 69\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 70\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.00e-47 / 28\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e8e-36 / 81\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eC. sativus\u003c/em\u003e UDP-glucose-dependent flavonoid UGT703B1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eKJ381079\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e7e-92 / 41\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.0 / 89\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.0 / 68\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1e-10 / 34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eWhile only a few of these genes were annotated in \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e with a low identity percentage, meaningful values for all the candidate genes were observed in \u003cem\u003eC. chrysanthus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, and \u003cem\u003eC. korolkowii\u003c/em\u003e (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Transcripts of these three species showed a considerable similarity to the candidate saffron genes. When considering the identity percentage of above 70%, the highest number of reads with the highest identity percentage were observed in \u003cem\u003eC. korolkowii.\u003c/em\u003e In fact, the transcripts of \u003cem\u003eC. korolkowii\u003c/em\u003e were easily aligned to the main genes involved in the apocarotenoid and flavonoid pathways, including CCD2, UGTs and ALDHs.\u003c/p\u003e \u003cp\u003eAmong the assembled transcripts, we observed potential isoforms of some candidate genes including a flavonoid glucosyltransferase (GT45), UDP-glucose-dependent flavonoid glucosyltransferase (UGT703B1), carotenoid cleavage dioxygenase 2 (CCD2), and beta-carotene hydroxylase in \u003cem\u003eC. chrysanthus\u003c/em\u003e, glucosyltransferase 2 (GLT2) in \u003cem\u003eC. chrysanthus, C. graveolens\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e, aldehyde dehydrogenase 2B4 (ALDH2B4) and glucosyltransferase (kaempferol) in \u003cem\u003eC. chrysanthus\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e (Table S2?).\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eAs more research is conducted on the potential health benefits of saffron value-added metabolites, the demand for these compounds is likely to increase. Finding alternative sources of these metabolites could help ensure a stable supply for the market. While historically the stigmas of some other crocuses like \u003cem\u003eC. cartwrightianus\u003c/em\u003e, the known ancestor of saffron [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e], were used as wild saffron, only \u003cem\u003eC. sativus\u003c/em\u003e stigmas contain the most considerable amount of these metabolites. Since previous investigations confirmed the presence of value-added metabolites in saffron tepals, tepals of other crocuses are a potential source worth exploring. In the present work, using the targeted metabolomics and transcriptomics approaches, four yellow-tepal crocus species were investigated: \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e, \u003cem\u003eC. korolkowii\u003c/em\u003e and \u003cem\u003eC. chrysanthus\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eBy examining the saffron apocarotenoid pathway, we could gain a deeper understanding of the genes involved in producing saffron-like metabolites. The presence of various types of crocins was proved in all four species. However, major crocins of \u003cem\u003eC. sativus\u003c/em\u003e (\u003cem\u003et\u003c/em\u003e-C-2, \u003cem\u003et-\u003c/em\u003eC-3, \u003cem\u003et-\u003c/em\u003eC-4) were only detected in \u003cem\u003eC. chrysanthus\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e3\u003c/span\u003e), which is also supported by our transcriptomic analysis (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). We analyzed three key enzymes of the apocarotenoid pathway, CCDs, ALDHs, and UGTs. We observed that CCD2, the enzyme catalyzing the first committed step of crocin biosynthesis, is present in the tepals from all of the four yellow-tepal crocuses. We also identified transcripts for several ALDH homologs including ALDH2B4, ALDH2C4, ALDH3I1, ALDH5F1, ALDH6B2, ALDH7B4 (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003eA, Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe presence of CCD2, ALDH and crocin-related UGTs transcripts was confirmed in tepals of \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003esissectus\u003c/em\u003e and \u003cem\u003eC. graveolens\u003c/em\u003e, despite no saffron-like crocins being detected in their LC-PDA analysis.\u003c/p\u003e \u003cp\u003eThe absence of detected stigma-like crocins in LC-PDA analysis prompted us to explore a combination of three different hypotheses to understand the underlying reason. The post-harvest degradation of crocins could be considered as the first explanation for the absence of crocins. Since all transcripts of the crocin pathway genes were observed in \u003cem\u003eC. graveolens\u003c/em\u003e, crocins might have become degraded subsequently. The effect of environmental factors on crocins were investigated by HPLC-DAD-MS [\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e] where light and temperature were introduced as the main elements that could degrade crocins after only one week. Due to the detection of other types of crocins in all four species, as well as their detection in saffron stigmas, this hypothesis is not likely.\u003c/p\u003e \u003cp\u003eWe considered the possibility that other compounds and CCD enzymes, rather than those directly involved in the synthesis of saffron-like crocins, may play a more important role in yellow tepal color. Besides crocins, there are various apocarotenoids such as β-cyclocytral and β-ionone which are CCD cleavage products that provide specific aromas and colors ranging from yellow to red in fruits and flowers [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e]. CCD2 belongs to the larger family of CCD enzymes which also includes CCD1, CCD4, CCD7, and CCD8 [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. As it is shown in our phylogenetic analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003eA), candidate genes for CsCCD and CCDs of other crocin-producing plants clustered separately while BdCCD1 clustered in CsCCD. CsCCDs clustered into two groups led by CCD1 and CCD4. The role of CCD1 and CCD4 in providing unique aromas and colors through carotenoid degradation was demonstrated in different plants [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. For instance, in \u003cem\u003eMedicago truncatula\u003c/em\u003e CCD1 is responsible for yellow color [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e, \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e], while in petals of chrysanthemum, CmCCD4a results in colorless compounds in white flowers [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e]. In contrast, the reddish color in citrus fruits is attributed to CitCCD4 [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e]. Moreover, previous studies validated the presence of CCD1a-, CCD1b- and CCD4a/b-encoding genes in the stigmas of \u003cem\u003eC. sativus\u003c/em\u003e and in tepals of \u003cem\u003eC. chrysanthus\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e [\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. CsCCD1 and CsCCD4 convert β-carotene into β-ionone and β-cyclocitral [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e] that play significant roles in the synthesis of aroma and flavor.\u003c/p\u003e \u003cp\u003eIn our study, CCD1 and CCD4 transcripts were also observed in all the studied yellow-tepal species. The presence of these transcripts could contribute to their yellow color along with the presence of crocins (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003e). What is evident is that crocins, although present, are found in low abundance on a dry weight basis compared to saffron stigmas. However, it is essential to consider them as potential commercial sources due to their significantly higher mass compared to saffron stigmas.\u003c/p\u003e \u003cp\u003eThe third and the most probable hypothesis refers to the glycosylation of plant secondary metabolites that is done by UGTs displaying the well-known PSPG (plant secondary-product-glycosyltransferase) box [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. While more than a hundred UGT-encoding genes exist in each plant genome [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e], the saffron-related ones were identified recently [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e] that are shown in phylogenetic tree (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). While UGT74AD1 (also named GL2) allows the formation of crocins with one and two glucose molecules, UGT91P3 does not use crocetin as a substrate and generate crocins with more than two glucose molecules [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. Studies demonstrated the expression patterns of the six identified genes (UGT91K2, UGT91K3, UGT91P3, UGT91P4, UGT91P5 and UGT91P6) of the UGT91 subfamily in several tissues of \u003cem\u003eC. sativus\u003c/em\u003e [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. While only UGT91P3 showed higher expression levels in the stigma, higher expression levels of UGT91K3 were detected in leaves [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e], suggesting the involvement of the latter in the synthesis of other compounds rather than crocins, which are absent in the leaves. As it is shown in the phylogenetic tree, these genes are grouped together (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). UGT91K2, UGT91K3 and UGT91P3 had expression patterns similar to UGT74AD1[\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. Based on the \u003cspan type=\"SmallCaps\" class=\"SmallCaps\" name=\"Emphasis\"\u003eBlast\u003c/span\u003e analysis of the crocin-related UGTs, we observed a higher identity percentage among UGT74AD1 (GL2), UGT91K2, UGT91K3, UGT91P3, UGT91P4, UGT91P5 and UGT91P in \u003cem\u003eC. korolkowii\u003c/em\u003e (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). It could potentially confirm the observed variety of crocins at the metabolome level. In contrast, a lower number of these genes with a high identity percentage were found in \u003cem\u003eC. chrysanthus\u003c/em\u003e and the decreasing trend continued in \u003cem\u003eC. graveolens\u003c/em\u003e and reached to almost non-significant ones in \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e. This observation could also suggest that these species are increasingly distantly related to saffron.\u003c/p\u003e \u003cp\u003eOur metabolome results provide evidence explaining why we still find crocins in yellow tepals, but these differ significantly from those found in the stigmas. The main reason behind these differences could be attributed to variations in the substrate affinity and catalytic properties of the putative crocetin UGTs. Notably, all the detected crocins in the tepals exhibited higher polarity compared to those present in saffron stigmas, leading to their earlier appearance in the LC chromatogram. This higher polarity suggests a greater degree of glycosylation. Investigating glucose availability in the tissues, as glucose is used for glycosylating the crocetin backbone, could be a compelling aspect to explore further. This might shed light on the discrepancy between tepals and stigmas in terms of crocin composition and abundance. This is in line with another study on some yellow-tepal crocuses showing that the crocins in tepals possess higher glycosylation levels in comparison to the crocins in \u003cem\u003eC. sativus\u003c/em\u003e stigmas [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. While crocins in saffron stigmas contain up to five and six glucose molecules [\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e], crocins in yellow-tepal crocuses include up to eight sugar molecules [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. In fact, these highly glycosylated crocins demonstrate differences in retention times indicating different arrangements of the glucose molecules on the ends of the crocetin [\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e]. Moreover, previous studies showed that these crocins are also present in the stigmas of spring crocuses and because of their absence in the stigma of autumn crocuses could be considered a distinguishing factor between spring and autumn crocuses [\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe apocarotenoid HTCC is also glycosylated by UGTs to produce picrocrocin, the precursor of safranal. UGT709G1 catalyzes the HTCC glucosyltransferase reaction [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e] Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Not HTCC, neither picrocrocin nor safranal were detected in yellow-tepals samples (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Nevertheless, at the transcriptome level, the UGT709G1 was identified in almost all samples. While the final enzyme in safranal production is still unknown [\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e, \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e], we could only check HTCC glucosyltransferase that eventually produces picrocrocin. Considering the unknown safranal pathway, several possibilities could be contemplated to explain the absence of safranal despite the presence of UGT709G1. Furthermore, it has been reported that the presence of picrocrocin and safranal in stigmas is greatly influenced by the developmental stage of the flowers [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. Conducting a time-course analysis of tepals could potentially shed light on the ability of crocus tepals to produce these monoterpenoids. Alternatively, the substrates of UGT709G1 might differ in a species- and/or tissue-specific manner, leading to variations in the production of certain compounds.\u003c/p\u003e \u003cp\u003eTo provide a comprehensive view of the presence of saffron value-added metabolites in yellow-tepal crocuses, the presence of saffron stigma-like kaempferols was also investigated. Our results confirmed the presence in \u003cem\u003eC. korolkowii\u003c/em\u003e of two kaempferol derivatives also identified in \u003cem\u003eC. sativus\u003c/em\u003e stigmas. Alongside with kaempferol glycosides, which are the sole flavonoids in saffron stigmas, a higher diversity of flavonoids, such as anthocyanins, quercetin- and flavone-related compounds (luteolin, tricin, acacetin, apigenin, and scutellarein) are found in saffron tepals. [\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e, \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. The detection of multiple compounds with maximum UV-Vis absorption values at 320 nm in all tepal samples, supports the presence of a higher diversity of flavonoids in the yellow-tepal crocuses. Analysis of \u003cem\u003eC. chrysanthus\u003c/em\u003e tepals has proved the existence of kaempferol \u003cem\u003e3-O\u003c/em\u003e glucoside, [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. Other flavonoids such as luteolin-derivatives in \u003cem\u003eC. corsicus\u003c/em\u003e and \u003cem\u003eC. minimus\u003c/em\u003e, and tricin in \u003cem\u003eC. heuffelianus\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e have been reported [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. In our study, we found flavonoid glucosyltransferase transcripts with the highest to the lowest E-value and identity percentage in \u003cem\u003eC. korolkowii\u003c/em\u003e, \u003cem\u003eC. chrysanthus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e and \u003cem\u003eC. flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e, respectively.\u003c/p\u003e \u003cp\u003eUltimately, tepals of \u003cem\u003eC. korolkowii\u003c/em\u003e with three proven types of saffron-like crocins and two kaempferol-derivatives showed the highest similarity to saffron stigmas. It is important to highlight that additional crocin species were found in yellow tepals which despite differing from the ones in saffron, support the metabolic capability of these tissues and species to synthesize crocins. The identified crocins differ in the decorations of the crocetin backbone, likely displaying high glycosylation levels. Although the crocins found in the tepals are different from those in saffron stigmas, they could still be a potential source of pigments that meet industrial demands. The identification of additional CCD transcripts in the yellow crocuses may also be involved in the synthesis of non-crocin pigments, which could also be of commercial interest. Moreover, our results support the presence of diverse flavonoid in tepals, which could be exploited as a source of natural products for nutritional and medicinal purposes.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003ePrevious observations noted that saffron tepals contain metabolites similar to the ones in stigma, leading to recent exploration of tepals as a potential compound source. Our study focused on the metabolomic and transcriptomic analysis of saffron stigma-like apocarotenoids and kaempferol-derivatives in the tepals of four common yellow-tepal crocuses.\u003c/p\u003e \u003cp\u003eAmong these crocuses, the metabolites in \u003cem\u003eC. korolkowii\u003c/em\u003e tepals displayed the highest similarity to those in saffron stigmas, containig three verified stigma-like crocins and two kaempferol-derivatives. However, most crocins in yellow tepals differed structurally from those in saffron, featuring variations in crocetin backbone decorations and increased glycosylation. Despite these distinctions, the identified crocins in tepals remain a potential source for industrial pigments.\u003c/p\u003e \u003cp\u003eMoreover, our results indicated a diverse array of flavonoids in tepals, presenting an additional opportunity for exploration. These flavonoids hold promise as a source for the development of natural products with applications in nutrition and medicine. In summary, while crocins in yellow-tepal crocus tepals differ from those in saffron stigmas, they still hold potential for industrial pigments. Additionally, the presence of diverse flavonoids in tepals expands the potential applications of these crocuses, offering opportunities for the development of natural products with nutritional and medicinal benefits.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003ePlant materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eStigmas of the autumn-flowering saffron crocus (\u003cem\u003eC. sativus\u003c/em\u003e) along with tepal samples of four species each of the yellow, spring-flowering species \u003cem\u003eC. flavus\u0026nbsp;\u003c/em\u003esubsp.\u003cem\u003e\u0026nbsp;dissectus\u003c/em\u003e, \u003cem\u003eC. graveolens\u003c/em\u003e and \u003cem\u003eC. korolkowii\u003c/em\u003e were grown and collected in IPK Gatersleben greenhouses while\u003cem\u003e\u0026nbsp;C. chrysanthus\u003c/em\u003e was provided by Vladimir Randjelovic from the University of Ni\u0026scaron;. Samples for metabolite and transcriptome analyses were collected with two different methods. For metabolic analysis, tepals and stigmas were fast-dried in silica gel and stored dry till their analysis. For the transcriptome analyses, tissue samples of tepals and stigmas were stabilized in RNA\u003cem\u003elater\u003c/em\u003e (Qiagen, Germany) to preserve RNA integrity during storage.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMetabolite extraction and analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSoluble semi-polar metabolites were extracted from tepal and stigma samples. All extraction steps were performed at room temperature under darkness. Firstly, samples were transferred into 2-ml Eppendorf tubes and weighted (5 to 15 mg). These were added with 100 \u0026micro;l of 50% (v/v) methanol per mg of material. Samples were grinded using a Precellys homogeneizer (Bertin Instruments, France) using 1-1.2 mm-diameter zirconium silicate beads (M\u0026uuml;hlmeier GmbH, Germany) with four cycles of 20 s each, at a frequency of 6000 Hz. Samples were vortexed and extracted through continuous agitation at 500 rpm and room temperature for 20 min under dark conditions (Eppendorf ThermoMixer, Germany). Samples were centrifuged at 17,000 g for 5 min and the supernatant was transferred into a new tube. The pellet was re-extracted two more times, each by adding 100 \u0026micro;l of 50% (v/v) methanol per mg of material, followed by vortexing, shaking (20 min, 500 rpm), and centrifugation. In each step, the supernatants were combined with the first one. Extracts were covered with Argon and stored at -20 \u0026deg;C prior to analysis.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSemi-polar metabolites were analyzed by reversed phase ultra-performance liquid chromatography coupled with photodiode array detection (RP-UPLC-PDA), using an Acquity UPLC system (Waters, Germany) equipped with an Acquity UPLC PDA e\u0026lambda; detector (Waters, Germany). Extracts were first centrifuged at 17,000 \u003cem\u003eg\u003c/em\u003e for 10 min, and 100 \u0026micro;l of supernatant were transferred into glass vials for analysis. Sample injection volumes of 5 \u0026micro;l were employed for UPLC analysis, using a partial-loop with needle-overfill (PLNO) injection mode with a 10 \u0026micro;l loop. Compounds were separated in an Acquity UPLC BEH phenyl column (130 \u0026Aring;, 2.1 x 100 mm, 1.7 \u0026mu;m; Waters, Germany) combined with an Acquity UPLC BEH phenyl VanGuard pre-column (130 \u0026Aring;, 2.1 x 5 mm, 1.7 \u0026mu;m; Waters, Germany) using the following gradient: from 10 to 36% of solvent B for the first 3.9 min, isocratic hold from 3.9 to 5 min, from 36 to 97% B from 5 to 6 min, and an isocratic hold from 6 to 7 min to clean the column; after each run, the column was equilibrated to the starting conditions (10% B); solvent A was LC-MS grade water (CHEMSOLUTE, Th. Geyer, Germany) with 0.5% (v/v) formic acid, and solvent B was LC-MS grade acetonitrile (CHEMSOLUTE, Th. Geyer, Germany) with 0.5% (v/v) formic acid. The column temperature was maintained at 35\u0026deg;C and the flow was set to 500 \u0026mu;l min\u003csup\u003e-1\u003c/sup\u003e. PDA-detection was performed in a range between 210 and 800 nm, at a resolution of 1.2 nm and a sampling rate of 20 points s\u003csup\u003e-1\u003c/sup\u003e. Within the major \u003cem\u003eCrocus\u0026nbsp;\u003c/em\u003esemi-polar metabolites, crocins, picrocrocin, and flavonoids and safranal, were detected at 440, 250, and 320 nm, respectively. The identity of crocin (\u003cem\u003etrans\u003c/em\u003e-crocetin-digentibiosyl ester), picrocrocin, and safranal, was confirmed by retention time and UV absorption spectra with commercial standards. Processing and analysis of the acquired PDA spectra was done with the Empower 3 software (Waters, Germany).\u003c/p\u003e\n\u003cp\u003eTo generate a reliable reference composition, the extracts from \u003cem\u003eC. sativus\u003c/em\u003e stigmas were further analyzed via ESI-UHR-QTOF-MS (ElectroSpray Ionization-Ultra-High-Resolution-Quadrupole Time of Flight-Mass Spectrometry) by coupling a maXis Impact ESI-QTOF MS (Bruker Daltonik GmbH; Germany) to the RP-UPLC-PDA system. The MS analyses were performed in positive and negative ionization modes. Small molecules (\u003cem\u003ei.e.,\u003c/em\u003e safranal, picrocrocin, flavonoids; 50-1000 \u003cem\u003em/z\u003c/em\u003e) were analyzed using the MS1 and MS/MS settings for barley and sunflower phenylpropanoids described in Garibay-Hern\u0026aacute;ndez \u003cem\u003eet al.\u003c/em\u003e 2021. Large molecules (\u003cem\u003ei.e.,\u003c/em\u003e crocins; 50-1500 \u003cem\u003em/z\u003c/em\u003e) were analyzed with the MS1 settings in positive ionization mode for anthocyanins as described in Garibay-Hern\u0026aacute;ndez \u003cem\u003eet al.\u003c/em\u003e 2021. The MS/MS analysis of large molecules in positive mode was performed in auto MS/MS using CID (Collision-Induced Dissociation) with the following settings: absolute area threshold: 5000 counts; exclusion activation: 15 spectra; exclusion release: 30\u0026thinsp;s; collision energy values (z = 1, 2, 3; isolation mass = 500; width = 8): 35, 25, 20 eV; collision energy values (z = 1, 2, 3; isolation mass = 1000; width = 10): 50, 40, 35\u0026thinsp;eV. The analysis of large molecules in negative mode was done with the following MS1 settings: 50\u0026ndash;1000\u0026thinsp;\u003cem\u003em/z\u003c/em\u003e; capillary voltage: 3.5 kV; nebulizer: 3 bar; dry gas: 8 l min\u003csup\u003e\u0026minus;\u003c/sup\u003e\u003csup\u003e1\u003c/sup\u003e; dry temperature: 200 \u0026deg;C; hexapole RF (Ratio Frequency) voltage: 150 Vpp (V peak-to-peak); funnel 1 RF: 400 Vpp; funnel 2 RF: 400 Vpp; pre-pulse storage time: 8\u0026nbsp;\u0026mu;s; transfer time: 60\u0026thinsp;\u0026mu;s; low mass: 40\u0026thinsp;\u003cem\u003em/z\u003c/em\u003e; collision cell RF: 800 Vpp; collision energy: 8 eV. The MS/MS analysis of large molecules in negative mode was done in auto MS/MS using CID as follows: absolute area threshold: 5000 counts; exclusion activation: 2 spectra; exclusion release: 12\u0026thinsp;s; collision energy values (z = 1; isolation mass = 500; width = 6): 20 eV; collision energy values (z = 1; isolation mass = 1000; width = 8): 20\u0026thinsp;eV. The Compass HyStar 3.2 SR2 software (Bruker Daltonik GmbH, Germany) was used to operate and coordinate LC-PDA-MS data acquisition. The analysis of MS data was performed using the Compass DataAnalysis 4.4 SR1 package (Bruker Daltonik GmbH).\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRNA extraction and transcriptome sequencing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTotal RNA was extracted using the Qiagen RNeasy plant mini kit. Total RNA purity and concentration were determined by the Qubit 1 2.0 Flurometer (Life Technologies, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, USA). A minimum amount of 50 ng\u0026nbsp;\u0026mu;l\u003csup\u003e-1\u003c/sup\u003e RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using the llumina TruSeq RNA Sample Prep Kit v2 for Illumina sequencing platforms.\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTranscriptome assembly and annotation\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe library preparations were sequenced on an Illumina Novaseq 6000 platform and paired-end reads were generated. In the quality control step, raw reads of fastq format were first processed through trimmomatic\u0026nbsp;[28]. All the downstream analyses were based on clean data with high quality. The transcriptome was assembled using Trinity-v2.9.0\u0026nbsp;[29]\u0026nbsp;with min_kmer_cov set to 2 by default and all other parameters set default.\u003c/p\u003e\n\u003cp\u003eGenes linked to crocin and kaempferol biosynthesis in \u003cem\u003eC. sativus\u003c/em\u003e were selected from\u0026nbsp;[7, 17, 30\u0026ndash;33]. These genes are listed in Table 1 and were subsequently used as a database in TBlastN and potential orthologs of these genes in the Trinity transcriptome assemblies were searched using the TBlastN parameters with e-value cut-of 0.01. Genes were then selected as potential full-length orthologs if they contained an ORF that had an overlap with at least 90% of the reference gene. Shorter sequences were listed as partial genes (Table S2).\u003c/p\u003e\n\u003cp\u003ePhylogenetic trees for three main genes (CCDs, ALDHs and UGTs) are calculated using MEGA3 with ClustalW alignment along with maximum likelihood with bootstrap 500 to generate the trees. \u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eALDH: Aldehyde dehydrogenase\u003c/p\u003e\n\u003cp\u003eCHY-\u0026beta;: \u0026beta;-Carotene hydrolase\u003c/p\u003e\n\u003cp\u003eBLAST: Basic Local Alignment Search Tool\u003c/p\u003e\n\u003cp\u003eCCD: Carotenoid Cleavage Dioxygenase\u003c/p\u003e\n\u003cp\u003ec-C: Cis-crocetin\u003c/p\u003e\n\u003cp\u003eESI-UHR-QTOF-MS: ElectroSpray Ionization-Ultra-High-Resolution-Quadrupole Time of Flight-Mass Spectrometry\u003c/p\u003e\n\u003cp\u003eGT: Glycosyltransferase\u003c/p\u003e\n\u003cp\u003eHPLC-DAD-MS: High-Performance Liquid Chromatography coupled to Diode Array Detection and Mass Spectrometry\u003c/p\u003e\n\u003cp\u003eHTCC: 2,6,6-Trimethyl-4-hydroxy-1-carboxyaldehyde-1-cyclohexene\u003c/p\u003e\n\u003cp\u003eK1: Kaempferol 3-O-sophoroside 7-O-glucoside\u003c/p\u003e\n\u003cp\u003eK2: Kaempferol 3-O-sophoroside\u003c/p\u003e\n\u003cp\u003eLC: Liquid Chromatography\u003c/p\u003e\n\u003cp\u003eLC-PDA: Liquid Chromatography coupled to Photodiode Array\u003c/p\u003e\n\u003cp\u003eLCYB: Lycopene \u0026beta;-cyclase\u003c/p\u003e\n\u003cp\u003eLC-MS: Liquid Chromatography-Mass Spectrometry\u003c/p\u003e\n\u003cp\u003eMEP: Methylerythritol phosphate\u003c/p\u003e\n\u003cp\u003eMS: Mass Spectrometry\u003c/p\u003e\n\u003cp\u003eMVA: Mevalonate pathways\u003c/p\u003e\n\u003cp\u003eORF: Open Reading Frame\u003c/p\u003e\n\u003cp\u003ePDA: Photodiode Array Detection\u003c/p\u003e\n\u003cp\u003ePSY: Phytoene synthase\u003c/p\u003e\n\u003cp\u003ePSPG: Plant Secondary-Product-Glycosyltransferase\u003c/p\u003e\n\u003cp\u003eQTOF: Quadrupole Time of Flight\u003c/p\u003e\n\u003cp\u003eRNA: Ribonucleic Acid\u003c/p\u003e\n\u003cp\u003eRP-UPLC-PDA: Reversed Phase Ultra-Performance Liquid Chromatography coupled with Photodiode Array Detection\u003c/p\u003e\n\u003cp\u003et-C: Trans-crocetin\u003c/p\u003e\n\u003cp\u003eTBLASTN: Translated Basic Local Alignment Search Tool with Nucleotides\u003c/p\u003e\n\u003cp\u003eUGTs: Uridine Diphosphate Glucosyltransferases\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe RNA-seq data for this study can be accessed at the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra) under accession number PRJNA926329.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe extend our thanks to Elena Brueckner for her expertise in metabolic profiling, Dr. Reinhard M. Fritsch for providing \u003cem\u003eC\u003c/em\u003e.\u003cem\u003e\u0026nbsp;flavus\u003c/em\u003e subsp. \u003cem\u003edissectus\u003c/em\u003e samples, and Prof. Vladimir Randjelovic for sharing \u003cem\u003eC\u003c/em\u003e.\u003cem\u003e\u0026nbsp;chrysanthus\u003c/em\u003e specimens. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research was funded by Deutsche Forschungsgemeinschaft (DFG) for the financial support of part of this work.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eZN, FRB, HPM, and BU planned the experiment; ZN, SSKS, MHWS analysed of the transcriptome data; AGH performed the metabolite analysis. The first version of the manuscript was prepared by SSKS. All authors participated in reviewing and approving the final version of the paper.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eKazemi-Shahandashti S-S, Mann L, El-Nagish A, Harpke D, Nemati Z, Usadel B, et al. Ancient Artworks and Crocus Genetics Both Support Saffron\u0026rsquo;s Origin in Early Greece. Front Plant Sci. 2022;13:834416.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDiretto G, Ahrazem O, Rubio-Moraga \u0026Aacute;, Fiore A, Sevi F, Argando\u0026ntilde;a J, et al. UGT709G1: a novel uridine diphosphate glycosyltransferase involved in the biosynthesis of picrocrocin, the precursor of safranal in saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e). New Phytol. 2019;224:725\u0026ndash;40.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZeka K, Ruparelia KC, Continenza MA, Stagos D, Vegli\u0026ograve; F, Arroo RRJ. Petals of \u003cem\u003eCrocus sativus\u003c/em\u003e L. as a potential source of the antioxidants crocin and kaempferol. Fitoterapia. 2015;107:128\u0026ndash;34.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBisti S, Maccarone R, Falsini B. Saffron and retina: neuroprotection and pharmacokinetics. Vis Neurosci. 2014;31:355\u0026ndash;61.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAbdullaev FI. Cancer Chemopreventive and Tumoricidal Properties of Saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e L). Exp Biol Med (Maywood). 2002;227:20\u0026ndash;5.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRubio Moraga A, Ahrazem O, Rambla JL, Granell A, G\u0026oacute;mez G\u0026oacute;mez L. Crocins with high levels of sugar conjugation contribute to the yellow colours of early-spring flowering crocus tepals. PLoS ONE. 2013;8:e71946.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJi A, Jia J, Xu Z, Li Y, Bi W, Ren F, et al. Transcriptome-Guided Mining of Genes Involved in Crocin Biosynthesis. Front Plant Sci. 2017;8:518.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen Y, Zhang H, Li Y-X, Cai L, Huang J, Zhao C, et al. Crocin and geniposide profiles and radical scavenging activity of gardenia fruits (\u003cem\u003eGardenia jasminoides\u003c/em\u003e Ellis) from different cultivars and at the various stages of maturation. Fitoterapia. 2010;81:269\u0026ndash;73.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIshizuka F, Shimazawa M, Umigai N, Ogishima H, Nakamura S, Tsuruma K, et al. Crocetin, a carotenoid derivative, inhibits retinal ischemic damage in mice. Eur J Pharmacol. 2013;703:1\u0026ndash;10.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSong Y-N, Wang Y, Zheng Y-H, Liu T-L, Zhang C. Crocins: A comprehensive review of structural characteristics, pharmacokinetics and therapeutic effects. Fitoterapia. 2021;153:104969.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaba SA, Mohiuddin T, Basu S, Swarnkar MK, Malik AH, Wani ZA, et al. Comprehensive transcriptome analysis of \u003cem\u003eCrocus sativus\u003c/em\u003e for discovery and expression of genes involved in apocarotenoid biosynthesis. BMC Genomics. 2015;16:698.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eN\u0026oslash;rb\u0026aelig;k R, Brandt K, Nielsen JK, \u0026Oslash;rgaard M, Jacobsen N. Flower pigment composition of \u003cem\u003eCrocus\u003c/em\u003e species and cultivars used for a chemotaxonomic investigation. Biochem Syst Ecol. 2002;30:763\u0026ndash;91.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLee G-A, Choi K-C, Hwang K-A. Kaempferol, a phytoestrogen, suppressed triclosan-induced epithelial-mesenchymal transition and metastatic-related behaviors of MCF-7 breast cancer cells. Environ Toxicol Pharmacol. 2017;49:48\u0026ndash;57.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoraga AR, Rambla JL, Ahrazem O, Granell A, G\u0026oacute;mez-G\u0026oacute;mez L. Metabolite and target transcript analyses during \u003cem\u003eCrocus sativus\u003c/em\u003e stigma development. Phytochemistry. 2009;70:1009\u0026ndash;16.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaito K, Yonekura-Sakakibara K, Nakabayashi R, Higashi Y, Yamazaki M, Tohge T, et al. The flavonoid biosynthetic pathway in Arabidopsis: Structural and genetic diversity. Plant Physiol Biochem. 2013;72:21\u0026ndash;34.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDuan L, Ding W, Liu X, Cheng X, Cai J, Hua E, et al. Biosynthesis and engineering of kaempferol in \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e. Microb Cell Fact. 2017;16:165.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWinterhalter P, Straubinger M. Saffron-Renewed interest in an ancient spice. Food Reviews Int. 2000;16:39\u0026ndash;59.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRavishankar B, Shukla VJ. Indian systems of medicine: a brief profile. Afr J Tradit Complement Altern Med. 2007;4:319\u0026ndash;37.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHosseini A, Razavi BM, Hosseinzadeh H. Saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e) petal as a new pharmacological target: a review. Iran J Basic Med Sci. 2018;21:1091\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNaidis P, Lykidou S, Mischopoulou M, Vouvoudi E, Nikolaidis NF. Study of the dyeing properties of saffron and ultrafiltrated saffron powders, as colourants for natural and synthetic fibres. Color Technol. 2023;139:565\u0026ndash;77.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHarpke D, Carta A, Tomović G, Ranđelović V, Ranđelović N, Blattner FR, et al. Phylogeny, karyotype evolution and taxonomy of \u003cem\u003eCrocus\u003c/em\u003e series Verni (Iridaceae). Plant Syst Evol. 2015;301:309\u0026ndash;25.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHarpke D, Kerndorff H, Pasche E, Peruzzi L. Neotypification of the name \u003cem\u003eCrocus biflorus\u003c/em\u003e Mill. (Iridaceae) and its consequences in the taxonomy of the genus. Phytotaxa. 2016;260:131.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePeruzzi L. \u003cem\u003eCrocus heuffelianus\u003c/em\u003e (Iridaceae), a new record for the Italian flora. Phytotaxa. 2016;261:291.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRaca I, Harpke D, Shuka L, Ranđelović V. A new species of \u003cem\u003eCrocus\u003c/em\u003e ser. \u003cem\u003eVerni\u003c/em\u003e (Iridaceae) with 2 \u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;12 chromosomes from the Balkans. Plant Biosystems - Int J Dealing all Aspects Plant Biology. 2022;156:36\u0026ndash;42.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShuka D, Kerndorff H, Waminal N, Ranđelovič V, Harpke D. Two new species of \u003cem\u003eCrocus\u003c/em\u003e series \u003cem\u003eScardici\u003c/em\u003e from the Western Balkans. preprint. Evolutionary Biology; 2022.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMohtashami L, Amiri MS, Ramezani M, Emami SA, Simal-Gandara J. The genus \u003cem\u003eCrocus\u003c/em\u003e L.: A review of ethnobotanical uses, phytochemistry and pharmacology. Ind Crops Prod. 2021;171:113923.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhu C, Bai C, Sanahuja G, Yuan D, Farr\u0026eacute; G, Naqvi S, et al. The regulation of carotenoid pigmentation in flowers. Arch Biochem Biophys. 2010;504:132\u0026ndash;41.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114\u0026ndash;20.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGrabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644\u0026ndash;52.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaba SA, Ashraf N. Apocarotenoids of Crocus sativus L: From biosynthesis to pharmacology. Singapore: Springer Singapore; 2016.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAhrazem O, Diretto G, Argando\u0026ntilde;a J, Rubio-Moraga \u0026Aacute;, Julve JM, Orz\u0026aacute;ez D, et al. Evolutionarily distinct carotenoid cleavage dioxygenases are responsible for crocetin production in \u003cem\u003eBuddleja davidii\u003c/em\u003e. J Exp Bot. 2017;68:4663\u0026ndash;77.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDemurtas OC, Frusciante S, Ferrante P, Diretto G, Azad NH, Pietrella M, et al. Candidate Enzymes for Saffron Crocin Biosynthesis Are Localized in Multiple Cellular Compartments. Plant Physiol. 2018;177:990\u0026ndash;1006.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eL\u0026oacute;pez-Jimenez AJ, Frusciante S, Niza E, Ahrazem O, Rubio-Moraga \u0026Aacute;, Diretto G, et al. A New Glycosyltransferase Enzyme from Family 91, UGT91P3, Is Responsible for the Final Glucosylation Step of Crocins in Saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e L). Int J Mol Sci. 2021;22:8815.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHughes J, Hughes MA. Multiple secondary plant product UDP-glucose glucosyltransferase genes expressed in cassava (Manihot esculenta Crantz) cotyledons. DNA Seq. 1994;5:41\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNemati Z, Blattner FR, Kerndorff H, Erol O, Harpke D. Phylogeny of the saffron-crocus species group, \u003cem\u003eCrocus\u003c/em\u003e series \u003cem\u003eCrocus\u003c/em\u003e (Iridaceae). Mol Phylogenet Evol. 2018;127:891\u0026ndash;7.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNemati Z, Harpke D, Gemicioglu A, Kerndorff H, Blattner FR. Saffron (Crocus sativus) is an autotriploid that evolved in Attica (Greece) from wild \u003cem\u003eCrocus cartwrightianus\u003c/em\u003e. Mol Phylogenet Evol. 2019;136:14\u0026ndash;20.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSuchareau M, Bordes A, Lem\u0026eacute;e L. Improved quantification method of crocins in saffron extract using HPLC-DAD after qualification by HPLC-DAD-MS. Food Chem. 2021;362:130199.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54:733\u0026ndash;49.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTimmins JJB, Kroukamp H, Walker RSK, Pretorius IS, Paulsen IT. Comparative Evaluation of Secreted Plant Carotenoid Cleavage Dioxygenase 1 (CCD1) Enzymes in \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e. Fermentation. 2022;8:395.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOhmiya A, Kato M, Shimada T, Nashima K, Kishimoto S, Nagata M. Molecular Basis of Carotenoid Accumulation in Horticultural Crops. Hort J. 2019;88:135\u0026ndash;49.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKlingner A, Hundeshagen B, Kernebeck H, Bothe H. Localization of the yellow pigment formed in roots of gramineous plants colonized by arbuscular fungi. Protoplasma. 1995;185:50\u0026ndash;7.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFester T, Schmidt D, Lohse S, Walter MH, Giuliano G, Bramley PM, et al. Stimulation of carotenoid metabolism in arbuscular mycorrhizal roots. Planta. 2002;216:148\u0026ndash;54.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOhmiya A, Kishimoto S, Aida R, Yoshioka S, Sumitomo K. Carotenoid cleavage dioxygenase (CmCCD4a) contributes to white color formation in chrysanthemum petals. Plant Physiol. 2006;142:1193\u0026ndash;201.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa G, Zhang L, Matsuta A, Matsutani K, Yamawaki K, Yahata M, et al. Enzymatic formation of β-citraurin from β-cryptoxanthin and Zeaxanthin by carotenoid cleavage dioxygenase4 in the flavedo of citrus fruit. Plant Physiol. 2013;163:682\u0026ndash;95.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRubio A, Rambla JL, Santaella M, G\u0026oacute;mez MD, Orzaez D, Granell A, et al. Cytosolic and plastoglobule-targeted carotenoid dioxygenases from Crocus sativus are both involved in beta-ionone release. J Biol Chem. 2008;283:24816\u0026ndash;25.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLe Roy J, Huss B, Creach A, Hawkins S, Neutelings G. Glycosylation Is a Major Regulator of Phenylpropanoid Availability and Biological Activity in Plants. Front Plant Sci. 2016;7:735.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePfister S, Meyer P, Steck A, Pfander H. Isolation and Structure Elucidation of Carotenoid\u0026thinsp;\u0026ndash;\u0026thinsp;Glycosyl Esters in Gardenia Fruits (\u003cem\u003eGardenia jasminoides\u003c/em\u003e Ellis) and Saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e Linne). J Agric Food Chem. 1996;44:2612\u0026ndash;5.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDufresne C, Cormier F, Dorion S, Niggli UA, Pfister S, Pfander H. Glycosylation of encapsulated crocetin by a Crocus sativus L. cell culture. Enzym Microb Technol. 1999;24:453\u0026ndash;62.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTarantilis PA, Tsoupras G, Polissiou M. Determination of saffron (Crocus sativus L.) components in crude plant extract using high-performance liquid chromatography-UV-visible photodiode-array detection-mass spectrometry. J Chromatogr A. 1995;699:107\u0026ndash;18.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCastillo R, Fern\u0026aacute;ndez J-A, G\u0026oacute;mez-G\u0026oacute;mez L. Implications of carotenoid biosynthetic genes in apocarotenoid formation during the stigma development of \u003cem\u003eCrocus sativus\u003c/em\u003e and its closer relatives. Plant Physiol. 2005;139:674\u0026ndash;89.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKanakis CD, Daferera DJ, Tarantilis PA, Polissiou MG. Qualitative determination of volatile compounds and quantitative evaluation of safranal and 4-hydroxy-2,6,6-trimethyl-1-cyclohexene-1-carboxaldehyde (HTCC) in Greek saffron. J Agric Food Chem. 2004;52:4515\u0026ndash;21.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoras B, Loffredo L, Rey S. Quality assessment of saffron (\u003cem\u003eCrocus sativus\u003c/em\u003e L.) extracts via UHPLC-DAD-MS analysis and detection of adulteration using gardenia fruit extract (\u003cem\u003eGardenia jasminoides\u003c/em\u003e Ellis). Food Chem. 2018;257:325\u0026ndash;32.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi C-Y, Wu T-S. Constituents of the pollen of Crocus sativus L. and their tyrosinase inhibitory activity. Chem Pharm Bull (Tokyo). 2002;50:1305\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMykhailenko O, Kovalyov V, Goryacha O, Ivanauskas L, Georgiyants V. Biologically active compounds and pharmacological activities of species of the genus \u003cem\u003eCrocus\u003c/em\u003e: A review. Phytochemistry. 2019;162:56\u0026ndash;89.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZengin G, Aumeeruddy MZ, Diuzheva A, Jekő J, Czi\u0026aacute;ky Z, Yıldıztugay A, et al. A comprehensive appraisal on Crocus chrysanthus (Herb.) Herb. flower extracts with HPLC-MS/MS profiles, antioxidant and enzyme inhibitory properties. J Pharm Biomed Anal. 2019;164:581\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHarborne JB, Williams CA. 6-Hydroxyflavones and Other Flavonoids of \u003cem\u003eCrocus\u003c/em\u003e. Z f\u0026uuml;r Naturforschung C. 1984;39:18\u0026ndash;23.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-plant-biology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pbio","sideBox":"Learn more about [BMC Plant Biology](http://bmcplantbiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/pbio/default.aspx","title":"BMC Plant Biology","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Saffron-like compounds, Yellow-tepal crocuses, Metabolomics, Transcriptomics, Crocins, Flavonoids, Alternative sources","lastPublishedDoi":"10.21203/rs.3.rs-3948859/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-3948859/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground\u003c/strong\u003e:\u003c/p\u003e\n\u003cp\u003eThe increasing demand for saffron metabolites in various commercial industries, including medicine, food, cosmetics, and dyeing, is driven by the discovery of their diverse applications. Saffron, derived from \u003cem\u003eCrocus sativus\u003c/em\u003estigmas, is the most expensive spice, and there is a need to explore additional sources to meet global consumption demands. In this study, we focused on yellow-flowering crocuses and examined their tepals to identify saffron-like compounds.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResults\u003c/strong\u003e:\u003c/p\u003e\n\u003cp\u003eThrough metabolomic and transcriptomic approaches, our investigation provides valuable insights into the biosynthesis of compounds in yellow-tepal crocuses that are similar to those found in saffron. The results of our study support the potential use of yellow-tepal crocuses as a source of various crocins (crocetin glycosylated derivatives) and flavonoids.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConclusions\u003c/strong\u003e:\u003c/p\u003e\n\u003cp\u003eOur findings suggest that yellow-tepal crocuses have the potential to serve as a viable excessive source of some saffron metabolites. The identification of crocins and flavonoids in these crocuses highlights their suitability for meeting the demands of various industries that utilize saffron compounds. Further exploration and utilization of yellow-tepal crocuses could contribute to addressing the growing global demand for saffron-related products.\u003c/p\u003e","manuscriptTitle":"Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-02-16 04:35:18","doi":"10.21203/rs.3.rs-3948859/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-03-19T07:07:01+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-03-08T09:30:38+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"9c89c2bb-4a02-41ee-aa1d-c3676263bb50","date":"2024-02-28T23:16:39+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-02-28T13:45:23+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2024-02-14T18:46:49+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-02-14T18:43:10+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-02-14T18:42:23+00:00","index":"","fulltext":""},{"type":"submitted","content":"BMC Plant Biology","date":"2024-02-11T16:01:08+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-plant-biology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pbio","sideBox":"Learn more about [BMC Plant Biology](http://bmcplantbiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/pbio/default.aspx","title":"BMC Plant Biology","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"c67cd020-a7ee-4216-8105-08acf854ae1e","owner":[],"postedDate":"February 16th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-05-09T03:57:58+00:00","versionOfRecord":{"articleIdentity":"rs-3948859","link":"https://doi.org/10.1186/s12870-024-05036-1","journal":{"identity":"bmc-plant-biology","isVorOnly":false,"title":"BMC Plant Biology"},"publishedOn":"2024-05-07 03:57:57","publishedOnDateReadable":"May 7th, 2024"},"versionCreatedAt":"2024-02-16 04:35:18","video":"","vorDoi":"10.1186/s12870-024-05036-1","vorDoiUrl":"https://doi.org/10.1186/s12870-024-05036-1","workflowStages":[]},"version":"v1","identity":"rs-3948859","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-3948859","identity":"rs-3948859","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}