Search

Summary: 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, RsVRN1In‐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. [ABSTRACT FROM AUTHOR]

Radish (Raphanus sativus L.) is an important root vegetable crop that is easily infected by various pathogens that result in decreased yield and quality. Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes play vital roles in resisting pathogen infection in plants. However, the genome-wide characterization and functional roles of NBS-LRR genes remain largely unexplored in radish. Here, a total of 187 RsNBS-LRR genes were identified at the whole-genome level in radish, among which 80 RsNBS-LRR genes were unevenly distributed on nine radish chromosomes. Interestingly, 15 clusters containing 36 RsNBS-LRR genes occurred in eight chromosomes. RNA-Seq data showed that several RsNBS-LRR genes exhibited significant differential expression profiles in different radish tissues. Moreover, a range of cis-acting regulatory elements associated with ABA, MeJA, or SA were identified in the promoter region of some RsNBS-LRR genes. RT-qPCR analysis showed that the expression of a few RsNBS-LRR genes (e.g., RsNBS021 and RsNBS163) was significantly induced under Peronospora parasitica infection and/or ABA treatment, indicating that they might play critical roles in ABA-dependent defense resistance processes. These results could enhance our understanding of the evolutionary relationship of RsNBS-LRR genes and facilitate the genetic manipulation of disease resistance in radish breeding programs. [ABSTRACT FROM AUTHOR]

Melatonin (MT) is a tryptophan-derived natural product that plays a vital role in plant response to abiotic stresses, including heavy metals (HMs). However, it remains elusive how exogenous MT mediates lead (Pb) accumulation and detoxification at the methylation and transcriptional levels in radish. In this study, decreased Pb accumulation and increased antioxidant enzyme activity were detected under MT treatment in radish. Single-base resolution maps of DNA methylation under Pb stress (Pb200) and Pb plus MT treatment (Pb_50MT) were first generated. The genome-wide methylation level was increased under Pb stress, while an overall loss of DNA methylation was observed under MT treatment. The differentially methylated region (DMR)-associated genes between Pb_50MT and Pb200 were uniquely enriched in ion binding terms, including cation binding, iron ion binding, and transition metal ion binding. Hyper-DMRs between Pb200 and Control exhibited a decreasing trend of methylation under Pb_50MT treatment. A few critical upregulated antioxidant genes (e.g., RsAPX2, RsPOD52 and RsGST) exhibited decreased methylation levels under MT treatment, which enabled the radish plants to scavenge lead-induced reactive oxygen species (ROS) and decrease oxidative stress. Notably, several MT-induced HM transporter genes with low methylation (e.g., RsABCF5, RsYSL7 and RsHMT) and transcription factors (e.g., RsWRKY41 and RsMYB2) were involved in reducing Pb accumulation in radish roots. These findings could facilitate comprehensive elucidation of the molecular mechanism underlying MT-mediated Pb accumulation and detoxification in radish and other root vegetable crops. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 Cd2+ 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. [ABSTRACT FROM AUTHOR]

Summary: 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. [ABSTRACT FROM AUTHOR]

• A total of 28 RsZIP genes were identified and classified into 18 subgroups. • Four tandem and six WGD/segmental duplication events exist in radish genome. • Four RsZIP genes were induced at 2 h under Pb and Cd treatment. • RsIRT1d and RsIRT1e confer yeast Pb but no Cd tolerance. Radish is an important root vegetable crop susceptible to heavy metal (HM) stresses. Zinc- and iron-regulated transporter-like proteins play vital roles in transporting metal ion to enhance HM tolerance or decrease detoxification in plants. In this study, 27 out of 28 RsZIP genes were unevenly located on eight chromosomes with exception of chromosome 3 (Chr.3). The evolutionary pattern divided the whole RsZIP proteins into 18 subgroups. Furthermore, four tandem and six WGD/segmental duplication events were identified, while 24 pairs containing 16 RsZIP and 13 AtZIP genes were orthologous. Moreover, the expression levels of 17 RsZIP genes including RsZIP10b,c and RsIRT1a,b were significantly decreased under 200 mg/L CdCl 2 and 1000 mg/L Pb(NO 3) 2 treatment. RT-qPCR analysis indicated that four RsZIP genes (RsIRT1d, RsIRT1e, RsIAR1 and RsZIP10d) were rapidly induced at 2 h under Pb and Cd treatment, while both RsZIP1 and RsZIP10b were downregulated. The results of yeast functional complementation indicated that RsIRT1d and RsIRT1e confer yeast Pb but no Cd tolerance, suggesting that they might play an important role in the transport of metal ions. These results could provide valuable insights into the characteristics of ZIP family genes and facilitate further exploring the molecular mechanism underlying HM stress responses in radish and other root vegetable crops. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Key message : The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. Abstract: 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Radish (Raphanus sativus L.) taproot contains high concentrations of flavonoids, including anthocyanins (ATCs), in red-skinned genotypes. However, little information on the genetic regulation of ATC biosynthesis in radish is available. A genome-wide association study of radish red skin color was conducted using whole-genome sequencing data derived from 179 radish genotypes. The R2R3-MYB transcription factor production of anthocyanin pigment 2 (PAP2) gene was found in the region associated with a leading SNP located on chromosome 2. The amino acid sequence encoded by the RsPAP2 gene was different from those of the other published RsMYB genes responsible for the red skin color of radish. The overexpression of the RsPAP2 gene resulted in ATC accumulation in Arabidopsis and radish, which was accompanied by the upregulation of several ATC-related structural genes. RsPAP2 was found to bind the RsUFGT and RsTT8 promoters, as shown by a dual-luciferase reporter system and a yeast one-hybrid assay. The promoter activities of the RsANS, RsCHI, RsPAL, and RsUFGT genes could be strongly activated by coinfiltration with RsPAP2 and RsTT8. These findings showed the effectiveness of GWAS in identifying candidate genes in radish and demonstrated that RsPAP2 could (either directly or together with its cofactor RsTT8) regulate the transcript levels of ATC-related genes to promote ATC biosynthesis, facilitating the genetic enhancement of ATC contents and other related traits in radish. [ABSTRACT FROM AUTHOR]

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 ArabidopsisbZIP 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 RsbZIPs 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. [ABSTRACT FROM AUTHOR]

• The identified 65 Aux/IAA genes show an uneven chromosomal distribution in radish. • Six Aux/IAA genes exhibit specifically high expression level in radish root tissues. • RsIAA33 gene might involve in promoting taproot thickening and development in radish. 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. [ABSTRACT FROM AUTHOR]

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. • Totally 247 AP2-ERF genes were firstly genome-widely identified in radish. • The AP2/ERF superfamily could be classified into five groups. • Tandem duplications play critical role in the expansion of RsAP2/ERF genes. • A proportion of AP2/ERF genes were differentially expressed under abiotic stresses. [ABSTRACT FROM AUTHOR]