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.

Black rot disease, caused by Xanthomonas campestris pv. campestris (Xcc), results in significant yield losses in Brassica oleracea crops worldwide. To find black rot disease-resistant cabbage lines, we carried out pathogenicity assays using the scissor-clipping method in 94 different B. oleracea lines. By comparing the lesion areas, we selected a relatively resistant line, Black rot Resistance 155 (BR155), and a highly susceptible line, SC31. We compared the two cabbage lines for the Xcc-induced expression pattern of 13 defense-related genes. Among them, the Xcc-induced expression level of PR1 and antioxidant-related genes (SOD, POD, APX, Trx H, and CHI) were more than two times higher in BR155 than SC31. Nitroblue tetrazolium (NBT) and diaminobenzidine tetrahydrochloride (DAB) staining analysis showed that BR155 accumulated less Xcc-induced reactive oxygen species (ROS) than did SC31. In addition, 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging assays showed that BR155 had higher antioxidant activity than SC31. This study, focused on the defense responses of cabbage during the early biotrophic stage of infection, indicated that Xcc-induced ROS might play a role in black rot disease development. We suggest that non-enzymatic antioxidants are important, particularly in the early defense mechanisms of cabbage against Xcc.

The unwanted detachment of organs such as flowers, leaves, and fruits from the main body of a plant (abscission) has significant effects on agricultural practice. Both timely and precise regulation of organ abscission from a plant is crucial as it influences the agricultural yield. The tomato (Solanum lycopersicum) has become a model system for research on organ abscission. Here, we characterized four tomato natural abscission variants named jointless (j), functionally impaired jointless (fij), functionally impaired jointless like (fij like), and normal joint (NJ), based on their cellular features within the flower abscission zones (AZ). Using eight INFLORESCENCE DEFICIENT IN ABSCISSION (SlIDA) genes and eight HAESA genes (SlHAE) identified in the genome sequence of tomato, we analyzed the pattern of gene expression during flower abscission. The AZ-specific expression for three tomato abscission polygalacturonases (SlTAPGs) in the development of flower AZ, and the progression of abscission validated our natural abscission system. Compared to that of j, fij, and fij like variants, the AZ-specific expression for SlIDA, SlIDL2, SlIDL3, SlIDL4, and SlIDL5 in the NJ largely corelated and increased with the process of abscission. Of eight SlHAE genes examined, the expression for SlHSL6 and SlHSL7 were found to be AZ-specific and increased as abscission progressed in the NJ variant. Unlike the result of gene expression obtained from natural abscission system, an in silico analysis of transcriptional binding sites uncovered that SlIDA genes (SlIDA, SlIDL6, and SlIDL7) are predominantly under the control of environmental stress, while most of the SlHSL genes are affiliated with the broader context in developmental processes and stress responses. Our result presents the potential bimodal transcriptional regulation of the tomato IDA-HAE module associated with flower abscission in tomatoes.

Leaf senescence is a developmental process induced by various molecular and environmental stimuli that may affect crop yield. The dark-induced leaf senescence-91 (DLS-91) plants displayed rapid leaf senescence, dramatically decreased chlorophyll contents, low photochemical efficiencies, and upregulation of the senescence-associated marker gene BrSAG12-1. To understand DLS molecular mechanism, we examined transcriptomic changes in DLS-91 and control line DLS-42 following 0, 1, and 4 days of dark treatment (DDT) stages. We identified 501, 446, and 456 DEGs, of which 16.7%, 17.2%, and 14.4% encoded TFs, in samples from the three stages. qRT-PCR validation of 16 genes, namely, 7 MADS, 6 NAC, and 3 WRKY, suggested that BrAGL8-1, BrAGL15-1, and BrWRKY70-1 contribute to the rapid leaf senescence of DLS-91 before (0 DDT) and after (1 and 4 DDT) dark treatment, whereas BrNAC046-2, BrNAC029-2/BrNAP, and BrNAC092-1/ORE1 TFs may regulate this process at a later stage (4 DDT). In-silico analysis of cis-acting regulatory elements of BrAGL8-1, BrAGL42-1, BrNAC029-2, BrNAC092-1, and BrWRKY70-3 of B. rapa provides insight into the regulation of these genes. Our study has uncovered several AGL-MADS, WRKY, and NAC TFs potentially worthy of further study to understand the underlying mechanism of rapid DLS in DLS-91.