Subclade of flavin-monooxygenases involved in aliphatic glucosinolate biosynthesis

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.

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

Unravelling Glucoraphanin and Glucoerucin Metabolism across Broccoli Sprout Development: Insights from Metabolite and Transcriptome Analysis

Variations in the concentration of glucoraphanin (GRA) and glucoerucin (GER), as well as the corresponding breakdown products, isothiocyanates (ITCs) and nitriles, were investigated during the growth of broccoli sprouts. The concentrations of GRA and GER decreased sharply from 33.66 µmol/g to 11.48 µmol/g and 12.98 µmol/g to 8.23 µmol/g, respectively, after seed germination. From the third to the seventh day, both GRA and GER were maintained as relatively stable. The highest concentrations of sulforaphane (17.16 µmol/g) and erucin (12.26 µmol/g) were observed on the first day. Hereafter, the concentrations of nitrile hydrolyzed from GRA or GER were higher than those of the corresponding ITCs. Moreover, the ratio of sulforaphane to sulforaphane nitrile decreased from 1.35 to 0.164 from 1 d to 5 d, with a similar trend exhibited for erucin/erucin nitrile after 2 d. RNA-seq analysis showed that BolMYB28 and BolCYP83A1, involved in aliphatic glucosinolate (GSL) biosynthesis, remained largely unexpressed until the third day. In contrast, the genes operating within the GSL-myrosinase hydrolysis pathway were highly expressed right from the beginning, with their expression levels increasing significantly after the third day. Additionally, we identified two BolESPs and six BolNSPs that might play important roles in promoting the production of nitriles during the development of broccoli sprouts.

Flavin-containing monooxygenases FMOGS-OXs integrate flowering transition and salt tolerance in Arabidopsis thaliana

Salt stress substantially leads to flowering delay. The regulation of salt-induced late flowering has been studied at the transcriptional and protein levels; however, the involvement of secondary metabolites has rarely been investigated. Here, we report that FMOGS-OXs (EC 1.14.13.237), the enzymes that catalyze the biosynthesis of glucosinolates (GSLs), promote flowering transition in Arabidopsis thaliana. It has been reported that WRKY75 is a positive regulator, and MAF4 is a negative regulator of flowering transition. The products of FMOGS-OXs, methylsulfinylalkyl GSLs (MS GSLs), facilitate flowering by inducing WRKY75 and repressing the MAS-MAF4 module. We further show that the degradation of MS GSLs is involved in salt-induced late flowering and salt tolerance. Salt stress induces the expression of myrosinase genes, resulting in the degradation of MS GSLs, thereby relieving the promotion of WRKY75 and inhibition of MAF4, leading to delayed flowering. In addition, the degradation products derived from MS GSLs enhance salt tolerance. Previous studies have revealed that FMOGS-OXs exhibit alternative catalytic activity to form trimethylamine N-oxide (TMAO) under salt stress, which activates multiple stress-related genes to promote salt tolerance. Therefore, FMOGS-OXs integrate flowering transition and salt tolerance in various ways. Our study shed light on the functional diversity of GSLs and established a connection between flowering transition, salt resistance, and GSL metabolism.

Aphid Resistance Segregates Independently of Cardenolide and Glucosinolate Content in an Erysimum cheiranthoides (Wormseed Wallflower) F2 Population

Plants in the genus Erysimum produce both glucosinolates and cardenolides as a defense mechanism against herbivory. Two natural isolates of Erysimum cheiranthoides (wormseed wallflower) differed in their glucosinolate content, cardenolide content, and their resistance to Myzus persicae (green peach aphid), a broad generalist herbivore. Both classes of defensive metabolites were produced constitutively and were not further induced by aphid feeding. To investigate the relative importance of glucosinolates and cardenolides in E. cheiranthoides defense, we generated an improved genome assembly, genetic map, and segregating F2 population. The genotypic and phenotypic analysis of the F2 plants identified quantitative trait loci, which affected glucosinolates and cardenolides, but not the aphid resistance. The abundance of most glucosinolates and cardenolides was positively correlated in the F2 population, indicating that similar processes regulate their biosynthesis and accumulation. Aphid reproduction was positively correlated with glucosinolate content. Although the overall cardenolide content had little effect on aphid growth and survival, there was a negative correlation between aphid reproduction and helveticoside abundance. However, this variation in defensive metabolites could not explain the differences in aphid growth on the two parental lines, suggesting that processes other than the abundance of glucosinolates and cardenolides have a predominant effect on aphid resistance in E. cheiranthoides.

Aphid resistance segregates independently of cardiac glycoside and glucosinolate content in an Erysimum cheiranthoides (wormseed wallflower) F2 population

Plants in the genus Erysimum produce both glucosinolates and cardiac glycosides as defense against herbivory. Two natural isolates of Erysimum cheiranthoides (wormseed wallflower) differed in their glucosinolate content, cardiac glycoside content, and resistance to Myzus persicae (green peach aphid), a broad generalist herbivore. Both classes of defensive metabolites were produced constitutively and were not induced further by aphid feeding. To investigate the relative importance of glucosinolates and cardiac glycosides in E. cheiranthoides defense, we generated an improved genome assembly, genetic map, and segregating F2 population. Genotypic and phenotypic analysis of the F2 plants identified quantitative trait loci affecting glucosinolates and cardiac glycosides, but not aphid resistance. The abundance of most glucosinolates and cardiac glycosides was positively correlated in the F2 population, indicating that similar processes regulate their biosynthesis and accumulation. Aphid reproduction was positively correlated with glucosinolate content. Although overall cardiac glycoside content had little effect on aphid growth and survival, there was a negative correlation between aphid reproduction and helveticoside abundance. However, this variation in defensive metabolites could not explain the differences in aphid growth on the two parental lines, suggesting that processes other than the abundance of glucosinolates and cardiac glycosides have a predominant effect on aphid resistance in E. cheiranthoides.

Large-scale identification of novel transcriptional regulators of the aliphatic glucosinolate pathway in Arabidopsis

Aliphatic glucosinolates are a large group of plant secondary metabolites characteristic of Brassicaceae, including the model plant Arabidopsis. The diverse and complex degradation products of aliphatic glucosinolates contribute to plant responses to herbivory, pathogen attack, and environmental stresses. Most of the biosynthesis genes in the aliphatic glucosinolate pathway have been cloned in Arabidopsis, and the research focus has recently shifted to the regulatory mechanisms controlling aliphatic glucosinolate accumulation. Up till now, more than 40 transcriptional regulators have been identified as regulating the aliphatic glucosinolate pathway, but many more novel regulators likely remain to be discovered based on research evidence over the past decade. In the current study, we took a systemic approach to functionally test 155 candidate transcription factors in Arabidopsis identified by yeast one-hybrid assay, and successfully validated at least 30 novel regulators that could significantly influence the accumulation of aliphatic glucosinolates in our experimental set-up. We also showed that the regulators of the aliphatic glucosinolate pathway have balanced positive and negative effects, and glucosinolate metabolism and plant development can be coordinated. Our work is the largest scale effort so far to validate transcriptional regulators of a plant secondary metabolism pathway, and provides new insights into how the highly diverse plant secondary metabolism is regulated at the transcriptional level.

Glucosinolate Catabolism Maintains Glucosinolate Profiles and Transport in Sulfur-Starved Arabidopsis

Glucosinolates (GSLs) are sulfur (S)-rich specialized metabolites present in Brassicales order plants. Our previous study found that GSL can function as a S source in Arabidopsis seedlings via its catabolism catalyzed by two β-glucosidases (BGLUs), BGLU28 and BGLU30. However, as GSL profiles in plants vary among growth stages and organs, the potential contribution of BGLU28/30-dependent GSL catabolism at the reproductive growth stage needs verification. Thus, in this study, we assessed growth, metabolic and transcriptional phenotypes of mature bglu28/30 double mutants grown under different S conditions. Our results showed that compared to wild-type plants grown under -S, mature bglu28/30 mutants displayed impaired growth and accumulated increased levels of GSL in their reproductive organs and rosette leaves of before-bolting plants. In contrast, the levels of primary S-containing metabolites, glutathione and cysteine decreased in their mature seeds. Furthermore, the transport of GSL from rosette leaves to the reproductive organs was stimulated in the bglu28/30 mutants under -S. Transcriptome analysis revealed that genes related to other biological processes, such as ethylene response, defense response and plant response to heat, responded differentially to -S in the bglu28/30 mutants. Altogether, these findings broadened our understanding of the roles of BGLU28/30-dependent GSL catabolism in plant adaptation to nutrient stress.

Quantitative Time-Course Analysis of Osmotic and Salt Stress in Arabidopsis thaliana Using Short Gradient Multi-CV FAIMSpro BoxCar DIA

A major limitation when undertaking quantitative proteomic time-course experimentation is the tradeoff between depth-of-analysis and speed-of-analysis. In high complexity and high dynamic range sample types, such as plant extracts, balance between resolution and time is especially apparent. To address this, we evaluate multiple compensation voltage (CV) high field asymmetric waveform ion mobility spectrometry (FAIMSpro) settings using the latest label-free single-shot Orbitrap-based DIA acquisition workflows for their ability to deeply quantify the Arabidopsis thaliana seedling proteome. Using a BoxCarDIA acquisition workflow with a -30 -50 -70 CV FAIMSpro setting, we were able to consistently quantify >5000 Arabidopsis seedling proteins over a 21-min gradient, facilitating the analysis of ∼42 samples per day. Utilizing this acquisition approach, we then quantified proteome-level changes occurring in Arabidopsis seedling shoots and roots over 24 h of salt and osmotic stress, to identify early and late stress response proteins and reveal stress response overlaps. Here, we successfully quantify >6400 shoot and >8500 root protein groups, respectively, quantifying nearly ∼9700 unique protein groups in total across the study. Collectively, we pioneer a short gradient, multi-CV FAIMSpro BoxCarDIA acquisition workflow that represents an exciting new analysis approach for undertaking quantitative proteomic time-course experimentation in plants.

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

Quantitative disease resistance (QDR) controls the association of the light leaf spot pathogen with Brassica napus; four QDR loci that were in linkage disequilibrium and eight gene expression markers were identified. Quantitative disease resistance (QDR) can provide durable control of pathogens in crops in contrast to resistance (R) gene-mediated resistance which can break down due to pathogen evolution. QDR is therefore a desirable trait in crop improvement, but little is known about the causative genes, and so it is difficult to incorporate into breeding programmes. Light leaf spot, caused by Pyrenopeziza brassicae, is an important disease of oilseed rape (canola, Brassica napus). To identify new QDR gene loci, we used a high-throughput screening pathosystem with P. brassicae on 195 lines of B. napus combined with an association transcriptomics platform. We show that all resistance against P. brassicae was associated with QDR and not R gene-mediated. We used genome-wide association analysis with an improved B. napus population structure to reveal four gene loci significantly (P = 0.0001) associated with QDR in regions showing linkage disequilibrium. On chromosome A09, enhanced resistance was associated with heterozygosity for a cytochrome P450 gene co-localising with a previously described locus for seed glucosinolate content. In addition, eight significant gene expression markers with a false discovery rate of 0.001 were associated with QDR against P. brassicae. For seven of these, expression was positively correlated with resistance, whereas for one, a HXXXD-type acyl-transferase, negative correlation indicated a potential susceptibility gene. The study identifies novel QDR loci for susceptibility and resistance, including novel cryptic QDR genes associated with heterozygosity, that will inform future crop improvement.

Formation, immunomodulatory activities, and enhancement of glucosinolates and sulforaphane in broccoli sprouts: a review for maximizing the health benefits to human

Broccoli sprouts have been considered as functional foods which have received increasing attention because they have been highly prized for glucosinolates, phenolics, and vitamins in particular glucosinolates. One of hydrolysates-sulforaphane from glucoraphanin is positively associated with the attenuation of inflammatory, which could reduce diabetes, cardiovascular and cancer risk. In recent decades, the great interest in natural bioactive components especially for sulforaphane promotes numerous researchers to investigate the methods to enhance glucoraphanin levels in broccoli sprouts and evaluate the immunomodulatory activities of sulforaphane. Therefore, glucosinolates profiles are different in broccoli sprouts varied with genotypes and inducers. Physicochemical, biological elicitors, and storage conditions were widely studied to promote the accumulation of glucosinolates and sulforaphane in broccoli sprouts. These inducers would stimulate the biosynthesis pathway gene expression and enzyme activities of glucosinolates and sulforaphane to increase the concentration in broccoli sprouts. The immunomodulatory activity of sulforaphane was summarized to be a new therapy for diseases with immune dysregulation. The perspective of this review served as a potential reference for customers and industries by application of broccoli sprouts as a functional food and clinical medicine.

Genomic Survey of Flavin Monooxygenases in Wild and Cultivated Rice Provides Insight into Evolution and Functional Diversities

The flavin monooxygenase (FMO) enzyme was discovered in mammalian liver cells that convert a carcinogenic compound, N-N'-dimethylaniline, into a non-carcinogenic compound, N-oxide. Since then, many FMOs have been reported in animal systems for their primary role in the detoxification of xenobiotic compounds. In plants, this family has diverged to perform varied functions like pathogen defense, auxin biosynthesis, and S-oxygenation of compounds. Only a few members of this family, primarily those involved in auxin biosynthesis, have been functionally characterized in plant species. Thus, the present study aims to identify all the members of the FMO family in 10 different wild and cultivated Oryza species. Genome-wide analysis of the FMO family in different Oryza species reveals that each species has multiple FMO members in its genome and that this family is conserved throughout evolution. Taking clues from its role in pathogen defense and its possible function in ROS scavenging, we have also assessed the involvement of this family in abiotic stresses. A detailed in silico expression analysis of the FMO family in Oryza sativa subsp. japonica revealed that only a subset of genes responds to different abiotic stresses. This is supported by the experimental validation of a few selected genes using qRT-PCR in stress-sensitive Oryza sativa subsp. indica and stress-sensitive wild rice Oryza nivara. The identification and comprehensive in silico analysis of FMO genes from different Oryza species carried out in this study will serve as the foundation for further structural and functional studies of FMO genes in rice as well as other crop types.

Genome-Wide Association Study of Glucosinolate Metabolites (mGWAS) in Brassica napus L

Glucosinolates (GSLs) are secondary plant metabolites that are enriched in rapeseed and related Brassica species, and they play important roles in defense due to their anti-nutritive and toxic properties. Here, we conducted a genome-wide association study of six glucosinolate metabolites (mGWAS) in rapeseed, including three aliphatic glucosinolates (m145 gluconapin, m150 glucobrassicanapin and m151 progoitrin), one aromatic glucosinolate (m157 gluconasturtiin) and two indole glucosinolates (m165 indolylmethyl glucosinolate and m172 4-hydroxyglucobrassicin), respectively. We identified 113 candidate intervals significantly associated with these six glucosinolate metabolites. In the genomic regions linked to the mGWAS peaks, 187 candidate genes involved in glucosinolate biosynthesis (e.g., BnaMAM1, BnaGGP1, BnaSUR1 and BnaMYB51) and novel genes (e.g., BnaMYB44, BnaERF025, BnaE2FC, BnaNAC102 and BnaDREB1D) were predicted based on the mGWAS, combined with analysis of differentially expressed genes. Our results provide insight into the genetic basis of glucosinolate biosynthesis in rapeseed and should facilitate marker-based breeding for improved seed quality in Brassica species.

Natural variation and artificial selection at the BnaC2.MYB28 locus modulate Brassica napus seed glucosinolate

The degradation products of glucosinolates (GSLs) greatly lower the nutritional value of rapeseed (Brassica napus) meal; thus, reduction of seed GSL content (SGC) has become an important objective of rapeseed breeding. In our previous study, we finely mapped a major QTL (qGSL-C2) for SGC to a 49-kb collinear region on B. rapa chromosome A2. Here, we experimentally validated that BnaC2.MYB28, encoding an R2R3-MYB transcription factor, is the causal gene of qGSL-C2. BnaC2.MYB28 is a nucleus-localized protein mainly expressed in vegetative tissues. Knockout of BnaC2.MYB28 in the high-SGC parent G120 reduced SGC to a value lower than that in the low-SGC parent ZY50, while overexpression of BnaC2.MYB28 in both parental lines (G120 and ZY50) led to extremely high SGC, indicating that BnaC2.MYB28 acts as a positive regulator of SGC in both parents. Molecular characterization revealed that BnaC2.MYB28 forms a homodimer and specifically interacts with BnaMYC3. Moreover, BnaC2.MYB28 can directly activate the expression of GSL biosynthesis genes. Differential expression abundance resulting from the polymorphic promoter sequences, in combination with the different capability in activating downstream genes involved in aliphatic GSL biosynthesis, caused the functional divergence of BnaC2.MYB28 in SGC regulation between the parents. Natural variation of BnaC2.MYB28 was highly associated with SGC in natural germplasm and has undergone artificial selection in modern low-GSL breeding. This study provides important insights into the core function of BnaC2.MYB28 in regulating SGC and a promising strategy for manipulating SGC in rapeseed.

Identification of a regiospecific S-oxygenase for the production of marasmin in traditional medicinal plant Tulbaghia violacea

Marasmin [S-(methylthiomethyl)-L-cysteine-4-oxide] is a pharmaceutically valuable sulfur-containing compound produced by the traditional medicinal plant, Tulbaghia violacea. Here, we report the identification of an S-oxygenase, TvMAS1, that produces marasmin from its corresponding sulfide, S-(methylthiomethyl)-L-cysteine. The amino acid sequence of TvMAS1 showed high sequence similarity to known flavin-containing S-oxygenating monooxygenases in plants. Recombinant TvMAS1 catalyzed regiospecific S-oxygenation at S4 of S-(methylthiomethyl)-L-cysteine to yield marasmin, with an apparent K _m value of 0.55 mM. TvMAS1 mRNA accumulated with S-(methylthiomethyl)-L-cysteine and marasmin in various organs of T. violacea. Our findings suggest that TvMAS1 catalyzes the S-oxygenation reaction during the last step of marasmin biosynthesis in T. violacea.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

IQD1 Involvement in Hormonal Signaling and General Defense Responses Against Botrytis cinerea

IQ Domain 1 (IQD1) is a novel Arabidopsis thaliana calmodulin-binding protein, which was found to be a positive regulator of glucosinolate (GS) accumulation and plant defense responses against insects. We demonstrate here that the IQD1 overexpressing line (IQD1 ^OXP ) was also more resistant also to the necrotrophic fungus Botrytis cinerea, whereas an IQD1 knockout line (iqd1-1) was much more sensitive. Furthermore, we showed that IQD1 is up-regulated by jasmonic acid (JA) and downregulated by salicylic acid (SA). A comparison of whole transcriptome expression between iqd1-1 and wild type plants revealed a substantial downregulation of genes involved in plant defense and hormone regulation. Further examination revealed a marked reduction of SA and increases in the levels of ethylene, JA and abscisic acid response genes in the iqd1-1 line. Moreover, quantification of SA, JA, and abscisic acids in IQD1 ^OXP and iqd1-1 lines relative to the wild type, showed a significant reduction in endogenous JA levels in the knockout line, simultaneously with increased SA levels. Relations between IQD1 ^OXP and mutants defective in plant-hormone response indicated that IQD1 cannot rescue the absence of NPR1 or impaired SA accumulation in the NahG line. IQD1 cannot rescue ein2 or eto1 mutations connected to the ethylene pathway involved in both defense responses against B. cinerea and in regulating GS accumulation. Furthermore, IQD1cannot rescue the aos, coi1 or jar1mutations, all involved in the defense response against B. cinerea and it depends on JAR1 to control indole glucosinolate accumulation. We also found that in the B. cinerea, which infected the iqd1-1 mutant, the most abundant upregulated group of proteins is involved in the degradation of complex carbohydrates, as correlated with the sensitivity of this mutant. In summary, our results suggest that IQD1 is an important A. thaliana defensive protein against B. cinerea that is integrated into several important pathways, such as those involved in plant defense and hormone responses.

The ease and complexity of identifying and using specialized metabolites for crop engineering

Plants produce a broad variety of specialized metabolites with distinct biological activities and potential applications. Despite this potential, most biosynthetic pathways governing specialized metabolite production remain largely unresolved across the plant kingdom. The rapid advancement of genetics and biochemical tools has enhanced our ability to identify plant specialized metabolic pathways. Further advancements in transgenic technology and synthetic biology approaches have extended this to a desire to design new pathways or move existing pathways into new systems to address long-running difficulties in crop systems. This includes improving abiotic and biotic stress resistance, boosting nutritional content, etc. In this review, we assess the potential and limitations for (1) identifying specialized metabolic pathways in plants with multi-omics tools and (2) using these enzymes in synthetic biology or crop engineering. The goal of these topics is to highlight areas of research that may need further investment to enhance the successful application of synthetic biology for exploiting the myriad of specialized metabolic pathways.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.)

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Analysis of Glucosinolate Content and Metabolism Related Genes in Different Parts of Chinese Flowering Cabbage

Glucosinolates (GSLs) are important secondary metabolites that play important defensive roles in cruciferous plants. Chinese flowering cabbage, one of the most common vegetable crops, is rich in GSLs and thus has the potential to reduce the risk of cancer in humans. Many genes that are involved in GSL biosynthesis and metabolism have been identified in the model plant Arabidopsis thaliana; however, few studies investigated the genes related to GSL biosynthesis and metabolism in Chinese flowering cabbage. In the present study, the GSL composition and content in three different organs of Chinese flowering cabbage (leaf, stalk, and flower bud) were determined. Our results showed that the total GSL content in flower buds was significantly higher than in stalks and leaves, and aliphatic GSLs were the most abundant GSL type. To understand the molecular mechanisms underlying the variations of GSL content, we analyzed the expression of genes encoding enzymes involved in GSL biosynthesis and transport in different tissues of Chinese flowering cabbage using RNA sequencing; the expression levels of most genes were found to be consistent with the pattern of total GSL content. Correlation and consistency analysis of differentially expressed genes from different organs with the GSL content revealed that seven genes (Bra029966, Bra012640, Bra016787, Bra011761, Bra006830, Bra011759, and Bra029248) were positively correlated with GSL content. These findings provide a molecular basis for further elucidating GSL biosynthesis and transport in Chinese flowering cabbage.

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.

Transcriptomic and metabolic analyses revealed the modulatory effect of vernalization on glucosinolate metabolism in radish (Raphanus sativus L.)

Vernalization is the process by which long-term cold like winter triggers transition to flowering in plants. Many biennial and perennial plants including Brassicaceae family plants require vernalization for floral transition. Not only floral transition, but dynamic physiological and metabolic changes might also take place during vernalization. However, vernalization-mediated metabolic change is merely investigated so far. One of secondary metabolites found in Brassiceceae family plants is glucosinolates (GSLs). GSLs provides defense against pathogens and herbivores attack in plants and also exhibits inhibitory activity against human cancer cell. Profiles of GSLs are highly modulated by different environmental stresses in Brassciaceae family plants. To grasp the effect of vernalization on GSLs metabolic dynamics in radish (Raphanus sativus L.), we performed transcriptomic and metabolic analysis during vernalization in radish. Through transcriptome analysis, we found many GSLs metabolic genes were significantly down-regulated by vernalization in radish plants. Ultra-High Performance Liquid Chromatography analysis also revealed that GSLs compounds were substantially reduced in vernalized radish samples compared to non-vernalized radish samples. Furthermore, we found that repressive histone modification (i.e. H3K27me3) is involved in the modulation of GSLs metabolism via epigenetic suppression of Glucoraphasatin Synthase 1 (GRS1) during vernalization in radish. This study revealed that GSLs metabolism is modulated by vernalization, suggestive of a newly identified target of vernalization in radish.

Identification, analysis, and modeling of the YUCCA protein family genome-wide in Coffea canephora

Auxin is involved in almost every aspect of plant growth and development, from embryogenesis to senescence. Indole-3-acetic acid (IAA) is the main known natural auxin that is synthesized by enzymes tryptophan aminotransferase of arabidopsis (TAA) and YUCCA (YUC) of the flavin-containing monooxygenases family (FMO) from one of the tryptophan-dependent pathways. Genome-wide identification and comprehensive analysis of the YUC-protein family have been conducted in Coffea canephora in the present study. A total of 10 members CcYUC gene family were identified in C. canephora. Phylogenetic analysis revealed that the CcYUC protein family is evolutionarily conserved, and they consist of four groups. In contrast, bioinformatic analysis predicted a hydrophobic transmembrane helix (TMH) for one CcYUC (YUC10) member only. Isoelectric point (pI), molecular mass (Ms), signal peptide, subcellular localization, and phosphorylation sites were predicted for CcYUC proteins. YUC enzymes require the prosthetic group flavin adenine dinucleotide (FAD) and the cofactor nicotinamide adenine dinucleotide phosphate (NADPH) for their enzymatic activity. Therefore, we include the molecular docking for CcYUC2-FAD-NADPH-IPyA and yucasin, which is a specific inhibitor for YUC activity. The docking results showed FAD and NADPH binding at the big and small domain sites, respectively, in CcYUC2. IPyA binds very close to FAD along the big domain, and yucasin competes for the same site as IPA, blocking IAA production. Furthermore, in silico point mutations affect the stability of the CcYUC2-4 proteins.

Influence of formic acid treatment on the proteome of the ectoparasite Varroa destructor

The ectoparasite Varroa destructor Anderson and Trueman is the most important parasites of the western honey bee, Apis mellifera L. The most widely currently used treatment uses formic acid (FA), but the understanding of its effects on V. destructor is limited. In order to understand the mechanism of action of FA, its effect on Varroa mites was investigated using proteomic analysis by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS). V. destructor was collected from honey bee colonies with natural mite infestation before and 24 h after the initiation of FA treatment and subjected to proteome analysis. A total of 2637 proteins were identified. Quantitative analysis of differentially expressed candidate proteins (fold change ≥ 1.5; p ≤ 0.05) revealed 205 differentially expressed proteins: 91 were induced and 114 repressed in the FA-treated group compared to the untreated control group. Impaired protein synthesis accompanied by increased protein and amino acid degradation suggest an imbalance in proteostasis. Signs of oxidative stress included significant dysregulation of candidate proteins of mitochondrial cellular respiration, increased endocytosis, and induction of heat shock proteins. Furthermore, an increased concentration of several candidate proteins associated with detoxification was observed. These results suggest dysregulated cellular respiration triggered by FA treatment as well as an increase in cellular defense mechanisms, including induced heat shock proteins and detoxification enzymes.

Structural Studies of Aliphatic Glucosinolate Chain-Elongation Enzymes

Plants evolved specialized metabolic pathways through gene duplication and functional divergence of enzymes involved in primary metabolism. The results of this process are varied pathways that produce an array of natural products useful to both plants and humans. In plants, glucosinolates are a diverse class of natural products. Glucosinolate function stems from their hydrolysis products, which are responsible for the strong flavors of Brassicales plants, such as mustard, and serve as plant defense molecules by repelling insects, fighting fungal infections, and discouraging herbivory. Additionally, certain hydrolysis products such as isothiocyanates can potentially serve as cancer prevention agents in humans. The breadth of glucosinolate function is a result of its great structural diversity, which comes from the use of aliphatic, aromatic and indole amino acids as precursors and elongation of some side chains by up to nine carbons, which, after the formation of the core glucosinolate structure, can undergo further chemical modifications. Aliphatic methionine-derived glucosinolates are the most abundant form of these compounds. Although both elongation and chemical modification of amino acid side chains are important for aliphatic glucosinolate diversity, its elongation process has not been well described at the molecular level. Here, we summarize new insights on the iterative chain-elongation enzymes methylthioalkylmalate synthase (MAMS) and isopropylmalate dehydrogenase (IPMDH).

The scope of flavin-dependent reactions and processes in the model plant Arabidopsis thaliana

Flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are utilized as coenzymes in many biochemical reduction-oxidation reactions owing to the ability of the tricyclic isoalloxazine ring system to employ the oxidized, radical and reduced state. We have analyzed the genome of Arabidopsis thaliana to establish an inventory of genes encoding flavin-dependent enzymes (flavoenzymes) as a basis to explore the range of flavin-dependent biochemical reactions that occur in this model plant. Expectedly, flavoenzymes catalyze many pivotal reactions in primary catabolism, which are connected to the degradation of basic metabolites, such as fatty and amino acids as well as carbohydrates and purines. On the other hand, flavoenzymes play diverse roles in anabolic reactions most notably the biosynthesis of amino acids as well as the biosynthesis of pyrimidines and sterols. Importantly, the role of flavoenzymes goes much beyond these basic reactions and extends into pathways that are equally crucial for plant life, for example the production of natural products. In this context, we outline the participation of flavoenzymes in the biosynthesis and maintenance of cofactors, coenzymes and accessory plant pigments (e. g. carotenoids) as well as phytohormones. Moreover, several multigene families have emerged as important components of plant immunity, for example the family of berberine bridge enzyme-like enzymes, flavin-dependent monooxygenases and NADPH oxidases. Furthermore, the versatility of flavoenzymes is highlighted by their role in reactions leading to tRNA-modifications, chromatin regulation and cellular redox homeostasis. The favorable photochemical properties of the flavin chromophore are exploited by photoreceptors to govern crucial processes of plant adaptation and development. Finally, a sequence- and structure-based approach was undertaken to gain insight into the catalytic role of uncharacterized flavoenzymes indicating their involvement in unknown biochemical reactions and pathways in A. thaliana.

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.

Trimethylamine N-oxide is a new plant molecule that promotes abiotic stress tolerance

Trimethylamine N-oxide (TMAO) is a well-known naturally occurring osmolyte in animals that counteracts the effect of different denaturants related to environmental stress and has recently been associated with severe human chronic diseases. In plants, however, the presence of TMAO has not yet been reported. In this study, we demonstrate that plants contain endogenous levels of TMAO, that it is synthesized by flavin-containing monooxygenases, and that its levels increase in response to abiotic stress conditions. In addition, our results reveal that TMAO operates as a protective osmolyte in plants, promoting appropriate protein folding and as an activator of abiotic stress-induced gene expression. Consistent with these functions, we show that TMAO enhances plant adaptation to low temperatures, drought, and high salt. We have thus uncovered a previously unidentified plant molecule that positively regulates abiotic stress tolerance.

Neighbor GWAS: incorporating neighbor genotypic identity into genome-wide association studies of field herbivory

An increasing number of field studies have shown that the phenotype of an individual plant depends not only on its genotype but also on those of neighboring plants; however, this fact is not taken into consideration in genome-wide association studies (GWAS). Based on the Ising model of ferromagnetism, we incorporated neighbor genotypic identity into a regression model, named "Neighbor GWAS". Our simulations showed that the effective range of neighbor effects could be estimated using an observed phenotype when the proportion of phenotypic variation explained (PVE) by neighbor effects peaked. The spatial scale of the first nearest neighbors gave the maximum power to detect the causal variants responsible for neighbor effects, unless their effective range was too broad. However, if the effective range of the neighbor effects was broad and minor allele frequencies were low, there was collinearity between the self and neighbor effects. To suppress the false positive detection of neighbor effects, the fixed effect and variance components involved in the neighbor effects should be tested in comparison with a standard GWAS model. We applied neighbor GWAS to field herbivory data from 199 accessions of Arabidopsis thaliana and found that neighbor effects explained 8% more of the PVE of the observed damage than standard GWAS. The neighbor GWAS method provides a novel tool that could facilitate the analysis of complex traits in spatially structured environments and is available as an R package at CRAN ( https://cran.rproject.org/package=rNeighborGWAS ).

Flavoprotein monooxygenases: Versatile biocatalysts

Flavoprotein monooxygenases (FPMOs) are single- or two-component enzymes that catalyze a diverse set of chemo-, regio- and enantioselective oxyfunctionalization reactions. In this review, we describe how FPMOs have evolved from model enzymes in mechanistic flavoprotein research to biotechnologically relevant catalysts that can be applied for the sustainable production of valuable chemicals. After a historical account of the development of the FPMO field, we explain the FPMO classification system, which is primarily based on protein structural properties and electron donor specificities. We then summarize the most appealing reactions catalyzed by each group with a focus on the different types of oxygenation chemistries. Wherever relevant, we report engineering strategies that have been used to improve the robustness and applicability of FPMOs.

Abscisic acid-mediated induction of FLAVIN-CONTAINING MONOOXYGENASE 2 leads to reduced accumulation of methylthioalkyl glucosinolates in Arabidopsis thaliana

Side-chain modification contributes to the structural diversity of aliphatic glucosinolates (GSLs), a class of sulfur-containing secondary metabolites found in Brassicales. The first step in side-chain modification of aliphatic GSLs is the S-oxygenation of the methylthioalkyl (MT) moiety to the methylsulfinylalkyl (MS) moiety. This reaction is catalyzed by flavin-containing monooxygenase (FMO_GS-OX), which is encoded by seven genes in Arabidopsis thaliana. Therefore, the regulation of FMO_GS-OX gene expression is key to controlling side-chain structural diversity. In this study, we demonstrated that the expression of FMO_GS-OX2 and FMO_GS-OX4 was induced by glucose treatment, independent of MYB28/29 and MYC2/3/4, the transcription factors that positively regulate aliphatic GSL biosynthesis. Glucose treatment of the abscisic acid (ABA)-related mutants indicated that glucose-triggered upregulation of FMO_GS-OX2 and FMO_GS-OX4 was partially regulated by ABA through the key negative regulators ABI1 and ABI2, and the positive regulator SnRK2, but not via the transcription factor ABI5. In wild-type plants, glucose treatment drastically reduced the accumulation of 4-methylthiobutyl (4MT) GSL, whereas a decrease in 4MT GSL was not observed in the fmo_gs-ox2, abi1-1, abi2-1, aba2-1, or aba3-1 mutants. This result indicated that the decreased accumulation of 4MT GSL by glucose treatment was attributed to upregulation of FMO_GS-OX2 via the ABA signaling pathway.

Evolutionary changes in the glucosinolate biosynthetic capacity in species representing Capsella, Camelina and Neslia genera

Glucosinolates are unique thioglucosides that evolved in the order Brassicales. These compounds function in plant adaptation to the environment, including combating plant pathogens, herbivore deterrence and abiotic stress tolerance. In line with their defensive functions glucosinolates usually accumulate constitutively in relatively high amounts in all tissues of Brassicaceae plants. Here we performed glucosinolate analysis in different organs of selected species representing Capsella, Camelina and Neslia genera, which similarly as the model plant Arabidopsis thaliana belong to the Camelineae tribe. We also identified orthologs of A. thaliana glucosinolate biosynthetic genes in the published genomes of some of the investigated species. Subsequent gene expression and phylogenetic analyses enabled us an insight into the evolutionary changes in the transcription of these genes and in the sequences of respective proteins that occurred within the Camelineae tribe. Our results indicated that glucosinolates are highly abundant in siliques and roots of the investigated species but hardly, if at all, produced in leaves. In addition to this unusual tissular distribution we revealed reduced structural diversity of methionine-derived aliphatic glucosinolates (AGs) with elevated accumulation of rare long chain AGs. This preference seems to correlate with evolutionary changes in genes encoding methylthioalkylmalate synthases that are responsible for the elongation of AG side chains. Finally, our results indicate that the biosynthetic pathway for tryptophan-derived indolic glucosinolates likely lost its main functions in immunity and resistance towards sucking insects and is on its evolutionary route to be shut off in the investigated species.

A chromosome-level reference genome of non-heading Chinese cabbage [Brassica campestris (syn. Brassica rapa) ssp. chinensis]

Non-heading Chinese cabbage (NHCC) is an important leafy vegetable cultivated worldwide. Here, we report the first high-quality, chromosome-level genome of NHCC001 based on PacBio, Hi-C, and Illumina sequencing data. The assembled NHCC001 genome is 405.33 Mb in size with a contig N50 of 2.83 Mb and a scaffold N50 of 38.13 Mb. Approximately 53% of the assembled genome is composed of repetitive sequences, among which long terminal repeats (LTRs, 20.42% of the genome) are the most abundant. Using Hi-C data, 97.9% (396.83 Mb) of the sequences were assigned to 10 pseudochromosomes. Genome assessment showed that this B. rapa NHCC001 genome assembly is of better quality than other currently available B. rapa assemblies and that it contains 48,158 protein-coding genes, 99.56% of which are annotated in at least one functional database. Comparative genomic analysis confirmed that B. rapa NHCC001 underwent a whole-genome triplication (WGT) event shared with other Brassica species that occurred after the WGD events shared with Arabidopsis. Genes related to ascorbic acid metabolism showed little variation among the three B. rapa subspecies. The numbers of genes involved in glucosinolate biosynthesis and catabolism were higher in NHCC001 than in Chiifu and Z1, due primarily to tandem duplication. The newly assembled genome will provide an important resource for research on B. rapa, especially B. rapa ssp. chinensis.

Glucosinolates: Natural Occurrence, Biosynthesis, Accessibility, Isolation, Structures, and Biological Activities

Glucosinolates (GSLs) are secondary plant metabolites abundantly found in plant order Brassicales. GSLs are constituted by an S-β-d-glucopyrano unit anomerically connected to O-sulfated (Z)-thiohydroximate moiety. The side-chain of the O-sulfate thiohydroximate moiety, which is derived from a different amino acid, contributes to the diversity of natural GSL, with more than 130 structures identified and validated to this day. Both the structural diversity of GSL and their biological implication in plants have been biochemically studied. Although chemical syntheses of GSL have been devised to give access to these secondary metabolites, direct extraction from biomass remains the conventional method to isolate natural GSL. While intact GSLs are biologically inactive, various products, including isothiocyanates, nitriles, epithionitriles, and cyanides obtained through their hydrolysis of GSLs, exhibit many different biological activities, among which several therapeutic benefits have been suggested. This article reviews natural occurrence, accessibility via chemical, synthetic biochemical pathways of GSL, and the current methodology of extraction, purification, and characterization. Structural information, including the most recent classification of GSL, and their stability and storage conditions will also be discussed. The biological perspective will also be explored to demonstrate the importance of these prominent metabolites.

flasher, a novel mutation in a glucosinolate modifying enzyme, conditions changes in plant architecture and hormone homeostasis

Meristem function is underpinned by numerous genes that affect hormone levels, ultimately controlling phyllotaxy, the transition to flowering and general growth properties. Class I KNOX genes are major contributors to this process, promoting cytokinin biosynthesis but repressing gibberellin production to condition a replication competent state. We identified a suppressor mutant of the KNOX1 mutant brevipedicellus (bp) that we termed flasher (fsh), which promotes stem and pedicel elongation, suppresses early senescence, and negatively affects reproductive development. Map-based cloning and complementation tests revealed that fsh is due to an E40K change in the flavin monooxygenase GS-OX5, a gene encoding a glucosinolate (GSL) modifying enzyme. In vitro enzymatic assays revealed that fsh poorly converts substrate to product, yet the levels of several GSLs are higher in the suppressor line, implicating FSH in feedback control of GSL flux. FSH is expressed predominantly in the vasculature in patterns that do not significantly overlap those of BP, implying a non-cell autonomous mode of meristem control via one or more GSL metabolites. Hormone analyses revealed that cytokinin levels are low in bp, but fsh restores cytokinin levels to near normal by activating cytokinin biosynthesis genes. In addition, jasmonate levels in the fsh suppressor are significantly lower than in bp, which is likely due to elevated expression of JA inactivating genes. These observations suggest the involvement of the GSL pathway in generating one or more negative effectors of growth that influence inflorescence architecture and fecundity by altering the balance of hormonal regulators.

The Flavoproteome of the Model Plant Arabidopsis thaliana

Flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are essential cofactors for enzymes, which catalyze a broad spectrum of vital reactions. This paper intends to compile all potential FAD/FMN-binding proteins encoded by the genome of Arabidopsis thaliana. Several computational approaches were applied to group the entire flavoproteome according to (i) different catalytic reactions in enzyme classes, (ii) the localization in subcellular compartments, (iii) different protein families and subclasses, and (iv) their classification to structural properties. Subsequently, the physiological significance of several of the larger flavoprotein families was highlighted. It is conclusive that plants, such as Arabidopsis thaliana, use many flavoenzymes for plant-specific and pivotal metabolic activities during development and for signal transduction pathways in response to biotic and abiotic stress. Thereby, often two up to several homologous genes are found encoding proteins with high protein similarity. It is proposed that these gene families for flavoproteins reflect presumably their need for differential transcriptional control or the expression of similar proteins with modified flavin-binding properties or catalytic activities.

A Comprehensive Gene Inventory for Glucosinolate Biosynthetic Pathway in Arabidopsis thaliana

Glucosinolates (GSLs) are plant secondary metabolites comprising sulfur and nitrogen mainly found in plants from the order of Brassicales, such as broccoli, cabbage, and Arabidopsis thaliana. The activated forms of GSL play important roles in fighting against pathogens and have health benefits to humans. The increasing amount of data on A. thaliana generated from various omics technologies can be investigated more deeply in search of new genes or compounds involved in GSL biosynthesis and metabolism. This review describes a comprehensive inventory of A. thaliana GSLs identified from published literature and databases such as KNApSAcK, KEGG, and AraCyc. A total of 113 GSL genes encoding for 23 transcription components, 85 enzymes, and five protein transporters were experimentally characterized in the past two decades. Continuous efforts are still on going to identify all molecules related to the production of GSLs. A manually curated database known as SuCCombase (http://plant-scc.org) was developed to serve as a comprehensive GSL inventory. Realizing lack of information on the regulation of GSL biosynthesis and degradation mechanisms, this review also includes relevant information and their connections with crosstalk among various factors, such as light, sulfur metabolism, and nitrogen metabolism, not only in A. thaliana but also in other crucifers.

Structure and function of a flavin-dependent S-monooxygenase from garlic (Allium sativum)

Allicin is a component of the characteristic smell and flavor of garlic (Allium sativum). A flavin-containing monooxygenase (FMO) produced by A. sativum (AsFMO) was previously proposed to oxidize S-allyl-l-cysteine (SAC) to alliin, an allicin precursor. Here, we present a kinetic and structural characterization of AsFMO that suggests a possible contradiction to this proposal. Results of steady-state kinetic analyses revealed that AsFMO exhibited negligible activity with SAC; however, the enzyme was highly active with l-cysteine, N-acetyl-l-cysteine, and allyl mercaptan. We found that allyl mercaptan with NADPH was the preferred substrate-cofactor combination. Rapid-reaction kinetic analyses showed that NADPH binds tightly (K_D of ∼2 μm) to AsFMO and that the hydride transfer occurs with pro-R stereospecificity. We detected the formation of a long-wavelength band when AsFMO was reduced by NADPH, probably representing the formation of a charge-transfer complex. In the absence of substrate, the reduced enzyme, in complex with NADP⁺, reacted with oxygen and formed an intermediate with a spectrum characteristic of C4a-hydroperoxyflavin, which decays several orders of magnitude more slowly than the k _cat The presence of substrate enhanced C4a-hydroperoxyflavin formation and, upon hydroxylation, oxidation occurred with a rate constant similar to the k _cat The structure of AsFMO complexed with FAD at 2.08-Å resolution features two domains for binding of FAD and NADPH, representative of class B flavin monooxygenases. These biochemical and structural results are consistent with AsFMO being an S-monooxygenase involved in allicin biosynthesis through direct formation of sulfenic acid and not SAC oxidation.

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome-wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue-specific associated loci were related to total leaf and seed GSL content. Thirty-six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double-low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Plants are often attacked by insects and other herbivores. As a result, they have evolved to defend themselves by producing many different chemicals that are toxic to these pests. As producing each chemical costs energy, individual plants often only produce one type of chemical that is targeted towards their main herbivore. Related species of plants often use the same type of chemical defense so, if a particular herbivore gains the ability to cope with this chemical, it may rapidly become an important pest for the whole plant family.

To escape this threat, some plants have gained the ability to produce more than one type of chemical defense. Wallflowers, for example, are a group of plants in the mustard family that produce two types of toxic chemicals: mustard oils, which are common in most plants in this family; and cardenolides, which are an innovation of the wallflowers, and which are otherwise found only in distantly related plants such as foxglove and milkweed. The combination of these two chemical defenses within the same plant may have allowed the wallflowers to escape attacks from their main herbivores and may explain why the number of wallflower species rapidly increased within the last two million years.

Züst et al. have now studied the diversity of mustard oils and cardenolides present in many different species of wallflower. This analysis revealed that almost all of the tested wallflower species produced high amounts of both chemical defenses, while only one species lacked the ability to produce cardenolides. The levels of mustard oils had no relation to the levels of cardenolides in the tested species, which suggests that the regulation of these two defenses is not linked. Furthermore, Züst et al. found that closely related wallflower species produced more similar cardenolides, but less similar mustard oils, to each other. This suggests that mustard oils and cardenolides have evolved independently in wallflowers and have distinct roles in the defense against different herbivores.

The evolution of insect resistance to pesticides and other toxins is an important concern for agriculture. Applying multiple toxins to crops at the same time is an important strategy to slow the evolution of resistance in the pests. The findings of Züst et al. describe a system in which plants have naturally evolved an equivalent strategy to escape their main herbivores. Understanding how plants produce multiple chemical defenses, and the costs involved, may help efforts to breed crop species that are more resistant to herbivores and require fewer applications of pesticides.

The “Green” FMOs: Diversity, Functionality and Application of Plant Flavoproteins

Flavin-dependent monooxygenases (FMOs) are ancient enzymes present in all kingdoms of life. FMOs typically catalyze the incorporation of an oxygen atom from molecular oxygen into small molecules. To date, the majority of functional characterization studies have been performed on mammalian, fungal and bacterial FMOs, showing that they play fundamental roles in drug and xenobiotic metabolism. By contrast, our understanding of FMOs across the plant kingdom is very limited, despite plants possessing far greater FMO diversity compared to both bacteria and other multicellular organisms. Here, we review the progress of plant FMO research, with a focus on FMO diversity and functionality. Significantly, of the FMOs characterized to date, they all perform oxygenation reactions that are crucial steps within hormone metabolism, pathogen resistance, signaling and chemical defense. This demonstrates the fundamental role FMOs have within plant metabolism, and presents significant opportunities for future research pursuits and downstream applications.

Expression profiles of glucosinolate biosynthetic genes in turnip (Brassica rapa var. rapa) at different developmental stages and effect of transformed flavin-containing monooxygenase genes on hairy root glucosinolate content

BACKGROUND

Glucosinolates (GSLs) are secondary metabolites, mainly existing in Brassica vegetables. Their breakdown products have health benefits and contribute to the distinctive taste of these vegetables. Because of their high value, there is a lot of interest in developing breeding strategies to increase the content of beneficial GSLs in Brassica species. GSLs are synthesized from certain amino acids and their biological roles depend largely on the structure of their side chains. Flavin-containing monooxygenase (FMO_GS-OX ) genes are involved in the synthesis of these side chains. To better understand GSL biosynthesis, we sequenced the transcriptomes of turnip (Brassica rapa var. rapa) tubers at four developmental stages (S1-S4) and determined their GSL content.

RESULTS

The total GSL content was high at the early stage (S1) of tuber development and increased up to S3, then decreased at S4. We detected 61 differentially expressed genes, including five FMO_GS-OX genes, that were related for GSL biosynthesis among the four developmental stages. Most of these genes were highly expressed at stages S1 to S3, but their expression was much lower at S4. We estimated the effect of the five FMO_GS-OX genes on GSL content by overexpressing them in turnip hairy roots and found that the amount of aliphatic GSLs increased significantly in the transgenic plants.

CONCLUSION

The transcriptome data and characterization of genes involved in GSL biosynthesis, particularly the FMO_GS-OX genes, will be valuable for improving the yield of beneficial GSLs in turnip and other Brassica crops. © 2019 Society of Chemical Industry.

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.

QTL and candidate genes associated with leaf anion concentrations in response to phosphate supply in Arabidopsis thaliana

BACKGROUND

Phosphorus is often present naturally in the soil as inorganic phosphate, Pi, which bio-availability is limited in many ecosystems due to low soil solubility and mobility. Plants respond to low Pi with a Pi Starvation Response, involving Pi sensing and long-distance signalling. There is extensive cross-talk between Pi homeostasis mechanisms and the homeostasis mechanism for other anions in response to Pi availability.

RESULTS

Recombinant Inbred Line (RIL) and Genome Wide Association (GWA) mapping populations, derived from or composed of natural accessions of Arabidopsis thaliana, were grown under sufficient and deficient Pi supply. Significant treatment effects were found for all traits and significant genotype x treatment interactions for the leaf Pi and sulphate concentrations. Using the RIL/QTL population, we identified 24 QTLs for leaf concentrations of Pi and other anions, including a major QTL for leaf sulphate concentration (SUL2) mapped to the bottom of chromosome (Chr) 1. GWA mapping found 188 SNPs to be associated with the measured traits, corresponding to 152 genes. One of these SNPs, associated with leaf Pi concentration, mapped to PP2A-1, a gene encoding an isoform of the catalytic subunit of a protein phosphatase 2A. Of two additional SNPs, associated with phosphate use efficiency (PUE), one mapped to AT5G49780, encoding a leucine-rich repeat protein kinase involved in signal transduction, and the other to SIZ1, a gene encoding a SUMO E3 ligase, and a known regulator of P starvation-dependent responses. One SNP associated with leaf sulphate concentration was found in SULTR2;1, encoding a sulphate transporter, known to enhance sulphate translocation from root to shoot under P deficiency. Finally, one SNP was mapped to FMO GS-OX4, a gene encoding glucosinolate S-oxygenase involved in glucosinolate biosynthesis, which located within the confidence interval of the SUL2 locus.

CONCLUSION

We identified several candidate genes with known functions related to anion homeostasis in response to Pi availability. Further molecular studies are needed to confirm and validate these candidate genes and understand their roles in examined traits. Such knowledge will contribute to future breeding for improved crop PUE .

Plant Networks as Traits and Hypotheses: Moving Beyond Description

Biology relies on the central thesis that the genes in an organism encode molecular mechanisms that combine with stimuli and raw materials from the environment to create a final phenotypic expression representative of the genomic programming. While conceptually simple, the genotype-to-phenotype linkage in a eukaryotic organism relies on the interactions of thousands of genes and an environment with a potentially unknowable level of complexity. Modern biology has moved to the use of networks in systems biology to try to simplify this complexity to decode how an organism's genome works. Previously, biological networks were basic ways to organize, simplify, and analyze data. However, recent advances are allowing networks to move beyond description and become phenotypes or hypotheses in their own right. This review discusses these efforts, like mapping responses across biological scales, including relationships among cellular entities, and the direct use of networks as traits or hypotheses.

Mechanisms underlying the enhanced biomass and abiotic stress tolerance phenotype of an Arabidopsis MIOX over-expresser

Myo-inositol oxygenase (MIOX) is the first enzyme in the inositol route to ascorbate (L-ascorbic acid, AsA, vitamin C). We have previously shown that Arabidopsis plants constitutively expressing MIOX have elevated foliar AsA content and displayed enhanced growth rate, biomass accumulation, and increased tolerance to multiple abiotic stresses. In this work, we used a combination of transcriptomics, chromatography, microscopy, and physiological measurements to gain a deeper understanding of the underlying mechanisms mediating the phenotype of the AtMIOX4 line. Transcriptomic analysis revealed increased expression of genes involved in auxin synthesis, hydrolysis, transport, and metabolism, which are supported by elevated auxin levels both in vitro and in vivo, and confirmed by assays demonstrating their effect on epidermal cell elongation in the AtMIOX4 over-expressers. Additionally, we detected up-regulation of transcripts involved in photosynthesis and this was validated by increased efficiency of the photosystem II and proton motive force. We also found increased expression of amylase leading to higher intracellular glucose levels. Multiple gene families conferring plants tolerance/expressed in response to cold, water limitation, and heat stresses were found to be elevated in the AtMIOX4 line. Interestingly, the high AsA plants also displayed up-regulation of transcripts and hormones involved in defense including jasmonates, defensin, glucosinolates, and transcription factors that are known to be important for biotic stress tolerance. These results overall indicate that elevated levels of auxin and glucose, and enhanced photosynthetic efficiency in combination with up-regulation of abiotic stresses response genes underly the higher growth rate and abiotic stresses tolerance phenotype of the AtMIOX4 over-expressers.

S-Alk(en)ylcysteine sulfoxides in the genus Allium: proposed biosynthesis, chemical conversion, and bioactivities

S-Alk(en)ylcysteine sulfoxides are sulfur-containing natural products characteristic of the genus Allium. Both the flavor and medicinal properties of Allium plants are attributed to a wide variety of sulfur-containing compounds that are generated from S-alk(en)ylcysteine sulfoxides. Previous radiotracer experiments proposed that S-alk(en)ylcysteine sulfoxides are biosynthesized from glutathione. The recent identification of γ-glutamyl transpeptidases and a flavin-containing S-oxygenase involved in the biosynthesis of S-allylcysteine sulfoxide (alliin) in garlic (Allium sativum) provided insights into the reaction order of deglutamylation and S-oxygenation together with the localization of the biosynthesis, although the rest of the enzymes in the pathway still await discovery. In intact plants, S-alk(en)ylcysteine sulfoxides are stored in the cytosol of storage mesophyll cells. During tissue damage, the vacuolar enzyme alliinase contacts and hydrolyzes S-alk(en)ylcysteine sulfoxides to produce the corresponding sulfenic acids, which are further converted into various sulfur-containing bioactive compounds mainly via spontaneous reactions. The formed sulfur-containing compounds exhibit bioactivities related to pathogen defense, the prevention and alleviation of cancer and cardiovascular diseases, and neuroprotection. This review summarizes the current understanding of the occurrence, biosynthesis, and alliinase-triggered chemical conversion of S-alk(en)ylcysteine sulfoxides in Allium plants as well as the impact of S-alk(en)ylcysteine sulfoxides and their derivatives on medicinal, food, and agricultural sciences.

Evolution of Glucosinolate Diversity via Whole-Genome Duplications, Gene Rearrangements, and Substrate Promiscuity

Over several decades, glucosinolates have become a model system for the study of specialized metabolic diversity in plants. The near-complete identification of biosynthetic enzymes, regulators, and transporters has provided support for the role of gene duplication and subsequent changes in gene expression, protein function, and substrate specificity as the evolutionary bases of glucosinolate diversity. Here, we provide examples of how whole-genome duplications, gene rearrangements, and substrate promiscuity potentiated the evolution of glucosinolate biosynthetic enzymes, regulators, and transporters by natural selection. This in turn may have led to the repeated evolution of glucosinolate metabolism and diversity in higher plants.

Transcriptome reveals the gene expression patterns of sulforaphane metabolism in broccoli florets

Sulforaphane is a new and effective anti-cancer component that is abundant in broccoli. In the past few years, the patterns of variability in glucosinolate content and its regulation in A. thaliana have been described in detail. However, the diversity of glucosinolate and sulforaphane contents in different organs during vegetative and reproductive stages has not been clearly explained. In this paper, we firstly investigated the transcriptome profiles of the developing buds and leaves at bolting stage of broccoli (B52) to further assess the gene expression patterns involved in sulforaphane synthesis. The CYP79F1 gene, as well as nine other genes related to glucorahpanin biosynthesis, MAM1, MAM3, St5b-2, FMO GS-OX1, MY, AOP2, AOP3, ESP and ESM1 were selected by digital gene expression analysis and were validated by quantitative real-time PCR (qRT-PCR). Meanwhile, the compositions of glucosinolates and sulforaphane were detected for correlation analysis with related genes. Finally the RNA sequencing libraries generated 147 957 344 clean reads, and 8 539 unigene assemblies were produced. In digital result, only CYP79F1, in the glucoraphanin pathway, was up-regulated in young buds but absent from the other organs, which was consistent with the highest level of sulforaphane content being in this organ compared to mature buds, buds one day before flowering, flowers and leaves. The sequencing results also presented that auxin and cytokinin might affect glucoraphanin accumulation. The study revealed that up-regulated expression of CYP79F1 plays a fundamental and direct role in sulforaphane production in inflorescences. Two genes of MAM1 and St5b-2 could up-regulated glucoraphanin generation. Synergistic expression of MAM1, MAM3, St5b-2, FMO GS-OX1, MY, ESP and ESM1 was found in sulforaphane metabolism. This study will be beneficial for understanding the diversity of sulforaphane in broccoli organs.