Search

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1 In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish., (© 2023 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Auxin/indole-3-acetic acid (Aux/IAA) genes encode short lived nuclear proteins that cooperated with auxin or auxin response factor (ARF), which are involved in plant growth and developmental processes. However, it's still ambiguous how the Aux/IAA genes regulate the process governing taproot thickening in radish. Herein, 65 Aux/IAA genes were identified from the radish genome. Gene duplication analysis showed that two pairs of tandem duplication and 17 (27%) segmental duplication events were identified among Aux/IAA family genes in radish. Transcriptomic analysis revealed that most of Aux/IAA genes (52/65) exhibited differential expression pattern in different root tissues, and six root-specific genes were highly expressed in root cortex, cambium, xylem, and root tip in radish. RT-qPCR analysis showed that the expression level of RsIAA33 was the highest at cortex splitting stage (CSS), and early expanding stage (ES). Furthermore, amiRNA-mediated gene silencing of RsIAA33 indicated that it could inhibit the reproductive growth, thus promoting taproot thickening and development. These results would provide valuable information for elucidating the molecular function of Aux/IAA genes involved in taproot thickening in radish., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Background: Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size., Results: A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis., Conclusion: This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

Cadmium (Cd) is an environmental pollutant that causes health hazard to living organisms. Melatonin (MT) has emerged as a ubiquitous pleiotropic molecule capable of coordinating heavy metal (HM) stresses in plants. However, it remains unclear how melatonin mediates Cd homeostasis and detoxification at transcriptional and/or post-transcriptional levels in radish. Herein, the activities of five key antioxidant enzymes were increased, while root and shoot Cd contents were dramatically decreased by melatonin. A combined small RNA and transcriptome sequencing analysis showed that 14 differentially expressed microRNAs (DEMs) and 966 differentially expressed genes (DEGs) were shared between the Cd and Cd + MT conditions. In all, 23 and ten correlated miRNA-DEG pairs were identified in Con vs. Cd and Con vs. Cd + MT comparisons, respectively. Several DEGs encoding yellow stripe 1-like (YSL), heavy metal ATPases (HMA), and ATP-binding cassette (ABC) transporters were involved in Cd transportation and sequestration in radish. Root exposure to Cd 2+ induced several specific signaling molecules, which consequently trigger some HM chelators, transporters, and antioxidants to achieve reactive oxygen species (ROS) scavenging and detoxification and eliminate Cd toxicity in radish plants. Notably, transgenic analysis revealed that overexpression of the RsMT1 (Metallothionein 1) gene could enhance Cd tolerance of tobacco plants, indicating that the exogenous melatonin confers Cd tolerance, which might be attributable to melatonin-mediated upregulation of RsMT1 gene in radish plants. These results could contribute to dissecting the molecular basis governing melatonin-mediated Cd stress response in plants and pave the way for high-efficient genetically engineering low-Cd-content cultivars in radish breeding programs., (© 2020 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

High-density genetic map is a valuable tool for exploring novel genomic information, quantitative trait locus (QTL) mapping and gene discovery of economically agronomic traits in plant species. However, high-resolution genetic map applied to tag QTLs associated with important traits and to investigate genomic features underlying recombination landscape in radish (Raphanus sativus) remains largely unexplored. In this study, an ultra-high-density genetic map with 378 738 SNPs covering 1306.8 cM in nine radish linkage groups (LGs) was developed by a whole-genome sequencing-based approach. A total of 18 QTLs for 11 horticulture traits were detected. The map-based cloning data indicated that the R2R3-MYB transcription factor RsMYB90 was a crucial candidate gene determining the taproot skin colour. Comparative genomics analysis among radish, Brassica rapa and B. oleracea genome revealed several genomic rearrangements existed in the radish genome. The highly uneven distribution of recombination was observed across the nine radish chromosomes. Totally, 504 recombination hot regions (RHRs) were enriched near gene promoters and terminators. The recombination rate in RHRs was positively correlated with the density of SNPs and gene, and GC content, respectively. Functional annotation indicated that genes within RHRs were mainly involved in metabolic process and binding. Three QTLs for three traits were found in the RHRs. The results provide novel insights into the radish genome evolution and recombination landscape, and facilitate the development of effective strategies for molecular breeding by targeting and dissecting important traits in radish., (© 2019 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Basic leucine zipper (bZIP) transcription factors play crucial roles in various abiotic stress responses as well as anthocyanin accumulation. Anthocyanins are most abundant in colorful skin radish, which exhibit strong antioxidant activity that offers benefits for human health. Here, a total of 135 bZIP-encoding genes were identified from radish genome. Synteny analysis showed that 104 radish and 63 Arabidopsis bZIP genes were orthologous. Transcriptome analysis revealed that 10 RsbZIP genes exhibited high-expression levels in radish taproot (RPKM>10). Specifically, RsbZIP010 exhibited down-regulated expression under Cd, Cr and Pb stresses, whereas RsbZIP031 and RsbZIP059 showed significant down-regulation under heat and salt stresses, respectively. RT-qPCR analysis indicated that RsbZIP011 and RsbZIP102 were significantly up-regulated in the tissues of radish with high anthocyanin contents. Furthermore, the promoter sequences of 39 anthocyanin-related genes were found to contain G-box or ACE-box elements that could be recognized by bZIP family members. Taken together, several RsbZIP s might be served as critical regulators in radish taproot under Cd, Cr, Pb, heat and salt stresses. RsbZIP011 and RsbZIP102 were the potential participants in anthocyanin biosynthesis pathway of radish. These results facilitate further investigation on functional characterization of bZIP genes in response to abiotic stress and anthocyanin biosynthesis in radish.

Main conclusion Among 247 RsAP2/ERF identified, the majority of the 21 representatives were preferably expressed under drought and heat while suppressed under heavy metals, indicating their potential roles in abiotic stress responses and tolerance. APETALA2/Ethylene-Responsive factor (AP2/ERF) transcription factor (TF) is one of the largest gene families in plants that play a fundamental role in growth and development as well as biotic and/or abiotic stresses responses. Although AP2/ERFs have been extensively characterized in many plant species, little is known about this family in radish, which is an important root vegetable with various medicinal properties. The available genome provides valuable opportunity to identify and characterize the global information on AP2/ERF TFs in radish. In this study, a total of 247 ERF family genes were identified from the radish genome, and sequence alignment and phylogenetic analyses classified the AP2/ERF superfamily into five groups (AP2, ERF, DREB, RAV and soloist). Motif analysis showed that other than AP2/ERF domains, other conserved regions were selectively distributed among different clades in the phylogenetic tree. Chromosome location analysis showed that tandem duplication may result in the expansion of RsAP2/ERF gene family. The RT-qPCR analysis confirmed that a proportion of AP2/ERF genes were preferably expressed under drought and heat stresses, whereas they were suppressed under the ABA and heavy metal stresses. These results provided valuable information for further evolutionary and functional characterization of RsAP2/ERF genes, and contributed to genetic improvement of stress tolerances in radish and other root vegetable crops., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: Abiotic stresses due to climate change pose a great threat to crop production. Heat shock transcription factors (HSFs) are vital regulators that play key roles in protecting plants against various abiotic stresses. Therefore, the identification and characterization of HSFs is imperative to dissect the mechanism responsible for plant stress responses. Although the HSF gene family has been extensively studied in several plant species, its characterization, evolutionary history and expression patterns in the radish (Raphanus sativus L.) remain limited., Results: In this study, 33 RsHSF genes were obtained from the radish genome, which were classified into three main groups based on HSF protein domain structure. Chromosomal localization analysis revealed that 28 of 33 RsHSF genes were located on nine chromosomes, and 10 duplicated RsHSF genes were grouped into eight gene pairs by whole genome duplication (WGD). Moreover, there were 23 or 9 pairs of orthologous HSFs were identified between radish and Arabidopsis or rice, respectively. Comparative analysis revealed a close relationship among radish, Chinese cabbage and Arabidopsis. RNA-seq data showed that eight RsHSF genes including RsHSF-03, were highly expressed in the leaf, root, cortex, cambium and xylem, indicating that these genes might be involved in plant growth and development. Further, quantitative real-time polymerase chain reaction (RT-qPCR) indicated that the expression patterns of 12 RsHSF genes varied upon exposure to different abiotic stresses including heat, salt, and heavy metals. These results indicated that the RsHSFs may be involved in abiotic stress response., Conclusions: These results could provide fundamental insights into the characteristics and evolution of the HSF family and facilitate further dissection of the molecular mechanism responsible for radish abiotic stress responses.

Polyploidy is an important evolutionary factor in most land plant lineages which possess more than two complete sets of chromosomes. Radish (Raphanus sativus L.) is an economically annual/biennial root vegetable crop worldwide. However, the expression patterns of duplicated homologs involved in the autopolyploidization remains unclear. In present study, the autotetraploid radish plants (2n = 4x = 36) were produced with colchicine and exhibited an increase in the size of flowers, leaves, stomata and pollen grains. The differential gene expression (DGE) profiling was performed to investigate the differences in gene expression patterns between diploid and its corresponding autotetraploid by RNA-Sequencing (RNA-Seq). Totally, 483 up-regulated differentially expressed genes (DEGs) and 408 down-regulated DEGs were detected in diploid and autotetraploid radishes, which majorly involved in the pathways of hormones, photosynthesis and stress response. Moreover, the xyloglucan endotransglucosylase/hydrolase (XTH) and pectin methylesterases (PME) family members related to cell enlargement and cell wall construction were found to be enriched in GO enrichment analysis, of which XTH family members enriched in "apoplast" and "cell wall" terms, while PME family members enriched in "cell wall" term. Reverse-transcription quantitative PCR (RT-qPCR) analysis indicated that the expression profile of DEGs were consistent with results from the RNA-Seq analysis. The DEGs involved in cell wall construction and auxin metabolism were predicted to be associated with organs size increase of autotetraploid radishes in the present study. These results could provide valuable information for elucidating the molecular mechanism underlying polyploidization and facilitating further genetic improvements of important traits in radish breeding programs.

Key Message: The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. The WRKY is among the largest transcription factors (TFs) associated with multiple biological activities for plant survival, including control response mechanisms against abiotic stresses such as heat, salinity, and heavy metals. Radish is an important root vegetable crop and therefore characterization and expression pattern investigation of WRKY transcription factors in radish is imperative. In the present study, 126 putative WRKY genes were retrieved from radish genome database. Protein sequence and annotation scrutiny confirmed that RsWRKY proteins possessed highly conserved domains and zinc finger motif. Based on phylogenetic analysis results, RsWRKYs candidate genes were divided into three groups (Group I, II and III) with the number 31, 74, and 20, respectively. Additionally, gene structure analysis revealed that intron-exon patterns of the WRKY genes are highly conserved in radish. Linkage map analysis indicated that RsWRKY genes were distributed with varying densities over nine linkage groups. Further, RT-qPCR analysis illustrated the significant variation of 36 RsWRKY genes under one or more abiotic stress treatments, implicating that they might be stress-responsive genes. In total, 126 WRKY TFs were identified from the R. sativus genome wherein, 35 of them showed abiotic stress-induced expression patterns. These results provide a genome-wide characterization of RsWRKY TFs and baseline for further functional dissection and molecular evolution investigation, specifically for improving abiotic stress resistances with an ultimate goal of increasing yield and quality of radish.

Cadmium (Cd) is a widespread heavy metal of particular concern with respect to the environment and human health. Although intensive studies have been conducted on Cd-exposed transcriptome profiling, little systematic proteome information is available on the molecular mechanism of Cd stress response in radish. In this study, the radish root proteome under Cd stress was investigated using a quantitative multiplexed proteomics approach. Seedlings were grown in nutrient solution without Cd (control) or with 10 or 50 μM CdCl2 for 12 h (Cd10 and Cd50, respectively). In total, 91 up- and 66 down-regulated proteins were identified in the control vs Cd10 comparison, while 340 up- and 286 down-regulated proteins were identified in the control vs Cd50 comparison. Functional annotation indicated that these differentially expressed proteins (DEPs) were mainly involved in carbohydrate and energy metabolism, stress and defense and signal transduction processes. Correlation analysis showed that 33 DEPs matched with their transcripts, indicating a relatively low correlation between transcript and protein levels under Cd stress. Quantitative real-time PCR evidenced the expression patterns of 12 genes encoding their corresponding DEPs. In particular, several pivotal proteins associated with carbohydrate metabolism, ROS scavenging, cell transport and signal transduction were involved in the coordinated regulatory network of the Cd stress response in radish. Root exposure to Cd2+ activated several key signaling molecules and metal-containing transcription factors, and subsequently some Cd-responsive functional genes were mediated to reduce Cd toxicity and re-establish redox homeostasis in radish. This is a first report on comprehensive proteomic characterization of Cd-exposed root proteomes in radish. These findings could facilitate unraveling of the molecular mechanism underlying the Cd stress response in radish and provide fundamental insights into the development of genetically engineered low-Cd-content radish cultivars., (© The Author 2017. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Radish is an important root vegetable crop with high nutritional, economic, and medicinal value. Lignin is an important secondary metabolite possessing a great effect on plant growth and product quality. To date, lignin biosynthesis-related genes have been identified in some important plant species. However, little information on characterization of critical genes involved in plant lignin biosynthesis is available in radish. In this study, a total of 71,148 transcripts sequences were obtained from radish root, of which 66 assembled unigenes and ten candidate genes were identified to be involved in lignin monolignol biosynthesis. Full-length cDNA sequences of seven randomly selected genes were isolated and sequenced from radish root, and the assembled unigenes covered more than 80% of their corresponding cDNA sequences. Moreover, the lignin content gradually accumulated in leaf during the developmental stages, and it increased from pre-cortex to cortex splitting stage, followed by a decrease at thickening stage and then increased at mature stage in root. RT-qPCR analysis revealed that all these genes except RsF5H exhibited relatively low expression level in root at thickening stage. The expression profiles of Rs4CL5, RsCCoAOMT1, and RsCOMT genes were consistent with the changes of root lignin content, implying that these candidate genes may play important roles in lignin formation in radish root. These findings would provide valuable information for identification of lignin biosynthesis-related genes and facilitate dissection of molecular mechanism underlying lignin biosynthesis in radish and other root vegetable crops.

Background: Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish., Results: In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions., Conclusions: These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.