BSA-seq integrated with transcriptomics and metabolomics revealing the candidate genes associated with safflower colors and flavonoid glycosides biosynthesis

Abstract

Safflower is featured with time-honored medical and economic values and developing into diverse phenotypic and genetic variations. In this study, to explore the critical genes associated with color phenotypes and flavonoid derivatives biosynthesis of safflower, BSA-seq, conflated with transcriptomic and metabolic methods were performed in two extreme colors (yellow and white) in the population of ‘ZHH0119’ and ‘XHH007.’ After crossing two parent plants reciprocally, the F₄ generation of two accessions were used to construct near-isogenic gene pools for the two extreme traits of yellow and white safflower. BSA-seq results located five QTLs regions on chromosomes 2, 8, 9, 10, and 12 including 6 CtPALs, 3 CtC4Hs, 2 Ct4CLs, 1 CtCHS, 32 CtUGTs, and 70 CtCYPs, which tied to the yellow color phenotype of safflower. Through transcriptome analysis of two accessions and at different flowering stages, 1 CtPAL, 5 CtC4Hs, 4 CtCHSs, 3 CtCHIs, 3 CtFLSs, 48 CtUGTs, 51 CtCYPs, and 75 transcription factors were revealed as significantly upregulated in the yellow accession compared to the white. Integrated analysis identified eight CtUGTs ( CtUGT50-57) that exhibited significant positive correlations with chalcone glycosides of yellow safflower. Based on functional characterization, CtUGT52 was found to boost Hydroxysafflor yellow A (HSYA) content in yellow safflower which possessing substrate promiscuity (chalones, flavonols, and flavonoids) and catalytic promiscuity (flavonols and flavonoids), revealing its vital role in the HSYA biosynthesis through transgenic overexpression. Combining catalytic mechanism verification of CtUGT52 towards phloretin, kaempferol, and luteolin, our study to some extent, elucidated the modification function of CtUGTs for flavonoid aglycones in the flavonoid biosynthesis pathway of safflower.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (81973421) and the National Key R&D Program of China (2019YFC1711100).

Author contributions

H.W. and B.H. conducted the experiments and data analysis, wrote the manuscript. S.Q. contributed to the protein function verification experiments. Y.G. and X.D. guided quantitative analysis and provided valuable advice. M.G. conceived this project, reviewed, and edited. All authors provided final approval for publication.

Data availability

The transcriptome data have been uploaded on NCBI under BioProject: PRJNA1262321.

Conflicts of interest statement

The authors declare that they have no conflicts of interest.

References

[1] Zafari M, Ebadi A, Jahanbakhsh S, et al. Safflower (Carthamus tinctorius) biochemical properties, yield, and oil content affected by 24-epibrassinosteroid and genotype under drought stress . J Agric Food Chem. 2020; 68: 6040-7.

[2] Cheng H, Yang CL, Ge PL, et al. Genetic diversity, clinical uses, and phytochemical and pharmacological properties of safflower (Carthamus tinctorius L.): an important medicinal plant . Front Pharmacol. 2024; 15: 1374680.

[3] Wang YQ, Zhang SS, Ni HL, et al. Autophagy is involved in the neuroprotective effect of nicotiflorin. J Ethnopharmacol. 2021; 278: 114279.

[4] Yan Z, Alimu R, Wan J, et al. Composition of major quinochalcone hydroxysafflor yellow A and anhydrosafflor yellow B is associated with colour of safflower (Carthamus tinctorius) during colour-transition but not with overall antioxidant capacity: a study on 144 cultivars . Food Res Int. 2022; 162: 112098.

[5] He BX, Zhang YJ, Wang LN, et al. Both two transcripts promoting the accumulation of the flavonoid profiles in overexpressed transgenic safflower. Front Plant Sci. 2022; 13: 833811.

[6] Chen Y, Li MF, Wen JY, et al. Pharmacological activities of safflower yellow and its clinical applications. Evid Based Complement Alternat Med. 2022; 2022: 1-24.

[7] Wang Y, Li XY, Deng FF, et al. Hydroxy-safflower yellow a alleviates osteoporosis in ovariectomized rat model by inhibiting carbonic anhydrase 2 activity. Front Pharmacol. 2021; 12: 734539.

[8] Rong JJ, Li CY, Zhang Q, et al. Hydroxysafflor yellow A inhibits endothelial cell ferroptosis in diabetic atherosclerosis mice by regulating miR-429/SLC7A11. Pharm Biol. 2023; 61: 404-15.

[9] Qin ED, Zhan W, Yan XC, et al. Genetic diversity of main agronomic triats in 482 safflower accessions. Chin J Oil Crop Sci. 2022; 44: 78-86.

[10] Wang YQ, Jin JQ, Zhang R, et al. Association analysis of BSA-seq, BSR-seq, and RNA-seq reveals key genes involved in purple leaf formation in a tea population (Camellia sinensis) . Hortic Res. 2024; 11: uhae191.

[11] Wang ZL, Wang HT, Chang GW, et al. Elucidation of the biosynthetic pathway of hydroxysafflor yellow A. Nat Commun. 2025; 16: 4489.

[12] Yin X, Xiang Y, Huang FQ, et al. Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis . Plant Biotechnol J. 2023; 21: 2209-23.

[13] Sardouei-Nasab S, Nemati Z, Mohammadi-Nejad G, et al. Phylogenomic investigation of safflower (Carthamus tinctorius) and related species using genotyping-by-sequencing (GBS) . Sci Rep. 2023; 13: 6212.

[14] Yildiz M, Altaf MT, Baloch FS, et al. Assessment of genetic diversity among 131 safflower (Carthamus tinctorius L.) accessions using peroxidase gene polymorphism (POGP) markers . Mol Biol Rep. 2022; 49: 6531-9.

[15] Singh KN, Rawat S, Kumar K, et al. Identification of significant marker-trait associations for Fusarium wilt resistance in a genetically diverse core collection of safflower using AFLP and SSR markers. J Appl Genet. 2022; 63: 447-62.

[16] Wu Z, Liu H, Zhan W, et al. The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis . Plant Biotechnol J. 2021; 19: 1725-42.

[17] Wang X, Han LQ, Li J, et al. Next-generation bulked segregant analysis for Breeding 4.0. Cell Rep. 2023; 42: 113039.

[18] Azam M, Zhang S, Huai Y, et al. Identification of genes for seed isoflavones based on bulk segregant analysis sequencing in soybean natural population. Theor Appl Genet. 2023; 136: 13.

[19] Zhai RN, Huang AH, Mo RX, et al. SNP-based bulk segregant analysis revealed disease resistance QTLs associated with northern corn leaf blight in maize. Front Genet. 2022; 13: 1038948.

[20] Li D, Ye G, Li J, et al. High light triggers flavonoid and polysaccharide synthesis through DoHY5-dependent signaling in Dendrobium officinale . Plant J. 2023; 115: 1114-33.

[21] Zhang HJ, Li ZX, Zhou S, et al. A fungal NRPS-PKS enzyme catalyses the formation of the flavonoid naringenin. Nat Commun. 2022; 13: 6361.

[22] Yin YC, Zhang XD, Gao ZQ, et al. The research progress of chalcone isomerase (CHI) in plants. Mol Biotechnol. 2018; 61: 32-52.

[23] Liu WX, Feng Y, Yu S, et al. The flavonoid biosynthesis network in plants. Int J Mol Sci. 2021; 22: 12824.

[24] Qi SY, He BX, Wang HT, et al. A muti-substrate flavonol O-glucosyltransferases from safflower. Molecules. 2023; 28: 7613.

[25] Xu XY, Xia M, Han Y, et al. Highly promiscuous flavonoid Di-O-glycosyltransferases from Carthamus tinctorius L. Molecules. 2024; 29: 604.

[26] Ren CX, Xi ZQ, Xian B, et al. Identification and characterization of CtUGT3 as the key player of astragalin biosynthesis in Carthamus tinctorius L. J Agric Food Chem. 2023; 71: 16221-32.

[27] Xi Z, Li Y, Liu S, et al. Functional analysis and molecular characterization of UGT95A2, a specialized glycosyltransferase for flavonoid 3′-O-glycosylation in Carthamus tinctorius L. Plant J. 2025; 122: e70213.

[28] Kaashyap M, Kaur S, Ford R, et al. Comprehensive transcriptomic analysis of two RIL parents with contrasting salt responsiveness identifies polyadenylated and non-polyadenylated flower lncRNAs in chickpea. Plant Biotechnol J. 2022; 20: 1402-16.

[29] Xian B, Zhou Y, Hu Y, et al. Genome-wide screen and multi-omics analysis reveal OGT1 participate in the biosynthesis of safflower flavonoid glycosides. Hortic Res. 2024; 11: uhae261.

[30] Wang MY, Ji QS, Lai B, et al. Structure-function and engineering of plant UDP-glycosyltransferase. Comput Struct Biotechnol J. 2023; 21: 5358-71.

[31] Guo D, Gao Y, Liu F, et al. Integrating molecular characterization and metabolites profile revealed CtCHI1’s significant role in Carthamus tinctorius L. BMC Plant Biol. 2019; 19: 376.