Three genes encoding AOP2, a protein involved in aliphatic glucosinolate biosynthesis, are differentially expressed in Brassica rapa

Metabolomic and Transcriptomic Profiles in Diverse Brassica oleracea Crops Provide Insights into the Genetic Regulation of Glucosinolate Profiles

Glucosinolates (GSLs) are plant secondary metabolites commonly found in the cruciferous vegetables of the Brassicaceae family, offering health benefits to humans and defense against pathogens and pests to plants. In this study, we investigated 23 GSL compounds' relative abundance in four tissues of five different Brassica oleracea morphotypes. Using the five corresponding high-quality B. oleracea genome assemblies, we identified 183 GSL-related genes and analyzed their expression with mRNA-Seq data. GSL abundance and composition varied strongly, among both tissues and morphotypes, accompanied by different gene expression patterns. Interestingly, broccoli exhibited a nonfunctional AOP2 gene due to a conserved 2OG-FeII_Oxy domain loss, explaining the unique accumulation of two health-promoting GSLs. Additionally, transposable element (TE) insertions were found to affect the gene structure of MAM3 genes. Our findings deepen the understanding of GSL variation and genetic regulation in B. oleracea morphotypes, providing valuable insights for breeding with tailored GSL profiles in these crops.

Chromosome-scale reference genome of broccoli (Brassica oleracea var. italica Plenck) provides insights into glucosinolate biosynthesis

Broccoli (Brassica oleracea var. italica Plenck) is an important vegetable crop, as it is rich in health-beneficial glucosinolates (GSLs). However, the genetic basis of the GSL diversity in Brassicaceae remains unclear. Here we report a chromosome-level genome assembly of broccoli generated using PacBio HiFi reads and Hi-C technology. The final genome assembly is 613.79 Mb in size, with a contig N50 of 14.70 Mb. The GSL profile and content analysis of different B. oleracea varieties, combined with a phylogenetic tree analysis, sequence alignment, and the construction of a 3D model of the methylthioalkylmalate synthase 1 (MAM1) protein, revealed that the gene copy number and amino acid sequence variation both contributed to the diversity of GSL biosynthesis in B. oleracea. The overexpression of BoMAM1 (BolI0108790) in broccoli resulted in high accumulation and a high ratio of C4-GSLs, demonstrating that BoMAM1 is the key enzyme in C4-GSL biosynthesis. These results provide valuable insights for future genetic studies and nutritive component applications of Brassica crops.

CRISPR/Cas9-mediated BoaAOP2s editing alters aliphatic glucosinolate side-chain metabolic flux and increases the glucoraphanin content in Chinese kale

Glucoraphanin (GRA) is an aliphatic glucosinolate (GSL), and its hydrolysis product has powerful anticancer activity. ALKENYL HYDROXALKYL PRODUCING 2 (AOP2) gene, encodes a 2-oxoglutarate-dependent dioxygenase, which can catalyze GRA to form gluconapin (GNA). However, GRA only present in trace amounts in Chinese kale. To increase the content of GRA in Chinese kale, three copies of BoaAOP2 were isolated and edited using CRISPR/Cas9 system. The content of GRA was 11.71- to 41.29-fold (0.082-0.289 μmol g^-1 FW) higher in T1 generation of boaaop2 mutants than in wild-type plants, and this was accompanied by an increase in the GRA/GNA ratio and reductions in the content of GNA and total aliphatic GSLs. BoaAOP2.1 is an effective gene for the alkenylation of aliphatic GSLs in Chinese kale. Overall, targeted editing of CRISPR/Cas9-mediated BoaAOP2s altered aliphatic GSL side-chain metabolic flux and enhanced the GRA content in Chinese kale, suggesting that metabolic engineering of BoaAOP2s has huge potential in improving nutritional quality of Chinese kale.

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

Glucosinolate Biosynthetic Genes of Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis

Cabbage (Brassica oleracea var. capitata) is a vegetable rich in glucosinolates (GSLs) that have proven health benefits. To gain insights into the synthesis of GSLs in cabbage, we systematically analyzed GSLs biosynthetic genes (GBGs) in the entire cabbage genome. In total, 193 cabbage GBGs were identified, which were homologous to 106 GBGs in Arabidopsis thaliana. Most GBGs in cabbage have undergone negative selection. Many homologous GBGs in cabbage and Chinese cabbage differed in expression patterns indicating the unique functions of these homologous GBGs. Spraying five exogenous hormones significantly altered expression levels of GBGs in cabbage. For example, MeJA significantly upregulated side chain extension genes BoIPMILSU1-1 and BoBCAT-3-1, and the expression of core structure construction genes BoCYP83A1 and BoST5C-1, while ETH significantly repressed the expression of side chain extension genes such as BoIPMILSU1-1, BoCYP79B2-1, and BoMAMI-1, and some transcription factors, namely BoMYB28-1, BoMYB34-1, BoMYB76-1, BoCYP79B2-1, and BoMAMI-1. Phylogenetically, the CYP83 family and CYP79B and CYP79F subfamilies may only be involved in GSL synthesis in cruciferous plants. Our unprecedented identification and analysis of GBGs in cabbage at the genome-wide level lays a foundation for the regulation of GSLs synthesis through gene editing and overexpression.

Transcriptome and QTL mapping analyses of major QTL genes controlling glucosinolate contents in vegetable- and oilseed-type Brassica rapa plants

Glucosinolates (GSLs) are secondary metabolites providing defense against pathogens and herbivores in plants, and anti-carcinogenic activity against human cancer cells. Profiles of GSLs vary greatly among members of genus Brassica. In this study, we found that a reference line of Chinese cabbage (B. rapa ssp. pekinensis), 'Chiifu' contains significantly lower amounts of total GSLs than the oilseed-type B. rapa (B. rapa ssp. trilocularis) line 'LP08'. This study aimed to identify the key regulators of the high accumulation of GSLs in Brassica rapa plants using transcriptomic and linkage mapping approaches. Comparative transcriptome analysis showed that, in total, 8,276 and 9,878 genes were differentially expressed between 'Chiifu' and 'LP08' under light and dark conditions, respectively. Among 162 B. rapa GSL pathway genes, 79 were related to GSL metabolism under light conditions. We also performed QTL analysis using a single nucleotide polymorphism-based linkage map constructed using 151 F₅ individuals derived from a cross between the 'Chiifu' and 'LP08' inbred lines. Two major QTL peaks were successfully identified on chromosome 3 using high-performance liquid chromatography to obtain GSL profiles from 97 F₅ recombinant inbred lines. The MYB-domain transcription factor gene BrMYB28.1 (Bra012961) was found in the highest QTL peak region. The second highest peak was located near the 2-oxoacid-dependent dioxygenase gene BrGSL-OH.1 (Bra022920). This study identified major genes responsible for differing profiles of GSLs between 'Chiifu' and 'LP08'. Thus, our study provides molecular insights into differences in GSL profiles between vegetative- and oilseed-type B. rapa plants.

Transcriptome analysis of genes related to glucoraphanin and sulforaphane synthesis in methyl jasmonate treated broccoli (Brassica oleracea var. italica) hairy roots

Hairy roots obtained by infecting broccoli (Brassica oleracea var. italica) leaves with Agrobacterium rhizogenes (ATCC15834) have the characteristics of phytohormone autonomy, genetic stability and can produce a large amount of the anti-cancer substance Sulforaphane (SF) and the biosynthetic precursor Glucoraphanin (GRA). Under the induction of the exogenous signaling molecule methyl jasmonate (MeJA), the production of SF in broccoli hairy roots was significantly increased. However, the molecular mechanism of GRA and SF synthesis in hairy roots of broccoli treated with MeJA has not been reported. In this study, according to the yield of GRA and SF, the best concentration of MeJA treatment for hairy roots of broccoli was selected. After 18 days of growth, broccoli hairy roots were treated with 10 mmol L^-1 MeJA for 0, 3, 6, 9 and 12 h. Compared with 0 h, the yield of GRA and SF increased under other treatments. The highest yield of GRA and SF occurred at 9 h, which were 2.22-fold and 1.74-fold higher than those at 0 h. Brassica oleracea var. botrytis was used as reference genome, and 5,757 differentially expressed genes (DEG) were observed at 0, 3, 6, 9 and 12 h under 10 mmol L^-1 MeJA treatment, of which 4,673 were down-regulated and 1084 were up-regulated. The key genes regulating GRA synthesis, CYP79F1, CYP83A1, UGT74B1, FMOGS-OX5 and GSL-OH, were up-regulated at 0 and 3 h, and down-regulated the rest of the time; BCAT2 was up-regulated at 6, 9, 12 h, and at 0, 3 h expression was down-regulated, transcription factors MYB28 and MYB29 were down-regulated by exogenous MeJA treatment. A pathway of GRA biosynthesis and transformation pathways in MeJA-treated broccoli hairy roots was simulated and the molecular mechanism of GRA biosynthesis and SF accumulation in broccoli hairy roots under MeJA treatment was revealed.

Mapping and Validation of BrGOLDEN: A Dominant Gene Regulating Carotenoid Accumulation in Brassica rapa

In plants, the accumulation of carotenoids can maintain the balance of the photosystem and improve crop nutritional quality. Therefore, the molecular mechanisms underlying carotenoid synthesis and accumulation should be further explored. In this study, carotenoid accumulation differed significantly among parental Brassica rapa. Genetic analysis was carried out using the golden inner leaf ‘1900264′ line and the light−yellow inner leaf ‘1900262′ line, showing that the golden inner leaf phenotype was controlled by a single dominant gene. Using bulked−segregant analysis sequencing, BraA09g007080.3C encoding the ORANGE protein was selected as a candidate gene. Sequence alignment revealed that a 4.67 kb long terminal repeat insertion in the third exon of the BrGOLDEN resulted in three alternatively spliced transcripts. The spatiotemporal expression results indicated that BrGOLDEN might regulate the expression levels of carotenoid−synthesis−related genes. After transforming BrGOLDEN into Arabidopsis thaliana, the seed−derived callus showed that BrGOLDEN_Ins and BrGOLDEN_Del lines presented a yellow color and the BrGOLDEN_Ldel line presented a transparent phenotype. In addition, using the yeast two−hybrid assay, BrGOLDEN_Ins, BrGOLDEN_Ldel, and Brgolden_wt exhibited strong interactions with BrPSY1, but BrGOLDEN_Del did not interact with BrPSY1 in the split−ubiquitin membrane system. In the secondary and 3D structure analysis, BrGOLDEN_Del was shown to have lost the PNFPSFIPFLPPL sequences at the 125 amino acid position, which resulted in the α−helices of BrGOLDEN_Del being disrupted, restricting the formation of the 3D structure and affecting the functions of the protein. These findings may provide new insights into the regulation of carotenoid synthesis in B. rapa.

Glucosinolates and Biotic Stress Tolerance in Brassicaceae with Emphasis on Cabbage: A Review

Glucosinolates (GSLs) and GSL-associated genes are receiving increasing attention from molecular biologists due to their multifunctional properties. GSLs are secondary metabolites considered to be highly active in most Brassica species. Their importance has motivated the discovery and functional analysis of the GSLs and GSL hydrolysis products involved in disease development in brassicas and other plants. Comprehensive knowledge of the GSL content of Brassica species and the molecular details of GSL-related genes will help elucidate the molecular control of this plant defense system. This report provides an overview of the current status of knowledge on GSLs, GSL biosynthesis, as well as hydrolysis related genes, and GSL hydrolysis products that regulate fungal, bacterial, and insect resistance in cabbage and other brassicas.

Transcriptome analysis of melatonin regulating the transformation of glucoraphanin to sulforaphane in broccoli hairy roots

Sulforaphane (SF) is one of the most effective natural products in preventing and fighting cancer, found in cruciferous plants. In this study, broccoli hairy roots grown for 20 d were used as the experimental material, and it was treated with 500 μmol/L melatonin (MT) for 0, 12 and 32 h to explore the effect of MT on the conversion of glucoraphanin (GRA) to SF. Results showed that the yields of GRA and SF were the largest under MT treatment for 12 h, which were 1.53 and 1.93-fold, respectively, compared to 0 h. However, Myrosinases activity was the highest under MT treatment for 32 h, which was 1.42-fold compared to that of the 0 h. The differential expression of key genes involved in GRA conversion to SF in broccoli hairy roots was identified transcriptome sequencing, and the path of the transformation from GRA to SF was simulated, which provided a theoretical basis for establishing an efficient transformation system from GRA to SF.

Transcriptomics analysis of genes induced by melatonin related to glucosinolates synthesis in broccoli hairy roots

Glucoraphanin (GRA) is found in the seeds and vegetative organs of broccoli (Brassica oleracea L. var. italica Planch) as the precursor of anti-carcinogen sulforaphane (SF). The yield of GRA obtained from these materials is weak and the cost is high. In recent years, the production of plant secondary metabolites by large-scale hairy roots culture in vitro has succeeded in some species. Melatonin (MT) is a natural hormone which existed in numerous organisms. Studies have demonstrated that MT can improve the synthesis of secondary metabolites in plants. At present, it has not been reported that MT regulates the biosynthesis of glucoraphanin in broccoli hairy roots. In this study, the broccoli hairy roots that grew for 20 d were respectively treated by 500 µM MT for 0, 6, 12, 20 and 32. To explore the reason of changes in secondary metabolites and reveal the biosynthetic pathway of glucoraphanin at transcriptional level. Compared with 0 h, the yield of GRA under other treatments was increased, and the overall trend was firstly increased and then decreased. The total yield of GRA reached the highest at 12 h, which was 1.22-fold of 0 h. Then, the genome of broccoli as the reference, a total of 13234 differentially expressed genes (DEGs) were identified in broccoli hairy roots under treatment with 500 µM MT for 0, 6, 12, 20 and 32 h, respectively. Among these DEGs, 6266 (47.35%) were upregulated and 6968 (52.65%) were downregulated. It was found that the pathway of 'Glucosinolates biosynthesis (ko00966)' was enriched in the 16th place by Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of the upregulated DEGs. The expression of key genes in the GRA biosynthesis pathway was upregulated at all time points, and a deduced GRA biosynthesis pathway map was constructed for reference.

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa

This study was conducted to investigate doubled haploid (DH) lines produced between high GSL (HGSL) Brassica rapa ssp. trilocularis (yellow sarson) and low GSL (LGSL) B. rapa ssp. chinensis (pak choi) parents. In total, 161 DH lines were generated. GSL content of HGSL DH lines ranged from 44.12 to 57.04 μmol·g⁻¹·dry weight (dw), which is within the level of high GSL B. rapa ssp. trilocularis (47.46 to 59.56 μmol g⁻¹ dw). We resequenced five of the HGSL DH lines and three of the LGSL DH lines. Recombination blocks were formed between the parental and DH lines with 108,328 single-nucleotide polymorphisms in all chromosomes. In the measured GSL, gluconapin occurred as the major substrate in HGSL DH lines. Among the HGSL DH lines, BrYSP_DH005 had glucoraphanin levels approximately 12-fold higher than those of the HGSL mother plant. The hydrolysis capacity of GSL was analyzed in HGSL DH lines with a Korean pak choi cultivar as a control. Bioactive compounds, such as 3-butenyl isothiocyanate, 4-pentenyl isothiocyanate, 2-phenethyl isothiocyanate, and sulforaphane, were present in the HGSL DH lines at 3-fold to 6.3-fold higher levels compared to the commercial cultivar. The selected HGSL DH lines, resequencing data, and SNP identification were utilized for genome-assisted selection to develop elite GSL-enriched cultivars and the industrial production of potential anti-cancerous metabolites such as gluconapin and glucoraphanin.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

Divergence of three BRX homoeologs in Brassica rapa and its effect on leaf morphology

The leafy head characteristic is a special phenotype of Chinese cabbage resulting from artificial selection during domestication and breeding. BREVIS RADIX (BRX) has been suggested to control root elongation, shoot growth, and tiller angle in Arabidopsis and rice. In Brassica rapa, three BrBRX homoeologs have been identified, but only BrBRX.1 and BrBRX.2 were found to be under selection in leaf-heading accessions, indicating their functional diversification in leafy head formation. Here, we show that these three BrBRX genes belong to a plant-specific BRX gene family but that they have significantly diverged from other BRX-like members on the basis of different phylogenetic classifications, motif compositions and expression patterns. Moreover, although the expression of these three BrBRX genes differed, compared with BrBRX.3, BrBRX.1, and BrBRX.2 displayed similar expression patterns. Arabidopsis mutant complementation studies showed that only BrBRX.1 could rescue the brx root phenotype, whereas BrBRX.2 and BrBRX.3 could not. However, overexpression of each of the three BrBRX genes in Arabidopsis resulted in similar pleiotropic leaf phenotypes, including epinastic leaf morphology, with an increase in leaf number and leaf petiole length and a reduction in leaf angle. These leaf traits are associated with leafy head formation. Further testing of a SNP (T/C) in BrBRX.2 confirmed that this allele in the heading accessions was strongly associated with the leaf-heading trait of B. rapa. Our results revealed that all three BrBRX genes may be involved in the leaf-heading trait, but they may have functionally diverged on the basis of their differential expression.

Glucosinolate Profile and Glucosinolate Biosynthesis and Breakdown Gene Expression Manifested by Black Rot Disease Infection in Cabbage

Cabbage (Brassica oleracea var. capitata) is an economically important crop in the family Brassicaceae. Black rot disease is a top ranked cabbage disease, which is caused by Xanthomonas campestris pv. campestris (Xcc) and may reduce 50% crop loss. Therefore, we need a clear understanding of black rot disease resistance for sustainable disease management. The secondary metabolites, like Glucosinolate (GSL) presents in Brassica species, which plays a potential role in the defense mechanism against pathogens. However, there is little known about GSL-regulated resistance mechanisms and GSL biosynthesis and the breakdown related gene expression after black rot disease infection in cabbage. In this study, relative expression of 43 biosynthetic and breakdown related GSLs were estimated in the black rot resistant and susceptible cabbage lines after Xcc inoculation. Ten different types of GSL from both aliphatic and indolic groups were identified in the contrasting cabbage lines by HPLC analysis, which included six aliphatic and four indolic compounds. In the resistant line, nine genes (MYB122-Bol026204, MYB34-Bol017062, AOP2-Bo9g006240, ST5c-Bol030757, CYP81F1-Bol017376, CYP81F2-Bol012237, CYP81F4-Bol032712, CYP81F4-Bol032714 and PEN2-Bol030092) showed consistent expression patterns. Pearson's correlation coefficient showed positive and significant association between aliphatic GSL compounds and expression values of ST5c-Bol030757 and AOP2-Bo9g006240 genes as well as between indolic GSL compounds and the expression of MYB34-Bol017062, MYB122-Bol026204, CYP81F2-Bol012237, CYP81F4-Bol032712 and CYP81F4-Bol032714 genes. This study helps in understanding the role of GSL biosynthesis and breakdown related genes for resistance against black rot pathogen in cabbage, which could be further confirmed through functional characterization either by overexpression or knock-out mutation.

Differential expression of major genes involved in the biosynthesis of aliphatic glucosinolates in intergeneric Baemoochae (Brassicaceae) and its parents during development

Thus study found the temporal and spatial relationship between production of aliphatic glucosinolate compounds and the expression profile of glucosinolate-related genes during growth and development in radish, Chinese cabbage, and their intergeneric hybrid baemoochae plants. Glucosinolates (GSLs) are one of major bioactive compounds in Brassicaceae plants. GSLs play a role in defense against microbes as well as chemo-preventative activity against cancer, which draw attentions from plant scientists. We investigated the temporal relationship between production of aliphatic Glucosinolate (GSLs) compounds and the expression profile of GSL related genes during growth and development in radish, Chinese cabbage, and their intergeneric hybrid, baemoochae. Over the complete life cycle, Glucoraphasatin (GRH) and glucoraphanin (GRE) predominated in radish, whereas gluconapin (GNP), glucobrassicanapin (GBN), and glucoraphanin (GRA) abounded in Chinese cabbage. Baemoochae contained intermediate levels of all GSLs studied, indicating inheritance from both radish and Chinese cabbage. Expression patterns of BCAT4, CYP79F1, CYP83A1, UGT74B1, GRS1, FMOgs-ox1, and AOP2 genes showed a correlation to their corresponding encoded proteins in radish, Chinese cabbage, and baemoochae. Interestingly, there is a sharp change in gene expression pattern involved in side chain modification, particularly GRS1, FMOgs-ox1, and AOP2, among these plants during the vegetative and reproductive stage. For instance, the GRS1 was strongly expressed during leaf development, while both of FMOgs-ox1 and AOP2 was manifested high in floral tissues. Furthermore, expression of GRS1 gene which is responsible for GRH production was predominantly expressed in leaf tissues of radish and baemoochae, whereas it was only slightly detected in Chinese cabbage root tissue, explaining why radish has an abundance of GRH compared to other Brassica plants. Altogether, our comprehensive and comparative data proved that aliphatic GSLs biosynthesis is dynamically and precisely regulated in a tissue- and development-dependent manner in Brassicaceae family members.

Overexpression of the MYB29 transcription factor affects aliphatic glucosinolate synthesis in Brassica oleracea

KEY MESSAGE

Overexpression of BoMYB29 gene up-regulates the aliphatic glucosinolate pathway in Brassica oleracea plants increasing the production of the anti-cancer metabolite glucoraphanin, and the toxic and pungent sinigrin. Isothiocyanates, the bio-active hydrolysis products of glucosinolates, naturally produced by several Brassicaceae species, play an important role in human health and agriculture. This study aims at correlating the content of aliphatic glucosinolates to the expression of genes involved in their synthesis in Brassica oleracea, and perform functional analysis of BoMYB29 gene. To this purpose, three genotypes were used: a sprouting broccoli, a cabbage, and a wild genotype (Winspit), a high glucosinolate containing accession. Winspit showed the highest transcript level of BoMYB28, BoMYB29 and BoAOP2 genes, and BoAOP2 expression was positively correlated with that of the two MYB genes. Further analyses of the aliphatic glucosinolates also showed a positive correlation between the expression of BoAOP2 and the production of sinigrin and gluconapin in Winspit. The Winspit BoMYB29 CDS was cloned and overexpressed in Winspit and in the DH AG1012 line. Overexpressing Winspit plants produced higher quantities of alkenyl glucosinolates, such as sinigrin. Conversely, the DH AG1012 transformants showed a higher production of methylsulphinylalkyl glucosinolates, including glucoraphanin, and, despite an up-regulation of the aliphatic glucosinolate genes, no increase in alkenyl glucosinolates. The latter may be explained by the absence of a functional AOP2 gene in DH AG1012. Nevertheless, an extract of DH AG1012 lines overexpressing BoMYB29 provided a chemoprotective effect on human colon cells. This work exemplifies how the genetic diversity of B. oleracea may be used by breeders to select for higher expression of transcription factors for glucosinolate biosynthesis to improve its natural, health-promoting properties.

Comparative transcriptome analyses revealed different heat stress responses in high- and low-GS Brassica alboglabra sprouts

Background

Chinese kale (Brassica alboglabra) contains high nutritional elements and functional molecules, especially anticarcinogenic and antioxidant glucosinolates (GS), which was highly affected by environment temperature. To investigate the link of GS biosynthesis with heat stress response in Chinese kale, global transcription profiles of high-GS line (HG), low-GS line (LG), high-GS line under heat stress (HGT) and low-GS line under heat stress (LGT) were analyzed.

Results

Based on three biological replicates of each RNA sequencing data, 3901, 4062 and 2396 differentially expressed genes in HG vs HGT, LG vs LGT and HGT vs LGT were obtained, respectively. GO annotation, KEGG pathway analysis and a comprehensive analysis of DEGs showed a strong correlation between the GS biosynthesis and heat stress response. It was noticed that 11 differentially expressed genes tied to the GS biosynthesis were down-regulated, 23 heat shock transcription factors and 61 heat shock proteins were up-regulated upon the heat treatment. Another two Chinese kale varieties Cuibao and Shunbao with high- and low- GS content respectively, were used to validate the relationship of GS content and heat-response, and the results showed that high-GS content variety were more thermotolerant than the low-GS content one although GS significantly decreased in both varieties under heat stress. In addition, HSP100/ClpB, HSP90, HSP70 and sHSPs were differentially expressed in high- and low-GS varieties. Notably, HSP90 and sHSPs showed an obviously early response to heat stress than other related genes.

Conclusion

The higher heat resistance of high-GS Chinese kale and the sharp decrease of glucosinolate content under heat stress indicated a strong relationship of GS accumulation and heat stress response. Combined with the previous report on the low expression of HSP90 at elevated temperatures in GS-deficient mutant TU8 of Arabidopsis, the differential expression pattern of HSP90 in high- and low- GS varieties and its early heat response implied it might be a key regulator in GS metabolism and heat-resistance in Chinese kale.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5652-y) contains supplementary material, which is available to authorized users.

Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus

Flowering time is an important agronomic trait that is highly influenced by the environment. To elucidate the genetic mechanism of flowering time in rapeseed (Brassica napus L.), a genome-wide QTL analysis was performed in a doubled haploid population grown in winter, semi-winter and spring ecological conditions. Fifty-five consensus QTLs were identified after combining phenotype and genomic data, including 12 environment-stable QTLs and 43 environment-specific QTLs. Importantly, six major QTLs for flowering time were identified, of which two were considered environment-specific QTLs in spring ecological condition and four were considered environment-stable QTLs in winter and semi-winter ecological conditions. Through QTL comparison, 18 QTLs were colocalized with QTLs from six other published studies. Combining the candidate genes with their functional annotation, in 49 of 55 consensus QTLs, 151 candidate genes in B. napus corresponding to 95 homologous genes in Arabidopsis thaliana related to flowering were identified, including BnaC03g32910D (CO), BnaA02g12130D (FT) and BnaA03g13630D (FLC). Most of the candidate genes were involved in different flowering regulatory pathways. Based on re-sequencing and differences in sequence annotation between the two parents, we found that regions containing some candidate genes have numerous non-frameshift InDels and many non- synonymous mutations, which might directly lead to gene functional variation. Flowering time was negativly correlated with seed yield and thousand seed weight based on a QTL comparison of flowering time and seed yield traits, which has implications in breeding new early-maturing varieties of B. napus. Moreover, a putative flowering regulatory network was constructed, including the photoperiod, circadian clock, vernalization, autonomous and gibberellin pathways. Multiple copies of genes led to functional difference among the different copies of homologous genes, which also increased the complexity of the flowering regulatory networks. Taken together, the present results not only provide new insights into the genetic regulatory network underlying the control of flowering time but also improve our understanding of flowering time regulatory pathways in rapeseed.

Functional Characterization of BoaMYB51s as Central Regulators of Indole Glucosinolate Biosynthesis in Brassica oleracea var. alboglabra Bailey

R2R3-MYB transcription factor MYB51 is known to control indole glucosinolate (indole GSL) biosynthesis in Arabidopsis. Here, two copies of BoaMYB51 have been isolated in Chinese kale (Brassica oleracea var. alboglabra Bailey), designated BoaMYB51.1 and BoaMYB51.2, which exhibit overlapping but distinct expression levels among different organs and respond to signaling molecules in a similar pattern. It has been demonstrated a structural and functional conservation between BoaMYB51s and AtMYB51 by phylogenetic analysis, complementation studies and transient expression assay. To further investigate the transcriptional mechanism, we identified the transcriptional activation domain (TAD) and putative interacting proteins of BoaMYB51s by means of yeast (Saccharomyces cerevisiae) two hybrid. Using tobacco (Nicotiana benthamiana) transient expression assay, we confirmed that the carboxy-end is required for transcriptional activation activity of BoaMYB51s. In addition, several BoaMYB51-interacting proteins have been identified by yeast two-hybrid screening. These results provide important insights into the molecular mechanisms by which MYB51 transcriptionally regulates indole GSL biosynthesis.

A naturally occurring variation in the BrMAM-3 gene is associated with aliphatic glucosinolate accumulation in Brassica rapa leaves

Glucosinolate profiles significantly vary among Brassica rapa genotypes. However, the molecular basis of these variations is largely unknown. In this study, we investigated a major quantitative trait locus (QTL) controlling aliphatic glucosinolate accumulation in B. rapa leaves. The QTL, which encompasses three tandem MAM genes and two MYB genes, was detected in two BC₂DH populations. Among the five-candidate genes, only the expression level of BrMAM-3 (Bra013007) was significantly correlated with the accumulation of aliphatic glucosinolates in B. rapa leaves. We identified a naturally occurring insertion within exon 1 of BrMAM-3, which is predicted to be a loss-of-function mutation, as confirmed by qRT-PCR. We determined that the loss of function was associated with the low glucosinolate content in B. rapa accessions. Furthermore, overexpressing the BrMAM-3 gene resulted in an increase in total aliphatic glucosinolates in Arabidopsis transgenic lines. Our study provides insights into the molecular mechanism underlying the accumulation of aliphatic glucosinolates in B. rapa leaves, thereby facilitating in the manipulation of total aliphatic glucosinolate content in Brassica crops.

De novo Transcriptome Assembly of Chinese Kale and Global Expression Analysis of Genes Involved in Glucosinolate Metabolism in Multiple Tissues

Chinese kale, a vegetable of the cruciferous family, is a popular crop in southern China and Southeast Asia due to its high glucosinolate content and nutritional qualities. However, there is little research on the molecular genetics and genes involved in glucosinolate metabolism and its regulation in Chinese kale. In this study, we sequenced and characterized the transcriptomes and expression profiles of genes expressed in 11 tissues of Chinese kale. A total of 216 million 150-bp clean reads were generated using RNA-sequencing technology. From the sequences, 98,180 unigenes were assembled for the whole plant, and 49,582~98,423 unigenes were assembled for each tissue. Blast analysis indicated that a total of 80,688 (82.18%) unigenes exhibited similarity to known proteins. The functional annotation and classification tools used in this study suggested that genes principally expressed in Chinese kale, were mostly involved in fundamental processes, such as cellular and molecular functions, the signal transduction, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in various tissues of Chinese kale. A large number of candidate genes involved in glucosinolate metabolism and its regulation were identified, and the expression patterns of these genes were analyzed. We found that most of the genes involved in glucosinolate biosynthesis were highly expressed in the root, petiole, and in senescent leaves. The expression patterns of ten glucosinolate biosynthetic genes from RNA-seq were validated by quantitative RT-PCR in different tissues. These results provided an initial and global overview of Chinese kale gene functions and expression activities in different tissues.

Enriching Glucoraphanin in Brassica rapa Through Replacement of BrAOP2.2/BrAOP2.3 with Non-functional Genes

Sulforaphane, the hydrolytic product of glucoraphanin glucosinolate, is a potent anticarcinogen that reduces the risk of several human cancers. However, in most B. rapa vegetables, glucoraphanin is undetectable or only present in trace amounts, since the glucoraphanin that is present is converted to gluconapin by three functional BrAOP2 genes. In this study, to enrich beneficial glucoraphanin content in B. rapa, the functional BrAOP2 alleles were replaced by non-functional counterparts through marker-assisted backcrossing (MAB). We identified non-functional mutations of two BrAOP2 genes from B. rapa. The backcross progenies with introgression of both non-functional braop2.2 and braop2.3 alleles significantly increased the glucoraphanin content by 18 times relative to the recurrent parent. In contrast, replacement or introgression of single non-functional braop2.2 or braop2.3 locus did not change glucoraphanin content. Our results suggest that replacement of these two functional BrAOP2 genes with non-functional alleles has the potential for producing improved Brassica crops with enriched beneficial glucoraphanin content.

Fire and Brimstone: Molecular Interactions between Sulfur and Glucosinolate Biosynthesis in Model and Crop Brassicaceae

Glucosinolates (GSLs) represent one of the most widely studied classes of plant secondary metabolite, and have a wide range of biological activities. Their unique properties also affect livestock and human health, and have been harnessed for food and other end-uses. Since GSLs are sulfur (S)-rich there are many lines of evidence suggesting that plant S status plays a key role in determining plant GSL content. However, there is still a need to establish a detailed knowledge of the distribution and remobilization of S and GSLs throughout the development of Brassica crops, and to represent this in terms of primary and secondary sources and sinks. The increased genome complexity, gene duplication and divergence within brassicas, together with their ontogenetic plasticity during crop development, appear to have a marked effect on the regulation of S and GSLs. Here, we review the current understanding of inorganic S (sulfate) assimilation into organic S forms, including GSLs and their precursors, the intracellular and inter-organ transport of inorganic and organic S forms, and the accumulation of GSLs in specific tissues. We present this in the context of overlapping sources and sinks, transport processes, signaling molecules and their associated molecular interactions. Our analysis builds on recent insights into the molecular regulation of sulfate uptake and transport by different transporters, transcription factors and miRNAs, and the role that these may play in GSL biosynthesis. We develop a provisional model describing the key processes that could be targeted in crop breeding programs focused on modifying GSL content.

The pangenome of an agronomically important crop plant Brassica oleracea

There is an increasing awareness that as a result of structural variation, a reference sequence representing a genome of a single individual is unable to capture all of the gene repertoire found in the species. A large number of genes affected by presence/absence and copy number variation suggest that it may contribute to phenotypic and agronomic trait diversity. Here we show by analysis of the Brassica oleracea pangenome that nearly 20% of genes are affected by presence/absence variation. Several genes displaying presence/absence variation are annotated with functions related to major agronomic traits, including disease resistance, flowering time, glucosinolate metabolism and vitamin biosynthesis.

Brassica oleracea is a single species that includes diverse crops such as cabbage, broccoli and Brussels sprouts. Here, the authors identify genes not captured in existing B. oleracea reference genomes by the assembly of a pangenome and show variations in gene content that may be related to important agronomic traits