Combined widely targeted metabolomics and transcriptomics analysis reveals differentially accumulated metabolites and the underlying molecular bases in fleshy taproots of distinct radish genotypes

Radish is an important taproot crop with medicinal and edible uses that is cultivated worldwide. However, the differences in metabolites and the underlying molecular bases among different radish types remain largely unknown. In the present study, a combined analysis of liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) and RNA-Seq data was conducted to uncover important differentially accumulated metabolites (DAMs) among radish accessions with green, white and red taproot flesh colours. A total of 657 metabolites were identified and 138 DAMs were commonly present in red vs. green and red vs. white accessions. Red accessions were rich in phenolic compounds, while green and white accessions had more amino acids. Additionally, 41 metabolites and 98 genes encoding 37 enzymes were enriched in the shikimate and phenolic biosynthesis pathways. CHS is the rate-limiting enzyme determining flavonoid differences among accessions. A total of 119 candidate genes might contribute to red accession-specific accumulated metabolites. Specifically, one gene cluster consisting of 16 genes, including one RsMYB1, RsMYBL2, RsTT8, RsDFR, RsANS, Rs4CL3, RsSCPL10, Rs3AT1 and RsSAP2 gene, two RsTT19 and RsWRKY44 genes and three RsUGT genes, might be involved in anthocyanin accumulation in red radish fleshy taproots. More importantly, an InDel marker was developed based on an RsMYB1 promoter mutation, and the accuracy reached 95.9% when it was used to select red-fleshed radishes. This study provides comprehensive insights into the metabolite differences and underlying molecular mechanisms in fleshy taproots among different radish genotypes and will be beneficial for the genetic improvement of radish nutritional quality.