CALMODULIN6 negatively regulates cold tolerance by attenuating ICE1-dependent stress responses in tomato.

Chilling temperatures induce an increase in cytoplasmic calcium (Ca2+) ions to transmit cold signals, but the precise role of Calmodulins (CaMs), a type of Ca2+ sensor, in plant tolerance to cold stress remains elusive. In this study, we characterized a tomato (Solanum lycopersicum) CaM gene, CALMODULIN6 (CaM6), which responds to cold stimulus. Overexpressing CaM6 increased tomato sensitivity to cold stress whereas silencing CaM6 resulted in a cold-insensitive phenotype. We showed that CaM6 interacts with Inducer of CBF expression 1 (ICE1) in a Ca2+-independent process and ICE1 contributes to cold tolerance in tomato plants. By integrating RNA-sequencing (RNA-seq) and chromatin immunoprecipitation-sequencing (ChIP-seq) assays, we revealed that ICE1 directly altered the expression of 76 downstream cold-responsive (COR) genes that potentially confer cold tolerance to tomato plants. Moreover, the physical interaction of CaM6 with ICE1 attenuated ICE1 transcriptional activity during cold stress. These findings reveal that CaM6 attenuates the cold tolerance of tomato plants by suppressing ICE1-dependent COR gene expression. We propose a CaM6/ICE1 module in which ICE1 is epistatic to CaM6 under cold stress. Our study sheds light on the mechanism of plant response to cold stress and reveals CaM6 is involved in the regulation of ICE1.

Genome-wide characterization of cytochrome P450 genes reveals the potential roles in fruit ripening and response to cold stress in tomato.

Plant cytochrome P450 (CYP) superfamily, the largest enzyme metabolism family, has been identified in many species and plays a vital role in plant development and stress response via secondary metabolite biosynthesis. A comprehensive identification and functional investigation of CYPs in tomato plants would contribute to deeper understanding of their biological significance. In this study, 268 tomato CYP genes were identified and found to be unevenly located on 12 chromosomes. Based on the phylogenetic analysis, these 268 SlCYPs were classed into two distinct clades (A-type and non-A-type) and nine clans, including 48 families. Moreover, 67 tandem and 22 WGD (whole genome duplication)/segmental duplication events were detected, of which 12 SlCYP genes experienced both WGD/segmental and tandem duplication events, indicating that tandem duplication plays a major role in the expansion of the SlCYP family. Besides, 48 pairs containing 41 SlCYP and 44 AtCYP genes were orthologous, while 216 orthologous pairs were obtained between tomato and potato. The expression level of all SlCYP genes in tomato tissues at different development stages was analyzed, and most expressed SlCYPs showed a tissue-specific pattern. Meanwhile, 143 differentially expressed SlCYPs were identified under cold stress. Furthermore, the RT-qPCR results indicated that SlCYPs may be involved in fruit ripening and cold tolerance in tomato seedlings. These findings provide valuable insights into the evolutionary relationships and functional characteristics of SlCYPs, which can be utilized for further investigation of fruit metabolic pathways and cold tolerance in tomato.

The miR164a-NAM3 module confers cold tolerance by inducing ethylene production in tomato.

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.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation.

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.

SEC1-C3H39 module fine-tunes cold tolerance by mediating its target mRNA degradation in tomato.

Plants adapt to cold stress at the physiological and biochemical levels, thus enabling them to maintain growth and development. However, the molecular mechanism of fine-tuning cold signals remains largely unknown. We addressed the function of SlSEC1-SlC3H39 module in cold tolerance by using SlSEC1 and SlC3H39 knockout and overexpression tomato lines. A tandem CCCH zinc-finger protein SlC3H39 negatively modulates cold tolerance in tomato. SlC3H39 binds to AU-rich elements in the 3'-untranslated region (UTR) to induce mRNA degradation and regulates gene expression post-transcriptionally. We further validate that SlC3H39 participates in post-transcriptional regulation of a variety of cold-responsive genes. An O-linked N-acetylglucosamine transferase SlSEC1 physically interacts with SlC3H39 proteins and negatively regulates cold tolerance in tomato. Further study shows that SlSEC1 is essential for SlC3H39 protein stability and maintains SlC3H39 function in cold tolerance. Genetic analysis shows that SlC3H39 is epistatic to SlSEC1 in cold tolerance. The findings indicate that SlC3H39 negatively modulates plant cold tolerance through post-transcriptional regulation by binding to cold-responding mRNA 3'-UTR and reducing those transcripts. SlSEC1 promotes the O-GlcNAclation status of SlC3H39 and maintains SlC3H39 function in cold tolerance. Taken together, we propose a SlSEC1-SlC3H39 module, which allows plants to balance defense responses and growth processes.

SlMPK1- and SlMPK2-mediated SlBBX17 phosphorylation positively regulates CBF-dependent cold tolerance in tomato.

B-box (BBX) proteins are an important class of zinc finger transcription factors that play a critical role in plant growth and stress response. However, the mechanisms of how BBX proteins participate in the cold response in tomato remain unclear. Here, using approaches of reverse genetics, biochemical and molecular biology we characterized a BBX transcription factor, SlBBX17, which positively regulates cold tolerance in tomato (Solanum lycopersicum). Overexpressing SlBBX17 enhanced C-repeat binding factor (CBF)-dependent cold tolerance in tomato plants, whereas silencing SlBBX17 increased plant susceptibility to cold stress. Crucially, the positive role of SlBBX17 in CBF-dependent cold tolerance was dependent on ELONGATED HYPOCOTYL5 (HY5). SlBBX17 physically interacted with SlHY5 to directly promote the protein stability of SlHY5 and subsequently increased the transcriptional activity of SlHY5 on SlCBF genes under cold stress. Further experiments showed that cold-activated mitogen-activated protein kinases, SlMPK1 and SlMPK2, also physically interact with and phosphorylate SlBBX17 to enhance the interaction between SlBBX17 and SlHY5, leading to enhanced CBF-dependent cold tolerance. Collectively, the study unveiled a mechanistic framework by which SlMPK1/2-SlBBX17-SlHY5 regulated transcription of SlCBFs to enhance cold tolerance, thereby shedding light on the molecular mechanisms of how plants respond to cold stress via multiple transcription factors.

SlIAA23-SlARF6 module controls arbuscular mycorrhizal symbiosis by regulating strigolactone biosynthesis in tomato.

Auxins are a class of phytohormones with roles involved in the establishment and maintenance of the arbuscular mycorrhizal symbiosis (AMS). Auxin response factors (ARFs) and Auxin/Indole-acetic acids (AUX/IAAs), as two transcription factors of the auxin signaling pathway, coregulate the transcription of auxin response genes. However, the interrelation and regulatory mechanism of ARFs and AUX/IAAs in regulating AMS are still unclear. In this study, we found that the content of auxin in tomato roots increased sharply and revealed the importance of the auxin signaling pathway in the early stage of AMS. Notably, SlARF6 was found to play a negative role in AMF colonization. Silencing SlARF6 significantly increased the expression of AM-marker genes, as well as AMF-induced phosphorus uptake. SlIAA23 could interact with SlARF6 in vivo and in vitro, and promoted the AMS and phosphorus uptake. Interestingly, SlARF6 and SlIAA23 played a contrary role in strigolactone (SL) synthesis and accumulation in AMF-colonized roots of tomato plants. SlARF6 could directly bind to the AuxRE motif of the SlCCD8 promoter and inhibited its transcription, however, this effect was attenuated by SlIAA23 through interaction with SlARF6. Our results suggest that SlIAA23-SlARF6 coregulated tomato-AMS via an SL-dependent pathway, thus affecting phosphorus uptake in tomato plants.

Tomato SlMYB15 transcription factor targeted by sly-miR156e-3p positively regulates ABA-mediated cold tolerance.

Cold is a common abiotic stress that seriously affects plant growth and development. MYB transcription factors are regulatory molecules that play important roles in various biological processes. We have previously demonstrated that SlMYB15 positively regulates cold tolerance in tomato. However, the underlying mechanism of SlMYB15-induced cold tolerance remains largely unexplored. Here, cold-induced SlMYB15 was found to be targeted by Solanum lycopersicum (sly)-miR156e-3p, which was decreased by cold stimulus in tomato. Tomato plants overexpressing sly-MIR156e-3p displayed significant enhancement in susceptibility to cold stress, while silencing of sly-miR156e-3p by an artificial microRNA interference strategy caused tomato plants to be more tolerant to cold. Moreover, both overexpression of SlMYB15 and silencing of sly-miR156e-3p increased the accumulation of ABA. SlMYB15 directly binds to the promoter regions of ABA biosynthesis and signalling genes, SlNCED1 and SlABF4, resulting in enhanced cold tolerance. Further experiments showed that SlMYB15 and sly-miR156e-3p also coordinated the cold tolerance of tomato via the reactive oxygen species (ROS) signalling pathway, as reflected by the increased expression of SlRBOH1, enhanced H2O2 and O2•-accumulation, and amplified activity of antioxidant enzymes in SlMYB15-overexpressing and sly-miR156e-3p-silenced plants. Taken together, our results demonstrate that SlMYB15 targeted by sly-miR156e-3p confers higher survivability to cold stress via ABA and ROS signals. This study provides valuable information for breeding improved crop cultivars better equipped with cold tolerance.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.).

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.

An ultra-high-density genetic map provides insights into genome synteny, recombination landscape and taproot skin colour in radish (Raphanus sativus L.).

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.

Melatonin confers cadmium tolerance by modulating critical heavy metal chelators and transporters in radish plants.

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.

Genome-wide characterization and evolutionary analysis of heat shock transcription factors (HSFs) to reveal their potential role under abiotic stresses in radish (Raphanus sativus L.).

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.

Transcriptome-based gene expression profiling of diploid radish (Raphanus sativus L.) and the corresponding autotetraploid.

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.

Genome- and Transcriptome-Wide Characterization of bZIP Gene Family Identifies Potential Members Involved in Abiotic Stress Response and Anthocyanin Biosynthesis in Radish (Raphanus sativus L.).

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.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.).

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.

Dissecting Root Proteome Changes Reveals New Insight into Cadmium Stress Response in Radish (Raphanus sativus L.).

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.

Identification of critical genes associated with lignin biosynthesis in radish (Raphanus sativus L.) by de novo transcriptome sequencing.

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.

Genome-wide characterization of the WRKY gene family in radish (Raphanus sativus L.) reveals its critical functions under different abiotic stresses.

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.

Genome-wide identification of long non-coding RNAs and their potential functions in radish response to salt stress.

Long non-coding RNAs (lncRNAs) are increasingly recognized as cis- and trans-acting regulators of protein-coding genes in plants, particularly in response to abiotic stressors. Among these stressors, high soil salinity poses a significant challenge to crop productivity. Radish (Raphanus sativus L.) is a prominent root vegetable crop that exhibits moderate susceptibility to salt stress, particularly during the seedling stage. Nevertheless, the precise regulatory mechanisms through which lncRNAs contribute to salt response in radish remain largely unexplored. In this study, we performed genome-wide identification of lncRNAs using strand-specific RNA sequencing on radish fleshy root samples subjected to varying time points of salinity treatment. A total of 7,709 novel lncRNAs were identified, with 363 of them displaying significant differential expression in response to salt application. Furthermore, through target gene prediction, 5,006 cis- and 5,983 trans-target genes were obtained for the differentially expressed lncRNAs. The predicted target genes of these salt-responsive lncRNAs exhibited strong associations with various plant defense mechanisms, including signal perception and transduction, transcription regulation, ion homeostasis, osmoregulation, reactive oxygen species scavenging, photosynthesis, phytohormone regulation, and kinase activity. Notably, this study represents the first comprehensive genome-wide analysis of salt-responsive lncRNAs in radish, to the best of our knowledge. These findings provide a basis for future functional analysis of lncRNAs implicated in the defense response of radish against high salinity, which will aid in further understanding the regulatory mechanisms underlying radish response to salt stress.

BRASSINAZOLE RESISTANT 1 Mediates Brassinosteroid-Induced Calvin Cycle to Promote Photosynthesis in Tomato.

Calvin cycle is a sequence of enzymatic reactions that assimilate atmospheric CO2 in photosynthesis. Multiple components are known to participate in the induction or suppression of the Calvin cycle but the mechanism of its regulation by phytohormones is still unclear. Brassinosteroids (BRs) are steroid phytohormones that promote photosynthesis and crop yields. In this study, we study the role of BRs in regulating Calvin cycle genes to further understand the regulation of the Calvin cycle by phytohormones in tomatoes. BRs and their signal effector BRASSINAZOLE RESISTANT 1 (BZR1) can enhance the Calvin cycle activity and improve the photosynthetic ability. BRs increased the accumulation of dephosphorylated form of BZR1 by 94% and induced an 88-126% increase in the transcription of key genes in Calvin cycle FBA1, RCA1, FBP5, and PGK1. BZR1 activated the transcription of these Calvin cycle genes by directly binding to their promoters. Moreover, silencing these Calvin cycle genes impaired 24-epibrassinolide (EBR)-induced enhancement of photosynthetic rate, the quantum efficiency of PSII, and V c,max and J max . Taken together, these results strongly suggest that BRs regulate the Calvin cycle in a BZR1-dependent manner in tomatoes. BRs that mediate coordinated regulation of photosynthetic genes are potential targets for increasing crop yields.