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 uptake and accumulation at transcriptional and/or posttranscriptional level in radish. Herein, the activities of five key antioxidant enzymes were increased, while root and shoot Cd content 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 Cd and Cd+MT conditions. In all, 18 and 8 correlated miRNA-DEG pairs were identified in Con vs. Cd and Con vs. Cd+MT comparison, respectively. Several DEGs encoding YSL (yellow stripe 1-like), HMA (heavy metal ATPases) and ABC (ATP-binding cassette) transporters were involved in Cd chelation and sequestration in radish. Root exposure to Cd 2+ induced excess reactive oxygen species (ROS) and some specific signaling molecules including CaMs, CDPKs and MAPKs, which consequently trigger several metal-binding peptides and metal transporter genes to establish redox homeostasis and eliminate Cd toxicity in radish plants. Notably, transgenic analysis revealed that over-expression of RsMT1 gene could enhance Cd tolerance of tobacco plants, indicating that the exogenous melatonin confers Cd tolerance might attribute to melatonin-mediated up-regulation of RsMT1 gene in radish plants. These results could contribute to clarify molecular basis of melatonin-mediated Cd stress response in radish, and pave the way for high-efficient genetically engineering low-Cd-content cultivars in radish breeding programs.

The fleshy taproot of radish is an important storage organ determining its yield and quality. Taproot thickening is a complex developmental process in radish. However, the molecular mechanisms governing this process remain unclear at the proteome level. In this study, a comparative proteomic analysis was performed to analyze the proteome changes at three developmental stages of taproot thickening using iTRAQ approach. In total, 1862 differentially expressed proteins (DEPs) were identified from 6342 high-confidence proteins, among which 256 up-regulated proteins displayed overlapped accumulation in S1 (pre-cortex splitting stage) vs. S2 (cortex splitting stage) and S1 vs. S3 (expanding stage) pairs, whereas 122 up-regulated proteins displayed overlapped accumulation in S1 vs. S3 and S2 vs. S3 pairs. Gene Ontology (GO) and pathway enrichment analysis showed that these DEPs were mainly involved in several processes such as “starch and sucrose metabolism”, “plant hormone signal transduction”, and “biosynthesis of secondary metabolites”. A high concordance existed between iTRAQ and RT-qPCR at the mRNA expression levels. Furthermore, association analysis showed that 187, 181, and 96 DEPs were matched with their corresponding differentially expressed genes (DEGs) in S1 vs. S2, S1 vs. S3, and S2 vs. S3 comparison, respectively. Notably, several functional proteins including cell division cycle 5-like protein (CDC5), expansin B1 (EXPB1), and xyloglucan endotransglucosylase/hydrolase protein 24 (XTH24) were responsible for cell division and expansion during radish taproot thickening process. These results could facilitate a better understanding of the molecular mechanism underlying taproot thickening, and provide valuable information for the identification of critical genes/proteins responsible for taproot thickening in root vegetable crops.

SummaryHigh‐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.

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 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.

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.

SUMMARY Because of a high sensitivity to cold, both the yield and quality of tomato (Solanum lycopersicum L.) are severely restricted by cold stress. The NAC transcription factor (TF) family has been characterized as an important player in plant growth, development, and the stress response, but the role of NAC TFs in cold stress and their interaction with other post‐transcriptional regulators such as microRNAs in cold tolerance remains elusive. Here, we demonstrated that SlNAM3, the predicted target of Sl‐miR164a/b‐5p, improved cold tolerance as indicated by a higher maximum quantum efficiency of photosystem II (Fv/Fm), lower relative electrolyte leakage, and less wilting in SlNAM3‐overexpression plants compared to wild‐type. Further genetic and molecular confirmation revealed that Sl‐miR164a/b‐5p functioned upstream of SlNAM3 by inhibiting the expression of the latter, thus playing a negative role in cold tolerance. Interestingly, this role is partially mediated by an ethylene‐dependent pathway because either Sl‐miR164a/b‐5p silencing or SlNAM3 overexpression improved cold tolerance in the transgenic lines by promoting ethylene production. Moreover, silencing of the ethylene synthesis genes, SlACS1A, SlACS1B, SlACO1, and SlACO4, resulted in a significant decrease in cold tolerance. Further experiments demonstrated that NAM3 activates SlACS1A, SlACS1B, SlACO1, and SlACO4 transcription by directly binding to their promoters. Taken together, the present study identified the miR164a‐NAM3 module conferring cold tolerance in tomato plants via the direct regulation of SlACS1A, SlACS1B, SlACO1, and SlACO4 expression to induce ethylene synthesis.

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.

No abstract