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.

Fine mapping and candidate gene analysis of CRA3.7 conferring clubroot resistance in Brassica rapa.

KEY MESSAGE: In this study, we fine-mapped a clubroot resistance gene CRA3.7 in Chinese cabbage and developed its closely linked marker syau-InDel3008 for marker-assisted selection in CR cultivars breeding. Chinese cabbage is an important leafy vegetable rich in many nutrients widely grown in China. Clubroot disease caused by an obligate biotrophic pathogen Plasmodiophora brassicae was rapidly spread and challenged to Chinese cabbage production. A clubroot resistance (CR) gene, CRA3.7, was mapped on chromosome A03 of Brassica rapa. A Chinese cabbage line 'CR510', which harbor homozygous resistance locus CRA3.7 was selected from a BC4F3 family. 'CR510' was crossed with a clubroot susceptible Chinese cabbage inbred line '59-1'. Total 51 recombinant plants were identified from an F2 population including 3000 individuals. These recombinants were selfed and the clubroot resistance of F2/3 families was evaluated. Finally, a clubroot resistance gene CRA3.7 was fine-mapped to an interval of approximately 386 kb between marker syau-InDel3024 and syau-InDel3008. According to the reference genome, total 54 genes including five encoding the TIR-NBS-LRR proteins was annotated in the fine-mapped region. Further, nine candidate's gene expression in parental lines at 7, 14 and 21 days after inoculation of P. brassicae were evaluated. Bra019376, Bra019401, Bra019403 and Bra019410 are highly expressed in 'CR510' than '59-1'. Gene sequence of Bra019410 from 'CR510' was cloned and identified different from CRa. Therefore, Bra019376, Bra019401, Bra019403 and Bra019410 are the most likely candidates for CRA3.7. Our research provides a valuable germplasm resource against P. brassicae Pb3 and CRA3.7 closely linked marker for marker-assisted selection in CR cultivars breeding.

Genome-wide identification and characterization of NBS-encoding genes in Raphanus sativus L. and their roles related to Fusarium oxysporum resistance.

BACKGROUND: The nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are important for plant development and disease resistance. Although genome-wide studies of NBS-encoding genes have been performed in several species, the evolution, structure, expression, and function of these genes remain unknown in radish (Raphanus sativus L.). A recently released draft R. sativus L. reference genome has facilitated the genome-wide identification and characterization of NBS-encoding genes in radish. RESULTS: A total of 225 NBS-encoding genes were identified in the radish genome based on the essential NB-ARC domain through HMM search and Pfam database, with 202 mapped onto nine chromosomes and the remaining 23 localized on different scaffolds. According to a gene structure analysis, we identified 99 NBS-LRR-type genes and 126 partial NBS-encoding genes. Additionally, 80 and 19 genes respectively encoded an N-terminal Toll/interleukin-like domain and a coiled-coil domain. Furthermore, 72% of the 202 NBS-encoding genes were grouped in 48 clusters distributed in 24 crucifer blocks on chromosomes. The U block on chromosomes R02, R04, and R08 had the most NBS-encoding genes (48), followed by the R (24), D (23), E (23), and F (17) blocks. These clusters were mostly homogeneous, containing NBS-encoding genes derived from a recent common ancestor. Tandem (15 events) and segmental (20 events) duplications were revealed in the NBS family. Comparative evolutionary analyses of orthologous genes among Arabidopsis thaliana, Brassica rapa, and Brassica oleracea reflected the importance of the NBS-LRR gene family during evolution. Moreover, examinations of cis-elements identified 70 major elements involved in responses to methyl jasmonate, abscisic acid, auxin, and salicylic acid. According to RNA-seq expression analyses, 75 NBS-encoding genes contributed to the resistance of radish to Fusarium wilt. A quantitative real-time PCR analysis revealed that RsTNL03 (Rs093020) and RsTNL09 (Rs042580) expression positively regulates radish resistance to Fusarium oxysporum, in contrast to the negative regulatory role for RsTNL06 (Rs053740). CONCLUSIONS: The NBS-encoding gene structures, tandem and segmental duplications, synteny, and expression profiles in radish were elucidated for the first time and compared with those of other Brassicaceae family members (A. thaliana, B. oleracea, and B. rapa) to clarify the evolution of the NBS gene family. These results may be useful for functionally characterizing NBS-encoding genes in radish.

QTL mapping for Fusarium wilt resistance based on the whole-genome resequencing and their association with functional genes in Raphanus sativus.

KEY MESSAGE: Two major QTL associated with resistance to Fusarium wilt (FW) were identified using whole-genome resequencing. Sequence variations and gene expression level differences suggest that TIR-NBS and LRR-RLK are candidate genes associated with FW-resistance. Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. raphani is an important disease in radish, leading to severe decrease in yield and quality. YR4 as a novel genetic source to resistant to FW was confirmed through screening with five pathogen isolates. We have generated F2 and F2:3 populations segregated with FW resistance using YR4 and YR18 inbred lines. The disease symptom was evaluated in F2:3 population (n = 180) in three independent studies over two years. We identified 4 QTL including the two major QTL (FoRsR7.159A and FoRsR9.359A). FoRsR7.159A and FoRsR9.359A were detected in three replicated experiments. FoRsR7.159A was delimited to the 2.18-Mb physical interval on chromosome R07, with a high LOD value (5.17-12.84) and explained phenotypic variation (9.34%-27.97%). The FoRsR9.359A represented relatively low LOD value (3.38-4.52) and explained phenotypic variation (6.24%-8.82%). On the basis of the re-sequencing data for the parental lines, we identified five putative resistance-related genes and 13 unknown genes with sequence variations at the gene and protein levels. A semi-quantitative RT-PCR analysis revealed that Rs382940 (TIR-NBS) and Rs382200 (RLK) were expressed only in 'YR4' from 0 to 6 days after the inoculation. Moreover, Rs382950 (TIR-NBS-LRR) was more highly expressed in 'YR4' from 3 to 6 days after the inoculation. These three genes might be important for FW-resistance in radish. We identified several markers based on these potential candidate genes. The marker set should be useful for breeding system to introduce the FW resistance loci from 'YR4' to improve tolerance to FW.

Genome-Wide Identification, Evolution, and Comparative Analysis of B-Box Genes in Brassica rapa, B. oleracea, and B. napus and Their Expression Profiling in B. rapa in Response to Multiple Hormones and Abiotic Stresses.

The B-box zinc-finger transcription factors are important for plant growth, development, and various physiological processes such as photomorphogenesis, light signaling, and flowering, as well as for several biotic and abiotic stress responses. However, there is relatively little information available regarding Brassica B-box genes and their expression. In this study, we identified 51, 52, and 101 non-redundant genes encoding B-box proteins in Brassica rapa (BrBBX genes), B. oleracea (BoBBX genes), and B. napus (BnBBX genes), respectively. A whole-genome identification, characterization, and evolutionary analysis (synteny and orthology) of the B-box gene families in the diploid species B. rapa (A genome) and B. oleracea (C genome) and in the allotetraploid species B. napus (AC genome) revealed segmental duplications were the major contributors to the expansion of the BrassicaBBX gene families. The BrassicaBBX genes were classified into five subgroups according to phylogenetic relationships, gene structures, and conserved domains. Light-responsive cis-regulatory elements were detected in many of the BBX gene promoters. Additionally, BrBBX expression profiles in different tissues and in response to various abiotic stresses (heat, cold, salt, and drought) or hormones (abscisic acid, methyl jasmonate, and gibberellic acid) were analyzed by qRT-PCR. The data indicated that many B-box genes (e.g., BrBBX13, BrBBX15, and BrBBX17) may contribute to plant development and growth as well as abiotic stress tolerance. Overall, the identified BBX genes may be useful as functional genetic markers for multiple stress responses and plant developmental processes.

Fine-mapping of a major QTL (Fwr1) for fusarium wilt resistance in radish.

KEY MESSAGE: A major radish QTL (Fwr1) for fusarium wilt resistance was fine-mapped. Sequence and expression analyses suggest that a gene encoding a serine/arginine-rich protein kinase is a candidate gene for Fwr1. Fusarium wilt resistance locus 1 (Fwr1) is a major quantitative trait locus (QTL) mediating the resistance of radish inbred line 'B2' to Fusarium oxysporum, which is responsible for fusarium wilt. We previously detected Fwr1 on radish linkage group 3 (i.e., chromosome 5). In this study, a high-resolution genetic map of the Fwr1 locus was constructed by analyzing 354 recombinant F2 plants derived from a cross between 'B2' and '835', the latter of which is susceptible to fusarium wilt. The Fwr1 QTL was fine-mapped to a 139.8-kb region between markers FM82 and FM87 in the middle part of chromosome 5. Fifteen candidate genes were predicted in this region based on a sequence comparison with the 'WK10039' radish reference genome. Additionally, we examined the time-course expression patterns of these predicted genes following an infection by the fusarium wilt pathogen. The ORF4 expression level was significantly higher in the resistant 'B2' plants than in the susceptible '835' plants. The ORF4 sequence was predicted to encode a serine/arginine-rich protein kinase and includes SNPs that result in nonsynonymous mutations, which may have important functional consequences. This study reveals a novel gene responsible for fusarium wilt resistance in radish. Further analyses of this gene may elucidate the molecular mechanisms underlying the fusarium wilt resistance of plants.

Quantitative Trait Locus Mapping of Clubroot Resistance and Plasmodiophora brassicae Pathotype Banglim-Specific Marker Development in Brassica rapa.

Clubroot resistance is an economically important trait in Brassicaceae crops. Although many quantitative trait loci (QTLs) for clubroot resistance have been identified in Brassica, disease-related damage continues to occur owing to differences in host variety and constant pathogen variation. Here, we investigated the inheritance of clubroot resistance in a double haploid population developed by crossing clubroot resistant and susceptible lines "09CR500" and "09CR501", respectively. The resistance of "09CR500" to Plasmodiophora brassicae pathotype "Banglim" was controlled as a single dominant gene, with the segregation of resistance and susceptibility being nearly 1:1. PbBrA08Banglim was identified as having a logarithm of odds value of 7.9-74.8, and a phenotypic variance of 26.0-97.1% with flanking marker "09CR.11390652" in A08. After aligning QTL regions to the B. rapa reference genome, 11 genes were selected as candidates. PbBrA08Banglim was located near Crr1, CRs, and Rcr9 loci, but differences were validated by marker analysis, gene structural variations, and gene expression levels, as well as phenotypic responses to the pathotype. Genotyping using the "09CR.11390652" marker accurately distinguished the Banglim-resistance phenotypes in the double haploid population. Thus, the developed marker will be useful in Brassica breeding programs, marker-assisted selection, and gene pyramiding to identify and develop resistant cultivars.

Genetic and physiological analyses of root cracking in radish (Raphanus sativus L.).

KEY MESSAGE: A major QTL conferring tolerance to radish (Raphanus sativus) root cracking was mapped for the first time and two calcium regulatory genes were identified that positively associated with the cracking phenomenon. Root cracking is a severe physiological disorder that significantly decreases the yield and commercial value of radish. The genetic and physiological mechanisms underlying this root cracking disorder have not been characterized. In this study, quantitative trait loci (QTLs) putatively associated with radish root cracking were mapped. Ten QTLs were distributed in six linkage groups, among these QTLs, 'RCr1' in LG1 was detected over 3 consecutive years and was considered to be a major QTL for root cracking. The QTL 'RCr1' was responsible for 4.47-18.11% of variance in the root cracking phenotype. We subsequently identified two candidate genes, RsANNAT and RsCDPK. Both genes encode proteins involved in calcium binding, ion transport, and Ca2+ signal transduction, which are important for regulating plant development and adaptations to the environment. These genes were co-localized to the major QTL region. Additionally, we analyzed physiological changes (i.e., root firmness, cell wall content, and cell-wall-bound calcium content) in two parental lines during different developmental stages. Moreover, we observed that the RsANNAT and RsCDPK expression levels are positively correlated with Ca2+ contents in the roots of the cracking-tolerant '835' cultivar. Thus, these genes may influence root cracking. The data provided herein may support the useful information to understand root cracking behavior in radish and may enable breeders to develop new cultivars exhibiting increased tolerance to root and fruit cracking.

RsGLK2.1-RsNF-YA9a module positively regulates the chlorophyll biosynthesis by activating RsHEMA2 in green taproot of radish.

Radish (Raphanus sativus L.) is an economically important and widely cultivated root vegetable crop. The coloration of the green skin and green flesh is an important trait influencing the nutrition and flavor quality in fruit radish. GOLDEN2-LIKEs (GLKs) play critically important roles in plastid development and chlorophyll biosynthesis in plants. However, the molecular mechanism underlying chlorophyll biosynthesis still remain elusive in green fruit radish taproot. Herein, the RsGLK2.1 gene exhibited higher expression level in taproot with a green skin (GS) and green flesh (GF) than that in taproot of the white or red radish genotypes. RsGLK2.1 is a nuclear transcription factor that has intrinsic transcriptional activation activity. Overexpression of RsGLK2.1 increased the total chlorophyll content of 20.68%-45.84% in radish leaves. Knockout of the RsGLK2.1 gene via CRISPR/Cas9 technology resulted in a significant decrease in the chlorophyll content. Overexpression of the RsGLK2.1 gene could restore the phenotype of the glk1glk2 mutant Arabidopsis. RsGLK2.1 was participated in regulating the chlorophyll biosynthesis by directly binding to the promoter of RsHEMA2 and activating its transcription. The interaction of RsNF-YA9a with RsGLK2.1 increased the transcriptional activity of the downstream gene RsHEMA2 under the light condition rather than the dark condition, indicating that both of them regulate the chlorophyll biosynthesis in a light-dependent manner of radish. Overall, these results provided insights into the molecular framework of the RsGLK2.1-RsNF-YA9a module, and could facilitate dissecting the regulatory mechanism underlying chlorophyll biosynthesis in green taproot of radish, and genetic improvement of quality traits in fruit radish breeding programs.

RsERF40 contributes to cold stress tolerance and cell expansion of taproot in radish (Raphanus sativus L.).

The growth and development of taproots are inhibited by cold stress in radish (Raphanus sativus L.). Ethylene-responsive element binding factors (ERF) are key participators in the cold stress response and growth regulation of plants. However, the function of ERF genes in cold tolerance and root development in radish remains elusive. Here, we showed that the secondary growth of radish taproots was inhibited by cold stress. Comparative transcriptome analysis demonstrated that the RsERF40 gene is an important regulator of the cold stress response and root growth regulation. The cold tolerance of transgenic Arabidopsis plants overexpressing the RsERF40 gene was significantly improved. Overexpressing RsERF40 in the cold-sensitive radish genotype and silencing RsERF40 in the cold-tolerant radish genotype indicated that RsERF40 was beneficial for alleviating oxidative damage under cold stress in radish. Transgenic Arabidopsis seedlings showed an increase in the elongation and radial growth of dark-grown roots. RT-qPCR analysis showed that the expression of the cold-related genes (CORs) RsCOR78 and RsCOR413PM1 and the cell wall strengthening-related genes RsCESA6 and RsEXPB3 was upregulated in transgenic Arabidopsis seedlings. Yeast one-hybrid (Y1H) and dual-luciferase reporter assays (DLA) revealed that RsERF40 directly regulates RsCOR78, RsCOR413PM1, RsCESA6 and RsEXPB3 expression, illustrating that RsERF40 enhances cold tolerance and taproot growth by modulating osmotic adjustment and cell wall mechanical strength in radish. In this study, the RsERF40-regulon was firstly found to be a new cold response pathway independent of the CBF-COR pathway conferring cold stress tolerance with increasing radish taproot growth. These results provided novel insight into the molecular mechanism underlying cold stress response and would facilitate the genetic improvement of cold tolerance in radish and other root vegetable crops.

RsCDF3, a member of Cycling Dof Factors, positively regulates cold tolerance via auto-regulation and repressing two RsRbohs transcription in radish (Raphanus sativus L.).

Radish is one of the most economical root vegetable crops worldwide. Cold stress dramatically impedes radish taproot formation and development as well as reduces its yield and quality. Although the Cycling Dof Factors (CDFs) play crucial roles in plant growth, development and abiotic stress responses, how CDF TFs mediate the regulatory network of cold stress response remains largely unexplored in radish. Herein, a total of nine RsCDF genes were identified from the radish genome. Among them, the RsCDF3 exhibited obviously up-regulated expression under cold stress, especially at 12 h and 24 h. RsCDF3 was localized to the nucleus and displayed dramatic cold-induced promoter activity in tobacco leaves. Moreover, overexpression of RsCDF3 significantly enhanced cold tolerance of radish plants, whereas its knock-down plants exhibited the opposite phenotype. Interestingly, both in vitro and in vivo assays indicated that the RsCDF3 repressed the transcription of RsRbohA and RsRbohC via directly binding to their promoters, which contributed to maintaining the cellular homeostasis of reactive oxygen species (ROS) production and scavenging in radish. In addition, the RsCDF3 bound to its own promoter to mediate its transcription, thereby forming an autoregulatory feedback loop to cooperatively trigger RsRbohs-dependent cold tolerance. Together, we revealed a novel RsCDF3-RsRbohs module to promote the cold tolerance in radish plants. These findings would facilitate unveiling the molecular mechanism governing RsCDF3-mediated cold stress response in radish.

[Study on plant regeneration from somatic embryos of vulnerable medicinal plant Glehnia littoralis].

OBJECTIVE: To study the cause of the seeds dormancy of Glehnia littoralis in vitro and to establish plant regeneration methods via somatic embryos. METHOD: The effects of endosperm and exogenous hormone on the seed dormancy breaking of G. littoralis and the effect of hormone concentration on embryonic callus induction and plant regeneration via somatic embryos were observed, RESULTS: The germination rate of the seeds with 1/3 endosperm was the highest which achieved 31%. TDZ, 6-BA and GA3 treatment could not break seed dormancy but easily lead to abnormal seedlings. Embryogenic callus induction rates was up to 57% on MS supplemented with 1.0 mg x L(-1) 2,4-D. After 20 days culture, embryogenic calli were transferred to MS medium and cotyledonary embryos were formed in 40 days. The regenerated plants were obtained in 20 days. CONCLUSION: An effective system of plant regeneration of G. littoralis was established in this study.

Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection.

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection on Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at seven days after the P. brassicae inoculation, whereas the expression levels of the starch degradation-related genes increased at 14 days after the inoculation. Then genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization as indicated by the histological observation and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

Transcriptome Sequencing and Chemical Analysis Reveal the Formation Mechanism of White Florets in Carthamus tinctorius L.

Carthamus tinctorius L. (safflower), an economic crop and herb, has been extensively studied for its diverse chemical constituents and pharmacological effects, but the mechanism of safflower pigments (SP) leading to different colors of florets has not been clarified. In the present study, we compared the contents of SP in two varieties of safflower with white and red florets, named Xinhonghua No. 7 (WXHH) and Yunhong No. 2 (RYH). The results showed the contents of SP in RYH were higher than WXHH. To investigate genes related to SP, we obtained six cDNA libraries of florets from the two varieties by transcriptome sequencing. A total of 225,008 unigenes were assembled and 40 unigenes related to safflower pigment biosynthesis were annotated, including 7 unigenes of phenylalanine ammonia-lyase (PAL), 20 unigenes of 4-coumarate-CoA ligase (4CL), 1 unigene of trans-cinnamate 4-monooxygenase (C4H), 7 unigenes of chalcone synthase (CHS), 4 unigenes of chalcone isomerase (CHI), and 1 unigene of flavanone 3-hydroxylase (F3H). Based on expression levels we selected 16 differentially expressed unigenes (DEGs) and tested them using reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR), which was consistent with the sequencing results. Consequently, we speculated that in WXHH, 3 PALs, 3 4CLs, 1 C4H, 1 CHS, and 1 CHI, which were down-regulated, and 1 F3H, which was up-regulated, may play a key role in the formation of white florets.

Isolation and characterization of fusarium wilt resistance gene analogs in radish.

The resistance gene analog (RGA)-based marker strategy is an effective supplement for current marker reservoir of radish disease-resistance breeding. In this study, we identified RGAs based on the conserved nucleotide-binding site (NBS) and S-receptor-like kinase (SRLK) domains. A total of 68 NBS-RGAs and 46 SRLK-RGAs were isolated from two FW-resistant radish inbred lines, B2 and YR31, and one susceptible line, YR15. A BLASTx search revealed that the NBS-RGAs contained six conserved motifs (i.e., P loop, RNBS-A, Kinase-2, RNBS-B, RNBS-C, and GLPL) and the SRLK-RGAs, contained two conserved motifs (i.e., G-type lectin and PAN-AP). A phylogenetic analysis indicated that the NBS-RGAs could be separated into two classes (i.e., toll/interleukin receptor and coiled-coil types), with six subgroups, and the SRLK-RGAs were divided into three subgroups. Moreover, we designed RGA-specific markers from data-mining approach in radish databases. Based on marker analysis, 24 radish inbred lines were clustered into five main groups with a similarity index of 0.44 and showing genetic diversity with resistance variation in those radish inbred lines. The development of RGA-specific primers would be valuable for marker-assisted selection during the breeding of disease-resistant radish cultivars.