Raphanus has undergone a lengthy evolutionary process and has rich diversity. However, the inter- and intraspecific phylogenetic relationships and genetic diversity of this genus are not well understood. Through SSR-sequencing and multi-analysis of 939 wild, semi-wild and cultivated accessions, we discovered that the European wild radish (EWR) population is separated from cultivated radishes and has a higher genetic diversity. Frequent intraspecific genetic exchanges occurred in the whole cultivated radish (WCR) population; there was considerable genetic differentiation within the European cultivated radish (ECR) population, which could drive radish diversity formation. Among the ECR subpopulations, European primitive cultivated radishes (EPCRs) with higher genetic diversity are most closely related to the EWR population and exhibit a gene flow with rat-tail radishes (RTRs) and black radishes (BRs)/oil radishes (ORs). Among Asian cultivated radishes (ACRs), Chinese big radishes (CBRs) with a relatively high diversity are furthest from the EWR population, and most Japanese/Korean big radishes (JKBRs) are close to CBR accessions, except for a few old Japanese landraces that are closer to the EPCR. The CBR and JKBR accessions are independent of RTR accessions; however, phylogenetic analysis indicates that the RTR is sister to the clade of CBR (including JWR), which suggests that the RTR may share the most recent common ancestry with CBRs and JWRs. In addition, Japanese wild radishes (JWRs), (namely, R. sativus forma raphanistroides) are mainly scattered between CBRs and EPCRs in PCoA analysis. Moreover, JWRs have a strong gene exchange with the JKBR, OR and RTR subpopulations. American wild radishes (AWRs) are closely related to European wild and cultivated radishes, and have a gene flow with European small radishes (ESRs), suggesting that the AWR developed from natural hybridization between the EWR and the ESR. Overall, this demonstrates that Europe was the origin center of the radish, and that Europe, South Asia and East Asia appear to have been three independent domestication centers. The EPCR, AWR and JWR, as semi-wild populations, might have played indispensable transitional roles in radish evolution. Our study provides new perspectives into the origin, evolution and genetic diversity of Raphanus and facilitates the conservation and exploitation of radish germplasm resources.

Radish (Raphanus sativus L.) is rich in specific glucosinolates (GSLs), which benefit human health and special flavor formation. Although the basic GSLs metabolic pathway in Brassicaceae plants is clear, the regulating mechanism for specific glucosinolates content in radish fleshy taproots is not well understood. In this study, we discovered that there was a significant difference in the GSLs profiles and the content of various GSLs components. Glucoraphasatin (GRH) is the most predominant GSL in radish taproots of different genotypes as assessed by HPLC analysis. Further, we compared the taproot transcriptomes of three radish genotypes with high and low GSLs content by employing RNA-seq. Totally, we identified forty-one differentially expressed genes related to GSLs metabolism. Among them, thirteen genes (RsBCAT4, RsIPMDH1, RsMAM1a, RsMAM1b, RsCYP79F1, RsGSTF9, RsGGP1, RsSUR1, RsUGT74C1, RsST5b, RsAPK1, RsGSL-OH, and RsMYB28) were significantly higher co-expressed in the high content genotypes than in low content genotype. Notably, correlation analysis indicated that the expression level of RsMYB28, as an R2R3 transcription factor directly regulating aliphatic glucosinolate biosynthesis, was positively correlated with the GRH content. Co-expression network showed that RsMYB28 probably positively regulated the expression of the above genes, particularly RsSUR1, and consequently the synthesis of GRH. Moreover, the molecular mechanism of the accumulation of this 4-carbon (4C) GSL in radish taproots was explored. This study provides new perspectives on the GSLs accumulation mechanism and genetic improvements in radish taproots.