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.

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.

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.

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.

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.

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.

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.

Melatonin-induced DNA demethylation of metal transporters and antioxidant genes alleviates lead stress in radish plants.

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.

Genome-wide identification of AUX/IAA in radish and functional characterization of RsIAA33 gene during taproot thickening.

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.

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.

A genome-wide association study uncovers a critical role of the RsPAP2 gene in red-skinned Raphanus sativus L.

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.

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.

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.

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.

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.

Genome-wide characterization of the AP2/ERF gene family in radish (Raphanus sativus L.): Unveiling evolution and patterns in response to abiotic stresses.

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.

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.

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.