Barbarea vulgaris is a wild cruciferous plant and include two distinct types: the G- and P-types named after their glabrous and pubescent leaves, respectively. The types differ significantly in resistance to a range of insects and diseases as well as glucosinolates and other chemical defenses. A high-density linkage map was needed for further progress to be made in the molecular research of this plant. We performed restriction site-associated DNA sequencing (RAD-Seq) on an F2 population generated from G- and P-type B. vulgaris. A total of 1545 SNP markers were mapped and ordered in eight linkage groups, which represents the highest density linkage map to date for the crucifer tribe Cardamineae. A total of 722 previously published genome contigs (50.2 Mb, 30% of the total length) can be anchored to this high density genetic map, an improvement compared to a previously published map (431 anchored contigs, 38.7 Mb, 23% of the assembly genome). Most of these (572 contigs, 31.2 Mb) were newly anchored to the map, representing a significant improvement. On the basis of the present high-density genetic map, 37 QTL were detected for eleven traits, each QTL explaining 2.9–71.3% of the phenotype variation. QTL of glucosinolates, leaf size and color traits were in most cases overlapping, possibly implying a functional connection. This high-density linkage map and the QTL obtained in this study will be useful for further understanding of the genetic of the B. vulgaris and molecular basis of these traits, many of which are shared in the related crop watercress.

Barbarea vulgaris contains two genotypes: the glabrous type (G-type), which confers resistance to the diamondback moth (DBM) and other insect pests, and the pubescent type (P-type), which is susceptible to the DBM. Herein, the transcriptomes of P-type B. vulgaris before and after DBM infestation were subjected to Illumina (Solexa) pyrosequencing and comparative analysis. 5.0 gigabase pairs of clean nucleotides were generated. Non-redundant unigenes (33,721) were assembled and 94.1 % of them were annotated. Compared with our previous G-type transcriptome, the expression patterns of many insect responsive genes, including those related to secondary metabolism, phytohormones and transcription factors, which were significantly induced by DBM in G-type plants, were less sensitive to DBM infestation in P-type plants. The genes of the triterpenoid saponin pathway were identified in both G- and P-type plants. The upstream genes of the pathway showed similar expression patterns between the two genotypes. However, gene expression for two downstream enzymes, the glucosyl transferase (UGT73C11) and an oxidosqualene cyclase (OSC), were significantly upregulated in the P-type compared with the G-type plant. The homologous genes from P- and G-type plants were detected by BLAST unigenes with a cutoff level E-value < e−10. 12,980 gene families containing 26,793 P-type and 36,944 G-type unigenes were shared by the two types of B. vulgaris. 38,397 single nucleotide polymorphisms (SNPs) were found in 9,452 orthologous genes between the P- and G-type plants. We also detected 5,105 simple sequence repeats (SSRs) in the B. vulgaris transcriptome, comprising mono-nucleotide-repeats (2,477; 48.5 %) and triple-nucleotide-repeats (1,590; 31.1 %). Of these, 1,657 SSRs displayed polymorphisms between the P- and G-type. Consequently, 913 SSR primer pairs were designed with a resolution of more than two nucleotides. We randomly chose 30 SSRs to detect the genetic diversity of 32 Barbarea germplasms. The distance tree showed that these accessions were clearly divided into groups, with the G-type grouping with available Western and Central European B. vulgaris accessions in contrast to the P-type accession, B. stricta and B. verna. These data represent useful information for pest-resistance gene mining and for the investigation of the molecular basis of plant-pest interactions.

Background: Sugar not only is an important biomacromolecule that plays important roles in plant growth, development, and biotic and abiotic stress tolerance but also provides a skeleton for other macromolecules, such as proteins and nucleic acids. Sugar transporter proteins (STPs) play essential roles in plant sugar transport and ultimately affect the abovementioned life processes. However, the evolutionary dynamics of this important gene family in Brassicaceae crops are still largely unknown, and the functional differentiation of radish STP genes remains unclear. Results: In the present study, a comparative genomic study of STP genes in five representative Brassicaceae crops was conducted, and a total of 25, 25, 28, 36 and 49 STP genes were individually identified in Raphanus sativus (Rs), Brassica oleracea (Bo), B. rapa (Br), B. napus (Bn) and B. juncea (Bj), which were divided into four clades by phylogenetic analysis. The number of STP genes was no direct correlation with genome size and the total number of coding genes in Brassicaceae crops, and their physical and chemical properties showed no significant difference. Expression analysis showed that radish STP genes play vital roles not only in flower and seedpod development but also under heavy metal (cadmium, chromium and lead), NaCl and PEG-6000 stresses, Agrobacterium tumefaciens infection, and exogenous sugar treatment. RsSTP13.2 was significantly upregulated in the resistant radish cultivar by A. tumefaciens infection and induced by heavy metal, NaCl and PEG-6000 stress, indicating that it is involved in resistance to both biotic and abiotic stress in radish. Conclusions: The present study provides insights into the evolutionary patterns of the STP gene family in Brassicaceae genomes and provides a theoretical basis for future functional analysis of STP genes in Brassicaceae crops.

Radish (Raphanus sativus L.) is rich in specific glucosinolates (GSLs), which benefit human health and special flavor formation. Although the basic GSLs metabolic pathway in Brassicaceae plants is clear, the regulating mechanism for specific glucosinolates content in radish fleshy taproots is not well understood. In this study, we discovered that there was a significant difference in the GSLs profiles and the content of various GSLs components. Glucoraphasatin (GRH) is the most predominant GSL in radish taproots of different genotypes as assessed by HPLC analysis. Further, we compared the taproot transcriptomes of three radish genotypes with high and low GSLs content by employing RNA-seq. Totally, we identified forty-one differentially expressed genes related to GSLs metabolism. Among them, thirteen genes (RsBCAT4, RsIPMDH1, RsMAM1a, RsMAM1b, RsCYP79F1, RsGSTF9, RsGGP1, RsSUR1, RsUGT74C1, RsST5b, RsAPK1, RsGSL-OH, and RsMYB28) were significantly higher co-expressed in the high content genotypes than in low content genotype. Notably, correlation analysis indicated that the expression level of RsMYB28, as an R2R3 transcription factor directly regulating aliphatic glucosinolate biosynthesis, was positively correlated with the GRH content. Co-expression network showed that RsMYB28 probably positively regulated the expression of the above genes, particularly RsSUR1, and consequently the synthesis of GRH. Moreover, the molecular mechanism of the accumulation of this 4-carbon (4C) GSL in radish taproots was explored. This study provides new perspectives on the GSLs accumulation mechanism and genetic improvements in radish taproots.